Natural underlying mtDNA heteroplasmy as a potential source of intra-person hiPSC variability.

Functional variability among human clones of induced pluripotent stem cells (hiPSCs) remains a limitation in assembling high-quality biorepositories. Beyond inter-person variability, the root cause of intra-person variability remains unknown. Mitochondria guide the required transition from oxidative to glycolytic metabolism in nuclear reprogramming. Moreover, mitochondria have their own genome (mitochondrial DNA [mtDNA]). Herein, we performed mtDNA next-generation sequencing (NGS) on 84 hiPSC clones derived from a cohort of 19 individuals, including mitochondrial and non-mitochondrial patients. The analysis of mtDNA variants showed that low levels of potentially pathogenic mutations in the original fibroblasts are revealed through nuclear reprogramming, generating mutant hiPSCs with a detrimental effect in their differentiated progeny. Specifically, hiPSC-derived cardiomyocytes with expanded mtDNA mutations non-related with any described human disease, showed impaired mitochondrial respiration, being a potential cause of intra-person hiPSC variability. We propose mtDNA NGS as a new selection criterion to ensure hiPSC quality for drug discovery and regenerative medicine.

Recessive TAF1A mutations reveal ribosomopathy in siblings with end-stage pediatric dilated cardiomyopathy.

Non-ischemic dilated cardiomyopathy (DCM) has been recognized as a heritable disorder for over 25 years, yet clinical genetic testing is non-diagnostic in >50% of patients, underscoring the ongoing need for DCM gene discovery. Here, whole exome sequencing uncovered a novel molecular basis for idiopathic end-stage heart failure in two sisters who underwent cardiac transplantation at three years of age. Compound heterozygous recessive mutations in TAF1A, encoding an RNA polymerase I complex protein, were associated with marked fibrosis of explanted hearts and gene-specific nucleolar segregation defects in cardiomyocytes, indicative of impaired ribosomal RNA synthesis. Knockout of the homologous gene in zebrafish recapitulated a heart failure phenotype with pericardial edema, decreased ventricular systolic function, and embryonic mortality. These findings expand the clinical spectrum of ribosomopathies to include pediatric DCM.

Noncardiac genetic predisposition in sudden infant death syndrome.

PURPOSE: Sudden infant death syndrome (SIDS) is the commonest cause of sudden death of an infant; however, the genetic basis remains poorly understood. We aimed to identify noncardiac genes underpinning SIDS and determine their prevalence compared with ethnically matched controls. METHODS: Using exome sequencing we assessed the yield of ultrarare nonsynonymous variants (minor allele frequency [MAF] ≤0.00005, dominant model; MAF ≤0.01, recessive model) in 278 European SIDS cases (62% male; average age =2.7 ± 2 months) versus 973 European controls across 61 noncardiac SIDS-susceptibility genes. The variants were classified according to American College of Medical Genetics and Genomics criteria. Case-control, gene-collapsing analysis was performed in eight candidate biological pathways previously implicated in SIDS pathogenesis. RESULTS: Overall 43/278 SIDS cases harbored an ultrarare single-nucleotide variant compared with 114/973 controls (15.5 vs. 11.7%, p=0.10). Only 2/61 noncardiac genes were significantly overrepresented in cases compared with controls (ECE1, 3/278 [1%] vs. 1/973 [0.1%] p=0.036; SLC6A4, 2/278 [0.7%] vs. 1/973 [0.1%] p=0.049). There was no difference in yield of pathogenic or likely pathogenic variants between cases and controls (1/278 [0.36%] vs. 4/973 [0.41%]; p=1.0). Gene-collapsing analysis did not identify any specific biological pathways to be significantly associated with SIDS. CONCLUSIONS: A monogenic basis for SIDS amongst the previously implicated noncardiac genes and their encoded biological pathways is negligible.

A Comprehensive Approach to Sequence-oriented IsomiR annotation (CASMIR): demonstration with IsomiR profiling in colorectal neoplasia.

