Genetic Contribution to End-Stage Cardiomyopathy Requiring Heart Transplantation.

BACKGROUND: Many cardiovascular disorders propel the development of advanced heart failure that necessitates cardiac transplantation. When treatable causes are excluded, studies to define causes are often abandoned, resulting in a diagnosis of end-stage idiopathic cardiomyopathy. We studied whether DNA sequence analyses could identify unrecognized causes of end-stage nonischemic cardiomyopathy requiring heart transplantation and whether the prevalence of genetic causes differed from ambulatory cardiomyopathy cases. METHODS: We performed whole exome and genome sequencing of 122 explanted hearts from 101 adult and 21 pediatric patients with idiopathic cardiomyopathy from a single center. Data were analyzed for pathogenic/likely pathogenic variants in nuclear and mitochondrial genomes and assessed for nonhuman microbial sequences. The frequency of damaging genetic variants was compared among cardiomyopathy cohorts with different clinical severity. RESULTS: Fifty-four samples (44.3%) had pathogenic/likely pathogenic cardiomyopathy gene variants. The frequency of pathogenic variants was similar in pediatric (42.9%) and adult (43.6%) samples, but the distribution of mutated genes differed (P=8.30×10-4). The prevalence of causal genetic variants was significantly higher in end-stage than in previously reported ambulatory adult dilated cardiomyopathy cases (P<0.001). Among remaining samples with unexplained causes, no damaging mitochondrial variants were identified, but 28 samples contained parvovirus genome sequences, including 2 samples with 6- to 9-fold higher levels than the overall mean levels in other samples. CONCLUSIONS: Pathogenic variants and viral myocarditis were identified in 45.9% of patients with unexplained end-stage cardiomyopathy. Damaging gene variants are significantly more frequent among transplant compared with patients with ambulatory cardiomyopathy. Genetic analyses can help define cause of end-stage cardiomyopathy to guide management and risk stratification of patients and family members.

Contribution of Previously Unrecognized RNA Splice-Altering Variants to Congenital Heart Disease.

BACKGROUND: Known genetic causes of congenital heart disease (CHD) explain <40% of CHD cases, and interpreting the clinical significance of variants with uncertain functional impact remains challenging. We aim to improve diagnostic classification of variants in patients with CHD by assessing the impact of noncanonical splice region variants on RNA splicing. METHODS: We tested de novo variants from trio studies of 2649 CHD probands and their parents, as well as rare (allele frequency, <2×10-6) variants from 4472 CHD probands in the Pediatric Cardiac Genetics Consortium through a combined computational and in vitro approach. RESULTS: We identified 53 de novo and 74 rare variants in CHD cases that alter splicing and thus are loss of function. Of these, 77 variants are in known dominant, recessive, and candidate CHD genes, including KMT2D and RBFOX2. In 1 case, we confirmed the variant's predicted impact on RNA splicing in RNA transcripts from the proband's cardiac tissue. Two probands were found to have 2 loss-of-function variants for recessive CHD genes HECTD1 and DYNC2H1. In addition, SpliceAI-a predictive algorithm for altered RNA splicing-has a positive predictive value of ≈93% in our cohort. CONCLUSIONS: Through assessment of RNA splicing, we identified a new loss-of-function variant within a CHD gene in 78 probands, of whom 69 (1.5%; n=4472) did not have a previously established genetic explanation for CHD. Identification of splice-altering variants improves diagnostic classification and genetic diagnoses for CHD. REGISTRATION: URL: https://clinicaltrials.gov; Unique identifier: NCT01196182.

An ancient founder mutation located between ROBO1 and ROBO2 is responsible for increased microtia risk in Amerindigenous populations.

