Dual transcriptional activator and repressor roles of TBX20 regulate adult cardiac structure and function

Role of TBX20 Truncating Variants in Dilated Cardiomyopathy and Left Ventricular Noncompaction

BACKGROUND

Less than 40% of patients with dilated cardiomyopathy (DCM) have a pathogenic/likely pathogenic genetic variant identified. TBX20 has been linked to congenital heart defects; although an association with left ventricular noncompaction (LVNC) and DCM has been proposed, it is still considered a gene with limited evidence for these phenotypes. This study sought to investigate the association between the TBX20 truncating variant (TBX20tv) and DCM/LVNC.

METHODS

TBX20 was sequenced by next-generation sequencing in 7463 unrelated probands with a diagnosis of DCM or LVNC, 22 773 probands of an internal comparison group (hypertrophic cardiomyopathy, channelopathies, or aortic diseases), and 124 098 external controls (individuals from the gnomAD database). Enrichment of TBX20tv in DCM/LVNC was calculated, cosegregation was determined in selected families, and clinical characteristics and outcomes were analyzed in carriers.

RESULTS

TBX20tv was enriched in DCM/LVNC (24/7463; 0.32%) compared with internal (1/22 773; 0.004%) and external comparison groups (4/124 098; 0.003%), with odds ratios of 73.23 (95% CI, 9.90-541.45; P<0.0001) and 99.76 (95% CI, 34.60-287.62; P<0.0001), respectively. TBX20tv was cosegregated with DCM/LVNC phenotype in 21 families for a combined logarythm of the odds score of 4.53 (strong linkage). Among 57 individuals with TBX20tv (49.1% men; mean age, 35.9±20.8 years), 41 (71.9%) exhibited DCM/LVNC, of whom 14 (34.1%) had also congenital heart defects. After a median follow-up of 6.9 (95% CI, 25-75:3.6-14.5) years, 9.7% of patients with DCM/LVNC had end-stage heart failure events and 4.8% experienced malignant ventricular arrhythmias.

CONCLUSIONS

TBX20tv is associated with DCM/LVNC; congenital heart defect is also present in around one-third of cases. TBX20tv-associated DCM/LVNC is characterized by a nonaggressive phenotype, with a low incidence of major cardiovascular events. TBX20 should be considered a definitive gene for DCM and LVNC and routinely included in genetic testing panels for these phenotypes.

Transcriptome analysis of human hypertrophic cardiomyopathy reveals inhibited cardiac development pathways in children

The epidemiological, etiological, and clinical characteristics vary greatly between pediatric (P-HCM) and adult (A-HCM) hypertrophic cardiomyopathy (HCM) patients, and the understanding of the heterogeneous pathogenesis mechanisms is insufficient to date. In this study, we aimed to comprehensively assess the respective transcriptome signatures and uncover the essential differences in gene expression patterns among A-HCM and P-HCM. The transcriptome data of adults were collected from public data (GSE89714), and novel pediatric data were first obtained by RNA sequencing from 14 P-HCM and 9 infantile donor heart samples. Our study demonstrates the common signatures of myofilament or protein synthesis and calcium ion regulation pathways in HCM. Mitochondrial function is specifically dysregulated in A-HCM, whereas the inhibition of cardiac developing networks typifies P-HCM. These findings not only distinguish the transcriptome characteristics in children and adults with HCM but also reveal the potential mechanism of the higher incidence of septal defects in P-HCM patients.

The Tbx20-TLE interaction is essential for the maintenance of the second heart field

T-box transcription factor 20 (Tbx20) plays a multifaceted role in cardiac morphogenesis and controls a broad gene regulatory network. However, the mechanism by which Tbx20 activates and represses target genes in a tissue-specific and temporal manner remains unclear. Studies show that Tbx20 directly interacts with the Transducin-like Enhancer of Split (TLE) family of proteins to mediate transcriptional repression. However, a function for the Tbx20-TLE transcriptional repression complex during heart development has yet to be established. We created a mouse model with a two amino acid substitution in the Tbx20 EH1 domain, thereby disrupting the Tbx20-TLE interaction. Disruption of this interaction impaired crucial morphogenic events, including cardiac looping and chamber formation. Transcriptional profiling of Tbx20EH1Mut hearts and analysis of putative direct targets revealed misexpression of the retinoic acid pathway and cardiac progenitor genes. Further, we show that altered cardiac progenitor development and function contribute to the severe cardiac defects in our model. Our studies indicate that TLE-mediated repression is a primary mechanism by which Tbx20 controls gene expression.

Discovery of TBX20 as a Novel Gene Underlying Atrial Fibrillation

Atrial fibrillation (AF), the most prevalent type of sustained cardiac dysrhythmia globally, confers strikingly enhanced risks for cognitive dysfunction, stroke, chronic cardiac failure, and sudden cardiovascular demise. Aggregating studies underscore the crucial roles of inherited determinants in the occurrence and perpetuation of AF. However, due to conspicuous genetic heterogeneity, the inherited defects accounting for AF remain largely indefinite. Here, via whole-genome genotyping with genetic markers and a linkage assay in a family suffering from AF, a new AF-causative locus was located at human chromosome 7p14.2-p14.3, a ~4.89 cM (~4.43-Mb) interval between the markers D7S526 and D7S2250. An exome-wide sequencing assay unveiled that, at the defined locus, the mutation in the TBX20 gene, NM_001077653.2: c.695A>G; p.(His232Arg), was solely co-segregated with AF in the family. Additionally, a Sanger sequencing assay of TBX20 in another family suffering from AF uncovered a novel mutation, NM_001077653.2: c.862G>C; p.(Asp288His). Neither of the two mutations were observed in 600 unrelated control individuals. Functional investigations demonstrated that the two mutations both significantly reduced the transactivation of the target gene KCNH2 (a well-established AF-causing gene) and the ability to bind the promoter of KCNH2, while they had no effect on the nuclear distribution of TBX20. Conclusively, these findings reveal a new AF-causative locus at human chromosome 7p14.2-p14.3 and strongly indicate TBX20 as a novel AF-predisposing gene, shedding light on the mechanism underlying AF and suggesting clinical significance for the allele-specific treatment of AF patients.

Age-dependent Lamin changes induce cardiac dysfunction via dysregulation of cardiac transcriptional programs

As we age, structural changes contribute to progressive decline in organ function, which in the heart act through poorly characterized mechanisms. Taking advantage of the short lifespan and conserved cardiac proteome of the fruit fly, we found that cardiomyocytes exhibit progressive loss of Lamin C (mammalian Lamin A/C homologue) with age, coincident with decreasing nuclear size and increasing nuclear stiffness. Premature genetic reduction of Lamin C phenocopies aging's effects on the nucleus, and subsequently decreases heart contractility and sarcomere organization. Surprisingly, Lamin C reduction downregulates myogenic transcription factors and cytoskeletal regulators, possibly via reduced chromatin accessibility. Subsequently, we find a role for cardiac transcription factors in regulating adult heart contractility and show that maintenance of Lamin C, and cardiac transcription factor expression, prevents age-dependent cardiac decline. Our findings are conserved in aged non-human primates and mice, demonstrating that age-dependent nuclear remodeling is a major mechanism contributing to cardiac dysfunction.

PRDM16 Is a Compact Myocardium-Enriched Transcription Factor Required to Maintain Compact Myocardial Cardiomyocyte Identity in Left Ventricle

BACKGROUND

Left ventricular noncompaction cardiomyopathy (LVNC) was discovered half a century ago as a cardiomyopathy with excessive trabeculation and a thin ventricular wall. In the decades since, numerous studies have demonstrated that LVNC primarily has an effect on left ventricles (LVs) and is often associated with LV dilation and dysfunction. However, in part because of the lack of suitable mouse models that faithfully mirror the selective LV vulnerability in patients, mechanisms underlying the susceptibility of LVs to dilation and dysfunction in LVNC remain unknown. Genetic studies have revealed that deletions and mutations in PRDM16 (PR domain-containing 16) cause LVNC, but previous conditional Prdm16 knockout mouse models do not mirror the LVNC phenotype in patients, and the underlying molecular mechanisms by which PRDM16 deficiency causes LVNC are still unclear.

METHODS

Prdm16 cardiomyocyte-specific knockout (Prdm16cKO) mice were generated and analyzed for cardiac phenotypes. RNA sequencing and chromatin immunoprecipitation deep sequencing were performed to identify direct transcriptional targets of PRDM16 in cardiomyocytes. Single-cell RNA sequencing in combination with spatial transcriptomics was used to determine cardiomyocyte identity at the single-cell level.

RESULTS

