Variants in a cis-regulatory element of TBX1 in conotruncal heart defect patients impair GATA6-mediated transactivation

Novel pathogenic GATA6 variant associated with congenital heart disease, diabetes mellitus and necrotizing enterocolitis

BACKGROUND

Pathogenic GATA6 variants have been associated with congenital heart disease (CHD) and a spectrum of extracardiac abnormalities, including pancreatic agenesis, congenital diaphragmatic hernia, and developmental delay. However, the comprehensive genotype-phenotype correlation of pathogenic GATA6 variation in humans remains to be fully understood.

METHODS

Exome sequencing was performed in a family where four members had CHD. In vitro functional analysis of the GATA6 variant was performed using immunofluorescence, western blot, and dual-luciferase reporter assay.

RESULTS

A novel, heterozygous missense variant in GATA6 (c.1403 G > A; p.Cys468Tyr) segregated with affected members in a family with CHD, including three with persistent truncus arteriosus. In addition, one member had childhood onset diabetes mellitus (DM), and another had necrotizing enterocolitis (NEC) with intestinal perforation. The p.Cys468Tyr variant was located in the c-terminal zinc finger domain encoded by exon 4. The mutant protein demonstrated an abnormal nuclear localization pattern with protein aggregation and decreased transcriptional activity.

CONCLUSIONS

We report a novel, familial GATA6 likely pathogenic variant associated with CHD, DM, and NEC with intestinal perforation. These findings expand the phenotypic spectrum of pathologic GATA6 variation to include intestinal abnormalities.

IMPACT

Exome sequencing identified a novel heterozygous GATA6 variant (p.Cys468Tyr) that segregated in a family with CHD including persistent truncus arteriosus, atrial septal defects and bicuspid aortic valve. Additionally, affected members displayed extracardiac findings including childhood-onset diabetes mellitus, and uniquely, necrotizing enterocolitis with intestinal perforation in the first four days of life. In vitro functional assays demonstrated that GATA6 p.Cys468Tyr variant leads to cellular localization defects and decreased transactivation activity. This work supports the importance of GATA6 as a causative gene for CHD and expands the phenotypic spectrum of pathogenic GATA6 variation, highlighting neonatal intestinal perforation as a novel extracardiac phenotype.

Establishment of the Myeloid TBX-Code Reveals Aberrant Expression of T-Box Gene TBX1 in Chronic Myeloid Leukemia

T-box genes encode transcription factors, which control developmental processes and promote cancer if deregulated. Recently, we described the lymphoid TBX-code, which collates T-box gene activities in normal lymphopoiesis, enabling identification of members deregulated in lymphoid malignancies. Here, we have extended this analysis to cover myelopoiesis, compiling the myeloid TBX-code and, thus, highlighting which of these genes might be deregulated in myeloid tumor types. We analyzed public T-box gene expression datasets bioinformatically for normal and malignant cells. Candidate T-box-gene-expressing model cell lines were identified and examined by RQ-PCR, Western Blotting, genomic profiling, and siRNA-mediated knockdown combined with RNA-seq analysis and live-cell imaging. The established myeloid TBX-code comprised 10 T-box genes, including progenitor-cell-restricted TBX1. Accordingly, we detected aberrant expression of TBX1 in 10% of stem/progenitor-cell-derived chronic myeloid leukemia (CML) patients. The classic CML cell line K-562 expressed TBX1 at high levels and served as a model to identify TBX1 activators, including transcription factor GATA1 and genomic amplification of the TBX1 locus at 22q11; inhibitors, including BCR::ABL1 fusion and downregulated GNAI2, as well as BMP, FGF2, and WNT signaling; and the target genes CDKN1A, MIR17HG, NAV1, and TMEM38A. The establishment of the myeloid TBX-code permitted identification of aberrant TBX1 expression in subsets of CML patients and cell lines. TBX1 forms an integral part of an oncogenic regulatory network impacting proliferation, survival, and differentiation. Thus, the data spotlight novel diagnostic markers and potential therapeutic targets for this malignancy.

Prioritizing cardiovascular disease-associated variants altering NKX2-5 and TBX5 binding through an integrative computational approach

