The SIX Family of Transcription Factors: Common Themes Integrating Developmental and Cancer Biology

Tumors and their microenvironments: Learning from pediatric brain pathologies

Early clues to tumors and their microenvironments come from embryonic development. Here we review the literature and consider whether the embryonic brain and its pathologies can serve as a better model. Among embryonic organs, the brain is the most heterogenous and complex, with multiple lineages leading to wide spectrum of cell states and types. Its dysregulation promotes neurodevelopmental brain pathologies and pediatric tumors. Embryonic brain pathologies point to the crucial importance of spatial heterogeneity over time, akin to the tumor microenvironment. Tumors dedifferentiate through genetic mutations and epigenetic modulations; embryonic brains differentiate through epigenetic modulations. Our innovative review proposes learning developmental brain pathologies to target tumor evolution-and vice versa. We describe ways through which tumor pharmacology can learn from embryonic brains and their pathologies, and how learning tumor, and its microenvironment, can benefit targeting neurodevelopmental pathologies. Examples include pediatric low-grade versus high-grade brain tumors as in rhabdomyosarcomas and gliomas.

A population-based analysis of the molecular landscape of glioma in adolescents and young adults reveals insights into gliomagenesis

Gliomas are a major cause of cancer-related deaths in adolescents and young adults (AYAs; ages 15-39 years). Different molecular alterations drive gliomas in children and adults, leading to distinct biology and clinical consequences, but the implications of pediatric- versus adult-type alterations in AYAs are unknown. Our population-based analysis of 1,456 clinically and molecularly characterized gliomas in patients aged 0-39 years addresses this gap. Pediatric-type alterations were found in 31% of AYA gliomas and conferred superior outcomes compared to adult-type alterations. AYA low-grade gliomas with specific RAS-MAPK alterations exhibited senescence, tended to arise in different locations and were associated with superior outcomes compared to gliomas in children, suggesting different cellular origins. Hemispheric IDH-mutant, BRAF p.V600E and FGFR-altered gliomas were associated with the risk of malignant transformation, having worse outcomes with increased age. These insights into gliomagenesis may provide a rationale for earlier intervention for certain tumors to disrupt the typical behavior, leading to improved outcomes.

Sine oculis homeobox homolog family function in gastrointestinal cancer: Progression and comprehensive analysis

The sine oculis homeobox homolog (SIX) family, a group of transcription factors characterized by a conserved DNA-binding homology domain, plays a critical role in orchestrating embryonic development and organogenesis across various organisms, including humans. Comprising six distinct members, from SIX1 to SIX6, each member contributes uniquely to the development and differentiation of diverse tissues and organs, underscoring the versatility of the SIX family. Dysregulation or mutations in SIX genes have been implicated in a spectrum of developmental disorders, as well as in tumor initiation and progression, highlighting their pivotal role in maintaining normal developmental trajectories and cellular functions. Efforts to target the transcriptional complex of the SIX gene family have emerged as a promising strategy to inhibit tumor development. While the development of inhibitors targeting this gene family is still in its early stages, the significant potential of such interventions holds promise for future therapeutic advances. Therefore, this review aimed to comprehensively explore the advancements in understanding the SIX family within gastrointestinal cancers, focusing on its critical role in normal organ development and its implications in gastrointestinal cancers, including gastric, pancreatic, colorectal cancer, and hepatocellular carcinomas. In conclusion, this review deepened the understanding of the functional roles of the SIX family and explored the potential of utilizing this gene family for the diagnosis, prognosis, and treatment of gastrointestinal cancers.

Familial Oculoauriculovertebral Spectrum: A Genomic Investigation of Autosomal Dominant Inheritance

OBJECTIVE

Oculoauriculovertebral spectrum (OAVS) encompasses abnormalities on derivatives from the first and second pharyngeal arches including macrostomia, hemifacial microsomia, micrognathia, preauricular tags, ocular, and vertebral anomalies. We present genetic findings on a 3-generation family affected with macrostomia, preauricular tags and ptosis following an autosomal dominant pattern.

