Learning epistatic polygenic phenotypes with Boolean interactions.

Detecting epistatic drivers of human phenotypes is a considerable challenge. Traditional approaches use regression to sequentially test multiplicative interaction terms involving pairs of genetic variants. For higher-order interactions and genome-wide large-scale data, this strategy is computationally intractable. Moreover, multiplicative terms used in regression modeling may not capture the form of biological interactions. Building on the Predictability, Computability, Stability (PCS) framework, we introduce the epiTree pipeline to extract higher-order interactions from genomic data using tree-based models. The epiTree pipeline first selects a set of variants derived from tissue-specific estimates of gene expression. Next, it uses iterative random forests (iRF) to search training data for candidate Boolean interactions (pairwise and higher-order). We derive significance tests for interactions, based on a stabilized likelihood ratio test, by simulating Boolean tree-structured null (no epistasis) and alternative (epistasis) distributions on hold-out test data. Finally, our pipeline computes PCS epistasis p-values that probabilisticly quantify improvement in prediction accuracy via bootstrap sampling on the test set. We validate the epiTree pipeline in two case studies using data from the UK Biobank: predicting red hair and multiple sclerosis (MS). In the case of predicting red hair, epiTree recovers known epistatic interactions surrounding MC1R and novel interactions, representing non-linearities not captured by logistic regression models. In the case of predicting MS, a more complex phenotype than red hair, epiTree rankings prioritize novel interactions surrounding HLA-DRB1, a variant previously associated with MS in several populations. Taken together, these results highlight the potential for epiTree rankings to help reduce the design space for follow up experiments.

Genome-wide association studies highlight novel risk loci for septal defects and left-sided congenital heart defects.

BACKGROUND: Congenital heart defects (CHD) are structural defects of the heart affecting approximately 1% of newborns. They exhibit low penetrance and non-Mendelian patterns of inheritance as varied and complex traits. While genetic factors are known to play an important role in the development of CHD, the specific genetics remain unknown for the majority of patients. To elucidate the underlying genetic risk, we performed a genome wide association study (GWAS) of CHDs in general and specific CHD subgroups using the FinnGen Release 10 (R10) (N > 393,000), followed by functional fine-mapping through eQTL and co-localization analyses using the GTEx database. RESULTS: We discovered three genome-wide significant loci associated with general CHD. Two of them were located in chromosome 17: 17q21.32 (rs2316327, intronic: LRRC37A2, Odds ratio (OR) [95% Confidence Interval (CI)] = 1.17[1.12-1.23], p = 1.5 × 10-9) and 17q25.3 (rs1293973611, nearest: BAHCC1, OR[95%CI] = 4.48[2.80-7.17], p = 7.0 × 10-10), respectively, and in addition to general CHD, the rs1293973611 locus was associated with the septal defect subtype. The third locus was in band 1p21.2 (rs35046143, nearest: PALMD, OR[95%CI] = 1.15[1.09-1.21], p = 7.1 × 10-9), and it was associated with general CHD and left-sided lesions. In the subgroup analysis, two additional loci were associated with septal defects (rs75230966 and rs6824295), and one with left-sided lesions (rs1305393195). In the eQTL analysis the variants rs2316327 (general CHD), and rs75230966 (septal defects) both located in 17q21.32 (with a LD r2 of 0.41) were both predicted to significantly associate with the expression of WNT9B in the atrial appendage tissue category. This effect was further confirmed by co-localization analysis, which also implicated WNT3 expression in the atrial appendage. A meta-analysis of general CHD together with the UK Biobank (combined N = 881,678) provided a different genome-wide significant locus in LRRC37A2; rs16941382 (OR[95%CI] = 1.15[1.11-1.20], p = 1.5 × 10-9) which is in significant LD with rs2316327. CONCLUSIONS: Our results of general CHD and different CHD subcategories identified a complex risk locus on chromosome 17 near BAHCC1 and LRRC37A2, interacting with the genes WNT9B, WNT3 and MYL4, may constitute potential novel CHD risk associated loci, warranting future experimental tests to determine their role.

