The 5'-flanking sequence and regulatory elements of the cystatin S gene.

The gene encoding rat cystatin S (Cys S), a salivary gland-specific secretory protein, has CAAT and TATA boxes upstream of the inititation codon (Cox and Shaw, 1992), and contains regions that resemble those of other hormonally responsive eukaryotic genes. The 5'-flanking sequence of the rat Cys S gene has a potential CREB/AP-1 binding site (Rupp et al., 1990; Trejo et al., 1992), two potential glucocorticoid responsive elements (GREs, Drouin et al., 1989), and a possible GR/PR (glucocorticoid/progesterone) responsive element (Forman and Samuels, 1990). One of these potential GREs is adjacent to a potential AP-2 binding site, and another is typical of the glucocorticoid and progesterone receptor binding site. In this report, we have identified three regions in the 5'-flanking region of the Cys S gene that are found in salivary gland-specific genes (Ting et al., 1992) with a GT-rich region located between conserved elements II and III. Transfection experiments described in this paper suggest that a 281-bp DNA fragment from the Cys S gene promoter region with conserved elements II and III, the GT-rich region, and a possible GR/PR responsive element contains a negative regulatory element. In addition, our experiments suggest that the GT-rich region by itself is acting as a positive regulatory element.     

New Phytologist 140 (1998), 363–383

In the November 1998 issue of New Phytologist , we published the Tansley review 'Gibberellins: regulating genes and germination' by Sian Ritchie and Simon Gilroy ( New Phytol. (1998) 140 , 363–383). Since its publication, it has come to our attention that text associated with Fig. 4 was omitted during production. The correct figure is reprinted here in full.
We apologise to the author and to our readers for this mistake.
Figure 4. Promoter sequences of various genes expressed in the cereal aleurone and shown to be regulated by GA. The position of each sequence is indicated relative to the start codon. Regions identified as being involved in regulation of the genes are highlighted, as are similar regions in other genes. Sites at which protein has been shown to bind are also indicated. ( a ) Barley Amy 32b (Sutcliff et al ., 1993; Whittier et al ., 1987); wheat Amy 2/54 (Huttley et al ., 1992; Rushton et al ., 1992; Rushton et al ., 1995); barley Amy 46 (Khursheed & Rogers, 1988); barley Amy 2/p155 (Knox et al ., 1987); barley aleurain (Whittier et al ., 1987); barley β-glucanase II (Wolf, 1992); wheat cathepsin B-like (Cejudo et al ., 1992); rice ubiquitin-conjugating enzyme (Chen et al ., 1995). ( b ). Wheat Amy 1/18 (Rushton et al ., 1992); barley Amy pHV 19 (Jacobsen & Close, 1991; Gubler & Jacobsen, 1992)/ Amy 1 / 6-4 (Khursheed & Rogers, 1988; Rogers, Lanahan & Rogers 1994); rice OSamy-a / Amy 3c (Ou-Lee et al ., 1988; Sutcliff et al ., 1991; Yu et al ., 1992; Goldman et al ., 1994); rice Amy 3B (Sutcliffe et al ., 1991); rice OSamy-c (Kim et al ., 1992; Kim & Wu, 1992; Tanida et al ., 1994); rice Amy 1A (Huang et al ., 1990; Itoh et al ., 1995).
Figure 4 ( b ). For legend see facing page.     

Polymorphisms and evolutionary history of retrotransposon insertions in rice promoters

