Increased endothelial sclerostin caused by elevated DSCAM mediates multiple trisomy 21 phenotypes.

Trisomy 21 (T21), a recurrent aneuploidy occurring in 1:800 births, predisposes to congenital heart disease (CHD) and multiple extracardiac phenotypes. Despite a definitive genetic etiology, the mechanisms by which T21 perturbs development and homeostasis remain poorly understood. We compared the transcriptome of CHD tissues from 49 patients with T21 and 226 with euploid CHD (eCHD). We resolved cell lineages that misexpressed T21 transcripts by cardiac single-nucleus RNA sequencing and RNA in situ hybridization. Compared with eCHD samples, T21 samples had increased chr21 gene expression; 11-fold-greater levels (P = 1.2 × 10-8) of SOST (chr17), encoding the Wnt inhibitor sclerostin; and 1.4-fold-higher levels (P = 8.7 × 10-8) of the SOST transcriptional activator ZNF467 (chr7). Euploid and T21 cardiac endothelial cells coexpressed SOST and ZNF467; however, T21 endothelial cells expressed 6.9-fold more SOST than euploid endothelial cells (P = 2.7 × 10-27). Wnt pathway genes were downregulated in T21 endothelial cells. Expression of DSCAM, residing within the chr21 CHD critical region, correlated with SOST (P = 1.9 × 10-5) and ZNF467 (P = 2.9 × 10-4). Deletion of DSCAM from T21 endothelial cells derived from human induced pluripotent stem cells diminished sclerostin secretion. As Wnt signaling is critical for atrioventricular canal formation, bone health, and pulmonary vascular homeostasis, we concluded that T21-mediated increased sclerostin levels would inappropriately inhibit Wnt activities and promote Down syndrome phenotypes. These findings imply therapeutic potential for anti-sclerostin antibodies in T21.

Endocardial primary cilia and blood flow are required for regulation of EndoMT during endocardial cushion development.

Blood flow is critical for heart valve formation, and cellular mechanosensors are essential to translate flow into transcriptional regulation of development. Here, we identify a role for primary cilia in vivo in the spatial regulation of cushion formation, the first stage of valve development, by regionally controlling endothelial to mesenchymal transition (EndoMT) via modulation of Kruppel-like Factor 4 (Klf4) . We find that high shear stress intracardiac regions decrease endocardial ciliation over cushion development, correlating with KLF4 downregulation and EndoMT progression. Mouse embryos constitutively lacking cilia exhibit a blood-flow dependent accumulation of KLF4 in these regions, independent of upstream left-right abnormalities, resulting in impaired cushion cellularization. snRNA-seq revealed that cilia KO endocardium fails to progress to late-EndoMT, retains endothelial markers and has reduced EndoMT/mesenchymal genes that KLF4 antagonizes. Together, these data identify a mechanosensory role for endocardial primary cilia in cushion development through regional regulation of KLF4.

Functional dissection of human cardiac enhancers and noncoding de novo variants in congenital heart disease.

Rare coding mutations cause ∼45% of congenital heart disease (CHD). Noncoding mutations that perturb cis-regulatory elements (CREs) likely contribute to the remaining cases, but their identification has been problematic. Using a lentiviral massively parallel reporter assay (lentiMPRA) in human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs), we functionally evaluated 6,590 noncoding de novo variants (ncDNVs) prioritized from the whole-genome sequencing of 750 CHD trios. A total of 403 ncDNVs substantially affected cardiac CRE activity. A majority increased enhancer activity, often at regions with undetectable reference sequence activity. Of ten DNVs tested by introduction into their native genomic context, four altered the expression of neighboring genes and iPSC-CM transcriptional state. To prioritize future DNVs for functional testing, we used the MPRA data to develop a regression model, EpiCard. Analysis of an independent CHD cohort by EpiCard found enrichment of DNVs. Together, we developed a scalable system to measure the effect of ncDNVs on CRE activity and deployed it to systematically assess the contribution of ncDNVs to CHD.

Direct reprogramming of non-limb fibroblasts to cells with properties of limb progenitors.

The early limb bud consists of mesenchymal limb progenitors derived from the lateral plate mesoderm (LPM). The LPM also gives rise to the mesodermal components of the flank and neck. However, the cells at these other levels cannot produce the variety of cell types found in the limb. Taking advantage of a direct reprogramming approach, we find a set of factors (Prdm16, Zbtb16, and Lin28a) normally expressed in the early limb bud and capable of imparting limb progenitor-like properties to mouse non-limb fibroblasts. The reprogrammed cells show similar gene expression profiles and can differentiate into similar cell types as endogenous limb progenitors. The further addition of Lin41 potentiates the proliferation of the reprogrammed cells. These results suggest that these same four factors may play pivotal roles in the specification of endogenous limb progenitors.

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.