Cardiomyocyte-specific ablation of Prdm16 in mice caused LV-specific dilation and dysfunction, as well as biventricular noncompaction, which fully recapitulated LVNC in patients. PRDM16 functioned mechanistically as a compact myocardium-enriched transcription factor that activated compact myocardial genes while repressing trabecular myocardial genes in LV compact myocardium. Consequently, Prdm16cKO LV compact myocardial cardiomyocytes shifted from their normal transcriptomic identity to a transcriptional signature resembling trabecular myocardial cardiomyocytes or neurons. Chamber-specific transcriptional regulation by PRDM16 was attributable in part to its cooperation with LV-enriched transcription factors Tbx5 and Hand1.

CONCLUSIONS

These results demonstrate that disruption of proper specification of compact cardiomyocytes may play a key role in the pathogenesis of LVNC. They also shed light on underlying mechanisms of the LV-restricted transcriptional program governing LV chamber growth and maturation, providing a tangible explanation for the susceptibility of LV in a subset of LVNC cardiomyopathies.

Understanding the molecular basis of cardiomyopathy

Inherited cardiomyopathies are a major cause of mortality and morbidity worldwide and can be caused by mutations in a wide range of proteins located in different cellular compartments. The present review is based on Dr. Ju Chen's 2021 Robert M. Berne Distinguished Lectureship of the American Physiological Society Cardiovascular Section, in which he provided an overview of the current knowledge on the cardiomyopathy-associated proteins that have been studied in his laboratory. The review provides a general summary of the proteins in different compartments of cardiomyocytes associated with cardiomyopathies, with specific focus on the proteins that have been studied in Dr. Chen's laboratory.

The Role of Tbx20 in Cardiovascular Development and Function

Tbx20 is a member of the Tbx1 subfamily of T-box-containing genes and is known to play a variety of fundamental roles in cardiovascular development and homeostasis as well as cardiac remodeling in response to pathophysiological stresses. Mutations in TBX20 are widely associated with the complex spectrum of congenital heart defects (CHDs) in humans, which includes defects in chamber septation, chamber growth, and valvulogenesis. In addition, genetic variants of TBX20 have been found to be associated with dilated cardiomyopathy and heart arrhythmia. This broad spectrum of cardiac morphogenetic and functional defects is likely due to its broad expression pattern in multiple cardiogenic cell lineages and its critical regulation of transcriptional networks during cardiac development. In this review, we summarize recent findings in our general understanding of the role of Tbx20 in regulating several important aspects of cardiac development and homeostasis and heart function.

Genomic Copy Number Variation Study of Nine Macaca Species Provides New Insights into Their Genetic Divergence, Adaptation, and Biomedical Application

Copy number variation (CNV) can promote phenotypic diversification and adaptive evolution. However, the genomic architecture of CNVs among Macaca species remains scarcely reported, and the roles of CNVs in adaptation and evolution of macaques have not been well addressed. Here, we identified and characterized 1,479 genome-wide hetero-specific CNVs across nine Macaca species with bioinformatic methods, along with 26 CNV-dense regions and dozens of lineage-specific CNVs. The genes intersecting CNVs were overrepresented in nutritional metabolism, xenobiotics/drug metabolism, and immune-related pathways. Population-level transcriptome data showed that nearly 46% of CNV genes were differentially expressed across populations and also mainly consisted of metabolic and immune-related genes, which implied the role of CNVs in environmental adaptation of Macaca. Several CNVs overlapping drug metabolism genes were verified with genomic quantitative polymerase chain reaction, suggesting that these macaques may have different drug metabolism features. The CNV-dense regions, including 15 first reported here, represent unstable genomic segments in macaques where biological innovation may evolve. Twelve gains and 40 losses specific to the Barbary macaque contain genes with essential roles in energy homeostasis and immunity defense, inferring the genetic basis of its unique distribution in North Africa. Our study not only elucidated the genetic diversity across Macaca species from the perspective of structural variation but also provided suggestive evidence for the role of CNVs in adaptation and genome evolution. Additionally, our findings provide new insights into the application of diverse macaques to drug study.

Nuclear receptor corepressors in intellectual disability and autism

Autism spectrum disorder (ASD) is characterized by neurocognitive dysfunctions, such as impaired social interaction and language learning. Gene-environment interactions have a pivotal role in ASD pathogenesis. Nuclear receptor corepressors (NCORs) are transcription co-regulators physically associated with histone deacetylases (HDACs) and many known players in ASD etiology such as transducin β-like 1 X-linked receptor 1 and methyl-CpG binding protein 2. The epigenome-modifying NCOR complex is sensitive to many ASD risk factors, including HDAC inhibitor valproic acid and a variety of endocrine factors, xenobiotic chemicals, or metabolites that can directly bind to multiple nuclear receptors. Here, we review recent studies of NCORs in neurocognition using animal models and human genetics approaches. We discuss functional interplays between NCORs and other known players in ASD etiology. It is conceivable that the NCOR complex may bridge the in utero environmental risk factors of ASD with epigenetic remodeling and can serve as a converging point for many gene-environment interactions in the pathogenesis of ASD and intellectual disability.

Tbx20 Induction Promotes Zebrafish Heart Regeneration by Inducing Cardiomyocyte Dedifferentiation and Endocardial Expansion

Heart regeneration requires replenishment of lost cardiomyocytes (CMs) and cells of the endocardial lining. However, the signaling regulation and transcriptional control of myocardial dedifferentiation and endocardial activation are incompletely understood during cardiac regeneration. Here, we report that T-Box Transcription Factor 20 (Tbx20) is induced rapidly in the myocardial wound edge in response to various sources of cardiac damages in zebrafish. Inducing Tbx20 specifically in the adult myocardium promotes injury-induced CM proliferation through CM dedifferentiation, leading to loss of CM cellular contacts and re-expression of cardiac embryonic or fetal gene programs. Unexpectedly, we identify that myocardial Tbx20 induction activates the endocardium at the injury site with enhanced endocardial cell extension and proliferation, where it induces the endocardial Bone morphogenetic protein 6 (Bmp6) signaling. Pharmacologically inactivating endocardial Bmp6 signaling reduces expression of its targets, Id1 and Id2b, attenuating the increased endocardial regeneration in tbx20-overexpressing hearts. Altogether, our study demonstrates that Tbx20 induction promotes adult heart regeneration by inducing cardiomyocyte dedifferentiation as well as non-cell-autonomously enhancing endocardial cell regeneration.

Genomic analyses implicate noncoding de novo variants in congenital heart disease

A genetic etiology is identified for one-third of patients with congenital heart disease (CHD), with 8% of cases attributable to coding de novo variants (DNVs). To assess the contribution of noncoding DNVs to CHD, we compared genome sequences from 749 CHD probands and their parents with those from 1,611 unaffected trios. Neural network prediction of noncoding DNV transcriptional impact identified a burden of DNVs in individuals with CHD (n = 2,238 DNVs) compared to controls (n = 4,177; P = 8.7 × 10-4). Independent analyses of enhancers showed an excess of DNVs in associated genes (27 genes versus 3.7 expected, P = 1 × 10-5). We observed significant overlap between these transcription-based approaches (odds ratio (OR) = 2.5, 95% confidence interval (CI) 1.1-5.0, P = 5.4 × 10-3). CHD DNVs altered transcription levels in 5 of 31 enhancers assayed. Finally, we observed a DNV burden in RNA-binding-protein regulatory sites (OR = 1.13, 95% CI 1.1-1.2, P = 8.8 × 10-5). Our findings demonstrate an enrichment of potentially disruptive regulatory noncoding DNVs in a fraction of CHD at least as high as that observed for damaging coding DNVs.

Endodermal Maternal Transcription Factors Establish Super-Enhancers during Zygotic Genome Activation

Elucidation of the sequence of events underlying the dynamic interaction between transcription factors and chromatin states is essential. Maternal transcription factors function at the top of the regulatory hierarchy to specify the primary germ layers at the onset of zygotic genome activation (ZGA). We focus on the formation of endoderm progenitor cells and examine the interactions between maternal transcription factors and chromatin state changes underlying the cell specification process. Endoderm-specific factors Otx1 and Vegt together with Foxh1 orchestrate endoderm formation by coordinated binding to select regulatory regions. These interactions occur before the deposition of enhancer histone marks around the regulatory regions, and these TFs recruit RNA polymerase II, regulate enhancer activity, and establish super-enhancers associated with important endodermal genes. Therefore, maternal transcription factors Otx1, Vegt, and Foxh1 combinatorially regulate the activity of super-enhancers, which in turn activate key lineage-specifying genes during ZGA.

