Adaptation to environmental temperature is a major determinant of molecular evolutionary rates in archaea

Methods to infer the ancestral conditions of life are commonly based on geological and paleontological analyses. Recently, several studies used genome sequences to gain information about past ecological conditions taking advantage of the property that the G+C and amino acid contents of bacterial and archaeal ribosomal DNA genes and proteins, respectively, are strongly influenced by the environmental temperature. The adaptation to optimal growth temperature (OGT) since the Last Universal Common Ancestor (LUCA) over the universal tree of life was examined, and it was concluded that LUCA was likely to have been a mesophilic organism and that a parallel adaptation to high temperature occurred independently along the two lineages leading to the ancestors of Bacteria on one side and of Archaea and Eukarya on the other side. Here, we focus on Archaea to gain a precise view of the adaptation to OGT over time in this domain. It has been often proposed on the basis of indirect evidence that the last archaeal common ancestor was a hyperthermophilic organism. Moreover, many results showed the influence of environmental temperature on the evolutionary dynamics of archaeal genomes: Thermophilic organisms generally display lower evolutionary rates than mesophiles. However, to our knowledge, no study tried to explain the differences of evolutionary rates for the entire archaeal domain and to investigate the evolution of substitution rates over time. A comprehensive archaeal phylogeny and a non homogeneous model of the molecular evolutionary process allowed us to estimate ancestral base and amino acid compositions and OGTs at each internal node of the archaeal phylogenetic tree. The last archaeal common ancestor is predicted to have been hyperthermophilic and adaptations to cooler environments can be observed for extant mesophilic species. Furthermore, mesophilic species present both long branches and high variation of nucleotide and amino acid compositions since the last archaeal common ancestor. The increase of substitution rates observed in mesophilic lineages along all their branches can be interpreted as an ongoing adaptation to colder temperatures and to new metabolisms. We conclude that environmental temperature is a major factor that governs evolutionary rates in Archaea.     

Global variation in G+C content along vertebrate genome DNA. Possible correlation with chromosome band structures

The global, rather than local, variation in G+C content along the nuclear DNA sequences of various organisms was studied using GenBank sequence data. When long DNA sequences of the genomes of Escherichia coli and Saccharomyces cerevisiae were examined, the levels of their G+C content (G+C%) were found to be within a narrow range around that of the whole genome. The G+C% levels for sequences of vertebrate genomes, however, were found to cover a wide range, showing that their genome is a mosaic of sequences with different G+C% levels, in each of which the sequence is fairly homogeneous in its G+C% for a very long distance. Through surveying a human genetic map and GenBank DNA sequences, the global variations in G+C% along the human genome DNA were found to be correlated with chromosome band structures.     

Molecular cloning of DNA from specific chromosomal regions by microdissection and sequence-independent amplification of DNA

A method that permits the in vitro amplification and cloning of DNA dissected from specific regions of a chromosome and does not require prior knowledge of the DNA sequence is described. DNA from several different chromosomal loci in the Drosophila melanogaster genome has been isolated by this method. Although the procedure was developed to permit the isolation of DNA sequences from serial sections of a single microdissected polytene chromosome, it should be useful for obtaining DNA clones from specific regions of the nonpolytene chromosomes of other organisms as well.     

Analysis of invariant sequences in 266 complete genomes

To date, the complete genome sequences of more than 250 organisms have been determined. This information can now be used to determine whether there exist any invariant sequences that are conserved among all organisms, from bacteria to plants, animals, and humans. The existence of invariant sequences would strongly suggest that these sequences have been inherited unchanged from the last common ancestor of all life, and that they have essential functions. We have developed a new software program to identify invariant sequences conserved among the currently sequenced genomes and applied this analysis to the complete genome sequences of 266 organisms. We have identified 3 invariant DNA sequences longer than or equal to 11 bp and 6 invariant amino acid sequences longer than or equal to 6 aa. The longest invariant DNA sequence, AAGTCGTACAAGGT (15 bp), was found in the 16S/18S rRNA gene. Two 8 aa sequences, GHVDHGKT in IF2 and EF-Tu and DTPGHVDF in EF-G, were the longest invariant amino acid sequences detected. These sequences could be essential elements from the genome of the last common ancestor and may have remained unchanged throughout evolution.     

Isolating promoters of multigene family members from the polyploid sugarcane genome by PCR-based walking in BAC DNA

The availability of a wider range of promoters for regulated expression in valuable transgenic crops would benefit functional genomics studies and current biotechnology programs aimed at improved productivity. Polymerase chain reaction (PCR)-based genome walking techniques are commonly used to isolate promoters or 5' flanking genomic regions adjacent to known cDNA sequences in genomes that are not yet completely sequenced. However, these techniques are problematic when applied directly to DNA isolated from crops with highly complex and large genomes. An adaptor ligation-mediated PCR-based BAC genome walking method is described here for the efficient isolation of promoters of multigene family members, such as the putative defense and fiber biosynthesis DIRIGENT genes that are abundant in the stem of sugarcane, a species with a highly polyploid genome. The advantage of this method is the efficient and specific amplification of the target promoter using BAC genomic DNA as template for the adaptor ligation-mediated PCR walking.     