DESIGN

We generated whole-genome sequencing data for the proband, affected father, and unaffected paternal grandmother followed by Sanger sequencing on 23 family members for the top candidate gene mutations. We performed parent and sibling-based transmission disequilibrium tests (TDTs) and burden analysis via a penalized linear mixed model, for segregation and mutation burden, respectively. Next, via bioinformatic tools we predicted protein function, mutation pathogenicity, and pathway enrichment to investigate the biological relevance of mutations identified.

RESULTS

Rare missense mutations in SIX1, KDR/VEGFR2, and PDGFRA showed the best segregation with the OAVS phenotypes in this family. When considering any of the 3 OAVS phenotypes as an outcome, SIX1 had the strongest associations in parent-TDTs and sib-TDTs (P = 0.025, P = 0.052) (unadjusted P-values). Burden analysis identified SIX1 (RC = 0.87) and PDGFRA (RC = 0.98) strongly associated with OAVS severity. Using phenotype-specific outcomes, sib-TDTs identified SIX1 with uni- or bilateral ptosis (P = 0.049) and ear tags (P = 0.01), and PDGFRA and KDR/VEGFR2 with ear tags (both P < 0.01).

CONCLUSION

SIX1, PDGFRA, and KDR/VEGFR2 are strongly associated to OAVS phenotypes. SIX1 has been previously associated with OAVS ear malformations and is co-expressed with EYA1 during ear development. Efforts to strengthen the genotype-phenotype co-relation underlying the OAVS are key to discover etiology, family counseling, and prevention.

SIX1 aggravates the progression of spinal cord injury in mice by promoting M1 polarization of microglia

Inflammation aggravates secondary damage following spinal cord injury (SCI). M1 microglia induce inflammation and exert neurotoxic effects, whereas M2 microglia exert anti-inflammatory and neuroprotective effects. The sine oculis homeobox (SIX) gene family consists of six members, including sine oculis homeobox homolog 1 (SIX1)-SIX6. SIX1 is expressed in microglia and promotes inflammation. This study aimed to evaluate the role and underlying mechanisms of SIX1 in microglia polarization in vitro (LPS-treated mouse microglia; BV2 cells) and in vivo (a mouse model of SCI). SIX1 expression was increased in the microglia of mice with SCI. SIX1 was positively correlated with the M1 microglia marker inducible nitric oxide synthase (iNOS) and negatively correlated with the M2 microglia marker arginase 1 (Arg1) in mice with SCI. Knockdown of SIX1 promoted functional recovery by enhancing M2 microglia polarization in mice with SCI. The transcription, expression, and activity of enhancer of zeste homolog 2 (EZH2) were decreased in LPS-stimulated BV2 cells. Downregulation of EZH2 promoted SIX1 expression in LPS-treated BV2 cells by inhibiting the methylation of the SIX1 promoter. SIX1 enhanced the transcription of vascular endothelial growth factor-C (VEGF-C) in LPS-stimulated BV2 cells with downregulated EZH2. VEGF-C promoted M1 polarization and inhibited M2 polarization in BV2 cells by binding to vascular endothelial growth factor receptor 3 (VEGFR3). Overall, the results suggest that SIX1 promotes M1 polarization of microglia following SCI by upregulating the VEGF-C/VEGFR3 axis, whereas the blockade of SIX1 can improve the recovery of locomotor function following SCI, demonstrating a novel strategy for the treatment of SCI.

SIX transcription factors are necessary for the activation of DUX4 expression in facioscapulohumeral muscular dystrophy

BACKGROUND

