Conserved non-coding sequences and transcriptional regulation

Single-Base Resolution Map of Evolutionary Constraints and Annotation of Conserved Elements across Major Grass Genomes

Conserved noncoding sequences (CNSs) are evolutionarily conserved DNA sequences that do not encode proteins but may have potential regulatory roles in gene expression. CNS in crop genomes could be linked to many important agronomic traits and ecological adaptations. Compared with the relatively mature exon annotation protocols, efficient methods are lacking to predict the location of noncoding sequences in the plant genomes. We implemented a computational pipeline that is tailored to the comparisons of plant genomes, yielding a large number of conserved sequences using rice genome as the reference. In this study, we used 17 published grass genomes, along with five monocot genomes as well as the basal angiosperm genome of Amborella trichopoda. Genome alignments among these genomes suggest that at least 12.05% of the rice genome appears to be evolving under constraints in the Poaceae lineage, with close to half of the evolutionarily constrained sequences located outside protein-coding regions. We found evidence for purifying selection acting on the conserved sequences by analyzing segregating SNPs within the rice population. Furthermore, we found that known functional motifs were significantly enriched within CNS, with many motifs associated with the preferred binding of ubiquitous transcription factors. The conserved elements that we have curated are accessible through our public database and the JBrowse server. In-depth functional annotations and evolutionary dynamics of the identified conserved sequences provide a solid foundation for studying gene regulation, genome evolution, as well as to inform gene isolation for cereal biologists.

Differences in enhancer activity in mouse and zebrafish reporter assays are often associated with changes in gene expression

Background

Phenotypic evolution in animals is thought to be driven in large part by differences in gene expression patterns, which can result from sequence changes in cis-regulatory elements (cis-changes) or from changes in the expression pattern or function of transcription factors (trans-changes). While isolated examples of trans-changes have been identified, the scale of their overall contribution to regulatory and phenotypic evolution remains unclear.

Results

Here, we attempt to examine the prevalence of trans-effects and their potential impact on gene expression patterns in vertebrate evolution by comparing the function of identical human tissue-specific enhancer sequences in two highly divergent vertebrate model systems, mouse and zebrafish. Among 47 human conserved non-coding elements (CNEs) tested in transgenic mouse embryos and in stable zebrafish lines, at least one species-specific expression domain was observed in the majority (83%) of cases, and 36% presented dramatically different expression patterns between the two species. Although some of these discrepancies may be due to the use of different transgenesis systems in mouse and zebrafish, in some instances we found an association between differences in enhancer activity and changes in the endogenous gene expression patterns between mouse and zebrafish, suggesting a potential role for trans-changes in the evolution of gene expression.

Conclusions

In total, our results: (i) serve as a cautionary tale for studies investigating the role of human enhancers in different model organisms, and (ii) suggest that changes in the trans environment may play a significant role in the evolution of gene expression in vertebrates.

Foxa2 regulates the expression of Nato3 in the floor plate by a novel evolutionarily conserved promoter

The development of the neural tube into a complex central nervous system involves morphological, cellular and molecular changes, all of which are tightly regulated. The floor plate (FP) is a critical organizing center located at the ventral-most midline of the neural tube. FP cells regulate dorsoventral patterning, differentiation and axon guidance by secreting morphogens. Here we show that the bHLH transcription factor Nato3 (Ferd3l) is specifically expressed in the spinal FP of chick and mouse embryos. Using in ovo electroporation to understand the regulation of the FP-specific expression of Nato3, we have identified an evolutionarily conserved 204 bp genomic region, which is necessary and sufficient to drive expression to the chick FP. This promoter contains two Foxa2-binding sites, which are highly conserved among distant phyla. The two sites can bind Foxa2 in vitro, and are necessary for the expression in the FP in vivo. Gain and loss of Foxa2 function in vivo further emphasize its role in Nato3 promoter activity. Thus, our data suggest that Nato3 is a direct target of Foxa2, a transcription activator and effector of Sonic hedgehog, the hallmark regulator of FP induction and spinal cord development. The identification of the FP-specific promoter is an important step towards a better understanding of the molecular mechanisms through which Nato3 transcription is regulated and for uncovering its function during nervous system development. Moreover, the promoter provides us with a powerful tool for conditional genetic manipulations in the FP.

An intronic enhancer driven by NF-κB contributes to transcriptional regulation of peptidylarginine deiminase type I gene in human keratinocytes