Microtia is a congenital malformation that encompasses mild hypoplasia to complete loss of the external ear, or pinna. Although the contribution of genetic variation and environmental factors to microtia remains elusive, Amerindigenous populations have the highest reported incidence. Here, using both transmission disequilibrium tests and association studies in microtia trios (parents and affected child) and microtia cohorts enrolled in Latin America, we map an ∼10-kb microtia locus (odds ratio = 4.7; P = 6.78e-18) to the intergenic region between Roundabout 1 (ROBO1) and Roundabout 2 (ROBO2) (chr3: 78546526 to 78555137). While alleles at the microtia locus significantly increase the risk of microtia, their penetrance is low (<1%). We demonstrate that the microtia locus contains a polymorphic complex repeat element that is expanded in affected individuals. The locus is located near a chromatin loop region that regulates ROBO1 and ROBO2 expression in induced pluripotent stem cell­derived neural crest cells. Furthermore, we use single nuclear RNA sequencing to demonstrate ROBO1 and ROBO2 expression in both fibroblasts and chondrocytes of the mature human pinna. Because the microtia allele is enriched in Amerindigenous populations and is shared by some East Asian subjects with craniofacial malformations, we propose that both populations share a mutation that arose in a common ancestor prior to the ancient migration of Eurasian populations into the Americas and that the high incidence of microtia among Amerindigenous populations reflects the population bottleneck that occurred during the migration out of Eurasia.

Neither cardiac mitochondrial DNA variation nor copy number contribute to congenital heart disease risk.

The well-established manifestation of mitochondrial mutations in functional cardiac disease (e.g., mitochondrial cardiomyopathy) prompted the hypothesis that mitochondrial DNA (mtDNA) sequence and/or copy number (mtDNAcn) variation contribute to cardiac defects in congenital heart disease (CHD). MtDNAcns were calculated and rare, non-synonymous mtDNA mutations were identified in 1,837 CHD-affected proband-parent trios, 116 CHD-affected singletons, and 114 paired cardiovascular tissue/blood samples. The variant allele fraction (VAF) of heteroplasmic variants in mitochondrial RNA from 257 CHD cardiovascular tissue samples was also calculated. On average, mtDNA from blood had 0.14 rare variants and 52.9 mtDNA copies per nuclear genome per proband. No variation with parental age at proband birth or CHD-affected proband age was seen. mtDNAcns in valve/vessel tissue (320 ± 70) were lower than in atrial tissue (1,080 ± 320, p = 6.8E-21), which were lower than in ventricle tissue (1,340 ± 280, p = 1.4E-4). The frequency of rare variants in CHD-affected individual DNA was indistinguishable from the frequency in an unaffected cohort, and proband mtDNAcns did not vary from those of CHD cohort parents. In both the CHD and the comparison cohorts, mtDNAcns were significantly correlated between mother-child, father-child, and mother-father. mtDNAcns among people with European (mean = 52.0), African (53.0), and Asian haplogroups (53.5) were calculated and were significantly different for European and Asian haplogroups (p = 2.6E-3). Variant heteroplasmic fraction (HF) in blood correlated well with paired cardiovascular tissue HF (r = 0.975) and RNA VAF (r = 0.953), which suggests blood HF is a reasonable proxy for HF in heart tissue. We conclude that mtDNA mutations and mtDNAcns are unlikely to contribute significantly to CHD risk.

Contribution of Noncanonical Splice Variants to TTN Truncating Variant Cardiomyopathy.

BACKGROUND: Heterozygous TTN truncating variants cause 10% to 20% of idiopathic dilated cardiomyopathy (DCM). Although variants which disrupt canonical splice signals (ie, invariant dinucleotide of the splice donor site, invariant dinucleotide of the splice acceptor site) at exon-intron junctions are readily recognized as TTN truncating variants, the effects of other nearby sequence variations on splicing and their contribution to disease is uncertain. METHODS: Rare variants of unknown significance located in the splice regions of highly expressed TTN exons from 203 DCM cases, 3329 normal subjects, and clinical variant databases were identified. The effects of these variants on splicing were assessed using an in vitro splice assay. RESULTS: Splice-altering variants of unknown significance were enriched in DCM cases over controls and present in 2% of DCM patients (P=0.002). Application of this method to clinical variant databases demonstrated 20% of similar variants of unknown significance in TTN splice regions affect splicing. Noncanonical splice-altering variants were most frequently located at position +5 of the donor site (P=4.4×107) and position -3 of the acceptor site (P=0.002). SpliceAI, an emerging in silico prediction tool, had a high positive predictive value (86%-95%) but poor sensitivity (15%-50%) for the detection of splice-altering variants. Alternate exons spliced out of most TTN transcripts frequently lacked the consensus base at +5 donor and -3 acceptor positions. CONCLUSIONS: Noncanonical splice-altering variants in TTN explain 1-2% of DCM and offer a 10-20% increase in the diagnostic power of TTN sequencing in this disease. These data suggest rules that may improve efforts to detect splice-altering variants in other genes and may explain the low percent splicing observed for many alternate TTN exons.

