Integrative analysis of gene expression, DNA methylation, physiological traits, and genetic variation in human skeletal muscle

We integrate comeasured gene expression and DNA methylation (DNAme) in 265 human skeletal muscle biopsies from the FUSION study with >7 million genetic variants and eight physiological traits: height, waist, weight, waist-hip ratio, body mass index, fasting serum insulin, fasting plasma glucose, and type 2 diabetes. We find hundreds of genes and DNAme sites associated with fasting insulin, waist, and body mass index, as well as thousands of DNAme sites associated with gene expression (eQTM). We find that controlling for heterogeneity in tissue/muscle fiber type reduces the number of physiological trait associations, and that long-range eQTMs (>1 Mb) are reduced when controlling for tissue/muscle fiber type or latent factors. We map genetic regulators (quantitative trait loci; QTLs) of expression (eQTLs) and DNAme (mQTLs). Using Mendelian randomization (MR) and mediation techniques, we leverage these genetic maps to predict 213 causal relationships between expression and DNAme, approximately two-thirds of which predict methylation to causally influence expression. We use MR to integrate FUSION mQTLs, FUSION eQTLs, and GTEx eQTLs for 48 tissues with genetic associations for 534 diseases and quantitative traits. We identify hundreds of genes and thousands of DNAme sites that may drive the reported disease/quantitative trait genetic associations. We identify 300 gene expression MR associations that are present in both FUSION and GTEx skeletal muscle and that show stronger evidence of MR association in skeletal muscle than other tissues, which may partially reflect differences in power across tissues. As one example, we find that increased RXRA muscle expression may decrease lean tissue mass.Peer reviewe

An RBP4 promoter polymorphism increases risk of type 2 diabetes

Aims/hypothesis: Retinol-binding protein 4 (RBP4), originally known for retinol transport, was recently identified as an adipokine affecting insulin resistance. The RBP4 -803GA promoter polymorphism influences binding of hepatic nuclear factor 1α and is associated with type 2 diabetes in case-control studies. We hypothesised that the RBP4 -803GA polymorphism increases type 2 diabetes risk at a population-based level. In addition, information on retinol intake and plasma vitamin A levels enabled us to explore the possible underlying mechanism. Methods: In the Rotterdam Study, a prospective, population-based, follow-up study, the -803GA polymorphism was genotyped. In Cox proportional hazards models, associations of the -803GA polymorphism and retinol intake with type 2 diabetes risk were examined. Moreover, the interaction of the polymorphism with retinol intake on type 2 diabetes risk was assessed. In a subgroup of participants the association of the polymorphism and vitamin A plasma levels was investigated. Results: Homozygous carriers of the -803A allele had increased risk of type 2 diabetes (HR 1.83; 95% CI 1.26-2.66). Retinol intake was not associated with type 2 diabetes risk and showed no interaction with the RBP4 -803GA polymorphism. Furthermore, there was no significant association of the polymorphism with plasma vitamin A levels. Conclusions/interpretation: Our results provide evidence that homozygosity for the RBP4 -803A allele is associated with increased risk of type 2 diabetes in the Rotterdam population. This relationship was not clearly explained by retinol intake and vitamin A plasma levels. Therefore, we cannot differentiate between a retinol-dependent or -independent mechanism of this RBP4 variant

Interactions between genetic variation and cellular environment in skeletal muscle gene expression

From whole organisms to individual cells, responses to environmental conditions are influenced by genetic makeup, where the effect of genetic variation on a trait depends on the environmental context. RNA-sequencing quantifies gene expression as a molecular trait, and is capable of capturing both genetic and environmental effects. In this study, we explore opportunities of using allele-specific expression (ASE) to discover cis-acting genotype-environment interactions (GxE)-genetic effects on gene expression that depend on an environmental condition. Treating 17 common, clinical traits as approximations of the cellular environment of 267 skeletal muscle biopsies, we identify 10 candidate environmental response expression quantitative trait loci (reQTLs) across 6 traits (12 unique gene-environment trait pairs; 10% FDR per trait) including sex, systolic blood pressure, and low-density lipoprotein cholesterol. Although using ASE is in principle a promising approach to detect GxE effects, replication of such signals can be challenging as validation requires harmonization of environmental traits across cohorts and a sufficient sampling of heterozygotes for a transcribed SNP. Comprehensive discovery and replication will require large human transcriptome datasets, or the integration of multiple transcribed SNPs, coupled with standardized clinical phenotyping.Peer reviewe