Peptidylarginine deiminases (PADs) catalyze the conversion of protein-bound arginine to citrulline residues. In human epidermis, where filaggrin is the main deiminated protein, three PADs are detected with specific patterns of expression depending on the keratinocyte (KC) differentiation state. Previous characterizations of the PAD-encoding gene promoters have shown that proximal regulation alone is not sufficient to explain this specificity of expression. In this work, we describe an evolutionarily highly conserved nucleotide segment located in the first intron of the PAD1 gene (PADI1). Luciferase reporter assays showed that it enhances the activity of the PADI1 promoter, in a calcium- and orientation-independent manner. Mutation of a putative NF-κB cis-element markedly reduced its enhancer activity, which also confirmed its potential regulatory function. Chromatin immunoprecipitation assays evidenced the binding of both p65 and p50 NF-κB subunits to the cis-element, and RNA interference inhibition assays confirmed that NF-κB contributes to the PADI1 transcriptional control. Furthermore, the intronic enhancer and promoter of PADI1 potentially interact through chromatin looping, as indicated by chromosome conformation capture assays. Our findings provide evidence that an NF-κB-mediated signaling pathway is involved in PADI1 regulation in human epidermal KCs.

The importance of being cis: evolution of orthologous fish and mammalian enhancer activity

Conserved noncoding elements (CNEs) in vertebrate genomes often act as developmental enhancers, but a critical issue is how well orthologous CNE sequences retain the same activity in their respective species, a characteristic important for generalization of model organism studies. To quantify how well CNE enhancer activity has been preserved, we compared the anatomy-specific activities of 41 zebra fish CNEs in zebra fish embryos with the activities of orthologous human CNEs in mouse embryos. We found that 13/41 (∼30%) of the orthologous CNE pairs exhibit conserved positive activity in zebra fish and mouse. Conserved positive activity is only weakly associated with either sequence conservation or the absence of bases undergoing accelerated evolution. A stronger effect is that disparate activity is associated with transcription factor binding site divergence. To distinguish the contributions of cis- versus trans-regulatory changes, we analyzed 13 CNEs in a three-way experimental comparison: human CNE tested in zebra fish, human CNE tested in mouse, and orthologous zebra fish CNE tested in zebra fish. Both cis- and trans-changes affect a significant fraction of CNEs, although human and zebra fish sequences exhibit disparate activity in zebra fish (indicating cis regulatory changes) twice as often as human sequences show disparate activity when tested in mouse and zebra fish (indicating trans regulatory changes). In all four cases where the zebra fish and human CNE display a similar expression pattern in zebra fish, the human CNE also displays a similar expression pattern in mouse. This suggests that the endogenous enhancer activity of ∼30% of human CNEs can be determined from experiments in zebra fish alone, and to identify these CNEs, both the zebra fish and the human sequences should be tested.

Midline signaling and evolution of the forebrain in chordates: a focus on the lamprey Hedgehog case

Lampreys are agnathans (vertebrates without jaws). They occupy a key phylogenetic position in the emergence of novelties and in the diversification of morphology at the dawn of vertebrates. We have used lampreys to investigate the possibility that embryonic midline signaling systems have been a driving force for the evolution of the forebrain in vertebrates. We have focused on Sonic Hedgehog/Hedgehog (Shh/Hh) signaling. In this article, we first review and summarize our recent work on the comparative analysis of embryonic expression patterns for Shh/Hh, together with Fgf8 (fibroblast growth factor 8) and Wnt (wingless-Int) pathway components, in the embryonic lamprey forebrain. Comparison with nonvertebrate chordates on one hand, and jawed vertebrates on the other hand, shows that these morphogens/growth factors acquired new expression domains in the most rostral part of the neural tube in lampreys compared to nonvertebrate chordates, and in jawed vertebrates compared to lampreys. These data are consistent with the idea that changes in Shh, Fgf8 or Wnt signaling in the course of evolution have been instrumental for the emergence and diversification of the telencephalon, a part of the forebrain that is unique to vertebrates. We have then used comparative genomics on Shh/Hh loci to identify commonalities and differences in noncoding regulatory sequences across species and phyla. Conserved noncoding elements (CNEs) can be detected in lamprey Hh introns, even though they display unique structural features and need adjustments of parameters used for in silico alignments to be detected, because of lamprey-specific properties of the genome. The data also show conservation of a ventral midline enhancer located in Shh/Hh intron 2 of all chordates, the very species which possess a notochord and a floor plate, but not in earlier emerged deuterostomes or protostomes. These findings exemplify how the Shh/Hh locus is one of the best loci to study genome evolution with regards to developmental events.

A single-nucleotide deletion in the POMP 5' UTR causes a transcriptional switch and altered epidermal proteasome distribution in KLICK genodermatosis

