Substitution rate comparisons between grasses and palms: synonymous rate differences at the nuclear gene Adh parallel rate differences at the plastid gene rbcL.

A number of studies have noted that nucleotide substitution rates at the chloroplast-encoded rbcL locus violate the molecular clock principle. Substitution rate variation at this plastid gene is particularly pronounced between palms and grasses; for example, a previous study estimated that substitution rates in rbcL sequences are approximately 5-fold faster in grasses than in palms. To determine whether a proportionate change in substitution rates also occurs in plant nuclear genes, we characterized nucleotide substitution rates in palm and grass sequences for the nuclear gene Adh. In this article, we report that palm sequences evolve at a rate of 2.61 x 10(-9) substitution per synonymous site per year, a rate which is slower than most plant nuclear genes. Grass Adh sequences evolve approximately 2.5-fold faster than palms at synonymous sites. Thus, synonymous rates in nuclear Adh genes show a marked decrease in palms relative to grasses, paralleling the pattern found at the plastid rbcL locus. This shared pattern indicates that synonymous rates are correlated between a nuclear and a plastid gene. Remarkably, nonsynonymous rates do not show this correlation. Nonsynonymous rates vary between two duplicated grass Adh loci, and nonsynonymous rates at the palm Adh locus are not markedly reduced relative to grasses.

Substitution rate calibration of small subunit ribosomal RNA identifies chlorarachniophyte endosymbionts as remnants of green algae.

Chlorarachniophytes are amoeboid algae with chlorophyll a and b containing plastids that are surrounded by four membranes instead of two as in plants and green algae. These extra membranes form important support for the hypothesis that chlorarachniophytes have acquired their plastids by the ingestion of another eukaryotic plastid-containing alga. Chlorarachniophytes also contain a small nucleus-like structure called the nucleomorph situated between the two inner and the two outer membranes surrounding the plastid. This nucleomorph is a remnant of the endosymbiont's nucleus and encodes, among other molecules, small subunit ribosomal RNA. Previous phylogenetic analyses on the basis of this molecule provided unexpected and contradictory evidence for the origin of the chlorarachniophyte endosymbiont. We developed a new method for measuring the substitution rates of the individual nucleotides of small subunit ribosomal RNA. From the resulting substitution rate distribution, we derived an equation that gives a more realistic relationship between sequence dissimilarity and evolutionary distance than equations previously available. Phylogenetic trees constructed on the basis of evolutionary distances computed by this new method clearly situate the chlorarachniophyte nucleomorphs among the green algae. Moreover, this relationship is confirmed by transversion analysis of the Chlorarachnion plastid small subunit ribosomal RNA.

Malaria’s Eve: Evidence of a recent population bottleneck throughout the world populations of Plasmodium falciparum

We have analyzed DNA sequences from world-wide geographic strains of Plasmodium falciparum and found a complete absence of synonymous DNA polymorphism at 10 gene loci. We hypothesize that all extant world populations of the parasite have recently derived (within several thousand years) from a single ancestral strain. The upper limit of the 95% confidence interval for the time when this most recent common ancestor lived is between 24,500 and 57,500 years ago (depending on different estimates of the nucleotide substitution rate); the actual time is likely to be much more recent. The recent origin of the P. falciparum populations could have resulted from either a demographic sweep (P. falciparum has only recently spread throughout the world from a small geographically confined population) or a selective sweep (one strain favored by natural selection has recently replaced all others). The selective sweep hypothesis requires that populations of P. falciparum be effectively clonal, despite the obligate sexual stage of the parasite life cycle. A demographic sweep that started several thousand years ago is consistent with worldwide climatic changes ensuing the last glaciation, increased anthropophilia of the mosquito vectors, and the spread of agriculture. P. falciparum may have rapidly spread from its African tropical origins to the tropical and subtropical regions of the world only within the last 6,000 years. The recent origin of the world-wide P. falciparum populations may account for its virulence, as the most malignant of human malarial parasites.

Determining divergence times with a protein clock: Update and reevaluation