How do maternal transcription factors interact with chromatin regions to coordinate the endodermal gene regulatory program? Paraiso et al. demonstrate that combinatorial binding of maternal Otx1, Vegt, and Foxh1 to select enhancers and super-enhancers in the genome controls endodermal cell fate specification during zygotic gene activation.

A Soybean bZIP Transcription Factor GmbZIP19 Confers Multiple Biotic and Abiotic Stress Responses in Plant

The basic leucine zipper (bZIP) is a plant-specific transcription factor family that plays crucial roles in response to biotic and abiotic stresses. However, little is known about the function of bZIP genes in soybean. In this study, we isolated a bZIP gene, GmbZIP19, from soybean. A subcellular localization study of GmbZIP19 revealed its nucleus localization. We showed that GmbZIP19 expression was significantly induced by ABA (abscisic acid), JA (jasmonic acid) and SA (salicylic acid), but reduced under salt and drought stress conditions. Further, GmbZIP19 overexpression Arabidopsis lines showed increased resistance to S. sclerotiorum and Pseudomonas syringae associated with upregulated ABA-, JA-, ETH- (ethephon-)and SA-induced marker genes expression, but exhibited sensitivity to salt and drought stresses in association with destroyed stomatal closure and downregulated the salt and drought stresses marker genes' expression. We generated a soybean transient GmbZIP19 overexpression line, performed a Chromatin immunoprecipitation assay and found that GmbZIP19 bound to promoters of ABA-, JA-, ETH-, and SA-induced marker genes in soybean. The yeast one-hybrid verified the combination. The current study suggested that GmbZIP19 is a positive regulator of pathogen resistance and a negative regulator of salt and drought stress tolerance.

The selector genes midline and H15 control ventral leg pattern by both inhibiting Dpp signaling and specifying ventral fate

mid and H15 encode Tbx20 transcription factors that specify ventral pattern in the Drosophila leg. We find that there are at least two pathways for mid and H15 specification of ventral fate. In the first pathway, mid and H15 negatively regulate Dpp, the dorsal signal in leg development. mid and H15 block the dorsalizing effects of Dpp signaling in the ventral leg. In loss- and gain-of-function experiments in imaginal discs, we show that mid and H15 block the accumulation of phospho-Mad, the activated form of the Drosophila pSmad1/5 homolog. In a second pathway, we find mid and H15 must also directly promote ventral fate because simultaneously blocking Dpp signaling in mid H15 mutants does not rescue the ventral to dorsal transformation in most ventral leg structures. We show that mid and H15 act as transcriptional repressors in ventral leg development. The two genes repress the Dpp target gene Dad, the laterally expressed gene Upd, and the mid VLE enhancer. This repression depends on the eh1 domain, a binding site for the Groucho co-repressor, and is likely direct because Mid localizes to target gene enhancers in PCR-ChIP assays. A mid allele mutant for the repressing domain (eh1), mid^eh1, was found to be compromised in gain-of-function assays and in rescue of mid H15 loss-of-function. We propose that mid and H15 specify ventral fate through inhibition of Dpp signaling and through coordinating the repression of genes in the ventral leg.

Tbx20 Is Required in Mid-Gestation Cardiomyocytes and Plays a Central Role in Atrial Development

Supplemental Digital Content is available in the text.

Rationale:

Mutations in the transcription factor TBX20 (T-box 20) are associated with congenital heart disease. Germline ablation of Tbx20 results in abnormal heart development and embryonic lethality by embryonic day 9.5. Because Tbx20 is expressed in multiple cell lineages required for myocardial development, including pharyngeal endoderm, cardiogenic mesoderm, endocardium, and myocardium, the cell type–specific requirement for TBX20 in early myocardial development remains to be explored.

Objective:

Here, we investigated roles of TBX20 in midgestation cardiomyocytes for heart development.

Methods and Results:

Ablation of Tbx20 from developing cardiomyocytes using a doxycycline inducible cTnTCre transgene led to embryonic lethality. The circumference of developing ventricular and atrial chambers, and in particular that of prospective left atrium, was significantly reduced in Tbx20 conditional knockout mutants. Cell cycle analysis demonstrated reduced proliferation of Tbx20 mutant cardiomyocytes and their arrest at the G1-S phase transition. Genome-wide transcriptome analysis of mutant cardiomyocytes revealed differential expression of multiple genes critical for cell cycle regulation. Moreover, atrial and ventricular gene programs seemed to be aberrantly regulated. Putative direct TBX20 targets were identified using TBX20 ChIP-Seq (chromatin immunoprecipitation with high throughput sequencing) from embryonic heart and included key cell cycle genes and atrial and ventricular specific genes. Notably, TBX20 bound a conserved enhancer for a gene key to atrial development and identity, COUP-TFII/Nr2f2 (chicken ovalbumin upstream promoter transcription factor 2/nuclear receptor subfamily 2, group F, member 2). This enhancer interacted with the NR2F2 promoter in human cardiomyocytes and conferred atrial specific gene expression in a transgenic mouse in a TBX20-dependent manner.

Conclusions:

Myocardial TBX20 directly regulates a subset of genes required for fetal cardiomyocyte proliferation, including those required for the G1-S transition. TBX20 also directly downregulates progenitor-specific genes and, in addition to regulating genes that specify chamber versus nonchamber myocardium, directly activates genes required for establishment or maintenance of atrial and ventricular identity. TBX20 plays a previously unappreciated key role in atrial development through direct regulation of an evolutionarily conserved COUPT-FII enhancer.

T-Box20 inhibits osteogenic differentiation in adipose-derived human mesenchymal stem cells: the role of T-Box20 on osteogenesis

Background

Skeletal development and its cellular function are regulated by various transcription factors. The T-box (Tbx) family of transcription factors have critical roles in cellular differentiation as well as heart and limbs organogenesis. These factors possess activator and/or repressor domains to modify the expression of target genes. Despite the obvious effects of Tbx20 on heart development, its impact on bone development is still unknown.

Methods

To investigate the consequence by forced Tbx20 expression in the osteogenic differentiation of human mesenchymal stem cells derived from adipose tissue (Ad-MSCs), these cells were transduced with a bicistronic lentiviral vector encoding Tbx20 and an enhanced green fluorescent protein.

Results

Tbx20 gene delivery system suppressed the osteogenic differentiation of Ad-MSCs, as indicated by reduction in alkaline phosphatase activity and Alizarin Red S staining. Consistently, reverse transcription-polymerase chain reaction analyses showed that Tbx20 gain-of-function reduced the expression levels of osteoblast marker genes in osteo-inductive Ad-MSCs cultures. Accordingly, Tbx20 negatively affected osteogenesis through modulating expression of key factors involved in this process.

Conclusion

The present study suggests that Tbx20 could inhibit osteogenic differentiation in adipose-derived human mesenchymal stem cells.

Transcriptomic Signatures of Postnatal and Adult Intrinsically Photosensitive Ganglion Cells

Intrinsically photosensitive retinal ganglion cells (ipRGCs) are rare mammalian photoreceptors essential for non-image-forming vision functions, such as circadian photoentrainment and the pupillary light reflex. They comprise multiple subtypes distinguishable by morphology, physiology, projections, and levels of expression of melanopsin (Opn4), their photopigment. The molecular programs that distinguish ipRGCs from other ganglion cells and ipRGC subtypes from one another remain elusive. Here, we present comprehensive gene expression profiles of early postnatal and adult mouse ipRGCs purified from two lines of reporter mice that mark different sets of ipRGC subtypes. We find dozens of novel genes highly enriched in ipRGCs. We reveal that Rasgrp1 and Tbx20 are selectively expressed in subsets of ipRGCs, though these molecularly defined groups imperfectly match established ipRGC subtypes. We demonstrate that the ipRGCs regulating circadian photoentrainment are diverse at the molecular level. Our findings reveal unexpected complexity in gene expression patterns across mammalian ipRGC subtypes.

Dynamics of cardiomyocyte transcriptome and chromatin landscape demarcates key events of heart development

Organogenesis involves dynamic regulation of gene transcription and complex multipathway interactions. Despite our knowledge of key factors regulating various steps of heart morphogenesis, considerable challenges in understanding its mechanism still exist because little is known about their downstream targets and interactive regulatory network. To better understand transcriptional regulatory mechanism driving heart development and the consequences of its disruption in vivo, we performed time-series analyses of the transcriptome and genome-wide chromatin accessibility in isolated cardiomyocytes (CMs) from wild-type zebrafish embryos at developmental stages corresponding to heart tube morphogenesis, looping, and maturation. We identified genetic regulatory modules driving crucial events of heart development that contained key cardiac TFs and are associated with open chromatin regions enriched for DNA sequence motifs belonging to the family of the corresponding TFs. Loss of function of cardiac TFs Gata5, Tbx5a, and Hand2 affected the cardiac regulatory networks and caused global changes in chromatin accessibility profile, indicating their role in heart development. Among regions with differential chromatin accessibility in mutants were highly conserved noncoding elements that represent putative enhancers driving heart development. The most prominent gene expression changes, which correlated with chromatin accessibility modifications within their proximal promoter regions, occurred between heart tube morphogenesis and looping, and were associated with metabolic shift and hematopoietic/cardiac fate switch during CM maturation. Our results revealed the dynamic regulatory landscape throughout heart development and identified interactive molecular networks driving key events of heart morphogenesis.

