Genetic Dissection of Hypertrophic Cardiomyopathy with Myocardial RNA-Seq

Identifying biomarkers of ginseng medicines with different natures on heart failure

ETHNOPHARMACOLOGICAL RELEVANCE

The nature of Chinese medicine is a unique index to measure its efficacy. Generally, treating the hot syndrome with cold nature medicine and vice versa. Ginseng medicines, a renowned Chinese medicine known for its qi tonifying action, encompasses various herbal materials such as ginseng, red ginseng, and black ginseng (GS, RG, and BG, respectively), ginseng leaves (GL), and American ginseng (AG), which exhibited different natures, thought contained similar ginsenosides. This traditional effect of GS and RG "reinvigorate the pulse for relieving qi depletion". It is closely linked to anti-heart failure (HF), HF is a clinical manifestation of deficiency of "heart-qi". However, the elucidation of the mechanism underlying the anti-HF effects of ginseng medicines with different natures remains a significant challenge.

AIM OF THE STUDY

To elucidate pharmacological mechanisms underlying the effect of ginseng medicines on HF, and to identify biomarkers associated with their various natures. Furthermore, it provides the basis for the different applications of ginseng medicines with various natures.

MATERIALS AND METHODS

This study established a rat model of HF induced by isoproterenol (ISO) combined with a specific diet. Four representative hot/cold herbs were selected as compared references for the medicine natures. The divergent effects of these herbs on the HF model were investigated by analyzing RNA-seq data to identify genes expressed differentially. Additionally, pathways associated with medicine natures were obtained using KEGG. Furthermore, UPLC-QqQ-MS/MS, as well as ELISA, were used to measure indexes associated with the nervous system, energy metabolisms, and endocrinology systems, such as BNP, CK, IL-1, T3, T4, cAMP, cGMP, AD, adrenal hormones (DOC, CORT, and COR), progestogens (pregnenolone, P, 17-OH-PR, and 17-OH-P), androgens (DHEA, A4, and T), and estrogens hormones (E2).

RESULTS

All ginseng medicines demonstrated varying levels of efficacy in alleviating HF and GS exhibited a significant protective effect on HF. The ginseng medicines with qi tonifying primarily achieve their effect by enhancing the levels of adrenal hormones (DOC, CORT, and COR), T4, elevation of cAMP/cGMP, and activation of AchE. Warm nature qi tonifying ginseng medicines increased the levels of 17-OH-PR and P while decreasing 17-OH-P and the ratio of E2/T. On the other hand, cold nature qi tonifying ginseng medicines decreased the levels of A4 and T while increasing the ratio of E2/T.

CONCLUSION

Overall, the effects of warm nature ginseng medicines are stronger on HF compared to cold nature ginseng medicines. Our research firstly reported that the E2/T ratio, progestogens (17-OH-PR, 17-OH-P, and P), and androgens (A4 and T) have been identified as significant biomarkers for discerning the mechanism differences of ginseng medicines with differences natures in treatment of HF.

Genome-wide enhancer-associated tandem repeats are expanded in cardiomyopathy

BACKGROUND

Cardiomyopathy is a clinically and genetically heterogeneous heart condition that can lead to heart failure and sudden cardiac death in childhood. While it has a strong genetic basis, the genetic aetiology for over 50% of cardiomyopathy cases remains unknown.

METHODS

In this study, we analyse the characteristics of tandem repeats from genome sequence data of unrelated individuals diagnosed with cardiomyopathy from Canada and the United Kingdom (n = 1216) and compare them to those found in the general population. We perform burden analysis to identify genomic and epigenomic features that are impacted by rare tandem repeat expansions (TREs), and enrichment analysis to identify functional pathways that are involved in the TRE-associated genes in cardiomyopathy. We use Oxford Nanopore targeted long-read sequencing to validate repeat size and methylation status of one of the most recurrent TREs. We also compare the TRE-associated genes to those that are dysregulated in the heart tissues of individuals with cardiomyopathy.

