Tianxiao Huan

Epigenetic Signatures of Cigarette Smoking

Background—DNA methylation leaves a long-term signature of smoking exposure and is one potential mechanism by which tobacco exposure predisposes to adverse health outcomes, such as cancers, osteoporosis, lung, and cardiovascular disorders. Methods and Results—To comprehensively determine the association between cigarette smoking and DNA methylation, we conducted a meta-analysis of genome-wide DNA methylation assessed using the Illumina BeadChip 450K array on 15 907 blood-derived DNA samples from participants in 16 cohorts (including 2433 current, 6518 former, and 6956 never smokers). Comparing current versus never smokers, 2623 cytosine–phosphate–guanine sites (CpGs), annotated to 1405 genes, were statistically significantly differentially methylated at Bonferroni threshold of P<1×10−7 (18 760 CpGs at false discovery rate <0.05). Genes annotated to these CpGs were enriched for associations with several smoking-related traits in genome-wide studies including pulmonary function, cancers, inflammatory diseases, and heart disease. Comparing former versus never smokers, 185 of the CpGs that differed between current and never smokers were significant P<1×10−7 (2623 CpGs at false discovery rate <0.05), indicating a pattern of persistent altered methylation, with attenuation, after smoking cessation. Transcriptomic integration identified effects on gene expression at many differentially methylated CpGs. Conclusions—Cigarette smoking has a broad impact on genome-wide methylation that, at many loci, persists many years after smoking cessation. Many of the differentially methylated genes were novel genes with respect to biological effects of smoking and might represent therapeutic targets for prevention or treatment of tobacco-related diseases. Methylation at these sites could also serve as sensitive and stable biomarkers of lifetime exposure to tobacco smoke.

A Systems Biology Framework Identifies Molecular Underpinnings of Coronary Heart Disease

Objective—Genetic approaches have identified numerous loci associated with coronary heart disease (CHD). The molecular mechanisms underlying CHD gene–disease associations, however, remain unclear. We hypothesized that genetic variants with both strong and subtle effects drive gene subnetworks that in turn affect CHD. Approach and Results—We surveyed CHD-associated molecular interactions by constructing coexpression networks using whole blood gene expression profiles from 188 CHD cases and 188 age- and sex-matched controls. Twenty-four coexpression modules were identified, including 1 case-specific and 1 control-specific differential module (DM). The DMs were enriched for genes involved in B-cell activation, immune response, and ion transport. By integrating the DMs with gene expression–associated single-nucleotide polymorphisms and with results of genome-wide association studies of CHD and its risk factors, the control-specific DM was implicated as CHD causal based on its significant enrichment for both CHD and lipid expression–associated single-nucleotide polymorphisms. This causal DM was further integrated with tissue-specific Bayesian networks and protein–protein interaction networks to identify regulatory key driver genes. Multitissue key drivers (SPIB and TNFRSF13C) and tissue-specific key drivers (eg, EBF1) were identified. Conclusions—Our network-driven integrative analysis not only identified CHD-related genes, but also defined network structure that sheds light on the molecular interactions of genes associated with CHD risk.

Meta-analysis identifies common and rare variants influencing blood pressure and overlapping with metabolic trait loci

Meta-analyses of association results for blood pressure using exome-centric single-variant and gene-based tests identified 31 new loci in a discovery stage among 146,562 individuals, with follow-up and meta-analysis in 180,726 additional individuals (total n = 327,288). These blood pressure–associated loci are enriched for known variants for cardiometabolic traits. Associations were also observed for the aggregation of rare and low-frequency missense variants in three genes, NPR1, DBH, and PTPMT1. In addition, blood pressure associations at 39 previously reported loci were confirmed. The identified variants implicate biological pathways related to cardiometabolic traits, vascular function, and development. Several new variants are inferred to have roles in transcription or as hubs in protein–protein interaction networks. Genetic risk scores constructed from the identified variants were strongly associated with coronary disease and myocardial infarction. This large collection of blood pressure–associated loci suggests new therapeutic strategies for hypertension, emphasizing a link with cardiometabolic risk.

Lymphocyte adaptor protein LNK deficiency exacerbates hypertension and end-organ inflammation

The lymphocyte adaptor protein LNK (also known as SH2B3) is primarily expressed in hematopoietic and endothelial cells, where it functions as a negative regulator of cytokine signaling and cell proliferation. Single-nucleotide polymorphisms in the gene encoding LNK are associated with autoimmune and cardiovascular disorders; however, it is not known how LNK contributes to hypertension. Here, we determined that loss of LNK exacerbates angiotensin II-induced (Ang II-induced) hypertension and the associated renal and vascular dysfunction. At baseline, kidneys from Lnk-/- mice exhibited greater levels of inflammation, oxidative stress, and glomerular injury compared with WT animals, and these parameters were further exacerbated by Ang II infusion. Aortas from Lnk-/- mice exhibited enhanced inflammation, reduced nitric oxide levels, and impaired endothelial-dependent relaxation. Bone marrow transplantation studies demonstrated that loss of LNK in hematopoietic cells is primarily responsible for the observed renal and vascular inflammation and predisposition to hypertension. Ang II infusion increased IFN-γ-producing CD8+ T cells in the spleen and kidneys of Lnk-/- mice compared with WT mice. Moreover, IFN-γ deficiency resulted in blunted hypertension in response to Ang II infusion. Together, these results suggest that LNK is a potential therapeutic target for hypertension and its associated renal and vascular sequela.

