Abstract 439: Trajectory Analysis of Left Ventricular Dimensions From Biobank Data Uncovers Novel Genetic Associations

Background: Approximately 6 million adults in the United States have heart failure. The progression of heart failure is variable arising from differences in sex, age, genetic background including ancestry, and medication response. Many population-based genetic studies of heart failure have been cross-sectional in nature, failing to gain additional power from longitudinal analyses. As heart failure is known to change over time, using longitudinal data trajectories as a quantitative trait will increase power in genome wide association studies (GWAS).

Methods: We used the electronic health record in a racially and ethnically diverse medical biobank from a single, metropolitan US center. We used whole genome data from 896 unrelated participants analyzed, including 494 who had at least 1 electrocardiogram and 324 who had more than 1 echocardiogram (average of 3 observations per person). A mixture model based semiparametric latent growth curve model was used to cluster outcome measures used for genome-wide analyses.

Results: GWAS on the trajectory probability of QTc interval identified significant associations with variants in regulatory regions proximal to the WLS gene, which encodes Wntless, a Wnt ligand secretion mediator. WLS was previously associated with QTc and myocardial infarction, thus confirming the power of the method. GWAS on the trajectory probability of left ventricular diameter (LVIDd) identified significant associations with variants in regulatory regions near MYO10 , which encodes unconventional Myosin-10. MYO10 was previously associated with obesity and metabolic syndrome.

Conclusions: This is the first study to show an association with variants in or near MYO10 and left ventricular dimension changes over time. Further, we found that using trajectory probabilities can provide increased power to find novel associations with longitudinal data. This reduces the need for larger cohorts, and increases yield from smaller, well-phenotyped cohorts, such as those found in biobanks. This approach should be useful in the study of rare diseases and underrepresented populations.

Pathogenic and Uncertain Genetic Variants Have Clinical Cardiac Correlates in Diverse Biobank Participants

Background

Genome sequencing coupled with electronic heath record data can uncover medically important genetic variation. Interpretation of rare genetic variation and its role in mediating cardiovascular phenotypes is confounded by variants of uncertain significance.

Methods and Results

We analyzed the whole genome sequence of 900 racially and ethnically diverse biobank participants selected from a single US center. Participants were equally divided among European, African, Hispanic, and mixed races/ethnicities. We evaluated the American College of Medical Genetics and Genomics medically actionable list of 59 genes, focusing on the cardiac genes. Variation was interpreted using the most recent reports in ClinVar, a database of medically relevant human variation. We identified 19 individuals with pathogenic or likely pathogenic variants in cardiac actionable genes (2%) and found evidence of related clinical correlates in the electronic health record. Participants of African ancestry, compared with those of European ancestry, had more variants of uncertain significance in the medically actionable genes including the 30 cardiac actionable genes, even when normalized to total variant count per person. Longitudinal measures of left ventricle size from ≈400 biobank participants (1723 patient‐years) were correlated with genetic findings. The presence of ≥1 uncertain variant in the actionable cardiac genes and a cardiomyopathy diagnosis correlated with increased left ventricular internal diameter in diastole and in systole. In particular, MYBPC3 was identified as a gene with excess variants of uncertain significance.

Conclusions

These data indicate that a subset of uncertain genetic variants may confer risk and should not be considered benign.

Association of Cardiomyopathy With MYBPC3 D389V and MYBPC3Δ25bpIntronic Deletion in South Asian Descendants

Importance

The genetic variant MYBPC3Δ25bp occurs in 4% of South Asian descendants, with an estimated 100 million carriers worldwide. MYBPC3 Δ25bp has been linked to cardiomyopathy and heart failure. However, the high prevalence of MYBPC3Δ25bp suggests that other stressors act in concert with MYBPC3Δ25bp.

Objective

To determine whether there are additional genetic factors that contribute to the cardiomyopathic expression of MYBPC3Δ25bp.

Design, Setting, andParticipants

South Asian individuals living in the United States were screened for MYBPC3Δ25bp, and a subgroup was clinically evaluated using electrocardiograms and echocardiograms at Loyola University, Chicago, Illinois, between January 2015 and July 2016.

