Phylogenetic footprinting reveals a nuclear protein which binds to silencer sequences in the human gamma and epsilon globin genes

Extensive binding of uncharacterized human transcription factors to genomic dark matter

Most of the human genome is thought to be non-functional, and includes large segments often referred to as "dark matter" DNA. The genome also encodes hundreds of putative and poorly characterized transcription factors (TFs). We determined genomic binding locations of 166 uncharacterized human TFs in living cells. Nearly half of them associated strongly with known regulatory regions such as promoters and enhancers, often at conserved motif matches and co-localizing with each other. Surprisingly, the other half often associated with genomic dark matter, at largely unique sites, via intrinsic sequence recognition. Dozens of these, which we term "Dark TFs", mainly bind within regions of closed chromatin. Dark TF binding sites are enriched for transposable elements, and are rarely under purifying selection. Some Dark TFs are KZNFs, which contain the repressive KRAB domain, but many are not: the Dark TFs also include known or potential pioneer TFs. Compiled literature information supports that the Dark TFs exert diverse functions ranging from early development to tumor suppression. Thus, our results sheds light on a large fraction of previously uncharacterized human TFs and their unappreciated activities within the dark matter genome.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state regulatory potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbor distinctive transcription factor binding motifs that are similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we show that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

A chromosome-scale epigenetic map of the Hydra genome reveals conserved regulators of cell state

The epithelial and interstitial stem cells of the freshwater polyp Hydra are the best-characterized stem cell systems in any cnidarian, providing valuable insight into cell type evolution and the origin of stemness in animals. However, little is known about the transcriptional regulatory mechanisms that determine how these stem cells are maintained and how they give rise to their diverse differentiated progeny. To address such questions, a thorough understanding of transcriptional regulation in Hydra is needed. To this end, we generated extensive new resources for characterizing transcriptional regulation in Hydra, including new genome assemblies for Hydra oligactis and the AEP strain of Hydra vulgaris, an updated whole-animal single-cell RNA-seq atlas, and genome-wide maps of chromatin interactions, chromatin accessibility, sequence conservation, and histone modifications. These data revealed the existence of large kilobase-scale chromatin interaction domains in the Hydra genome that contain transcriptionally coregulated genes. We also uncovered the transcriptomic profiles of two previously molecularly uncharacterized cell types: isorhiza-type nematocytes and somatic gonad ectoderm. Finally, we identified novel candidate regulators of cell type-specific transcription, several of which have likely been conserved at least since the divergence of Hydra and the jellyfish Clytia hemisphaerica more than 400 million years ago.

Analysis of a new begomovirus unveils a composite element conserved in the CP gene promoters of several Geminiviridae genera: Clues to comprehend the complex regulation of late genes

A novel bipartite begomovirus, Blechum interveinal chlorosis virus (BleICV), was characterized at the genome level. Comparative analyses revealed that BleICV coat protein (CP) gene promoter is highly divergent from the equivalent region of other begomoviruses (BGVs), with the single exception of Tomato chino La Paz virus (ToChLPV) with which it shares a 23-bp phylogenetic footprint exhibiting dyad symmetry. Systematic examination of the homologous CP promoter segment of 132 New World BGVs revealed the existence of a quasi-palindromic DNA segment displaying a strongly conserved ACTT-(N7)-AAGT core. The spacer sequence between the palindromic motifs is constant in length, but its sequence is highly variable among viral species, presenting a relaxed consensus (TT)GGKCCCY, which is similar to the Conserved Late Element or CLE (GTGGTCCC), a putative TrAP-responsive element. The homologous CP promoter region of Old World BGVs exhibited a distinct organization, with the putative TATA-box overlapping the left half of the ACTT-N7 composite element. Similar CP promoter sequences, dubbed "TATA-associated composite element" or TACE, were found in viruses belonging to different Geminiviridae genera, hence hinting unsuspected evolutionary relationships among those lineages. To get cues about the TACE function, the regulatory function of the CLE was explored in distinct experimental systems. Transgenic tobacco plants harboring a GUS reporter gene driven by a promoter composed by CLE multimers expressed high beta-glucuronidase activity in absence of viral factors, and that expression was increased by begomovirus infection. On the other hand, the TrAP-responsiveness of a truncated CP promoter of Tomato golden mosaic virus (TGMV) was abolished by site-directed mutation of the only CLE present in it, whereas the artificial addition of one CLE to the -125 truncated promoter strongly enhanced the transactivation level in tobacco protoplasts. These results indicate that the CLE is a TrAP-responsive element, hence providing valuable clues to interpret the recurrent association of the CLE with the TACE. On the basis of the aforesaid direct evidences and the insights afforded by the extensive comparative analysis of BleICV CP promoter, we propose that the TACE might be involved in the TrAP-mediated derepression of CP gene in vascular tissues.

MYOD and HAND transcription factors have conserved recognition sites in mTOR promoter: insights from in silico analysis

mTOR regulates multiple cellular processes that are critical for proper maintenance of cell growth and development. However, mechanisms and factors responsible for transcriptional regulation of mTOR are partially known. To identify different transcription factor binding sites in promoter region of mTOR, we performed in silico phylogenetic foot printing analysis of diverse set of human orthologs. Phylogenetic tree for the orthologs was generated to establish the evolutionary relationships among them. Conserved binding sites among the species were predicted by tool MEME. The predicted conserved sites were further analyzed for binding of transcription factors by MatInspector program. Predicted TFs were then integrated with known physical interactions and coexpression data to decipher the important transcriptional regulators of mTOR signaling. Our study suggests that motifs AGGCGGG (+ 15 to + 21) and GGCGGC (+ 60 to + 65) are highly conserved across the species and are recognition sequence for HAND and MYOD transcription factors, respectively. Also these two transcription factors show direct physical interaction in protein-protein interaction map, indicating their regulatory role on expression of mTOR for control of myogenesis. Our study provides novel clues on differential regulation of mTOR under diverse environmental conditions.

Original Research: Stable expression of miR-34a mediates fetal hemoglobin induction in K562 cells

Sickle cell anemia is a common genetic disorder caused by a point mutation in the sixth codon of the β-globin gene affecting people of African descent worldwide. A wide variety of clinical phenotypes ranging from mild to severe symptoms and complications occur due to hemoglobin S polymerization, red blood cell sickling, and vaso-occlusion. Research efforts are ongoing to develop strategies of fetal hemoglobin (HbF; α2γ2) induction to inhibit sickle hemoglobin polymerization and improve clinical outcomes. Insights have been gained from investigating mutations in the β-globin locus or transcription factors involved in the mechanisms of hemoglobin switching. Recent efforts to expand molecular targets that modulate γ-globin expression involve microRNAs that work through posttranscriptional gene regulation. Therefore, the goal of our study was to identify novel microRNA genes involved in fetal hemoglobin expression. Using in silico analysis, we identified a miR-34a binding site in the γ-globin mRNA which was tested for functional relevance. Stable expression of the shMIMIC miR-34a lentivirus vector increased fetal hemoglobin levels in single cell K562 clones consistent with silencing of a γ-globin gene repressor. Furthermore, miR-34a promoted cell differentiation supported by increased expression of KLF1, glycophorin A, and the erythropoietin receptor. Western blot analysis of known negative regulators of γ-globin including YY1, histone deacetylase 1, and STAT3, which are regulated by miR-34a showed no change in YY1 and histone deacetylase 1 levels; however, total- and phosphorylated-STAT3 levels were decreased in single cell miR-34a K562 clones. These data support a mechanism of fetal hemoglobin activation by miR-34a involving STAT3 gene silencing.