Retrotransposons are ubiquitous in higher plant genomes. The presence or absence of retrotransposons in whole genome and high throughput genomic sequence (HTGS) from cultivated and wild rice was investigated to understand the organization and evolution of retrotransposon insertions in promoter regions. Approximately half of the Oryza sativa subsp. japonica 'Nipponbare' promoters with retrotransposons conserved in Oryza sativa subsp. indica '93-11' and four wild rice species showed higher sequence conservation in retrotransposon than nonretrotransposon regions. We further investigated, in detail, the evolutionary dynamics of five retrotransposons in the promoter regions of 95 rice genotypes. Our data suggest that four of five insertions (Rp2-Rp5) occurred in the ancestor of AA genome, while the other insertion (Rp1) predates the ancestral divergence of Oryza officinalis (CC genome). Four retrotransposons (Rp2-Rp5) were present in 52% (Rp2), 29% (Rp3), 53% (Rp4), and 43% (Rp5) of the rice genotypes with AA genome type, and the fifth retrotransposon (Rp1) was present in 95% of the rice genotypes with AA, BBCC, or CC genome types. Furthermore, most of these retrotransposons were found to evolve slower than flanking promoter regions, suggesting a role in promoter function for regulating downstream genes.     

Tricross : using dot-plots in sequence-id space to detect uncataloged intergenic features.