Main Outcomes and Measures

Next-generation sequencing of 174 cardiovascular disease genes was applied to identify additional modifying gene mutations and correlate genotype-phenotype parameters. Cardiomyocytes derived from human-induced pluripotent stem cells were established and examined to assess the role of MYBPC3Δ25bp.

Results

In this genotype-phenotype study, individuals of South Asian descent living in the United States from both sexes (36.23% female) with a mean population age of 48.92 years (range, 18-84 years) were recruited. Genetic screening of 2401 US South Asian individuals found an MYBPC3Δ25bpcarrier frequency of 6%. A higher frequency of missense TTN variation was found in MYBPC3Δ25bp carriers compared with noncarriers, identifying distinct genetic backgrounds within the MYBPC3Δ25bp carrier group. Strikingly, 9.6% of MYBPC3Δ25bp carriers also had a novel MYBPC3 variant, D389V. Family studies documented D389V was in tandem on the same allele as MYBPC3Δ25bp, and D389V was only seen in the presence of MYBPC3Δ25bp. In contrast to MYBPC3Δ25bp, MYBPC3Δ25bp/D389V was associated with hyperdynamic left ventricular performance (mean [SEM] left ventricular ejection fraction, 66.7 [0.7%]; left ventricular fractional shortening, 36.6 [0.6%]; P < .03) and stem cell-derived cardiomyocytes exhibited cellular hypertrophy with abnormal Ca2+ transients.

Conclusions and Relevance

MYBPC3Δ25bp/D389V is associated with hyperdynamic features, which are an early finding in hypertrophic cardiomyopathy and thought to reflect an unfavorable energetic state. These findings support that a subset of MYBPC3Δ25bp carriers, those with D389V, account for the increased risk attributed to MYBPC3Δ25bp.

Harmonizing Clinical Sequencing and Interpretation for the eMERGE III Network

The advancement of precision medicine requires new methods to coordinate and deliver genetic data from heterogeneous sources to physicians and patients. The eMERGE III Network enrolled >25,000 participants from biobank and prospective cohorts of predominantly healthy individuals for clinical genetic testing to determine clinically actionable findings. The network developed protocols linking together the 11 participant collection sites and 2 clinical genetic testing laboratories. DNA capture panels targeting 109 genes were used for testing of DNA and sample collection, data generation, interpretation, reporting, delivery, and storage were each harmonized. A compliant and secure network enabled ongoing review and reconciliation of clinical interpretations, while maintaining communication and data sharing between clinicians and investigators. A total of 202 individuals had positive diagnostic findings relevant to the indication for testing and 1,294 had additional/secondary findings of medical significance deemed to be returnable, establishing data return rates for other testing endeavors. This study accomplished integration of structured genomic results into multiple electronic health record (EHR) systems, setting the stage for clinical decision support to enable genomic medicine. Further, the established processes enable different sequencing sites to harmonize technical and interpretive aspects of sequencing tests, a critical achievement toward global standardization of genomic testing. The eMERGE protocols and tools are available for widespread dissemination.

Abstract 469: Genomic Context Predicts Dilated but Not Hypertrophic Cardiomyopathy

Rationale: Cardiomyopathies are a major cause of heart failure and have a strong heritable component.

Objective: Determine the role of both common and rare nonsynonymous genetic variation in hypertrophic (HCM) and dilated familial cardiomyopathy (DCM).

Methods and Results: Whole genome sequencing was used to determine common and rare nonsynonymous genetic variation in familial cases of HCM (n=56) or DCM (n=71). Variation was evaluated in 102 cardiomyopathy genes routinely assayed in clinical and commercial gene testing panels. We used cardiac gene expression data from GTEx (Genotype-Tissue Expression database) to define additional genes expressed in the heart. The number of nonsynonymous single nucleotide variants (nsSNVs), the majority of which were missense variants, was correlated with echocardiographic measurements. Principal component analysis (PCA) of left ventricular measures separated HCM and DCM. Regression of the first principal component using all nonsynonymous single nucleotide variants (nsSNVs) in cardiomyopathy genes showed the number of nsSNVs predicted DCM but not HCM. DCM probability in the cohort significantly increased as the number of cardiomyopathy gene nsSNVs increased (p<0.02). The increase in nsSNVs in cardiomyopathy genes significantly associated with reduced left ventricular ejection fraction and increased left ventricular diameter in DCM. Resampling methods identified genes with deviant cumulative allele frequencies, identifying potential modifier genes for cardiomyopathy.