BACKGROUND: MicroRNA (miRNA) profiling is an important step in studying biological associations and identifying marker candidates. miRNA exists in isoforms, called isomiRs, which may exhibit distinct properties. With conventional profiling methods, limitations in assay and analysis platforms may compromise isomiR interrogation. RESULTS: We introduce a comprehensive approach to sequence-oriented isomiR annotation (CASMIR) to allow unbiased identification of global isomiRs from small RNA sequencing data. In this approach, small RNA reads are maintained as independent sequences instead of being summarized under miRNA names. IsomiR features are identified through step-wise local alignment against canonical forms and precursor sequences. Through customizing the reference database, CASMIR is applicable to isomiR annotation across species. To demonstrate its application, we investigated isomiR profiles in normal and neoplastic human colorectal epithelia. We also ran miRDeep2, a popular miRNA analysis algorithm to validate isomiRs annotated by CASMIR. With CASMIR, specific and biologically relevant isomiR patterns could be identified. We note that specific isomiRs are often more abundant than their canonical forms. We identify isomiRs that are commonly up-regulated in both colorectal cancer and advanced adenoma, and illustrate advantages in targeting isomiRs as potential biomarkers over canonical forms. CONCLUSIONS: Studying miRNAs at the isomiR level could reveal new insight into miRNA biology and inform assay design for specific isomiRs. CASMIR facilitates comprehensive annotation of isomiR features in small RNA sequencing data for isomiR profiling and differential expression analysis.

Exome-Wide Rare Variant Analyses in Sudden Infant Death Syndrome.

OBJECTIVE: To determine whether a monogenic basis explains sudden infant death syndrome (SIDS) using an exome-wide focus. STUDY DESIGN: A cohort of 427 unrelated cases of SIDS (257 male; average age = 2.7 ± 1.9 months) underwent whole-exome sequencing. Exome-wide rare variant analyses were carried out with 278 SIDS cases of European ancestry (173 male; average age = 2.7 ± 1.98 months) and 973 ethnic-matched controls based on 6 genetic models. Ingenuity Pathway Analysis also was performed. The cohort was collected in collaboration with coroners, medical examiners, and pathologists by St George's University of London, United Kingdom, and Mayo Clinic, Rochester, Minnesota. Whole-exome sequencing was performed at the Genomic Laboratory, Kings College London, United Kingdom, or Mayo Clinic's Medical Genome Facility, Rochester, Minnesota. RESULTS: Although no exome-wide significant (P < 2.5 × 10-6) difference in burden of ultra-rare variants was detected for any gene, 405 genes had a greater prevalence (P < .05) of ultra-rare nonsynonymous variants among cases with 17 genes at P < .005. Some of these potentially overrepresented genes may represent biologically plausible novel candidate genes for a monogenic basis for a portion of patients with SIDS. The top canonical pathway identified was glucocorticoid biosynthesis (P = .01). CONCLUSIONS: The lack of exome-wide significant genetic associations indicates an extreme heterogeneity of etiologies underlying SIDS. Our approach to understanding the genetic mechanisms of SIDS has far reaching implications for the SIDS research community as a whole and may catalyze new evidence-based SIDS research across multiple disciplines. Perturbations in glucocorticoid biosynthesis may represent a novel SIDS-associated biological pathway for future SIDS investigative research.

Enhancer of Zeste Homolog 2 Inhibition Stimulates Bone Formation and Mitigates Bone Loss Caused by Ovariectomy in Skeletally Mature Mice.

