Common variation neighbouring micro-RNA 22 is associated with increased left ventricular mass

Single-cell reconstruction and mutation enrichment analysis identifies dysregulated cardiomyocyte and endothelial cells in congenital heart disease

Congenital heart disease (CHD) is one of the most prevalent neonatal congenital anomalies. To catalog the putative candidate CHD risk genes, we collected 16,349 variants [single-nucleotide variants (SNVs) and Indels] impacting 8,308 genes in 3,166 CHD cases for a comprehensive meta-analysis. Using American College of Medical Genetics (ACMG) guidelines, we excluded the 0.1% of benign/likely benign variants and the resulting dataset consisted of 83% predicted loss of function variants and 17% missense variants. Seventeen percent were de novo variants. A stepwise analysis identified 90 variant-enriched CHD genes, of which six (GPATCH1, NYNRIN, TCLD2, CEP95, MAP3K19, and TTC36) were novel candidate CHD genes. Single-cell transcriptome cluster reconstruction analysis on six CHD tissues and four controls revealed upregulation of the top 10 frequently mutated genes primarily in cardiomyocytes. NOTCH1 (highest number of variants) and MYH6 (highest number of recurrent variants) expression was elevated in endocardial cells and cardiomyocytes, respectively, and 60% of these gene variants were associated with tetralogy of Fallot and coarctation of the aorta, respectively. Pseudobulk analysis using the single-cell transcriptome revealed significant (P < 0.05) upregulation of both NOTCH1 (endocardial cells) and MYH6 (cardiomyocytes) in the control heart data. We observed nine different subpopulations of CHD heart cardiomyocytes of which only four were observed in the control heart. This is the first comprehensive meta-analysis combining genomics and CHD single-cell transcriptomics, identifying the most frequently mutated CHD genes, and demonstrating CHD gene heterogeneity, suggesting that multiple genes contribute to the phenotypic heterogeneity of CHD. Cardiomyocytes and endocardial cells are identified as major CHD-related cell types.NEW & NOTEWORTHY Congential heart disease (CHD) is one of the most prevalent neonatal congenital anomalies. We present a comprehensive analysis combining genomics and CHD single-cell transcriptome. Our study identifies 90 potential candidate CHD risk genes of which 6 are novel. The risk genes have heterogenous expression suggestive of multiple genes contributing to the phenotypic heterogeneity of CHD. Cardiomyocytes and endocardial cells are identified as major CHD-related cell types.

Identification of key genes for hypertrophic cardiomyopathy using integrated network analysis of differential lncRNA and gene expression

Objective

Hypertrophic cardiomyopathy (HCM) is a complex heterogeneous heart disease. Recent reports found that long non-coding RNAs (lncRNAs) play an important role in the progression of cardiovascular diseases. The present study aimed to identify the novel lncRNAs and messenger RNAs (mRNAs) and determine the key pathways involved in HCM.

Methods

The lncRNA and mRNA sequencing datasets of GSE68316 and GSE130036 were downloaded from the Gene Expression Omnibus (GEO) database. An integrated co-expression network analysis was conducted to identify differentially expressed lncRNAs and differentially expressed mRNAs in patients with HCM. Then, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were explored to identify the biological functions and signaling pathways of the co-expression network. Protein–protein interaction (PPI) and hub gene networks were constructed by using Cytoscape software. Plasma samples of patients with HCM and the GSE89714 dataset were used to validate the bioinformatics results.

Results