MOTIVATION: The process of determining the functional sequence content of an organism is confounded by several factors. Large protein coding sequences are relatively easy to find by statistical methods. Smaller proteins however may escape detection due to their size falling below some arbitrary researcher-defined minimum cutoff, or the inability to precisely define a promoter, or translational start (Delcher et al., Nucleic Acids Res., 27, 4636-4641, 1999). Promoter and regulatory sequences themselves are difficult to define due to a significant amount of allowable sequence variation, as well as a probable lack of any completely accurate whole-organismal gene catalogs to date. Finally, certain genes coding functional RNAs may have insufficient structural or sequence constraints to be detectable by normal sequence structure/pattern searching methods (Eddy and Rivas, Bioinformatics, 16, 583-605, 2000). In those cases where there are multiple closely related organisms that have been sequenced, there is additional information that may be used in the investigation of sequence content-that being the possible conserved nature of functional sequences between the organisms. We present a method for the utilization of this conserved information to detect genes and other potentially functional sequences that may be missed by standard ORF-calling, RNA finding, and pattern matching software. The tricross programs produce a multi-way cross comparison of three sets of sequences, determine which are conserved in all three sets, and produce a graphical (Virtual Reality Modelling Language-VRML; (ISO/IEC 14772-1: 1997, VDC), 1997) representation as well as alignments of all sequence triples found. The software can also be applied to a pair of sequence sets, though the noise in the results increases. RESULTS: Tricross has been used to examine the intergenic-sequence content of the three archaeal Pyrococcus genomes to determine the most highly related sequences remaining between the annotated protein and RNA coding sequences. Set to relatively stringent similarity requirements for the search, tricross found 101 intergenic sequences conserved among the three organisms. Interestingly, 29 of these appear to contain members of a family of small RNA molecules (Kiss-Laszlo et al., EMBO J., 17, 797-807, 1998) only recently discovered in the Archaea (Armbruster, OSU, Diss., 1988; Omer et al., Science, 288, 517-522, 2000; Gaspin et al., J. Mol. Biol., 297, 895-906, 2000). While some of the remaining 72 appear to be individual highly conserved promoter sequences, others have no currently known biological significance. Although originally developed to facilitate the examination of intergenic sequences, none of the tricross logic is inherently specific to intergenic sequences. The software can also be applied to gene sequences, and has been used to produce inter-genomic gene order dot-plots for Haemophilus influenzae (Fleischmann et al., Science, 269, 496-512, 1995) versus H.ducreyi (unpublished data), and Neisseria meningiditis Z2491 (serogroup A) (Parkhill et al., Nature, 404, 502-506, 2000) versus Neisseria meningiditis Z58 (serogroup B) (Tettelin et al., Science, 287, 1809-1815, 2000) versus Neisseria gonorrhoeae (Lewis et al., http://micro-gen.ouhsc.edu/, 2000). AVAILABILITY: The tricross software package is available from http://www.biosci.ohio-state.edu/~ray/bioinformatics/tricross.html. CONTACT: ray@biosci.ohio-state.edu; daniels.7@osu.edu; munsonr@pediatrics.ohio-state.edu Supplementary information: Additional data from the cross-genomic comparisons examined in the discussion section are linked from http://www.biosci.ohio-state.edu/~ray/bioinformatics/tricross.html.     

Rapid genome evolution revealed by comparative sequence analysis of orthologous regions from four triticeae genomes

Bread wheat (Triticum aestivum) is an allohexaploid species, consisting of three subgenomes (A, B, and D). To study the molecular evolution of these closely related genomes, we compared the sequence of a 307-kb physical contig covering the high molecular weight (HMW)-glutenin locus from the A genome of durum wheat (Triticum turgidum, AABB) with the orthologous regions from the B genome of the same wheat and the D genome of the diploid wheat Aegilops tauschii (Anderson et al., 2003; Kong et al., 2004). Although gene colinearity appears to be retained, four out of six genes including the two paralogous HMW-glutenin genes are disrupted in the orthologous region of the A genome. Mechanisms involved in gene disruption in the A genome include retroelement insertions, sequence deletions, and mutations causing in-frame stop codons in the coding sequences. Comparative sequence analysis also revealed that sequences in the colinear intergenic regions of these different genomes were generally not conserved. The rapid genome evolution in these regions is attributable mainly to the large number of retrotransposon insertions that occurred after the divergence of the three wheat genomes. Our comparative studies indicate that the B genome diverged prior to the separation of the A and D genomes. Furthermore, sequence comparison of two distinct types of allelic variations at the HMW-glutenin loci in the A genomes of different hexaploid wheat cultivars with the A genome locus of durum wheat indicates that hexaploid wheat may have more than one tetraploid ancestor.     

Isolation and analysis of cDNAs encoding small GTP-binding proteins of Arabidopsis thaliana.

We previously isolated a DNA fragment from Arabidopsis thaliana homologous to the mammalian ras gene and named it ara [Matsui et al., Gene 76 (1989) 313-319]. Screening of cDNA clones homologous to ara in A. thaliana resulted in the isolation of four homologous genes. The products of these genes, ARA-2, ARA-3, ARA-4 and ARA-5, showed conservation of amino acids (aa) in four regions, all of which are present in small GTP-binding proteins, and are important for GTPase/GTP-binding activities. These products were highly homologous to those of the YPT genes of Saccharomyces cerevisiae and the ypt gene of Schizosaccharomyces pombe in the regions around aa 45, which is thought to be the site interacting with effector molecules. The products of these four genes showed characteristic aa sequence at their C termini, Cys-Cys-Xaa-Xaa. Another characteristic of this family is presence of Ser in place of Gly in the first conserved region (Gly12 of mammalian GTP-binding Ras protein).     

Adaptive evolution of newly emerged micro-RNA genes in Drosophila

How often micro-RNA (miRNA) genes emerged and how fast theyevolved soon after their emergence are some of the central questionsin the evolution of miRNAs. Because most known miRNA genes areancient and highly conserved, these questions can be best answeredby identifying newly emerged miRNA genes. Among the 78 miRNAgenes in Drosophila reported before 2007, only 5 are confirmedto be newly emerged in the genus (although many more can befound in the newly reported data set; e.g., Ruby et al. 2007;Stark et al. 2007; Lu et al. 2008). These new miRNA genes haveundergone numerous changes, even in the normally invariant maturesequences. Four of them (the miR-310/311/312/313 cluster, denotedmiR-310s) were duplicated from other conserved miRNA genes.The fifth one (miR-303) appears to be a very young gene, originatingde novo from a non-miRNA sequence recently. We sequenced these5 miRNA genes and their neighboring regions from a worldwidecollection of Drosophila melanogaster lines. The levels of divergenceand polymorphism in these miRNA genes, vis-à-vis thoseof the neighboring DNA sequences, suggest that these 5 genesare evolving adaptively. Furthermore, the polymorphism patternof miR-310s in D. melanogaster is indicative of hitchhikingunder positive selection. Thus, a large number of adaptive changesover a long period of time may be essential for the evolutionof newly emerged miRNA genes.     

Overexpression of rice OSH genes induces ectopic shoots on leaf sheaths of transgenic rice plants

Five rice homeobox (OSH) genes were overexpressed under the control of the cauliflower mosaic virus 35S promoter or the rice actin gene promoter in transgenic rice plants. Almost all of the transgenic plants showed abnormal phenotypes, which could be classified into three types according to their severity. Plants with the most severe phenotype formed only green organs, with many shoot apices on their adaxial sides. Plants with an intermediate phenotype formed bladeless leaves with normally developed leaf sheaths. Plants with a mild phenotype formed normal leaf sheaths and blades, but lacked ligules and showed diffusion of the blade-sheath boundary. The leaf structure of this phenotype was similar to that of dominant maize mutants, such as Kn1, Rs1, Lg3, and Lg4. Based on these phenotypes, we suggest that ectopic expression of the rice OSH genes interferes with the development of leaf blades and maintains leaves in less differentiated states. These results are discussed in relation to the leaf maturation schedule hypothesis (M. Freeling et al., 1992, BioEssays 14, 227-236).     

Detailed analysis of the mouse H-2Kb promoter: enhancer-like sequences and their role in the regulation of class I gene expression

Sequencing and deletion analyses of the H-2Kb promoter have suggested that several regions may be important for expression and regulation of this gene. Two of these regions are conserved inside the promoter of several genes coding for classical transplantation antigens, but not in the promoter of class I genes located in the Qa region. They display enhancer-like activity in cells that express H-2 genes, but show some tissue specificity in that they function very poorly in undifferentiated embryonal carcinoma cells in which H-2 genes are not expressed. They also have been shown not to be the target of the adenovirus-12 induced repression of class I gene expression recently demonstrated by Schrier et al. The promoter of the beta 2-microglobulin gene also contains a sequence with enhancer-like activity, but shares no homology with the H-2Kb promoter region.     

Analysis of regulatory elements and genes required for carbon tetrachloride degradation in Pseudomonas stutzeri strain KC

Previously, we described the generation and initial characterization of four Tn5 mutants of Pseudomonas stutzeri strain KC with impaired ability to degrade carbon tetrachloride (Sepúlveda-Torres et al., 1999). In this study, we show cloning and sequencing of an 8.3 kbp region in which all four transposons were located. This fragment encodes eight potential genes and is located in the central part of the 25 kbp fragment recently identified by Lewis et al. (2000) and shown by them to be sufficient to confer carbon tetrachloride transformation capability upon other pseudomonads. The four transposon insertion mutants mapped in ORF's F and I designated by Lewis et al. (2000). This is consistent with the results by Lewis et al. (2000) that orfFis required for carbon tetrachloride degradation. We further established that orfl is required for CCl4 degradation since the three mutants in this ORF were unable to degrade carbon tetrachloride. We present our analysis of the gene and protein sequences from the 8.3 kbp region and propose a tentative model for the role of different genes in the synthesis and activity of pyridine-2,6-bis(thiocarboxylate) (PDTC), the secreted factor responsible for carbon tetrachloride dechlorination. We also found a putative promoter that overlaps with a Fur-box-like sequence in the region upstream of mutated genes. To test this putative promoter region and Fur-box, we generated and ligated DNA fragments containing wild-type and mutant Fur-boxes to a lacZ reporter. The wild-type fragment showed promoter activity that is regulated by the concentration of iron in the medium. Finally, we screened a selection of Pseudomonas strains, including P. putida DSMZ 3601--a strain known to produce PDTC--for the presence of the genes characterized in this study. None of the strains tested positive, suggesting that Pseudomonas stutzeri strain KC may possess a distinct biosynthetic pathway for PDTC production.     

Sequence of the gene (pheA) encoding phenol monooxygenase from Pseudomonas sp. EST1001: expression in Escherichia coli and Pseudomonas putida.

The plasmid pEST1412 contains the genes, pheA and pheB, encoding phenol monooxygenase (PMO) and catechol 1,2-dioxygenase (C12]), respectively. Thse were originally cloned from the plasmid DNA of Pseudomonas sp. EST1001 [Kivisaar et al., Plasmid 24 (1990) 25-36]. Although pheA and pheB are cotranscribed using the promoter sequences derived from Tn4652 and the level of expression of C120 activities from pEST1412 was equal both in Escherichia coli and in Pseudomonas putida, the level of PMO activity measured in the cell-free extracts of E. coli was lower than that in P. putida. The nucleotide sequence of the 2.0-kb PstI-HindIII fragment of pEST1412 carrying pheA was determined. A 1821-bp ORF was found in this DNA. The structural gene (tfdB) encoding 2,4-dichlorophenol hydroxylase from pJP4 has been sequenced [Perkins et al., J. Bacteriol. 172 (1990) 2351-2359]. Comparison of the deduced amino acid sequences of tfdB and pheA revealed highly conserved regions in the protein products of these genes.     