AAV9:PKP2 improves heart function and survival in a Pkp2-deficient mouse model of arrhythmogenic right ventricular cardiomyopathy.

BACKGROUND: Arrhythmogenic right ventricular cardiomyopathy (ARVC) is a familial cardiac disease associated with ventricular arrhythmias and an increased risk of sudden cardiac death. Currently, there are no approved treatments that address the underlying genetic cause of this disease, representing a significant unmet need. Mutations in Plakophilin-2 (PKP2), encoding a desmosomal protein, account for approximately 40% of ARVC cases and result in reduced gene expression. METHODS: Our goal is to examine the feasibility and the efficacy of adeno-associated virus 9 (AAV9)-mediated restoration of PKP2 expression in a cardiac specific knock-out mouse model of Pkp2. RESULTS: We show that a single dose of AAV9:PKP2 gene delivery prevents disease development before the onset of cardiomyopathy and attenuates disease progression after overt cardiomyopathy. Restoration of PKP2 expression leads to a significant extension of lifespan by restoring cellular structures of desmosomes and gap junctions, preventing or halting decline in left ventricular ejection fraction, preventing or reversing dilation of the right ventricle, ameliorating ventricular arrhythmia event frequency and severity, and preventing adverse fibrotic remodeling. RNA sequencing analyses show that restoration of PKP2 expression leads to highly coordinated and durable correction of PKP2-associated transcriptional networks beyond desmosomes, revealing a broad spectrum of biological perturbances behind ARVC disease etiology. CONCLUSIONS: We identify fundamental mechanisms of PKP2-associated ARVC beyond disruption of desmosome function. The observed PKP2 dose-function relationship indicates that cardiac-selective AAV9:PKP2 gene therapy may be a promising therapeutic approach to treat ARVC patients with PKP2 mutations.

Arrhythmogenic right ventricular cardiomyopathy (ARVC) is a heart disease that leads to abnormal heartbeats and a higher risk of sudden cardiac death. ARVC is often caused by changes in a gene called PKP2, that then makes less PKP2 protein. PKP2 protein is important for the normal structure and function of the heart. Human ARVC characteristics can be mimicked in a mouse model missing this gene. Given no therapeutic option, our goal was to test if adding a working copy of PKP2 gene in the heart of this mouse model, using a technique called gene therapy that can deliver genes to cells, could improve heart function. Here, we show that a single dose of PKP2 gene therapy can improve heart function and heartbeats as well as extend lifespan in mice. PKP2 gene therapy may be a promising approach to treat ARVC patients with PKP2 mutations.

Targeting HDAC6 to treat heart failure with preserved ejection fraction in mice.

Heart failure with preserved ejection fraction (HFpEF) poses therapeutic challenges due to the limited treatment options. Building upon our previous research that demonstrates the efficacy of histone deacetylase 6 (HDAC6) inhibition in a genetic cardiomyopathy model, we investigate HDAC6's role in HFpEF due to their shared mechanisms of inflammation and metabolism. Here, we show that inhibiting HDAC6 with TYA-018 effectively reverses established heart failure and its associated symptoms in male HFpEF mouse models. Additionally, in male mice lacking Hdac6 gene, HFpEF progression is delayed and they are resistant to TYA-018's effects. The efficacy of TYA-018 is comparable to a sodium-glucose cotransporter 2 (SGLT2) inhibitor, and the combination shows enhanced effects. Mechanistically, TYA-018 restores gene expression related to hypertrophy, fibrosis, and mitochondrial energy production in HFpEF heart tissues. Furthermore, TYA-018 also inhibits activation of human cardiac fibroblasts and enhances mitochondrial respiratory capacity in cardiomyocytes. In this work, our findings show that HDAC6 impacts on heart pathophysiology and is a promising target for HFpEF treatment.

Epistasis regulates genetic control of cardiac hypertrophy.