Genomic approaches towards finding cis-regulatory modules in animals

Differential gene expression is the fundamental mechanism underlying animal development and cell differentiation. However, it is a challenge to identify comprehensively and accurately the DNA sequences that are required to regulate gene expression: namely, cis-regulatory modules (CRMs). Three major features, either singly or in combination, are used to predict CRMs: clusters of transcription factor binding site motifs, non-coding DNA that is under evolutionary constraint and biochemical marks associated with CRMs, such as histone modifications and protein occupancy. The validation rates for predictions indicate that identifying diagnostic biochemical marks is the most reliable method, and understanding is enhanced by the analysis of motifs and conservation patterns within those predicted CRMs.

Validation of skeletal muscle cis-regulatory module predictions reveals nucleotide composition bias in functional enhancers

We performed a genome-wide scan for muscle-specific cis-regulatory modules (CRMs) using three computational prediction programs. Based on the predictions, 339 candidate CRMs were tested in cell culture with NIH3T3 fibroblasts and C2C12 myoblasts for capacity to direct selective reporter gene expression to differentiated C2C12 myotubes. A subset of 19 CRMs validated as functional in the assay. The rate of predictive success reveals striking limitations of computational regulatory sequence analysis methods for CRM discovery. Motif-based methods performed no better than predictions based only on sequence conservation. Analysis of the properties of the functional sequences relative to inactive sequences identifies nucleotide sequence composition can be an important characteristic to incorporate in future methods for improved predictive specificity. Muscle-related TFBSs predicted within the functional sequences display greater sequence conservation than non-TFBS flanking regions. Comparison with recent MyoD and histone modification ChIP-Seq data supports the validity of the functional regions.

For efficient identification of genomic sequences responsible for regulating gene expression, a number of computer programs have been developed for automatic annotation of these regulatory regions. We searched for potential regulatory regions responsible for controlling the expression of skeletal muscle-specific genes using these programs, and validated the predictions in a popular cell culture model for muscle. We were able to identify 19 previously uncharacterized regulatory regions for muscle genes. The accuracy of the predictions made by these programs leaves much to be desired, leading us to conclude that other signals in addition to the sequence information will be required to achieve sufficient predictive power for genome annotation. Genomic regions with confirmed regulatory function were compared against non-functional sequences, revealing sequence conservation, composition and chromatin modification properties as important signals in determining regulatory region functionality.

Identification of two conserved cis-acting elements, MYCS and P1BS, involved in the regulation of mycorrhiza-activated phosphate transporters in eudicot species

• In this study, six putative promoter regions of phosphate transporter Pht1;3, Pht1;4 and Pht1;5 genes were isolated from eggplant and tobacco using the inverse polymerase chain reaction (iPCR). The isolated sequences show evolutionary conservation and divergence within/between the two groups of Pht1;3 and Pht1;4/Pht1;5. • Histochemical analyses showed that all six promoter fragments were sufficient to drive β-glucuronidase (GUS) expression specifically in arbuscular mycorrhizal (AM) tobacco roots and were confined to distinct cells containing AM fungal structures (arbuscules or intracellular hyphae). • A series of promoter truncation and mutation analyses combined with phylogenetic footprinting of these promoters revealed that at least two cis-regulatory elements--the mycorrhiza transcription factor binding sequence (MYCS) first identified in this study and P1BS--mediated the transcriptional activation of the AM-mediated inorganic phosphate (Pi) transporter genes. Deletion or partial mutation of either of the two motifs in the promoters could cause a remarkable decrease, or even complete absence, of the promoter activity. • Our results propose that uptake of inorganic phosphate (Pi) by AM fungi is regulated, at least partially, in an MYCS- and P1BS-dependent manner in eudicot species. Our finding offers new insights into the molecular mechanisms underlying the coordination between the AM and the Pi signalling pathways.

Modeling the evolution of regulatory elements by simultaneous detection and alignment with phylogenetic pair HMMs

The computational detection of regulatory elements in DNA is a difficult but important problem impacting our progress in understanding the complex nature of eukaryotic gene regulation. Attempts to utilize cross-species conservation for this task have been hampered both by evolutionary changes of functional sites and poor performance of general-purpose alignment programs when applied to non-coding sequence. We describe a new and flexible framework for modeling binding site evolution in multiple related genomes, based on phylogenetic pair hidden Markov models which explicitly model the gain and loss of binding sites along a phylogeny. We demonstrate the value of this framework for both the alignment of regulatory regions and the inference of precise binding-site locations within those regions. As the underlying formalism is a stochastic, generative model, it can also be used to simulate the evolution of regulatory elements. Our implementation is scalable in terms of numbers of species and sequence lengths and can produce alignments and binding-site predictions with accuracy rivaling or exceeding current systems that specialize in only alignment or only binding-site prediction. We demonstrate the validity and power of various model components on extensive simulations of realistic sequence data and apply a specific model to study Drosophila enhancers in as many as ten related genomes and in the presence of gain and loss of binding sites. Different models and modeling assumptions can be easily specified, thus providing an invaluable tool for the exploration of biological hypotheses that can drive improvements in our understanding of the mechanisms and evolution of gene regulation.

Characterization, expression and function of DORMANCY ASSOCIATED MADS-BOX genes from leafy spurge

DORMANCY ASSOCIATED MADS-BOX (DAM) genes are related to AGAMOUS-LIKE 24 and SHORT VEGETATIVE PHASE genes of arabidopsis and are differentially regulated coordinately with endodormancy induction and release in buds of several perennial plant species. DAM genes were first shown to directly impact endodormancy in peach where a deletion of a series of DAM resulted in loss of endodormancy induction. We have cloned and characterized several MADS box genes from the model perennial weed leafy spurge. Leafy spurge DAM genes are preferentially expressed in shoot tips and buds in response to cold temperatures and day length in a manner that is relative to the level of endodormancy induced by various environmental conditions. Over-expression of one DAM gene in arabidopsis delays flowering. Additionally, we show that at least one DAM gene is differentially regulated by chromatin remodeling. Comparisons of the DAM gene promoters between poplar and leafy spurge have identified several conserved sequences that may be important for their expression patterns in response to dormancy-inducing stimuli.

Isolation and phylogenetic footprinting analysis of the 5'-regulatory region of the floral homeotic gene OrcPI from Orchis italica (Orchidaceae)