Molecular microbial diversity of an agricultural soil in Wisconsin.

A culture-independent survey of the soil microbial diversity in a clover-grass pasture in southern Wisconsin was conducted by sequence analysis of a universal clone library of genes coding for small-subunit rRNA (rDNA). A rapid and efficient method for extraction of DNA from soils which resulted in highly purified DNA with minimal shearing was developed. Universal small-subunit-rRNA primers were used to amplify DNA extracted from the pasture soil. The PCR products were cloned into pGEM-T, and either hypervariable or conserved regions were sequenced. The relationships of 124 sequences to those of cultured organisms of known phylogeny were determined. Of the 124 clones sequenced, 98.4% were from the domain Bacteria. Two of the rDNA sequences were derived from eukaryotic organelles. Two of the 124 sequences were of nuclear origin, one being fungal and the other a plant sequence. No sequences of the domain Archaea were found. Within the domain, Bacteria, three kingdoms were highly represented: the Proteobacteria (16.1%), the Cytophaga-Flexibacter-Bacteroides group (21.8%), and the low G+C-content gram-positive group (21.8%). Some kingdoms, such as the Thermotogales, the green nonsulfur group, Fusobacteria, and the Spirochaetes, were absent. A large number of the sequences (39.4%) were distributed among several clades that are not among the major taxa described by Olsen et al. (G.J. Olsen, C.R. Woese, and R. Overbeek, J. Bacteriol., 176:1-6, 1994). From the alignments of the sequence data, distance matrices were calculated to display the enormous microbial diversity found in this soil in two ways, as phylogenetic trees and as multidimensional-scaling plots.     

Fidelity of Thermococcus litoralis DNA polymerase (Vent) in PCR determined by denaturing gradient gel electrophoresis.

DNA synthesis fidelities of two thermostable DNA polymerases, Thermus aquaticus (Taq) and Thermococcus litoralis (Tli, also known as Vent), and a non-thermostable enzyme, a modified T7 DNA polymerase (Sequenase), were determined by analyzing polymerase chain reaction (PCR) products using denaturing gradient gel electrophoresis (DGGE). The error rates were 4.4, 8.9, and 2.4 x 10(-5) errors/bp for modified T7, Taq, and Tli polymerase, respectively. Reducing the nucleotide triphosphate concentration for Tli polymerase during PCR did not alter the fidelity. The ability of DGGE to detect a mutant present at several percent in a wild type population is related to the polymerase fidelity. To examine the sensitivity of mutant detection, human genomic DNA containing a 1% fraction of a known base pair substitution mutant was PCR-amplified with the three enzymes using primers that flank the mutant sequence. The PCR products were analyzed by DGGE. The signal from the mutant present at 1% was visible in the samples amplified with modified T7 and Tli polymerase, but the higher error rate of Taq polymerase did not permit visualization of the signal in DNA amplified with Taq polymerase.     

Real-time PCR-based method for the estimation of genome sizes

The fast and reliable estimation of the genome sizes of various species would allow for a systematic analysis of many organisms and could reveal insights into evolutionary processes. Many methods for the estimation of genome sizes have already been described. The classical methods are based on the determination of the phosphate content in the DNA backbone of total DNA isolated from a defined number of cells or on reassociation kinetics of high molecular weight genomic DNA (c(0)t assay). More recent techniques employ DNA-specific fluorescent dyes in flow cytometry analysis, image analysis or absorption cytometry after Feulgen staining. The method presented here is based on the absolute quantification of genetic elements in a known amount (mass) of genomic DNA by real-time quantitative PCR. The method was evaluated on three different eukaryotic species, Saccharomyces cerevisiae (12.1 Mb), Xiphophorus maculatus (550 Mb) and Homo sapiens sapiens (2.9 Gb), and found to be fast, highly accurate and reliable.     

A simple and efficient method to determine the terminal sequences of restriction fragments containing known sequences.

We present an improvement of the inverse PCR method for the determination of end sequences of restriction fragments containing unknown DNA sequences flanked by known segments. In this approach, a short "bridge" DNA is inserted during the self-ligation step of the inverse PCR technique. This bridge DNA acts as primer annealing sites for amplification and subsequent direct sequencing. Successive PCR amplifications enable selective amplification of the unknown sequences from a complex mixture. Unlike previously described methods, our method does not require special materials, such as synthetic adapters or biotinylated primers that must be prepared each time to adapt the target. Furthermore, no complex steps such as dephosphorylation or purification are needed. Our method can save time and reduce the cost of cloning unknown sequences; it is ideal for routine, rapid gene walking. We applied this method to a GC-rich bacterial genome and succeeded in determining the end sequences of a 4.5-kb fragment.     