Facioscapulohumeral muscular dystrophy (FSHD) is a common and progressive muscle wasting disease that is characterized by muscle weakness often first noticed in the face, the shoulder girdle and upper arms before progressing to the lower limb muscles. FSHD is caused by the misexpression of the Double Homeobox 4 (DUX4) transcription factor in skeletal muscle. While epigenetic derepression of D4Z4 macrosatellite repeats underlies DUX4 misexpression, our understanding of the complex transcriptional activation of DUX4 is incomplete.

METHODS

To identify potential DUX4-regulatory factors, we used small interfering RNAs (siRNAs) to knockdown SIX family transcription factors (SIX1, 2, 4, 5) in patient-derived FSHD1 and FSHD2 myoblasts that were differentiated to form multinucleated myotubes. Quantitative real-time polymerase chain reaction was used to measure changes in DUX4 mRNA, DUX4 target gene expression and myogenic markers. Staining for SIX1 and SIX2 with specific antibodies was performed in FSHD myoblasts and myotubes. To assess reciprocal effects of DUX4 on SIX1, 2, and 4 expression, we utilized a doxycycline-inducible DUX4 myoblast cell line.

RESULT

We show that SIX1, 2 and 4 transcription factors, regulators of embryonic development, muscle differentiation, regeneration and homeostasis, are necessary for myogenic differentiation-dependent DUX4 expression in FSHD muscle cells. Using siRNA, we demonstrate SIX1, SIX2, and SIX4 to be critical factors involved in the induction of DUX4 transcription in differentiating FSHD myotubes in vitro. siRNA dual knockdown of SIX1 and SIX2 resulted in a ~ 98% decrease of DUX4 and DUX4 target genes, suggesting that SIX1 and SIX2 are the most critical in promoting DUX4 expression. Importantly, we show that DUX4 downregulates SIX RNA levels, suggesting negative feedback regulation.

CONCLUSIONS

In this study, we identified a family of developmental regulators that promote aberrant DUX4 expression in FSHD1 and FSHD2 differentiating muscle cells. Our findings highlight the critical involvement of SIX transcription factors (SIX1, 2, 4) in the pathogenesis of FSHD by serving as necessary factors that function in the promotion of DUX4 expression following epigenetic derepression of the D4Z4 repeats.

EBV infection alters DNA methylation in primary human colon cells: A path to inflammation and carcinogenesis?

Epstein-Barr Virus (EBV) is associated with several types of human cancers, and changes in DNA methylation are reported to contribute to viral-driven carcinogenesis, particularly in cancers of epithelial origin. In a previous study, we demonstrated that EBV infects human primary colonic cells (HCoEpC) and replicates within these cells, leading to pro-inflammatory and pro-tumorigenic effects. Notably, these effects were mostly prevented by inhibiting viral replication with PAA. Interestingly, the EBV-induced effects correlated with the upregulation of DNMT1 and were counteracted by pretreating cells with 5-AZA, suggesting a role for DNA hypermethylation. Building on this background, the current study investigates the methylation changes induced by EBV infection in HCoEpC, both in the presence and absence of PAA, or ERK1/2 and STAT3 inhibitors, pathways known to be activated by EBV and involved in the dysregulation of methylation in tumor cells. The genome-wide methylation analysis conducted in this study allowed us to identify several biological processes and genes affected by these epigenetic changes, providing insights into the possible underlying mechanisms leading to the pathological effects induced by EBV. Specifically, we found that the virus induced significant methylation changes, with hypermethylation being more prevalent than hypomethylation. Several genes involved in embryogenesis, carcinogenesis, and inflammation were affected.

Candidate SNP Markers Significantly Altering the Affinity of the TATA-Binding Protein for the Promoters of Human Genes Associated with Primary Open-Angle Glaucoma

