Precise sequence complementarity between yeast chromosome ends and two classes of just-subtelomeric sequences

The terminal regions (last 20 kb) of Saccharomyces cerevisiae chromosomes universally contain blocks of precise sequence similarity to other chromosome terminal regions. The left and right terminal regions are distinct in the sense that the sequence similarities between them are reverse complements. Direct sequence similarity occurs between the left terminal regions and also between the right terminal regions, but not between any left ends and right ends. With minor exceptions the relationships range from 80% to 100% match within blocks. The regions of similarity are composites of familiar and unfamiliar repeated sequences as well as what could be considered “single-copy” (or better “two-copy”) sequences. All terminal regions were compared with all other chromosomes, forward and reverse complement, and 768 comparisons are diagrammed. It appears there has been an extensive history of sequence exchange or copying between terminal regions. The subtelomeric sequences fall into two classes. Seventeen of the chromosome ends terminate with the Y′ repeat, while 15 end with the 800-nt “X2” repeats just adjacent to the telomerase simple repeats. The just-subterminal repeats are very similar to each other except that chromosome 1 right end is more divergent.

Characteristic enrichment of DNA repeats in different genomes

Using computer programs developed for this purpose, we searched for various repeated sequences including inverted, direct tandem, and homopurine–homopyrimidine mirror repeats in various prokaryotes, eukaryotes, and an archaebacterium. Comparison of observed frequencies with expectations revealed that in bacterial genomes and organelles the frequency of different repeats is either random or enriched for inverted and/or direct tandem repeats. By contrast, in all eukaryotic genomes studied, we observed an overrepresentation of all repeats, especially homopurine–homopyrimidine mirror repeats. Analysis of the genomic distribution of all abundant repeats showed that they are virtually excluded from coding sequences. Unexpectedly, the frequencies of abundant repeats normalized for their expectations were almost perfect exponential functions of their size, and for a given repeat this function was indistinguishable between different genomes.

Mutational inactivation of the proapoptotic gene BAX confers selective advantage during tumor clonal evolution

A remarkable instability at simple repeated sequences characterizes gastrointestinal cancer of the microsatellite mutator phenotype (MMP). Mutations in the DNA mismatch repair gene family underlie the MMP, a landmark for hereditary nonpolyposis colorectal cancer. These tumors define a distinctive pathway for carcinogenesis because they display a particular spectrum of mutated cancer genes containing target repeats for mismatch repair deficiency. One such gene is BAX, a proapoptotic member of the Bcl-2 family of proteins, which plays a key role in programmed cell death. More than half of colon and gastric cancers of the MMP contain BAX frameshifts in a (G)8 mononucleotide tract. However, the functional significance of these mutations in tumor progression has not been established. Here we show that inactivation of the wild-type BAX allele by de novo frameshift mutations confers a strong advantage during tumor clonal evolution. Tumor subclones with only mutant alleles frequently appeared after inoculation into nude mice of single-cell clones of colon tumor cell lines with normal alleles. In contrast, no clones of BAX-expressing cells were found after inoculation of homozygous cell clones without wild-type BAX. These results support the interpretation that BAX inactivation contributes to tumor progression by providing a survival advantage. In this context, survival analyses show that BAX mutations are indicators of poor prognosis for both colon and gastric cancer of the MMP.

RIGS (repeat-induced gene silencing) in Arabidopsis is transcriptional and alters chromatin configuration.

We have previously reported repeat-induced gene silencing (RIGS) in Arabidopsis, in which transgene expression may be silenced epigenetically when repeated sequences are present. Among an allelic series of lines comprising a primary transformant and various recombinant progeny carrying different numbers of drug resistance gene copies at the same locus, silencing was found to depend strictly on repeated sequences and to correlate with an absence of steady-state mRNA. We now report characterization, in nuclei isolated from the same transgenic lines, of gene expression by nuclear run-on assay and of chromatin structure by nuclease protection assay. We find that silencing is correlated with absence of run-on transcripts, indicating that expression is silenced at the level of transcription. We find further that silencing is also correlated with increased resistance to both DNase I and micrococcal nuclease, indicating that the silenced state reflects a change in chromatin configuration. We propose that silencing results when a locally paired region of homologous repeated nucleotide sequences is flanked by unpaired heterologous DNA, which leads chromatin to adopt a local configuration that is difficult to transcribe, and possibly akin to heterochromatin.