A recent study of the divergence times of the major groups of organisms as gauged by amino acid sequence comparison has been expanded and the data have been reanalyzed with a distance measure that corrects for both constraints on amino acid interchange and variation in substitution rate at different sites. Beyond that, the availability of complete genome sequences for several eubacteria and an archaebacterium has had a great impact on the interpretation of certain aspects of the data. Thus, the majority of the archaebacterial sequences are not consistent with currently accepted views of the Tree of Life which cluster the archaebacteria with eukaryotes. Instead, they are either outliers or mixed in with eubacterial orthologs. The simplest resolution of the problem is to postulate that many of these sequences were carried into eukaryotes by early eubacterial endosymbionts about 2 billion years ago, only very shortly after or even coincident with the divergence of eukaryotes and archaebacteria. The strong resemblances of these same enzymes among the major eubacterial groups suggest that the cyanobacteria and Gram-positive and Gram-negative eubacteria also diverged at about this same time, whereas the much greater differences between archaebacterial and eubacterial sequences indicate these two groups may have diverged between 3 and 4 billion years ago.

Replication infidelity during a single cycle of Ty1 retrotransposition.

Retroviruses undergo a high frequency of genetic alterations during the process of copying their RNA genomes. However, little is known about the replication fidelity of other elements that transpose via reverse transcription of an RNA intermediate. The complete sequence of 29 independently integrated copies of the yeast retrotransposon Ty1 (173,043 nt) was determined, and the mutation rate during a single cycle of replication was calculated. The observed base substitution rate of 2.5 x 10(-5) bp per replication cycle suggests that this intracellular element can mutate as rapidly as retroviruses. The pattern and distribution of errors in the Ty1 genome is nonrandom and provides clues to potential in vivo molecular mechanisms of reverse transcriptase-mediated error generation, including heterogeneous RNase H cleavage of Ty1 RNA, addition of terminal nontemplated bases, and transient dislocation and realignment of primer-templates. Overall, analysis of errors generated during Ty1 replication underscores the utility of a genetically tractable model system for the study of reverse transcriptase fidelity.

Mutation in bone morphogenetic protein receptor-IB is associated with increased ovulation rate in Booroola Mérino ewes

Ewes from the Booroola strain of Australian Mérino sheep are characterized by high ovulation rate and litter size. This phenotype is due to the action of the FecBB allele of a major gene named FecB, as determined by statistical analysis of phenotypic data. By genetic analysis of 31 informative half-sib families from heterozygous sires, we showed that the FecB locus is situated in the region of ovine chromosome 6 corresponding to the human chromosome 4q22–23 that contains the bone morphogenetic protein receptor IB (BMPR-IB) gene encoding a member of the transforming growth factor-β (TGF-β) receptor family. A nonconservative substitution (Q249R) in the BMPR-IB coding sequence was found to be associated fully with the hyperprolificacy phenotype of Booroola ewes. In vitro, ovarian granulosa cells from FecBB/FecBB ewes were less responsive than granulosa cells from FecB+/FecB+ ewes to the inhibitory effect on steroidogenesis of GDF-5 and BMP-4, natural ligands of BMPR-IB. It is suggested that in FecBB/FecBB ewes, BMPR-IB would be inactivated partially, leading to an advanced differentiation of granulosa cells and an advanced maturation of ovulatory follicles.

Rates of nucleotide substitution in sexual and anciently asexual rotifers

The class Bdelloidea of the phylum Rotifera is the largest well studied eukaryotic taxon in which males and meiosis are unknown, and the only one for which these indications of ancient asexuality are supported by cytological and molecular genetic evidence. We estimated the rates of synonymous and nonsynonymous substitutions in the hsp82 heat shock gene in bdelloids and in facultatively sexual rotifers of the class Monogononta, employing distance based and maximum likelihood methods. Relative-rate tests, using acanthocephalan rotifers as an outgroup, showed slightly higher rates of nonsynonymous substitution and slightly lower rates of synonymous substitution in bdelloids as compared with monogononts. The opposite trend, however, was seen in intraclass pairwise comparisons. If, as it seems, bdelloids have evolved asexually, an equality of bdelloid and monogonont substitution rates would suggest that the maintenance of sexual reproduction in monogononts is not attributable to an effect of sexual reproduction in limiting the load of deleterious nucleotide substitutions.

Estimation of evolutionary distances under stationary and nonstationary models of nucleotide substitution

