Conserved sequences and the evolution of gene regulatory signals

Evaluation of conserved and ultra-conserved non-genic sequences in chromosome 15q15-linked periodic catatonia

Conserved and ultra-conserved non-genic sequence elements (CNGs, UCEs) between human and other mammalian genomes seem to constitute a heterogeneous group of functional sequences which likely have important biological function. To determine whether variation in CNGs and UCEs contributes to risk for the schizophrenic subphenotype of periodic catatonia (according to K. Leonhard; OMIM 605419), we evaluated non-coding elements at a critical 7.35 Mb interval on chromosome 15q15 in 8 unrelated cases with periodic catatonia (derived from pedigrees compatible with linkage to chromosome 15q15) and 8 controls, followed by association studies in a cohort of 510 cases and controls. Among 65 CNGs (≥100 bp, 100% identity; human-mouse comparison), 7 CNGs matched criteria for UCE (≥200  bp, 100% identity). A hot spot of 62/65 CNGs (95%) appeared at the MEIS2 locus, which implicates functional importance of associated (ultra-)conserved elements to this early developmental gene, which is present in the human fetal neocortex and associated with metabolic side effects to antipsychotic drugs. Further CNGs were identified at the PLCB2 and DLL4 locus or located intergenic between TYRO3 and MAPKBP1. Automated sequencing revealed genetic variation in 12.3% of CNGs, but frequencies were low (MAF: 0.06-0.4) in cases. Three variants located inside CNGs/UCEs were found in cases only. In a case-control association study we could not confirm a significant association of these three CNG-variants with periodic catatonia. Our results suggest genetic variation in (ultra-)conserved non-genic sequence elements which might alter functional properties. The identified variants are genetically not associated with the phenotype of periodic catatonia.

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.

Resources and transgenesis techniques for functional genomics in Xenopus

Recent developments in genomic resources and high-throughput transgenesis techniques have allowed Xenopus to 'metamorphose' from a classic model for embryology to a leading-edge experimental system for functional genomics. This process has incorporated the fast-breeding diploid frog, Xenopus tropicalis, as a new model-system for vertebrate genomics and genetics. Sequencing of the X. tropicalis genome is nearly complete, and its comparison with mammalian sequences offers a reliable guide for the genome-wide prediction of cis-regulatory elements. Unique cDNA sets have been generated for both X. tropicalis and X. laevis, which have facilitated non-redundant, systematic gene expression screening and comprehensive gene expression analysis. A variety of transgenesis techniques are available for both X. laevis and X. tropicalis, and the appropriate procedure may be chosen depending on the purpose for which it is required. Effective use of these resources and techniques will help to reveal the overall picture of the complex wiring of gene regulatory networks that control vertebrate development.

A publish-subscribe model of genetic networks

We present a simple model of genetic regulatory networks in which regulatory connections among genes are mediated by a limited number of signaling molecules. Each gene in our model produces (publishes) a single gene product, which regulates the expression of other genes by binding to regulatory regions that correspond (subscribe) to that product. We explore the consequences of this publish-subscribe model of regulation for the properties of single networks and for the evolution of populations of networks. Degree distributions of randomly constructed networks, particularly multimodal in-degree distributions, which depend on the length of the regulatory sequences and the number of possible gene products, differed from simpler Boolean NK models. In simulated evolution of populations of networks, single mutations in regulatory or coding regions resulted in multiple changes in regulatory connections among genes, or alternatively in neutral change that had no effect on phenotype. This resulted in remarkable evolvability in both number and length of attractors, leading to evolved networks far beyond the expectation of these measures based on random distributions. Surprisingly, this rapid evolution was not accompanied by changes in degree distribution; degree distribution in the evolved networks was not substantially different from that of randomly generated networks. The publish-subscribe model also allows exogenous gene products to create an environment, which may be noisy or stable, in which dynamic behavior occurs. In simulations, networks were able to evolve moderate levels of both mutational and environmental robustness.