FINDINGS

We demonstrate that tandem repeats that are rarely expanded in the general population are predominantly expanded in cardiomyopathy. We find that rare TREs are disproportionately present in constrained genes near transcriptional start sites, have high GC content, and frequently overlap active enhancer H3K27ac marks, where expansion-related DNA methylation may reduce gene expression. We demonstrate the gene silencing effect of expanded CGG tandem repeats in DIP2B through promoter hypermethylation. We show that the enhancer-associated loci are found in genes that are highly expressed in human cardiomyocytes and are differentially expressed in the left ventricle of the heart in individuals with cardiomyopathy.

INTERPRETATION

Our findings highlight the underrecognized contribution of rare tandem repeat expansions to the risk of cardiomyopathy and suggest that rare TREs contribute to ∼4% of cardiomyopathy risk.

FUNDING

Government of Ontario (RKCY), The Canadian Institutes of Health Research PJT 175329 (RKCY), The Azrieli Foundation (RKCY), SickKids Catalyst Scholar in Genetics (RKCY), The University of Toronto McLaughlin Centre (RKCY, SM), Ted Rogers Centre for Heart Research (SM), Data Sciences Institute at the University of Toronto (SM), The Canadian Institutes of Health Research PJT 175034 (SM), The Canadian Institutes of Health Research ENP 161429 under the frame of ERA PerMed (SM, RL), Heart and Stroke Foundation of Ontario & Robert M Freedom Chair in Cardiovascular Science (SM), Bitove Family Professorship of Adult Congenital Heart Disease (EO), Canada Foundation for Innovation (SWS, JR), Canada Research Chair (PS), Genome Canada (PS, JR), The Canadian Institutes of Health Research (PS).

Transcriptome analysis of human hypertrophic cardiomyopathy reveals inhibited cardiac development pathways in children

The epidemiological, etiological, and clinical characteristics vary greatly between pediatric (P-HCM) and adult (A-HCM) hypertrophic cardiomyopathy (HCM) patients, and the understanding of the heterogeneous pathogenesis mechanisms is insufficient to date. In this study, we aimed to comprehensively assess the respective transcriptome signatures and uncover the essential differences in gene expression patterns among A-HCM and P-HCM. The transcriptome data of adults were collected from public data (GSE89714), and novel pediatric data were first obtained by RNA sequencing from 14 P-HCM and 9 infantile donor heart samples. Our study demonstrates the common signatures of myofilament or protein synthesis and calcium ion regulation pathways in HCM. Mitochondrial function is specifically dysregulated in A-HCM, whereas the inhibition of cardiac developing networks typifies P-HCM. These findings not only distinguish the transcriptome characteristics in children and adults with HCM but also reveal the potential mechanism of the higher incidence of septal defects in P-HCM patients.

Feline myocardial transcriptome in health and in hypertrophic cardiomyopathy-A translational animal model for human disease

Hypertrophic cardiomyopathy (HCM) is the most common heart disease in cats, characterized by primary left ventricular hypertrophy. Feline HCM closely resembles human HCM and is suggested as translational animal model for the human disease. A genetic cause is established in humans and suspected for cats, but little is known about the gene expression and pathways involved in the pathogenesis of HCM. To investigate the myocardial transcriptome changes in HCM, RNA sequencing was conducted on left ventricle (LV) and left atrium (LA) samples of healthy cats and cats with HCM (each n = 5; 20 samples). Ingenuity Pathway Analysis was used to determine functional pathways, regulators, and networks. Distinct gene expression profiles were identified in the LV and LA of the feline healthy and HCM myocardium. Analysis of differentially expressed mRNAs (>2 fold; FDR < 0.01) found chamber-specific (LV vs. LA) expression in both healthy and HCM groups, with higher transcriptional activity in the LA. Genes that contribute to the distinct structure and function of each chamber in health and HCM were identified in the regional comparison. The gene expression profiles of HCM compared to healthy hearts revealed disease related genes, including THBS4 and KLHL33 (LV), FAM177B and THRSP (LA), the latter 3 have not been reported for the myocardium so far, as the top differently expressed genes in the HCM heart. Differently expressed genes and functional pathways found in the HCM heart are associated with cardiac remodeling and fibrosis, inflammation, microvascular changes, calcium signaling and cardiac metabolism, with some regional differences. RhoGDI-RhoGTPase signaling, integrin and ILK signaling pathways, the LXR/RXR pathway in the LA, and the PPARα/RXRα, HIF1α and CXCR4 pathways in the LV might be of particular importance in the HCM disease process. This study identified region-specific myocardial gene transcription patterns as well as novel genes and pathways associated with HCM.

