We performed genome-wide meta-analysis of lipid traits on three samples of Mexican and Mexican American ancestry comprising 4,383 individuals, and followed up significant and highly suggestive associations in three additional Hispanic samples comprising 7,876 individuals. Genome-wide significant signals were observed in or near CELSR2, ZNF259/APOA5, KANK2/DOCK6 and NCAN/MAU2 for total cholesterol, LPL, ABCA1, ZNF259/APOA5, LIPC and CETP for HDL cholesterol, CELSR2, APOB and NCAN/MAU2 for LDL cholesterol, and GCKR, TRIB1, ZNF259/APOA5 and NCAN/MAU2 for triglycerides. Linkage disequilibrium and conditional analyses indicate that signals observed at ABCA1 and LIPC for HDL cholesterol and NCAN/MAU2 for triglycerides are independent of previously reported lead SNP associations. Analyses of lead SNPs from the European Global Lipids Genetics Consortium (GLGC) dataset in our Hispanic samples show remarkable concordance of direction of effects as well as strong correlation in effect sizes. A meta-analysis of the European GLGC and our Hispanic datasets identified five novel regions reaching genome-wide significance: two for total cholesterol (FN1 and SAMM50), two for HDL cholesterol (LOC100996634 and COPB1) and one for LDL cholesterol (LINC00324/CTC1/PFAS). The top meta-analysis signals were found to be enriched for SNPs associated with gene expression in a tissue-specific fashion, suggesting an enrichment of tissue-specific function in lipid-associated loci.

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.

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.

Polychlorinated biphenyls (PCBs) are a group of man-made environmental pollutants which accumulate in humans with adverse health effects. To date, very little effort has been devoted to the study of the metabolism of PCBs on a genome-wide level.

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.

Understanding the genetic basis of airflow obstruction and smoking behaviour is key to determining the pathophysiology of chronic obstructive pulmonary disease (COPD). We used UK Biobank data to study the genetic causes of smoking behaviour and lung health.

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.

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.

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV 1), forced vital capacity (FVC) and the ratio of FEV 1 to FVC (FEV 1/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P<2·8x10 -7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs ( SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.

Transcriptome imputation has become a popular method for integrating genotype data with publicly available expression data to investigate the potentially causal role of genes in complex traits. Here, we compare three approaches (PrediXcan, MetaXcan and FUSION) via application to genome-wide association study (GWAS) data for Crohn's disease and type 1 diabetes from the Wellcome Trust Case Control Consortium. We investigate: (i) how the results of each approach compare with each other and with those of standard GWAS analysis; and (ii) how variants in the models used by the prediction tools compare with variants previously reported as eQTLs. We find that all approaches produce highly correlated results when applied to the same GWAS data, although for a subset of genes, mostly in the major histocompatibility complex, the approaches strongly disagree. We also observe that most associations detected by these methods occur near known GWAS risk loci. PrediXcan and MetaXcan's models for predicting expression more consistently recapitulate known effects of genotype on expression, suggesting they are more robust than FUSION. Application of these transcriptome imputation approaches to summary statistics from meta-analyses in Crohn's disease and type 1 diabetes detects 53 significant expression-Crohn's disease associations and 154 significant expression-type 1 diabetes associations, providing insight into biology underlying these diseases. We conclude that while current implementations of transcriptome imputation typically detect fewer associations than GWAS, they nonetheless provide an interesting way of interpreting association signals to identify potentially causal genes, and that PrediXcan and MetaXcan generally produce more reliable results than FUSION.