A promoter interaction map for cardiovascular disease genetics

Over 500 genetic loci have been associated with risk of cardiovascular diseases (CVDs); however, most loci are located in gene-distal non-coding regions and their target genes are not known. Here, we generated high-resolution promoter capture Hi-C (PCHi-C) maps in human induced pluripotent stem cells (iPSCs) and iPSC-derived cardiomyocytes (CMs) to provide a resource for identifying and prioritizing the functional targets of CVD associations. We validate these maps by demonstrating that promoters preferentially contact distal sequences enriched for tissue-specific transcription factor motifs and are enriched for chromatin marks that correlate with dynamic changes in gene expression. Using the CM PCHi-C map, we linked 1999 CVD-associated SNPs to 347 target genes. Remarkably, more than 90% of SNP-target gene interactions did not involve the nearest gene, while 40% of SNPs interacted with at least two genes, demonstrating the importance of considering long-range chromatin interactions when interpreting functional targets of disease loci.

Genome-wide use of high- and low-affinity Tbrain transcription factor binding sites during echinoderm development

Sea stars and sea urchins are model systems for interrogating the types of deep evolutionary changes that have restructured developmental gene regulatory networks (GRNs). Although cis-regulatory DNA evolution is likely the predominant mechanism of change, it was recently shown that Tbrain, a Tbox transcription factor protein, has evolved a changed preference for a low-affinity, secondary binding motif. The primary, high-affinity motif is conserved. To date, however, no genome-wide comparisons have been performed to provide an unbiased assessment of the evolution of GRNs between these taxa, and no study has attempted to determine the interplay between transcription factor binding motif evolution and GRN topology. The study here measures genome-wide binding of Tbrain orthologs by using ChIP-sequencing and associates these orthologs with putative target genes to assess global function. Targets of both factors are enriched for other regulatory genes, although nonoverlapping sets of functional enrichments in the two datasets suggest a much diverged function. The number of low-affinity binding motifs is significantly depressed in sea urchins compared with sea star, but both motif types are associated with genes from a range of functional categories. Only a small fraction (∼10%) of genes are predicted to be orthologous targets. Collectively, these data indicate that Tbr has evolved significantly different developmental roles in these echinoderms and that the targets and the binding motifs in associated cis-regulatory sequences are dispersed throughout the hierarchy of the GRN, rather than being biased toward terminal process or discrete functional blocks, which suggests extensive evolutionary tinkering.

Formation of a TBX20-CASZ1 protein complex is protective against dilated cardiomyopathy and critical for cardiac homeostasis

By the age of 40, one in five adults without symptoms of cardiovascular disease are at risk for developing congestive heart failure. Within this population, dilated cardiomyopathy (DCM) remains one of the leading causes of disease and death, with nearly half of cases genetically determined. Though genetic and high throughput sequencing-based approaches have identified sporadic and inherited mutations in a multitude of genes implicated in cardiomyopathy, how combinations of asymptomatic mutations lead to cardiac failure remains a mystery. Since a number of studies have implicated mutations of the transcription factor TBX20 in congenital heart diseases, we investigated the underlying mechanisms, using an unbiased systems-based screen to identify novel, cardiac-specific binding partners. We demonstrated that TBX20 physically and genetically interacts with the essential transcription factor CASZ1. This interaction is required for survival, as mice heterozygous for both Tbx20 and Casz1 die post-natally as a result of DCM. A Tbx20 mutation associated with human familial DCM sterically interferes with the TBX20-CASZ1 interaction and provides a physical basis for how this human mutation disrupts normal cardiac function. Finally, we employed quantitative proteomic analyses to define the molecular pathways mis-regulated upon disruption of this novel complex. Collectively, our proteomic, biochemical, genetic, and structural studies suggest that the physical interaction between TBX20 and CASZ1 is required for cardiac homeostasis, and further, that reduction or loss of this critical interaction leads to DCM. This work provides strong evidence that DCM can be inherited through a digenic mechanism.

A molecular understanding of cardiomyocyte development is an essential goal for improving clinical approaches to CHD. While TBX20 is an essential transcription factor for heart development and its disease relevance is well established, many fundamental questions remain about the mechanism of TBX20 function. Principle among these is how TBX20 mutations associated with adult dilated cardiomyopathy circumvent (DCM) the essential embryonic requirement for TBX20 in heart development. Here we report using an integrated approach that TBX20 complexes with the cardiac transcription factor CASZ1 in vivo. We confirmed TBX20 and CASZ1 interact biochemically and genetically, and show mice heterozygous for both Tbx20 and Casz1 die, beginning at 4 to 8 weeks post birth, exhibiting hallmarks of DCM. Interestingly, the human mutant TBX20^F256I bypasses the early essential requirement for TBX20 but leads to DCM. We report here that TBX20^F256I disrupts the TBX20-CASZ1 interaction, ascribing clinical relevance to this protein complex. Further, by using quantitative proteomics we have identified the molecular pathways altered in TBX20-CASZ1-mediated DCM. Together, these results identify a novel interaction between TBX20 and CASZ1 that is essential for maintaining cardiac homeostasis and imply that DCM can be inherited through a digenic mechanism.

Transcriptional regulation of FOXP3 requires integrated activation of both promoter and CNS regions in tumor-induced CD8+ Treg cells

T-regulatory cells are an upsurge in the tumor microenvironment and induce immune-evasion. CD4+ Treg cells are well characterized whereas the role of CD8+ Tregs in cancer has recently started to crease attention. Here, we report an augmentation CD8+FOXP3+ Tregs in breast tumor microenvironment. FOXP3, the lineage-specific transcription factor, is a dominant regulator of Treg cell development and function. FOXP3 is induced preferentially by divergent signaling in CD4+ Treg cells. But how FOXP3 is induced and maintained in tumor-CD8+ Tregs is the Cinderella of the investigation. We observed that RUNX3, a CD8+ lineage-specific transcription factor, binds at the FOXP3-promoter to induce its transcription. In addition to promoter activation, involvement of cis-elements CNS1 and CNS2 in the transcriptional regulation of FOXP3 was also evident in these cells. SMAD3 binds to CNS1 region and acts as transcription inducer, whereas GATA3 plays a temporal role in the FOXP3 transcription by differential chromatin modification in CNS regions. In CNS1 region, GATA3 acts as a repressor for FOXP3 in naïve CD8+ T cells. Whereas in CD8+ Tregs, GATA3 binds directly at CNS2 region and persuaded the maintenance of FOXP3. Therefore, the intervention of these concerted transcriptional machinery may have a therapeutic potential in immunotherapy of cancer.

Coordinated control of terminal differentiation and restriction of cellular plasticity

The acquisition of a specific cellular identity is usually paralleled by a restriction of cellular plasticity. Whether and how these two processes are coordinated is poorly understood. Transcription factors called terminal selectors activate identity-specific effector genes during neuronal differentiation to define the structural and functional properties of a neuron. To study restriction of plasticity, we ectopically expressed C. elegans CHE-1, a terminal selector of ASE sensory neuron identity. In undifferentiated cells, ectopic expression of CHE-1 results in activation of ASE neuron type-specific effector genes. Once cells differentiate, their plasticity is restricted and ectopic expression of CHE-1 no longer results in activation of ASE effector genes. In striking contrast, removal of the respective terminal selectors of other sensory, inter-, or motor neuron types now enables ectopically expressed CHE-1 to activate its ASE-specific effector genes, indicating that terminal selectors not only activate effector gene batteries but also control the restriction of cellular plasticity. Terminal selectors mediate this restriction at least partially by organizing chromatin. The chromatin structure of a CHE-1 target locus is less compact in neurons that lack their resident terminal selector and genetic epistasis studies with H3K9 methyltransferases suggest that this chromatin modification acts downstream of a terminal selector to restrict plasticity. Taken together, terminal selectors activate identity-specific genes and make non-identity-defining genes less accessible, thereby serving as a checkpoint to coordinate identity specification with restriction of cellular plasticity.

DOI:http://dx.doi.org/10.7554/eLife.24100.001