Expression levels of FAS are regulated through an evolutionary conserved element in intron 2, which modulates cystic fibrosis disease severity

We have analyzed frequent naturally occurring variants in the autogene FAS in two independent cystic fibrosis (CF) patient populations. Analysis of FAS expression levels from intestinal epithelial biopsies from 16 unrelated F508del-CFTR homozygotes showed a correlation between FAS intron 2 SNP rs7901656 and signals for Affymetrix GeneChip U133 Plus 2.0 probeset 204781_s_at consistent with a dominant model (P=0.0009). Genotype and haplotype analysis at six informative SNPs spanning the FAS gene locus was carried out on 37 nuclear families representing extreme clinical phenotypes that were selected from the European CF Twin and Sibling Study population of more than 300 affected sibling pairs. Case-control comparison of the haplotype composed of rs2296603-rs7901656-rs1571019 encompassing intron 2 of FAS reached significance (P=0.0246). Comparative phylogenetic analysis and functional annotation of the FAS intron 2 sequence revealed a conserved non-coding sequence surrounding rs7901656 and predicted binding sites for four transcription factors whereby the binding site of c-Rel is altered by rs7901656. Taken together, these findings from two independent CF patient cohorts indicate that allelic variants within FAS intron 2 alter FAS gene expression and that these functional variants modulate the manifestation of CF disease.

Computational methods to dissect cis-regulatory transcriptional networks

The formation of diverse cell types from an invariant set of genes is governed by biochemical and molecular processes that regulate gene activity. A complete understanding of the regulatory mechanisms of gene expression is the major function of genomics. Computational genomics is a rapidly emerging area for deciphering the regulation of metazoan genes as well as interpreting the results of high-throughput screening. The integration of computer science with biology has expedited molecular modelling and processing of large-scale data inputs such as microarrays, analysis of genomes, transcriptomes and proteomes. Many bioinformaticians have developed various algorithms for predicting transcriptional regulatory mechanisms from the sequence, gene expression and interaction data. This review contains compiled information of various computational methods adopted to dissect gene expression pathways.

Comparative genomic sequence and expression analyses of Medicago truncatula and alfalfa subspecies falcata COLD-ACCLIMATION-SPECIFIC genes

In Arabidopsis (Arabidopsis thaliana) the low-temperature induction of genes encoding the C-REPEAT BINDING FACTOR (CBF) transcriptional activators is a key step in cold acclimation. CBFs in turn activate a battery of downstream genes known as the CBF regulon, which collectively act to increase tolerance to low temperatures. Fundamental questions are: What determines the size and scope of the CBF regulon, and is this is a major determinant of the low-temperature tolerance capacity of individual plant species? Here we have begun to address these questions through comparative analyses of Medicago truncatula and Medicago sativa subsp. falcata. M. truncatula survived to -4 degrees C but did not cold acclimate, whereas Medicago falcata cold acclimated and survived -14 degrees C. Both species possessed low-temperature-induced CBFs but differed in the expression of the COLD-ACCLIMATION-SPECIFIC (CAS) genes, which are candidate CBF targets. M. falcata CAS30 was robustly cold-responsive whereas the MtCAS31 homolog was not. M. falcata also possessed additional CAS30 homologs in comparison to the single CAS31 gene in M. truncatula. MfCAS30 possessed multiple pairs of closely spaced C-REPEAT/DEHYDRATION RESPONSIVE ELEMENT (CRT/DRE) motifs, the cognate CBF binding site in its upstream region whereas MtCAS31 lacked one CRT/DRE partner of the two proximal partner pairs. CAS genes also shared a promoter structure comprising modules proximal and distal to the coding sequence. CAS15, highly cold-responsive in both species, harbored numerous CRT/DRE motifs, but only in the distal module. However, fusion of the MtCAS15 promoter, including the distal module, to a reporter gene did not result in low-temperature responsiveness in stably transformed Arabidopsis. In contrast, both MtCAS31 and MfCAS30 promoter fusions were low-temperature responsive, although the MfCAS31 fusion was less robust than the MfCAS30 fusion. From these studies we conclude that CAS genes harbor CRT/DRE motifs, their proximity to one another is likely key to regulatory output in Medicago, and they may be located kilobases distal to the transcriptional start site. We hypothesize that these differences in CRT/DRE copy numbers in CAS30/CAS31 upstream regions combined with differences in gene copy numbers may be a factor in determining differences in low-temperature tolerance between M. truncatula and M. falcata.