Cardiovascular diseases (CVDs) are the leading cause of death worldwide and are heavily influenced by genetic factors. Genome-wide association studies have mapped >90% of CVD-associated variants within the noncoding genome, which can alter the function of regulatory proteins, such as transcription factors (TFs). However, due to the overwhelming number of single-nucleotide polymorphisms (SNPs) (>500,000) in genome-wide association studies, prioritizing variants for in vitro analysis remains challenging. In this work, we implemented a computational approach that considers support vector machine (SVM)-based TF binding site classification and cardiac expression quantitative trait loci (eQTL) analysis to identify and prioritize potential CVD-causing SNPs. We identified 1535 CVD-associated SNPs within TF footprints and putative cardiac enhancers plus 14,218 variants in linkage disequilibrium with genotype-dependent gene expression in cardiac tissues. Using ChIP-seq data from two cardiac TFs (NKX2-5 and TBX5) in human-induced pluripotent stem cell-derived cardiomyocytes, we trained a large-scale gapped k-mer SVM model to identify CVD-associated SNPs that altered NKX2-5 and TBX5 binding. The model was tested by scoring human heart TF genomic footprints within putative enhancers and measuring in vitro binding through electrophoretic mobility shift assay. Five variants predicted to alter NKX2-5 (rs59310144, rs6715570, and rs61872084) and TBX5 (rs7612445 and rs7790964) binding were prioritized for in vitro validation based on the magnitude of the predicted change in binding and are in cardiac tissue eQTLs. All five variants altered NKX2-5 and TBX5 DNA binding. We present a bioinformatic approach that considers tissue-specific eQTL analysis and SVM-based TF binding site classification to prioritize CVD-associated variants for in vitro analysis.

An Anterior Second Heart Field Enhancer Regulates the Gene Regulatory Network of the Cardiac Outflow Tract

BACKGROUND

Conotruncal defects due to developmental abnormalities of the outflow tract (OFT) are an important cause of cyanotic congenital heart disease. Dysregulation of transcriptional programs tuned by NKX2-5 (NK2 homeobox 5), GATA6 (GATA binding protein 6), and TBX1 (T-box transcription factor 1) have been implicated in abnormal OFT morphogenesis. However, there remains no consensus on how these transcriptional programs function in a unified gene regulatory network within the OFT.

METHODS

We generated mice harboring a 226-nucleotide deletion of a highly conserved cardiac enhancer containing 2 GATA-binding sites located ≈9.4 kb upstream of the transcription start site of Nkx2-5 (Nkx2-5∆enh) using CRISPR-Cas9 gene editing and assessed phenotypes. Cardiac defects in Nkx2-5∆enh/∆enh mice were structurally characterized using histology and scanning electron microscopy, and physiologically assessed using electrocardiography, echocardiography, and optical mapping. Transcriptome analyses were performed using RNA sequencing and single-cell RNA sequencing data sets. Endogenous GATA6 interaction with and activity on the NKX2-5 enhancer was studied using chromatin immunoprecipitation sequencing and transposase-accessible chromatin sequencing in human induced pluripotent stem cell-derived cardiomyocytes.

RESULTS

Nkx2-5∆enh/∆enh mice recapitulated cyanotic conotruncal defects seen in patients with NKX2-5, GATA6, and TBX1 mutations. Nkx2-5∆enh/∆enh mice also exhibited defects in right Purkinje fiber network formation, resulting in right bundle-branch block. Enhancer deletion reduced embryonic Nkx2-5 expression selectively in the right ventricle and OFT of mutant hearts, indicating that enhancer activity is localized to the anterior second heart field. Transcriptional profiling of the mutant OFT revealed downregulation of important genes involved in OFT rotation and septation, such as Tbx1, Pitx2, and Sema3c. Endogenous GATA6 interacted with the highly conserved enhancer in human induced pluripotent stem cell-derived cardiomyocytes and in wild-type mouse hearts. We found critical dose dependency of cardiac enhancer accessibility on GATA6 gene dosage in human induced pluripotent stem cell-derived cardiomyocytes.

CONCLUSIONS

Our results using human and mouse models reveal an essential gene regulatory network of the OFT that requires an anterior second heart field enhancer to link GATA6 with NKX2-5-dependent rotation and septation gene programs.

Prioritizing Cardiovascular Disease-Associated Variants Altering NKX2-5 Binding through an Integrative Computational Approach

Cardiovascular diseases (CVDs) are the leading cause of death worldwide and are heavily influenced by genetic factors. Genome-wide association studies (GWAS) have mapped > 90% of CVD-associated variants within the non-coding genome, which can alter the function of regulatory proteins, like transcription factors (TFs). However, due to the overwhelming number of GWAS single nucleotide polymorphisms (SNPs) (>500,000), prioritizing variants for in vitro analysis remains challenging. In this work, we implemented a computational approach that considers support vector machine (SVM)-based TF binding site classification and cardiac expression quantitative trait loci (eQTL) analysis to identify and prioritize potential CVD-causing SNPs. We identified 1,535 CVD-associated SNPs that occur within human heart footprints/enhancers and 9,309 variants in linkage disequilibrium (LD) with differential gene expression profiles in cardiac tissue. Using hiPSC-CM ChIP-seq data from NKX2-5 and TBX5, two cardiac TFs essential for proper heart development, we trained a large-scale gapped k-mer SVM (LS-GKM-SVM) predictive model that can identify binding sites altered by CVD-associated SNPs. The computational predictive model was tested by scoring human heart footprints and enhancers in vitro through electrophoretic mobility shift assay (EMSA). Three variants (rs59310144, rs6715570, and rs61872084) were prioritized for in vitro validation based on their eQTL in cardiac tissue and LS-GKM-SVM prediction to alter NKX2-5 DNA binding. All three variants altered NKX2-5 DNA binding. In summary, we present a bioinformatic approach that considers tissue-specific eQTL analysis and SVM-based TF binding site classification to prioritize CVD-associated variants for in vitro experimental analysis.

