NIMA-related kinase 9 regulates the phosphorylation of the essential myosin light chain in the heart.

To adapt to changing hemodynamic demands, regulatory mechanisms modulate actin-myosin-kinetics by calcium-dependent and -independent mechanisms. We investigate the posttranslational modification of human essential myosin light chain (ELC) and identify NIMA-related kinase 9 (NEK9) to interact with ELC. NEK9 is highly expressed in the heart and the interaction with ELC is calcium-dependent. Silencing of NEK9 results in blunting of calcium-dependent ELC-phosphorylation. CRISPR/Cas9-mediated disruption of NEK9 leads to cardiomyopathy in zebrafish. Binding to ELC is mediated via the protein kinase domain of NEK9. A causal relationship between NEK9 activity and ELC-phosphorylation is demonstrated by genetic sensitizing in-vivo. Finally, we observe significantly upregulated ELC-phosphorylation in dilated cardiomyopathy patients and provide a unique map of human ELC-phosphorylation-sites. In summary, NEK9-mediated ELC-phosphorylation is a calcium-dependent regulatory system mediating cardiac contraction and inotropy.

Multi-omics assessment of dilated cardiomyopathy using non-negative matrix factorization.

Dilated cardiomyopathy (DCM), a myocardial disease, is heterogeneous and often results in heart failure and sudden cardiac death. Unavailability of cardiac tissue has hindered the comprehensive exploration of gene regulatory networks and nodal players in DCM. In this study, we carried out integrated analysis of transcriptome and methylome data using non-negative matrix factorization from a cohort of DCM patients to uncover underlying latent factors and covarying features between whole-transcriptome and epigenome omics datasets from tissue biopsies of living patients. DNA methylation data from Infinium HM450 and mRNA Illumina sequencing of n = 33 DCM and n = 24 control probands were filtered, analyzed and used as input for matrix factorization using R NMF package. Mann-Whitney U test showed 4 out of 5 latent factors are significantly different between DCM and control probands (P<0.05). Characterization of top 10% features driving each latent factor showed a significant enrichment of biological processes known to be involved in DCM pathogenesis, including immune response (P = 3.97E-21), nucleic acid binding (P = 1.42E-18), extracellular matrix (P = 9.23E-14) and myofibrillar structure (P = 8.46E-12). Correlation network analysis revealed interaction of important sarcomeric genes like Nebulin, Tropomyosin alpha-3 and ERC-protein 2 with CpG methylation of ATPase Phospholipid Transporting 11A0, Solute Carrier Family 12 Member 7 and Leucine Rich Repeat Containing 14B, all with significant P values associated with correlation coefficients >0.7. Using matrix factorization, multi-omics data derived from human tissue samples can be integrated and novel interactions can be identified. Hypothesis generating nature of such analysis could help to better understand the pathophysiology of complex traits such as DCM.

Doxorubicin induces cardiotoxicity in a pluripotent stem cell model of aggressive B cell lymphoma cancer patients.

Cancer therapies with anthracyclines have been shown to induce cardiovascular complications. The aims of this study were to establish an in vitro induced pluripotent stem cell model (iPSC) of anthracycline-induced cardiotoxicity (ACT) from patients with an aggressive form of B-cell lymphoma and to examine whether doxorubicin (DOX)-treated ACT-iPSC cardiomyocytes (CM) can recapitulate the clinical features exhibited by patients, and thus help uncover a DOX-dependent pathomechanism. ACT-iPSC CM generated from individuals with CD20+ B-cell lymphoma who had received high doses of DOX and suffered cardiac dysfunction were studied and compared to control-iPSC CM from cancer survivors without cardiac symptoms. In cellular studies, ACT-iPSC CM were persistently more susceptible to DOX toxicity including augmented disorganized myofilament structure, changed mitochondrial shape, and increased apoptotic events. Consistently, ACT-iPSC CM and cardiac fibroblasts isolated from fibrotic human ACT myocardium exhibited higher DOX-dependent reactive oxygen species. In functional studies, Ca2+ transient amplitude of ACT-iPSC CM was reduced compared to control cells, and diastolic sarcoplasmic reticulum Ca2+ leak was DOX-dependently increased. This could be explained by overactive CaMKIIδ in ACT CM. Together with DOX-dependent augmented proarrhythmic cellular triggers and prolonged action potentials in ACT CM, this suggests a cellular link to arrhythmogenic events and contractile dysfunction especially found in ACT engineered human myocardium. CamKIIδ inhibition prevented proarrhythmic triggers in ACT. In contrast, control CM upregulated SERCA2a expression in a DOX-dependent manner, possibly to avoid heart failure conditions. In conclusion, we developed the first human patient-specific stem cell model of DOX-induced cardiac dysfunction from patients with B-cell lymphoma. Our results suggest that DOX-induced stress resulted in arrhythmogenic events associated with contractile dysfunction and finally in heart failure after persistent stress activation in ACT patients.

