Genome-wide association studies (GWASs) have identified more than 150 loci associated with blood lipid and cholesterol levels; however, the functional and molecular mechanisms for many associations are unknown. We examined the functional regulatory effects of candidate variants at the GALNT2 locus associated with high-density lipoprotein cholesterol (HDL-C). Fine-mapping and conditional analyses in the METSIM study identified a single locus harboring 25 noncoding variants (r(2) > 0.7 with the lead GWAS variants) strongly associated with total cholesterol in medium-sized HDL (e.g., rs17315646, p = 3.5 × 10(-12)). We used luciferase reporter assays in HepG2 cells to test all 25 variants for allelic differences in regulatory enhancer activity. rs2281721 showed allelic differences in transcriptional activity (75-fold [T] versus 27-fold [C] more than the empty-vector control), as did a separate 780-bp segment containing rs4846913, rs2144300, and rs6143660 (49-fold [AT(-) haplotype] versus 16-fold [CC(+) haplotype] more). Using electrophoretic mobility shift assays, we observed differential CEBPB binding to rs4846913, and we confirmed this binding in a native chromatin context by performing chromatin-immunoprecipitation (ChIP) assays in HepG2 and Huh-7 cell lines of differing genotypes. Additionally, sequence reads in HepG2 DNase-I-hypersensitivity and CEBPB ChIP-seq signals spanning rs4846913 showed significant allelic imbalance. Allelic-expression-imbalance assays performed with RNA from primary human hepatocyte samples and expression-quantitative-trait-locus (eQTL) data in human subcutaneous adipose tissue samples confirmed that alleles associated with increased HDL-C are associated with a modest increase in GALNT2 expression. Together, these data suggest that at least rs4846913 and rs2281721 play key roles in influencing GALNT2 expression at this HDL-C locus.

Genome-wide association studies of gene expression traits and other cellular phenotypes have successfully identified links between genetic variation and biological processes. The majority of discoveries have uncovered cis-expression quantitative trait locus (eQTL) effects via mass univariate testing of SNPs against gene expression in single tissues. Here we present a Bayesian method for multiple-tissue experiments focusing on uncovering gene networks linked to genetic variation. Our method decomposes the 3D array (or tensor) of gene expression measurements into a set of latent components. We identify sparse gene networks that can then be tested for association against genetic variation across the genome. We apply our method to a data set of 845 individuals from the TwinsUK cohort with gene expression measured via RNA-seq analysis in adipose, lymphoblastoid cell lines (LCLs) and skin. We uncover several gene networks with a genetic basis and clear biological and statistical significance. Extensions of this approach will allow integration of different omics, environmental and phenotypic data sets.

Glycated haemoglobin (HbA1c) is recommended as an additional tool to glucose-based measures (fasting plasma glucose [FPG] and 2-hour plasma glucose [2PG] during oral glucose tolerance test [OGTT]) for the diagnosis of diabetes; however, its use in sub-Saharan African populations is not established. We assessed prevalence estimates and the diagnosis and detection of diabetes based on OGTT, FPG, and HbA1c in an urban black South African population.