Photodynamic augmentation of oncolytic virus therapy for central nervous system malignancies.

Oncolytic viruses (OVs) have emerged as a clinical therapeutic modality potentially effective for cancers that evade conventional therapies, including central nervous system malignancies. Rationally designed combinatorial strategies can augment the efficacy of OVs by boosting tumor-selective cytotoxicity and modulating the tumor microenvironment (TME). Photodynamic therapy (PDT) of cancer not only mediates direct neoplastic cell death but also primes the TME to sensitize the tumor to secondary therapies, allowing for the combination of two potentially synergistic therapies with broader targets. Here, we created G47Δ-KR, clinical oncolytic herpes simplex virus G47Δ that expresses photosensitizer protein KillerRed (KR). Optical properties and cytotoxic effects of G47Δ-KR infection followed by amber LED illumination (peak wavelength: 585-595 nm) were examined in human glioblastoma (GBM) and malignant meningioma (MM) models in vitro. G47Δ-KR infection of tumor cells mediated KR expression that was activated by LED and produced reactive oxygen species, leading to cell death that was more robust than G47Δ-KR without light. In vivo, we tested photodynamic-oncolytic virus (PD-OV) therapy employing intratumoral injection of G47Δ-KR followed by laser light tumor irradiation (wavelength: 585 nm) in GBM and MM xenografts. PD-OV therapy was feasible in these models and resulted in potent anti-tumor effects that were superior to G47Δ-KR alone (without laser light) or laser light alone. RNA sequencing analysis of post-treatment tumor samples revealed PD-OV therapy-induced increases in TME infiltration of variable immune cell types. This study thus demonstrated the proof-of-concept that G47Δ-KR enables PD-OV therapy for neuro-oncological malignancies and warrants further research to advance potential clinical translation.

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.

Tbx5 maintains atrial identity by regulating an atrial enhancer network.

Understanding how the atrial and ventricular chambers of the heart maintain their distinct identity is a prerequisite for treating chamber-specific diseases. Here, we selectively inactivated the transcription factor Tbx5 in the atrial working myocardium of the neonatal mouse heart to show that it is required to maintain atrial identity. Atrial Tbx5 inactivation downregulated highly chamber specific genes such as Myl7 and Nppa , and conversely, increased the expression of ventricular identity genes including Myl2 . Using combined single nucleus transcriptome and open chromatin profiling, we assessed genomic accessibility changes underlying the altered atrial identity expression program, identifying 1846 genomic loci with greater accessibility in control atrial cardiomyocytes compared to KO aCMs. 69% of the control-enriched ATAC regions were bound by TBX5, demonstrating a role for TBX5 in maintaining atrial genomic accessibility. These regions were associated with genes that had higher expression in control aCMs compared to KO aCMs, suggesting they act as TBX5-dependent enhancers. We tested this hypothesis by analyzing enhancer chromatin looping using HiChIP and found 510 chromatin loops that were sensitive to TBX5 dosage. Of the loops enriched in control aCMs, 73.7% contained anchors in control-enriched ATAC regions. Together, these data demonstrate a genomic role for TBX5 in maintaining the atrial gene expression program by binding to atrial enhancers and preserving tissue-specific chromatin architecture of atrial enhancers.

Prevalence and Clinical Consequences of Multiple Pathogenic Variants in Dilated Cardiomyopathy.

BACKGROUND: Dilated cardiomyopathy (DCM) was considered a monogenetic disease that can be caused by over 60 genes. Evidence suggests that the combination of multiple pathogenic variants leads to greater disease severity and earlier onset. So far, not much is known about the prevalence and disease course of multiple pathogenic variants in patients with DCM. To gain insight into these knowledge gaps, we (1) systematically collected clinical information from a well-characterized DCM cohort and (2) created a mouse model. METHODS: Complete cardiac phenotyping and genotyping was performed in 685 patients with consecutive DCM. Compound heterozygous digenic (LMNA [lamin]/titin deletion A-band) with monogenic (LMNA/wild-type) and wild-type/wild-type mice were created and phenotypically followed over time. RESULTS: One hundred thirty-one likely pathogenic/pathogenic (LP/P) variants in robust DCM-associated genes were found in 685 patients with DCM (19.1%) genotyped for the robust genes. Three of the 131 patients had a second LP/P variant (2.3%). These 3 patients had a comparable disease onset, disease severity, and clinical course to patients with DCM with one LP/P. The LMNA/Titin deletion A-band mice had no functional differences compared with the LMNA/wild-type mice after 40 weeks of follow-up, although RNA-sequencing suggests increased cardiac stress and sarcomere insufficiency in the LMNA/Titin deletion A-band mice. CONCLUSIONS: In this study population, 2.3% of patients with DCM with one LP/P also have a second LP/P in a different gene. Although the second LP/P does not seem to influence the disease course of DCM in patients and mice, the finding of a second LP/P can be of importance to their relatives.