The combinatorial effect of genetic variants is often assumed to be additive. Although genetic variation can clearly interact non-additively, methods to uncover epistatic relationships remain in their infancy. We develop low-signal signed iterative random forests to elucidate the complex genetic architecture of cardiac hypertrophy. We derive deep learning-based estimates of left ventricular mass from the cardiac MRI scans of 29,661 individuals enrolled in the UK Biobank. We report epistatic genetic variation including variants close to CCDC141, IGF1R, TTN, and TNKS. Several loci not prioritized by univariate genome-wide association analysis are identified. Functional genomic and integrative enrichment analyses reveal a complex gene regulatory network in which genes mapped from these loci share biological processes and myogenic regulatory factors. Through a network analysis of transcriptomic data from 313 explanted human hearts, we show that these interactions are preserved at the level of the cardiac transcriptome. We assess causality of epistatic effects via RNA silencing of gene-gene interactions in human induced pluripotent stem cell-derived cardiomyocytes. Finally, single-cell morphology analysis using a novel high-throughput microfluidic system shows that cardiomyocyte hypertrophy is non-additively modifiable by specific pairwise interactions between CCDC141 and both TTN and IGF1R. Our results expand the scope of genetic regulation of cardiac structure to epistasis.

Epistasis regulates genetic control of cardiac hypertrophy.

The combinatorial effect of genetic variants is often assumed to be additive. Although genetic variation can clearly interact non-additively, methods to uncover epistatic relationships remain in their infancy. We develop low-signal signed iterative random forests to elucidate the complex genetic architecture of cardiac hypertrophy. We derive deep learning-based estimates of left ventricular mass from the cardiac MRI scans of 29,661 individuals enrolled in the UK Biobank. We report epistatic genetic variation including variants close to CCDC141, IGF1R, TTN, and TNKS. Several loci not prioritized by univariate genome-wide association analysis are identified. Functional genomic and integrative enrichment analyses reveal a complex gene regulatory network in which genes mapped from these loci share biological processes and myogenic regulatory factors. Through a network analysis of transcriptomic data from 313 explanted human hearts, we show that these interactions are preserved at the level of the cardiac transcriptome. We assess causality of epistatic effects via RNA silencing of gene-gene interactions in human induced pluripotent stem cell-derived cardiomyocytes. Finally, single-cell morphology analysis using a novel high-throughput microfluidic system shows that cardiomyocyte hypertrophy is non-additively modifiable by specific pairwise interactions between CCDC141 and both TTN and IGF1R. Our results expand the scope of genetic regulation of cardiac structure to epistasis.

Rare variants in CAPN2 increase risk for isolated hypoplastic left heart syndrome.

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart defect (CHD) characterized by hypoplasia of the left ventricle and aorta along with stenosis or atresia of the aortic and mitral valves. HLHS represents only ∼4%-8% of all CHDs but accounts for ∼25% of deaths. HLHS is an isolated defect (i.e., iHLHS) in 70% of families, the vast majority of which are simplex. Despite intense investigation, the genetic basis of iHLHS remains largely unknown. We performed exome sequencing on 331 families with iHLHS aggregated from four independent cohorts. A Mendelian-model-based analysis demonstrated that iHLHS was not due to single, large-effect alleles in genes previously reported to underlie iHLHS or CHD in >90% of families in this cohort. Gene-based association testing identified increased risk for iHLHS associated with variation in CAPN2 (p = 1.8 × 10-5), encoding a protein involved in functional adhesion. Functional validation studies in a vertebrate animal model (Xenopus laevis) confirmed CAPN2 is essential for cardiac ventricle morphogenesis and that in vivo loss of calpain function causes hypoplastic ventricle phenotypes and suggest that human CAPN2707C>T and CAPN21112C>T variants, each found in multiple individuals with iHLHS, are hypomorphic alleles. Collectively, our findings show that iHLHS is typically not a Mendelian condition, demonstrate that CAPN2 variants increase risk of iHLHS, and identify a novel pathway involved in HLHS pathogenesis.

Genetic Determinants of the Interventricular Septum Are Linked to Ventricular Septal Defects and Hypertrophic Cardiomyopathy.