Xenopus borealis and Xenopus laevis 28S ribosomal DNA and the complete 40S ribosomal precursor RNA coding units of both species.

We have determined the nucleotide sequence of Xenopus borealis 28S ribosomal DNA (rDNA) and have revised the sequence of Xenopus laevis 28S rDNA (Ware et al., Nucl. Acids Res. 11, 7795-7817 (1983)). In the regions encoding the conserved structural core of 28S rRNA (2490 nucleotides) there are only four differences between the two species, each difference being a base substitution. In the variable regions, also called eukaryotic expansion segments (ca. 1630 nucleotides) there are some 61 differences, due to substitutions, mini-insertions and mini-deletions. Thus, evolutionary divergence in the variable regions has been at least 20-fold more rapid than in the conserved core. A search for intraspecies sequence variation has revealed minimal heterogeneity in X. laevis and none in X. borealis. At three out of four sites where heterogeneity was found in X. laevis (all in variable regions) the minority variant corresponded to the standard form in X. borealis. Intraspecies heterogeneity and interspecies divergence in the 28S variable regions are much less extensive than in the transcribed spacers. The 28S sequences are from the same clones that were used previously for sequencing the 18S genes and transcribed spacers. The complete sequences of the 40S precursor regions of the two reference clones are given.     

Cloning and characterization of a third type of human alpha-amylase gene, AMY2B

