FLT4 causes developmental disorders of the cardiovascular and lymphovascular systems via pleiotropic molecular mechanisms.

AIMS: Rare, deleterious genetic variants in FLT4 are associated with Tetralogy of Fallot (TOF), the most common cyanotic congenital heart disease (CHD). Distinct genetic variants in FLT4 are also an established cause of Milroy disease, the most prevalent form of primary hereditary lymphoedema. Phenotypic features of these two conditions are non-overlapping, implying pleiotropic cellular mechanisms during development. METHODS AND RESULTS: Here, we show that FLT4 variants identified in TOF patients, when expressed in primary human endothelial cells, cause aggregation of FLT4 protein in the perinuclear endoplasmic reticulum, activating proteostatic and metabolic signalling, whereas lymphoedema-associated FLT4 variants and wildtype FLT4 do not. FLT4 TOF variants display characteristic gene expression profiles in key developmental signalling pathways, revealing a role for FLT4 in cardiogenesis distinct from its role in lymphatic development. Inhibition of proteostatic signalling abrogates these effects, identifying potential avenues for therapeutic intervention. Depletion of flt4 in zebrafish caused cardiac phenotypes of reduced heart size and altered heart looping. These phenotypes were rescued with coinjection of wildtype human FLT4 mRNA, but incompletely or not at all by mRNA harbouring FLT4 TOF variants. CONCLUSIONS: Taken together, we identify a pathogenic mechanism for FLT4 variants predisposing to TOF that is distinct from the known dominant negative mechanism of Milroy-causative variants. FLT4 variants give rise to conditions of the two circulatory subdivisions of the vascular system via distinct developmental pleiotropic molecular mechanisms. TRANSLATIONAL PERSPECTIVE: Proteostatic dysfunction, if confirmed as a mechanism of CHD pathogenesis for other predisposing genes, may identify pathways to therapeutic interventions. Distinguishing mechanistically how variants in FLT4 give rise to CHD may have potential to individualise genetic counselling in affected families.

Clinical risk associated with COVID-19 among 86000 patients with congenital heart disease.

OBJECTIVE: To determine the magnitude of any excess risk of mortality and hospitalisation due to COVID-19 infection in patients with congenital heart disease (CHD) in the UK healthcare system. METHODS: Matched case-control study within the Clinical Practice Research Datalink study of anonymised general practice records in the National Health Service in England. Patients with CHD were stratified for disease severity according to the European Society of Cardiology guidelines. Presence of a positive COVID-19 test, hospitalisation with a diagnosis of COVID-19 and COVID-19-related mortality were compared in case and control groups. RESULTS: 86 441 patients with CHD and 335 839 controls were studied. Of patients with a positive COVID-19 test, patients with CHD were more likely than controls to be hospitalised (22.4% vs 14.5%; OR=1.77 (95% CI 1.60 to 1.96); p=2.11e-28) and suffer COVID-19-related death (6.1% vs 3.8%; OR=1.60 (95% CI 1.35 to 1.89); p=7.00e-08). The excess risk of COVID-19 hospitalisation and death rose with increasing physiological severity of CHD (presence of pulmonary vascular disease and/or cyanosis), rather than anatomical complexity. CONCLUSIONS: In this study of the COVID-19 pandemic experience, using population health records in over 86000 patients with CHD in England, patients with CHD with COVID-19 were at around 50-75% higher risk of hospitalisation and mortality compared with matched controls with COVID-19. We provide the first primary care-derived estimates for COVID-19 hospitalisation and case-fatality rates in patients with CHD. Some factors predictive of worse COVID-19 outcome in general populations (such as non-white ethnic group), and other CHD-specific comorbidities (such as pulmonary hypertension), influenced outcomes among patients with CHD.

Data on cardiac lncRNA STX18-AS1 expression in developing human hearts and function during in vitro hESC-cardiomyocyte differentiation.