Different modes of enhancer-specific regulation by Runt and Even-skipped during Drosophila segmentation

Expression of the Drosophila slp1 gene depends on nonadditive interactions between two cis-regulatory enhancers. These enhancers are repressed by preventing either Pol II recruitment or release of promoter-proximal paused Pol II in a manner that is both enhancer and transcription factor specific and can account for their nonadditive interaction.

The initial metameric expression of the Drosophila sloppy paired 1 (slp1) gene is controlled by two distinct cis-regulatory DNA elements that interact in a nonadditive manner to integrate inputs from transcription factors encoded by the pair-rule segmentation genes. We performed chromatin immunoprecipitation on reporter genes containing these elements in different embryonic genotypes to investigate the mechanism of their regulation. The distal early stripe element (DESE) mediates both activation and repression by Runt. We find that the differential response of DESE to Runt is due to an inhibitory effect of Fushi tarazu (Ftz) on P-TEFb recruitment and the regulation of RNA polymerase II (Pol II) pausing. The proximal early stripe element (PESE) is also repressed by Runt, but in this case, Runt prevents PESE-dependent Pol II recruitment and preinitiation complex (PIC) assembly. PESE is also repressed by Even-skipped (Eve), but, of interest, this repression involves regulation of P-TEFb recruitment and promoter-proximal Pol II pausing. These results demonstrate that the mode of slp1 repression by Runt is enhancer specific, whereas the mode of repression of the slp1 PESE enhancer is transcription factor specific. We propose a model based on these differential regulatory interactions that accounts for the nonadditive interactions between the PESE and DESE enhancers during Drosophila segmentation.

Tbx20 controls the expression of the KCNH2 gene and of hERG channels

Long QT syndrome (LQTS) exhibits great phenotype variability among family members carrying the same mutation, which can be partially attributed to genetic factors. We functionally analyzed the KCNH2 (encoding for Kv11.1 or hERG channels) and TBX20 (encoding for the transcription factor Tbx20) variants found by next-generation sequencing in two siblings with LQTS in a Spanish family of African ancestry. Affected relatives harbor a heterozygous mutation in KCNH2 that encodes for p.T152HfsX180 Kv11.1 (hERG). This peptide, by itself, failed to generate any current when transfected into Chinese hamster ovary (CHO) cells but, surprisingly, exerted "chaperone-like" effects over native hERG channels in both CHO cells and mouse atrial-derived HL-1 cells. Therefore, heterozygous transfection of native (WT) and p.T152HfsX180 hERG channels generated a current that was indistinguishable from that generated by WT channels alone. Some affected relatives also harbor the p.R311C mutation in Tbx20. In human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs), Tbx20 enhanced human KCNH2 gene expression and hERG currents (I_hERG) and shortened action-potential duration (APD). However, Tbx20 did not modify the expression or activity of any other channel involved in ventricular repolarization. Conversely, p.R311C Tbx20 did not increase KCNH2 expression in hiPSC-CMs, which led to decreased I_hERG and increased APD. Our results suggest that Tbx20 controls the expression of hERG channels responsible for the rapid component of the delayed rectifier current. On the contrary, p.R311C Tbx20 specifically disables the Tbx20 protranscriptional activity over KCNH2 Therefore, TBX20 can be considered a KCNH2-modifying gene.

TBX5: A Key Regulator of Heart Development

TBX5 is a member of the T-box transcription factor family and is primarily known for its role in cardiac and forelimb development. Human patients with dominant mutations in TBX5 are characterized by Holt-Oram syndrome, and show defects of the cardiac septa, cardiac conduction system, and the anterior forelimb. The range of cardiac defects associated with TBX5 mutations in humans suggests multiple roles for the transcription factor in cardiac development and function. Animal models demonstrate similar defects and have provided a useful platform for investigating the roles of TBX5 during embryonic development. During early cardiac development, TBX5 appears to act primarily as a transcriptional activator of genes associated with cardiomyocyte maturation and upstream of morphological signals for septation. During later cardiac development, TBX5 is required for patterning of the cardiac conduction system and maintenance of mature cardiomyocyte function. A comprehensive understanding of the integral roles of TBX5 throughout cardiac development and adult life will be critical for understanding human cardiac morphology and function.

iPSC-derived cardiomyocytes reveal abnormal TGF-β signalling in left ventricular non-compaction cardiomyopathy

Left ventricular non-compaction (LVNC) is the third most prevalent cardiomyopathy in children and its pathogenesis has been associated with the developmental defect of the embryonic myocardium. We show that patient-specific induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) generated from LVNC patients carrying a mutation in the cardiac transcription factor TBX20 recapitulate a key aspect of the pathological phenotype at the single-cell level and this was associated with perturbed transforming growth factor beta (TGF-β) signalling. LVNC iPSC-CMs have decreased proliferative capacity due to abnormal activation of TGF-β signalling. TBX20 regulates the expression of TGF-β signalling modifiers including one known to be a genetic cause of LVNC, PRDM16, and genome editing of PRDM16 caused proliferation defects in iPSC-CMs. Inhibition of TGF-β signalling and genome correction of the TBX20 mutation were sufficient to reverse the disease phenotype. Our study demonstrates that iPSC-CMs are a useful tool for the exploration of pathological mechanisms underlying poorly understood cardiomyopathies including LVNC.

T-Box Genes in the Kidney and Urinary Tract

T-box (Tbx) genes encode an ancient group of transcription factors that play important roles in patterning, specification, proliferation, and differentiation programs in vertebrate organogenesis. This is testified by severe organ malformation syndromes in mice homozygous for engineered null alleles of specific T-box genes and by the large number of human inherited organ-specific diseases that have been linked to mutations in these genes. One of the organ systems that has not been associated with loss of specific T-box gene function in human disease for long is the excretory system. However, this has changed with the finding that mutations in TBX18, a member of a vertebrate-specific subgroup within the Tbx1-subfamily of T-box transcription factor genes, cause congenital anomalies of the kidney and urinary tract, predominantly hydroureter and ureteropelvic junction obstruction. Gene expression analyses, loss-of-function studies, and lineage tracing in the mouse suggest a primary role for this transcription factor in specifying the ureteric mesenchyme in the common anlage of the kidney, the ureter, and the bladder. We review the function of Tbx18 in ureterogenesis and discuss the body of evidence that Tbx18 and other members of the T-box gene family, namely, Tbx1, Tbx2, Tbx3, and Tbx20, play additional roles in development and homeostasis of other components of the excretory system in vertebrates.

Probing chromatin landscape reveals roles of endocardial TBX20 in septation

Mutations in the T-box transcription factor TBX20 are associated with multiple forms of congenital heart defects, including cardiac septal abnormalities, but our understanding of the contributions of endocardial TBX20 to heart development remains incomplete. Here, we investigated how TBX20 interacts with endocardial gene networks to drive the mesenchymal and myocardial movements that are essential for outflow tract and atrioventricular septation. Selective ablation of Tbx20 in murine endocardial lineages reduced the expression of extracellular matrix and cell migration genes that are critical for septation. Using the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), we identified accessible chromatin within endocardial lineages and intersected these data with TBX20 ChIP-seq and chromatin loop maps to determine that TBX20 binds a conserved long-range enhancer to regulate versican (Vcan) expression. We also observed reduced Vcan expression in Tbx20-deficient mice, supporting a direct role for TBX20 in Vcan regulation. Further, we show that the Vcan enhancer drove reporter gene expression in endocardial lineages in a TBX20-binding site-dependent manner. This work illuminates gene networks that interact with TBX20 to orchestrate cardiac septation and provides insight into the chromatin landscape of endocardial lineages during septation.

Lack of Genetic Interaction between Tbx18 and Tbx2/Tbx20 in Mouse Epicardial Development

The epicardium, the outermost layer of the heart, is an essential source of cells and signals for the formation of the cardiac fibrous skeleton and the coronary vasculature, and for the maturation of the myocardium during embryonic development. The molecular factors that control epicardial mobilization and differentiation, and direct the epicardial-myocardial cross-talk are, however, insufficiently understood. The T-box transcription factor gene Tbx18 is specifically expressed in the epicardium of vertebrate embryos. Loss of Tbx18 is dispensable for epicardial development, but may influence coronary vessel maturation. In contrast, over-expression of an activator version of TBX18 severely impairs epicardial development by premature differentiation of epicardial cells into SMCs indicating a potential redundancy of Tbx18 with other repressors of the T-box gene family. Here, we show that Tbx2 and Tbx20 are co-expressed with Tbx18 at different stages of epicardial development. Using a conditional gene targeting approach we find that neither the epicardial loss of Tbx2 nor the combined loss of Tbx2 and Tbx18 affects epicardial development. Similarly, we observed that the heterozygous loss of Tbx20 with and without additional loss of Tbx18 does not impact on epicardial integrity and mobilization in mouse embryos. Thus, Tbx18 does not function redundantly with Tbx2 or Tbx20 in epicardial development.