The nucleotide sequences of regulatory elements from homologous genes can be strongly divergent. Phylogenetic footprinting, a comparative analysis of noncoding regions, can detect putative transcription factor binding sites (TFBSs) shared among the regulatory regions of 2 or more homologous genes. These conserved motifs have the potential to serve the same regulatory function in distantly related taxa. We isolated the 5'-noncoding region of the OrcPI gene, a MADS-box transcription factor involved in flower development in Orchis italica, using the thermal asymmetric interlaced polymerase chain reaction technique. This region (comprising 1352 bp) induced transient beta-glucuronidase expression in the petal tissue of white Rosa hybrida flowers and represents the 5'-regulatory sequence of the OrcPI gene. Phylogenetic footprinting analysis detected conserved regions within the 5'-regulatory sequence of OrcPI and the homologous regions of Oryza sativa, Lilium regale, and Arabidopsis thaliana. Some of these sequences are known TFBSs described in databases of plant regulatory elements. Nucleotide sequence data reported are available in the DDBJ/EMBL/GenBank databases under the following accession numbers: AF198055 promoter region of the PISTILLATA (PI) gene of A. thaliana; AB094985 cDNA of OrcPI (PI/GLOBOSA [PI/GLO] homologue) of O. italica; AB378089 5'-regulatory region of the OrcPI gene of O. italica; AP008211 putative promoter region of OSMADS2 (PI/GLO homologue) of O. sativa; AP008207 putative promoter region of OSMADS4 (PI/GLO homologue) of O. sativa; and AB158292 putative promoter region of the PI/GLO homologue of L. regale.

Binding sites for ETS family of transcription factors dominate the promoter regions of differentially expressed genes in abdominal aortic aneurysms

BACKGROUND

Previously, we identified 3274 distinct differentially expressed genes in abdominal aortic aneurysm (AAA) tissue compared with nonaneurysmal controls. As transcriptional control is responsible for these expression changes, we sought to find common transcriptional elements in the promoter regions of the differentially expressed genes.

METHODS AND RESULTS

We analyzed the up- and downregulated gene sets with Whole Genome rVISTA to determine the transcription factor (TF) binding sites overrepresented in the 5-kb promoter regions of the 3274 genes. The downregulated gene set yielded 144 TF binding sites that were overrepresented in the subset when compared with the entire genome. In contrast, the upregulated gene set yielded only 13 distinct overrepresented TF binding sites. Interestingly, as classified by TRANSFAC, 8 of the 13 TFs binding to these regions belong to the ETS family. Additionally, nuclear factor kB and its subunits p50 and p65 showed enrichment. Immunohistochemical analyses of 10 TFs from the upregulated set showed 9 to be present in AAA tissue. Based on gene ontology analysis of biological process categories of the upregulated target genes of enriched TFs, 10 TFs had enrichment in immune system process among their target genes.

CONCLUSIONS

Our genome-wide analysis provides further evidence of ETS and nuclear factor kB involvement in AAA. Additionally, our results provide novel insight for future studies aiming to dissect the pathogenesis of AAA and have uncovered potential therapeutic targets for AAA prevention.

Assessing phylogenetic motif models for predicting transcription factor binding sites

MOTIVATION

A variety of algorithms have been developed to predict transcription factor binding sites (TFBSs) within the genome by exploiting the evolutionary information implicit in multiple alignments of the genomes of related species. One such approach uses an extension of the standard position-specific motif model that incorporates phylogenetic information via a phylogenetic tree and a model of evolution. However, these phylogenetic motif models (PMMs) have never been rigorously benchmarked in order to determine whether they lead to better prediction of TFBSs than obtained using simple position weight matrix scanning.

RESULTS

We evaluate three PMM-based prediction algorithms, each of which uses a different treatment of gapped alignments, and we compare their prediction accuracy with that of a non-phylogenetic motif scanning approach. Surprisingly, all of these algorithms appear to be inferior to simple motif scanning, when accuracy is measured using a gold standard of validated yeast TFBSs. However, the PMM scanners perform much better than simple motif scanning when we abandon the gold standard and consider the number of statistically significant sites predicted, using column-shuffled 'random' motifs to measure significance. These results suggest that the common practice of measuring the accuracy of binding site predictors using collections of known sites may be dangerously misleading since such collections may be missing 'weak' sites, which are exactly the type of sites needed to discriminate among predictors. We then extend our previous theoretical model of the statistical power of PMM-based prediction algorithms to allow for loss of binding sites during evolution, and show that it gives a more accurate upper bound on scanner accuracy. Finally, utilizing our theoretical model, we introduce a new method for predicting the number of real binding sites in a genome. The results suggest that the number of true sites for a yeast TF is in general several times greater than the number of known sites listed in the Saccharomyces cerevisiae Database (SCPD). Among the three scanning algorithms that we test, the MONKEY algorithm has the highest accuracy for predicting yeast TFBSs.

Tuning in to the signals: noncoding sequence conservation in vertebrate genomes

Aligning and comparing genomic sequences enables the identification of conserved sequence signatures and can enrich for coding and noncoding functional regions. In vertebrates, the comparison of human and rodent genomes and the comparison of evolutionarily distant genomes, such as human and pufferfish, have identified specific sets of 'ultraconserved' sequence elements associated with the control of early development. However, is this just the tip of a 'conservation iceberg' or do these sequences represent a specific class of regulatory element? Studies on the zebrafish phox2b gene region and the ENCODE project suggest that many regulatory elements are not highly conserved, posing intriguing questions about the relationship between noncoding sequence conservation and function and the evolution of regulatory sequences.

Evolutionary rates and patterns for human transcription factor binding sites derived from repetitive DNA

BACKGROUND

The majority of human non-protein-coding DNA is made up of repetitive sequences, mainly transposable elements (TEs). It is becoming increasingly apparent that many of these repetitive DNA sequence elements encode gene regulatory functions. This fact has important evolutionary implications, since repetitive DNA is the most dynamic part of the genome. We set out to assess the evolutionary rate and pattern of experimentally characterized human transcription factor binding sites (TFBS) that are derived from repetitive versus non-repetitive DNA to test whether repeat-derived TFBS are in fact rapidly evolving. We also evaluated the position-specific patterns of variation among TFBS to look for signs of functional constraint on TFBS derived from repetitive and non-repetitive DNA.

RESULTS

We found numerous experimentally characterized TFBS in the human genome, 7-10% of all mapped sites, which are derived from repetitive DNA sequences including simple sequence repeats (SSRs) and TEs. TE-derived TFBS sequences are far less conserved between species than TFBS derived from SSRs and non-repetitive DNA. Despite their rapid evolution, several lines of evidence indicate that TE-derived TFBS are functionally constrained. First of all, ancient TE families, such as MIR and L2, are enriched for TFBS relative to younger families like Alu and L1. Secondly, functionally important positions in TE-derived TFBS, specifically those residues thought to physically interact with their cognate protein binding factors (TF), are more evolutionarily conserved than adjacent TFBS positions. Finally, TE-derived TFBS show position-specific patterns of sequence variation that are highly distinct from random patterns and similar to the variation seen for non-repeat derived sequences of the same TFBS.

CONCLUSION