A total of 1,426 differentially expressed long non-coding RNAs (lncRNAs) and 1,715 differentially expressed mRNAs were obtained from GSE68316, of which 965 lncRNAs and 896 mRNAs were upregulated and 461 lncRNAs and 819 mRNAs were downregulated. A total of 469 differentially expressed lncRNAs and 2,407 differentially expressed mRNAs were screened from GSE130036, of which 183 lncRNAs and 1,283 mRNAs were upregulated and 286 lncRNAs and 1,124 mRNAs were downregulated. A co-expression network was constructed and contained 30 differentially expressed lncRNAs and 63 differentially expressed mRNAs, which were primarily involved in ‘G-protein beta/gamma-subunit complex binding,' ‘polyubiquitin modification-dependent protein binding,' ‘Apelin signaling pathway,' and ‘Wnt signaling pathway.' The 10 hub genes in the upregulated network [G Protein Subunit Alpha I2 (GNAI2), G Protein Subunit Alpha I1 (GNAI1), G Protein Subunit Alpha I3 (GNAI3), G Protein Subunit Gamma 2 (GNG2), G Protein Subunit Beta 1 (GNB1), G Protein Subunit Gamma 13 (GNG13), G Protein Subunit Gamma Transducin 1 (GNGT1), G Protein Subunit Gamma 12 (GNG12), AKT Serine/Threonine Kinase 1 (AKT1) and GNAS Complex Locus (GNAS)] and the 10 hub genes in the downregulated network [Nucleotide-Binding Oligomerization Domain Containing Protein 2 (NOD2), Receptor-Interacting Serine/Threonine Kinase 2 (RIPK2), Nucleotide-Binding Oligomerization Domain Containing Protein 1 (NOD1), Mitochondrial Antiviral Signaling Protein (MAVS), Autophagy Related 16-Like 1 (ATG16L1), Interferon Induced With Helicase C Domain 1 (IFIH1), Autophagy Related 5 (ATG5), TANK-Binding Kinase 1 (TBK1), Caspase Recruitment Domain Family Member 9 (CARD9), and von Willebrand factor (VWF)] were screened using cytoHubba. The expression of LA16c-312E8.2 and RP5-1160K1.3 in the plasma of patients with HCM was elevated, and the expression of the MIR22 host gene (MIR22HG) was decreased, which was consistent with our analysis, while the expression of LINC00324 and Small Nucleolar RNA Host Gene 12 (SNHG12) was not significantly different between the two groups. Verification analyses performed on GSE89714 showed the upregulated mRNAs of Chloride Voltage-Gated Channel 7 (CLCN7), N-Acetylglucosamine-1-Phosphate Transferase Subunit Gamma (GNPTG), Unk Like Zinc Finger (UNKL), Adenosine Monophosphate Deaminase 2 (AMPD2), GNAI3, WD Repeat Domain 81 (WDR81), and Serpin Family F Member 1 (SERPINF1) and downregulated mRNAs of TATA-Box Binding Protein Associated Factor 12 (TAF12) co-expressed with five crucial lncRNAs. Moreover, GNAI2, GNAI3, GNG12, and vWF were upregulated and GNAS was downregulated in the top 10 hub genes of upregulated and downregulated PPI networks.

Conclusion

These findings from integrative biological analysis of lncRNA-mRNA co-expression networks explored the key genes and pathways and provide new insights into the understanding of the mechanism and discovering new therapeutic targets for HCM. Three differentially expressed pivotal lncRNAs (LA16c-312E8.2, RP5-1160K1.3, and MIR22HG) in the co-expression network may serve as biomarkers and intervention targets for the diagnosis and treatment of HCM.

Sequencing of Linkage Region on Chromosome 12p11 Identifies PKP2 as a Candidate Gene for Left Ventricular Mass in Dominican Families

Increased left ventricular mass (LVM) is an intermediate phenotype for cardiovascular disease (CVD) and a predictor of stroke. Using families from the Dominican Republic, we have previously shown LVM to be heritable and found evidence for linkage to chromosome 12p11. Our current study aimed to further characterize the QTL by sequencing the 1 LOD unit down region in 10 families from the Dominican Republic with evidence for linkage to LVM. Within this region, we tested 5477 common variants [CVs; minor allele frequency (MAF) ≥5%] using the Quantitative Transmission-Disequilibrium Test (QTDT). Gene-based analyses were performed to test rare variants (RVs; MAF < 5%) in 181 genes using the family-based sequence kernel association test. A sample of 618 unrelated Dominicans from the Northern Manhattan Study (NOMAS) and 12 Dominican families with Exome Array data were used for replication analyses. The most strongly associated CV with evidence for replication was rs1046116 (Discovery families P = 9.0 × 10-4; NOMAS P = 0.03; replication families P = 0.46), a missense variant in PKP2 In nonsynonymous RV analyses, PKP2 was one of the most strongly associated genes (P = 0.05) with suggestive evidence for replication in NOMAS (P = 0.05). PKP2 encodes the plakophilin 2 protein and is a desmosomal gene implicated in arrythmogenic right ventricular cardiomyopathy and recently in arrhythmogenic left ventricular cardiomyopathy, which makes PKP2 an excellent candidate gene for LVM. In conclusion, sequencing of our previously reported QTL identified common and rare variants within PKP2 to be associated with LVM. Future studies are necessary to elucidate the role these variants play in influencing LVM.

The T Allele of rs8075977 in the 5'-Flanking Region of the PEDF Gene Is Associated with Reduced Risk of Coronary Artery Disease in Elderly Chinese Men