A novel conditional mouse model for Nkx2-5 reveals transcriptional regulation of cardiac ion channels

Nkx2-5 is one of the master regulators of cardiac development, homeostasis and disease. This transcription factor has been previously associated with a suite of cardiac congenital malformations and impairment of electrical activity. When disease causative mutations in transcription factors are considered, NKX2-5 gene dysfunction is the most common abnormality found in patients. Here we describe a novel mouse model and subsequent implications of Nkx2-5 loss for aspects of myocardial electrical activity. In this work we have engineered a new Nkx2-5 conditional knockout mouse in which flox sites flank the entire Nkx2-5 locus, and validated this line for the study of heart development, differentiation and disease using a full deletion strategy. While our homozygous knockout mice show typical embryonic malformations previously described for the lack of the Nkx2-5 gene, hearts of heterozygous adult mice show moderate morphological and functional abnormalities that are sufficient to sustain blood supply demands under homeostatic conditions. This study further reveals intriguing aspects of Nkx2-5 function in the control of cardiac electrical activity. Using a combination of mouse genetics, biochemistry, molecular and cell biology, we demonstrate that Nkx2-5 regulates the gene encoding Kcnh2 channel and others, shedding light on potential mechanisms generating electrical abnormalities observed in patients bearing NKX2-5 dysfunction and opening opportunities to the study of novel therapeutic targets for anti-arrhythmogenic therapies.

Complex Interdependence Regulates Heterotypic Transcription Factor Distribution and Coordinates Cardiogenesis

Transcription factors (TFs) are thought to function with partners to achieve specificity and precise quantitative outputs. In the developing heart, heterotypic TF interactions, such as between the T-box TF TBX5 and the homeodomain TF NKX2-5, have been proposed as a mechanism for human congenital heart defects. We report extensive and complex interdependent genomic occupancy of TBX5, NKX2-5, and the zinc finger TF GATA4 coordinately controlling cardiac gene expression, differentiation, and morphogenesis. Interdependent binding serves not only to co-regulate gene expression but also to prevent TFs from distributing to ectopic loci and activate lineage-inappropriate genes. We define preferential motif arrangements for TBX5 and NKX2-5 cooperative binding sites, supported at the atomic level by their co-crystal structure bound to DNA, revealing a direct interaction between the two factors and induced DNA bending. Complex interdependent binding mechanisms reveal tightly regulated TF genomic distribution and define a combinatorial logic for heterotypic TF regulation of differentiation.

Overexpression of Tbx20 in Adult Cardiomyocytes Promotes Proliferation and Improves Cardiac Function After Myocardial Infarction

BACKGROUND

Adult mammalian cardiomyocytes (CMs) have the potential to proliferate, but this is not sufficient to generate adequate CMs after myocardial infarction (MI). The transcription factor Tbx20 is required for CM proliferation during development and adult CM homeostasis. The ability of Tbx20 overexpression (Tbx20(OE)) to promote adult CM proliferation and to improve cardiac function after MI was examined.

METHODS AND RESULTS

Tbx20(OE) was induced specifically in adult mouse differentiated CMs. Increased CM proliferation and fetal-like characteristics were found in Tbx20(OE) hearts compared with controls without causing pathology 4 weeks after Tbx20(OE) at baseline. Moreover, Tbx20(OE) in adult CM after MI significantly improved survival, cardiac function, and infarct size 4 weeks after MI. Improved cardiac repair, as indicated by increased CM proliferation and capillary density, was observed in the MI border zone of Tbx20(OE) hearts compared with controls. Expression of proliferation activator (cyclin D1, E1, and IGF1) and fetal contractile protein (ssTNI, βMHC) mRNA was increased whereas negative cell-cycle regulators (p21, Meis1) were decreased in Tbx20(OE) hearts compared with controls under both baseline and MI conditions. Tbx20(OE) in adult hearts activates multiple proproliferation pathways, including Akt, YAP and BMP. Interestingly, p21, Meis1, and a novel cell-cycle inhibitory gene, Btg2, are directly bound and repressed by Tbx20 with induction of proliferation in neonatal CM.

CONCLUSIONS

Tbx20(OE), specifically in adult CM, activates multiple cardiac proliferative pathways, directly represses cell-cycle inhibitory genes p21, Meis1, and Btg2, promotes adult CM proliferation; and preserves cardiac performance after MI.

Deep sequencing-generated modules demonstrate coherent expression patterns for various cardiac diseases

As sequencing technology rapidly develops, gene annotations have also become increasingly sophisticated with incorporation of information regarding the temporal-spatial context of alternative splicing patterns, developmental stages, and tissue specificity. The present study aimed to identify the heart-enriched genes based on next-generation sequencing data and to investigate the gene modules demonstrating coherent expression patterns for various cardiac disease-related perturbations. Seven gene modules, including 382 heart-enriched genes, were identified. At least two modules containing differentially expressed genes were experimentally confirmed to be highly sensitive to various cardiac diseases. Transcription factors regulating the gene modules were then analyzed based on knowledgebase information; the expression of eight transcription factors changed significantly during pressure-overload cardiac hypertrophy, suggesting possible regulation of the modules by the identified transcription factors. Collectively, our results contribute to the classification of heart-enriched genes and their modules and would aid in identification of the transcription factors involved in cardiac pathogenesis in the future.

Sequential Binding of MEIS1 and NKX2-5 on the Popdc2 Gene: A Mechanism for Spatiotemporal Regulation of Enhancers during Cardiogenesis

The homeobox transcription factors NKX2-5 and MEIS1 are essential for vertebrate heart development and normal physiology of the adult heart. We show that, during cardiac differentiation, the two transcription factors have partially overlapping expression patterns, with the result that as cardiac progenitors from the anterior heart field differentiate and migrate into the cardiac outflow tract, they sequentially experience high levels of MEIS1 and then increasing levels of NKX2-5. Using the Popdc2 gene as an example, we also show that a significant proportion of target genes for NKX2-5 contain a binding motif recognized by NKX2-5, which overlaps with a binding site for MEIS1. Binding of the two factors to such overlapping sites is mutually exclusive, and this provides a simple regulatory mechanism for spatial and temporal synchronization of a common pool of targets between NKX2-5 and MEIS1.

•

NKX2-5 shares a DNA-binding site with MEIS1

•

MEIS1 and NKX2-5 successively bind a Popdc2 enhancer

•

Successive binding by MEIS1 and NKX2-5 is a general mechanism of regulation

•

NKX2-5 represses fast troponin isoforms in the atria

Negative regulation of ABA signaling by WRKY33 is critical for Arabidopsis immunity towards Botrytis cinerea 2100

The Arabidopsis mutant wrky33 is highly susceptible to Botrytis cinerea. We identified &gt;1680 Botrytis-induced WRKY33 binding sites associated with 1576 Arabidopsis genes. Transcriptional profiling defined 318 functional direct target genes at 14 hr post inoculation. Comparative analyses revealed that WRKY33 possesses dual functionality acting either as a repressor or as an activator in a promoter-context dependent manner. We confirmed known WRKY33 targets involved in hormone signaling and phytoalexin biosynthesis, but also uncovered a novel negative role of abscisic acid (ABA) in resistance towards B. cinerea 2100. The ABA biosynthesis genes NCED3 and NCED5 were identified as direct targets required for WRKY33-mediated resistance. Loss-of-WRKY33 function resulted in elevated ABA levels and genetic studies confirmed that WRKY33 acts upstream of NCED3/NCED5 to negatively regulate ABA biosynthesis. This study provides the first detailed view of the genome-wide contribution of a specific plant transcription factor in modulating the transcriptional network associated with plant immunity.

Negative regulation of ABA signaling by WRKY33 is critical for Arabidopsis immunity towards Botrytis cinerea 2100

The Arabidopsis mutant wrky33 is highly susceptible to Botrytis cinerea. We identified &gt;1680 Botrytis-induced WRKY33 binding sites associated with 1576 Arabidopsis genes. Transcriptional profiling defined 318 functional direct target genes at 14 hr post inoculation. Comparative analyses revealed that WRKY33 possesses dual functionality acting either as a repressor or as an activator in a promoter-context dependent manner. We confirmed known WRKY33 targets involved in hormone signaling and phytoalexin biosynthesis, but also uncovered a novel negative role of abscisic acid (ABA) in resistance towards B. cinerea 2100. The ABA biosynthesis genes NCED3 and NCED5 were identified as direct targets required for WRKY33-mediated resistance. Loss-of-WRKY33 function resulted in elevated ABA levels and genetic studies confirmed that WRKY33 acts upstream of NCED3/NCED5 to negatively regulate ABA biosynthesis. This study provides the first detailed view of the genome-wide contribution of a specific plant transcription factor in modulating the transcriptional network associated with plant immunity.