The abundance of experimentally characterized human TFBS that are derived from repetitive DNA speaks to the substantial regulatory effects that this class of sequence has on the human genome. The unique evolutionary properties of repeat-derived TFBS are perhaps even more intriguing. TE-derived TFBS in particular, while clearly functionally constrained, evolve extremely rapidly relative to non-repeat derived sites. Such rapidly evolving TFBS are likely to confer species-specific regulatory phenotypes, i.e. divergent expression patterns, on the human evolutionary lineage. This result has practical implications with respect to the widespread use of evolutionary conservation as a surrogate for functionally relevant non-coding DNA. Most TE-derived TFBS would be missed using the kinds of sequence conservation-based screens, such as phylogenetic footprinting, that are used to help characterize non-coding DNA. Thus, the very TFBS that are most likely to yield human-specific characteristics will be neglected by the comparative genomic techniques that are currently de rigeur for the identification of novel regulatory sites.

Phylogenetic footprints in fern spore- and seed-specific gene promoters

Spermatophyte seed-storage proteins have descended from a group of proteins involved in cellular desiccation/hydration processes. Conserved protein structures are found across all plant phyla and in the fungi and Archaea. We investigated whether conservation in the coding region sequence is paralleled by common gene regulatory processes. Seed- and spore-specific gene promoters of three phylogenetically diverse plants were analysed by transient and transgenic expression in Arabidopsis thaliana and tobacco. The transcription factors FUS3 and ABI3, which are central regulators of seed maturation processes, interact with cis-motifs of seed-specific promoters from distantly related plants. The promoter of a fern spore-specific gene encoding a seed-storage globulin-like protein exhibits strong seed-specific activity in both Arabidopsis and tobacco. The existence of phylogenetic footprints indicates good conservation of regulatory pathways controlling gene expression in fern spores and in gymnosperm and angiosperm seeds, reflecting the concerted evolution of coding and regulatory regions.

Comparative analysis of orthologous cellulose synthase promoters from Arabidopsis, Populus and Eucalyptus: evidence of conserved regulatory elements in angiosperms

* The cellulose synthase (CesA) gene family encodes the catalytic subunits of a large protein complex responsible for the deposition of cellulose into plant cell walls. Early in vascular plant evolution, the gene family diverged into distinct members with conserved structures and functions (e.g. primary or secondary cell wall biosynthesis). Although the functions and expression domains of CesA genes have been extensively studied in plants, little is known about transcriptional regulation and promoter evolution in this gene family. * Here, comparative sequence analysis of orthologous CesA promoters from three angiosperm genera, Arabidopsis, Populus and Eucalyptus, was performed to identify putative cis-regulatory sequences. The promoter sequences of groups of Arabidopsis genes that are co-expressed with the primary or secondary cell wall-related CesA genes were also analyzed. * Reporter gene analysis of newly isolated promoter regions of six E. grandis CesA genes in Arabidopsis revealed the conserved functionality of the promoter sequences. Comparative sequence analysis identified 71 conserved sequence motifs, of which 66 were significantly over-represented in either primary or secondary wall-associated promoters. * The presence of conserved cis-regulatory elements in the evolutionary distant CesA promoters of Arabidopsis, Populus and Eucalyptus suggests an ancient transcriptional network regulating cellulose biosynthesis in vascular plants.

Time scale for cyclostome evolution inferred with a phylogenetic diagnosis of hagfish and lamprey cDNA sequences

The Cyclostomata consists of the two orders Myxiniformes (hagfishes) and Petromyzoniformes (lampreys), and its monophyly has been unequivocally supported by recent molecular phylogenetic studies. Under this updated vertebrate phylogeny, we performed in silico evolutionary analyses using currently available cDNA sequences of cyclostomes. We first calculated the GC-content at four-fold degenerate sites (GC(4)), which revealed that an extremely high GC-content is shared by all the lamprey species we surveyed, whereas no striking pattern in GC-content was observed in any of the hagfish species surveyed. We then estimated the timing of diversification in cyclostome evolution using nucleotide and amino acid sequences. We obtained divergence times of 470-390 million years ago (Mya) in the Ordovician-Silurian-Devonian Periods for the interordinal split between Myxiniformes and Petromyzoniformes; 90-60 Mya in the Cretaceous-Tertiary Periods for the split between the two hagfish subfamilies, Myxininae and Eptatretinae; 280-220 Mya in the Permian-Triassic Periods for the split between the two lamprey subfamilies, Geotriinae and Petromyzoninae; and 30-10 Mya in the Tertiary Period for the split between the two lamprey genera, Petromyzon and Lethenteron. This evolutionary configuration indicates that Myxiniformes and Petromyzoniformes diverged shortly after the common ancestor of cyclostomes split from the future gnathostome lineage. Our results also suggest that intra-subfamilial diversification in hagfish and lamprey lineages (especially those distributed in the northern hemisphere) occurred in the Cretaceous or Tertiary Periods.

Identification of regulatory elements by gene family footprinting and in vivo analysis

Gene families of recently duplicated but subsequently diverged genes provide an unique opportunity for comparative analysis of regulatory elements. We have studied the human SPRR gene family of small proline rich proteins involved in barrier function of stratified squamous epithelia. These genes are all expressed in normal human keratinocytes, but respond differently to environmental insults. Comparisons of the functional promoter regions allows the rapid identification of both conserved and of novel regulatory elements that appeared after gene duplication. Competitive electrophoretic mobility shift assays can be used to confirm their presence. Here we show the power of gene family footprinting by the identification of two novel elements in the SPRR3 promoter, not present in SPRR1A and SPRR2A. One of these elements binds a protein similar to GAAP-1, a pro-apoptotic activator of IRF-1 and p53. In vivo analysis shows that this element functions as an inhibitor of SPRR3 transcription. The second novel element functions as an activator of promoter activity and is characterized by its A/T rich sequence. The latter interacting protein indeed binds through contacts in the minor groove, and strikingly, depends on the presence of calcium for DNA interaction.

Enhancer identification through comparative genomics

With the availability of genomic sequence from numerous vertebrates, a paradigm shift has occurred in the identification of distant-acting gene regulatory elements. In contrast to traditional gene-centric studies in which investigators randomly scanned genomic fragments that flank genes of interest in functional assays, the modern approach begins electronically with publicly available comparative sequence datasets that provide investigators with prioritized lists of putative functional sequences based on their evolutionary conservation. However, although a large number of tools and resources are now available, application of comparative genomic approaches remains far from trivial. In particular, it requires users to dynamically consider the species and methods for comparison depending on the specific biological question under investigation. While there is currently no single general rule to this end, it is clear that when applied appropriately, comparative genomic approaches exponentially increase our power in generating biological hypotheses for subsequent experimental testing. It is anticipated that cardiac-related genes and the identification of their distant-acting transcriptional enhancers are particularly poised to benefit from these modern capabilities.

Computational biology: toward deciphering gene regulatory information in mammalian genomes