Mechanisms of Congenital Heart Disease Caused by NAA15 Haploinsufficiency.

[Figure: see text].

EM-mosaic detects mosaic point mutations that contribute to congenital heart disease.

BACKGROUND: The contribution of somatic mosaicism, or genetic mutations arising after oocyte fertilization, to congenital heart disease (CHD) is not well understood. Further, the relationship between mosaicism in blood and cardiovascular tissue has not been determined. METHODS: We developed a new computational method, EM-mosaic (Expectation-Maximization-based detection of mosaicism), to analyze mosaicism in exome sequences derived primarily from blood DNA of 2530 CHD proband-parent trios. To optimize this method, we measured mosaic detection power as a function of sequencing depth. In parallel, we analyzed our cohort using MosaicHunter, a Bayesian genotyping algorithm-based mosaic detection tool, and compared the two methods. The accuracy of these mosaic variant detection algorithms was assessed using an independent resequencing method. We then applied both methods to detect mosaicism in cardiac tissue-derived exome sequences of 66 participants for which matched blood and heart tissue was available. RESULTS: EM-mosaic detected 326 mosaic mutations in blood and/or cardiac tissue DNA. Of the 309 detected in blood DNA, 85/97 (88%) tested were independently confirmed, while 7/17 (41%) candidates of 17 detected in cardiac tissue were confirmed. MosaicHunter detected an additional 64 mosaics, of which 23/46 (50%) among 58 candidates from blood and 4/6 (67%) of 6 candidates from cardiac tissue confirmed. Twenty-five mosaic variants altered CHD-risk genes, affecting 1% of our cohort. Of these 25, 22/22 candidates tested were confirmed. Variants predicted as damaging had higher variant allele fraction than benign variants, suggesting a role in CHD. The estimated true frequency of mosaic variants above 10% mosaicism was 0.14/person in blood and 0.21/person in cardiac tissue. Analysis of 66 individuals with matched cardiac tissue available revealed both tissue-specific and shared mosaicism, with shared mosaics generally having higher allele fraction. CONCLUSIONS: We estimate that ~ 1% of CHD probands have a mosaic variant detectable in blood that could contribute to cardiac malformations, particularly those damaging variants with relatively higher allele fraction. Although blood is a readily available DNA source, cardiac tissues analyzed contributed ~ 5% of somatic mosaic variants identified, indicating the value of tissue mosaicism analyses.

SarcTrack.

RATIONALE: Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) in combination with CRISPR/Cas9 genome editing provide unparalleled opportunities to study cardiac biology and disease. However, sarcomeres, the fundamental units of myocyte contraction, are immature and nonlinear in hiPSC-CMs, which technically challenge accurate functional interrogation of contractile parameters in beating cells. Furthermore, existing analysis methods are relatively low-throughput, indirectly assess contractility, or only assess well-aligned sarcomeres found in mature cardiac tissues. OBJECTIVE: We aimed to develop an analysis platform that directly, rapidly, and automatically tracks sarcomeres in beating cardiomyocytes. The platform should assess sarcomere content, contraction and relaxation parameters, and beat rate. METHODS AND RESULTS: We developed SarcTrack, a MatLab software that monitors fluorescently tagged sarcomeres in hiPSC-CMs. The algorithm determines sarcomere content, sarcomere length, and returns rates of sarcomere contraction and relaxation. By rapid measurement of hundreds of sarcomeres in each hiPSC-CM, SarcTrack provides large data sets for robust statistical analyses of multiple contractile parameters. We validated SarcTrack by analyzing drug-treated hiPSC-CMs, confirming the contractility effects of compounds that directly activate (CK-1827452) or inhibit (MYK-461) myosin molecules or indirectly alter contractility (verapamil and propranolol). SarcTrack analysis of hiPSC-CMs carrying a heterozygous truncation variant in the myosin-binding protein C ( MYBPC3) gene, which causes hypertrophic cardiomyopathy, recapitulated seminal disease phenotypes including cardiac hypercontractility and diminished relaxation, abnormalities that normalized with MYK-461 treatment. CONCLUSIONS: SarcTrack provides a direct and efficient method to quantitatively assess sarcomere function. By improving existing contractility analysis methods and overcoming technical challenges associated with functional evaluation of hiPSC-CMs, SarcTrack enhances translational prospects for sarcomere-regulating therapeutics and accelerates interrogation of human cardiac genetic variants.