Reducing the expression of the Indy (I'm Not Dead Yet) gene in lower organisms extends life span by mechanisms resembling caloric restriction. Similarly, deletion of the mammalian homolog, mIndy (Slc13a5), encoding for a plasma membrane tricarboxylate transporter, protects from aging- and diet-induced adiposity and insulin resistance in mice. The organ specific contribution to this phenotype is unknown. We examined the impact of selective inducible hepatic knockdown of mIndy on whole body lipid and glucose metabolism using 2'-O-methoxyethyl chimeric anti-sense oligonucleotides (ASOs) in high-fat fed rats. 4-week treatment with 2'-O-methoxyethyl chimeric ASO reduced mIndy mRNA expression by 91% (P=0.001) compared to control ASO. Besides similar body weights between both groups, mIndy-ASO treatment lead to a 74% reduction in fasting plasma insulin concentrations as well as a 35% reduction in plasma triglycerides. Moreover, hepatic triglyceride content was significantly reduced by the knockdown of mIndy, likely mediating a trend to decreased basal rates of endogenous glucose production as well as an increased suppression of hepatic glucose production by 25% during a hyperinsulinemic-euglycemic clamp. Together, these data suggest that inducible liver-selective reduction of mIndy in rats is able to ameliorate hepatic steatosis and insulin resistance, conditions occurring with high calorie diets and during aging.

The genetic architecture of common traits, including the number, frequency, and effect sizes of inherited variants that contribute to individual risk, has been long debated. Genome-wide association studies have identified scores of common variants associated with type 2 diabetes, but in aggregate, these explain only a fraction of the heritability of this disease. Here, to test the hypothesis that lower-frequency variants explain much of the remainder, the GoT2D and T2D-GENES consortia performed whole-genome sequencing in 2,657 European individuals with and without diabetes, and exome sequencing in 12,940 individuals from five ancestry groups. To increase statistical power, we expanded the sample size via genotyping and imputation in a further 111,548 subjects. Variants associated with type 2 diabetes after sequencing were overwhelmingly common and most fell within regions previously identified by genome-wide association studies. Comprehensive enumeration of sequence variation is necessary to identify functional alleles that provide important clues to disease pathophysiology, but large-scale sequencing does not support the idea that lower-frequency variants have a major role in predisposition to type 2 diabetes.

Type 2 diabetes (T2D) results from the combined effects of genetic and environmental factors on multiple tissues over time. Of the >100 variants associated with T2D and related traits in genome-wide association studies (GWAS), >90% occur in non-coding regions, suggesting a strong regulatory component to T2D risk. Here to understand how T2D status, metabolic traits and genetic variation influence gene expression, we analyse skeletal muscle biopsies from 271 well-phenotyped Finnish participants with glucose tolerance ranging from normal to newly diagnosed T2D. We perform high-depth strand-specific mRNA-sequencing and dense genotyping. Computational integration of these data with epigenome data, including ATAC-seq on skeletal muscle, and transcriptome data across diverse tissues reveals that the tissue-specific genetic regulatory architecture of skeletal muscle is highly enriched in muscle stretch/super enhancers, including some that overlap T2D GWAS variants. In one such example, T2D risk alleles residing in a muscle stretch/super enhancer are linked to increased expression and alternative splicing of muscle-specific isoforms of ANK1.

Directed differentiation of stem cells offers a scalable solution to the need for human cell models recapitulating islet biology and T2D pathogenesis. We profiled mRNA expression at 6 stages of an induced pluripotent stem cell (iPSC) model of endocrine pancreas development from 2 donors, and characterized the distinct transcriptomic profiles associated with each stage. Established regulators of endodermal lineage commitment, such as SOX17 (log2 fold change [FC] compared to iPSCs = 14.2, p-value = 4.9 × 10(-5)) and the pancreatic agenesis gene GATA6 (log2 FC = 12.1, p-value = 8.6 × 10(-5)), showed transcriptional variation consistent with their known developmental roles. However, these analyses highlighted many other genes with stage-specific expression patterns, some of which may be novel drivers or markers of islet development. For example, the leptin receptor gene, LEPR, was most highly expressed in published data from in vivo-matured cells compared to our endocrine pancreas-like cells (log2 FC = 5.5, p-value = 2.0 × 10(-12)), suggesting a role for the leptin pathway in the maturation process. Endocrine pancreas-like cells showed significant stage-selective expression of adult islet genes, including INS, ABCC8, and GLP1R, and enrichment of relevant GO-terms (e.g. "insulin secretion"; odds ratio = 4.2, p-value = 1.9 × 10(-3)): however, principal component analysis indicated that in vitro-differentiated cells were more immature than adult islets. Integration of the stage-specific expression information with genetic data from T2D genome-wide association studies revealed that 46 of 82 T2D-associated loci harbor genes present in at least one developmental stage, facilitating refinement of potential effector transcripts. Together, these data show that expression profiling in an iPSC islet development model can further understanding of islet biology and T2D pathogenesis.

Adolescence is a sensitive period for weight gain and risky health behaviors, such as smoking. Genome-wide association studies (GWAS) have identified loci contributing to adult body mass index (BMI). Evidence suggests that many of these loci have a larger influence on adolescent BMI. However, few studies have examined interactions between smoking and obesity susceptibility loci on BMI. This study investigates the interaction of current smoking and established BMI SNPs on adolescent BMI. Using data from the National Longitudinal Study of Adolescent to Adult Health, a nationally-representative, prospective cohort of the US school-based population in grades 7 to 12 (12-20 years of age) in 1994-95 who have been followed into adulthood (Wave II 1996; ages 12-21, Wave III; ages 18-27), we assessed (in 2014) interactions of 40 BMI-related SNPs and smoking status with percent of the CDC/NCHS 2000 median BMI (%MBMI) in European Americans (n = 5075), African Americans (n = 1744) and Hispanic Americans (n = 1294).

Hundreds of genes reside in structurally complex, poorly understood regions of the human genome. One such region contains the three amylase genes (AMY2B, AMY2A and AMY1) responsible for digesting starch into sugar. Copy number of AMY1 is reported to be the largest genomic influence on obesity, although genome-wide association studies for obesity have found this locus unremarkable. Using whole-genome sequence analysis, droplet digital PCR and genome mapping, we identified eight common structural haplotypes of the amylase locus that suggest its mutational history. We found that the AMY1 copy number in an individual's genome is generally even (rather than odd) and partially correlates with nearby SNPs, which do not associate with body mass index (BMI). We measured amylase gene copy number in 1,000 obese or lean Estonians and in 2 other cohorts totaling ∼3,500 individuals. We had 99% power to detect the lower bound of the reported effects on BMI, yet found no association.

Adrenal steroid hormones, which regulate a plethora of physiological functions, are produced via tightly controlled pathways. Investigations of these pathways, based on experimental data, can be facilitated by computational modeling for calculations of metabolic rate alterations. We therefore used a model system, based on mass balance and mass reaction equations, to kinetically evaluate adrenal steroidogenesis in human adrenal cortex-derived NCI H295R cells. For this purpose a panel of 10 steroids was measured by liquid chromatographic-tandem mass spectrometry. Time-dependent changes in cell incubate concentrations of steroids - including cortisol, aldosterone, dehydroepiandrosterone and their precursors - were measured after incubation with angiotensin II, forskolin and abiraterone. Model parameters were estimated based on experimental data using weighted least square fitting. Time-dependent angiotensin II- and forskolin-induced changes were observed for incubate concentrations of precursor steroids with peaks that preceded maximal increases in aldosterone and cortisol. Inhibition of 17-alpha-hydroxylase/17,20-lyase with abiraterone resulted in increases in upstream precursor steroids and decreases in downstream products. Derived model parameters, including rate constants of enzymatic processes, appropriately quantified observed and expected changes in metabolic pathways at multiple conversion steps. Our data demonstrate limitations of single time point measurements and the importance of assessing pathway dynamics in studies of adrenal cortical cell line steroidogenesis. Our analysis provides a framework for evaluation of steroidogenesis in adrenal cortical cell culture systems and demonstrates that computational modeling-derived estimates of kinetic parameters are an effective tool for describing perturbations in associated metabolic pathways.

DNA methylation has a great potential for understanding the aetiology of common complex traits such as Type 2 diabetes (T2D). Here we perform genome-wide methylated DNA immunoprecipitation sequencing (MeDIP-seq) in whole-blood-derived DNA from 27 monozygotic twin pairs and follow up results with replication and integrated omics analyses. We identify predominately hypermethylated T2D-related differentially methylated regions (DMRs) and replicate the top signals in 42 unrelated T2D cases and 221 controls. The strongest signal is in the promoter of the MALT1 gene, involved in insulin and glycaemic pathways, and related to taurocholate levels in blood. Integrating the DNA methylome findings with T2D GWAS meta-analysis results reveals a strong enrichment for DMRs in T2D-susceptibility loci. We also detect signals specific to T2D-discordant twins in the GPR61 and PRKCB genes. These replicated T2D associations reflect both likely causal and consequential pathways of the disease. The analysis indicates how an integrated genomics and epigenomics approach, utilizing an MZ twin design, can provide pathogenic insights as well as potential drug targets and biomarkers for T2D and other complex traits.

Genome wide association studies (GWAS) for fasting glucose (FG) and insulin (FI) have identified common variant signals which explain 4.8% and 1.2% of trait variance, respectively. It is hypothesized that low-frequency and rare variants could contribute substantially to unexplained genetic variance. To test this, we analyzed exome-array data from up to 33,231 non-diabetic individuals of European ancestry. We found exome-wide significant (P<5×10-7) evidence for two loci not previously highlighted by common variant GWAS: GLP1R (p.Ala316Thr, minor allele frequency (MAF)=1.5%) influencing FG levels, and URB2 (p.Glu594Val, MAF = 0.1%) influencing FI levels. Coding variant associations can highlight potential effector genes at (non-coding) GWAS signals. At the G6PC2/ABCB11 locus, we identified multiple coding variants in G6PC2 (p.Val219Leu, p.His177Tyr, and p.Tyr207Ser) influencing FG levels, conditionally independent of each other and the non-coding GWAS signal. In vitro assays demonstrate that these associated coding alleles result in reduced protein abundance via proteasomal degradation, establishing G6PC2 as an effector gene at this locus. Reconciliation of single-variant associations and functional effects was only possible when haplotype phase was considered. In contrast to earlier reports suggesting that, paradoxically, glucose-raising alleles at this locus are protective against type 2 diabetes (T2D), the p.Val219Leu G6PC2 variant displayed a modest but directionally consistent association with T2D risk. Coding variant associations for glycemic traits in GWAS signals highlight PCSK1, RREB1, and ZHX3 as likely effector transcripts. These coding variant association signals do not have a major impact on the trait variance explained, but they do provide valuable biological insights.

Inflammation in the central nervous system (CNS) and disruption of its immune privilege are major contributors to the pathogenesis of multiple sclerosis (MS) and of its rodent counterpart, experimental autoimmune encephalomyelitis (EAE). We have previously identified developmental endothelial locus-1 (Del-1) as an endogenous anti-inflammatory factor, which inhibits integrin-dependent leukocyte adhesion. Here we show that Del-1 contributes to the immune privilege status of the CNS. Intriguingly, Del-1 expression decreased in chronic-active MS lesions and in the inflamed CNS in the course of EAE. Del-1-deficiency was associated with increased EAE severity, accompanied by increased demyelination and axonal loss. As compared with control mice, Del-1(-/-) mice displayed enhanced disruption of the blood-brain barrier and increased infiltration of neutrophil granulocytes in the spinal cord in the course of EAE, accompanied by elevated levels of inflammatory cytokines, including interleukin-17 (IL-17). The augmented levels of IL-17 in Del-1-deficiency derived predominantly from infiltrated CD8(+) T cells. Increased EAE severity and neutrophil infiltration because of Del-1-deficiency was reversed in mice lacking both Del-1 and IL-17 receptor, indicating a crucial role for the IL-17/neutrophil inflammatory axis in EAE pathogenesis in Del-1(-/-) mice. Strikingly, systemic administration of Del-1-Fc ameliorated clinical relapse in relapsing-remitting EAE. Therefore, Del-1 is an endogenous homeostatic factor in the CNS protecting from neuroinflammation and demyelination. Our findings provide mechanistic underpinnings for the previous implication of Del-1 as a candidate MS susceptibility gene and suggest that Del-1-centered therapeutic approaches may be beneficial in neuroinflammatory and demyelinating disorders.

Background: Neuropathic pain is an increasingly prevalent condition and has a major impact on health and quality of life. However, the risk factors for the development and maintenance of neuropathic pain are poorly understood. Clinical, genetic and psychosocial factors all contribute to chronic pain, but their interactions have not been studied in large cohorts. The DOLORisk study aims to study these factors. Protocol: Multicentre cross-sectional and longitudinal cohorts covering the main causes leading to neuropathic pain (e.g. diabetes, surgery, chemotherapy, traumatic injury), as well as rare conditions, follow a common protocol for phenotyping of the participants. This core protocol correlates answers given by the participants on a set of questionnaires with the results of their genetic analyses. A smaller number of participants undergo deeper phenotyping procedures, including neurological examination, nerve conduction studies, threshold tracking, quantitative sensory testing, conditioned pain modulation and electroencephalography. Ethics and dissemination: All studies have been approved by their regional ethics committees as required by national law. Results are disseminated through the DOLORisk website, scientific meetings, open-access publications, and in partnership with patient organisations. Strengths and limitations: Large cohorts covering many possible triggers for neuropathic painMulti-disciplinary approach to study the interaction of clinical, psychosocial and genetic risk factorsHigh comparability of the data across centres thanks to harmonised protocolsOne limitation is that the length of the questionnaires might reduce the response rate and quality of responses of participants.

Non-secretor status due to homozygosity for the common FUT2 variant c.461G>A (p.Trp154∗) is associated with either risk for autoimmune diseases or protection against viral diarrhea and HIV. We determined the role of FUT2 in otitis media susceptibility by obtaining DNA samples from 609 multi-ethnic families and simplex case subjects with otitis media. Exome and Sanger sequencing, linkage analysis, and Fisher exact and transmission disequilibrium tests (TDT) were performed. The common FUT2 c.604C>T (p.Arg202∗) variant co-segregates with otitis media in a Filipino pedigree (LOD = 4.0). Additionally, a rare variant, c.412C>T (p.Arg138Cys), is associated with recurrent/chronic otitis media in European-American children (p = 1.2 × 10-5) and US trios (TDT p = 0.01). The c.461G>A (p.Trp154∗) variant was also over-transmitted in US trios (TDT p = 0.01) and was associated with shifts in middle ear microbiota composition (PERMANOVA p < 10-7) and increased biodiversity. When all missense and nonsense variants identified in multi-ethnic US trios with CADD > 20 were combined, FUT2 variants were over-transmitted in trios (TDT p = 0.001). Fut2 is transiently upregulated in mouse middle ear after inoculation with non-typeable Haemophilus influenzae. Four FUT2 variants-namely p.Ala104Val, p.Arg138Cys, p.Trp154∗, and p.Arg202∗-reduced A antigen in mutant-transfected COS-7 cells, while the nonsense variants also reduced FUT2 protein levels. Common and rare FUT2 variants confer susceptibility to otitis media, likely by modifying the middle ear microbiome through regulation of A antigen levels in epithelial cells. Our families demonstrate marked intra-familial genetic heterogeneity, suggesting that multiple combinations of common and rare variants plus environmental factors influence the individual otitis media phenotype as a complex trait.

Many genetic loci affect circulating lipid levels, but it remains unknown whether lifestyle factors, such as physical activity, modify these genetic effects. To identify lipid loci interacting with physical activity, we performed genome-wide analyses of circulating HDL cholesterol, LDL cholesterol, and triglyceride levels in up to 120,979 individuals of European, African, Asian, Hispanic, and Brazilian ancestry, with follow-up of suggestive associations in an additional 131,012 individuals. We find four loci, in/near CLASP1, LHX1, SNTA1, and CNTNAP2, that are associated with circulating lipid levels through interaction with physical activity; higher levels of physical activity enhance the HDL cholesterol-increasing effects of the CLASP1, LHX1, and SNTA1 loci and attenuate the LDL cholesterol-increasing effect of the CNTNAP2 locus. The CLASP1, LHX1, and SNTA1 regions harbor genes linked to muscle function and lipid metabolism. Our results elucidate the role of physical activity interactions in the genetic contribution to blood lipid levels.

The concepts of allostatic load and overload, i. e., a dramatic increase in the allostatic load that predisposes to disease, have been extensively described in the literature. Here, we show that rats engaging in active offensive response (AOR) behavioral strategies to chronic predator scent stress (PSS) display less anxiety behavior and lower plasma cortisol levels vs. rats engaging in passive defensive response (PDR) behavioral strategies to chronic PSS. In the same chronic PSS paradigm, AOR rats also have higher lactate and lower glutamate levels in amygdala but not in control-region hippocampus vs. PDR rats. The implications of these findings for regulation of allostatic and stress responses, and post-traumatic stress disorder (PTSD) are discussed.

To explore novel genetic loci for diabetic nephropathy, we performed genome-wide association studies (GWAS) for diabetic nephropathy in Japanese patients with type 2 diabetes. We analyzed the association of 5,768,242 single nucleotide polymorphisms (SNPs) in Japanese patients with type 2 diabetes, 2,380 nephropathy cases and 5,234 controls. We further performed GWAS for diabetic nephropathy using independent Japanese patients with type 2 diabetes, 429 cases and 358 controls and the results of these two GWAS were combined with an inverse variance meta-analysis (stage-1), followed by a de novo genotyping for the candidate SNP loci (p < 1.0 × 10(-4)) in an independent case-control study (Stage-2; 1,213 cases and 1,298 controls). After integrating stage-1 and stage-2 data, we identified one SNP locus, significantly associated with diabetic nephropathy; rs56094641 in FTO, P = 7.74 × 10(-10). We further examined the association of rs56094641 with diabetic nephropathy in independent Japanese patients with type 2 diabetes (902 cases and 1,221 controls), and found that the association of this locus with diabetic nephropathy remained significant after integrating all association data (P = 7.62 × 10(-10)). We have identified FTO locus as a novel locus for conferring susceptibility to diabetic nephropathy in Japanese patients with type 2 diabetes.

To investigate the genetic basis of type 2 diabetes (T2D) to high resolution, the GoT2D and T2D-GENES consortia catalogued variation from whole-genome sequencing of 2,657 European individuals and exome sequencing of 12,940 individuals of multiple ancestries. Over 27M SNPs, indels, and structural variants were identified, including 99% of low-frequency (minor allele frequency [MAF] 0.1-5%) non-coding variants in the whole-genome sequenced individuals and 99.7% of low-frequency coding variants in the whole-exome sequenced individuals. Each variant was tested for association with T2D in the sequenced individuals, and, to increase power, most were tested in larger numbers of individuals (>80% of low-frequency coding variants in ~82 K Europeans via the exome chip, and ~90% of low-frequency non-coding variants in ~44 K Europeans via genotype imputation). The variants, genotypes, and association statistics from these analyses provide the largest reference to date of human genetic information relevant to T2D, for use in activities such as T2D-focused genotype imputation, functional characterization of variants or genes, and other novel analyses to detect associations between sequence variation and T2D.

Human islet transplantation has been hampered by donor cell death associated with the islet preparation procedure before transplantation. Regulated necrosis pathways are biochemically and morphologically distinct from apoptosis. Recently, ferroptosis was identified as a non-apoptotic form of iron-dependent regulated necrosis implicated in various pathological conditions. Mediators of islet oxidative stress, including glutathione peroxidase-4 (GPX4), have been identified as inhibitors of ferroptosis, and mechanisms that affect GPX4 function can impact islet function and viability. Ferroptosis has not been investigated directly in human islets, and its relevance in islet transplantation remains unknown. Herein, we sought to determine whether in vitro human islet viability and function is compromised in the presence of two distinct ferroptosis-inducing agents (FIA), erastin or RSL3, and whether these effects could be rescued with ferroptosis inhibitors, ferrostatin-1 (Fer-1), or desferrioxamine (DFO). Viability, as assessed by lactate dehydrogenase (LDH) release, revealed significant death in erastin- and RSL3-treated islets, 20.3% ± 3.8 and 24.4% ± 2.5, 24 h post culture, respectively. These effects were ameliorated in islets pre-treated with Fer-1 or the iron chelator, desferrioxamine (DFO). Stimulation index, a marker of islet function revealed a significant reduction in function in erastin-treated islets (control 1.97 ± 0.13 vs. 50 μM erastin 1.32 ± 0.1) (p < 0.05). Fer-1 and DFO pre-treatment alone did not augment islet viability or function. Pre-treatment of islets with erastin or Fer-1 did not impact in vivo engraftment in an immunodeficient mouse transplant model. Our data reveal that islets are indeed susceptible to ferroptosis in vitro, and induction of this novel cell death modality leads to compromised islet function, which can be recoverable in the presence of the ferroptosis inhibitors. The in vivo impact of this pathway in islet transplantation remains elusive given the constraints of our study, but warrants continued investigation.