Cardiomyocyte infection by Trypanosoma cruzi promotes innate immune response and glycolysis activation.

Introduction: Chagas cardiomyopathy, a disease caused by Trypanosoma cruzi (T. cruzi) infection, is a major contributor to heart failure in Latin America. There are significant gaps in our understanding of the mechanism for infection of human cardiomyocytes, the pathways activated during the acute phase of the disease, and the molecular changes that lead to the progression of cardiomyopathy. Methods: To investigate the effects of T. cruzi on human cardiomyocytes during infection, we infected induced pluripotent stem cell-derived cardiomyocytes (iPSC-CM) with the parasite and analyzed cellular, molecular, and metabolic responses at 3 hours, 24 hours, and 48 hours post infection (hpi) using transcriptomics (RNAseq), proteomics (LC-MS), and metabolomics (GC-MS and Seahorse) analyses. Results: Analyses of multiomic data revealed that cardiomyocyte infection caused a rapid increase in genes and proteins related to activation innate and adaptive immune systems and pathways, including alpha and gamma interferons, HIF-1α signaling, and glycolysis. These responses resemble prototypic responses observed in pathogen-activated immune cells. Infection also caused an activation of glycolysis that was dependent on HIF-1α signaling. Using gene editing and pharmacological inhibitors, we found that T. cruzi uptake was mediated in part by the glucose-facilitated transporter GLUT4 and that the attenuation of glycolysis, HIF-1α activation, or GLUT4 expression decreased T. cruzi infection. In contrast, pre-activation of pro-inflammatory immune responses with LPS resulted in increased infection rates. Conclusion: These findings suggest that T. cruzi exploits a HIF-1α-dependent, cardiomyocyte-intrinsic stress-response activation of glycolysis to promote intracellular infection and replication. These chronic immuno-metabolic responses by cardiomyocytes promote dysfunction, cell death, and the emergence of cardiomyopathy.

Efficient in vivo genome editing prevents hypertrophic cardiomyopathy in mice.

Dominant missense pathogenic variants in cardiac myosin heavy chain cause hypertrophic cardiomyopathy (HCM), a currently incurable disorder that increases risk for stroke, heart failure and sudden cardiac death. In this study, we assessed two different genetic therapies-an adenine base editor (ABE8e) and a potent Cas9 nuclease delivered by AAV9-to prevent disease in mice carrying the heterozygous HCM pathogenic variant myosin R403Q. One dose of dual-AAV9 vectors, each carrying one half of RNA-guided ABE8e, corrected the pathogenic variant in ≥70% of ventricular cardiomyocytes and maintained durable, normal cardiac structure and function. An additional dose provided more editing in the atria but also increased bystander editing. AAV9 delivery of RNA-guided Cas9 nuclease effectively inactivated the pathogenic allele, albeit with dose-dependent toxicities, necessitating a narrow therapeutic window to maintain health. These preclinical studies demonstrate considerable potential for single-dose genetic therapies to correct or silence pathogenic variants and prevent the development of HCM.

Damaging variants in FOXI3 cause microtia and craniofacial microsomia.

PURPOSE: Craniofacial microsomia (CFM) represents a spectrum of craniofacial malformations, ranging from isolated microtia with or without aural atresia to underdevelopment of the mandible, maxilla, orbit, facial soft tissue, and/or facial nerve. The genetic causes of CFM remain largely unknown. METHODS: We performed genome sequencing and linkage analysis in patients and families with microtia and CFM of unknown genetic etiology. The functional consequences of damaging missense variants were evaluated through expression of wild-type and mutant proteins in vitro. RESULTS: We studied a 5-generation kindred with microtia, identifying a missense variant in FOXI3 (p.Arg236Trp) as the cause of disease (logarithm of the odds = 3.33). We subsequently identified 6 individuals from 3 additional kindreds with microtia-CFM spectrum phenotypes harboring damaging variants in FOXI3, a regulator of ectodermal and neural crest development. Missense variants in the nuclear localization sequence were identified in cases with isolated microtia with aural atresia and found to affect subcellular localization of FOXI3. Loss of function variants were found in patients with microtia and mandibular hypoplasia (CFM), suggesting dosage sensitivity of FOXI3. CONCLUSION: Damaging variants in FOXI3 are the second most frequent genetic cause of CFM, causing 1% of all cases, including 13% of familial cases in our cohort.