Negative regulation of ABA signaling by WRKY33 is critical for Arabidopsis immunity towards Botrytis cinerea 2100

The Arabidopsis mutant wrky33 is highly susceptible to Botrytis cinerea. We identified &gt;1680 Botrytis-induced WRKY33 binding sites associated with 1576 Arabidopsis genes. Transcriptional profiling defined 318 functional direct target genes at 14 hr post inoculation. Comparative analyses revealed that WRKY33 possesses dual functionality acting either as a repressor or as an activator in a promoter-context dependent manner. We confirmed known WRKY33 targets involved in hormone signaling and phytoalexin biosynthesis, but also uncovered a novel negative role of abscisic acid (ABA) in resistance towards B. cinerea 2100. The ABA biosynthesis genes NCED3 and NCED5 were identified as direct targets required for WRKY33-mediated resistance. Loss-of-WRKY33 function resulted in elevated ABA levels and genetic studies confirmed that WRKY33 acts upstream of NCED3/NCED5 to negatively regulate ABA biosynthesis. This study provides the first detailed view of the genome-wide contribution of a specific plant transcription factor in modulating the transcriptional network associated with plant immunity.

VISIONET: intuitive visualisation of overlapping transcription factor networks, with applications in cardiogenic gene discovery

BACKGROUND

Existing de novo software platforms have largely overlooked a valuable resource, the expertise of the intended biologist users. Typical data representations such as long gene lists, or highly dense and overlapping transcription factor networks often hinder biologists from relating these results to their expertise.

RESULTS

VISIONET, a streamlined visualisation tool built from experimental needs, enables biologists to transform large and dense overlapping transcription factor networks into sparse human-readable graphs via numerically filtering. The VISIONET interface allows users without a computing background to interactively explore and filter their data, and empowers them to apply their specialist knowledge on far more complex and substantial data sets than is currently possible. Applying VISIONET to the Tbx20-Gata4 transcription factor network led to the discovery and validation of Aldh1a2, an essential developmental gene associated with various important cardiac disorders, as a healthy adult cardiac fibroblast gene co-regulated by cardiogenic transcription factors Gata4 and Tbx20.

CONCLUSIONS

We demonstrate with experimental validations the utility of VISIONET for expertise-driven gene discovery that opens new experimental directions that would not otherwise have been identified.

Novel mutations in the transcriptional activator domain of the human TBX20 in patients with atrial septal defect

Background. The relevance of TBX20 gene in heart development has been demonstrated in many animal models, but there are few works that try to elucidate the effect of TBX20 mutations in human congenital heart diseases. In these studies, all missense mutations associated with atrial septal defect (ASD) were found in the DNA-binding T-box domain, none in the transcriptional activator domain. Methods. We search for TBX20 mutations in a group of patients with ASD or ventricular septal defect (VSD) using the High Resolution Melting (HRM) method and DNA sequencing. Results. We report three missense mutations (Y309D, T370O, and M395R) within the transcriptional activator domain of human TBX20 that were associated with ASD. Conclusions. This is the first association of TBX20 transcriptional activator domain missense mutations with ASD. These findings could have implications for diagnosis, genetic screening, and patient follow-up.

Right ventricular long noncoding RNA expression in human heart failure

The expression of long noncoding RNAs (lncRNAs) in human heart failure (HF) has not been widely studied. Using RNA sequencing (RNA-Seq), we compared lncRNA expression in 22 explanted human HF hearts with lncRNA expression in 5 unused donor human hearts. We used Cufflinks to identify isoforms and DESeq to identify differentially expressed genes. We identified the noncoding RNAs by cross-reference to Ensembl release 73 (Genome Reference Consortium human genome build 37) and explored possible functional roles using a variety of online tools. In HF hearts, RNA-Seq identified 84,793 total messenger RNA coding and noncoding different transcripts, including 13,019 protein-coding genes, 2,085 total lncRNA genes, and 1,064 pseudogenes. By Ensembl noncoding RNA categories, there were 48 lncRNAs, 27 pseudogenes, and 30 antisense RNAs for a total of 105 differentially expressed lncRNAs in HF hearts. Compared with donor hearts, HF hearts exhibited differential expression of 7.7% of protein-coding genes, 3.7% of lncRNAs (including pseudogenes), and 2.5% of pseudogenes. There were not consistent correlations between antisense lncRNAs and parent genes and between pseudogenes and parent genes, implying differential regulation of expression. Exploratory in silico functional analyses using online tools suggested a variety of possible lncRNA regulatory roles. By providing a comprehensive profile of right ventricular polyadenylated messenger RNA transcriptome in HF, RNA-Seq provides an inventory of differentially expressed lncRNAs, including antisense transcripts and pseudogenes, for future mechanistic study.

Right ventricular myocardial biomarkers in human heart failure

BACKGROUND

Right ventricular (RV) dysfunction contributes to mortality in chronic heart failure (HF). However, the molecular mechanisms of RV failure remain poorly understood, and RV myocardial biomarkers have yet to be developed.

METHODS AND RESULTS

We performed RNA sequencing (RNA-seq) on 22 explanted human HF RVs and 5 unused donor human heart RVs (DON RV) and compared results to those recently reported from 16 explanted human LVs We used Bowtie-Tophat for transcript alignment and transcriptome assembly, DESeq for identification of differentially expressed genes (DEGs) and Ingenuity for exploration of gene ontologies. In the HF RV, RNA-seq identified 130,790 total RNA transcripts including 13,272 protein coding genes, 10,831 long non-coding RNA genes and 8,605 pseudogenes. There were 800-1000 DEGs between DON and HF RV comparison groups with differences concentrated in cytoskeletal, basement membrane, extracellular matrix (ECM), inflammatory mediator, hemostasis, membrane transport and transcription factor genes, lncRNAs and pseudogenes. In an unbiased approach, the top 10 DEGs SERPINA3, SERPINA5, LCN6, LCN10, STEAP4, AKR1C1, STAC2, SPARCL1, VSIG4 and F8 exhibited no overlap in read counts between DON and HF RVs, high sensitivities, specificities, predictive values and areas under the receiver operating characteristic curves. STEAP4, SPARCL1 and VSIG4 were differentially expressed between RVs and LVs, supporting their roles as RV-specific myocardial biomarkers.

CONCLUSIONS

Unbiased, comprehensive profiling of the RV transcriptome by RNA-seq suggests structural changes and abnormalities in inflammatory processes and yields specific, novel HF RV vs HF LV myocardial biomarkers not previously identified by more limited transcriptome profiling approaches.

Modular evolution of DNA-binding preference of a Tbrain transcription factor provides a mechanism for modifying gene regulatory networks

Gene regulatory networks (GRNs) describe the progression of transcriptional states that take a single-celled zygote to a multicellular organism. It is well documented that GRNs can evolve extensively through mutations to cis-regulatory modules (CRMs). Transcription factor proteins that bind these CRMs may also evolve to produce novelty. Coding changes are considered to be rarer, however, because transcription factors are multifunctional and hence are more constrained to evolve in ways that will not produce widespread detrimental effects. Recent technological advances have unearthed a surprising variation in DNA-binding abilities, such that individual transcription factors may recognize both a preferred primary motif and an additional secondary motif. This provides a source of modularity in function. Here, we demonstrate that orthologous transcription factors can also evolve a changed preference for a secondary binding motif, thereby offering an unexplored mechanism for GRN evolution. Using protein-binding microarray, surface plasmon resonance, and in vivo reporter assays, we demonstrate an important difference in DNA-binding preference between Tbrain protein orthologs in two species of echinoderms, the sea star, Patiria miniata, and the sea urchin, Strongylocentrotus purpuratus. Although both orthologs recognize the same primary motif, only the sea star Tbr also has a secondary binding motif. Our in vivo assays demonstrate that this difference may allow for greater evolutionary change in timing of regulatory control. This uncovers a layer of transcription factor binding divergence that could exist for many pairs of orthologs. We hypothesize that this divergence provides modularity that allows orthologous transcription factors to evolve novel roles in GRNs through modification of binding to secondary sites.

Functional networks of nucleocytoplasmic transport-related genes differentiate ischemic and dilated cardiomyopathies. A new therapeutic opportunity