The anchorage function of CipA (CelL), a scaffolding protein of the Clostridium thermocellum cellulosome.

Enzymatic cellulose degradation is a heterogeneous reaction requiring binding of soluble cellulase molecules to the solid substrate. Based on our studies of the cellulase complex of Clostridium thermocellum (the cellulosome), we have previously proposed that such binding can be brought about by a special "anchorage subunit." In this "anchor-enzyme" model, CipA (a major subunit of the cellulosome) enhances the activity of CelS (the most abundant catalytic subunit of the cellulosome) by anchoring it to the cellulose surface. We have subsequently reported that CelS contains a conserved duplicated sequence at its C terminus and that CipA contains nine repeated sequences with a cellulose binding domain (CBD) in between the second and third repeats. In this work, we reexamined the anchor-enzyme mechanism by using recombinant CelS (rCelS) and various CipA domains, CBD, R3 (the repeat next to CBD), and CBD/R3, expressed in Escherichia coli. As analyzed by non-denaturing gel electrophoresis, rCelS, through its conserved duplicated sequence, formed a stable complex with R3 or CBD/R3 but not with CBD. Although R3 or CBD alone did not affect the binding of rCelS to cellulose, such binding was dependent on CBD/R3, indicating the anchorage role of CBD/R3. Such anchorage apparently increased the rCelS activity toward crystalline cellulose. These results substantiate the proposed anchor-enzyme model and the expected roles of individual CipA domains and the conserved duplicated sequence of CelS.

Evolution of gilled mushrooms and puffballs inferred from ribosomal DNA sequences

Homobasidiomycete fungi display many complex fruiting body morphologies, including mushrooms and puffballs, but their anatomical simplicity has confounded efforts to understand the evolution of these forms. We performed a comprehensive phylogenetic analysis of homobasidiomycetes, using sequences from nuclear and mitochondrial ribosomal DNA, with an emphasis on understanding evolutionary relationships of gilled mushrooms and puffballs. Parsimony-based optimization of character states on our phylogenetic trees suggested that strikingly similar gilled mushrooms evolved at least six times, from morphologically diverse precursors. Approximately 87% of gilled mushrooms are in a single lineage, which we call the “euagarics.” Recently discovered 90 million-year-old fossil mushrooms are probably euagarics, suggesting that (i) the origin of this clade must have occurred no later than the mid-Cretaceous and (ii) the gilled mushroom morphology has been maintained in certain lineages for tens of millions of years. Puffballs and other forms with enclosed spore-bearing structures (Gasteromycetes) evolved at least four times. Derivation of Gasteromycetes from forms with exposed spore-bearing structures (Hymenomycetes) is correlated with repeated loss of forcible spore discharge (ballistospory). Diverse fruiting body forms and spore dispersal mechanisms have evolved among Gasteromycetes. Nevertheless, it appears that Hymenomycetes have never been secondarily derived from Gasteromycetes, which suggests that the loss of ballistospory has constrained evolution in these lineages.

Thermal adaptation analyzed by comparison of protein sequences from mesophilic and extremely thermophilic Methanococcus species

The genome sequence of the extremely thermophilic archaeon Methanococcus jannaschii provides a wealth of data on proteins from a thermophile. In this paper, sequences of 115 proteins from M. jannaschii are compared with their homologs from mesophilic Methanococcus species. Although the growth temperatures of the mesophiles are about 50°C below that of M. jannaschii, their genomic G+C contents are nearly identical. The properties most correlated with the proteins of the thermophile include higher residue volume, higher residue hydrophobicity, more charged amino acids (especially Glu, Arg, and Lys), and fewer uncharged polar residues (Ser, Thr, Asn, and Gln). These are recurring themes, with all trends applying to 83–92% of the proteins for which complete sequences were available. Nearly all of the amino acid replacements most significantly correlated with the temperature change are the same relatively conservative changes observed in all proteins, but in the case of the mesophile/thermophile comparison there is a directional bias. We identify 26 specific pairs of amino acids with a statistically significant (P < 0.01) preferred direction of replacement.

Effective inhibition of influenza virus production in cultured cells by external guide sequences and ribonuclease P