Genome-wide B1 retrotransposon binds the transcription factors dioxin receptor and Slug and regulates gene expression in vivo

Alterations in tissue-specific gene expression greatly affect cell function. Transcription factors (TFs) interact with cis-acting binding sites in noncoding enhancer promoter regions. Transposable elements (TEs) are abundant and similarly represented among mammalian genomes. TEs are important in gene regulation, but their function is not well understood. We have characterized a TE containing functional TF-binding sites for the carcinogen-activated dioxin receptor xenobiotic responsive element (XRE) and the epithelial-mesenchymal transition regulator Slug (Slug site). A Mus promoter database was scanned for XREs to predict coregulation with other TFs. We identified an overrepresented (1,398 genes) B1 retrotransposon containing XRE and Slug sites within 35 bp of each other (designated as B1-X35S). This B1-X35S retrotransposon differed from classic B1s by the presence of the Slug site and by its differential nucleotide conservation outside the X35S region. Phylogenetically, B1-X35S appeared recently in evolution, close to the B1-B subfamily. Comparative gene expression in 61 mouse tissues revealed that B1-X35S-containing genes had lower median expression levels than those with canonical B1 TEs, suggesting a repressive role for X35S. Indeed, X35S was functional and able to bind aryl hydrocarbon (dioxin) receptor (AhR) and Slug and, importantly, to repress cis-reporter genes. Moreover, AhR and Slug were recruited to X35S in vivo and repressed the endogenous expression of X35S-containing genes. Our results demonstrate the existence of a widely present B1 subfamily in the mouse. Because AhR and Slug are relevant in tumor development and differentiation, X35S may represent a genome-wide regulatory mechanism and a tool to modulate gene expression.

Activation-dependent intrachromosomal interactions formed by the TNF gene promoter and two distal enhancers

Here we provide a mechanism for specific, efficient transcription of the TNF gene and, potentially, other genes residing within multigene loci. We identify and characterize highly conserved noncoding elements flanking the TNF gene, which undergo activation-dependent intrachromosomal interactions. These elements, hypersensitive site (HSS)-9 and HSS+3 (9 kb upstream and 3 kb downstream of the TNF gene, respectively), contain DNase I hypersensitive sites in naive, T helper 1, and T helper 2 primary T cells. Both HSS-9 and HSS+3 inducibly associate with acetylated histones, indicative of chromatin remodeling, bind the transcription factor nuclear factor of activated T cells (NFAT)p in vitro and in vivo, and function as enhancers of NFAT-dependent transactivation mediated by the TNF promoter. Using the chromosome conformation capture assay, we demonstrate that upon T cell activation intrachromosomal looping occurs in the TNF locus. HSS-9 and HSS+3 each associate with the TNF promoter and with each other, circularizing the TNF gene and bringing NFAT-containing nucleoprotein complexes into close proximity. TNF gene regulation thus reveals a mode of intrachromosomal interaction that combines a looped gene topology with interactions between enhancers and a gene promoter.

Lef1 regulates COX-2 transcription in chondrocytes