Mitochondria and G-quadruplex evolution: an intertwined relationship

G-quadruplexes (G4s) are non-canonical structures formed in guanine (G)-rich sequences through stacked G tetrads by Hoogsteen hydrogen bonding. Several studies have demonstrated the existence of G4s in the genome of various organisms, including humans, and have proposed that G4s have a regulatory role in various cellular functions. However, little is known regarding the dissemination of G4s in mitochondria. In this review, we report the observation that the number of potential G4-forming sequences in the mitochondrial genome increases with the evolutionary complexity of different species, suggesting that G4s have a beneficial role in higher-order organisms. We also discuss the possible function of G4s in mitochondrial (mt)DNA and long noncoding (lnc)RNA and their role in various biological processes.     

A method for site-directed transplacement of in vitro altered DNA sequences in Rhizobium

Summary A method that directs transplacement of in vitro altered DNA sequences to substitute the corresponding wild-type DNA sequences at the original location in the Rhizobium genome is described. The NPTII gene of transposon Tn5 which confers resistance to several antibiotics in a wide variety of organisms is used as a selectable marker on the vector. To generate DNA fragment substitution, it is not essential to place a selectable genetic marker directly linked to the in vitro altered region of the DNA fragment. The method utilizes the vector pSUP201 that replicates in E. coli but neither in Rhizobium nor in many other host systems. DNA transplacement occurs by the integration of the in vitro altered DNA fragment at its homologous site in the chromosome along with the vector. A duplicated target sequence is generated, which is then followed by a homologous recombination event between the duplicated DNA sequences leading to the excision of the vector carrying the selectable marker together with one of the duplicated copies of DNA. The excised DNA fragment is subsequently lost due to the inability of the vector pSUP201 to replicate in Rhizobium. Using this method, we have transplaced a nif DNA fragment containing a large deletion of approximately 9 Kb into its homologous site in the genome of cowpea Rhizobium IRc78. The IRc78 strain carrying the deletion forms nodules on host-plants but is unable to reduce atmospheric nitrogen. The method allows precise alterations of the genomic DNA in an organism that is genetically not well characterized, provided a vector carrying a selectable marker can be transferred.Dedicated to Professor Georg Melchers to celebrate his 50-year association with the journal     

Maximum-likelihood estimation of phylogeny from DNA sequences when substitution rates differ over sites

Felsenstein's maximum-likelihood approach for inferring phylogeny from DNA sequences assumes that the rate of nucleotide substitution is constant over different nucleotide sites. This assumption is sometimes unrealistic, as has been revealed by analysis of real sequence data. In the present paper Felsenstein's method is extended to the case where substitution rates over sites are described by the gamma distribution. A numerical example is presented to show that the method fits the data better than do previous models.      

Applications of retrotransposons as genetic tools in plant biology

Retrotransposons are mobile genetic elements that accomplish transposition via an RNA intermediate that is reverse transcribed before integration into a new location within the host genome. They are ubiquitous in eukaryotic organisms and constitute a major portion of the nuclear genome (often more than half of the total DNA) in plants. Furthermore, they are dispersed as interspersed repetitive sequences throughout most of the length of all host chromosomes. These unique properties of retrotransposons have been exploited as genetic tools for plant genome analysis. Major applications are in determining phylogeny and genetic diversity and in the functional analyses of genes in plants. Here, recent advances in molecular markers, gene tagging and functional genomics technologies using plant retrotransposons are described.     

Thermus aquaticus DNA polymerase I mutants with altered fidelity. Interacting mutations in the O-helix

Phe(667) in the conserved O-helix of Thermus aquaticus (Taq) DNA polymerase I (pol I) is known to be important for discrimination against dideoxy-NTPs. We show here that Phe(667) is also important for base selection fidelity. In a forward mutation assay at high polymerase concentration, wild type pol I catalyzed frequent A --> T and G --> T transversions and -1 frameshifts at nonreiterated sites involving loss of a purine immediately downstream of a pyrimidine. The mutants F667L and A661E,I665T,F667L exhibited large decreases in A --> T and G --> T transversions, and the triple mutant displayed reduction in the aforementioned -1 frameshifts as well. Kinetic analysis showed that the F667L and A661E,I665T,F667L polymerases discriminated against synthesis of A:A mispairs more effectively and catalyzed less extension of A:A mispairs than the wild type enzyme. These data indicate that Phe(667) functions in maintaining the error frequency and spectrum, and the catalytic efficiency, of wild type pol I. We also found that the strong general mutator activity conferred by the single A661E substitution was entirely suppressed in the A661E, I665T,F667L polymerase, exemplifying how interactions among O-helix residues can contribute to fidelity. We discuss the mutator and anti-mutator mutations in light of recently obtained three-dimensional structures of T. aquaticus pol I.     

