A meta-analysis of gene expression signatures of blood pressure and hypertension

Genome-wide association studies (GWAS) have uncovered numerous genetic variants (SNPs) that are associated with blood pressure (BP). Genetic variants may lead to BP changes by acting on intermediate molecular phenotypes such as coded protein sequence or gene expression, which in turn affect BP variability. Therefore, characterizing genes whose expression is associated with BP may reveal cellular processes involved in BP regulation and uncover how transcripts mediate genetic and environmental effects on BP variability. A meta-analysis of results from six studies of global gene expression profiles of BP and hypertension in whole blood was performed in 7017 individuals who were not receiving antihypertensive drug treatment. We identified 34 genes that were differentially expressed in relation to BP (Bonferroni-corrected p<0.05). Among these genes, FOS and PTGS2 have been previously reported to be involved in BP-related processes; the others are novel. The top BP signature genes in aggregate explain 5%-9% of inter-individual variance in BP. Of note, rs3184504 in SH2B3, which was also reported in GWAS to be associated with BP, was found to be a trans regulator of the expression of 6 of the transcripts we found to be associated with BP (FOS, MYADM, PP1R15A, TAGAP, S100A10, and FGBP2). Gene set enrichment analysis suggested that the BP-related global gene expression changes include genes involved in inflammatory response and apoptosis pathways. Our study provides new insights into molecular mechanisms underlying BP regulation, and suggests novel transcriptomic markers for the treatment and prevention of hypertension

Genome-wide association analyses identify 143 risk variants and putative regulatory mechanisms for type 2 diabetes

Type 2 diabetes (T2D) is a very common disease in humans. Here we conduct a meta-analysis of genome-wide association studies (GWAS) with ~16 million genetic variants in 62,892 T2D cases and 596,424 controls of European ancestry. We identify 139 common and 4 rare variants associated with T2D, 42 of which (39 common and 3 rare variants) are independent of the known variants. Integration of the gene expression data from blood (n = 14,115 and 2765) with the GWAS results identifies 33 putative functional genes for T2D, 3 of which were targeted by approved drugs. A further integration of DNA methylation (n = 1980) and epigenomic annotation data highlight 3 genes (CAMK1D, TP53INP1, and ATP5G1) with plausible regulatory mechanisms, whereby a genetic variant exerts an effect on T2D through epigenetic regulation of gene expression. Our study uncovers additional loci, proposes putative genetic regulatory mechanisms for T2D, and provides evidence of purifying selection for T2D-associated variants

A whole-blood transcriptome meta-analysis identifies gene expression signatures of cigarette smoking