Conclusions: DCM subjects carry a greater burden of nsSNVs in cardiomyopathy genes. This genomic burden translates to impaired systolic cardiac function in DCM. In contrast, nsSNV burden in cardiomyopathy genes did not correlate with the probability or manifestation of left ventricular measures in HCM. These findings support a complex inheritance for DCM where increased variation in cardiomyopathy genes creates a genetic background that predisposes to DCM and increased disease severity. The distinct genetic landscapes of HCM and DCM suggest that greater genetic variation in cardiac genes provokes unfavorable ventricular remodeling with reduced systolic function.

Prelamin A causes aberrant myonuclear arrangement and results in muscle fiber weakness

Physiological and premature aging are frequently associated with an accumulation of prelamin A, a precursor of lamin A, in the nuclear envelope of various cell types. Here, we aimed to underpin the hitherto unknown mechanisms by which prelamin A alters myonuclear organization and muscle fiber function. By experimentally studying membrane-permeabilized myofibers from various transgenic mouse lines, our results indicate that, in the presence of prelamin A, the abundance of nuclei and myosin content is markedly reduced within muscle fibers. This leads to a concept by which the remaining myonuclei are very distant from each other and are pushed to function beyond their maximum cytoplasmic capacity, ultimately inducing muscle fiber weakness.

Abstract 116: Mybphl is a Novel Myofilament Component Implicated in Arrhythmia and Dilated Cardiomyopathy

Dilated cardiomyopathy is under significant genetic influence. Using whole-genome sequencing, we found a premature stop variant (R255X) in the MYBPHL gene in a family with DCM and conduction system disease. A Mybphl null mouse model revealed systolic dysfunction with atrial and ventricular conduction system (VCS) abnormalities. We found MyBP-HL protein highly expressed in both human and mouse atria and in a subset of ventricular cardiomyocytes. We hypothesize that MyBP-HL regulates VCS function and loss of MyBP-HL causes ventricular arrhythmia and cardiac dysfunction. We now examined hearts and isolated VCMs from Mybphl WT, het, and null mice using immunofluorescence microscopy (IFM). Clusters of MyBP-HL+ ventricular cardiomyocytes were evident in WT hearts (~1/200,000), with only 10% as many MyBP-HL+ cardiomyocytes in het samples. In WT hearts, MyBP-HL+ cells localized in the right ventricular free wall, with significantly fewer MyBP-HL+ cells in the right ventricular free wall in het hearts. The VCS marker contactin-2 (Cntn2) was used to assess MyBP-HL co-localization in the VCS. MyBP-HL co-localized with Cntn2+ cardiomyocytes in the transition zone between the atria and ventricle near the AV node, as well as in a subset of Cntn2+ Purkinje fibers. MyBP-HL+ ventricular cardiomyocytes were also found outside the VCS, suggesting other roles for those cells in regulating ventricular function. To discern MyBP-HL myofilament localization, we used super-resolution structured illumination microscopy on atrial cardiomyocytes. MyBP-HL co-localized with cMyBP-C in the C-zone of the A-band but also occupied space in the inner portion of the A-band near the M-line. This localization was not altered in het mice; however, the ratio of MyBP-HL in the inner portion of the A-band was diminished in het mice compared to WT. This suggests that MyBP-HL may bind the light meromyosin coiled-coil tail that originates from the M-line. Following the discovery that loss of MyBP-HL is associated with DCM and arrhythmia, we identified that MyBP-HL associates with the VCS and reduced MyBP-HL levels reduce the amount of MyBP-HL positive VCMs. These data suggest that loss of MyBP-HL may impair the development and function of the VCS, resulting in arrhythmias.