This article presents data concerning STX18-AS1, a long noncoding RNA gene identified from a Genome-wide association study of Atrial Septal Defect (ASD). The data describes its expression patterns in human tissues and functions in regulating cardiomyocyte differentiation in vitro. STX18-AS1 is a lncRNA with a higher abundance in developing tissues, including hearts. Its transcription distribution within the embryonic hearts during key heart septation stages supports STX18-AS1's association with risk SNPs for ASD. The CRISPR stem cell pool in which STX18-AS1 was knocked down, showed reduced CM differentiation efficiency and lower expression of key cardiac transcriptional factors. This indicated its regulative role in supporting the lineage specification from cardiac mesoderm into cardiac progenitors and cardiomyocytes. These data can benefit the understanding of human embryonic heart developmental biology, and the time-course changes of cardiac transcriptional factors during in vitro cardiomyocyte differentiation from human embryonic stem cells.

Significantly increased risk of chronic obstructive pulmonary disease amongst adults with predominantly mild congenital heart disease.

Adults with congenital heart disease (CHD) face increased risk of various comorbid diseases. Previous work on lung dysfunction in this population has mainly focused on restrictive lung disease, in patients with severe CHD phenotypes. We examined the association of mild CHD with chronic obstructive pulmonary disease (COPD) in the UK Biobank (UKB). Electronic health records (EHR) were used to identify 3385 CHD cases and 479,765 healthy controls in UKB, before performing a case-control analysis over a 20-year study period for a total of > 9.5 M person-years of follow-up. Our analysis showed that UKB participants with CHD are at substantially greater risk of developing COPD than healthy controls (8.7% vs 3.1% prevalence, unadjusted OR 2.98, 95% CI 2.63, 3.36, P = 1.40e-53). Slightly increased rates of smoking were observed amongst CHD cases, however the association with COPD was shown to be robust to adjustment for smoking and other factors known to modulate COPD risk within a multivariable-adjusted Cox regression framework (fully adjusted HR 2.21, 95% CI 1.97, 2.48, P = 5.5e-41). Care for adults with CHD should aim to mitigate their increased risk of COPD, possibly via increased smoking cessation support.

Investigation of Copy Number Variation in South African Patients With Congenital Heart Defects.

BACKGROUND: Congenital heart disease (CHD) is a leading non-infectious cause of pediatric morbidity and mortality worldwide. Although the etiology of CHD is poorly understood, genetic factors including copy number variants (CNVs) contribute to the risk of CHD in individuals of European ancestry. The presence of rare CNVs in African CHD populations is unknown. This study aimed to identify pathogenic and likely pathogenic CNVs in South African patients with CHD. METHODS: Genotyping was performed on 90 patients with nonsyndromic CHD using the Affymetrix CytoScan HD platform. These data were used to identify large, rare CNVs in known CHD-associated genes and candidate genes. RESULTS: We identified eight CNVs overlapping known CHD-associated genes (GATA4, CRKL, TBX1, FLT4, B3GAT3, NSD1) in six patients. The analysis also revealed CNVs encompassing five candidate genes likely to play a role in the development of CHD (DGCR8, KDM2A, JARID2, FSTL1, CYFIP1) in five patients. One patient was found to have 47, XXY karyotype. We report a total discovery yield of 6.7%, with 5.6% of the cohort carrying pathogenic or likely pathogenic CNVs expected to cause the observed phenotypes. CONCLUSIONS: In this study, we show that chromosomal microarray is an effective technique for identifying CNVs in African patients diagnosed with CHD and have demonstrated results similar to previous CHD genetic studies in Europeans. Novel potential CHD genes were also identified, indicating the value of genetic studies of CHD in ancestrally diverse populations.

Joint analysis of functionally related genes yields further candidates associated with Tetralogy of Fallot.