Computational biology is a rapidly evolving area where methodologies from computer science, mathematics, and statistics are applied to address fundamental problems in biology. The study of gene regulatory information is a central problem in current computational biology. This article reviews recent development of statistical methods related to this field. Starting from microarray gene selection, we examine methods for finding transcription factor binding motifs and cis-regulatory modules in coregulated genes, and methods for utilizing information from cross-species comparisons and ChIP-chip experiments. The ultimate understanding of cis-regulatory logic in mammalian genomes may require the integration of information collected from all these steps.

Exploring conservation of transcription factor binding sites with CONREAL

Prediction of transcription factor binding sites (TFBS) is commonly used to formulate working hypotheses for experimental studies on gene regulation. Computational identification of functional TFBS is complicated because of short length and degeneracy of sequence motifs recognized by transcription factors. Information on conservation of predicted sites in orthologous sequences from different species (phylogenetic footprinting) can be used to distinguish potentially functional elements from background predictions. Results of phylogenetic footprinting may substantially depend on the algorithm used to construct an alignment of orthologous sequences, from which conservation of predicted TFBS is estimated. The CONREAL web server allows prediction and comparison of conserved TFBS based on AVID, BLASTZ, CONREAL, and LAGAN alignments. The web tool is particularly suited for the analysis of individual genes or genomic regions, although the underlying algorithm can also be used in high-throughput promoter analysis.

Phylogenetic analysis of the promoter region of the CD40L gene in primates and other mammals

CD40L is a type II membrane protein comprised of 261 amino acids. CD40L plays a crucial role in the immune system where it is primarily expressed on activated T cells and triggers immunoglobulin class switching. The genetic disease X-linked hypergammaglobulinemia (HIGM1, XHIGM or XHIM) is caused by mutations in the CD40L gene. Individuals with HIGM1 are susceptible to recurrent infections to pathogens and a relationship has been shown to exist with malaria [Sabeti, P., Usen, S., Farhadian, S., Jallow, M., Doherty, T., Newport, M., Pinder, M., Ward, R., Kwiatkowski, D., 2002a. CD40L association with protection from severe malaria. Genes Immun. 3, 286-291]. In this paper, we phylogenetically examine the promoter region of CD40L in primates and other mammals via phylogenetic shadowing. This analysis revealed several regions of the CD40L promoter that were highly constrained and thereby inferred to be functional. These constrained regions confirmed many known regulatory sites. In addition, a novel, highly constrained upstream region was also identified which had an NF-AT recognition motif. These analyses also showed that the different mammal groups do not share an exactly similar set of promoter binding sites and taxon-specific promoters have evolved.

Functional conservation of a root hair cell-specific cis-element in angiosperms with different root hair distribution patterns

Vascular plants develop distinctive root hair distribution patterns in the root epidermis, depending on the taxon. The three patterns, random (Type 1), asymmetrical cell division (Type 2), and positionally cued (Type 3), are controlled by different upstream fate-determining factors that mediate expression of root hair cell-specific genes for hair morphogenesis. Here, we address whether these root hair genes possess a common transcriptional regulatory module (cis-element) determining cell-type specificity despite differences in the final root hair pattern. We identified Arabidopsis thaliana expansinA7 (At EXPA7) orthologous (and paralogous) genes from diverse angiosperm species with different hair distribution patterns. The promoters of these genes contain conserved root hair-specific cis-elements (RHEs) that were functionally verified in the Type-3 Arabidopsis root. The promoter of At EXPA7 (Type-3 pattern) also showed hair cell-specific expression in the Type 2 rice (Oryza sativa) root. Root hair-specific genes other than EXPAs also carry functionally homologous RHEs in their promoters. The RHE core consensus was established by a multiple alignment of functionally characterized RHEs from different species and by high-resolution analysis of At EXPA7 RHE1. Our results suggest that this regulatory module of root hair-specific genes has been conserved across angiosperms despite the divergence of upstream fate-determining machinery.

Antiviral activity of CHO-SS cell-derived human omega interferon and other human interferons against HCV RNA replicons and related viruses

The fully glycosylated human omega interferon produced from CHO-SS cells (glycosylated IFN-omega) has been shown to be well-tolerated in man and to induce a sustained virologic response in patients infected with hepatitis C virus (HCV). We examined the antiviral activity of glycosylated IFN-omega and various human IFNs (IFN-alpha, -beta, -gamma and non-glycosylated bacterial (Escherichia coli) recombinant IFN-omega (non-glycosylated IFN-omega)) against HCV RNA replicons and several viruses related to HCV. Since none of the IFNs displayed cytotoxicity we compared their activities based on the effective concentration of the IFN that inhibited virus growth by 50% (EC50). Glycosylated IFN-omega was found to be the most potent antiviral agent of all the IFNs tested against bovine viral diarrhea virus (BVDV), yellow fever virus and West Nile virus. With HCV RNA replicons, non-glycosylated IFN-omega was comparable in activity to IFN-alpha while glycosylated IFN-omega was markedly more potent, indicating that glycosylation has an important effect on its activity. Drug combination analysis of glycosylated IFN-omega+ribavirin (RBV) in BVDV showed a synergy of antiviral effects similar to IFN-alpha+RBV, as well as a unique antagonism of RBV cytotoxic effects by glycosylated IFN-omega. Transcription factor (TF) profiling indicated that IFN-alpha or glycosylated IFN-omega treatment upregulated the same 17 TFs. IFN-alpha and glycosylated IFN-omega also upregulated 9 and 40 additional unique TFs, respectively. The differences in the expression of these TFs were modest, but statistically significantly different for eight of the TFs that were upregulated exclusively by glycosylated IFN-omega. The activation of these additional TFs by glycosylated IFN-omega might contribute to its high potency.

Transforming omics data into context: Bioinformatics on genomics and proteomics raw data

Differential gene expression analysis and proteomics have exerted significant impact on the elucidation of concerted cellular processes, as simultaneous measurement of hundreds to thousands of individual objects on the level of RNA and protein ensembles became technically feasible. The availability of such data sets has promised a profound understanding of phenomena on an aggregate level, expressed as the phenotypic response (observables) of cells, e.g., in the presence of drugs, or characterization of cells and tissue displaying distinct patho-physiological states. However, the step of transforming these data into context, i.e., linking distinct expression or abundance patterns with phenotypic observables - and furthermore enabling a sound biological interpretation on the level of reaction networks and concerted pathways, is still a major shortcoming. This finding is certainly based on the enormous complexity embedded in cellular reaction networks, but a variety of computational approaches have been developed over the last few years to overcome these issues. This review provides an overview on computational procedures for analysis of genomic and proteomic data introducing a sequential analysis workflow: Explorative statistics for deriving a first, from the purely statistical viewpoint, relevant candidate gene/protein list, followed by co-regulation and network analysis to biologically expand this core list toward functional networks and pathways. The review on these procedures is complemented by example applications tailored at identification of disease-associated proteins. Optimization of computational procedures involved, in conjunction with the continuous increase in additional biological data, clearly has the potential of boosting our understanding of processes on a cell-wide level.

Developmental changes in DNA methylation and covalent histone modifications of chromatin associated with the ε-, γ-, and β-globin gene promoters in Papio anubis