Primary open-angle glaucoma (POAG) is the most common form of glaucoma. This condition leads to optic nerve degeneration and eventually to blindness. Tobacco smoking, alcohol consumption, fast-food diets, obesity, heavy weight lifting, high-intensity physical exercises, and many other bad habits are lifestyle-related risk factors for POAG. By contrast, moderate-intensity aerobic exercise and the Mediterranean diet can alleviate POAG. In this work, we for the first time estimated the phylostratigraphic age indices (PAIs) of all 153 POAG-related human genes in the NCBI Gene Database. This allowed us to separate them into two groups: POAG-related genes that appeared before and after the phylum Chordata, that is, ophthalmologically speaking, before and after the camera-type eye evolved. Next, in the POAG-related genes' promoters, we in silico predicted all 3835 candidate SNP markers that significantly change the TATA-binding protein (TBP) affinity for these promoters and, through this molecular mechanism, the expression levels of these genes. Finally, we verified our results against five independent web services-PANTHER, DAVID, STRING, MetaScape, and GeneMANIA-as well as the ClinVar database. It was concluded that POAG is likely to be a symptom of the human self-domestication syndrome, a downside of being civilized.

Targeting non-coding RNAs as a promising biomarker in peritoneal metastasis: Background, mechanism, and therapeutic approach

Peritoneal metastasis (PM) pathophysiology is complex and not fully understood. PM, originating from gastrointestinal (GI) cancer, is a condition that significantly worsens patient prognosis due to its complex nature and limited treatment options. The non-coding RNAs (ncRNAs) have been shown to play pivotal roles in cancer biology, influencing tumorigenesis, progression, metastasis, and therapeutic resistance. Increasing evidence has demonstrated the regulatory functions of different classes of ncRNAs, including microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) in PM. Identifying biomarkers for early detection of PM is a crucial step towards improving patient outcomes, and how ncRNA profiles correlate with survival rates, response to therapy, and recurrence risks have raised much attention in recent years. Additionally, exploring innovative therapeutic approaches utilizing ncRNAs, such as targeted therapy and gene silencing, may offer new horizons in treating this dire condition. Recent advances in systemic treatments and the development of novel loco-regional therapies have opened doors to multimodal treatment approaches. Radical surgeries combined with hyperthermic intraperitoneal chemotherapy (HIPEC) have shown promising results, leading to extended patient survival. Current research is focused on the molecular characterization of PM, which is crucial for early detection and developing future therapeutic strategies. By summarizing the latest findings, this study underscores the transformative potential of ncRNAs in enhancing the diagnosis, prognosis, and treatment of PM in GI cancer, paving the way for more personalized and effective clinical strategies.

The homeodomain transcription factor Six3 regulates hypothalamic Pomc expression and its absence from POMC neurons induces hyperphagia and mild obesity in male mice

OBJECTIVE

Proopiomelanocortin (POMC) neurons release potent anorexigenic neuropeptides, which suppress food intake and enhance energy expenditure via melanocortin receptors. Although the importance of central melanocortin in physiological regulation is well established, the underlying genetic mechanisms that define the functional identity of melanocortin neurons and maintain hypothalamic Pomc expression remain to be fully determined. In this study, we investigate the functional significance of Six3, a transcriptional regulator notably expressed in embryonic and adult mouse POMC neurons, in the regulation of hypothalamic Pomc expression and downstream physiological consequences.

METHODS

We first evaluated the expression of Six3 in the developing and adult hypothalamus by double fluorescence in situ hybridization. Next, we assessed POMC immunoreactivity in mutant mice selectively lacking Six3 from Pomc-expressing neurons and quantified Pomc mRNA levels in a tamoxifen-inducible Six3 knockout mouse model activated at embryonic E9.5 days. We also determined glucose and insulin sensitivity, daily food intake, body composition and body weight in adult male and female mice lacking Six3 specifically from POMC neurons. Lastly, we assessed the physiological consequences of ablating Six3 from POMC neurons in adult mice.

RESULTS

Six3 and Pomc were co-expressed in mouse hypothalamic neurons during development and adulthood. Mouse embryos deficient in Six3 showed reduced Pomc expression in the developing hypothalamus. Targeted deletion of Six3 specifically from POMC neurons resulted in decreased hypothalamic Pomc expression, increased daily food intake, enhanced glucose sensitivity and mild obesity in male but not in female mice. Finally, conditional removal of Six3 from POMC neurons in adult mice led to a reduction in hypothalamic POMC immunoreactivity with no significant effects in body weight or food intake.

