Ultra high-resolution gene centric genomic structural analysis of a non-syndromic congenital heart defect, Tetralogy of Fallot

Identification and analysis of inflammation-related biomarkers in tetralogy of Fallot

Background

Studies have revealed that inflammatory response is relevant to the tetralogy of Fallot (TOF). However, there are no studies to systematically explore the role of the inflammation-related genes (IRGs) in TOF. Therefore, based on bioinformatics, we explored the biomarkers related to inflammation in TOF, laying a theoretical foundation for its in-depth study.

Methods

TOF-related datasets (GSE36761 and GSE35776) were downloaded from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) between TOF and control groups were identified in GSE36761. And DEGs between TOF and control groups were intersected with IRGs to obtain differentially expressed IRGs (DE-IRGs). Afterwards, the least absolute shrinkage and selection operator (LASSO) and random forest (RF) were utilized to identify the biomarkers. Next, immune analysis was carried out. The transcription factor (TF)-mRNA, lncRNA-miRNA-mRNA, and miRNA-single nucleotide polymorphism (SNP)-mRNA networks were created. Finally, the potential drugs targeting the biomarkers were predicted.

Results

There were 971 DEGs between TOF and control groups, and 29 DE-IRGs were gained through the intersection between DEGs and IRGs. Next, a total of five biomarkers (MARCO, CXCL6, F3, SLC7A2, and SLC7A1) were acquired via two machine learning algorithms. Infiltrating abundance of 18 immune cells was significantly different between TOF and control groups, such as activated B cells, neutrophil, CD56dim natural killer cells, etc. The TF-mRNA network contained 4 mRNAs, 31 TFs, and 33 edges, for instance, ELF1-CXCL6, CBX8-SLC7A2, ZNF423-SLC7A1, ZNF71-F3. The lncRNA-miRNA-mRNA network was created, containing 4 mRNAs, 4 miRNAs, and 228 lncRNAs. Afterwards, nine SNPs locations were identified in the miRNA-SNP-mRNA network. A total of 21 drugs were predicted, such as ornithine, lysine, arginine, etc.

Conclusions

Our findings detected five inflammation-related biomarkers (MARCO, CXCL6, F3, SLC7A2, and SLC7A1) for TOF, providing a scientific reference for further studies of TOF.

A CAMK2B variant associated with tetralogy of Fallot, developmental delay, and growth retardation

CAMK2B encodes the beta-subunit of calcium/calmodulin-dependent protein kinase II (CAMKII), which is expressed mainly in the brain. Variants of CAMK2A and CAMK2B cause neurodevelopmental disorders, and CAMK2B alterations have been described in at least 14 patients with intellectual disability and developmental delay. Here, we describe a novel CAMK2B variant in a patient with tetralogy of Fallot (TOF), developmental delay, and growth retardation. The patient was a 2-year-old female. She was delivered at 36 weeks 6 days gestational age by caesarean section due to non-reassuring fetal status, with birth weight of 1680 g (-2.0 SD), birth length of 43.5 cm (-1.5 SD), and occipital-frontal head circumference (OFC) of 29.4 cm (-1.7 SD). Growth retardation, microcephaly, developmental delay, tetralogy of Fallot, and dysmorphic features were present. The patient controlled her head position at four months, rolled at six months, sat at 13 months, crawled at 18 months, and walked with support at 21 months. She began speaking words at 2 years old. Her dysmorphic features included a wide face, broad forehead, puffy eyelids, broad nasal base, broad and prominent philtrum, pointed chin, full cheeks, and prominent ears. A de novo missense CAMK2B variant (NM_172079.2:c.895A > G (p.Lys299Glu) NC_000007.14:g.44241708T > C (hg38)) was identified by proband exome sequencing and confirmed by Sanger sequencing. The variant was located at an autoregulatory segment and highly conserved among species. This patient displayed many of the physical features of CAMK2B-related neurodevelopmental disorder (NDD), but the TOF present in the current case is not a feature of patients with the NDD. Since a de novo CAMK2B (p.Leu443Val) variant has previously been found in a cohort of TOF, we conclude that CAMK2B variants may be associated with this specific cardiac defect.

Mutations in genes related to myocyte contraction and ventricular septum development in non-syndromic tetralogy of Fallot

Objective