Although several genes involved in the development of Tetralogy of Fallot have been identified, no genetic diagnosis is available for the majority of patients. Low statistical power may have prevented the identification of further causative genes in gene-by-gene survey analyses. Thus, bigger samples and/or novel analytic approaches may be necessary. We studied if a joint analysis of groups of functionally related genes might be a useful alternative approach. Our reanalysis of whole-exome sequencing data identified 12 groups of genes that exceedingly contribute to the burden of Tetralogy of Fallot. Further analysis of those groups showed that genes with high-impact variants tend to interact with each other. Thus, our results strongly suggest that additional candidate genes may be found by studying the protein interaction network of known causative genes. Moreover, our results show that the joint analysis of functionally related genes can be a useful complementary approach to classical single-gene analyses.

Genetic analyses of circulating PUFA-derived mediators identifies heritable dihydroxyeicosatrienoic acid species.

Estimates of heritability are the first step in identifying a trait with substantial variation due to genetic factors. Large-scale genetic analyses can identify the DNA variants that influence the levels of circulating lipid species and the statistical technique Mendelian randomisation can use these DNA variants to address potential causality of these lipids in disease. We estimated the heritability of plasma eicosanoids, octadecanoids and docosanoids to identify those lipid species with substantial heritability. We analysed plasma lipid mediators in 31 White British families (196 participants) ascertained for high blood pressure and deeply clinically and biochemically phenotyped over a 25-year period. We found that the dihydroxyeicosatrienoic acid (DHET) species, 11,12-DHET and 14,15-DHET, products of arachidonic acid metabolism by cytochrome P450 (CYP) monooxygenase and soluble epoxide hydrolase (sEH), exhibited substantial heritability (h2 = 33%-37%; Padj<0.05). Identification of these two heritable bioactive lipid species allows for future large-scale, targeted, lipidomics-genomics analyses to address causality in cardiovascular and other diseases.

Rare GATA6 variants associated with risk of congenital heart disease phenotypes in 200,000 UK Biobank exomes.

Congenital heart disease (CHD) has a complex and largely uncharacterised genetic etiology. Using 200,000 UK Biobank (UKB) exomes, we assess the burden of ultra-rare, potentially pathogenic variants in the largest case/control cohort of predominantly mild CHD to date. We find an association with GATA6, a member of the GATA family of transcription factors that play an important role during heart development and has been linked with several CHD phenotypes previously. Several identified GATA6 variants are previously unreported and their roles in conferring risk to CHD warrants further study. We demonstrate that despite limitations regarding detailed familial phenotype information in large-scale biobank projects, through careful consideration of case and control cohorts it is possible to derive important associations.

A novel RNA-mediated mechanism causing down-regulation of insulating promoter interactions in human embryonic stem cells.

The genome-wide promoter interactome is primarily maintained and regulated by architectural proteins such as CTCF and cohesin. However, some studies suggest a role for non-coding RNAs (ncRNAs) in this process. We aimed to characterise the regulatory role of RNA-mediated promoter interactions in the control of gene expression. We integrated genome-wide datasets of RNA-chromatin and promoter-genome interactions in human embryonic stem cells (hESCs) to identify putative RNA-mediated promoter interactions. We discovered that CTCF sites were enriched in RNA-PIRs (promoter interacting regions co-localising with RNA-chromatin interaction sites) and genes interacting with RNA-PIRs containing CTCF sites showed higher expression levels. One of the long noncoding RNAs (lncRNAs) expressed in hESCs, Syntaxin 18-Antisense 1 (STX18-AS1), appeared to be involved in an insulating promoter interaction with the neighbouring gene, MSX1. By knocking down STX18-AS1, the MSX1 promoter-PIR interaction was intensified and the target gene (MSX1) expression was down-regulated. Conversely, reduced MSX1 promoter-PIR interactions, resulting from CRISPR-Cas9 deletion of the PIR, increased the expression of MSX1. We conclude that STX18-AS1 RNA antagonised local CTCF-mediated insulating promoter interactions to augment gene expression. Such down-regulation of the insulating promoter interactions by this novel mechanism may explain the higher expression of genes interacting with RNA-PIRs linked to CTCF sites.