Perturbations in skeletal development and bone degeneration may result in reduced bone mass and quality, leading to greater fracture risk. Bone loss is mitigated by bone protective therapies, but there is a clinical need for new bone-anabolic agents. Previous work has demonstrated that Ezh2 (enhancer of zeste homolog 2), a histone 3 lysine 27 (H3K27) methyltransferase, suppressed differentiation of osteogenic progenitors. Here, we investigated whether inhibition of Ezh2 can be leveraged for bone stimulatory applications. Pharmacologic inhibition and siRNA knockdown of Ezh2 enhanced osteogenic commitment of MC3T3 preosteoblasts. Next generation RNA sequencing of mRNAs and real time quantitative PCR profiling established that Ezh2 inactivation promotes expression of bone-related gene regulators and extracellular matrix proteins. Mechanistically, enhanced gene expression was linked to decreased H3K27 trimethylation (H3K27me3) near transcriptional start sites in genome-wide sequencing of chromatin immunoprecipitations assays. Administration of an Ezh2 inhibitor modestly increases bone density parameters of adult mice. Furthermore, Ezh2 inhibition also alleviated bone loss in an estrogen-deficient mammalian model for osteoporosis. Ezh2 inhibition enhanced expression of Wnt10b and Pth1r and increased the BMP-dependent phosphorylation of Smad1/5. Thus, these data suggest that inhibition of Ezh2 promotes paracrine signaling in osteoblasts and has bone-anabolic and osteoprotective potential in adults.

Homozygous/Compound Heterozygous Triadin Mutations Associated With Autosomal-Recessive Long-QT Syndrome and Pediatric Sudden Cardiac Arrest: Elucidation of the Triadin Knockout Syndrome.

BACKGROUND: Long-QT syndrome (LQTS) may result in syncope, seizures, or sudden cardiac arrest. Although 16 LQTS-susceptibility genes have been discovered, 20% to 25% of LQTS remains genetically elusive. METHODS AND RESULTS: We performed whole-exome sequencing child-parent trio analysis followed by recessive and sporadic inheritance modeling and disease-network candidate analysis gene ranking to identify a novel underlying genetic mechanism for LQTS. Subsequent mutational analysis of the candidate gene was performed with polymerase chain reaction, denaturing high-performance liquid chromatography, and DNA sequencing on a cohort of 33 additional unrelated patients with genetically elusive LQTS. After whole-exome sequencing and variant filtration, a homozygous p.D18fs*13 TRDN-encoded triadin frameshift mutation was discovered in a 10-year-old female patient with LQTS with a QTc of 500 milliseconds who experienced recurrent exertion-induced syncope/cardiac arrest beginning at 1 year of age. Subsequent mutational analysis of TRDN revealed either homozygous or compound heterozygous frameshift mutations in 4 of 33 unrelated cases of LQTS (12%). All 5 TRDN-null patients displayed extensive T-wave inversions in precordial leads V1 through V4, with either persistent or transient QT prolongation and severe disease expression of exercise-induced cardiac arrest in early childhood (≤3 years of age) and required aggressive therapy. The overall yield of TRDN mutations was significantly greater in patients ≤10 years of age (5 of 10, 50%) compared with older patients (0 of 24, 0%; P=0.0009). CONCLUSIONS: We identified TRDN as a novel underlying genetic basis for recessively inherited LQTS. All TRDN-null patients had strikingly similar phenotypes. Given the recurrent nature of potential lethal arrhythmias, patients fitting this phenotypic profile should undergo cardiac TRDN genetic testing.

TNNI3K mutation in familial syndrome of conduction system disease, atrial tachyarrhythmia and dilated cardiomyopathy.

