Alan Lai

Epigenome-Wide Association Study Identifies Cardiac Gene Patterning and a Novel Class of Biomarkers for Heart Failure

Background: Biochemical DNA modification resembles a crucial regulatory layer among genetic information, environmental factors, and the transcriptome. To identify epigenetic susceptibility regions and novel biomarkers linked to myocardial dysfunction and heart failure, we performed the first multi-omics study in myocardial tissue and blood of patients with dilated cardiomyopathy and controls. Methods: Infinium human methylation 450 was used for high-density epigenome-wide mapping of DNA methylation in left-ventricular biopsies and whole peripheral blood of living probands. RNA deep sequencing was performed on the same samples in parallel. Whole-genome sequencing of all patients allowed exclusion of promiscuous genotype-induced methylation calls. Results: In the screening stage, we detected 59 epigenetic loci that are significantly associated with dilated cardiomyopathy (false discovery corrected P⩽0.05), with 3 of them reaching epigenome-wide significance at P⩽5×10−8. Twenty-seven (46%) of these loci could be replicated in independent cohorts, underlining the role of epigenetic regulation of key cardiac transcription regulators. Using a staged multi-omics study design, we link a subset of 517 epigenetic loci with dilated cardiomyopathy and cardiac gene expression. Furthermore, we identified distinct epigenetic methylation patterns that are conserved across tissues, rendering these CpGs novel epigenetic biomarkers for heart failure. Conclusions: The present study provides to our knowledge the first epigenome-wide association study in living patients with heart failure using a multi-omics approach.

Clinical genetics and outcome of left ventricular non-compaction cardiomyopathy

Aims In this study, we aimed to clinically and genetically characterize LVNC patients and investigate the prevalence of variants in known and novel LVNC disease genes. Introduction Left ventricular non-compaction cardiomyopathy (LVNC) is an increasingly recognized cause of heart failure, arrhythmia, thromboembolism, and sudden cardiac death. We sought here to dissect its genetic causes, phenotypic presentation and outcome. Methods and results In our registry with follow-up of in the median 61 months, we analysed 95 LVNC patients (68 unrelated index patients and 27 affected relatives; definite familial LVNC = 23.5%) by cardiac phenotyping, molecular biomarkers and exome sequencing. Cardiovascular events were significantly more frequent in LVNC patients compared with an age-matched group of patients with non-ischaemic dilated cardiomyopathy (hazard ratio = 2.481, P = 0.002). Stringent genetic classification according to ACMG guidelines revealed that TTN, LMNA, and MYBPC3 are the most prevalent disease genes (13 patients are carrying a pathogenic truncating TTN variant, odds ratio = 40.7, Confidence interval = 21.6-76.6, P < 0.0001, percent spliced in 76-100%). We also identified novel candidate genes for LVNC. For RBM20, we were able to perform detailed familial, molecular and functional studies. We show that the novel variant p.R634L in the RS domain of RBM20 co-segregates with LVNC, leading to titin mis-splicing as revealed by RNA sequencing of heart tissue in mutation carriers, protein analysis, and functional splice-reporter assays. Conclusion Our data demonstrate that the clinical course of symptomatic LVNC can be severe. The identified pathogenic variants and distribution of disease genes-a titin-related pathomechanism is found in every fourth patient-should be considered in genetic counselling of patients. Pathogenic variants in the nuclear proteins Lamin A/C and RBM20 were associated with worse outcome.

Genomic structural variations lead to dysregulation of important coding and non-coding RNA species in dilated cardiomyopathy

The transcriptome needs to be tightly regulated by mechanisms that include transcription factors, enhancers, and repressors as well as non‐coding RNAs. Besides this dynamic regulation, a large part of phenotypic variability of eukaryotes is expressed through changes in gene transcription caused by genetic variation. In this study, we evaluate genome‐wide structural genomic variants (SVs) and their association with gene expression in the human heart. We detected 3,898 individual SVs affecting all classes of gene transcripts (e.g., mRNA, miRNA, lncRNA) and regulatory genomic regions (e.g., enhancer or TFBS). In a cohort of patients (n = 50) with dilated cardiomyopathy (DCM), 80,635 non‐protein‐coding elements of the genome are deleted or duplicated by SVs, containing 3,758 long non‐coding RNAs and 1,756 protein‐coding transcripts. 65.3% of the SV‐eQTLs do not harbor a significant SNV‐eQTL, and for the regions with both classes of association, we find similar effect sizes. In case of deleted protein‐coding exons, we find downregulation of the associated transcripts, duplication events, however, do not show significant changes over all events. In summary, we are first to describe the genomic variability associated with SVs in heart failure due to DCM and dissect their impact on the transcriptome. Overall, SVs explain up to 7.5% of the variation of cardiac gene expression, underlining the importance to study human myocardial gene expression in the context of the individual genome. This has immediate implications for studies on basic mechanisms of cardiac maladaptation, biomarkers, and (gene) therapeutic studies alike.

Personalized Computer Simulation of Diastolic Function in Heart Failure

The search for a parameter representing left ventricular relaxation from non-invasive and invasive diagnostic tools has been extensive, since heart failure (HF) with preserved ejection fraction (HF-pEF) is a global health problem. We explore here the feasibility using patient-specific cardiac computer modeling to capture diastolic parameters in patients suffering from different degrees of systolic HF. Fifty eight patients with idiopathic dilated cardiomyopathy have undergone thorough clinical evaluation, including cardiac magnetic resonance imaging (MRI), heart catheterization, echocardiography, and cardiac biomarker assessment. A previously-introduced framework for creating multi-scale patient-specific cardiac models has been applied on all these patients. Novel parameters, such as global stiffness factor and maximum left ventricular active stress, representing cardiac active and passive tissue properties have been computed for all patients. Invasive pressure measurements from heart catheterization were then used to evaluate ventricular relaxation using the time constant of isovolumic relaxation Tau (τ). Parameters from heart catheterization and the multi-scale model have been evaluated and compared to patient clinical presentation. The model parameter global stiffness factor, representing diastolic passive tissue properties, is correlated significantly across the patient population with τ. This study shows that multi-modal cardiac models can successfully capture diastolic (dys) function, a prerequisite for future clinical trials on HF-pEF.

The Role of Quality Control in Targeted Next-generation Sequencing Library Preparation

Next-generation sequencing (NGS) is getting routinely used in the diagnosis of hereditary diseases, such as human cardiomyopathies. Hence, it is of utter importance to secure high quality sequencing data, enabling the identification of disease-relevant mutations or the conclusion of negative test results. During the process of sample preparation, each protocol for target enrichment library preparation has its own requirements for quality control (QC); however, there is little evidence on the actual impact of these guidelines on resulting data quality. In this study, we analyzed the impact of QC during the diverse library preparation steps of Agilent SureSelect XT target enrichment and Illumina sequencing. We quantified the parameters for a cohort of around 600 samples, which include starting amount of DNA, amount of sheared DNA, smallest and largest fragment size of the starting DNA; amount of DNA after the pre-PCR, and smallest and largest fragment size of the resulting DNA; as well as the amount of the final library, the corresponding smallest and largest fragment size, and the number of detected variants. Intriguingly, there is a high tolerance for variations in all QC steps, meaning that within the boundaries proposed in the current study, a considerable variance at each step of QC can be well tolerated without compromising NGS quality.