Efficient exon skipping of SGCG mutations mediated by phosphorodiamidate morpholino oligomers

Exon skipping uses chemically modified antisense oligonucleotides to modulate RNA splicing. Therapeutically, exon skipping can bypass mutations and restore reading frame disruption by generating internally truncated, functional proteins to rescue the loss of native gene expression. Limb-girdle muscular dystrophy type 2C is caused by autosomal recessive mutations in the SGCG gene, which encodes the dystrophin-associated protein γ-sarcoglycan. The most common SGCG mutations disrupt the transcript reading frame abrogating γ-sarcoglycan protein expression. In order to treat most SGCG gene mutations, it is necessary to skip 4 exons in order to restore the SGCG transcript reading frame, creating an internally truncated protein referred to as Mini-Gamma. Using direct reprogramming of human cells with MyoD, myogenic cells were tested with 2 antisense oligonucleotide chemistries, 2'-O-methyl phosphorothioate oligonucleotides and vivo-phosphorodiamidate morpholino oligomers, to induce exon skipping. Treatment with vivo-phosphorodiamidate morpholino oligomers demonstrated efficient skipping of the targeted exons and corrected the mutant reading frame, resulting in the expression of a functional Mini-Gamma protein. Antisense-induced exon skipping of SGCG occurred in normal cells and those with multiple distinct SGCG mutations, including the most common 521ΔT mutation. These findings demonstrate a multiexon-skipping strategy applicable to the majority of limb-girdle muscular dystrophy 2C patients.

Novel nesprin-1 mutations associated with dilated cardiomyopathy cause nuclear envelope disruption and defects in myogenesis

Nesprins-1 and -2 are highly expressed in skeletal and cardiac muscle and together with SUN (Sad1p/UNC84)-domain containing proteins and lamin A/C form the LInker of Nucleoskeleton-and-Cytoskeleton (LINC) bridging complex at the nuclear envelope (NE). Mutations in nesprin-1/2 have previously been found in patients with autosomal dominant Emery–Dreifuss muscular dystrophy (EDMD) as well as dilated cardiomyopathy (DCM). In this study, three novel rare variants (R8272Q, S8381C and N8406K) in the C-terminus of the SYNE1 gene (nesprin-1) were identified in seven DCM patients by mutation screening. Expression of these mutants caused nuclear morphology defects and reduced lamin A/C and SUN2 staining at the NE. GST pull-down indicated that nesprin-1/lamin/SUN interactions were disrupted. Nesprin-1 mutations were also associated with augmented activation of the ERK pathway in vitro and in hearts in vivo. During C2C12 muscle cell differentiation, nesprin-1 levels are increased concomitantly with kinesin light chain (KLC-1/2) and immunoprecipitation and GST pull-down showed that these proteins interacted via a recently identified LEWD domain in the C-terminus of nesprin-1. Expression of nesprin-1 mutants in C2C12 cells caused defects in myoblast differentiation and fusion associated with dysregulation of myogenic transcription factors and disruption of the nesprin-1 and KLC-1/2 interaction at the outer nuclear membrane. Expression of nesprin-1α₂ WT and mutants in zebrafish embryos caused heart developmental defects that varied in severity. These findings support a role for nesprin-1 in myogenesis and muscle disease, and uncover a novel mechanism whereby disruption of the LINC complex may contribute to the pathogenesis of DCM.

Experimental Modeling Supports a Role for MyBP-HL as a Novel Myofilament Component in Arrhythmia and Dilated Cardiomyopathy

BACKGROUND

Cardiomyopathy and arrhythmias are under significant genetic influence. Here, we studied a family with dilated cardiomyopathy and associated conduction system disease in whom prior clinical cardiac gene panel testing was unrevealing.

METHODS

Whole-genome sequencing and induced pluripotent stem cells were used to examine a family with dilated cardiomyopathy and atrial and ventricular arrhythmias. We also characterized a mouse model with heterozygous and homozygous deletion of Mybphl.