Locus mapping has uncovered diverse etiologies for familial atrial fibrillation (AF), dilated cardiomyopathy (DCM), and mixed cardiac phenotype syndromes, yet the molecular basis for these disorders remains idiopathic in most cases. Whole-exome sequencing (WES) provides a powerful new tool for familial disease gene discovery. Here, synergistic application of these genomic strategies identified the pathogenic mutation in a familial syndrome of atrial tachyarrhythmia, conduction system disease (CSD), and DCM vulnerability. Seven members of a three-generation family exhibited the variably expressed phenotype, three of whom manifested CSD and clinically significant arrhythmia in childhood. Genome-wide linkage analysis mapped two equally plausible loci to chromosomes 1p3 and 13q12. Variants from WES of two affected cousins were filtered for rare, predicted-deleterious, positional variants, revealing an unreported heterozygous missense mutation disrupting the highly conserved kinase domain in TNNI3K. The G526D substitution in troponin I interacting kinase, with the most deleterious SIFT and Polyphen2 scores possible, resulted in abnormal peptide aggregation in vitro and in silico docking models predicted altered yet energetically favorable wild-type mutant dimerization. Ventricular tissue from a mutation carrier displayed histopathological hallmarks of DCM and reduced TNNI3K protein staining with unique amorphous nuclear and sarcoplasmic inclusions. In conclusion, mutation of TNNI3K, encoding a heart-specific kinase previously shown to modulate cardiac conduction and myocardial function in mice, underlies a familial syndrome of electrical and myopathic heart disease. The identified substitution causes a TNNI3K aggregation defect and protein deficiency, implicating a dominant-negative loss of function disease mechanism.

De novo RRAGC mutation activates mTORC1 signaling in syndromic fetal dilated cardiomyopathy.

Idiopathic dilated cardiomyopathy (DCM) is a heritable, genetically heterogeneous disorder with variable age-dependent penetrance. We sought to identify the genetic underpinnings of syndromic, sporadic DCM in a newborn female diagnosed in utero. Postnatal evaluation revealed ventricular dilation and systolic dysfunction, bilateral cataracts, and mild facial dysmorphisms. Comprehensive metabolic and genetic testing, including chromosomal microarray, mitochondrial DNA and targeted RASopathy gene sequencing, and clinical whole exome sequencing for known cardiomyopathy genes was non-diagnostic. Following exclusion of asymptomatic DCM in the parents, trio-based whole exome sequencing was carried out on a research basis, filtering for rare, predicted deleterious de novo and recessive variants. An unreported de novo S75Y mutation was discovered in RRAGC, encoding Ras-related GTP binding C, an essential GTPase in nutrient-activated mechanistic target of rapamycin complex 1 (mTORC1) signaling. In silico protein modeling and molecular dynamics simulation predicted the mutation to disrupt ligand interactions and increase the GDP-bound state. Overexpression of RagC(S75Y) rendered AD293 cells partially insensitive to amino acid deprivation, resulting in increased mTORC1 signaling compared to wild-type RagC. These findings implicate mTORC1 dysregulation through a gain-of-function mutation in RagC as a novel molecular basis for syndromic forms of pediatric heart failure, and expand genotype-phenotype correlation in RASopathy-related syndromes.

Compound heterozygous NOTCH1 mutations underlie impaired cardiogenesis in a patient with hypoplastic left heart syndrome.

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart defect (CHD) that necessitates staged, single ventricle surgical palliation. An increased frequency of bicuspid aortic valve (BAV) has been observed among relatives. We postulated number of mutant alleles as a molecular basis for variable CHD expression in an extended family comprised of an HLHS proband and four family members who underwent echocardiography and whole-genome sequencing (WGS). Dermal fibroblast-derived induced pluripotent stem cells (iPSC) were procured from the proband-parent trio and bioengineered into cardiomyocytes. Cardiac phenotyping revealed aortic valve atresia and a slit-like left ventricular cavity in the HLHS proband, isolated bicuspid pulmonary valve in his mother, BAV in a maternal 4° relative, and no CHD in his father or sister. Filtering of WGS for rare, functional variants that segregated with CHD and were compound heterozygous in the HLHS proband identified NOTCH1 as the sole candidate gene. An unreported missense mutation (P1964L) in the cytoplasmic domain, segregating with semilunar valve malformation, was maternally inherited and a rare missense mutation (P1256L) in the extracellular domain, clinically silent in the heterozygous state, was paternally inherited. Patient-specific iPSCs exhibited diminished transcript levels of NOTCH1 signaling pathway components, impaired myocardiogenesis, and a higher prevalence of heterogeneous myofilament organization. Extended, phenotypically characterized families enable WGS-derived variant filtering for plausible Mendelian modes of inheritance, a powerful strategy to discover molecular underpinnings of CHD. Identification of compound heterozygous NOTCH1 mutations and iPSC-based functional modeling implicate mutant allele burden and impaired myogenic potential as mechanisms for HLHS.