CONCLUSIONS

Altogether, our results demonstrate that Six3 plays an essential role in the early establishment of POMC neuron identity and the maintenance of physiological levels of hypothalamic Pomc expression. In addition, our study demonstrates that the functional significance of Six3 expression in POMC neurons is sexually dimorphic and age-dependent.

The Eyes Absent family: At the intersection of DNA repair, mitosis, and replication

The Eyes Absent family (EYA1-4) are a group of dual function proteins that act as both tyrosine phosphatases and transcriptional co-activators. EYA proteins play a vital role in development, but are also aberrantly overexpressed in cancers, where they often confer an oncogenic effect. Precisely how the EYAs impact cell biology is of growing interest, fuelled by the therapeutic potential of an expanding repertoire of EYA inhibitors. Recent functional studies suggest that the EYAs are important players in the regulation of genome maintenance pathways including DNA repair, mitosis, and DNA replication. While the characterized molecular mechanisms have predominantly been ascribed to EYA phosphatase activities, EYA co-transcriptional activity has also been found to impact the expression of genes that support these pathways. This indicates functional convergence of EYA phosphatase and co-transcriptional activities, highlighting the emerging importance of the EYA protein family at the intersection of genome maintenance mechanisms. In this review, we discuss recent progress in defining EYA protein substrates and transcriptional effects, specifically in the context of genome maintenance. We then outline future directions relevant to the field and discuss the clinical utility of EYA inhibitors.

Whole genome sequencing of a family with autosomal dominant features within the oculoauriculovertebral spectrum

Background

Oculoauriculovertebral Spectrum (OAVS) encompasses abnormalities on derivatives from the first and second pharyngeal arches including macrostomia, hemifacial microsomia, micrognathia, preauricular tags, ocular and vertebral anomalies. We present genetic findings on a three-generation family affected with macrostomia, preauricular tags and uni- or bilateral ptosis following an autosomal dominant pattern.

Methods

We generated whole genome sequencing data for the proband, affected parent and unaffected paternal grandparent followed by Sanger sequencing on 23 family members for the top 10 candidate genes: KCND2, PDGFRA, CASP9, NCOA3, WNT10A, SIX1, MTF1, KDR/VEGFR2, LRRK1, and TRIM2 We performed parent and sibling-based transmission disequilibrium tests and burden analysis via a penalized linear mixed model, for segregation and mutation burden respectively. Next, via bioinformatic tools we predicted protein function, mutation pathogenicity and pathway enrichment to investigate the biological relevance of mutations identified.

Results

Rare missense mutations in SIX1, KDR/VEGFR2, and PDGFRA showed the best segregation with the OAV phenotypes in this family. When considering any of the 3 OAVS phenotypes as an outcome, SIX1 had the strongest associations in parent-TDTs and sib-TDTs (p=0.025, p=0.052) (unadjusted p-values). Burden analysis identified SIX1 (RC=0.87) and PDGFRA (RC=0.98) strongly associated with OAVS severity. Using phenotype-specific outcomes, sib-TDTs identified SIX1 with uni- or bilateral ptosis (p=0.049) and ear tags (p=0.01), and PDGFRA and KDR/VEGFR2 with ear tags (both p<0.01).

Conclusion

SIX1, PDGFRA, and KDR/VEGFR2 are strongly associated to OAVS phenotypes. SIX1 has been previously associated with OAVS ear malformations and is co-expressed with EYA1 during ear development. Efforts to strengthen the genotype-phenotype co-relation underlying the OAVS are key to discover etiology, family counseling and prevention.

Reduced Nephron Endowment in Six2-TGCtg Mice Is Due to Six3 Misexpression by Aberrant Enhancer-Promoter Interactions in the Transgene