The polymerase (PB2) and nucleocapsid (NP) genes encoded by the genome of influenza virus are essential for replication of the virus. When synthetic genes that express RNAs for external guide sequences targeted to the mRNAs of the PB2 and NP genes are stably incorporated into mouse cells in tissue culture, infection of these cells with influenza virus is nonproductive. Endogenous RNase P cleaves the targeted influenza virus mRNAs when they are in a complex with the external guide sequences. Targeting two different mRNAs simultaneously inhibits viral particle production more efficiently than does targeting only one mRNA.

The ability of HIV type 1 to use CCR-3 as a coreceptor is controlled by envelope V1/V2 sequences acting in conjunction with a CCR-5 tropic V3 loop

Although infection by primary HIV type 1 (HIV-1) isolates normally requires the functional interaction of the viral envelope protein with both CD4 and the CCR-5 coreceptor, a subset of such isolates also are able to use the distinct CCR-3 receptor. By analyzing the ability of a series of wild-type and chimeric HIV-1 envelope proteins to mediate CCR-3-dependent infection, we have determined that CCR-3 tropism maps to the V1 and V2 variable region of envelope. Although substitution of the V1/V2 region of a CCR-3 tropic envelope into the context of a CCR-5 tropic envelope is both necessary and sufficient to confer CCR-3 tropism, this same substitution has no phenotypic effect when inserted into a CXCR-4 tropic HIV-1 envelope context. However, this latter chimera acquires both CCR-3 and CCR-5 tropism when a CCR-5 tropic V3 loop sequence also is introduced. These data demonstrate that the V1/2 region of envelope can, like the V3 loop region, encode a particular coreceptor requirement and suggest that a functional envelope:CCR-3 interaction may depend on the cooperative interaction of CCR-3 with both the V1/V2 and the V3 region of envelope.

Natural trans-splicing in carnitine octanoyltransferase pre-mRNAs in rat liver

Carnitine octanoyltransferase (COT) transports medium-chain fatty acids through the peroxisome. During isolation of a COT clone from a rat liver library, a cDNA in which exon 2 was repeated, was characterized. Reverse transcription-PCR amplifications of total RNAs from rat liver showed a three-band pattern. Sequencing of the fragments revealed that, in addition to the canonical exon organization, previously reported [Choi, S. J. et al. (1995) Biochim. Biophys. Acta 1264, 215–222], there were two other forms in which exon 2 or exons 2 and 3 were repeated. The possibility of this exonic repetition in the COT gene was ruled out by genomic Southern blot. To study the gene expression, we analyzed RNA transcripts by Northern blot after RNase H digestion of total RNA. Three different transcripts were observed. Splicing experiments also were carried out in vitro with different constructs that contain exon 2 plus the 5′ or the 3′ adjacent intron sequences. Our results indicate that accurate joining of two exons 2 occurs by a trans-splicing mechanism, confirming the potential of these structures for this process in nature. The trans-splicing can be explained by the presence of three exon-enhancer sequences in exon 2. Analysis by Western blot of the COT proteins by using specific antibodies showed that two proteins corresponding to the expected Mr are present in rat peroxisomes. This is the first time that a natural trans-splicing reaction has been demonstrated in mammalian cells.

Discrimination of the closely related A and D genomes of the hexaploid oat Avena sativa L.

A satellite DNA sequence, As120a, specific to the A-genome chromosomes in the hexaploid oat, Avena sativa L., was isolated by subcloning a fragment with internal tandem repeats from a plasmid, pAs120, that had been obtained from an Avena strigosa (As genome) genomic library. Southern and in situ hybridization showed that sequences with homology to sequences within pAs120 were dispersed throughout the genome of diploid (A and C genomes), tetraploid (AC genomes), and hexaploid (ACD genomes) Avena species. In contrast, sequences homologous to As120a were found in two A-genome species (A. strigosa and Avena longiglumis) and in the hexaploid A. sativa whereas this sequence was little amplified in the tetraploid Avena murphyi and was absent in the remaining A- and C-genome diploid species. In situ hybridization of pAs120a to hexaploid oat species revealed the distribution of elements of the As120a repeated family over both arms of 14 of 42 chromosomes of this species. By using double in situ hybridization with pAs120a and a C genome-specific probe, three sets of 14 chromosomes were revealed corresponding to the A, C, and D genomes of the hexaploid species. Simultaneous in situ hybridizations with pAs120a and ribosomal probes were used to assign the SAT chromosomes of hexaploid species to their correct genomes. This work reports a sequence able to distinguish between the closely related A and D genomes of hexaploid oats. This sequence offers new opportunities to analyze the relationships of Avena species and to explore the possible evolution of various polyploid oat species.