SLM2 Is A Novel Cardiac Splicing Factor Involved in Heart Failure due to Dilated Cardiomyopathy.

Alternative mRNA splicing is a fundamental process to increase the versatility of the genome. In humans, cardiac mRNA splicing is involved in the pathophysiology of heart failure. Mutations in the splicing factor RNA binding motif protein 20 (RBM20) cause severe forms of cardiomyopathy. To identify novel cardiomyopathy-associated splicing factors, RNA-seq and tissue-enrichment analyses were performed, which identified up-regulated expression of Sam68-Like mammalian protein 2 (SLM2) in the left ventricle of dilated cardiomyopathy (DCM) patients. In the human heart, SLM2 binds to important transcripts of sarcomere constituents, such as those encoding myosin light chain 2 (MYL2), troponin I3 (TNNI3), troponin T2 (TNNT2), tropomyosin 1/2 (TPM1/2), and titin (TTN). Mechanistically, SLM2 mediates intron retention, prevents exon exclusion, and thereby mediates alternative splicing of the mRNA regions encoding the variable proline-, glutamate-, valine-, and lysine-rich (PEVK) domain and another part of the I-band region of titin. In summary, SLM2 is a novel cardiac splicing regulator with essential functions for maintaining cardiomyocyte integrity by binding to and processing the mRNAs of essential cardiac constituents such as titin.

Energy Metabolites as Biomarkers in Ischemic and Dilated Cardiomyopathy.

With more than 25 million people affected, heart failure (HF) is a global threat. As energy production pathways are known to play a pivotal role in HF, we sought here to identify key metabolic changes in ischemic- and non-ischemic HF by using a multi-OMICS approach. Serum metabolites and mRNAseq and epigenetic DNA methylation profiles were analyzed from blood and left ventricular heart biopsy specimens of the same individuals. In total we collected serum from n = 82 patients with Dilated Cardiomyopathy (DCM) and n = 51 controls in the screening stage. We identified several metabolites involved in glycolysis and citric acid cycle to be elevated up to 5.7-fold in DCM (p = 1.7 × 10-6). Interestingly, cardiac mRNA and epigenetic changes of genes encoding rate-limiting enzymes of these pathways could also be found and validated in our second stage of metabolite assessment in n = 52 DCM, n = 39 ischemic HF and n = 57 controls. In conclusion, we identified a new set of metabolomic biomarkers for HF. We were able to identify underlying biological cascades that potentially represent suitable intervention targets.

microRNA neural networks improve diagnosis of acute coronary syndrome (ACS).