BACKGROUND: A large proportion of genetic risk remains unexplained for structural heart disease involving the interventricular septum (IVS) including hypertrophic cardiomyopathy and ventricular septal defects. This study sought to develop a reproducible proxy of IVS structure from standard medical imaging, discover novel genetic determinants of IVS structure, and relate these loci to diseases of the IVS, hypertrophic cardiomyopathy, and ventricular septal defect. METHODS: We estimated the cross-sectional area of the IVS from the 4-chamber view of cardiac magnetic resonance imaging in 32 219 individuals from the UK Biobank which was used as the basis of genome wide association studies and Mendelian randomization. RESULTS: Measures of IVS cross-sectional area at diastole were a strong proxy for the 3-dimensional volume of the IVS (Pearson r=0.814, P=0.004), and correlated with anthropometric measures, blood pressure, and diagnostic codes related to cardiovascular physiology. Seven loci with clear genomic consequence and relevance to cardiovascular biology were uncovered by genome wide association studies, most notably a single nucleotide polymorphism in an intron of CDKN1A (rs2376620; ß, 7.7 mm2 [95% CI, 5.8-11.0]; P=6.0×10-10), and a common inversion incorporating KANSL1 predicted to disrupt local chromatin structure (ß, 8.4 mm2 [95% CI, 6.3-10.9]; P=4.2×10-14). Mendelian randomization suggested that inheritance of larger IVS cross-sectional area at diastole was strongly associated with hypertrophic cardiomyopathy risk (pIVW=4.6×10-10) while inheritance of smaller IVS cross-sectional area at diastole was associated with risk for ventricular septal defect (pIVW=0.007). CONCLUSIONS: Automated estimates of cross-sectional area of the IVS supports discovery of novel loci related to cardiac development and Mendelian disease. Inheritance of genetic liability for either small or large IVS, appears to confer risk for ventricular septal defect or hypertrophic cardiomyopathy, respectively. These data suggest that a proportion of risk for structural and congenital heart disease can be localized to the common genetic determinants of size and shape of cardiovascular anatomy.

Oligogenic Architecture of Rare Noncoding Variants Distinguishes 4 Congenital Heart Disease Phenotypes.

BACKGROUND: Congenital heart disease (CHD) is highly heritable, but the power to identify inherited risk has been limited to analyses of common variants in small cohorts. METHODS: We performed reimputation of 4 CHD cohorts (n=55 342) to the TOPMed reference panel (freeze 5), permitting meta-analysis of 14 784 017 variants including 6 035 962 rare variants of high imputation quality as validated by whole genome sequencing. RESULTS: Meta-analysis identified 16 novel loci, including 12 rare variants, which displayed moderate or large effect sizes (median odds ratio, 3.02) for 4 separate CHD categories. Analyses of chromatin structure link 13 of the genome-wide significant loci to key genes in cardiac development; rs373447426 (minor allele frequency, 0.003 [odds ratio, 3.37 for Conotruncal heart disease]; P=1.49×10-8) is predicted to disrupt chromatin structure for 2 nearby genes BDH1 and DLG1 involved in Conotruncal development. A lead variant rs189203952 (minor allele frequency, 0.01 [odds ratio, 2.4 for left ventricular outflow tract obstruction]; P=1.46×10-8) is predicted to disrupt the binding sites of 4 transcription factors known to participate in cardiac development in the promoter of SPAG9. A tissue-specific model of chromatin conformation suggests that common variant rs78256848 (minor allele frequency, 0.11 [odds ratio, 1.4 for Conotruncal heart disease]; P=2.6×10-8) physically interacts with NCAM1 (PFDR=1.86×10-27), a neural adhesion molecule acting in cardiac development. Importantly, while each individual malformation displayed substantial heritability (observed h2 ranging from 0.26 for complex malformations to 0.37 for left ventricular outflow tract obstructive disease) the risk for different CHD malformations appeared to be separate, without genetic correlation measured by linkage disequilibrium score regression or regional colocalization. CONCLUSIONS: We describe a set of rare noncoding variants conferring significant risk for individual heart malformations which are linked to genes governing cardiac development. These results illustrate that the oligogenic basis of CHD and significant heritability may be linked to rare variants outside protein-coding regions conferring substantial risk for individual categories of cardiac malformation.