Global microRNA expression profiles in insulin target tissues in a spontaneous rat model of type 2 diabetes

AIMS/HYPOTHESIS: MicroRNAs regulate a broad range of biological mechanisms. To investigate the relationship between microRNA expression and type 2 diabetes, we compared global microRNA expression in insulin target tissues from three inbred rat strains that differ in diabetes susceptibility. METHODS: Using microarrays, we measured the expression of 283 microRNAs in adipose, liver and muscle tissue from hyperglycaemic (Goto-Kakizaki), intermediate glycaemic (Wistar Kyoto) and normoglycaemic (Brown Norway) rats (n = 5 for each strain). Expression was compared across strains and validated using quantitative RT-PCR. Furthermore, microRNA expression variation in adipose tissue was investigated in 3T3-L1 adipocytes exposed to hyperglycaemic conditions. RESULTS: We found 29 significantly differentiated microRNAs (p(adjusted) < 0.05): nine in adipose tissue, 18 in liver and two in muscle. Of these, five microRNAs had expression patterns that correlated with the strain-specific glycaemic phenotype. MiR-222 (p(adjusted) = 0.0005) and miR-27a (p(adjusted) = 0.006) were upregulated in adipose tissue; miR-195 (p(adjusted) = 0.006) and miR-103 (p(adjusted) = 0.04) were upregulated in liver; and miR-10b (p(adjusted) = 0.004) was downregulated in muscle. Exposure of 3T3-L1 adipocytes to increased glucose concentration upregulated the expression of miR-222 (p = 0.008), miR-27a (p = 0.02) and the previously reported miR-29a (p = 0.02). Predicted target genes of these differentially expressed microRNAs are involved in pathways relevant to type 2 diabetes. CONCLUSION: The expression patterns of miR-222, miR-27a, miR-195, miR-103 and miR-10b varied with hyperglycaemia, suggesting a role for these microRNAs in the pathophysiology of type 2 diabetes, as modelled by the Gyoto-Kakizaki rat. We observed similar patterns of expression of miR-222, miR-27a and miR-29a in adipocytes as a response to increased glucose levels, which supports our hypothesis that altered expression of microRNAs accompanies primary events related to the pathogenesis of type 2 diabetes

Discovery of rare variants associated with blood pressure regulation through meta-analysis of 1.3 million individuals

Correction: Volume53, Issue5 Page 762-762 DOI: 10.1038/s41588-021-00832-z Published MAY 2021Genetic studies of blood pressure (BP) to date have mainly analyzed common variants (minor allele frequency > 0.05). In a meta-analysis of up to similar to 1.3 million participants, we discovered 106 new BP-associated genomic regions and 87 rare (minor allele frequencyPeer reviewe

Interspecific Proteomic Comparisons Reveal Ash Phloem Genes Potentially Involved in Constitutive Resistance to the Emerald Ash Borer

The emerald ash borer (Agrilus planipennis) is an invasive wood-boring beetle that has killed millions of ash trees since its accidental introduction to North America. All North American ash species (Fraxinus spp.) that emerald ash borer has encountered so far are susceptible, while an Asian species, Manchurian ash (F. mandshurica), which shares an evolutionary history with emerald ash borer, is resistant. Phylogenetic evidence places North American black ash (F. nigra) and Manchurian ash in the same clade and section, yet black ash is highly susceptible to the emerald ash borer. This contrast provides an opportunity to compare the genetic traits of the two species and identify those with a potential role in defense/resistance. We used Difference Gel Electrophoresis (DIGE) to compare the phloem proteomes of resistant Manchurian to susceptible black, green, and white ash. Differentially expressed proteins associated with the resistant Manchurian ash when compared to the susceptible ash species were identified using nano-LC-MS/MS and putative identities assigned. Proteomic differences were strongly associated with the phylogenetic relationships among the four species. Proteins identified in Manchurian ash potentially associated with its resistance to emerald ash borer include a PR-10 protein, an aspartic protease, a phenylcoumaran benzylic ether reductase (PCBER), and a thylakoid-bound ascorbate peroxidase. Discovery of resistance-related proteins in Asian species will inform approaches in which resistance genes can be introgressed into North American ash species. The generation of resistant North American ash genotypes can be used in forest ecosystem restoration and urban plantings following the wake of the emerald ash borer invasion