The Human SIX Family of Homeobox Genes

ATAC-Seq for Assaying Chromatin Accessibility Protocol Using Echinoderm Embryos

Cis-regulatory elements (CREs) and transcription factors (TFs) associated with them determine temporal and spatial domains of gene expression. Therefore, identification of these CREs and TFs is crucial to elucidating transcriptional programs across taxa. With chromatin accessibility facilitating transcription factor access to DNA, the identification of regions of open chromatin sheds light both on the function of the regulatory elements and their evolution, thus allowing the recognition of potential CREs. Buenrostro and colleagues have developed a novel method for exploring chromatin accessibility: assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), which can be used for the purpose of identifying putative CREs. This method was shown to have considerable advantages when compared to traditional methods such as sequence conservation analyses or functional assays. Here we present the adaptation of the ATAC-seq method to echinoderm species and discuss how it can be used for CRE discovery.

Modeling Glaucoma: Retinal Ganglion Cells Generated from Induced Pluripotent Stem Cells of Patients with SIX6 Risk Allele Show Developmental Abnormalities

Glaucoma represents a group of multifactorial diseases with a unifying pathology of progressive retinal ganglion cell (RGC) degeneration, causing irreversible vision loss. To test the hypothesis that RGCs are intrinsically vulnerable in glaucoma, we have developed an in vitro model using the SIX6 risk allele carrying glaucoma patient-specific induced pluripotent stem cells (iPSCs) for generating functional RGCs. Here, we demonstrate that the efficiency of RGC generation by SIX6 risk allele iPSCs is significantly lower than iPSCs-derived from healthy, age- and sex-matched controls. The decrease in the number of RGC generation is accompanied by repressed developmental expression of RGC regulatory genes. The SIX6 risk allele RGCs display short and simple neurites, reduced expression of guidance molecules, and immature electrophysiological signature. In addition, these cells have higher expression of glaucoma-associated genes, CDKN2A and CDKN2B, suggesting an early onset of the disease phenotype. Consistent with the developmental abnormalities, the SIX6 risk allele RGCs display global dysregulation of genes which map on developmentally relevant biological processes for RGC differentiation and signaling pathways such as mammalian target of rapamycin that integrate diverse functions for differentiation, metabolism, and survival. The results suggest that SIX6 influences different stages of RGC differentiation and their survival; therefore, alteration in SIX6 function due to the risk allele may lead to cellular and molecular abnormalities. These abnormalities, if carried into adulthood, may make RGCs vulnerable in glaucoma. Stem Cells 2017;35:2239-2252.

Identification and characterization of variants and a novel 4 bp deletion in the regulatory region of SIX6, a risk factor for primary open-angle glaucoma

Background

Primary open‐angle glaucoma (POAG) is a complex disease of multigenic inheritance and the most common subtype of glaucoma. SIX6 encodes a transcription factor involved in retina, optic nerve, and pituitary development. Previous studies showed a genetic association between the SIX6 locus and POAG, identifying risk alleles. Whether these alleles are present also in the south Indian population is unclear.

Methods

To address this question, the SIX6 gene and an already characterized and highly conserved SIX6 enhancer (Ch14:60974427‐60974430) were sequenced in two south Indian cohorts, respectively, composed of 65/65 and 200/200 POAG cases/age‐matched controls. We next used Taqman‐based allelic discrimination assay to genotype a common variant (rs33912345: c.421A>C) and the rs1048372 SNP in two cohorts, respectively, composed of 557/387 and 590/448 POAG cases/age‐matched controls. An additional cohort of 153 POAG cases was subsequently recruited to assess the association of the rs33912345:c.421A>C and rs10483727 variants with more prominent changes in two POAG diagnostic parameters: retinal nerve fiber layer thickness and vertical cup/disc ratio, using spectral domain optical coherence tomography. The activity of the newly identified enhancer variants was assessed by transgenesis in zebrafish and luciferase assays.

Results

We identified two known rare and two common variants in the SIX6 locus and a novel 4 bp deletion in the analyzed enhancer. Contrary to previous studies, we could not establish a significant association between the rs10483727 and rs33912345:c.421A>C variants and PAOG in the south Indian ethnicity but patients carrying the corresponding C or T risk alleles exhibited a dose‐dependent reduction of the thickness of the retinal nerve fiber layer and a significant increase in the vertical cup/disc ratio. Transgenesis in zebrafish and luciferase assays demonstrated that the newly identified 4 bp deletion significantly reduced reporter expression in cells of the retinal ganglion and amacrine layers, where human SIX6 is expressed.

Conclusion

Altogether, our data further support the implication of SIX6 variants as POAG risk factors and implicates SIX6 haploinsufficiency in POAG pathogenesis.