Correction: Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease.

[This corrects the article DOI: 10.1371/journal.pgen.1009679.].

Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease.

Numerous genetic studies have established a role for rare genomic variants in Congenital Heart Disease (CHD) at the copy number variation (CNV) and de novo variant (DNV) level. To identify novel haploinsufficient CHD disease genes, we performed an integrative analysis of CNVs and DNVs identified in probands with CHD including cases with sporadic thoracic aortic aneurysm. We assembled CNV data from 7,958 cases and 14,082 controls and performed a gene-wise analysis of the burden of rare genomic deletions in cases versus controls. In addition, we performed variation rate testing for DNVs identified in 2,489 parent-offspring trios. Our analysis revealed 21 genes which were significantly affected by rare CNVs and/or DNVs in probands. Fourteen of these genes have previously been associated with CHD while the remaining genes (FEZ1, MYO16, ARID1B, NALCN, WAC, KDM5B and WHSC1) have only been associated in small cases series or show new associations with CHD. In addition, a systems level analysis revealed affected protein-protein interaction networks involved in Notch signaling pathway, heart morphogenesis, DNA repair and cilia/centrosome function. Taken together, this approach highlights the importance of re-analyzing existing datasets to strengthen disease association and identify novel disease genes and pathways.

Clinical Genetic Risk Variants Inform a Functional Protein Interaction Network for Tetralogy of Fallot.

BACKGROUND: Tetralogy of Fallot (TOF)-the most common cyanotic heart defect in newborns-has evidence of multiple genetic contributing factors. Identifying variants that are clinically relevant is essential to understand patient-specific disease susceptibility and outcomes and could contribute to delineating pathomechanisms. METHODS: Using a clinically driven strategy, we reanalyzed exome sequencing data from 811 probands with TOF, to identify rare loss-of-function and other likely pathogenic variants in genes associated with congenital heart disease. RESULTS: We confirmed a major contribution of likely pathogenic variants in FLT4 (VEGFR3 [vascular endothelial growth factor receptor 3]; n=14) and NOTCH1 (n=10) and identified 1 to 3 variants in each of 21 other genes, including ATRX, DLL4, EP300, GATA6, JAG1, NF1, PIK3CA, RAF1, RASA1, SMAD2, and TBX1. In addition, multiple loss-of-function variants provided support for 3 emerging congenital heart disease/TOF candidate genes: KDR (n=4), IQGAP1 (n=3), and GDF1 (n=8). In total, these variants were identified in 63 probands (7.8%). Using the 26 composite genes in a STRING protein interaction enrichment analysis revealed a biologically relevant network (P=3.3×10-16), with VEGFR2 (vascular endothelial growth factor receptor 2; KDR) and NOTCH1 (neurogenic locus notch homolog protein 1) representing central nodes. Variants associated with arrhythmias/sudden death and heart failure indicated factors that could influence long-term outcomes. CONCLUSIONS: The results are relevant to precision medicine for TOF. They suggest considerable clinical yield from genome-wide sequencing, with further evidence for KDR (VEGFR2) as a congenital heart disease/TOF gene and for VEGF (vascular endothelial growth factor) and Notch signaling as mechanisms in human disease. Harnessing the genetic heterogeneity of single gene defects could inform etiopathogenesis and help prioritize novel candidate genes for TOF.

Rare variants in KDR, encoding VEGF Receptor 2, are associated with tetralogy of Fallot.

