Identification of two evolutionarily conserved and functional regulatory elements in intron 2 of the human BRCA1 gene

Cross-species comparative genomics is a powerful strategy for identifying functional regulatory elements within noncoding DNA. In this paper, comparative analysis of human and mouse intronic sequences in the breast cancer susceptibility gene (BRCA1) revealed two evolutionarily conserved noncoding sequences (CNS) in intron 2, 5 kb downstream of the core BRCA1 promoter. The functionality of these elements was examined using homologous-recombination-based mutagenesis of reporter gene-tagged cosmids incorporating these regions and flanking sequences from the BRCA1 locus. This showed that CNS-1 and CNS-2 have differential transcriptional regulatory activity in epithelial cell lines. Mutation of CNS-1 significantly reduced reporter gene expression to 30% of control levels. Conversely mutation of CNS-2 increased expression to 200% of control levels. Regulation is at the level of transcription and shows promoter specificity. Both elements also specifically bind nuclear proteins in vitro. These studies demonstrate that the combination of comparative genomics and functional analysis is a successful strategy to identify novel regulatory elements and provide the first direct evidence that conserved noncoding sequences in BRCA1 regulate gene expression. (c) 2005 Elsevier Inc. All rights reserved.

Ultraconserved elements in the human genome

There are 481 segments longer than 200 base pairs (bp) that are absolutely conserved (100% identity with no insertions or deletions) between orthologous regions of the human, rat, and mouse genomes. Nearly all of these segments are also conserved in the chicken and dog genomes, with an average of 95 and 99% identity, respectively. Many are also significantly conserved in fish. These ultraconserved elements of the human genome are most often located either overlapping exons in genes involved in RNA processing or in introns or nearby genes involved in the regulation of transcription and development. Along with more than 5000 sequences of over 100 bp that are absolutely conserved among the three sequenced mammals, these represent a class of genetic elements whose functions and evolutionary origins are yet to be determined, but which are more highly conserved between these species than are proteins and appear to be essential for the ontogeny of mammals and other vertebrates.

Molecular analysis of the pathway for the synthesis of thiol tripeptides in the model legume Lotus japonicus

The thiol tripeptides, glutathione (GSH) and homoglutathione (hGSH), perform multiple roles in legumes, including protection against toxicity of free radicals and heavy metals. The three genes involved in the synthesis of GSH and hGSH in the model legume, Lotus japonicus, have been fully characterized and appear to be present as single copies in the genome. The gamma-glutamylcysteine synthetase (gammaecs) gene was mapped on the long arm of chromosome 4 (70.0 centimorgans [cM]) and consists of 15 exons, whereas the glutathione synthetase (gshs) and homoglutathione synthetase (hgshs) genes were mapped on the long arm of chromosome 1 (81.3 cM) and found to be arranged in tandem, with a separation of approximately 8 kb. Both genes consist of 12 exons of exactly the same size (except exon 1, which is similar). Two types of transcripts were detected for the gshs gene, which putatively encode proteins localized in the plastids and cytosol. Promoter regions contain cis-acting regulatory elements that may be involved in the plant's response to light, hormones, and stress. Determination of transcript levels, enzyme activities, and thiol contents in nodules, roots, and leaves revealed that gammaecs and hgshs are expressed in all three plant organs, whereas gshs is significantly functional only in nodules. This strongly suggests an important role of GSH in the rhizobia-legume symbiosis.

Isolation of an imaginal disc growth factor homologue from Pieris rapae and its expression following parasitization by Cotesia rubecula

Endoparasitoid insects introduce maternal factors into the body of their host at oviposition to suppress cellular defences for the protection of the developing parasitoid. We have shown that transient expression of polydnavirus genes from a hymenopteran parasitoid Cotesia rubecula (CrPDV) is responsible for the inactivation of hemocytes from the lepidopteran host Pieris rapae. Since the observed downregulation of CrPDV genes in infected host tissues is not due to cis-regulatory elements at the CrV1 gene locus, we speculated that the termination of CrPDV gene expression may be due to cellular inactivation caused by the CrV1-mediated immune suppression of infected tissues. To test this assumption, we isolated an imaginal disc growth factor (IDGF) that is expressed in fat body and hemocytes, the target of viral infection and expression of CrPDV genes. Time-course experiments showed that the level of P. rapae IDGF is not affected by parasitization and polydnavirus infection. However, the amount of highly expressed genes, such as storage proteins, arylphorin and lipophorin, are significantly reduced following parasitization. (C) 2004 Elsevier Ltd. All rights reserved.