Relationship Between Ascending Thoracic Aortic Diameter and Blood Pressure: A Mendelian Randomization Study.

BACKGROUND: Observational studies identified elevated blood pressure (BP) as a strong risk factor for thoracic aortic dilation, and BP reduction is the primary medical intervention recommended to prevent progression of aortic aneurysms. However, although BP may impact aortic dilation, aortic size may also impact BP. The causal relationship between BP and thoracic aortic size has not been reliably established. METHODS: Genome-wide association studies summary statistics were obtained for BP and ascending thoracic aortic diameter (AscAoD). Causal effects of BP on AscAoD were estimated using 2-sample Mendelian randomization using a range of pleiotropy-robust methods. RESULTS: Genetically predicted increased systolic BP, diastolic BP, and mean arterial pressure all significantly associate with higher AscAoD (systolic BP: ß estimate, 0.0041 mm/mm Hg [95% CI, 0.0008-0.0074]; P=0.02, diastolic BP: ß estimate, 0.0272 mm/mm Hg [95% CI, 0.0224-0.0320]; P<0.001, and mean arterial pressure: ß estimate, 0.0168 mm/mm Hg [95% CI, 0.0130-0.0206]; P<0.001). Genetically predicted pulse pressure, meanwhile, had an inverse association with AscAoD (ß estimate, -0.0155 mm/mm Hg [95% CI, -0.0213 to -0.0096]; P<0.001). Multivariable Mendelian randomization analyses showed that genetically predicted increased mean arterial pressure and reduced pulse pressure were independently associated with AscAoD. Bidirectional Mendelian randomization demonstrated that genetically predicted AscAoD was inversely associated with pulse pressure (ß estimate, -2.0721 mm Hg/mm [95% CI, -3.1137 to -1.0306]; P<0.001) and systolic BP (ß estimate, -1.2878 mm Hg/mm [95% CI, -2.3533 to -0.2224]; P=0.02), while directly associated with diastolic BP (0.8203 mm Hg/mm [95% CI, 0.2735-1.3672]; P=0.004). CONCLUSIONS: BP likely contributes causally to ascending thoracic aortic dilation. Increased AscAoD likely contributes to lower systolic BP and pulse pressure, but not diastolic BP, consistent with the hemodynamic consequences of a reduced aortic diameter.

Maternal first trimester metabolic profile in pregnancies with transposition of the great arteries.

BACKGROUND: Higher maternal body mass index (BMI) and abnormal glucose metabolism during early pregnancy are associated with congenital heart defects in the offspring, but the exact mechanisms are unknown. METHODS: We evaluated the association between maternal first trimester metabolic profile and transposition of the great arteries (TGA) in the offspring in a matched case-control study with 100 TGA mothers and 200 controls born in Finland during 2004-2014. Cases and controls were matched by birth year, child sex, and maternal age and BMI. Serum samples collected between 10- and 14-weeks of gestation were analyzed for 73 metabolic measures. Conditional logistic regression was used to assess the risk for TGA in the offspring, and a subgroup analysis among mothers with high BMI was conducted. RESULTS: Higher concentrations of four subtypes of extremely large very-low-density lipoprotein (VLDL) particles and one of large VLDL particles were observed in TGA mothers. This finding did not reach statistical significance after multiple testing correction. The pooled odds ratio (OR) of the all metabolic variables was slightly higher in TGA mothers in the subgroup with maternal BMI over 25 (OR 1.25) and significantly higher in the subgroup with maternal BMI over 30 (OR 1.95) compared to the original population (OR 1.18). CONCLUSIONS: Our findings indicate that an abnormal maternal early pregnancy metabolic profile might be associated with TGA in the offspring, especially in obese mothers. A trend indicating altered VLDL subtype composition in TGA pregnancies warrants further research.