Diversity, phylogeny and expression patterns of Pou and Six homeodomain transcription factors in hydrozoan jellyfish Craspedacusta sowerbyi

Formation of all metazoan bodies is controlled by a group of selector genes including homeobox genes, highly conserved across the entire animal kingdom. The homeobox genes from Pou and Six classes are key members of the regulation cascades determining development of sensory organs, nervous system, gonads and muscles. Besides using common bilaterian models, more attention has recently been targeted at the identification and characterization of these genes within the basal metazoan phyla. Cnidaria as a diploblastic sister group to bilateria with simple and yet specialized organs are suitable models for studies on the sensory organ origin and the associated role of homeobox genes. In this work, Pou and Six homeobox genes, together with a broad range of other sensory-specific transcription factors, were identified in the transcriptome of hydrozoan jellyfish Craspedacusta sowerbyi. Phylogenetic analyses of Pou and Six proteins revealed cnidarian-specific sequence motifs and contributed to the classification of individual factors. The majority of the Craspedacusta sowerbyi Pou and Six homeobox genes are predominantly expressed in statocysts, manubrium and nerve ring, the tissues with sensory and nervous activities. The described diversity and expression patterns of Pou and Six factors in hydrozoan jellyfish highlight their evolutionarily conserved functions. This study extends the knowledge of the cnidarian genome complexity and shows that the transcriptome of hydrozoan jellyfish is generally rich in homeodomain transcription factors employed in the regulation of sensory and nervous functions.

Comparative analysis of Six 3 and Six 6 distribution in the developing and adult mouse brain

Six 3 and Six 6 genes are two closely related members of the Six/sine oculis family of homeobox containing transcription factors. Their expression and function at early stages of embryonic development has been widely addressed in a variety of species. However, their mRNA distribution during late embryonic, postnatal, and adult brain barely has been analyzed. Here, we show that despite their initial overlap in the anterior neural plate, the expression of Six 3 and Six 6 progressively segregates to different regions during mammalian brain development, maintaining only few areas of partial overlap in the thalamic and hypothalamic regions. Six 3, but not Six 6, is additionally expressed in the olfactory bulb, cerebral cortex, hippocampus, midbrain, and cerebellum. These distinct patterns support the idea that Six 3 and Six 6 are differentially required during forebrain development.

Expression analysis of SIX3 and SIX6 in human tissues reveals differences in expression and a novel correlation between the expression of SIX3 and the genes encoding isocitrate dehyhrogenase and cadherin 18

SIX3 and SIX6 are transcription factors expressed during early stages of eye development. Limited expression data for SIX3 and SIX6 are available in the literature but, to date, there are no reports of the relative levels of expression of these genes throughout the human body and in adult tissues in particular. In this paper, we report extensive real-time quantitative PCR analyses of SIX3 and SIX6 expression in many different tissues of the adult human body, including ocular tissues, and a comparison of expression data with that of many other genes to identify similarity in expression. Using this powerful technique, we have detected a novel statistical correlation between the spatial distribution and the quantitative expression of SIX3 and 5 other transcripts including IDH1, the gene encoding the NADP(+)-dependent enzyme isocitrate dehydrogenase, and cadherin 18, type 2 (CDH14). Our data demonstrate that this novel technique can be used to generate hypotheses by comparison of gene expression profiles to identify possible interactions between genes or gene products.

Six3 and Six6 activity is modulated by members of the groucho family

Six3 and Six6 are two genes required for the specification and proliferation of the eye field in vertebrate embryos, suggesting that they might be the functional counterparts of the Drosophila gene sine oculis (so). Phylogenetic and functional analysis have however challenged this idea, raising the possibility that the molecular network in which Six3 and Six6 act may be different from that described for SO. To address this, we have performed yeast two-hybrid screens, using either Six3 or Six6 as a bait. In this paper, we report the results of the latter screen that led to the identification of TLE1 (a transcriptional repressor of the groucho family) and AES (a potential dominant negative form of TLE proteins) as cofactors for both SIX6 and SIX3. Biochemical and mutational analysis shows that the Six domain of both SIX3 and SIX6 strongly interact with the QD domain of TLE1 and AES, but that SIX3 also interacts with TLE proteins via the WDR domain. Tle1 and Aes are expressed in the developing eye of medaka fish (Oryzias latipes) embryos, overlapping with the distribution of both Six3 and Six6. Gain-of-function studies in medaka show a clear synergistic activity between SIX3/SIX6 and TLE1, which, on its own, can expand the eye field. Conversely, AES alone decreases the eye size and abrogates the phenotypic consequences of SIX3/6 over-expression. These data indicate that both Tle1 and Aes participate in the molecular network that control eye development and are consistent with the view that both Six3 and Six6 act in combination with either Tle1 and/or Aes.