BACKGROUND: Cardiac troponins are the preferred biomarkers of acute myocardial infarction. Despite superior sensitivity, serial testing of Troponins to identify patients suffering acute coronary syndromes is still required in many cases to overcome limited specificity. Moreover, unstable angina pectoris relies on reported symptoms in the troponin-negative group. In this study, we investigated genome-wide miRNA levels in a prospective cohort of patients with clinically suspected ACS and determined their diagnostic value by applying an in silico neural network. METHODS: PAXgene blood and serum samples were drawn and hsTnT was measured in patients at initial presentation to our Chest-Pain Unit. After clinical and diagnostic workup, patients were adjudicated by senior cardiologists in duty to their final diagnosis: STEMI, NSTEMI, unstable angina pectoris and non-ACS patients. ACS patients and a cohort of healthy controls underwent deep transcriptome sequencing. Machine learning was implemented to construct diagnostic miRNA classifiers. RESULTS: We developed a neural network model which incorporates 34 validated ACS miRNAs, showing excellent classification results. By further developing additional machine learning models and selecting the best miRNAs, we achieved an accuracy of 0.96 (95% CI 0.96-0.97), sensitivity of 0.95, specificity of 0.96 and AUC of 0.99. The one-point hsTnT value reached an accuracy of 0.89, sensitivity of 0.82, specificity of 0.96, and AUC of 0.96. CONCLUSIONS: Here we show the concept of neural network based biomarkers for ACS. This approach also opens the possibility to include multi-modal data points to further increase precision and perform classification of other ACS differential diagnoses.

Epigenetic Regulation of Alternative mRNA Splicing in Dilated Cardiomyopathy.

In recent years, the genetic architecture of dilated cardiomyopathy (DCM) has been more thoroughly elucidated. However, there is still insufficient knowledge on the modifiers and regulatory principles that lead to the failure of myocardial function. The current study investigates the association of epigenome-wide DNA methylation and alternative splicing, both of which are important regulatory principles in DCM. We analyzed screening and replication cohorts of cases and controls and identified distinct transcriptomic patterns in the myocardium that differ significantly, and we identified a strong association of intronic DNA methylation and flanking exons usage (p < 2 × 10-16). By combining differential exon usage (DEU) and differential methylation regions (DMR), we found a significant change of regulation in important sarcomeric and other DCM-associated pathways. Interestingly, inverse regulation of Titin antisense non-coding RNA transcript splicing and DNA methylation of a locus reciprocal to TTN substantiate these findings and indicate an additional role for non-protein-coding transcripts. In summary, this study highlights for the first time the close interrelationship between genetic imprinting by DNA methylation and the transport of this epigenetic information towards the dynamic mRNA splicing landscape. This expands our knowledge of the genome-environment interaction in DCM besides simple gene expression regulation.

m6A-mRNA methylation regulates cardiac gene expression and cellular growth.

Conceptually similar to modifications of DNA, mRNAs undergo chemical modifications, which can affect their activity, localization, and stability. The most prevalent internal modification in mRNA is the methylation of adenosine at the N6-position (m6A). This returns mRNA to a role as a central hub of information within the cell, serving as an information carrier, modifier, and attenuator for many biological processes. Still, the precise role of internal mRNA modifications such as m6A in human and murine-dilated cardiac tissue remains unknown. Transcriptome-wide mapping of m6A in mRNA allowed us to catalog m6A targets in human and murine hearts. Increased m6A methylation was found in human cardiomyopathy. Knockdown and overexpression of the m6A writer enzyme Mettl3 affected cell size and cellular remodeling both in vitro and in vivo. Our data suggest that mRNA methylation is highly dynamic in cardiomyocytes undergoing stress and that changes in the mRNA methylome regulate translational efficiency by affecting transcript stability. Once elucidated, manipulations of methylation of specific m6A sites could be a powerful approach to prevent worsening of cardiac function.

Fast and simple protein-alignment-guided assembly of orthologous gene families from microbiome sequencing reads.