Maternal and perinatal obesity induce bronchial obstruction and pulmonary hypertension via IL-6-FoxO1-axis in later life.

Obesity is a pre-disposing condition for chronic obstructive pulmonary disease, asthma, and pulmonary arterial hypertension. Accumulating evidence suggests that metabolic influences during development can determine chronic lung diseases (CLD). We demonstrate that maternal obesity causes early metabolic disorder in the offspring. Here, interleukin-6 induced bronchial and microvascular smooth muscle cell (SMC) hyperproliferation and increased airway and pulmonary vascular resistance. The key anti-proliferative transcription factor FoxO1 was inactivated via nuclear exclusion. These findings were confirmed using primary SMC treated with interleukin-6 and pharmacological FoxO1 inhibition as well as genetic FoxO1 ablation and constitutive activation. In vivo, we reproduced the structural and functional alterations in offspring of obese dams via the SMC-specific ablation of FoxO1. The reconstitution of FoxO1 using IL-6-deficient mice and pharmacological treatment did not protect against metabolic disorder but prevented SMC hyperproliferation. In human observational studies, childhood obesity was associated with reduced forced expiratory volume in 1 s/forced vital capacity ratio Z-score (used as proxy for lung function) and asthma. We conclude that the interleukin-6-FoxO1 pathway in SMC is a molecular mechanism by which perinatal obesity programs the bronchial and vascular structure and function, thereby driving CLD development. Thus, FoxO1 reconstitution provides a potential therapeutic option for preventing this metabolic programming of CLD.

High heritability of ascending aortic diameter and trans-ancestry prediction of thoracic aortic disease.

Enlargement of the aorta is an important risk factor for aortic aneurysm and dissection, a leading cause of morbidity in the developed world. Here we performed automated extraction of ascending aortic diameter from cardiac magnetic resonance images of 36,021 individuals from the UK Biobank, followed by genome-wide association. We identified lead variants across 41 loci, including genes related to cardiovascular development (HAND2, TBX20) and Mendelian forms of thoracic aortic disease (ELN, FBN1). A polygenic score significantly predicted prevalent risk of thoracic aortic aneurysm and the need for surgical intervention for patients with thoracic aneurysm across multiple ancestries within the UK Biobank, FinnGen, the Penn Medicine Biobank and the Million Veterans Program (MVP). Additionally, we highlight the primary causal role of blood pressure in reducing aortic dilation using Mendelian randomization. Overall, our findings provide a roadmap for using genetic determinants of human anatomy to understand cardiovascular development while improving prediction of diseases of the thoracic aorta.

Computational estimates of annular diameter reveal genetic determinants of mitral valve function and disease.

The fibrous annulus of the mitral valve plays an important role in valvular function and cardiac physiology, while normal variation in the size of cardiovascular anatomy may share a genetic link with common and rare disease. We derived automated estimates of mitral valve annular diameter in the 4-chamber view from 32,220 MRI images from the UK Biobank at ventricular systole and diastole as the basis for GWAS. Mitral annular dimensions corresponded to previously described anatomical norms, and GWAS inclusive of 4 population strata identified 10 loci, including possibly novel loci (GOSR2, ERBB4, MCTP2, MCPH1) and genes related to cardiac contractility (BAG3, TTN, RBFOX1). ATAC-Seq of primary mitral valve tissue localized multiple variants to regions of open chromatin in biologically relevant cell types and rs17608766 to an algorithmically predicted enhancer element in GOSR2. We observed strong genetic correlation with measures of contractility and mitral valve disease and clinical correlations with heart failure, cerebrovascular disease, and ventricular arrhythmias. Polygenic scoring of mitral valve annular diameter in systole was predictive of risk mitral valve prolapse across 4 cohorts. In summary, genetic and clinical studies of mitral valve annular diameter revealed genetic determinants of mitral valve biology, while highlighting clinical associations. Polygenic determinants of mitral valve annular diameter may represent an independent risk factor for mitral prolapse. Overall, computationally estimated phenotypes derived at scale from medical imaging represent an important substrate for genetic discovery and clinical risk prediction.