Understanding alternative splicing: Towards a cellular code

In violation of the 'one gene, one polypeptide' rule, alternative splicing allows individual genes to produce multiple protein isoforms - thereby playing a central part in generating complex proteomes. Alternative splicing also has a largely hidden function in quantitative gene control, by targeting RNAs for nonsense-mediated decay. Traditional gene-by-gene investigations of alternative splicing mechanisms are now being complemented by global approaches. These promise to reveal details of the nature and operation of cellular codes that are constituted by combinations of regulatory elements in pre-mRNA substrates and by cellular complements of splicing regulators, which together determine regulated splicing pathways.

Analysis of distinct tartrate-resistant acid phosphatase promoter regions in transgenic mice

The tartrate-resistant acid phosphatase (TRAP) is present in multiple tissues, including kidney, liver, lung, spleen, and bone. Recent study of (TRAP) gene expression has provided evidence for distinct promoters within the (TRAP) gene, suggesting that the gene has alternative, tissue-preferred mRNA transcripts. Examination of endogenous (TRAP) exon 1B and 1C mRNA transcripts revealed tissue-preferred transcript abundance with increased exon 1B transcripts detected in liver and kidney and increased exon 1C transcripts detected in bone and spleen. In this investigation, we have made transgenic mice that express a marker gene driven by two candidate promoters, designated BC and C, within the (TRAP) gene. The BC and C promoters are 2.2 and 1.6 kb, respectively, measured from the translation initiation site. Evaluation of BC transgenic lines demonstrated robust expression in multiple tissues. In contrast, significant transgene expression was not detected in C transgenic lines. Evaluation of transgene mRNAs in BC transgenic lines revealed that virtually all expression was in the form of B transcripts, suggesting that the tissue-preferred pattern of endogenous (TRAP) was not replicated in the BC transgenic line. Likewise, osteoclastogenic cultures from BC, but not C, transgenic bone marrow cells expressed the transgene following receptor activator of NFkappaB ligand/macrophage colony-stimulating factor stimulation. In conclusion, when compared with the 2.2-kb BC portion of the (TRAP) promoter region, the 1.6-kb C portion does not account for significant gene expression in vivo or in vitro; production of the bone- and spleen-preferred (TRAP) C transcript must depend on regulatory elements outside of the 2.2-kb promoter. As the majority of currently investigated transcription factors that influence transcriptional regulation of osteoclast gene expression bind within the 1.6-kb C portion of the (TRAP) promoter, it is likely that transcription binding sites outside of the 2.2-kb region will have profound effects on regulation of the gene in vivo and in vitro.

Assessing computational tools for the discovery of transcription factor binding sites

The prediction of regulatory elements is a problem where computational methods offer great hope. Over the past few years, numerous tools have become available for this task. The purpose of the current assessment is twofold: to provide some guidance to users regarding the accuracy of currently available tools in various settings, and to provide a benchmark of data sets for assessing future tools.

Evolutionary turnover of mammalian transcription start sites

Alignments of homologous genomic sequences are widely used to identify functional genetic elements and study their evolution. Most studies tacitly equate homology of functional elements with sequence homology. This assumption is violated by the phenomenon of turnover, in which functionally equivalent elements reside at locations that are nonorthologous at the sequence level. Turnover has been demonstrated previously for transcription-factor-binding sites. Here, we show that transcription start sites of equivalent genes do not always reside at equivalent locations in the human and mouse genomes. We also identify two types of partial turnover, illustrating evolutionary pathways that could lead to complete turnover. These findings suggest that the signals encoding transcription start sites are highly flexible and evolvable, and have cautionary implications for the use of sequence-level conservation to detect gene regulatory elements.

The prion protein gene: Identifying regulatory signals using marsupial sequence