PURPOSE: Rare genetic variants in KDR, encoding the vascular endothelial growth factor receptor 2 (VEGFR2), have been reported in patients with tetralogy of Fallot (TOF). However, their role in disease causality and pathogenesis remains unclear. METHODS: We conducted exome sequencing in a familial case of TOF and large-scale genetic studies, including burden testing, in >1,500 patients with TOF. We studied gene-targeted mice and conducted cell-based assays to explore the role of KDR genetic variation in the etiology of TOF. RESULTS: Exome sequencing in a family with two siblings affected by TOF revealed biallelic missense variants in KDR. Studies in knock-in mice and in HEK 293T cells identified embryonic lethality for one variant when occurring in the homozygous state, and a significantly reduced VEGFR2 phosphorylation for both variants. Rare variant burden analysis conducted in a set of 1,569 patients of European descent with TOF identified a 46-fold enrichment of protein-truncating variants (PTVs) in TOF cases compared to controls (P = 7 × 10-11). CONCLUSION: Rare KDR variants, in particular PTVs, strongly associate with TOF, likely in the setting of different inheritance patterns. Supported by genetic and in vivo and in vitro functional analysis, we propose loss-of-function of VEGFR2 as one of the mechanisms involved in the pathogenesis of TOF.

Circulating ceramides as biomarkers of cardiovascular disease: Evidence from phenotypic and genomic studies.

There is a need for new biomarkers of atherosclerotic cardiovascular disease (ACVD), the main cause of death globally. Ceramides, a class of potent bioactive lipid mediators, have signalling roles in apoptosis, cellular stress and inflammation. Recent studies have highlighted circulating ceramides as novel biomarkers of coronary artery disease, type-2 diabetes and insulin resistance. Ceramides are highly regulated by enzymatic reactions throughout the body in terms of their activity and metabolism, including production, degradation and transport. The genetic studies that have been completed to date on the main ceramide species found in circulation are described, highlighting the importance of DNA variants in genes involved in ceramide biosynthesis as key influencers of heritable, circulating ceramide levels. We also review studies of disease associations with ceramides and discuss mechanistic insights deriving from recent genomic studies. The signalling activities of ceramides in vascular inflammation and apoptosis, associations between circulating ceramides and coronary artery disease risk, type-2 diabetes and insulin resistance, and the potential importance of ceramides with regard to ACVD risk factors, such as blood pressure, lipoproteins and lifestyle factors, are also discussed.

Heritability and family-based GWAS analyses of the N-acyl ethanolamine and ceramide plasma lipidome.

Signalling lipids of the N-acyl ethanolamine (NAE) and ceramide (CER) classes have emerged as potential biomarkers of cardiovascular disease (CVD). We sought to establish the heritability of plasma NAEs (including the endocannabinoid anandamide) and CERs, to identify common DNA variants influencing the circulating concentrations of the heritable lipids, and assess causality of these lipids in CVD using 2-sample Mendelian randomization (2SMR). Nine NAEs and 16 CERs were analyzed in plasma samples from 999 members of 196 British Caucasian families, using targeted ultra-performance liquid chromatography with tandem mass spectrometry. All lipids were significantly heritable (h2 = 36-62%). A missense variant (rs324420) in the gene encoding the enzyme fatty acid amide hydrolase (FAAH), which degrades NAEs, associated at genome-wide association study (GWAS) significance (P < 5 × 10-8) with four NAEs (DHEA, PEA, LEA and VEA). For CERs, rs680379 in the SPTLC3 gene, which encodes a subunit of the rate-limiting enzyme in CER biosynthesis, associated with a range of species (e.g. CER[N(24)S(19)]; P = 4.82 × 10-27). We observed three novel associations between SNPs at the CD83, SGPP1 and DEGS1 loci, and plasma CER traits (P < 5 × 10-8). 2SMR in the CARDIoGRAMplusC4D cohorts (60 801 cases; 123 504 controls) and in the DIAGRAM cohort (26 488 cases; 83 964 controls), using the genetic instruments from our family-based GWAS, did not reveal association between genetically determined differences in CER levels and CVD or diabetes. Two of the novel GWAS loci, SGPP1 and DEGS1, suggested a casual association between CERs and a range of haematological phenotypes, through 2SMR in the UK Biobank, INTERVAL and UKBiLEVE cohorts (n = 110 000-350 000).

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.