Disruption of protein quality control of the human ether-à-go-go related gene K+ channel results in profound long QT syndrome.

BACKGROUND: Long QT syndrome (LQTS) is a hereditary disease that predisposes patients to life-threatening cardiac arrhythmias and sudden cardiac death. Our previous study of the human ether-à-go-go related gene (hERG)-encoded K+ channel (Kv11.1) supports an association between hERG and RING finger protein 207 (RNF207) variants in aggravating the onset and severity of LQTS, specifically T613M hERG (hERGT613M) and RNF207 frameshift (RNF207G603fs) mutations. However, the underlying mechanistic underpinning remains unknown. OBJECTIVE: The purpose of the present study was to test the role of RNF207 in the function of hERG-encoded K+ channel subunits. METHODS: Whole-cell patch-clamp experiments were performed in human embryonic kidney (HEK 293) cells and human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) together with immunofluorescent confocal and high resolution microscopy, auto-ubiquitinylation assays, and co-immunoprecipitation experiments to test the functional interactions between hERG and RNF207. RESULTS: Here, we demonstrated that RNF207 serves as an E3 ubiquitin ligase and targets misfolded hERGT613M proteins for degradation. RNF207G603fs exhibits decreased activity and hinders the normal degradation pathway; this increases the levels of hERGT613M subunits and their dominant-negative effect on the wild-type subunits, ultimately resulting in decreased current density. Similar findings are shown for hERGA614V, a known dominant-negative mutant subunit. Finally, the presence of RNF207G603fs with hERGT613M results in significantly prolonged action potential durations and reduced hERG current in human-induced pluripotent stem cell-derived cardiomyocytes. CONCLUSION: Our study establishes RNF207 as an interacting protein serving as a ubiquitin ligase for hERG-encoded K+ channel subunits. Normal function of RNF207 is critical for the quality control of hERG subunits and consequently cardiac repolarization. Moreover, our study provides evidence for protein quality control as a new paradigm in life-threatening cardiac arrhythmias in patients with LQTS.

Comprehensive Genetic Testing for Pediatric Hypertrophic Cardiomyopathy Reveals Clinical Management Opportunities and Syndromic Conditions.

Hypertrophic cardiomyopathy (HCM) has historically been diagnosed phenotypically. Through genetic testing, identification of a molecular diagnosis (MolDx) is increasingly common but the impact on pediatric patients is unknown. This was a retrospective study of next-generation sequencing data for 602 pediatric patients with a clinician-reported history of HCM. Diagnostic yield was stratified by gene and self-reported race/ethnicity. A MolDx of HCM was identified in 242 (40%) individuals. Sarcomeric genes were the highest yielding, but pathogenic and/or likely pathogenic (P/LP) variants in syndromic genes were found in 36% of individuals with a MolDx, often in patients without documented clinical suspicion for a genetic syndrome. Among all MolDx, 73% were in genes with established clinical management recommendations and 2.9% were in genes that conferred eligibility for clinical trial enrollment. Black patients were the least likely to receive a MolDx. In the current era, genetic testing can impact management of HCM, beyond diagnostics or prognostics, through disease-specific guidelines or clinical trial eligibility. Genetic testing frequently can help identify syndromes in patients for whom syndromes may not be suspected. These findings highlight the importance of pursuing broad genetic testing, independent of suspicion based on phenotype. Lower rates of MolDx in Black patients may contribute to health inequities. Further research is needed evaluating the genetics of HCM in underrepresented/underserved populations. Additionally, research related to the impact of genetic testing on clinical management of other diseases is warranted.

Single-cell transcriptomic landscape of cardiac neural crest cell derivatives during development.