Hypertrophic cardiomyopathy: an up-to-date snapshot of the clinical drug development pipeline

INTRODUCTION

Hypertrophic cardiomyopathy (HCM) is a complex cardiac disease with highly variable phenotypic expression and clinical course most often caused by sarcomeric gene mutations resulting in left ventricular hypertrophy, fibrosis, hypercontractility, and diastolic dysfunction. For almost 60 years, HCM has remained an orphan disease and still lacks a disease-specific treatment.

AREAS COVERED

This review summarizes recent preclinical and clinical trials with repurposed drugs and new emerging pharmacological and gene-based therapies for the treatment of HCM.

EXPERT OPINION

The off-label drugs routinely used alleviate symptoms but do not target the core pathophysiology of HCM or prevent or revert the phenotype. Recent advances in the genetics and pathophysiology of HCM led to the development of cardiac myosin adenosine triphosphatase inhibitors specifically directed to counteract the hypercontractility associated with HCM-causing mutations. Mavacamten, the first drug specifically developed for HCM successfully tested in a phase 3 trial, represents the major advance for the treatment of HCM. This opens new horizons for the development of novel drugs targeting HCM molecular substrates which hopefully modify the natural history of the disease. The role of current drugs in development and genetic-based approaches for the treatment of HCM are also discussed.

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.

CMYA5 is a novel interaction partner of FHL2 in cardiac myocytes

Four-and-a-half LIM domains protein 2 (FHL2) is an anti-hypertrophic adaptor protein that regulates cardiac myocyte signalling and function. Herein, we identified cardiomyopathy-associated 5 (CMYA5) as a novel FHL2 interaction partner in cardiac myocytes. In vitro pull-down assays demonstrated interaction between FHL2 and the N- and C-terminal regions of CMYA5. The interaction was verified in adult cardiac myocytes by proximity ligation assays. Immunofluorescence and confocal microscopy demonstrated co-localisation in the same subcellular compartment. The binding interface between FHL2 and CMYA5 was mapped by peptide arrays. Exposure of neonatal rat ventricular myocytes to a CMYA5 peptide covering one of the FHL2 interaction sites led to an increase in cell area at baseline, but a blunted response to chronic phenylephrine treatment. In contrast to wild-type hearts, loss or reduced FHL2 expression in Fhl2-targeted knockout mouse hearts or in a humanised mouse model of hypertrophic cardiomyopathy led to redistribution of CMYA5 into the perinuclear and intercalated disc region. Taken together, our results indicate a direct interaction of the two adaptor proteins FHL2 and CMYA5 in cardiac myocytes, which might impact subcellular compartmentation of CMYA5.