The baboon is a suitable and relevant animal model to study the mechanism of human globin gene switching. This investigation addresses the role of DNA methylation and histone coding in globin gene switching in the baboon, Papio anubis. Bisulfite sequencing and chromatin immunoprecipitation studies were performed in erythroid cells purified from fetuses of varying gestational ages and from adult bone marrow to analyze the manner that changes in DNA methylation of the epsilon-, gamma-, and beta-globin promoters and association of ac-H3, ac-H4, H3-dimeK4, H3-dimeK36, and H3-dimeK79 with the epsilon-, gamma-, and beta-globin promoters occur during development. Changes in DNA methylation of the epsilon- and gamma-globin gene promoters during transitional stages of globin gene switching were consistent with the stochastic model of methylation and a role of DNA methylation in gene silencing. Enrichment of ac-H3, ac-H4, and pol II at the promoters of developmentally active genes was observed, while the pattern of distribution of H3-dimeK4 and H3-dimeK79 suggests that these modifications are found near both currently and formerly active promoters. Enrichment of H3-dimeK36 at the silenced epsilon-globin gene promoter was observed. These studies demonstrate that coordinated epigenetic modifications in the chromatin structure of the beta-like globin gene promoters accompany the highly regulated changes in expression patterns of these genes during development.

Rate of promoter class turn-over in yeast evolution

BACKGROUND

Phylogenetic conservation at the DNA level is routinely used as evidence of molecular function, under the assumption that locations and sequences of functional DNA segments remain invariant in evolution. In particular, short DNA segments participating in initiation and regulation of transcription are often conserved between related species. However, transcription of a gene can evolve, and this evolution may involve changes of even such conservative DNA segments. Genes of yeast Saccharomyces have promoters of two classes, class 1 (TATA-containing) and class 2 (non-TATA-containing).

RESULTS

Comparison of upstream non-coding regions of orthologous genes from the five species of Saccharomyces sensu stricto group shows that among 212 genes which very likely have class 1 promoters in S. cerevisiae, 17 probably have class 2 promoters in one or more other species. Conversely, among 322 genes which very likely have class 2 promoters in S. cerevisiae, 44 probably have class 1 promoters in one or more other species. Also, for at least 2 genes from the set of 212 S. cerevisiae genes with class 1 promoters, the locations of the TATA consensus sequences are substantially different between the species.

CONCLUSION

Our results indicate that, in the course of yeast evolution, a promoter switches its class with the probability at least approximately 0.1 per time required for the accumulation of one nucleotide substitution at a non-coding site. Thus, key sequences involved in initiation of transcription evolve with substantial rates in yeast.

Genes involved in complex adaptive processes tend to have highly conserved upstream regions in mammalian genomes

BACKGROUND

Recent advances in genome sequencing suggest a remarkable conservation in gene content of mammalian organisms. The similarity in gene repertoire present in different organisms has increased interest in studying regulatory mechanisms of gene expression aimed at elucidating the differences in phenotypes. In particular, a proximal promoter region contains a large number of regulatory elements that control the expression of its downstream gene. Although many studies have focused on identification of these elements, a broader picture on the complexity of transcriptional regulation of different biological processes has not been addressed in mammals. The regulatory complexity may strongly correlate with gene function, as different evolutionary forces must act on the regulatory systems under different biological conditions. We investigate this hypothesis by comparing the conservation of promoters upstream of genes classified in different functional categories.

RESULTS

By conducting a rank correlation analysis between functional annotation and upstream sequence alignment scores obtained by human-mouse and human-dog comparison, we found a significantly greater conservation of the upstream sequence of genes involved in development, cell communication, neural functions and signaling processes than those involved in more basic processes shared with unicellular organisms such as metabolism and ribosomal function. This observation persists after controlling for G+C content. Considering conservation as a functional signature, we hypothesize a higher density of cis-regulatory elements upstream of genes participating in complex and adaptive processes.

CONCLUSION

We identified a class of functions that are associated with either high or low promoter conservation in mammals. We detected a significant tendency that points to complex and adaptive processes were associated with higher promoter conservation, despite the fact that they have emerged relatively recently during evolution. We described and contrasted several hypotheses that provide a deeper insight into how transcriptional complexity might have been emerged during evolution.

CREDO: a web-based tool for computational detection of conserved sequence motifs in noncoding sequences

SUMMARY

CREDO is a user-friendly, web-based tool that integrates the analysis and results of different algorithms widely used for the computational detection of conserved sequence motifs in noncoding sequences. It enables easy comparison of the individual results. CREDO offers intuitive interfaces for easy and rapid configuration of the applied algorithms and convenient views on the results in graphical and tabular formats.

AVAILABILITY

http://mips.gsf.de/proj/regulomips/credo.htm.

Control of foot differentiation in Hydra: Phylogenetic footprinting indicates interaction of head, bud and foot patterning systems

Homeodomain transcription factor CnNK-2 seems to play a major role in foot formation in Hydra. Recently, we reported in vitro evidence indicating that CnNK-2 has autoregulatory features and regulates expression of the morphogenetic peptide pedibin. We proposed that CnNK-2 and pedibin synergistically orchestrate foot differentiation processes. Here, we further analyzed the regulatory network controlling foot formation in Hydra. By phylogenetic footprinting we compared the CnNK-2 5'-flanking sequence from two closely related species, Hydra vulgaris and Hydra oligactis. Unexpectedly, we detected a highly conserved binding site for HNF-3beta, a vertebrate Forkhead transcription factor, in the CnNK-2 5'-flanking region. The Hydra HNF-3beta homolog budhead is predominantly expressed in the apical region of the body column and early during budding. Budhead is absent from tissue expressing CnNK-2 and thought to be involved in determining tissue for head differentiation. By electrophoretic mobility shift assays we demonstrate an in vitro interaction between recombinant budhead protein and the interspecific conserved HNF-3beta binding motif in the CnNK-2 5'-flanking region. Our results strengthen the view of CnNK-2 as an important regulator during foot patterning processes. Furtheron, they point to budhead as a candidate for a transcriptional regulator of CnNK-2 and to an interaction of foot and head patterning processes in Hydra on the molecular level.

Moving primate genomics beyond the chimpanzee genome

The comparative DNA sequence data that already exist on individual genomic loci depict the phylogenetic relationships of nearly all extant primate genera. Such a phylogenetic representation of the primates, validated by many sequenced primate genomes, and encompassing the full adaptive diversity of the order, is a prerequisite for identifying the genetic basis of humankind, and for testing the proposed human uniqueness of these traits. Some of these traits have been discovered recently, particularly in genes encoding proteins that are important for brain function.

Conservation genomics: applying whole genome studies to species conservation efforts