The migratory cardiac neural crest cells (CNCCs) contribute greatly to cardiovascular development. A thorough understanding of the cell lineages, developmental chronology, and transcriptomic states of CNCC derivatives during normal development is essential for deciphering the pathogenesis of CNCC-associated congenital anomalies. Here, we perform single-cell transcriptomic sequencing of 34,131 CNCC-derived cells in mouse hearts covering eight developmental stages between E10.5 and P7. We report the presence of CNCC-derived mural cells that comprise pericytes and microvascular smooth muscle cells (mVSMCs). Furthermore, we identify the transition from the CNCC-derived pericytes to mVSMCs and the key regulators over the transition. In addition, our data support that many CNCC derivatives had already committed or differentiated to a specific lineage when migrating into the heart. We explore the spatial distribution of some critical CNCC-derived subpopulations with single-molecule fluorescence in situ hybridization. Finally, we computationally reconstruct the differentiation path and regulatory dynamics of CNCC derivatives. Our study provides novel insights into the cell lineages, developmental chronology, and regulatory dynamics of CNCC derivatives during development.

Congenital heart disease risk loci identified by genome-wide association study in European patients.

Genetic factors undoubtedly affect the development of congenital heart disease (CHD) but still remain ill defined. We sought to identify genetic risk factors associated with CHD and to accomplish a functional analysis of SNP-carrying genes. We performed a genome-wide association study (GWAS) of 4034 White patients with CHD and 8486 healthy controls. One SNP on chromosome 5q22.2 reached genome-wide significance across all CHD phenotypes and was also indicative for septal defects. One region on chromosome 20p12.1 pointing to the MACROD2 locus identified 4 highly significant SNPs in patients with transposition of the great arteries (TGA). Three highly significant risk variants on chromosome 17q21.32 within the GOSR2 locus were detected in patients with anomalies of thoracic arteries and veins (ATAV). Genetic variants associated with ATAV are suggested to influence the expression of WNT3, and the variant rs870142 related to septal defects is proposed to influence the expression of MSX1. We analyzed the expression of all 4 genes during cardiac differentiation of human and murine induced pluripotent stem cells in vitro and by single-cell RNA-Seq analyses of developing murine and human hearts. Our data show that MACROD2, GOSR2, WNT3, and MSX1 play an essential functional role in heart development at the embryonic and newborn stages.

A phenome-wide association study of 26 mendelian genes reveals phenotypic expressivity of common and rare variants within the general population.

The clinical evaluation of a genetic syndrome relies upon recognition of a characteristic pattern of signs or symptoms to guide targeted genetic testing for confirmation of the diagnosis. However, individuals displaying a single phenotype of a complex syndrome may not meet criteria for clinical diagnosis or genetic testing. Here, we present a phenome-wide association study (PheWAS) approach to systematically explore the phenotypic expressivity of common and rare alleles in genes associated with four well-described syndromic diseases (Alagille (AS), Marfan (MS), DiGeorge (DS), and Noonan (NS) syndromes) in the general population. Using human phenotype ontology (HPO) terms, we systematically mapped 60 phenotypes related to AS, MS, DS and NS in 337,198 unrelated white British from the UK Biobank (UKBB) based on their hospital admission records, self-administrated questionnaires, and physiological measurements. We performed logistic regression adjusting for age, sex, and the first 5 genetic principal components, for each phenotype and each variant in the target genes (JAG1, NOTCH2 FBN1, PTPN1 and RAS-opathy genes, and genes in the 22q11.2 locus) and performed a gene burden test. Overall, we observed multiple phenotype-genotype correlations, such as the association between variation in JAG1, FBN1, PTPN11 and SOS2 with diastolic and systolic blood pressure; and pleiotropy among multiple variants in syndromic genes. For example, rs11066309 in PTPN11 was significantly associated with a lower body mass index, an increased risk of hypothyroidism and a smaller size for gestational age, all in concordance with NS-related phenotypes. Similarly, rs589668 in FBN1 was associated with an increase in body height and blood pressure, and a reduced body fat percentage as observed in Marfan syndrome. Our findings suggest that the spectrum of associations of common and rare variants in genes involved in syndromic diseases can be extended to individual phenotypes within the general population.