BACKGROUND: Microbiome sequencing projects typically collect tens of millions of short reads per sample. Depending on the goals of the project, the short reads can either be subjected to direct sequence analysis or be assembled into longer contigs. The assembly of whole genomes from metagenomic sequencing reads is a very difficult problem. However, for some questions, only specific genes of interest need to be assembled. This is then a gene-centric assembly where the goal is to assemble reads into contigs for a family of orthologous genes. METHODS: We present a new method for performing gene-centric assembly, called protein-alignment-guided assembly, and provide an implementation in our metagenome analysis tool MEGAN. Genes are assembled on the fly, based on the alignment of all reads against a protein reference database such as NCBI-nr. Specifically, the user selects a gene family based on a classification such as KEGG and all reads binned to that gene family are assembled. RESULTS: Using published synthetic community metagenome sequencing reads and a set of 41 gene families, we show that the performance of this approach compares favorably with that of full-featured assemblers and that of a recently published HMM-based gene-centric assembler, both in terms of the number of reference genes detected and of the percentage of reference sequence covered. CONCLUSIONS: Protein-alignment-guided assembly of orthologous gene families complements whole-metagenome assembly in a new and very useful way.

MEGAN Community Edition - Interactive Exploration and Analysis of Large-Scale Microbiome Sequencing Data.

There is increasing interest in employing shotgun sequencing, rather than amplicon sequencing, to analyze microbiome samples. Typical projects may involve hundreds of samples and billions of sequencing reads. The comparison of such samples against a protein reference database generates billions of alignments and the analysis of such data is computationally challenging. To address this, we have substantially rewritten and extended our widely-used microbiome analysis tool MEGAN so as to facilitate the interactive analysis of the taxonomic and functional content of very large microbiome datasets. Other new features include a functional classifier called InterPro2GO, gene-centric read assembly, principal coordinate analysis of taxonomy and function, and support for metadata. The new program is called MEGAN Community Edition (CE) and is open source. By integrating MEGAN CE with our high-throughput DNA-to-protein alignment tool DIAMOND and by providing a new program MeganServer that allows access to metagenome analysis files hosted on a server, we provide a straightforward, yet powerful and complete pipeline for the analysis of metagenome shotgun sequences. We illustrate how to perform a full-scale computational analysis of a metagenomic sequencing project, involving 12 samples and 800 million reads, in less than three days on a single server. All source code is available here: https://github.com/danielhuson/megan-ce.

Characterization of the Gut Microbial Community of Obese Patients Following a Weight-Loss Intervention Using Whole Metagenome Shotgun Sequencing.

BACKGROUND/OBJECTIVES: Cross-sectional studies suggested that obesity is promoted by the gut microbiota. However, longitudinal data on taxonomic and functional changes in the gut microbiota of obese patients are scarce. The aim of this work is to study microbiota changes in the course of weight loss therapy and the following year in obese individuals with or without co-morbidities, and to asses a possible predictive value of the gut microbiota with regard to weight loss maintenance. SUBJECTS/METHODS: Sixteen adult patients, who followed a 52-week weight-loss program comprising low calorie diet, exercise and behavioral therapy, were selected according to their weight-loss course. Over two years, anthropometric and metabolic parameters were assessed and microbiota from stool samples was functionally and taxonomically analyzed using DNA shotgun sequencing. RESULTS: Overall the microbiota responded to the dietetic and lifestyle intervention but tended to return to the initial situation both at the taxonomical and functional level at the end of the intervention after one year, except for an increase in Akkermansia abundance which remained stable over two years (12.7x103 counts, 95%CI: 322-25100 at month 0; 141x103 counts, 95%CI: 49-233x103 at month 24; p = 0.005). The Firmicutes/Bacteroidetes ratio was higher in obese subjects with metabolic syndrome (0.64, 95%CI: 0.34-0.95) than in the "healthy obese" (0.27, 95%CI: 0.08-0.45, p = 0.04). Participants, who succeeded in losing their weight consistently over the two years, had at baseline a microbiota enriched in Alistipes, Pseudoflavonifractor and enzymes of the oxidative phosphorylation pathway compared to patients who were less successful in weight reduction. CONCLUSIONS: Successful weight reduction in the obese is accompanied with increased Akkermansia numbers in feces. Metabolic co-morbidities are associated with a higher Firmicutes/Bacteroidetes ratio. Most interestingly, microbiota differences might allow discrimination between successful and unsuccessful weight loss prior to intervention.