PatternCNV: a versatile tool for detecting copy number changes from exome sequencing data.

MOTIVATION: Exome sequencing (exome-seq) data, which are typically used for calling exonic mutations, have also been utilized in detecting DNA copy number variations (CNVs). Despite the existence of several CNV detection tools, there is still a great need for a sensitive and an accurate CNV-calling algorithm with built-in QC steps, and does not require a paired reference for each sample. RESULTS: We developed a novel method named PatternCNV, which (i) accounts for the read coverage variations between exons while leveraging the consistencies of this variability across different samples; (ii) reduces alignment BAM files to WIG format and therefore greatly accelerates computation; (iii) incorporates multiple QC measures designed to identify outlier samples and batch effects; and (iv) provides a variety of visualization options including chromosome, gene and exon-level views of CNVs, along with a tabular summarization of the exon-level CNVs. Compared with other CNV-calling algorithms using data from a lymphoma exome-seq study, PatternCNV has higher sensitivity and specificity. AVAILABILITY AND IMPLEMENTATION: The software for PatternCNV is implemented using Perl and R, and can be used in Mac or Linux environments. Software and user manual are available at http://bioinformaticstools.mayo.edu/research/patterncnv/, and R package at https://github.com/topsoil/patternCNV/.

Exome sequencing establishes diagnosis of Alström syndrome in an infant presenting with non-syndromic dilated cardiomyopathy.

Idiopathic dilated cardiomyopathy is a heritable, genetically heterogeneous disorder characterized by progressive heart failure. Dilated cardiomyopathy typically exhibits autosomal dominant inheritance, yet frequently remains clinically silent until adulthood. We sought to discover the molecular basis of idiopathic, non-syndromic dilated cardiomyopathy in a one-month-old male presenting with severe heart failure. Previous comprehensive testing of blood, urine, and skin biopsy specimen was negative for metabolic, mitochondrial, storage, and infectious etiologies. Ophthalmologic examination was normal. Chromosomal microarray and commercial dilated cardiomyopathy gene panel testing failed to identify a causative mutation. Parental screening echocardiograms revealed no evidence of clinically silent dilated cardiomyopathy. Whole exome sequencing was carried out on the family trio on a research basis, filtering for rare, deleterious, recessive and de novo genetic variants. Pathogenic compound heterozygous truncating mutations were identified in ALMS1, diagnostic of Alström syndrome and prompting disclosure of genetic findings. Alström syndrome is a known cause for dilated cardiomyopathy in children yet delayed and mis-diagnosis are common owing to its rarity and age-dependent emergence of multisystem clinical manifestations. At six months of age the patient ultimately developed bilateral nystagmus and hyperopia, features characteristic of the syndrome. Early diagnosis is guiding clinical monitoring of other organ systems and allowing for presymptomatic intervention. Furthermore, recognition of recessive inheritance as the mechanism for sporadic disease has informed family planning. This case highlights a limitation of standard gene testing panels for pediatric dilated cardiomyopathy and exemplifies the potential for whole exome sequencing to solve a diagnostic dilemma and enable personalized care.

Histone deacetylase inhibition promotes osteoblast maturation by altering the histone H4 epigenome and reduces Akt phosphorylation.