Disease-associated non-coding variants alter NKX2-5 DNA-binding affinity

Genome-wide association studies (GWAS) have mapped over 90 % of disease- or trait-associated variants within the non-coding genome, like cis-regulatory elements (CREs). Non-coding single nucleotide polymorphisms (SNPs) are genomic variants that can change how DNA-binding regulatory proteins, like transcription factors (TFs), interact with the genome and regulate gene expression. NKX2-5 is a TF essential for proper heart development, and mutations affecting its function have been associated with congenital heart diseases (CHDs). However, establishing a causal mechanism between non-coding genomic variants and human disease remains challenging. To address this challenge, we identified 8475 SNPs predicted to alter NKX2-5 DNA-binding using a position weight matrix (PWM)-based predictive model. Five variants were prioritized for in vitro validation; four of them are associated with traits and diseases that impact cardiovascular health. The impact of these variants on NKX2-5 binding was evaluated with electrophoretic mobility shift assay (EMSA) using purified recombinant NKX2-5 homeodomain. Binding curves were constructed to determine changes in binding between variant and reference alleles. Variants rs7350789, rs7719885, rs747334, and rs3892630 increased binding affinity, whereas rs61216514 decreased binding by NKX2-5 when compared to the reference genome. Our findings suggest that differential TF-DNA binding affinity can be key in establishing a causal mechanism of pathogenic variants.

Evaluating High-Confidence Genes in Conotruncal Cardiac Defects by Gene Burden Analyses

Background In nonsyndromic conotruncal cardiac defects, the use of next-generation sequencing for clinical diagnosis is increasingly adopted, but gene-disease associations in research are only partially translated to diagnostic panels, suggesting a need for evidence-based consensus. Methods and Results In an exome data set of 245 patients with conotruncal cardiac defects, we performed burden analysis on a high-confidence congenital heart disease gene list (n=132) with rare (<0.01%) and ultrarare (absent in the Genome Aggregation Database) protein-altering variants. Overall, we confirmed an excess of rare variants compared with ethnicity-matched controls and identified 2 known genes (GATA6, NOTCH1) and 4 candidate genes supported by the literature (ANKRD11, DOCK6, NPHP4, and STRA6). Ultrarare variant analysis was performed in combination with 3 other published studies (n=1451) and identified 3 genes (FLT4, NOTCH1, TBX1) to be significant, whereas a subgroup analysis involving 391 Chinese subjects identified only GATA6 as significant. Conclusions We suggest that these significant genes in our rare and ultrarare burden analyses warrant prioritization for clinical testing implied for rare inherited and de novo variants. Additionally, associations on ClinVar for these genes were predominantly variants of uncertain significance. Therefore, a more stringent assessment of gene-disease associations in a larger and ethnically diverse cohort is required to be prudent for future curation of conotruncal cardiac defect genes.

Screening Method for 22q11 Deletion Syndrome Involving the Use of TaqMan qPCR for TBX1 in Patients with Conotruncal Congenital Heart Disease

22q11.2 deletion syndrome is a phenotypic spectrum that encompasses DiGeorge syndrome (OMIM: 188400) and velocardiofacial syndrome (OMIM: 192430). It is caused by a 1.5–3.0 Mb hemizygous deletion of locus 22q11.2, which leads to characteristic facies, conotruncal cardiovascular malformations, velopharyngeal insufficiency, T-lymphocyte dysfunction due to thymic aplasia, and parathyroid hypoplasia, and, less frequently, neurological manifestations such as delayed psychomotor development or schizophrenia. This study aimed to describe a screening method for the diagnosis of 22q11.2 deletion syndrome in patients with Conotruncal Congenital Heart Disease (CCHD), using qPCR to detect the copy number of the TBX1 gene in a single DNA sample. A total of 23 patients were included; 21 with a biallelic prediction of the TBX1 copy number gene and 2 with a monoallelic prediction who were suspected to be positive and subjected to MLPA confirmation. One patient (4.34%) with truncus arteriosus CCHD was confirmed to have 22q11.2 deletion syndrome. We propose this approach as a possible newborn screening method for 22q11.2 deletion syndrome in CCHD patients.