We have previously reported concerning the existence of a third type of human alpha-amylase gene, AMY3 [Emi et al., Gene 62 (1988) 229-235; Tomita et al., Gene 76 (1989) 11-18], which is expressed in a lung carcinoid tissue, and differs in nucleotide sequence from the two previously characterized human alpha-amylase genes coding for salivary and pancreatic isozymes, termed AMY1 and AMY2, respectively. Here, we rename this gene AMY2B to coincide with the designation by Gumucio et al. [Mol. Cell Biol. 8 (1988) 1197-1205] and describe its genetic properties as revealed by sequencing studies. It consists of ten major exons whose sequences are highly homologous to those of AMY1 and AMY2. Not only the exons, but also most of the introns seem to be highly conserved, as judged from physical mapping data. The AMY2B gene identified from mRNA in a lung carcinoid tissue has at least two additional untranslated exons in its 5' region; hence the promoter lies far upstream relative to the other two AMY genes.     

Isolation and sequencing of the 5′ end of the rat microtubule-associated protein (MAP1B)-encoding cDNA

We have isolated and sequenced the 5′ end of the cDNA encoding the rat microtubule-associated protein 1B (MAP1B). We found that this region is highly homologous to the corresponding regions of the human [Lien et al., 22 (1994) 273–280] and mouse [Noble et al., J. Cell Biol. 109 (1989) 3367–3376] MAPIB genes. The combination of the sequence that we are presenting with the previously published sequence [Zauner et al., Eur. J. Cell Biol. 57 (1992) 66–74], represents the complete rat MAP1B cDNA coding sequence.