Although the involvement of Wnt/beta-catenin signaling pathway in the regulation of COX-2 expression has been suggested, the underlying molecular mechanisms are still unclear. Here, we demonstrate that lymphoid enhancer-binding factor 1 (Lef1) in concert with beta-catenin is associated with the transcriptional regulation of the gene encoding COX-2 in mouse chondrocytes. Our results show that co-overexpression of Lef1 and beta-catenin synergistically upregulates COX-2 expression. Furthermore, decrease of Lef1 expression using Lef1 siRNA results in the downregulation of COX-2 expression. Using bioinformatic analysis, we identified a conserved Lef1-binding site that is mapped at the 3' region of the genomic COX-2 locus. In vivo and in vitro binding of Lef1 at the predicted binding site was proved using chromatin immunoprecipitation assays and electrophoretic mobility shift assays, respectively. Moreover, binding of Lef1 and beta-catenin to the Lef1-binding site transactivates COX-2 promoter activity. Our results indicate a pivotal role of Lef1 in the regulation of COX-2 transcription in arthritic chondrocytes.

Identification of eukaryotic promoter regulatory elements using nonhomologous random recombination

Understanding the regulatory logic of a eukaryotic promoter requires the elucidation of the regulatory elements within that promoter. Current experimental or computational methods to discover regulatory motifs within a promoter can be labor intensive and may miss redundant, unprecedented or weakly activating elements. We have developed an unbiased combinatorial approach to rapidly identify new upstream activating sequences (UASs) in a promoter. This approach couples nonhomologous random recombination with an in vivo screen to efficiently identify UASs and does not rely on preconceived hypotheses about promoter regulation or on similarity to known activating sequences. We validated this method using the unfolded protein response (UPR) in yeast and were able to identify both known and potentially novel UASs involved in the UPR. One of the new UASs discovered using this approach implicates Crz1 as a possible activator of Hac1, a transcription factor involved in the UPR. This method has several advantages over existing methods for UAS discovery including its speed, potential generality, sensitivity and lack of false positives and negatives.

Shuffling of cis-regulatory elements is a pervasive feature of the vertebrate lineage

Alignment of orthologous vertebrate loci reveals that a significant proportion of conserved cis-regulatory elements have undergone shuffling during evolution.

Background

All vertebrates share a remarkable degree of similarity in their development as well as in the basic functions of their cells. Despite this, attempts at unearthing genome-wide regulatory elements conserved throughout the vertebrate lineage using BLAST-like approaches have thus far detected noncoding conservation in only a few hundred genes, mostly associated with regulation of transcription and development.

Results

We used a unique combination of tools to obtain regional global-local alignments of orthologous loci. This approach takes into account shuffling of regulatory regions that are likely to occur over evolutionary distances greater than those separating mammalian genomes. This approach revealed one order of magnitude more vertebrate conserved elements than was previously reported in over 2,000 genes, including a high number of genes found in the membrane and extracellular regions. Our analysis revealed that 72% of the elements identified have undergone shuffling. We tested the ability of the elements identified to enhance transcription in zebrafish embryos and compared their activity with a set of control fragments. We found that more than 80% of the elements tested were able to enhance transcription significantly, prevalently in a tissue-restricted manner corresponding to the expression domain of the neighboring gene.

Conclusion

Our work elucidates the importance of shuffling in the detection of cis-regulatory elements. It also elucidates how similarities across the vertebrate lineage, which go well beyond development, can be explained not only within the realm of coding genes but also in that of the sequences that ultimately govern their expression.

Sequence analyses to study the evolutionary history and cis-regulatory elements of Hedgehog genes

Sequence analysis and comparative genomics are powerful tools to gain knowledge on multiple aspects of gene and protein regulation and function. These have been widely used to understand the evolutionary history and the biochemistry of Hedgehog (Hh) proteins, and the molecular control of Hedgehog gene expression. Here, we report on some of the methods available to retrieve protein and genomic sequences. We describe how protein sequence comparison can produce information on the evolutionary history of Hh proteins. Moreover, we describe the use of genomic sequence analysis including phylogenetic footprinting and transcription factor-binding site search tools, techniques that allow for the characterization of cis-regulatory elements of developmental genes such as the Hedgehog genes.