RESULTS

Whole-genome sequencing identified a premature stop codon, R255X, in the MYBPHL gene encoding MyBP-HL (myosin-binding protein-H like), a novel member of the myosin-binding protein family. MYBPHL was found to have high atrial expression with low ventricular expression. We determined that MyBP-HL protein was myofilament associated in the atria, and truncated MyBP-HL protein failed to incorporate into the myofilament. Human cell modeling demonstrated reduced expression from the mutant MYBPHL allele. Echocardiography of Mybphl heterozygous and null mouse hearts exhibited a 36% reduction in fractional shortening and an increased diastolic ventricular chamber size. Atria weight normalized to total heart weight was significantly increased in Mybphl heterozygous and null mice. Using a reporter system, we detected robust expression of Mybphl in the atria, and in discrete puncta throughout the right ventricular wall and septum, as well. Telemetric electrocardiogram recordings in Mybphl mice revealed cardiac conduction system abnormalities with aberrant atrioventricular conduction and an increased rate of arrhythmia in heterozygous and null mice.

CONCLUSIONS

The findings of reduced ventricular function and conduction system defects in Mybphl mice support that MYBPHL truncations may increase risk for human arrhythmias and cardiomyopathy.

Abstract 74: Whole Genome Sequencing as a Diagnostic Tool for Cardiomyopathy

Background: Congestive heart failure (CHF) is an increasing medical problem that is disabling and costly, and CHF currently affects more than 5 million Americans. A leading cause of CHF is cardiomyopathy, a disorder with a high heritable component. Genetic studies of familial cardiomyopathy have identified more than 50 responsible genes. Knowledge of the genetic underpinnings of disease provides prognostic information and can guide clinical decision- making. Clinical genetic testing employs gene panels for diagnosis, and typical panels include 20-100 genes. We assessed the viability of whole genome sequencing (WGS) to identify potentially pathogenic cardiomyopathy-associated genetic variants in established and novel genes.

Methods: WGS was conducted on 91 families (99 individuals) with cardiomyopathy. In 46/91 (50%) families, there was a prior history of unrevealing panel-based clinical genetic testing. WGS was performed using Illumina next generation sequencing. Variants were called and annotated using MegaSeq and MegaSeq2. MegaSeq2 refines the number of candidate variants from ~4 million to <30 using variant effect prediction, conservation at variant site, frequency, and expression data.

Results: MegaSeq2 was applied to 99 cardiomyopathy genomes and in 45/91 families (49%), a potentially pathogenic or likely pathogenic variant was identified. Here we describe identification and segregation analysis in 4 families. Two families have unique RBM20 variants that segregate with cardiomyopathy. One family was found to carry a premature stop in SCN5A and exhibits variable phenotypes including sudden death, arrhythmia and cardiomyopathy. The fourth family has cardiomyopathy that segregates with a BAG3 early termination variant.

Conclusions: WGS was able to identify potentially pathogenic variants in nearly half (45/91). While some of the variants identified were in novel genes, the majority of pathogenic variants were in established cardiomyopathy genes despite previous panel-based testing. More comprehensive analysis, like WGS, allows for re-evaluation of variation as additional genetic information becomes available.

Abstract 391: The Pathophysiological Role of MYBPHL in Dilated Cardiomyopathy

Background: Cardiomyopathy is a leading cause of heart failure and is highly heritable. One common form of cardiomyopathy is dilated cardiomyopathy (DCM), which currently has over 70 identified genes that have been described as causative for the disease. Genetic testing for DCM employs gene panels and has a sensitivity of mutation detection less than 50%, indicating that additional genes contribute to DCM. Here, we employed whole genome sequencing (WGS) in a family with DCM and heart block who had previously undergone unrevealing genetic testing. We identified a premature stop codon in the MYBPHL gene, a gene that has not previously been linked to DCM as a likely cause of DCM in this family. Myosin binding protein H Like (MyBP-HL) is a muscle-expressed protein bearing structural similarity to myosin binding protein C (MyBP-C), which is commonly mutated gene in cardiomyopathies.