The function of the prion protein gene (PRNP) and its normal product PrPC is elusive. We used comparative genomics as a strategy to understand the normal function of PRNP. As the reliability of comparisons increases with the number of species and increased evolutionary distance, we isolated and sequenced a 66.5 kb BAC containing the PRNP gene from a distantly related mammal, the model Australian marsupial Macropus eugenii (tammar wallaby). Marsupials are separated from eutherians such as human and mouse by roughly 180 million years of independent evolution. We found that tammar PRNP, like human PRNP, has two exons. Prion proteins encoded by the tammar wallaby and a distantly related marsupial, Monodelphis domestica (Brazilian opossum) PRNP contain proximal PrP repeats with a distinct, marsupial-specific composition and a variable number. Comparisons of tammar wallaby PRNP with PRNPs from human, mouse, bovine and ovine allowed us to identify non-coding gene regions conserved across the marsupial-eutherian evolutionary distance, which are candidates for regulatory regions. In the PRNP 3' UTR we found a conserved signal for nuclear-specific polyadenylation and the putative cytoplasmic polyadenylation element (CPE), indicating that post-transcriptional control of PRNP mRNA activity is important. Phylogenetic footprinting revealed conserved potential binding sites for the MZF-1 transcription factor in both upstream promoter and intron/intron 1, and for the MEF2, MyTI, Oct-1 and NFAT transcription factors in the intron(s). The presence of a conserved NFAT-binding site and CPE indicates involvement of PrPC in signal transduction and synaptic plasticity. (c) 2004 Elsevier B.V. All rights reserved.

Computational analysis of base composition pattern and promoter elements in the putative promoter regions in relation to expression profiles of 682 human genes on chromosome 22

The base composition pattern (BCP) in the putative promoter region (PPRs) up to 5 Kb lengths of 682 human genes on Chromosome 22 (Chr22) was examined. Two-dimensional (2D) and three-dimensional (3D) functions were designed to delineate the DNA base composition, with four major patterns identified. It is found that 17.6% genes include TATA box, 28.0% GC box, 18.9% CAAT box and 38.4% CpG islands, and approximately 10% genes have one of four putative initiator (Inr) motifs. The occurrence of the promoter elements is tightly associated with the base composition features in the promoter regions, and the associations of the base composition features with occurrence of the promoter elements in the promoter regions mediate tissue-wide expression of the genes in human. The occurrence of two or more promoter elements in the promoter regions is required for the medium- and wide-range expression profiles of the human genes on Chr22. Thus, the reported data shed light on the characteristics of the PPRs of the human genes on Chr22, which may improve our understanding of regulatory roles of the PPRs with occurrence of the promoter elements in gene expression.

Evolving noisy oscillatory dynamics in genetic regulatory networks

We introduce a genetic programming (GP) approach for evolving genetic networks that demonstrate desired dynamics when simulated as a discrete stochastic process. Our representation of genetic networks is based on a biochemical reaction model including key elements such as transcription, translation and post-translational modifications. The stochastic, reaction-based GP system is similar but not identical with algorithmic chemistries. We evolved genetic networks with noisy oscillatory dynamics. The results show the practicality of evolving particular dynamics in gene regulatory networks when modelled with intrinsic noise.

Multilevel modelling for inference of genetic regulatory networks

Time-course experiments with microarrays are often used to study dynamic biological systems and genetic regulatory networks (GRNs) that model how genes influence each other in cell-level development of organisms. The inference for GRNs provides important insights into the fundamental biological processes such as growth and is useful in disease diagnosis and genomic drug design. Due to the experimental design, multilevel data hierarchies are often present in time-course gene expression data. Most existing methods, however, ignore the dependency of the expression measurements over time and the correlation among gene expression profiles. Such independence assumptions violate regulatory interactions and can result in overlooking certain important subject effects and lead to spurious inference for regulatory networks or mechanisms. In this paper, a multilevel mixed-effects model is adopted to incorporate data hierarchies in the analysis of time-course data, where temporal and subject effects are both assumed to be random. The method starts with the clustering of genes by fitting the mixture model within the multilevel random-effects model framework using the expectation-maximization (EM) algorithm. The network of regulatory interactions is then determined by searching for regulatory control elements (activators and inhibitors) shared by the clusters of co-expressed genes, based on a time-lagged correlation coefficients measurement. The method is applied to two real time-course datasets from the budding yeast (Saccharomyces cerevisiae) genome. It is shown that the proposed method provides clusters of cell-cycle regulated genes that are supported by existing gene function annotations, and hence enables inference on regulatory interactions for the genetic network.