Coronary artery disease (CAD) is a multifactorial disease with a genetic component. Pigment epithelium-derived factor (PEDF) exerts anti-inflammatory, anti-oxidant, anti-thrombotic, and anti-angiogenic effects and thus has received increasing attention as a sensitive biomarker of atherosclerosis and CAD. To explore the potential association between PEDF single nucleotide polymorphisms (SNPs) and CAD, we performed this case-control study of consecutive elderly Chinese Han male patients (n = 416) and age-matched male controls (n = 528) without a history of CAD or electrocardiographic signs of CAD. The enrolled CAD patients (age ≥ 60 years) are not biologically related. A tag approach was used to examine 100% of common variations in the PEDF gene (r² ≥ 0.8, minor allele frequency > 0.1). PEDF tag SNPs (tSNPs) were selected using the HapMap Data-CHB which describes the common patterns of human DNA sequence variation and Tagger program. SNPs were genotyped using ligase detection reaction (LDR). Seven tSNPs (rs8075977, rs11658342, rs1136287, rs12603825, rs12453107, rs6828 and rs11078634) were selected. Among them, only one SNP, rs8075977 (C/T) located in the 5'-flanking region, showed the significant effect on the susceptibility to CAD. The frequency of its T allele was significantly higher in the controls (52.7%) than that in the CAD group (46.2%) (adjusted OR = 0.88, 95% CI: 0.80-0.96; P = 0.005). In conclusion, the T allele of rs8075977 in the 5'-flanking region of the PEDF gene may be protective for CAD. Conversely, the C allele at this variation site is associated with CAD in elderly Chinese Han men.

Genetic polymorphisms associated with heart failure: A literature review

Objective

To review possible associations reported between genetic variants and the risk, therapeutic response and prognosis of heart failure.

Methods

Electronic databases (PubMed, Web of Science and CNKI) were systematically searched for relevant papers, published between January 1995 and February 2015.

Results

Eighty-two articles covering 29 genes and 39 polymorphisms were identified.

Conclusion

Genetic association studies of heart failure have been highly controversial. There may be interaction or synergism of several genetic variants that together result in the ultimate pathological phenotype for heart failure.

RNA expression profiling of human iPSC-derived cardiomyocytes in a cardiac hypertrophy model

Cardiac hypertrophy is an independent risk factor for cardiovascular disease and heart failure. There is increasing evidence that microRNAs (miRNAs) play an important role in the regulation of messenger RNA (mRNA) and the pathogenesis of various cardiovascular diseases. However, the ability to comprehensively study cardiac hypertrophy on a gene regulatory level is impacted by the limited availability of human cardiomyocytes. Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) offer the opportunity for disease modeling. Here we utilize a previously established invitro model of cardiac hypertrophy to interrogate the regulatory mechanism associated with the cardiac disease process. We perform miRNA sequencing and mRNA expression analysis on endothelin 1 (ET-1) stimulated hiPSC-CMs to describe associated RNA expression profiles. MicroRNA sequencing revealed over 250 known and 34 predicted novel miRNAs to be differentially expressed between ET-1 stimulated and unstimulated control hiPSC-CMs. Messenger RNA expression analysis identified 731 probe sets with significant differential expression. Computational target prediction on significant differentially expressed miRNAs and mRNAs identified nearly 2000 target pairs. A principal component analysis approach comparing the invitro data with human myocardial biopsies detected overlapping expression changes between the invitro samples and myocardial biopsies with Left Ventricular Hypertrophy. These results provide further insights into the complex RNA regulatory mechanism associated with cardiac hypertrophy.

miR-22 in cardiac remodeling and disease

Regulation of gene expression during cardiac development and remodeling is very complicated, involving epigenetic, transcriptional, post-transcriptional, and translational regulation. Our understanding of the molecular mechanisms underlying cardiac remodeling is still far from complete. MicroRNAs are a class of small non-coding RNAs that have been shown to play critical roles in gene regulation in cardiovascular biology and disease. microRNA-22 (miR-22) is an evolutionally conserved miRNA that is highly expressed in the heart. Recent studies uncovered miR-22 as an important regulator for cardiac remodeling. miR-22 modulates the expression and function of genes involved in hypertrophic response, sarcomere reorganization, and metabolic program shift during cardiac remodeling. In this review, we will focus on the recent findings of miR-22 in cardiac remodeling and the therapeutic potential of this miRNA in the treatment of related defects resulting from adverse cardiac remodeling.