Circulating miR-499a-5p Is a Potential Biomarker of MYH7-Associated Hypertrophic Cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is the most common inherited myocardial disease with significant genetic and phenotypic heterogeneity. To search for novel biomarkers, which could increase the accuracy of HCM diagnosis and improve understanding of its phenotype formation, we analyzed the levels of circulating miRNAs—stable non-coding RNAs involved in post-transcriptional gene regulation. Performed high throughput sequencing of miRNAs in plasma of HCM patients and controls pinpointed miR-499a-5p as one of 35 miRNAs dysregulated in HCM. Further investigation on enlarged groups of individuals showed that its level was higher in carriers of pathogenic/likely pathogenic (P/LP) variants in MYH7 gene compared to controls (fold change, FC = 8.9; p < 0.0001). Just as important, carriers of variants in MYH7 gene were defined with higher miRNA levels than carriers of variants in the MYBPC3 gene (FC = 14.1; p = 0.0003) and other patients (FC = 4.1; p = 0.0008). The receiver operating characteristic analysis analysis showed the ability of miR-499a-5p to identify MYH7 variant carriers with the HCM phenotype with area under the curve value of 0.95 (95% confidence interval: 0.88−1.03, p = 0.0004); sensitivity and specificity were 0.86 and 0.91 (cut-off = 0.0014). Therefore, miR-499a-5p could serve as a circulating biomarker of HCM, caused by P/LP variants in MYH7 gene.

Whole genome sequencing delineates regulatory, copy number, and cryptic splice variants in early onset cardiomyopathy

Cardiomyopathy (CMP) is a heritable disorder. Over 50% of cases are gene-elusive on clinical gene panel testing. The contribution of variants in non-coding DNA elements that result in cryptic splicing and regulate gene expression has not been explored. We analyzed whole-genome sequencing (WGS) data in a discovery cohort of 209 pediatric CMP patients and 1953 independent replication genomes and exomes. We searched for protein-coding variants, and non-coding variants predicted to affect the function or expression of genes. Thirty-nine percent of cases harbored pathogenic coding variants in known CMP genes, and 5% harbored high-risk loss-of-function (LoF) variants in additional candidate CMP genes. Fifteen percent harbored high-risk regulatory variants in promoters and enhancers of CMP genes (odds ratio 2.25, p = 6.70 × 10-7 versus controls). Genes involved in α-dystroglycan glycosylation (FKTN, DTNA) and desmosomal signaling (DSC2, DSG2) were most highly enriched for regulatory variants (odds ratio 6.7-58.1). Functional effects were confirmed in patient myocardium and reporter assays in human cardiomyocytes, and in zebrafish CRISPR knockouts. We provide strong evidence for the genomic contribution of functionally active variants in new genes and in regulatory elements of known CMP genes to early onset CMP.

Cell cycle defects underlie childhood-onset cardiomyopathy associated with Noonan syndrome

Childhood-onset myocardial hypertrophy and cardiomyopathic changes are associated with significant morbidity and mortality in early life, particularly in patients with Noonan syndrome, a multisystemic genetic disorder caused by autosomal dominant mutations in genes of the Ras-MAPK pathway. Although the cardiomyopathy associated with Noonan syndrome (NS-CM) shares certain cardiac features with the hypertrophic cardiomyopathy caused by mutations in sarcomeric proteins (HCM), such as pathological myocardial remodeling, ventricular dysfunction, and increased risk for malignant arrhythmias, the clinical course of NS-CM significantly differs from HCM. This suggests a distinct pathophysiology that remains to be elucidated. Here, through analysis of sarcomeric myosin conformational states, histopathology, and gene expression in left ventricular myocardial tissue from NS-CM, HCM, and normal hearts complemented with disease modeling in cardiomyocytes differentiated from patient-derived PTPN11 N308S/+ induced pluripotent stem cells, we demonstrate distinct disease phenotypes between NS-CM and HCM and uncover cell cycle defects as a potential driver of NS-CM.

Combining whole exome sequencing with in silico analysis and clinical data to identify candidate variants in pediatric left ventricular noncompaction

BACKGROUND

Understanding the overall variant burden in pediatric patients with left ventricular noncompaction (LVNC) has clinical implications. Whole exome sequencing (WES) allows detection of coding variants in both candidate cardiomyopathy genes and those included on commercial panels. Other lines of evidence, including in silico analysis, are necessary to reduce the overwhelming number of variants to those most likely having a phenotypic impact.