Key Points

Aberrant enhancer–promoter interactions detected by Hi-C drive ectopic expression of Six3 in the Six2TGCtg line. Disruption of Six3 in the Six2TGCtg line restores nephron number, implicating SIX3 interference with SIX2 function in nephron progenitor cell renewal.

Background

Lifelong kidney function relies on the complement of nephrons generated during mammalian development from a mesenchymal nephron progenitor cell population. Low nephron endowment confers increased susceptibility to CKD. Reduced nephron numbers in the popular Six2TGC transgenic mouse line may be due to disruption of a regulatory gene at the integration site and/or ectopic expression of a gene(s) contained within the transgene.

Methods

Targeted locus amplification was performed to identify the integration site of the Six2TGC transgene. Genome-wide chromatin conformation capture (Hi-C) datasets were generated from nephron progenitor cells isolated from the Six2TGC +/tg mice, the Cited1 CreERT2/+ control mice, and the Six2TGC +/tg ; Tsc1 +/Flox mice that exhibited restored nephron number compared with Six2TGC +/tg mice. Modified transgenic mice lacking the C-terminal domain of Six3 were used to evaluate the mechanism of nephron number reduction in the Six2TGC +/tg mouse line.

Results

Targeted locus amplification revealed integration of the Six2TGC transgene within an intron of Cntnap5a on chr1, and Hi-C analysis mapped the precise integration of Six2TGC and Cited1 CreERT2 transgenes to chr1 and chr14, respectively. No changes in topology, accessibility, or expression were observed within the 50-megabase region centered on Cntnap5a in Six2TGC +/tg mice compared with control mice. By contrast, we identified an aberrant regulatory interaction between a Six2 distal enhancer and the Six3 promoter contained within the transgene. Increasing the Six2TGC tg to Six2 locus ratio or removing one Six2 allele in Six2TGC +/tg mice caused severe renal hypoplasia. Furthermore, clustered regularly interspaced short palindromic repeats disruption of Six3 within the transgene (Six2TGC ∆Six3CT ) restored nephron endowment to wild-type levels and abolished the stoichiometric effect.

Conclusions

These findings broadly demonstrate the utility of Hi-C data in mapping transgene integration sites and architecture. Data from genetic and biochemical studies together suggest that in Six2TGC kidneys, SIX3 interferes with SIX2 function in nephron progenitor cell renewal through its C-terminal domain.

Genetic architecture of trait variance in craniofacial morphology

The genetic architecture of trait variance has long been of interest in genetics and evolution. One of the earliest attempts to understand this architecture was presented in Lerner's Genetic Homeostasis (1954). Lerner proposed that heterozygotes should be better able to tolerate environmental perturbations because of functional differences between the alleles at a given locus, with each allele optimal for slightly different environments. This greater robustness to environmental variance, he argued, would result in smaller trait variance for heterozygotes. The evidence for Lerner's hypothesis has been inconclusive. To address this question using modern genomic methods, we mapped loci associated with differences in trait variance (vQTL) on 1,101 individuals from the F34 of an advanced intercross between LG/J and SM/J mice. We also mapped epistatic interactions for these vQTL in order to understand the influence of epistasis for the architecture of trait variance. We did not find evidence supporting Lerner's hypothesis, that heterozygotes tend to have smaller trait variances than homozygotes. We further show that the effects of most mapped loci on trait variance are produced by epistasis affecting trait means and that those epistatic effects account for about a half of the differences in genotypic-specific trait variances. Finally, we propose a model where the different interactions between the additive and dominance effects of the vQTL and their epistatic partners can explain Lerner's original observations but can also be extended to include other conditions where heterozygotes are not the least variable genotype.

The zebrafish paralog six2b is required for early proximal pronephros morphogenesis