Bone has remarkable regenerative capacity, but this ability diminishes during aging. Histone deacetylase inhibitors (HDIs) promote terminal osteoblast differentiation and extracellular matrix production in culture. The epigenetic events altered by HDIs in osteoblasts may hold clues for the development of new anabolic treatments for osteoporosis and other conditions of low bone mass. To assess how HDIs affect the epigenome of committed osteoblasts, MC3T3 cells were treated with suberoylanilide hydroxamic acid (SAHA) and subjected to microarray gene expression profiling and high-throughput ChIP-Seq analysis. As expected, SAHA induced differentiation and matrix calcification of osteoblasts in vitro. ChIP-Seq analysis revealed that SAHA increased histone H4 acetylation genome-wide and in differentially regulated genes, except for the 500 bp upstream of transcriptional start sites. Pathway analysis indicated that SAHA increased the expression of insulin signaling modulators, including Slc9a3r1. SAHA decreased phosphorylation of insulin receptor ß, Akt, and the Akt substrate FoxO1, resulting in FoxO1 stabilization. Thus, SAHA induces genome-wide H4 acetylation and modulates the insulin/Akt/FoxO1 signaling axis, whereas it promotes terminal osteoblast differentiation in vitro.

High-resolution molecular validation of self-renewal and spontaneous differentiation in clinical-grade adipose-tissue derived human mesenchymal stem cells.

Improving the effectiveness of adipose-tissue derived human mesenchymal stromal/stem cells (AMSCs) for skeletal therapies requires a detailed characterization of mechanisms supporting cell proliferation and multi-potency. We investigated the molecular phenotype of AMSCs that were either actively proliferating in platelet lysate or in a basal non-proliferative state. Flow cytometry combined with high-throughput RNA sequencing (RNASeq) and RT-qPCR analyses validate that AMSCs express classic mesenchymal cell surface markers (e.g., CD44, CD73/NT5E, CD90/THY1, and CD105/ENG). Expression of CD90 is selectively elevated at confluence. Self-renewing AMSCs express a standard cell cycle program that successively mediates DNA replication, chromatin packaging, cyto-architectural enlargement, and mitotic division. Confluent AMSCs preferentially express genes involved in extracellular matrix (ECM) formation and cellular communication. For example, cell cycle-related biomarkers (e.g., cyclins E2 and B2, transcription factor E2F1) and histone-related genes (e.g., H4, HINFP, NPAT) are elevated in proliferating AMSCs, while ECM genes are strongly upregulated (>10-fold) in quiescent AMSCs. AMSCs also express pluripotency genes (e.g., POU5F1, NANOG, KLF4) and early mesenchymal markers (e.g., NES, ACTA2) consistent with their multipotent phenotype. Strikingly, AMSCs modulate expression of WNT signaling components and switch production of WNT ligands (from WNT5A/WNT5B/WNT7B to WNT2/WNT2B), while upregulating WNT-related genes (WISP2, SFRP2, and SFRP4). Furthermore, post-proliferative AMSCs spontaneously express fibroblastic, osteogenic, chondrogenic, and adipogenic biomarkers when maintained in confluent cultures. Our findings validate the biological properties of self-renewing and multi-potent AMSCs by providing high-resolution quality control data that support their clinical versatility.

Identification of Rare Genetic Variants in Familial Spontaneous Coronary Artery Dissection and Evidence for Shared Biological Pathways.

Rare familial spontaneous coronary artery dissection (SCAD) kindreds implicate genetic disease predisposition and provide a unique opportunity for candidate gene discovery. Whole-genome sequencing was performed in fifteen probands with non-syndromic SCAD who had a relative with SCAD, eight of whom had a second relative with extra-coronary arteriopathy. Co-segregating variants and associated genes were prioritized by quantitative variant, gene, and disease-level metrics. Curated public databases were queried for functional relationships among encoded proteins. Fifty-four heterozygous coding variants in thirteen families co-segregated with disease and fulfilled primary filters of rarity, gene variation constraint, and predicted-deleterious protein effect. Secondary filters yielded 11 prioritized candidate genes in 12 families, with high arterial tissue expression (n = 7), high-confidence protein-level interactions with genes associated with SCAD previously (n = 10), and/or previous associations with connective tissue disorders and aortopathies (n = 3) or other vascular phenotypes in mice or humans (n = 11). High-confidence associations were identified among 10 familial SCAD candidate-gene-encoded proteins. A collagen-encoding gene was identified in five families, two with distinct variants in COL4A2. Familial SCAD is genetically heterogeneous, yet perturbations of extracellular matrix, cytoskeletal, and cell-cell adhesion proteins implicate common disease-susceptibility pathways. Incomplete penetrance and variable expression suggest genetic or environmental modifiers.