Heterogeneous dynamics of ionic liquids: A four-point time correlation function approach.

Many ionic liquids show behavior similar to that of glassy systems, e.g., large and long-lasted deviations from Gaussian dynamics and clustering of "mobile" and "immobile" groups of ions. Herein a time-dependent four-point density correlation function-typically used to characterize glassy systems-is implemented for the ionic liquids, choline acetate, and 1-butyl-3-methylimidazolium acetate. Dynamic correlation beyond the first ionic solvation shell on the time scale of nanoseconds is found in the ionic liquids, revealing the cooperative nature of ion motions. The traditional solvent, acetonitrile, on the other hand, shows a much shorter length-scale that decays after a few picoseconds.

Correction: A molecular dynamics study of the ionic liquid, choline acetate.

Correction for 'A molecular dynamics study of the ionic liquid, choline acetate' by Jon A. L. Willcox et al., Phys. Chem. Chem. Phys., 2016, 18, 14850-14858.

Molecular Dynamics Study of Water Flow across Multiple Layers of Pristine, Oxidized, and Mixed Regions of Graphene Oxide.

A molecular dynamics graphene oxide model is used to shed light on commonly overlooked features of graphene oxide membranes. The model features both perpendicular and parallel water flow across multiple sheets of pristine and/or oxidized graphene to simulate "brick-and-mortar" microstructures. Additionally, regions of pristine/oxidized graphene overlap that have thus far been overlooked in the literature are explored. Differences in orientational and hydrogen-bonding features between adjacent layers of water in this mixed region are found to be even more prominent than differences between pristine and oxidized channels. This region also shows lateral water flow in equilibrium simulations and orthogonal flow in non-equilibrium simulations significantly greater than those in the oxidized region, suggesting it may play a non-negligible role in the mechanism of water flow across graphene oxide membranes.

A molecular dynamics study of the ionic liquid, choline acetate.

Structural and dynamic properties of the ionic liquid (IL) choline acetate are studied using molecular dynamics (MD) simulations. The hydroxyl group of choline shows significant hydrogen-bonding interactions with the oxygen atoms of acetate. Nearly all choline cations are found to form a hydrogen bond with acetate anions at 400 K, while about 67% of cations participate in hydrogen-bonding interactions at 600 K. At 400 K, subdiffusive and prominent non-Gaussian behavior persist for t > 10 ns. At 600 K, the usual diffusion regime is obtained after a few hundred ps of subdiffusive behavior. Analysis of reorientational motions of acetate ions, particularly those of their short axes, indicates a high degree of dynamic heterogeneity, in agreement with previous work on different IL systems.

Bicephalic amphiphile architecture affects antibacterial activity.

A series of cationic amphiphiles, each with an aromatic core, was prepared and investigated for antimicrobial properties. The synthesized amphiphiles in this study are bicephalic (double-headed) in that they each possess two trimethylammonium head groups and a single linear alkoxy tail. Minimum inhibitory and minimum bactericidal concentrations of these amphiphiles were in the low micromolar range. Antimicrobial activities are highly sensitive to the chain length of the hydrophobic region, and modestly reliant on the relative positioning of the head groups on the aromatic core. These trends were more pronounced in time kill assays, wherein longer chain compounds required significantly shorter times to completely kill bacteria. Microscopy suggested that the mode of cell death was lysis. Strong inhibition was observed with both biscationic compounds and monocationic comparisons against Gram-positive bacteria; only biscationic amphiphiles maintained good activity versus the Gram-negative bacteria tested. These observations provide direction for future antimicrobial structural investigations.