Cigarette smoking is a leading modifiable cause of death worldwide. We hypothesized that cigarette smoking induces extensive transcriptomic changes that lead to target-organ damage and smoking-related diseases. We performed a metaanalysis of transcriptome-wide gene expression using whole blood-derived RNA from 10,233 participants of European ancestry in six cohorts (including 1421 current and 3955 former smokers) to identify associations between smoking and altered gene expression levels. At a false discovery rate (FDR) < 0.1, we identified 1270 differentially expressed genes in current vs. never smokers, and 39 genes in former vs. never smokers. Expression levels of 12 genes remained elevated up to 30 years after smoking cessation, suggesting that the molecular consequence of smoking may persist for decades. Gene ontology analysis revealed enrichment of smoking-related genes for activation of platelets and lymphocytes, immune response, and apoptosis. Many of the top smoking-related differentially expressed genes, including LRRN3 and GPR15, have DNA methylation loci in promoter regions that were recently reported to be hypomethylated among smokers. By linking differential gene expression with smoking-related disease phenotypes, we demonstrated that stroke and pulmonary function show enrichment for smoking-related gene expression signatures. Mediation analysis revealed the expression of several genes (e.g. ALAS2) to be putative mediators of the associations between smoking and inflammatory biomarkers (IL6 and C-re

Cell Specific eQTL Analysis without Sorting Cells

The functional consequences of trait associated SNPs are often investigated using expression quantitative trait locus (eQTL) mapping. While trait-associated variants may operate in a cell-type specific manner, eQTL datasets for such cell-types may not always be available. We performed a genome-environment interaction (GxE) meta-analysis on data from 5,683 samples to infer the cell type specificity of whole blood cis-eQTLs. We demonstrate that this method is able to predict neutrophil and lymphocyte specific cis-eQTLs and replicate these predictions in independent cell-type specific datasets. Finally, we show that SNPs associated with Crohn’s disease preferentially affect gene expression within neutrophils, including the archetypal NOD2 locus

Persistent Pesticides in Human Breast Milk and Cryptorchidism

INTRODUCTION: Prenatal exposure to some pesticides can adversely affect male reproductive health in animals. We investigated a possible human association between maternal exposure to 27 organochlorine compounds used as pesticides and cryptorchidism among male children. DESIGN: Within a prospective birth cohort, we performed a case–control study; 62 milk samples from mothers of cryptorchid boys and 68 from mothers of healthy boys were selected. Milk was collected as individual pools between 1 and 3 months postpartum and analyzed for 27 organochlorine pesticides. RESULTS: Eight organochlorine pesticides were measurable in all samples (medians; nanograms per gram lipid) for cases/controls: 1,1-dichloro-2,2-bis(4-chlorophenyl)ethylene (p,p′-DDE): 97.3/83.8; β-hexachlorocyclohexane (β-HCH): 13.6/12.3; hexachlorobenzene (HCB): 10.6/8.8; α -endosulfan: 7.0/6.7; oxychlordane: 4.5/4.1; 1,1,1-trichloro-2,2-bis(4-chlorophenyl)ethane (p,p′-DDT): 4.6/4.0; dieldrin: 4.1/3.1; cis-heptachloroepoxide (cis-HE): 2.5/2.2. Five compounds [octachlorostyrene (OCS); pentachlorobenzene, 1,1-dichloro-2,2-bis(4-chlorophenyl)ethane (p,p′-DDD); o,p′-DDT; mirex] were measurable in most samples (detection rates 90.8–99.2%) but in lower concentrations. For methoxychlor, cis-chlordane, pentachloroanisole (PCA), γ -HCH, 1,1-dichloro-2-(2-chlorophenyl)-2,2(4-chlorophenyl)ethane, trans-chlordane, α -HCH, and o,p′-DDE, both concentrations and detection rates were low (26.5–71.5%). Heptachlor, HCH (δ, ɛ ), aldrin, β-endosulfan and trans-heptachloroepoxide were detected at negligible concentrations and low detection rates and were not analyzed further. Seventeen of 21 organochlorine pesticides [p,p′-DDT, p,p′-DDE, p,p′-DDD, o,p′-DDT, HCH (α , β, γ ), HCB, PCA, α -endosulfan, cis-HE, chlordane (cis-, trans-) oxychlordane, methoxychlor, OCS, and dieldrin] were measured in higher median concentrations in case milk than in control milk. Apart from trans-chlordane (p = 0.012), there were no significant differences between cryptorchid and healthy boys for individual chemicals. However, combined statistical analysis of the eight most abundant persistent pesticides showed that pesticide levels in breast milk were significantly higher in boys with cryptorchidism (p = 0.032). CONCLUSION: The association between congenital cryptorchidism and some persistent pesticides in breast milk as a proxy for maternal exposure suggests that testicular descent in the fetus may be adversely affected

DNA Methylation Signatures of Chronic Low-Grade Inflammation Are Associated with Complex Diseases

Background: Chronic low-grade inflammation reflects a subclinical immune response implicated in the pathogenesis of complex diseases. Identifying genetic loci where DNA methylation is associated with chronic low-grade inflammation may reveal novel pathways or therapeutic targets for inflammation. Results: We performed a meta-analysis of epigenome-wide association studies (EWAS) of serum C-reactive protein (CRP), which is a sensitive marker of low-grade inflammation, in a large European population (n = 8863) and trans-ethnic replication in African Americans (n = 4111). We found differential methylation at 218 CpG sites to be associated with CRP (P \u3c 1.15 × 10–7) in the discovery panel of European ancestry and replicated (P \u3c 2.29 × 10–4) 58 CpG sites (45 unique loci) among African Americans. To further characterize the molecular and clinical relevance of the findings, we examined the association with gene expression, genetic sequence variants, and clinical outcomes. DNA methylation at nine (16%) CpG sites was associated with whole blood gene expression in cis (P \u3c 8.47 × 10–5), ten (17%) CpG sites were associated with a nearby genetic variant (P \u3c 2.50 × 10–3), and 51 (88%) were also associated with at least one related cardiometabolic entity (P \u3c 9.58 × 10–5). An additive weighted score of replicated CpG sites accounted for up to 6% inter-individual variation (R2) of age-adjusted and sex-adjusted CRP, independent of known CRP-related genetic variants. Conclusion: We have completed an EWAS of chronic low-grade inflammation and identified many novel genetic loci underlying inflammation that may serve as targets for the development of novel therapeutic interventions for inflammation

Genome-Wide Interaction Analyses of Serum Calcium on Ventricular Repolarization Time in 125 393 Participants

Background: Ventricular repolarization time (ECG QT and JT intervals) is associated with malignant arrhythmia. Genome-wide association studies have identified 230 independent loci for QT and JT; however, 50% of their heritability remains unexplained. Previous work supports a causal effect of lower serum calcium concentrations on longer ventricular repolarization time. We hypothesized calcium interactions with QT and JT variant associations could explain a proportion of the missing heritability. Methods and Results: We performed genome-wide calcium interaction analyses for QT and JT intervals. Participants were stratified by their calcium level relative to the study distribution (top or bottom 20%). We performed a 2-stage analysis (genome-wide discovery [N=62 532] and replication [N=59 861] of lead variants) and a single-stage genome-wide meta-analysis (N=122 393, [European ancestry N=117 581, African ancestry N=4812]). We also calculated 2-degrees of freedom joint main and interaction and 1-degree of freedom interaction P values. In 2-stage and single-stage analyses, 50 and 98 independent loci, respectively, were associated with either QT or JT intervals (2-degrees of freedom joint main and interaction P value &lt;5x10(-8)). No lead variant had a significant interaction result after correcting for multiple testing and sensitivity analyses provided similar findings. Two loci in the single-stage meta-analysis were not reported previously (SPPL2B and RFX6). Conclusions: We have found limited support for an interaction effect of serum calcium on QT and JT variant associations despite sample sizes with suitable power to detect relevant effects. Therefore, such effects are unlikely to explain a meaningful proportion of the heritability of QT and JT, and factors including rare variation and other environmental interactions need to be considered

GWAS on longitudinal growth traits reveals different genetic factors influencing infant, child, and adult BMI

Early childhood growth patterns are associated with adult health, yet the genetic factors and the developmental stages involved are not fully understood. Here, we combine genome-wide association studies with modeling of longitudinal growth traits to study the genetics of infant and child growth, followed by functional, pathway, genetic correlation, risk score, and colocalization analyses to determine how developmental timings, molecular pathways, and genetic determinants of these traits overlap with those of adult health. We found a robust overlap between the genetics of child and adult body mass index (BMI), with variants associated with adult BMI acting as early as 4 to 6 years old. However, we demonstrated a completely distinct genetic makeup for peak BMI during infancy, influenced by variation at the LEPR/LEPROT locus. These findings suggest that different genetic factors control infant and child BMI. In light of the obesity epidemic, these findings are important to inform the timing and targets of prevention strategies.</p

GWAS on longitudinal growth traits reveals different genetic factors influencing infant, child, and adult BMI

Early childhood growth patterns are associated with adult health, yet the genetic factors and the developmental stages involved are not fully understood. Here, we combine genome-wide association studies with modeling of longitudinal growth traits to study the genetics of infant and child growth, followed by functional, pathway, genetic correlation, risk score, and colocalization analyses to determine how developmental timings, molecular pathways, and genetic determinants of these traits overlap with those of adult health. We found a robust overlap between the genetics of child and adult body mass index (BMI), with variants associated with adult BMI acting as early as 4 to 6 years old. However, we demonstrated a completely distinct genetic makeup for peak BMI during infancy, influenced by variation at the LEPR/LEPROT locus. These findings suggest that different genetic factors control infant and child BMI. In light of the obesity epidemic, these findings are important to inform the timing and targets of prevention strategies

Large-scale cis- and trans-eQTL analyses identify thousands of genetic loci and polygenic scores that regulate blood gene expression

AbstractTrait-associated genetic variants affect complex phenotypes primarily via regulatory mechanisms on the transcriptome. To investigate the genetics of gene expression, we performed cis- and trans-expression quantitative trait locus (eQTL) analyses using blood-derived expression from 31,684 individuals through the eQTLGen Consortium. We detected cis-eQTL for 88% of genes, and these were replicable in numerous tissues. Distal trans-eQTL (detected for 37% of 10,317 trait-associated variants tested) showed lower replication rates, partially due to low replication power and confounding by cell type composition. However, replication analyses in single-cell RNA-seq data prioritized intracellular trans-eQTL. Trans-eQTL exerted their effects via several mechanisms, primarily through regulation by transcription factors. Expression of 13% of the genes correlated with polygenic scores for 1,263 phenotypes, pinpointing potential drivers for those traits. In summary, this work represents a large eQTL resource, and its results serve as a starting point for in-depth interpretation of complex phenotypes.Abstract Trait-associated genetic variants affect complex phenotypes primarily via regulatory mechanisms on the transcriptome. To investigate the genetics of gene expression, we performed cis- and trans-expression quantitative trait locus (eQTL) analyses using blood-derived expression from 31,684 individuals through the eQTLGen Consortium. We detected cis-eQTL for 88% of genes, and these were replicable in numerous tissues. Distal trans-eQTL (detected for 37% of 10,317 trait-associated variants tested) showed lower replication rates, partially due to low replication power and confounding by cell type composition. However, replication analyses in single-cell RNA-seq data prioritized intracellular trans-eQTL. Trans-eQTL exerted their effects via several mechanisms, primarily through regulation by transcription factors. Expression of 13% of the genes correlated with polygenic scores for 1,263 phenotypes, pinpointing potential drivers for those traits. In summary, this work represents a large eQTL resource, and its results serve as a starting point for in-depth interpretation of complex phenotypes