Tbx5 maintains atrial identity in post-natal cardiomyocytes by regulating an atrial-specific enhancer network.

Understanding how the atrial and ventricular heart chambers maintain distinct identities is a prerequisite for treating chamber-specific diseases. Here, we selectively knocked out (KO) the transcription factor Tbx5 in the atrial working myocardium to evaluate its requirement for atrial identity. Atrial Tbx5 inactivation downregulated atrial cardiomyocyte (aCM) selective gene expression. Using concurrent single nucleus transcriptome and open chromatin profiling, genomic accessibility differences were identified between control and Tbx5 KO aCMs, revealing that 69% of the control-enriched ATAC regions were bound by TBX5. Genes associated with these regions were downregulated in KO aCMs, suggesting they function as TBX5-dependent enhancers. Comparing enhancer chromatin looping using H3K27ac HiChIP identified 510 chromatin loops sensitive to TBX5 dosage, and 74.8% of control-enriched loops contained anchors in control-enriched ATAC regions. Together, these data demonstrate TBX5 maintains the atrial gene expression program by binding to and preserving the tissue-specific chromatin architecture of atrial enhancers.

Pathogenesis of Cardiomyopathy Caused by Variants in ALPK3, an Essential Pseudokinase in the Cardiomyocyte Nucleus and Sarcomere.

BACKGROUND: ALPK3 encodes α-kinase 3, a muscle-specific protein of unknown function. ALPK3 loss-of-function variants cause cardiomyopathy with distinctive clinical manifestations in both children and adults, but the molecular functions of ALPK3 remain poorly understood. METHODS: We explored the putative kinase activity of ALPK3 and the consequences of damaging variants using isogenic human induced pluripotent stem cell-derived cardiomyocytes, mice, and human patient tissues. RESULTS: Multiple sequence alignment of all human α-kinase domains and their orthologs revealed 4 conserved residues that were variant only in ALPK3, demonstrating evolutionary divergence of the ALPK3 α-kinase domain sequence. Phosphoproteomic evaluation of both ALPK3 kinase domain inhibition and overexpression failed to detect significant changes in catalytic activity, establishing ALPK3 as a pseudokinase. Investigations into alternative functions revealed that ALPK3 colocalized with myomesin proteins (MYOM1, MYOM2) at both the nuclear envelope and the sarcomere M-band. ALPK3 loss-of-function variants caused myomesin proteins to mislocalize and also dysregulated several additional M-band proteins involved in sarcomere protein turnover, which ultimately impaired cardiomyocyte structure and function. CONCLUSIONS: ALPK3 is an essential cardiac pseudokinase that inserts in the nuclear envelope and the sarcomere M-band. Loss of ALPK3 causes mislocalization of myomesins, critical force-buffering proteins in cardiomyocytes, and also dysregulates M-band proteins necessary for sarcomere protein turnover. We conclude that ALPK3 cardiomyopathy induces ventricular dilatation caused by insufficient myomesin-mediated force buffering and hypertrophy by impairment of sarcomere proteostasis.

Histone deacetylase inhibitors enhance oncolytic herpes simplex virus therapy for malignant meningioma.

Approximately 20% of meningiomas are not benign (higher grade) and tend to relapse after surgery and radiation therapy. Malignant (anaplastic) meningioma (MM) is a minor subset of high-grade meningioma that is lethal with no effective treatment options currently. Oncolytic herpes simplex virus (oHSV) is a powerful anti-cancer modality that induces both direct cell death and anti-tumor immunity, and has shown activity in preclinical models of MM. However, clinically meaningful efficacy will likely entail rational mechanistic combination approaches. We here show that epigenome modulator histone deacetylase inhibitors (HDACi) increase anti-cancer effects of oHSV in human MM models, IOMM-Lee (NF2 wild-type) and CH157 (NF2 mutant). Minimally toxic, sub-micromolar concentrations of pan-HDACi, Trichostatin A and Panobinostat, substantively increased the infectability and spread of oHSV G47Δ within MM cells in vitro, resulting in enhanced oHSV-mediated killing of target cells when infected at low multiplicity of infection (MOI). Transcriptomics analysis identified selective alteration of mRNA processing and splicing modules that might underlie the potent anti-MM effects of combining HDACi and oHSV. In vivo, HDACi treatment increased intratumoral oHSV replication and boosted the capacity of oHSV to control the growth of human MM xenografts. Thus, our work supports further translational development of the combination approach employing HDACi and oHSV for the treatment of MM.