Eighty percent of patients with a diagnosis of tetralogy of Fallot (TOF) do not have a known genetic etiology or syndrome. We sought to identify key molecular pathways and biological processes that are enriched in non-syndromic TOF, the most common form of cyanotic congenital heart disease, rather than single driver genes to elucidate the pathogenesis of this disease.

Methods

We undertook exome sequencing of 362 probands with non-syndromic TOF and their parents within the Pediatric Cardiac Genomics Consortium (PCGC). We identified rare (minor allele frequency <1 × 10-4), de novo variants to ascertain pathways and processes affected in this population to better understand TOF pathogenesis. Pathways and biological processes enriched in the PCGC TOF cohort were compared to 317 controls without heart defects (and their parents) from the Simons Foundation Autism Research Initiative (SFARI).

Results

A total of 120 variants in 117 genes were identified as most likely to be deleterious, with CHD7, CLUH, UNC13C, and WASHC5 identified in two probands each. Gene ontology analyses of these variants using multiple bioinformatic tools demonstrated significant enrichment in processes including cell cycle progression, chromatin remodeling, myocyte contraction and calcium transport, and development of the ventricular septum and ventricle. There was also a significant enrichment of target genes of SOX9, which is critical in second heart field development and whose loss results in membranous ventricular septal defects related to disruption of the proximal outlet septum. None of these processes was significantly enriched in the SFARI control cohort.

Conclusion

Innate molecular defects in cardiac progenitor cells and genes related to their viability and contractile function appear central to non-syndromic TOF pathogenesis. Future research utilizing our results is likely to have significant implications in stratification of TOF patients and delivery of personalized clinical care.

Identification of Novel Rare Copy Number Variants Associated with Sporadic Tetralogy of Fallot and Clinical Implications