Heart failure provokes alterations in the expression of nucleocytoplasmic transport-related genes. To elucidate the nucleocytoplasmic transport-linked functional network underlying the two major causes of heart failure, ischemic cardiomyopathy (ICM) and dilated cardiomyopathy (DCM), we examined global transcriptome profiles of left ventricular myocardium tissue samples from 31 patients (ICM, n = 10; DCM, n = 13) undergoing heart transplantation and control donors (CNT, n = 8) using RNA-Sequencing and GeneMANIA. Comparative profiling of ICM versus control and DCM versus control showed 1081 and 2440 differentially expressed genes, respectively (>1.29-fold; P<0.05). GeneMANIA revealed differentially regulated functional networks specific to ICM and DCM. In comparison with CNT, differential expression was seen in 9 and 12 nucleocytoplasmic transport-related genes in ICM and DCM groups, respectively. DDX3X, KPNA2, and PTK2B were related to ICM, while SMURF2, NUP153, IPO5, RANBP3, NOXA1, and RHOJ were involved in DCM pathogenesis. Furthermore, the two pathologies shared 6 altered genes: XPO1, ARL4, NFKB2, FHL3, RANBP2, and RHOU showing an identical trend in expression in both ICM and DCM. Notably, the core of the derived functional networks composed of nucleocytoplasmic transport-related genes (XPO1, RANBP2, NUP153, IPO5, KPNA2, and RANBP3) branched into several pathways with downregulated genes. Moreover, we identified genes whose expression levels correlated with left ventricular mass index and left ventricular function parameters in HF patients. Collectively, our study provides a clear distinction between the two pathologies at the transcriptome level and opens up new possibilities to search for appropriate therapeutic targets for ICM and DCM.

It is not the parts, but how they interact that determines the behaviour of circadian rhythms across scales and organisms

Biological rhythms, generated by feedback loops containing interacting genes, proteins and/or cells, time physiological processes in many organisms. While many of the components of the systems that generate biological rhythms have been identified, much less is known about the details of their interactions. Using examples from the circadian (daily) clock in three organisms, Neurospora, Drosophila and mouse, we show, with mathematical models of varying complexity, how interactions among (i) promoter sites, (ii) proteins forming complexes, and (iii) cells can have a drastic effect on timekeeping. Inspired by the identification of many transcription factors, for example as involved in the Neurospora circadian clock, that can both activate and repress, we show how these multiple actions can cause complex oscillatory patterns in a transcription–translation feedback loop (TTFL). Inspired by the timekeeping complex formed by the NMO–PER–TIM–SGG complex that regulates the negative TTFL in the Drosophila circadian clock, we show how the mechanism of complex formation can determine the prevalence of oscillations in a TTFL. Finally, we note that most mathematical models of intracellular clocks model a single cell, but compare with experimental data from collections of cells. We find that refitting the most detailed model of the mammalian circadian clock, so that the coupling between cells matches experimental data, yields different dynamics and makes an interesting prediction that also matches experimental data: individual cells are bistable, and network coupling removes this bistability and causes the network to be more robust to external perturbations. Taken together, we propose that the interactions between components in biological timekeeping systems are carefully tuned towards proper function. We also show how timekeeping can be controlled by novel mechanisms at different levels of organization.

A common genetic variant within SCN10A modulates cardiac SCN5A expression

Variants in SCN10A, which encodes a voltage-gated sodium channel, are associated with alterations of cardiac conduction parameters and the cardiac rhythm disorder Brugada syndrome; however, it is unclear how SCN10A variants promote dysfunctional cardiac conduction. Here we showed by high-resolution 4C-seq analysis of the Scn10a-Scn5a locus in murine heart tissue that a cardiac enhancer located in Scn10a, encompassing SCN10A functional variant rs6801957, interacts with the promoter of Scn5a, a sodium channel-encoding gene that is critical for cardiac conduction. We observed that SCN5A transcript levels were several orders of magnitude higher than SCN10A transcript levels in both adult human and mouse heart tissue. Analysis of BAC transgenic mouse strains harboring an engineered deletion of the enhancer within Scn10a revealed that the enhancer was essential for Scn5a expression in cardiac tissue. Furthermore, the common SCN10A variant rs6801957 modulated Scn5a expression in the heart. In humans, the SCN10A variant rs6801957, which correlated with slowed conduction, was associated with reduced SCN5A expression. These observations establish a genomic mechanism for how a common genetic variation at SCN10A influences cardiac physiology and predisposes to arrhythmia.

Large-scale quality analysis of published ChIP-seq data

ChIP-seq has become the primary method for identifying in vivo protein-DNA interactions on a genome-wide scale, with nearly 800 publications involving the technique appearing in PubMed as of December 2012. Individually and in aggregate, these data are an important and information-rich resource. However, uncertainties about data quality confound their use by the wider research community. Recently, the Encyclopedia of DNA Elements (ENCODE) project developed and applied metrics to objectively measure ChIP-seq data quality. The ENCODE quality analysis was useful for flagging datasets for closer inspection, eliminating or replacing poor data, and for driving changes in experimental pipelines. There had been no similarly systematic quality analysis of the large and disparate body of published ChIP-seq profiles. Here, we report a uniform analysis of vertebrate transcription factor ChIP-seq datasets in the Gene Expression Omnibus (GEO) repository as of April 1, 2012. The majority (55%) of datasets scored as being highly successful, but a substantial minority (20%) were of apparently poor quality, and another ∼25% were of intermediate quality. We discuss how different uses of ChIP-seq data are affected by specific aspects of data quality, and we highlight exceptional instances for which the metric values should not be taken at face value. Unexpectedly, we discovered that a significant subset of control datasets (i.e., no immunoprecipitation and mock immunoprecipitation samples) display an enrichment structure similar to successful ChIP-seq data. This can, in turn, affect peak calling and data interpretation. Published datasets identified here as high-quality comprise a large group that users can draw on for large-scale integrated analysis. In the future, ChIP-seq quality assessment similar to that used here could guide experimentalists at early stages in a study, provide useful input in the publication process, and be used to stratify ChIP-seq data for different community-wide uses.

Beyond the ENCODE project: using genomics and epigenomics strategies to study enhancer evolution

The complex expression patterns observed for many genes are often regulated by distal transcription enhancers. Changes in the nucleotide sequences of enhancers may therefore lead to changes in gene expression, representing a central mechanism by which organisms evolve. With the development of the experimental technique of chromatin immunoprecipitation (ChIP), in which discrete regions of the genome bound by specific proteins can be identified, it is now possible to identify transcription factor binding events (putative cis-regulatory elements) in entire genomes. Comparing protein-DNA binding maps allows us, for the first time, to attempt to identify regulatory differences and infer global patterns of change in gene expression across species. Here, we review studies that used genome-wide ChIP to study the evolution of enhancers. The trend is one of high divergence of cis-regulatory elements between species, possibly compensated by extensive creation and loss of regulatory elements and rewiring of their target genes. We speculate on the meaning of the differences observed and discuss that although ChIP experiments identify the biochemical event of protein-DNA interaction, it cannot determine whether the event results in a biological function, and therefore more studies are required to establish the effect of divergence of binding events on species-specific gene expression.

A Gro/TLE-NuRD corepressor complex facilitates Tbx20-dependent transcriptional repression

The cardiac transcription factor Tbx20 has a critical role in the proper morphogenetic development of the vertebrate heart, and its misregulation has been implicated in human congenital heart disease. Although it is established that Tbx20 exerts its function in the embryonic heart through positive and negative regulation of distinct gene programs, it is unclear how Tbx20 mediates proper transcriptional regulation of its target genes. Here, using a combinatorial proteomic and bioinformatic approach, we present the first characterization of Tbx20 transcriptional protein complexes. We have systematically investigated Tbx20 protein-protein interactions by immunoaffinity purification of tagged Tbx20 followed by proteomic analysis using GeLC-MS/MS, gene ontology classification, and functional network analysis. We demonstrate that Tbx20 is associated with a chromatin remodeling network composed of TLE/Groucho corepressors, members of the Nucleosome Remodeling and Deacetylase (NuRD) complex, the chromatin remodeling ATPases RUVBL1/RUVBL2, and the T-box repressor Tbx18. We determined that the interaction with TLE corepressors is mediated via an eh1 binding motif in Tbx20. Moreover, we demonstrated that ablation of this motif results in a failure to properly assemble the repression network and disrupts Tbx20 function in vivo. Importantly, we validated Tbx20-TLE interactions in the mouse embryonic heart, and identified developmental genes regulated by Tbx20-TLE binding, thereby confirming a primary role for a Tbx20-TLE repressor complex in embryonic heart development. Together, these studies suggest a model in which Tbx20 associates with a Gro/TLE-NuRD repressor complex to prevent inappropriate gene activation within the forming heart.