Efficiency of ligation-mediated PCR and TAIL-PCR methods for isolation of RbcS promoter sequences from green microalgae Ankistrodesmus convolutus

Isolation of promoter sequences from known gene sequences is a tedious task in genome-related research. An efficient method of obtaining the promoter sequences is necessary in order to successfully use targeted promoters for genetic manipulations. Here, efficiency and usefulness of two PCR-based methods, namely: ligation-mediated PCR and thermal asymmetric interlaced (TAIL) PCR, for isolation of promoter sequences of the ribulose-1,5-bisphosphate carboxylase/oxygenase small subunit (RbcS) gene from green microalgae Ankistrodesmus convolutus (A. convolutus) were evaluated. The results showed that the amplification efficiency of TAIL-PCR was higher than that of the ligation-mediated PCR method, i.e. the amplified promoter fragments of 1.2 and 0.8 kb in length or promoter sequences of 813 and 606 bp (after eliminating the unreadable sequences). The use of TAIL-PCR described here presents a low cost and efficient strategy for the isolation of promoter sequences of known genes, especially in GC-rich regions, and species with little or no available genome information such as A. convolutus.     

Pseudogenes, junk DNA, and the dynamics of Rickettsia genomes

Studies of neutrally evolving sequences suggest that differences in eukaryotic genome sizes result from different rates of DNA loss. However, very few pseudogenes have been identified in microbial species, and the processes whereby genes and genomes deteriorate in bacteria remain largely unresolved. The typhus-causing agent, Rickettsia prowazekii, is exceptional in that as much as 24% of its 1.1-Mb genome consists of noncoding DNA and pseudogenes. To test the hypothesis that the noncoding DNA in the R. prowazekii genome represents degraded remnants of ancestral genes, we systematically examined all of the identified pseudogenes and their flanking sequences in three additional Rickettsia species. Consistent with the hypothesis, we observe sequence similarities between genes and pseudogenes in one species and intergenic DNA in another species. We show that the frequencies and average sizes of deletions are larger than insertions in neutrally evolving pseudogene sequences. Our results suggest that inactivated genetic material in the Rickettsia genomes deteriorates spontaneously due to a mutation bias for deletions and that the noncoding sequences represent DNA in the final stages of this degenerative process.     

Multiple independent fusions of glucose-6-phosphate dehydrogenase with enzymes in the pentose phosphate pathway

Fusions of the first two enzymes in the pentose phosphate pathway, glucose-6-phosphate dehydrogenase (G6PD) and 6-phosphogluconolactonase (6PGL), have been previously described in two distant clades, chordates and species of the malarial parasite Plasmodium. We have analyzed genome and expressed sequence data from a variety of organisms to identify the origins of these gene fusion events. Based on the orientation of the domains and range of species in which homologs can be found, the fusions appear to have occurred independently, near the base of the metazoan and apicomplexan lineages. Only one of the two metazoan paralogs of G6PD is fused, showing that the fusion occurred after a duplication event, which we have traced back to an ancestor of choanoflagellates and metazoans. The Plasmodium genes are known to contain a functionally important insertion that is not seen in the other apicomplexan fusions, highlighting this as a unique characteristic of this group. Surprisingly, our search revealed two additional fusion events, one that combined 6PGL and G6PD in an ancestor of the protozoan parasites Trichomonas and Giardia, and another fusing G6PD with phosphogluconate dehydrogenase (6PGD) in a species of diatoms. This study extends the range of species known to contain fusions in the pentose phosphate pathway to many new medically and economically important organisms.     

Discrimination of reproductive forms of Thrips tabaci (Thysanoptera: Thripidae) by PCR with sequence specific primers

In agriculture, although it is important to identify species of pest insects, the morphological identification is often difficult. DNA genotyping is useful for the identification of species in morphologically indiscriminable species. Thrips tabaci (Lindeman) can be divided into two reproductive forms (arrhenotoky and thelytoky, each of which different in pesticide resistance) but morphological discrimination is not possible. Here, we establish a simple method to discriminate the strains based on their mitochondrial DNA sequences. Phylogenetic analysis including the T. tabaci and congeneric species provided ancestor sequences of each strain of T tabaci. Based on the ancestor sequences, we developed a primer set that include strain specific primers on sense strand and common primer on anti sense strand. Using this primer set, the strains of 196 individuals of T. tabaci were successfully assigned to each ofgenotypic forms. As the phylogeny and ancestor sequences were based on worldwide samples, this method will work well on most populations around the world.