The transcription factor Six2 plays a crucial role in maintaining self-renewing nephron progenitor cap mesenchyme (CM) during metanephric kidney development. In mouse and human, expression at single-cell resolution has detected Six2 in cells as they leave the CM pool and differentiate. The role Six2 may play in these cells as they differentiate remains unknown. Here, we took advantage of the zebrafish pronephric kidney which forms directly from intermediate mesoderm to test six2b function during pronephric tubule development and differentiation. Expression of six2b during early zebrafish development was consistent with a role in pronephros formation. Using morpholino knock-down and CRISPR/Cas9 mutagenesis, we show a functional role for six2b in the development of proximal elements of the pronephros. By 48 h post-fertilization, six2b morphants and mutants showed disrupted pronephric tubule morphogenesis. We observed a lower-than-expected frequency of phenotypes in six2b stable genetic mutants suggesting compensation. Supporting this, we detected increased expression of six2a in six2b stable mutant embryos. To further confirm six2b function, F0 crispant embryos were analyzed and displayed similar phenotypes as morphants and stable mutants. Together our data suggests a conserved role for Six2 during nephrogenesis and a role in the morphogenesis of the proximal tubule.

Reduced nephron endowment in the common Six2-TGC tg mouse line is due to Six3 misexpression by aberrant enhancer-promoter interactions in the transgene

Lifelong kidney function relies on the complement of nephrons generated during mammalian development from a mesenchymal nephron progenitor cell (NPC) population. Low nephron endowment confers increased susceptibility to chronic kidney disease. We asked whether reduced nephron numbers in the popular Six2TGC transgenic mouse line 1 was due to disruption of a regulatory gene at the integration site or to ectopic expression of a gene(s) contained within the transgene. Targeted locus amplification identified integration of the Six2TGC transgene within an intron of Cntnap5a on chr1. We generated Hi-C datasets from NPCs isolated from the Six2TGC tg/+ mice, the Cited1 CreERT2/+ control mice, and the Six2TGC tg/+ ; Tsc1 +/Flox,2 mice that exhibited restored nephron number compared with Six2TGC tg/+ mice, and mapped the precise integration of Six2TGC and Cited1 CreERT2 transgenes to chr1 and chr14, respectively. No changes in topology, accessibility, or expression were observed within the 50-megabase region centered on Cntnap5a in Six2TGC tg/+ mice compared with control mice. By contrast, we identified an aberrant regulatory interaction between a Six2 distal enhancer and the Six3 promoter contained within the transgene. Increasing the Six2TGC tg to Six2 locus ratio or removing one Six2 allele in Six2TGC tg/+ mice, caused severe renal hypoplasia. Furthermore, CRISPR disruption of Six3 within the transgene ( Six2TGC Δ Six3CT ) restored nephron endowment to wildtype levels and abolished the stoichiometric effect. Data from genetic and biochemical studies together suggest that in Six2TGC, SIX3 interferes with SIX2 function in NPC renewal through its C-terminal domain.

Significance

Using high-resolution chromatin conformation and accessibility datasets we mapped the integration site of two popular transgenes used in studies of nephron progenitor cells and kidney development. Aberrant enhancer-promoter interactions drive ectopic expression of Six3 in the Six2TGC tg line which was correlated with disruption of nephrogenesis. Disruption of Six3 within the transgene restored nephron numbers to control levels; further genetic and biochemical studies suggest that Six3 interferes with Six2 -mediated regulation of NPC renewal.

Patterns of transcription factor binding and epigenome at promoters allow interpretable predictability of multiple functions of non-coding and coding genes