Objective: Determine the physiological and pathophysiological role of Mybphl .

Results: RNA-seq and qPCR from mouse hearts revealed that Mybphl is highly expressed in the right and left atria with lower expression in the ventricle and virtually no expression in skeletal muscle. As MyBP-HL shares a high homology with the myofilament proteins cardiac myosin binding protein-C and H, we investigated if MyBP-HL is also myofilament-associated. We determined that MyBP-HL protein is myofilament-associated in the atria although not clearly so in ventricle. To assess the requirement of MyBP-HL in cardiac function, we used a mouse model with an insertional disruption of the Mybphl gene. These mice have deficits in in vivo cardiac function, with reduced fractional shortening. In addition, ECG recordings from the Mybphl null mice show conduction system abnormalities affecting atrioventricular conduction.

Conclusions: WGS identified a premature stop codon in MYBPHL in human DCM. A mouse model with a disrupted Mybphl gene showed similar pathophysiological features as the humans with reduced ventricular function and cardiac conduction system abnormalities. MyBP-HL is an important protein for normal cardiac function.

Genetic Variation in Cardiomyopathy and Cardiovascular Disorders

With the wider deployment of massively-parallel, next-generation sequencing, it is now possible to survey human genome data for research and clinical purposes. The reduced cost of producing short-read sequencing has now shifted the burden to data analysis. Analysis of genome sequencing remains challenged by the complexity of the human genome, including redundancy and the repetitive nature of genome elements and the large amount of variation in individual genomes. Public databases of human genome sequences greatly facilitate interpretation of common and rare genetic variation, although linking database sequence information to detailed clinical information is limited by privacy and practical issues. Genetic variation is a rich source of knowledge for cardiovascular disease because many, if not all, cardiovascular disorders are highly heritable. The role of rare genetic variation in predicting risk and complications of cardiovascular diseases has been well established for hypertrophic and dilated cardiomyopathy, where the number of genes that are linked to these disorders is growing. Bolstered by family data, where genetic variants segregate with disease, rare variation can be linked to specific genetic variation that offers profound diagnostic information. Understanding genetic variation in cardiomyopathy is likely to help stratify forms of heart failure and guide therapy. Ultimately, genetic variation may be amenable to gene correction and gene editing strategies.

Disruption of the lamin A and matrin-3 interaction by myopathic LMNA mutations

The nuclear face of the nuclear membrane is enriched with the intermediate filament protein lamin A. Mutations in LMNA, the gene encoding lamin A, lead to a diverse set of inherited conditions including myopathies that affect both the heart and skeletal muscle. To gain insight about lamin A protein interactions, binding proteins associated with the tail of lamin A were characterized. Of 130 nuclear proteins found associated with the lamin A tail, 17 (13%) were previously described lamin A binding partners. One protein not previously linked to lamin A, matrin-3, was selected for further study, because like LMNA mutations, matrin-3 has also been implicated in inherited myopathy. Matrin-3 binds RNA and DNA and is a nucleoplasmic protein originally identified from the insoluble nuclear fraction, referred to as the nuclear matrix. Anti-matrin-3 antibodies were found to co-immunoprecipitate lamin A, and the lamin-A binding domain was mapped to the carboxy-terminal half of matrin-3. Three-dimensional mapping of the lamin A-matrin-3 interface showed that the LMNA truncating mutation Δ303, which lacks the matrin-3 binding domain, was associated with an increased distance between lamin A and matrin-3. LMNA mutant cells are known to have altered biophysical properties and the matrin-3-lamin A interface is positioned to contribute to these defects.

The genetic landscape of cardiomyopathy and its role in heart failure

Heart failure is highly influenced by heritability, and nearly 100 genes link to familial cardiomyopathy. Despite the marked genetic diversity that underlies these complex cardiovascular phenotypes, several key genes and pathways have emerged. Hypertrophic cardiomyopathy is characterized by increased contractility and a greater energetic cost of cardiac output. Dilated cardiomyopathy is often triggered by mutations that disrupt the giant protein titin. The energetic consequences of these mutations offer molecular targets and opportunities for new drug development and gene correction therapies.