Studies of complete genomes are leading to a new understanding of the biology of mammals and providing ongoing insights into the fundamental aspects of the organization and evolution of biological systems. Comparison of primate genomes can identify aspects of their organization, regulation and function that appeared during the primate radiation, but without comparison to more evolutionarily distant mammals and other vertebrates, highly conserved aspects of genome architecture will not be accurately identified nor will the lineage-specific changes be identified as such. Many species of primates face risks of extinction; yet the knowledge of their genomes will provide a deeper understanding of primate adaptations, human origins, and provide the framework for discoveries anticipated to improve human medicine. The great apes, the closest relatives of the human species, are among the most vulnerable and most important for human medical studies. However, apes are not the only species whose genomic information will enrich humankind. Comparative genomic studies of endangered species can benefit conservation efforts on their behalf. Increased knowledge of genome makeup and variation in endangered species finds conservation application in population evaluation monitoring and management, understanding phylozoogeography, can enhance wildlife health management, identify risk factors for genetic disorders, and provide insights into demographic management of small populations in the wild and in captivity.

Systematic discovery of regulatory motifs in human promoters and 3' UTRs by comparison of several mammals

Comprehensive identification of all functional elements encoded in the human genome is a fundamental need in biomedical research. Here, we present a comparative analysis of the human, mouse, rat and dog genomes to create a systematic catalogue of common regulatory motifs in promoters and 3' untranslated regions (3' UTRs). The promoter analysis yields 174 candidate motifs, including most previously known transcription-factor binding sites and 105 new motifs. The 3'-UTR analysis yields 106 motifs likely to be involved in post-transcriptional regulation. Nearly one-half are associated with microRNAs (miRNAs), leading to the discovery of many new miRNA genes and their likely target genes. Our results suggest that previous estimates of the number of human miRNA genes were low, and that miRNAs regulate at least 20% of human genes. The overall results provide a systematic view of gene regulation in the human, which will be refined as additional mammalian genomes become available.

Computational prediction of conserved operons and phylogenetic footprinting of transcription regulatory elements in the metal-reducing bacterial family Geobacteraceae

Members of the family Geobacteraceae are an important group of microorganisms from the delta subdivision of Proteobacteria that couple the oxidation of organic compounds to metal reduction. In order to uncover transcription regulatory interactions in these organisms, we used computational methods to identify conserved operons and putative cis-regulatory transcription elements. We identified 26 putative operons with gene order and function conserved among two species of Geobacteraceae, Geobacter sulfurreducens and Geobacter metallireducens. Most of these operons were also conserved in Desulfovibrio vulgaris, an additional metal reducing organism from family Desulfovibrionaceae of the delta subdivision of Proteobacteria. The predicted conserved operons were investigated for the presence of transcription factor binding sites by two different methods, (i) comparison of non-coding regions in conserved operons, and (ii) neural network promoter prediction. Predicted motifs were screened to identify most likely transcription factor binding sites and ribosome-binding sites. We provide information on motifs in Geobacteraceae similar to known transcription factor binding sites in Escherichia coli, conserved motifs in other bacterial species, putative palindromic sites, and predicted ribosome-binding sites. These predictions will aid in further elucidation of regulatory networks of gene interactions in Geobacteraceae.

Comparative genomics

The genomes from three mammals (human, mouse, and rat), two worms, and several yeasts have been sequenced, and more genomes will be completed in the near future for comparison with those of the major model organisms. Scientists have used various methods to align and compare the sequenced genomes to address critical issues in genome function and evolution. This review covers some of the major new insights about gene content, gene regulation, and the fraction of mammalian genomes that are under purifying selection and presumed functional. We review the evolutionary processes that shape genomes, with particular attention to variation in rates within genomes and along different lineages. Internet resources for accessing and analyzing the treasure trove of sequence alignments and annotations are reviewed, and we discuss critical problems to address in new bioinformatic developments in comparative genomics.

PromoLign: a database for upstream region analysis and SNPs