Understanding the biological roles of all genes only through experimental methods is challenging. A computational approach with reliable interpretability is needed to infer the function of genes, particularly for non-coding RNAs. We have analyzed genomic features that are present across both coding and non-coding genes like transcription factor (TF) and cofactor ChIP-seq (823), histone modifications ChIP-seq (n = 621), cap analysis gene expression (CAGE) tags (n = 255), and DNase hypersensitivity profiles (n = 255) to predict ontology-based functions of genes. Our approach for gene function prediction was reliable (>90% balanced accuracy) for 486 gene-sets. PubMed abstract mining and CRISPR screens supported the inferred association of genes with biological functions, for which our method had high accuracy. Further analysis revealed that TF-binding patterns at promoters have high predictive strength for multiple functions. TF-binding patterns at the promoter add an unexplored dimension of explainable regulatory aspects of genes and their functions. Therefore, we performed a comprehensive analysis for the functional-specificity of TF-binding patterns at promoters and used them for clustering functions to reveal many latent groups of gene-sets involved in common major cellular processes. We also showed how our approach could be used to infer the functions of non-coding genes using the CRISPR screens of coding genes, which were validated using a long non-coding RNA CRISPR screen. Thus our results demonstrated the generality of our approach by using gene-sets from CRISPR screens. Overall, our approach opens an avenue for predicting the involvement of non-coding genes in various functions.

Chikungunya virus infection induces ultrastructural changes and impaired neuronal differentiation of human neurospheres

Chikungunya virus (CHIKV) is an arthropod-borne virus recently associated with large outbreaks in many parts of the world. Infection is typically manifested as a febrile and self-limited illness, characterized by joint pain and myalgia, albeit severe neurological manifestations are also reported. Although CHIKV is not recognized as a truly neurotropic virus, neurons, astrocytes, and oligodendrocytes are susceptible to infection in vitro. Here we employed a model of 3D cell culture to obtain neurospheres from ATRA/BNDF differentiated human neuroblastoma cells. We demonstrate that CHIKV is able to establish a productive infection, resulting in ultrastructural changes in cell morphology and impaired neuronal differentiation. Ultrastructural analysis of neurospheres infected with CHIKV during neuronal differentiation revealed diminished neuron dendrite formation, accumulation of viral particles associated with the plasma membrane, numerous cell vacuoles, and swollen mitochondria. Apoptotic cells were significantly increased at 72 h post-infection. Compared to Zika virus, a well-characterized neurotropic arbovirus, CHIKV infection resulted in a more discrete, albeit detectable upregulation of IL-6 levels. Finally, we found that CHIKV infection resulted in an altered profile expression, mainly downregulation, of a group of transcription factors named Hox genes. Altogether our findings highlight important features of CHIKV in the CNS, as well as the feasibility of neurospheres as robust experimental models that can support further studies for novel pharmacological interventions.

An atrial fibrillation-associated regulatory region modulates cardiac Tbx5 levels and arrhythmia susceptibility

Heart development and rhythm control are highly Tbx5 dosage-sensitive. TBX5 haploinsufficiency causes congenital conduction disorders, whereas increased expression levels of TBX5 in human heart samples has been associated with atrial fibrillation (AF). We deleted the conserved mouse orthologues of two independent AF-associated genomic regions in the Tbx5 locus, one intronic (RE(int)) and one downstream (RE(down)) of Tbx5. In both lines, we observed a modest (30%) increase of Tbx5 in the postnatal atria. To gain insight into the effects of slight dosage increase in vivo, we investigated the atrial transcriptional, epigenetic and electrophysiological properties of both lines. Increased atrial Tbx5 expression was associated with induction of genes involved in development, ion transport and conduction, with increased susceptibility to atrial arrhythmias, and increased action potential duration of atrial cardiomyocytes. We identified an AF-associated variant in the human RE(int) that increases its transcriptional activity. Expression of the AF-associated transcription factor Prrx1 was induced in Tbx5^RE(int)KO cardiomyocytes. We found that some of the transcriptional and functional changes in the atria caused by increased Tbx5 expression were normalized when reducing cardiac Prrx1 expression in Tbx5^RE(int)KO mice, indicating an interaction between these two AF genes. We conclude that modest increases in expression of dose-dependent transcription factors, caused by common regulatory variants, significantly impact on the cardiac gene regulatory network and disease susceptibility.