The Development of an Infrastructure to Facilitate the Use of Whole Genome Sequencing for Population Health.

The clinical use of genomic analysis has expanded rapidly resulting in an increased availability and utility of genomic information in clinical care. We have developed an infrastructure utilizing informatics tools and clinical processes to facilitate the use of whole genome sequencing data for population health management across the healthcare system. Our resulting framework scaled well to multiple clinical domains in both pediatric and adult care, although there were domain specific challenges that arose. Our infrastructure was complementary to existing clinical processes and well-received by care providers and patients. Informatics solutions were critical to the successful deployment and scaling of this program. Implementation of genomics at the scale of population health utilizes complicated technologies and processes that for many health systems are not supported by current information systems or in existing clinical workflows. To scale such a system requires a substantial clinical framework backed by informatics tools to facilitate the flow and management of data. Our work represents an early model that has been successful in scaling to 29 different genes with associated genetic conditions in four clinical domains. Work is ongoing to optimize informatics tools; and to identify best practices for translation to smaller healthcare systems.

Genetic Association Between Hypoplastic Left Heart Syndrome and Cardiomyopathies.

BACKGROUND: Hypoplastic left heart syndrome (HLHS) with risk of poor outcome has been linked to MYH6 variants, implicating overlap in genetic etiologies of structural and myopathic heart disease. METHODS: Whole genome sequencing was performed in 197 probands with HLHS, 43 family members, and 813 controls. Data were filtered for rare, segregating variants in 3 index families comprised of an HLHS proband and relative(s) with cardiomyopathy. Whole genome sequencing data from cases and controls were compared for rare variant burden across 56 cardiomyopathy genes utilizing a weighted burden test approach, accounting for multiple testing using a Bonferroni correction. RESULTS: A pathogenic MYBPC3 nonsense variant was identified in the first proband who underwent cardiac transplantation for diastolic heart failure, her father with left ventricular noncompaction, and 2 fourth-degree relatives with hypertrophic cardiomyopathy. A likely pathogenic RYR2 missense variant was identified in the second proband, a second-degree relative with aortic dilation, and a fourth-degree relative with dilated cardiomyopathy. A pathogenic RYR2 exon 3 in-frame deletion was identified in the third proband diagnosed with catecholaminergic polymorphic ventricular tachycardia and his father with left ventricular noncompaction and catecholaminergic polymorphic ventricular tachycardia. To further investigate HLHS-cardiomyopathy gene associations in cases versus controls, rare variant burden testing of 56 genes revealed enrichment in MYH6 (P=0.000068). Rare, predicted-damaging MYH6 variants were identified in 10% of probands in our cohort-4 with familial congenital heart disease, 4 with compound heterozygosity (3 with systolic ventricular dysfunction), and 4 with MYH6-FLNC synergistic heterozygosity. CONCLUSIONS: Whole genome sequencing in multiplex families, proband-parent trios, and case-control cohorts revealed defects in cardiomyopathy-associated genes in patients with HLHS, which may portend impaired functional reserve of the single-ventricle circulation.

Mitochondrial DNA: Hotspot for Potential Gene Modifiers Regulating Hypertrophic Cardiomyopathy.