METHODS

Five families, including five pediatric probands with LVNC, 5 other affected, and 10 unaffected family members, had WES performed, followed by bioinformatics filtering and Sanger sequencing. Review of the HGMD, variant classification by ACMG guidelines, and clinical information were used to further refine complex genotypes.

RESULTS

One nonsense and eleven missense variants were identified. In Family 1, affected siblings carried digenic heterozygous variants: E1350K-MYH7 and A276V-ANKRD1. The proband also carried heterozygous W143X-NRG1. Four affected members of Family 2 carried K184Q-MYH7 while unaffected members did not. In Family 3, homozygous A161T-MYH7 and heterozygous P4935T-OBSCN variants were identified in the proband with the latter being absent in his unaffected brother. In Family 4, proband's father and half-sibling have mild hypertrabeculation and carry T3796I-PLEC. The proband, carrying T3796I-PLEC and V2878A-OBSCN, demonstrated higher trabeculation burden. The proband in Family 5 carried four variants, R3247W-PLEC, C92Y-ERG, T1233M-NCOR2, and E54K-HIST1H4B. Application of ACMG criteria and clinical data revealed that W143X-NRG1, P4935T-OBSCN, and V2878A-OBSCN likely have no phenotypic role.

CONCLUSIONS

We report nine variants, including novel T3796I-PLEC and biallelic A161T-MYH7, likely contributing to phenotypes ranging from asymptomatic hypertrabeculation to severe LVNC with heart failure.

Case Report: Novel LIM domain-binding protein 3 (LDB3) mutations associated with hypertrophic cardiomyopathy family

Hypertrophic cardiomyopathy (HCM) is an autosomal dominant cardiomyopathy, which is one of the most common reasons for cardiac arrest in children or adolescents. It is characterized by ventricular hypertrophy (usually left ventricle), small ventricular cavity, and reduced ventricular diastolic compliance found by echocardiography in the absence of abnormal load (such as hypertension or aortic stenosis). HCM is usually caused by mutations in genes encoding sarcomere or sarcomere-related genes. Whole exome sequencing (WES) is performed to identify probable causative genes. Through WES, we identified LIM domain-binding protein 3 (LDB3) mutations (R547Q and P323S) respectively in an 11-year-old HCM girl and a 6-year-old HCM boy. Neural network analyses showed that the LDB3 (R547Q and P323S) mutation decreased its protein stability, with confidence scores of -0.9211 and -0.8967. The STRUM server also confirmed that the mutation decreased its protein stability. Thus, LDB3 mutation may be associated with heritable HCM. To our knowledge, this is the first time to report LDB3 heterozygous variants (R547Q and P323S) responsible for heritable HCM.

A Review on the Evolving Roles of MiRNA-Based Technologies in Diagnosing and Treating Heart Failure

MiRNA-regulated processes are pivotal in cardiovascular homeostasis and disease. These short non-coding RNAs have ideal properties that could be utilized as potential biomarkers; moreover, their functions as post-transcriptional regulators of mRNA make them interesting therapeutic targets. In this review, we summarize the current state of miRNA-based biomarkers in a variety of diseases leading to heart failure, as well as provide an outlook on developing miRNA-based therapies in the heart failure field.

The Emerging Role of Epigenetics in Therapeutic Targeting of Cardiomyopathies