The study of transcriptional regulation at the genomic level has been hindered by the lack of functional annotation in the putative regulatory regions. Phylogenetic footprinting, in which cross-species sequence alignment among orthologous genes is applied to locate conserved sequence blocks, is an effective strategy to attack this problem. Single nucleotide polymorphisms (SNPs) in transcription factor (TF) binding sites contribute to the heterogeneity of TF binding sites and might disrupt or enhance their regulatory activity. The correlation of SNPs with the TF sites will not only help in functional evaluation of SNPs, but will also help in the study of transcription regulation by focusing attention on specific TF sites. PromoLign (http://polly.wustl.edu/promolign/main.html) is an online database application that presents SNPs and TF binding profiles in the context of human-mouse orthologous sequence alignment with a hyperlinked graphical interface. PromoLign could be applied to a variety of SNPs and transcription related studies, including association genetics, population genetics, and pharmacogenetics.

Inhibition of β protein 1 expression enhances β-globin promoter activity and β-globin mRNA levels in the human erythroleukemia (K562) cell line

OBJECTIVE

In this paper, we report new observations related to the mechanism of the negative regulation of the important adult beta-globin gene in the erythroid cells at the embryonic-fetal stage of their development. We focused on the role of the silencer II region located upstream of the beta-globin gene, which along with its cognate binding protein BP1, negatively regulates beta-globin transcription.

MATERIALS AND METHODS

We prepared plasmid constructs containing the wild-type silencer II sequence, a mutated silencer II sequence, or a mutated control sequence in the beta-globin promoter 690-bp insert, which in turn was linked to an enhanced green fluorescent protein (EGFP) reporter gene. A human erythroleukemia cell line (K562) with embryonic-fetal phenotype was transfected with these EGFP constructs.

RESULTS

Flow cytometry and fluorescence digital imaging showed about threefold increase in the beta-globin promoter activity of the mutated silencer II construct. Introduction of a small interfering RNA (siRNA) complementary to BP1 into the cells caused a 75% decrease in BP1 expression and a simultaneous approximately 40% elevation of beta-globin promoter activity as well as an increase in beta-globin mRNA levels, as compared with controls. We detected no changes in the mRNA levels of positive regulators of hemoglobin transcription such as EKLF and GATA-1.

CONCLUSION

Our results support the involvement of BP1 in the mechanism of the negative regulation of beta-globin transcription. A better understanding of this mechanism may lay the groundwork for novel gene therapy approaches to inhibit the expression of abnormal structural variants of adult beta globin, such as sickle hemoglobin.

Genome sequence of the Brown Norway rat yields insights into mammalian evolution

The laboratory rat (Rattus norvegicus) is an indispensable tool in experimental medicine and drug development, having made inestimable contributions to human health. We report here the genome sequence of the Brown Norway (BN) rat strain. The sequence represents a high-quality 'draft' covering over 90% of the genome. The BN rat sequence is the third complete mammalian genome to be deciphered, and three-way comparisons with the human and mouse genomes resolve details of mammalian evolution. This first comprehensive analysis includes genes and proteins and their relation to human disease, repeated sequences, comparative genome-wide studies of mammalian orthologous chromosomal regions and rearrangement breakpoints, reconstruction of ancestral karyotypes and the events leading to existing species, rates of variation, and lineage-specific and lineage-independent evolutionary events such as expansion of gene families, orthology relations and protein evolution.

Comparative genomics

Comparing the genomes of two different species allow the exploration of a host of intriguing evolutionary and genetic questions

Genomic organization, promoter, amino acid sequence, chromosomal localization, and expression of the human gene for CORS-26 (collagenous repeat-containing sequence of 26-kDa protein)

The murine gene for CORS-26 is located on mouse chromosome 15A2 and its expression has been reported to be restricted to fibroblasts, cartilage and kidney. Here, the complete genomic organization of the corresponding human CORS-26 gene with exon/intron boundaries and exon-specific primer combinations is presented. Additionally, a 1.2 kb fragment of the TATA box-containing promoter region was cloned and analyzed for putative transcription factor binding sites. The deduced amino acid sequence is presented completely. Northern blot analysis using a human multiple-tissue cDNA panel demonstrated expression of human CORS-26 mRNA in colon and small intestine. Additionally, RT-PCR analysis revealed expression of CORS-26 mRNA in placenta, fibroblasts and white adipose tissue. The chromosomal localization of the human CORS-26 gene was mapped to human chromosome 5p by fluorescence in situ hybridization (FISH). In humans, chromosomal imbalances on chromosome 5p were reported to be involved in the pathogenesis of osteosarcoma. Therefore, a human bone tumor cDNA panel was investigated and a strong CORS-26 mRNA expression was found in osteosarcoma, chondroblastoma and giant cell tumor. The present data provide the basis for further investigation of CORS-26 gene regulation in the context of mesenchymal tissue development and in the pathogenesis of bone or skeletal disease.

A third isoform of cytochrome c oxidase subunit VIII is present in mammals

The terminal enzyme of the mitochondrial respiratory chain, cytochrome c oxidase (COX), contains three mitochondrial and ten nuclear encoded subunits in mammals. Three of the nuclear subunits (VIa, VIIa, and VIII) have muscle and non-muscle-specific isoforms, subunit IV contains a lung-specific isoform, and subunit VIb contains a testes-specific isoform. For subunit VIII, the smallest nuclear encoded COX polypeptide, we have now found a third gene (COX 8-3), which has been identified in human, lemur, rat, and mouse, suggesting that it is present in a broad range of Eutherian mammals. Sequence similarity and gene structure support the homology of COX8-3 to the other subunit VIII isoforms, indicating that all three are the product of gene duplications. COX VIII-3 protein is mitochondrially-targeted, as shown by a fluorescent COX VIII3/DsRed fusion protein. Both the mitochondrial targeting and its sequence conservation suggest that COXVIII-3 functions as part of the COX holoenzyme and could have a tissue-specific role, as is the case for the other two isoforms. Questions remain about where COX8-3 is predominantly expressed. However, detection of full-length cDNAs, lower levels of sequence divergence at the first and second codon positions compared to the third, and a conserved gene structure indicate that COX VIII-3 is an expressed gene whose origin dates to at least 91 million years ago.

Regulatory elements of the floral homeotic gene AGAMOUS identified by phylogenetic footprinting and shadowing

In Arabidopsis thaliana, cis-regulatory sequences of the floral homeotic gene AGAMOUS (AG) are located in the second intron. This 3-kb intron contains binding sites for two direct activators of AG, LEAFY (LFY) and WUSCHEL (WUS), along with other putative regulatory elements. We have used phylogenetic footprinting and the related technique of phylogenetic shadowing to identify putative cis-regulatory elements in this intron. Among 29 Brassicaceae species, several other motifs, but not the LFY and WUS binding sites identified previously, are largely invariant. Using reporter gene analyses, we tested six of these motifs and found that they are all functionally important for the activity of AG regulatory sequences in A. thaliana. Although there is little obvious sequence similarity outside the Brassicaceae, the intron from cucumber AG has at least partial activity in A. thaliana. Our studies underscore the value of the comparative approach as a tool that complements gene-by-gene promoter dissection but also demonstrate that sequence-based studies alone are insufficient for a complete identification of cis-regulatory sites.

The role of a retinoic acid response element in establishing the anterior neural expression border of Hoxd4 transgenes

The zebrafish hoxd4a locus was compared to its murine ortholog, Hoxd4. The sequence of regulatory elements, including a DR5 type retinoic acid response element (RARE) required for Hoxd4 neural enhancer activity, are highly conserved. Additionally, zebrafish and mouse neural enhancers function identically in transgenic mouse embryos. We tested whether sequence conservation reflects functional importance by altering the spacing and sequence of the RARE in the Hoxd4 neural enhancer. Stabilizing receptor-DNA interactions did not anteriorize transgene expression. By contrast, conversion of the RARE from a DR5 to a DR2 type element decreased receptor-DNA stability and posteriorized expression. Hence, the setting of the Hox anterior expression border is not a simple function of the affinity of retinoid receptors for their cognate element.

Comparative and functional analyses of LYL1 loci establish marsupial sequences as a model for phylogenetic footprinting

Comparative genomic sequence analysis is a powerful technique for identifying regulatory regions in genomic DNA. However, its utility largely depends on the evolutionary distances between the species involved. Here we describe the screening of a genomic BAC library from the stripe-faced dunnart, Sminthopsis macroura, formerly known as the narrow-footed marsupial mouse. We isolated a clone containing the LYL1 locus, completely sequenced the 60.6-kb insert, and compared it with orthologous human and mouse sequences. Noncoding homology was substantially reduced in the human/dunnart analysis compared with human/mouse, yet we could readily identify all promoters and exons. Human/mouse/dunnart alignments of the LYL1 candidate promoter allowed us to identify putative transcription factor binding sites, revealing a pattern highly reminiscent of critical regulatory regions of the LYL1 paralogue, SCL. This newly identified LYL1 promoter showed strong activity in myeloid progenitor cells and was bound in vivo by Fli1, Elf1, and Gata2-transcription factors all previously shown to bind to the SCL stem cell enhancer. This study represents the first large-scale comparative analysis involving marsupial genomic sequence and demonstrates that such comparisons provide a powerful approach to characterizing mammalian regulatory elements.

Globin gene switching in primates

Evolutionary approaches to the identification of DNA sequences required for transcription of the genes of the beta-globin cluster are reviewed. Sequence alignments of non-coding regions from widely divergent species revealed many conserved motifs (phylogenetic footprints) that were putative transcription factor binding sites and candidate binding proteins were identified. The differential timing of the prosimian and simian gamma-globin genes was analyzed by identifying base changes in the vicinity of the phylogenetic footprints. These differential phylogenetic footprints were shown to bind different nuclear factors, and the behavior of constructs with human or galago gamma-promoters in transgenic mice indicated that DNA motifs near the gamma-globin genes are sufficient to determine the developmental stage of expression. Locus control region alignments have identified many conserved sequence differences outside of the hypersensitive sites. Globin protein and mRNA expression profiles during embryological development in a series of catarrhine (Old World monkeys and apes) and platyrrhine (New World monkeys) primates have been determined. While all catarrhines examined to date have globin expression patterns that are highly similar to the well-established human switching behavior, platyrrhines have inactivated their gamma 1 genes by a variety of mechanisms, and have an earlier gamma-beta switch.