Hypertrophic cardiomyopathy (HCM) is a prevalent and untreatable cardiovascular disease with a highly complex clinical and genetic causation. HCM patients bearing similar sarcomeric mutations display variable clinical outcomes, implying the involvement of gene modifiers that regulate disease progression. As individuals exhibiting mutations in mitochondrial DNA (mtDNA) present cardiac phenotypes, the mitochondrial genome is a promising candidate to harbor gene modifiers of HCM. Herein, we sequenced the mtDNA of isogenic pluripotent stem cell-cardiomyocyte models of HCM focusing on two sarcomeric mutations. This approach was extended to unrelated patient families totaling 52 cell lines. By correlating cellular and clinical phenotypes with mtDNA sequencing, potentially HCM-protective or -aggravator mtDNA variants were identified. These novel mutations were mostly located in the non-coding control region of the mtDNA and did not overlap with those of other mitochondrial diseases. Analysis of unrelated patients highlighted family-specific mtDNA variants, while others were common in particular population haplogroups. Further validation of mtDNA variants as gene modifiers is warranted but limited by the technically challenging methods of editing the mitochondrial genome. Future molecular characterization of these mtDNA variants in the context of HCM may identify novel treatments and facilitate genetic screening in cardiomyopathy patients towards more efficient treatment options.

Patient-specific genomics and cross-species functional analysis implicate LRP2 in hypoplastic left heart syndrome.

Congenital heart diseases (CHDs), including hypoplastic left heart syndrome (HLHS), are genetically complex and poorly understood. Here, a multidisciplinary platform was established to functionally evaluate novel CHD gene candidates, based on whole-genome and iPSC RNA sequencing of a HLHS family-trio. Filtering for rare variants and altered expression in proband iPSCs prioritized 10 candidates. siRNA/RNAi-mediated knockdown in healthy human iPSC-derived cardiomyocytes (hiPSC-CM) and in developing Drosophila and zebrafish hearts revealed that LDL receptor-related protein LRP2 is required for cardiomyocyte proliferation and differentiation. Consistent with hypoplastic heart defects, compared to patents the proband's iPSC-CMs exhibited reduced proliferation. Interestingly, rare, predicted-damaging LRP2 variants were enriched in a HLHS cohort; however, understanding their contribution to HLHS requires further investigation. Collectively, we have established a multi-species high-throughput platform to rapidly evaluate candidate genes and their interactions during heart development, which are crucial first steps toward deciphering oligogenic underpinnings of CHDs, including hypoplastic left hearts.

Rare Missense Variants in TLN1 Are Associated With Familial and Sporadic Spontaneous Coronary Artery Dissection.

BACKGROUND: Spontaneous coronary artery dissection (SCAD) is an uncommon idiopathic disorder predominantly affecting young, otherwise healthy women. Rare familial cases reveal a genetic predisposition to disease. The aim of this study was to identify a novel susceptibility gene for SCAD. METHODS: Whole-exome sequencing was performed in a family comprised of 3 affected individuals and filtered to identify rare, predicted deleterious, segregating variants. Immunohistochemical staining was used to evaluate protein expression of the identified candidate gene. The prevalence and spectrum of rare (<0.1%) variants within binding domains was determined by next-generation sequencing or denaturing high-performance liquid chromatography in a sporadic SCAD cohort of 675 unrelated individuals. RESULTS: We identified a rare heterozygous missense variant within a highly conserved ß-integrin-binding domain of TLN1 segregating with familial SCAD. TLN1 encodes talin 1-a large cytoplasmic protein of the integrin adhesion complex that links the actin cytoskeleton and extracellular matrix. Consistent with high mRNA expression in arterial tissues, robust immunohistochemical staining of talin 1 was demonstrated in coronary arteries. Nine additional rare heterozygous missense variants in TLN1 were identified in 10 sporadic cases. Incomplete penetrance, suggesting genetic or environmental modifiers of this episodic disorder, was evident in the familial case and 5 individuals with sporadic SCAD from whom parental DNA was available. CONCLUSIONS: Our findings reveal TLN1 as a disease-associated gene in familial and sporadic SCAD and, together with abnormal vascular phenotypes reported in animal models of talin 1 disruption, implicate impaired structural integrity of the coronary artery cytoskeleton in SCAD susceptibility.