Cardiomyopathies (CMPs) are a heterogeneous group of myocardial diseases accountable for the majority of cases of heart failure (HF) and/or sudden cardiac death (SCD) worldwide. With the recent advances in genomics, the original classification of CMPs on the basis of morphological and functional criteria (dilated (DCM), hypertrophic (HCM), restrictive (RCM), and arrhythmogenic ventricular cardiomyopathy (AVC)) was further refined into genetic (inherited or familial) and acquired (non-inherited or secondary) forms. Despite substantial progress in the identification of novel CMP-associated genetic variations, as well as improved clinical recognition diagnoses, the functional consequences of these mutations and the exact details of the signaling pathways leading to hypertrophy, dilation, and/or contractile impairment remain elusive. To date, global research has mainly focused on the genetic factors underlying CMP pathogenesis. However, growing evidence shows that alterations in molecular mediators associated with the diagnosis of CMPs are not always correlated with genetic mutations, suggesting that additional mechanisms, such as epigenetics, may play a role in the onset or progression of CMPs. This review summarizes published findings of inherited CMPs with a specific focus on the potential role of epigenetic mechanisms in regulating these cardiac disorders.

Systemic Bioinformatic Analyses of Nuclear-Encoded Mitochondrial Genes in Hypertrophic Cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is an autosomal dominant disease and mitochondria plays a key role in the progression in HCM. Here, we analyzed the expression pattern of nuclear-encoded mitochondrial genes (NMGenes) in HCM and found that the expression of NMGenes was significantly changed. A total of 316 differentially expressed NMGenes (DE-NMGenes) were identified. Pathway enrichment analyses showed that energy metabolism-related pathways such as "pyruvate metabolism" and "fatty acid degradation" were dysregulated, which highlighted the importance of energy metabolism in HCM. Next, we constructed a protein-protein interaction network based on 316 DE-NMGenes and identified thirteen hubs. Then, a total of 17 TFs (transcription factors) were predicted to potentially regulate the expression of 316 DE-NMGenes according to iRegulon, among which 8 TFs were already found involved in pathological hypertrophy. The remaining TFs (like GATA1, GATA5, and NFYA) were good candidates for further experimental verification. Finally, a mouse model of transverse aortic constriction (TAC) was established to validate the genes and results showed that DDIT4, TKT, CLIC1, DDOST, and SNCA were all upregulated in TAC mice. The present study represents the first effort to evaluate the global expression pattern of NMGenes in HCM and provides innovative insight into the molecular mechanism of HCM.

Statistical integration of two omics datasets using GO2PLS

BACKGROUND

Nowadays, multiple omics data are measured on the same samples in the belief that these different omics datasets represent various aspects of the underlying biological systems. Integrating these omics datasets will facilitate the understanding of the systems. For this purpose, various methods have been proposed, such as Partial Least Squares (PLS), decomposing two datasets into joint and residual subspaces. Since omics data are heterogeneous, the joint components in PLS will contain variation specific to each dataset. To account for this, Two-way Orthogonal Partial Least Squares (O2PLS) captures the heterogeneity by introducing orthogonal subspaces and better estimates the joint subspaces. However, the latent components spanning the joint subspaces in O2PLS are linear combinations of all variables, while it might be of interest to identify a small subset relevant to the research question. To obtain sparsity, we extend O2PLS to Group Sparse O2PLS (GO2PLS) that utilizes biological information on group structures among variables and performs group selection in the joint subspace.

RESULTS

The simulation study showed that introducing sparsity improved the feature selection performance. Furthermore, incorporating group structures increased robustness of the feature selection procedure. GO2PLS performed optimally in terms of accuracy of joint score estimation, joint loading estimation, and feature selection. We applied GO2PLS to datasets from two studies: TwinsUK (a population study) and CVON-DOSIS (a small case-control study). In the first, we incorporated biological information on the group structures of the methylation CpG sites when integrating the methylation dataset with the IgG glycomics data. The targeted genes of the selected methylation groups turned out to be relevant to the immune system, in which the IgG glycans play important roles. In the second, we selected regulatory regions and transcripts that explained the covariance between regulomics and transcriptomics data. The corresponding genes of the selected features appeared to be relevant to heart muscle disease.

CONCLUSIONS

GO2PLS integrates two omics datasets to help understand the underlying system that involves both omics levels. It incorporates external group information and performs group selection, resulting in a small subset of features that best explain the relationship between two omics datasets for better interpretability.