Estimation of evolutionary distances has always been a major issue in the study of molecular evolution because evolutionary distances are required for estimating the rate of evolution in a gene, the divergence dates between genes or organisms, and the relationships among genes or organisms. Other closely related issues are the estimation of the pattern of nucleotide substitution, the estimation of the degree of rate variation among sites in a DNA sequence, and statistical testing of the molecular clock hypothesis. Mathematical treatments of these problems are considerably simplified by the assumption of a stationary process in which the nucleotide compositions of the sequences under study have remained approximately constant over time, and there now exist fairly extensive studies of stationary models of nucleotide substitution, although some problems remain to be solved. Nonstationary models are much more complex, but significant progress has been recently made by the development of the paralinear and LogDet distances. This paper reviews recent studies on the above issues and reports results on correcting the estimation bias of evolutionary distances, the estimation of the pattern of nucleotide substitution, and the estimation of rate variation among the sites in a sequence.

Novel human DNA alkyltransferases obtained by random substitution and genetic selection in bacteria.

DNA repair alkyltransferases protect organisms against the cytotoxic, mutagenic, and carcinogenic effects of alkylating agents by transferring alkyl adducts from DNA to an active cysteine on the protein, thereby restoring the native DNA structure. We used random sequence substitutions to gain structure-function information about the human O6-methylguanine-DNA methyltransferase (EC 2.1.1.63), as well as to create active mutants. Twelve codons surrounding but not including the active cysteine were replaced by a random nucleotide sequence, and the resulting random library was selected for the ability to provide alkyltransferase-deficient Escherichia coli with resistance to the methylating agent N-methyl-N'-nitro-N-nitrosoguanidine. Few amino acid changes were tolerated in this evolutionarily conserved region of the protein. One mutation, a valine to phenylalanine change at codon 139 (V139F), was found in 70% of the selected mutants; in fact, this mutant was selected much more frequently than the wild type. V139F provided alkyltransferase-deficient bacteria with greater protection than the wild-type protein against both the cytotoxic and mutagenic effects of N-methyl-N'-nitro-N-nitrosoguanidine, increasing the D37 over 4-fold and reducing the mutagenesis rate 2.7-5.5-fold. This mutant human alkyltransferase, or others similarly created and selected, could be used to protect bone marrow cells from the cytotoxic side effects of alkylation-based chemotherapeutic regimens.

A general additive distance with time-reversibility and rate variation among nucleotide sites.

As additivity is a very useful property for a distance measure, a general additive distance is proposed under the stationary time-reversible (SR) model of nucleotide substitution or, more generally, under the stationary, time-reversible, and rate variable (SRV) model, which allows rate variation among nucleotide sites. A method for estimating the mean distance and the sampling variance is developed. In addition, a method is developed for estimating the variance-covariance matrix of distances, which is useful for the statistical test of phylogenies and molecular clocks. Computer simulation shows (i) if the sequences are longer than, say, 1000 bp, the SR method is preferable to simpler methods; (ii) the SR method is robust against deviations from time-reversibility; (iii) when the rate varies among sites, the SRV method is much better than the SR method because the distance is seriously underestimated by the SR method; and (iv) our method for estimating the sampling variance is accurate for sequences longer than 500 bp. Finally, a test is constructed for testing whether DNA evolution follows a general Markovian model.

Internal ribosomal entry site substitution eliminates neurovirulence in intergeneric poliovirus recombinants.

Neuropathogenicity of poliovirus can be attenuated by mutations in the internal ribosomal entry site (IRES) within the 5' nontranslated region of its genome. The Sabin vaccine strains used in prevention of poliomyelitis carry such mutations in their IRES elements. In addition, mutations within the structural and nonstructural proteins of Sabin strains may equally contribute to the attenuation phenotype. Despite their effectiveness as vaccines, the Sabin strains retain a neuropathogenic potential in animal models for poliomyelitis and, at a very low rate, they can cause poliomyelitis in vaccine recipients. The elimination of the neurocytopathic phenotype was achieved through the exchange of the entire poliovirus IRES with its counterpart from human rhinovirus type 2 without affecting growth properties in nonneuronal cells. The attenuating effect of the human rhinovirus type 2 IRES within the context of a poliovirus genome has been mapped to the 3' portion of this genetic element.