Targeted analysis of whole genome sequence data to diagnose genetic cardiomyopathy

BACKGROUND

Cardiomyopathy is highly heritable but genetically diverse. At present, genetic testing for cardiomyopathy uses targeted sequencing to simultaneously assess the coding regions of >50 genes. New genes are routinely added to panels to improve the diagnostic yield. With the anticipated $1000 genome, it is expected that genetic testing will shift toward comprehensive genome sequencing accompanied by targeted gene analysis. Therefore, we assessed the reliability of whole genome sequencing and targeted analysis to identify cardiomyopathy variants in 11 subjects with cardiomyopathy.

METHODS AND RESULTS

Whole genome sequencing with an average of 37× coverage was combined with targeted analysis focused on 204 genes linked to cardiomyopathy. Genetic variants were scored using multiple prediction algorithms combined with frequency data from public databases. This pipeline yielded 1 to 14 potentially pathogenic variants per individual. Variants were further analyzed using clinical criteria and segregation analysis, where available. Three of 3 previously identified primary mutations were detected by this analysis. In 6 subjects for whom the primary mutation was previously unknown, we identified mutations that segregated with disease, had clinical correlates, and had additional pathological correlation to provide evidence for causality. For 2 subjects with previously known primary mutations, we identified additional variants that may act as modifiers of disease severity. In total, we identified the likely pathological mutation in 9 of 11 (82%) subjects.

CONCLUSIONS

These pilot data demonstrate that ≈30 to 40× coverage whole genome sequencing combined with targeted analysis is feasible and sensitive to identify rare variants in cardiomyopathy-associated genes.

Genetic profiling for risk reduction in human cardiovascular disease

Cardiovascular disease is a major health concern affecting over 80,000,000 people in the U.S. alone. Heart failure, cardiomyopathy, heart rhythm disorders, atherosclerosis and aneurysm formation have significant heritable contribution. Supported by familial aggregation and twin studies, these cardiovascular diseases are influenced by genetic variation. Family-based linkage studies and population-based genome-wide association studies (GWAS) have each identified genes and variants important for the pathogenesis of cardiovascular disease. The advent of next generation sequencing has ushered in a new era in the genetic diagnosis of cardiovascular disease, and this is especially evident when considering cardiomyopathy, a leading cause of heart failure. Cardiomyopathy is a genetically heterogeneous disorder characterized by morphologically abnormal heart with abnormal function. Genetic testing for cardiomyopathy employs gene panels, and these panels assess more than 50 genes simultaneously. Despite the large size of these panels, the sensitivity for detecting the primary genetic defect is still only approximately 50%. Recently, there has been a shift towards applying broader exome and/or genome sequencing to interrogate more of the genome to provide a genetic diagnosis for cardiomyopathy. Genetic mutations in cardiomyopathy offer the capacity to predict clinical outcome, including arrhythmia risk, and genetic diagnosis often provides an early window in which to institute therapy. This discussion is an overview as to how genomic data is shaping the current understanding and treatment of cardiovascular disease.

Supercomputing for the parallelization of whole genome analysis

MOTIVATION

The declining cost of generating DNA sequence is promoting an increase in whole genome sequencing, especially as applied to the human genome. Whole genome analysis requires the alignment and comparison of raw sequence data, and results in a computational bottleneck because of limited ability to analyze multiple genomes simultaneously.

RESULTS

We now adapted a Cray XE6 supercomputer to achieve the parallelization required for concurrent multiple genome analysis. This approach not only markedly speeds computational time but also results in increased usable sequence per genome. Relying on publically available software, the Cray XE6 has the capacity to align and call variants on 240 whole genomes in ∼50 h. Multisample variant calling is also accelerated.

AVAILABILITY AND IMPLEMENTATION

The MegaSeq workflow is designed to harness the size and memory of the Cray XE6, housed at Argonne National Laboratory, for whole genome analysis in a platform designed to better match current and emerging sequencing volume.