Multiplexed Single-Nucleus RNA Sequencing Using Lipid-Oligo Barcodes.

This protocol describes a robust pipeline for simultaneously analyzing multiple samples by single-nucleus (sn)RNA-seq. cDNA obtained from each single sample are labeled with the same lipid-coupled oligonucleotide barcode (10X Genomics). Nuclei from as many as 12 individual samples can be pooled together and simultaneously processed for cDNA library construction and subsequent DNA sequencing. While previous protocols using lipid-coupled oligonucleotide barcodes were optimized for analysis of samples consisting of viable cells, this protocol is optimized for analyses of quick-frozen cell samples. The protocol ensures efficient recovery of nuclei both by incorporating high sucrose buffered solutions and by including a tracking dye (trypan blue) during nuclei isolation. The protocol also describes a procedure for removing single nuclei 'artifacts' by removing cell debris prior to single nuclear fractionation. This protocol informs the use of computational tools for filtering poorly labeled nuclei and assigning sample identity using barcode unique molecular identifier (UMI) read counts percentages. The computational pipeline is applicable to either cultured or primary, fresh or frozen cells, regardless of their cell types and species. Overall, this protocol reduces batch effects and experimental costs while enhancing sample comparison. © 2022 Wiley Periodicals LLC. Basic Protocol 1: Labeling cells with lipid oligo barcodes and generating multiplexed single-nucleus RNA-seq libraries Basic Protocol 2: Bioinformatic deconvolution of the multiplexed snRNAseq libraries.

Pathogenic variants damage cell composition and single cell transcription in cardiomyopathies.

Pathogenic variants in genes that cause dilated cardiomyopathy (DCM) and arrhythmogenic cardiomyopathy (ACM) convey high risks for the development of heart failure through unknown mechanisms. Using single-nucleus RNA sequencing, we characterized the transcriptome of 880,000 nuclei from 18 control and 61 failing, nonischemic human hearts with pathogenic variants in DCM and ACM genes or idiopathic disease. We performed genotype-stratified analyses of the ventricular cell lineages and transcriptional states. The resultant DCM and ACM ventricular cell atlas demonstrated distinct right and left ventricular responses, highlighting genotype-associated pathways, intercellular interactions, and differential gene expression at single-cell resolution. Together, these data illuminate both shared and distinct cellular and molecular architectures of human heart failure and suggest candidate therapeutic targets.

Ablation of lysophosphatidic acid receptor 1 attenuates hypertrophic cardiomyopathy in a mouse model.

Myocardial fibrosis is a key pathologic feature of hypertrophic cardiomyopathy (HCM). However, the fibrotic pathways activated by HCM-causing sarcomere protein gene mutations are poorly defined. Because lysophosphatidic acid is a mediator of fibrosis in multiple organs and diseases, we tested the role of the lysophosphatidic acid pathway in HCM. Lysphosphatidic acid receptor 1 (LPAR1), a cell surface receptor, is required for lysophosphatidic acid mediation of fibrosis. We bred HCM mice carrying a pathogenic myosin heavy-chain variant (403+/-) with Lpar1-ablated mice to create mice carrying both genetic changes (403+/- LPAR1 -/-) and assessed development of cardiac hypertrophy and fibrosis. Compared with 403+/- LPAR1WT, 403+/- LPAR1 -/- mice developed significantly less hypertrophy and fibrosis. Single-nucleus RNA sequencing of left ventricular tissue demonstrated that Lpar1 was predominantly expressed by lymphatic endothelial cells (LECs) and cardiac fibroblasts. Lpar1 ablation reduced the population of LECs, confirmed by immunofluorescence staining of the LEC markers Lyve1 and Ccl21a and, by in situ hybridization, for Reln and Ccl21a. Lpar1 ablation also altered the distribution of fibroblast cell states. FB1 and FB2 fibroblasts decreased while FB0 and FB3 fibroblasts increased. Our findings indicate that Lpar1 is expressed predominantly by LECs and fibroblasts in the heart and is required for development of hypertrophy and fibrosis in an HCM mouse model. LPAR1 antagonism, including agents in clinical trials for other fibrotic diseases, may be beneficial for HCM.

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.