Tetralogy of Fallot (TOF) is the most common cyanotic congenital heart disease. Highly penetrant copy number variants (CNVs) and genes related to the etiology of TOF likely exist with differences among populations. We aimed to identify CNV contributions to sporadic TOF cases in Han Chinese. Genomic DNA was extracted from peripheral blood in 605 subjects (303 sporadic TOF and 302 unaffected Han Chinese [Control] from cardiac centers in China and analyzed by genome-wide association study (GWAS). The GWAS results were compared to existing Database of Genetic Variants. These CNVs were further validated by qPCR. Bioinformatics analyses were performed with Protein-Protein Interaction (PPI) network and KEGG pathway enrichment. Across all chromosomes 119 novel "TOF-specific CNVs" were identified with prevalence of CNVs of 21.5% in chromosomes 1-20 and 37.0% including Chr21/22. In chromosomes 1-20, CNVs on 11q25 (encompasses genes ACAD8, B3GAT1, GLB1L2, GLB1L3, IGSF9B, JAM3, LOC100128239, LOC283177, MIR4697, MIR4697HG, NCAPD3, OPCML, SPATA19, THYN1, and VPS26B) and 14q32.33 (encompasses genes THYN1, OPCML, and NCAPD3) encompass genes most likely to be associated with TOF. Specific CNVs found on the chromosome 21 (6.3%) and 22(11.9%) were also identified in details. PPI network analysis identified the genes covering the specific CNVs related to TOF and the signaling pathways. This study for first time identified novel TOF-specific CNVs in the Han Chinese with higher frequency than in Caucasians and with 11q25 and 14q32.33 not reported in TOF of Caucasians. These novel CNVs identify new candidate genes for TOF and provide new insights into genetic basis of TOF.

Identification of circRNA–miRNA–mRNA Regulatory Network and Crucial Signaling Pathway Axis Involved in Tetralogy of Fallot

Tetralogy of Fallot (TOF) is one of the most common cyanotic congenital heart diseases (CHD) worldwide; however, its pathogenesis remains unclear. Recent studies have shown that circular RNAs (circRNAs) act as “sponges” for microRNAs (miRNAs) to compete for endogenous RNA (ceRNA) and play important roles in regulating gene transcription and biological processes. However, the mechanism of ceRNA in TOF remains unclear. To explore the crucial regulatory connections and pathways of TOF, we obtained the human TOF gene, miRNA, and circRNA expression profiling datasets from the Gene Expression Omnibus (GEO) database. After data pretreatment, differentially expressed mRNAs (DEmRNAs), microRNAs (DEmiRNAs), and circRNAs (DEcircRNAs) were identified between the TOF and healthy groups, and a global triple ceRNA regulatory network, including circRNAs, miRNAs, and mRNAs based on the integrated data, was constructed. A functional enrichment analysis was performed on the Metascape website to explore the biological functions of the selected genes. Then, we constructed a protein-protein interaction (PPI) network and identified seven hub genes using the cytoHubba and MCODE plug-ins in the Cytoscape software, including BCL2L11, PIK3R1, SOCS3, OSMR, STAT3, RUNX3, and IL6R. Additionally, a circRNA–miRNA–hub gene subnetwork was established, and its enrichment analysis results indicated that the extrinsic apoptotic signaling pathway, JAK-STAT signaling pathway and PI3K-Akt signaling pathway may be involved in the pathogenesis of TOF. We further identified the hsa_circ_000601/hsa-miR-148a/BCL2L11 axis as a crucial signaling pathway axis from the subnetwork. This study provides a novel regulatory network for the pathogenesis of TOF, revealing the possible molecular mechanisms and crucial regulatory pathways that may provide new strategies for candidate diagnostic biomarkers or potential therapeutic targets for TOF.

Construction and investigation of a circRNA-associated ceRNA regulatory network in Tetralogy of Fallot

Background

As the most frequent type of cyanotic congenital heart disease (CHD), tetralogy of Fallot (TOF) has a relatively poor prognosis without corrective surgery. Circular RNAs (circRNAs) represent a novel class of endogenous noncoding RNAs that regulate target gene expression posttranscriptionally in heart development. Here, we investigated the potential role of the ceRNA network in the pathogenesis of TOF.

Methods

To identify circRNA expression profiles in TOF, microarrays were used to screen the differentially expressed circRNAs between 3 TOF and 3 control human myocardial tissue samples. Then, a dysregulated circRNA-associated ceRNA network was constructed using the established multistep screening strategy.

Results

In summary, a total of 276 differentially expressed circRNAs were identified, including 214 upregulated and 62 downregulated circRNAs in TOF samples. By constructing the circRNA-associated ceRNA network based on bioinformatics data, a total of 19 circRNAs, 9 miRNAs, and 34 mRNAs were further screened. Moreover, by enlarging the sample size, the qPCR results validated the positive correlations between hsa_circ_0007798 and HIF1A.

Conclusions

The findings in this study provide a comprehensive understanding of the ceRNA network involved in TOF biology, such as the hsa_circ_0007798/miR-199b-5p/HIF1A signalling axis, and may offer candidate diagnostic biomarkers or potential therapeutic targets for TOF. In addition, we propose that the ceRNA network regulates TOF progression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12872-021-02217-w.

The Needle in the Haystack-Searching for Genetic and Epigenetic Differences in Monozygotic Twins Discordant for Tetralogy of Fallot

Congenital heart defects (CHDs) are the most common birth defect in human with an incidence of almost 1% of all live births. Most cases have a multifactorial origin with both genetics and the environment playing a role in its development and progression. Adding an epigenetic component to this aspect is exemplified by monozygotic twins which share the same genetic background but have a different disease status. As a result, the interplay between the genetic, epigenetic and the environmental conditions might contribute to the etiology and phenotype. To date, the underlying causes of the majority of CHDs remain poorly understood. In this study, we performed genome-wide high-throughput sequencing to examine the genetic, structural genomic and epigenetic differences of two identical twin pairs discordant for Tetralogy of Fallot (TOF), representing the most common cyanotic form of CHDs. Our results show the almost identical genetic and structural genomic identity of the twins. In contrast, several epigenetic alterations could be observed given by DNA methylation changes in regulatory regions of known cardiac-relevant genes. Overall, this study provides first insights into the impact of genetic and especially epigenetic factors underlying monozygotic twins discordant for CHD like TOF.

Whole Exome Sequencing Reveals the Major Genetic Contributors to Nonsyndromic Tetralogy of Fallot

RATIONALE

Familial recurrence studies provide strong evidence for a genetic component to the predisposition to sporadic, nonsyndromic Tetralogy of Fallot (TOF), the most common cyanotic congenital heart disease phenotype. Rare genetic variants have been identified as important contributors to the risk of congenital heart disease, but relatively small numbers of TOF cases have been studied to date.

OBJECTIVE

We used whole exome sequencing to assess the prevalence of unique, deleterious variants in the largest cohort of nonsyndromic TOF patients reported to date.

METHODS AND RESULTS

Eight hundred twenty-nine TOF patients underwent whole exome sequencing. The presence of unique, deleterious variants was determined; defined by their absence in the Genome Aggregation Database and a scaled combined annotation-dependent depletion score of ≥20. The clustering of variants in 2 genes, NOTCH1 and FLT4, surpassed thresholds for genome-wide significance (assigned as P<5×10^-8) after correction for multiple comparisons. NOTCH1 was most frequently found to harbor unique, deleterious variants. Thirty-one changes were observed in 37 probands (4.5%; 95% CI, 3.2%-6.1%) and included 7 loss-of-function variants 22 missense variants and 2 in-frame indels. Sanger sequencing of the unaffected parents of 7 cases identified 5 de novo variants. Three NOTCH1 variants (p.G200R, p.C607Y, and p.N1875S) were subjected to functional evaluation, and 2 showed a reduction in Jagged1-induced NOTCH signaling. FLT4 variants were found in 2.4% (95% CI, 1.6%-3.8%) of TOF patients, with 21 patients harboring 22 unique, deleterious variants. The variants identified were distinct to those that cause the congenital lymphoedema syndrome Milroy disease. In addition to NOTCH1, FLT4 and the well-established TOF gene, TBX1, we identified potential association with variants in several other candidates, including RYR1, ZFPM1, CAMTA2, DLX6, and PCM1.

CONCLUSIONS

The NOTCH1 locus is the most frequent site of genetic variants predisposing to nonsyndromic TOF, followed by FLT4. Together, variants in these genes are found in almost 7% of TOF patients.

De novo large rare copy-number variations contribute to conotruncal heart disease in Chinese patients

Conotruncal heart anomalies (CTDs) are particularly prevalent congenital heart diseases (CHD) in Hong Kong. We surveyed large (>500 kb), rare (<1% frequency in controls) copy-number variations (CNVs) in Chinese patients with CTDs to identify potentially disease-causing variations. Adults who tested negative for 22q11.2 deletions were recruited from the adult CHD clinic in Hong Kong. Using a stringent calling criteria, high-confidence CNV calls were obtained, and a large control set comprising 3,987 Caucasian and 1,945 Singapore Chinese subjects was used to identify rare CNVs. Ten large rare CNVs were identified, and 3 in 108 individuals were confirmed to harbour de novo CNVs. All three patients were syndromic with a more complex phenotype, and each of these CNVs overlapped regions likely to be important in CHD. One was a 611 kb deletion at 17p13.3, telomeric to the Miller-Dieker syndrome (MDS) critical region, overlapping the NXN gene. Another was a 5 Mb deletion at 13q33.3, within a previously described critical region for CHD. A third CNV, previously unreported, was a large duplication at 2q22.3 overlapping the ZEB2 gene. The commonly reported 1q21.1 recurrent duplication was not observed in this Chinese cohort. We provide detailed phenotypic and genotypic descriptions of large rare genic CNVs that may represent CHD loci in the East Asian population. Larger samples of Chinese origin will be required to determine whether the genome-wide distribution differs from that found in predominantly European CHD cohorts.

Search for Rare Copy-Number Variants in Congenital Heart Defects Identifies Novel Candidate Genes and a Potential Role for FOXC1 in Patients With Coarctation of the Aorta

Background—

Congenital heart defects are the most frequent malformations among newborns and a frequent cause of morbidity and mortality. Although genetic variation contributes to congenital heart defects, their precise molecular bases remain unknown in the majority of patients.

Methods and Results—

We analyzed, by high-resolution array comparative genomic hybridization, 316 children with sporadic, nonsyndromic congenital heart defects, including 76 coarctation of the aorta, 159 transposition of the great arteries, and 81 tetralogy of Fallot, as well as their unaffected parents. We identified by array comparative genomic hybridization, and validated by quantitative real-time polymerase chain reaction, 71 rare de novo (n=8) or inherited (n=63) copy-number variants (CNVs; 50 duplications and 21 deletions) in patients. We identified 113 candidate genes for congenital heart defects within these CNVs, including BTRC , CHRNB3 , CSRP2BP , ERBB2 , ERMARD , GLIS3 , PLN , PTPRJ , RLN3 , and TCTE3 . No de novo CNVs were identified in patients with transposition of the great arteries in contrast to coarctation of the aorta and tetralogy of Fallot ( P =0.002; Fisher exact test). A search for transcription factor binding sites showed that 93% of the rare CNVs identified in patients with coarctation of the aorta contained at least 1 gene with FOXC1-binding sites. This significant enrichment ( P <0.0001; permutation test) was not observed for the CNVs identified in patients with transposition of the great arteries and tetralogy of Fallot. We hypothesize that these CNVs may alter the expression of genes regulated by FOXC1. Foxc1 belongs to the forkhead transcription factors family, which plays a critical role in cardiovascular development in mice.

Conclusions—

These data suggest that deregulation of FOXC1 or its downstream genes play a major role in the pathogenesis of coarctation of the aorta in humans.

A systematic variant screening in familial cases of congenital heart defects demonstrates the usefulness of molecular genetics in this field

The etiology of congenital heart defect (CHD) combines environmental and genetic factors. So far, there were studies reporting on the screening of a single gene on unselected CHD or on familial cases selected for specific CHD types. Our goal was to systematically screen a proband of familial cases of CHD on a set of genetic tests to evaluate the prevalence of disease-causing variant identification. A systematic screening of GATA4, NKX2-5, ZIC3 and Multiplex ligation-dependent probe amplification (MLPA) P311 Kit was setup on the proband of 154 families with at least two cases of non-syndromic CHD. Additionally, ELN screening was performed on families with supravalvular arterial stenosis. Twenty-two variants were found, but segregation analysis confirmed unambiguously the causality of 16 variants: GATA4 (1 ×), NKX2-5 (6 ×), ZIC3 (3 ×), MLPA (2 ×) and ELN (4 ×). Therefore, this approach was able to identify the causal variant in 10.4% of familial CHD cases. This study demonstrated the existence of a de novo variant even in familial CHD cases and the impact of CHD variants on adult cardiac condition even in the absence of CHD. This study showed that the systematic screening of genetic factors is useful in familial CHD cases with up to 10.4% elucidated cases. When successful, it drastically improved genetic counseling by discovering unaffected variant carriers who are at risk of transmitting their variant and are also exposed to develop cardiac complications during adulthood thus prompting long-term cardiac follow-up. This study provides an important baseline at dawning of the next-generation sequencing era.

Genetic causes of congenital diaphragmatic hernia

Congenital diaphragmatic hernia (CDH) is a moderately prevalent birth defect that, despite advances in neonatal care, is still a significant cause of infant death, and surviving patients have significant morbidity. The goal of ongoing research to elucidate the genetic causes of CDH is to develop better treatment and ultimately prevention. CDH is a complex developmental defect that is etiologically heterogeneous. This review summarizes the recurrent genetic causes of CDH including aneuploidies, chromosome copy number variants, and single gene mutations. It also discusses strategies for genetic evaluation and genetic counseling in an era of rapidly evolving technologies in clinical genetic diagnostics.

MicroRNA-421 Dysregulation is Associated with Tetralogy of Fallot

The importance of microRNAs for maintaining stability in the developing vertebrate heart has recently become apparent. In addition, there is a growing appreciation for the significance of microRNAs in developmental pathology, including the formation of congenital heart defects. We examined the expression of microRNAs in right ventricular (RV) myocardium from infants with idiopathic tetralogy of Fallot (TOF, without a 22q11.2 deletion), and found 61 microRNAs to be significantly changed in expression in myocardium from children with TOF compared to normally developing comparison subjects (O’Brien et al. 2012). Predicted targets of microRNAs with altered expression were enriched for gene networks that regulate cardiac development. We previously derived a list of 229 genes known to be critical to heart development, and found 44 had significantly changed expression in TOF myocardium relative to normally developing myocardium. These 44 genes had significant negative correlations with 33 microRNAs, each of which also had significantly changed expression. Here, we focus on miR-421, as it is significantly upregulated in RV tissue from infants with TOF; is predicted to interact with multiple members of cardiovascular regulatory pathways; and has been shown to regulate cell proliferation. We knocked down, and over expressed miR-421 in primary cells derived from the RV of infants with TOF, and infants with normally developing hearts, respectively. We found a significant inverse correlation between the expression of miR-421 and SOX4, a key regulator of the Notch pathway, which has been shown to be important for the cardiac outflow track. These findings suggest that the dysregulation of miR-421 warrants further investigation as a potential contributor to tetralogy of Fallot.