Gene Expression Signatures of Coronary Heart Disease

Objective—To identify transcriptomic biomarkers of coronary heart disease (CHD) in 188 cases with CHD and 188 age- and sex-matched controls who were participants in the Framingham Heart Study. Approach and Results—A total of 35 genes were differentially expressed in cases with CHD versus controls at false discovery rate<0.5, including GZMB, TMEM56, and GUK1. Cluster analysis revealed 3 gene clusters associated with CHD, 2 linked to increased erythrocyte production and a third to reduced natural killer and T cell activity in cases with CHD. Exon-level results corroborated and extended the gene-level results. Alternative splicing analysis suggested that GUK1 and 38 other genes were differentially spliced in cases with CHD versus controls. Gene Ontology analysis linked ubiquitination and T-cell–related pathways with CHD. Conclusions—Two bioinformatically defined groups of genes show consistent associations with CHD. Our findings are consistent with the hypotheses that hematopoesis is upregulated in CHD, possibly reflecting a compensatory mechanism, and that innate immune activity is disrupted in CHD or altered by its treatment. Transcriptomic signatures may be useful in identifying pathways associated with CHD and point toward novel therapeutic targets for its treatment and prevention.

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.

DNA methylation signatures of chronic low-grade inflammation are associated with complex diseases

BackgroundChronic 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.ResultsWe 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 < 1.15 × 10–7) in the discovery panel of European ancestry and replicated (P < 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 < 8.47 × 10–5), ten (17%) CpG sites were associated with a nearby genetic variant (P < 2.50 × 10–3), and 51 (88%) were also associated with at least one related cardiometabolic entity (P < 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.ConclusionWe 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.

Integrative network analysis reveals molecular mechanisms of blood pressure regulation.

Genome‐wide association studies (GWAS) have identified numerous loci associated with blood pressure (BP). The molecular mechanisms underlying BP regulation, however, remain unclear. We investigated BP‐associated molecular mechanisms by integrating BP GWAS with whole blood mRNA expression profiles in 3,679 individuals, using network approaches. BP transcriptomic signatures at the single‐gene and the coexpression network module levels were identified. Four coexpression modules were identified as potentially causal based on genetic inference because expression‐related SNPs for their corresponding genes demonstrated enrichment for BP GWAS signals. Genes from the four modules were further projected onto predefined molecular interaction networks, revealing key drivers. Gene subnetworks entailing molecular interactions between key drivers and BP‐related genes were uncovered. As proof‐of‐concept, we validated SH2B3, one of the top key drivers, using Sh2b3−/− mice. We found that a significant number of genes predicted to be regulated by SH2B3 in gene networks are perturbed in Sh2b3−/− mice, which demonstrate an exaggerated pressor response to angiotensin II infusion. Our findings may help to identify novel targets for the prevention or treatment of hypertension.

Identification of common genetic variants controlling transcript isoform variation in human whole blood

An understanding of the genetic variation underlying transcript splicing is essential to dissect the molecular mechanisms of common disease. The available evidence from splicing quantitative trait locus (sQTL) studies has been limited to small samples. We performed genome-wide screening to identify SNPs that might control mRNA splicing in whole blood collected from 5,257 Framingham Heart Study participants. We identified 572,333 cis sQTLs involving 2,650 unique genes. Many sQTL-associated genes (40%) undergo alternative splicing. Using the National Human Genome Research Institute (NHGRI) genome-wide association study (GWAS) catalog, we determined that 528 unique sQTLs were significantly enriched for 8,845 SNPs associated with traits in previous GWAS. In particular, we found 395 (4.5%) GWAS SNPs with evidence of cis sQTLs but not gene-level cis expression quantitative trait loci (eQTLs), suggesting that sQTL analysis could provide additional insights into the functional mechanism underlying GWAS results. Our findings provide an informative sQTL resource for further characterizing the potential functional roles of SNPs that control transcript isoforms relevant to common diseases.

Genome-wide Identification of microRNA Expression Quantitative Trait Loci

Identification of microRNA expression quantitative trait loci (miR-eQTL) can yield insights into regulatory mechanisms of microRNA transcription, and can help elucidate the role of microRNA as mediators of complex traits. Here we present a miR-eQTL mapping study of whole blood from 5239 individuals, and identify 5269 cis-miR-eQTLs for 76 mature microRNAs. Forty-nine percent of cis-miR-eQTLs are located 300–500kb upstream of their associated intergenic microRNAs, suggesting that distal regulatory elements may affect the interindividual variability in microRNA expression levels. We find that cis-miR-eQTLs are highly enriched for cis-mRNA-eQTLs and regulatory SNPs. Among 243 cis-miR-eQTLs that were reported to be associated with complex traits in prior genome-wide association studies, many cis-miR-eQTLs miRNAs display differential expression in relation to the corresponding trait (e.g., rs7115089, miR-125b-5p, and HDL cholesterol). Our study provides a roadmap for understanding the genetic basis of miRNA expression, and sheds light on miRNA involvement in a variety of complex traits.