Three functional luciferase domains in a single polypeptide chain

We report a unique case of a gene containing three homologous and contiguous repeat sequences, each of which, after excision, cloning, and expression in Escherichia coli, is shown to code for a peptide catalyzing the same reaction as the native protein, Gonyaulax polyedra luciferase (Mr = 137). This enzyme, which catalyzes the light-emitting oxidation of a linear tetrapyrrole (dinoflagellate luciferin), exhibits no sequence similarities to other luciferases in databases. Sequence analysis also reveals an unusual evolutionary feature of this gene: synonymous substitutions are strongly constrained in the central regions of each of the repeated coding sequences.

Direct isolation of human transcribed sequences from yeast artificial chromosomes through the application of RNA fingerprinting

The identification of cDNA clones from genomic regions known to contain human genes is usually the rate-limiting factor in positional cloning strategies. We demonstrate here that human genes present on yeast artificial chromosomes (YACs) are transcribed in yeast host cells. We have used the arbitrarily primed RNA (RAP) fingerprinting method to identify human-specific, transcribed sequences from YACs located in the 13q12 chromosome region. By comparing the RAP fingerprints generated using defined, arbitrary primers from various fragmented YACs, megaYACs, and host yeast, we were able to identify and map 20 products transcribed from the human YAC inserts. This method, therefore, permits the simultaneous isolation and mapping of novel expressed sequences directly from whole YACs.

Strained DNA is kinked by low concentrations of Zn2+

A novel atomic force microscope with a magnetically oscillated tip has provided unprecedented resolution of small DNA fragments spontaneously adsorbed to mica and imaged in situ in the presence of divalent ions. Kinks (localized bends of average angle 78°) were observed in axially strained minicircles consisting of tandemly repeated d(A)5 and d(GGGCC[C]) sequences. The frequency of kinks in identical minicircles increased 4-fold in the presence of 1 mM Zn2+ compared with 1 mM Mg2+. Kinking persisted in mixed Mg2+/Zn2+ electrolytes until the Zn2+ concentration dropped below 100 μM, indicating that this type of kinking may occur under physiological conditions. Kinking appears to replace intrinsic bending, and statistical analysis shows that kinks are not localized within any single sequence element. A surprisingly small free energy is associated with kink formation.

Two sterol regulatory element-like sequences mediate up-regulation of caveolin gene transcription in response to low density lipoprotein free cholesterol

Caveolae form the terminus for a major pathway of intracellular free cholesterol (FC) transport. Caveolin mRNA levels in confluent human skin fibroblasts were up-regulated following increased uptake of low density lipoprotein (LDL) FC. The increase induced by FC was not associated with detectable change in mRNA stability, indicating that caveolin mRNA levels were mediated at the level of gene transcription. A total of 924 bp of 5′ flanking region of the caveolin gene were cloned and sequenced. The promoter sequence included three G+C-rich potential sterol regulatory elements (SREs), a CAAT sequence and a Sp1 consensus sequence. Deletional mutagenesis of individual SRE-like sequences indicated that of these two (at −646 and −395 bp) were essential for the increased transcription rates mediated by LDL-FC, whereas the third was inconsequential. Gel shift analysis of protein binding from nuclear extracts to these caveolin promoter DNA sequences, together with DNase I footprinting, confirmed nucleoprotein binding to the SRE-like elements as part of the transcriptional response to LDL-FC. A supershift obtained with antibody to SRE-binding protein 1 (SPEBP-1) indicated that this protein binds at −395 bp. There was no reaction at −395 bp with anti-Sp1 antibody nor with either antibody at −646 bp. The cysteine protease inhibitor N-acetyl-leu-leu-norleucinal (ALLN), which inhibits SREBP catabolism, superinhibited caveolin mRNA levels regardless of LDL-FC. This finding suggests that SREBP inhibits caveolin gene transcription in contrast to its stimulating effect on other promoters. The findings of this study are consistent with the postulated role for caveolin as a regulator of cellular FC homeostasis in quiescent peripheral cells, and the coordinate regulation by SREBP of FC influx and efflux.