KLICK syndrome is a rare autosomal-recessive skin disorder characterized by palmoplantar keratoderma, linear hyperkeratotic papules, and ichthyosiform scaling. In order to establish the genetic cause of this disorder, we collected DNA samples from eight European probands. Using high-density genome-wide SNP analysis, we identified a 1.5 Mb homozygous candidate region on chromosome 13q. Sequence analysis of the ten annotated genes in the candidate region revealed homozygosity for a single-nucleotide deletion at position c.-95 in the proteasome maturation protein (POMP) gene, in all probands. The deletion is included in POMP transcript variants with long 5' untranslated regions (UTRs) and was associated with a marked increase of these transcript variants in keratinocytes from KLICK patients. POMP is a ubiquitously expressed protein and functions as a chaperone for proteasome maturation. Immunohistochemical analysis of skin biopsies from KLICK patients revealed an altered epidermal distribution of POMP, the proteasome subunit proteins alpha 7 and beta 5, and the ER stress marker CHOP. Our results suggest that KLICK syndrome is caused by a single-nucleotide deletion in the 5' UTR of POMP resulting in altered distribution of POMP in epidermis and a perturbed formation of the outermost layers of the skin. These findings imply that the proteasome has a prominent role in the terminal differentiation of human epidermis.

Neurogenin 1 (Neurog1) expression in the ventral neural tube is mediated by a distinct enhancer and preferentially marks ventral interneuron lineages

The bHLH transcription factor Neurog1 (Ngn1, Neurod3, neurogenin 1) is involved in neuronal differentiation and cell-type specification in distinct regions of the developing nervous system. Here, transgenic mouse models were developed that use a Bacterial Artificial Chromosome (BAC) containing 208kb flanking the Neurog1 gene to efficiently drive expression of GFP and Cre in all Neurog1 domains. Two characteristics of Neurog1 gene regulation were uncovered. First, a 4kb region previously shown to be sufficient for driving expression of a reporter gene to a subset of the Neurog1 pattern in the developing midbrain, hindbrain, and spinal cord is required uniformly for high levels of expression in all Neurog1 domains, even those not originally identified as being regulated by this region. Second, a 0.8 kb enhancer was identified that is sufficient to drive Neurog1-like expression specifically in the ventral neural tube. Furthermore, Neurog1 progenitor cells in the ventral neural tube are largely fated to interneuron lineages rather than to motoneurons. These studies provide new tools for directing tissue specific expression in the developing neural tube, define Neurog1 lineages in the spinal cord, and further define the complex genomic structure required for obtaining the correct levels and spatial restriction of the neuronal differentiation gene Neurog1.

Cis-regulation and chromosomal rearrangement of the fgf8 locus after the teleost/tetrapod split

The complex expression pattern of fibroblast growth factor 8 (Fgf8) and the cellular responses dependent on concentration of its mRNA in vertebrates suggest that Fgf8 should be tightly controlled at the transcriptional level. We found zebrafish conserved noncoding elements (CNEs) with pan-vertebrate as well as fish-specific orthologous sequences from across 200 kb of the zebrafish fgf8a genomic regulatory block to direct reporter expression in patterns consistent with the expression pattern of fgf8a. These included elements from inside the introns of the skin-specific slc2a15a and the ubiquitously expressed fbxw4 bystander genes. The fgf8a/fbxw4 gene pair, which has remained joined throughout three whole genome duplications in chordate evolution, is inverted in teleost genomes, but CNEs across both evolutionary breakpoints showed specific activity. While some CNEs directed highly reproducible expression patterns, others were subject to variation but showed, in a subset of transgenes, expression in the apical ectodermal ridge, the anterior boundaries of somites and the midbrain-hindbrain boundary, specific Fgf8 signaling domains, suggesting that their activity may be context specific. A human element with tetrapod-specific orthologous sequences directed reporter expression to the vasculature, possibly corresponding to a tetrapod innovation. We conclude that fgf8a transcriptional regulation employs pan-vertebrate and teleost-specific enhancers dispersed over three genes in the zebrafish genome.

Epigenetic regulation of cytokine gene expression in T lymphocytes

The developmental program of T helper and regulatory T cell lineage commitment is governed by both genetic and epigenetic mechanisms. The principal events, signaling pathways and the lineage determining factors involved have been extensively studied in the past ten years. Recent studies have elucidated the important role of chromatin remodeling and epigenetic changes for proper regulation of gene expression of lineage-specific cytokines. These include DNA methylation and histone modifications in epigenomic reprogramming during T helper cell development and effector T cell functions. This review discusses the basic epigenetic mechanisms and the role of transcription factors for the differential cytokine gene regulation in the T helper lymphocyte subsets.