Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects

A number of studies indicated that lineages of animals with high rates of mitochondrial (mt) gene rearrangement might have high rates of mt nucleotide substitution. We chose the hemipteroid assemblage and the Insecta to test the idea that rates of mt gene rearrangement and mt nucleotide substitution are correlated. For this purpose, we sequenced the mt genome of a lepidopsocid from the Psocoptera, the only order of hemipteroid insects for which an entire mtDNA sequence is not available. The mt genome of this lepidopsocid is circular, 16,924 bp long, and contains 37 genes and a putative control region; seven tRNA genes and a protein-coding gene in this genome have changed positions relative to the ancestral arrangement of mt genes of insects. We then compared the relative rates of nucleotide substitution among species from each of the four orders of hemipteroid insects and among the 20 insects whose mt genomes have been sequenced entirely. All comparisons among the hernipteroid insects showed that species with higher rates of gene rearrangement also had significantly higher rates of nucleotide substitution statistically than did species with lower rates of gene rearrangement. In comparisons among the 20 insects, where the mt genomes of the two species differed by more than five breakpoints, the more rearranged species always had a significantly higher rate of nucleotide substitution than the less rearranged species. However, in comparisons where the mt genomes of two species differed by five or less breakpoints, the more rearranged species did not always have a significantly higher rate of nucleotide substitution than the less rearranged species. We tested the statistical significance of the correlation between the rates of mt gene rearrangement and mt nucleotide substitution with nine pairs of insects that were phylogenetically independent from one 2 another. We found that the correlation was positive and statistically significant (R-2 = 0.73, P = 0.01; R-s = 0.67, P < 0.05). We propose that increased rates of nucleotide substitution may lead to increased rates of gene rearrangement in the mt genomes of insects.

Patterns of nucleotide substitution in angiosperm cpDNA trnL (UAA)-trnF (GAA) regions

Patterns of substitution in chloroplast encoded trnL_F regions were compared between species of Actaea (Ranunculales), Digitalis (Scrophulariales), Drosera (Caryophyllales), Panicoideae (Poales), the small chromosome species clade of Pelargonium (Geraniales), each representing a different order of flowering plants, and Huperzia (Lycopodiales). In total, the study included 265 taxa, each with > 900-bp sequences, totaling 0.24 Mb. Both pairwise and phylogeny-based comparisons were used to assess nucleotide substitution patterns. In all six groups, we found that transition/transversion ratios, as estimated by maximum likelihood on most-parsimonious trees, ranged between 0.8 and 1.0 for ingroups. These values occurred both at low sequence divergences, where substitutional saturation, i.e., multiple substitutions having occurred at the same (homologous) nucleotide position, was not expected, and at higher levels of divergence. This suggests that the angiosperm trnL-F regions evolve in a pattern different from that generally observed for nuclear and animal mtDNA (transitional/transversion ratio > or = 2). Transition/transversion ratios in the intron and the spacer region differed in all alignments compared, yet base compositions between the regions were highly similar in all six groups. A>-C transversions were significantly less frequent than the other four substitution types. This correlates with results from studies on fidelity mechanisms in DNA replication that predict A<->T and G<->C transversions to be least likely to occur. It therefore strengthens confidence in the link between mutation bias at the polymerase level and the actual fixation of substitutions as recorded on evolutionary trees, and concomitantly, in the neutrality of nucleotide substitutions as phylogenetic markers.

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.

Relationships among the three major lineages of the Acari (Arthropoda : Arachnida) inferred from small subunit rRNA: Paraphyly of the parasitiformes with respect to the opilioacariformes and relative rates of nucleotide substitution

We inferred phylogeny among the three major lineages of the Acari ( mites) from the small subunit rRNA gene. Our phylogeny indicates that the Opilioacariformes is the sister-group to the Ixodida+Holothyrida, not the Ixodida+Mesostigmata+Holothyrida, as previously thought. Support for this relationship increased when sites with the highest rates of nucleotide substitution, and thus the greatest potential for saturation with nucleotide substitutions, were removed. Indeed, the increase in support ( and resolution) was despite a 70% reduction in the number of parsimony-informative sites from 408 to 115. This shows that rather than 'noisy' sites having no impact on resolution of deep branches, 'noisy' sites have the potential to obscure phylogenetic relationships. The arrangement, Ixodida+Holothyrida+Opilioacariformes, however, may be an artefact of long-branch attraction since relative-rate tests showed that the Mesostigmata have significantly faster rates of nucleotide substitution than other parasitiform mites. Thus, the fast rates of nucleotide substitution of the Mesostigmata might have caused the Mesostigmata to be attracted to the outgroup in our trees. We tested the hypothesis that the high rate of nucleotide substitution in some mites was related to their short generation times. The Acari species that have high nucleotide substitution rates usually have short generation times; these mites also tend to be more active and thus have higher metabolic rates than other mites. Therefore, more than one factor may affect the rate of nucleotide substitution in these mites.

STUDIES ON THE PATTERNS OF NUCLEOTIDE AND AMINO ACID SUBSTITUTION (ELECTROPHORESIS, MUTATION RATE, SELECTION, BASE, COMPOSITION)

Theoretical and empirical studies were conducted on the pattern of nucleotide and amino acid substitution in evolution, taking into account the effects of mutation at the nucleotide level and purifying selection at the amino acid level. A theoretical model for predicting the evolutionary change in electrophoretic mobility of a protein was also developed by using information on the pattern of amino acid substitution. The specific problems studied and the main results obtained are as follows: (1) Estimation of the pattern of nucleotide substitution in DNA nuclear genomes. The pattern of point mutations and nucleotide substitutions among the four different nucleotides are inferred from the evolutionary changes of pseudogenes and functional genes, respectively. Both patterns are non-random, the rate of change varying considerably with nucleotide pair, and that in both cases transitions occur somewhat more frequently than transversions. In protein evolution, substitution occurs more often between amino acids with similar physico-chemical properties than between dissimilar amino acids. (2) Estimation of the pattern of nucleotide substitution in RNA genomes. The majority of mutations in retroviruses accumulate at the reverse transcription stage. Selection at the amino acid level is very weak, and almost non-existent between synonymous codons. The pattern of mutation is very different from that in DNA genomes. Nevertheless, the pattern of purifying selection at the amino acid level is similar to that in DNA genomes, although selection intensity is much weaker. (3) Evaluation of the determinants of molecular evolutionary rates in protein-coding genes. Based on rates of nucleotide substitution for mammalian genes, the rate of amino acid substitution of a protein is determined by its amino acid composition. The content of glycine is shown to correlate strongly and negatively with the rate of substitution. Empirical formulae, called indices of mutability, are developed in order to predict the rate of molecular evolution of a protein from data on its amino acid sequence. (4) Studies on the evolutionary patterns of electrophoretic mobility of proteins. A theoretical model was constructed that predicts the electric charge of a protein at any given pH and its isoelectric point from data on its primary and quaternary structures. Using this model, the evolutionary change in electrophoretic mobilities of different proteins and the expected amount of electrophoretically hidden genetic variation were studied. In the absence of selection for the pI value, proteins will on the average evolve toward a mildly basic pI. (Abstract shortened with permission of author.) ^

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.

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.

New Genotype of Dengue Type 3 Virus Circulating in Brazil and Colombia Showed a Close Relationship to Old Asian Viruses

Dengue type 3 genotype V viruses have been recently detected in Brazil and Colombia. In this study, we described another Brazilian isolate belonging to this genotype. Phylogenetic analysis including dengue type 3 viruses isolated worldwide showed that Brazilian and Colombian viruses were closely related to viruses isolated in Asia more than two decades ago. The characteristic evolutionary pattern of dengue type 3 virus cannot explain the close similarity of new circulating viruses with old viruses. Further studies are needed to confirm the origin of the new dengue type III genotype circulating in Brazil and Colombia.

Sensitivity to cisplatin-induced mutations and elevated chromosomal aberrations in lymphocytes from sickle cell disease patients

Sickle cell disease (SCD) is an inherited disorder caused by a single nucleotide substitution in the P-globin gene. The clinical heterogeneity observed in SCD patients has been attributed to environmental and genetic factors. The patients are subjected to increased oxidative stress, particularly during vaso-occlusive crises and acute chest pain. Another possible cause of oxidative stress in SCD is the high concentration of iron in the patients` plasma. The increase in oxidative stress could be a relevant risk factor for mutagenesis and carcinogenesis. Studies on the frequency of basal chromosomal aberrations in cultured lymphocytes from SCD patients have not been reported so far. In order to contribute to the understanding of the role of the different biomarkers and their relationship with the extremely variable clinical manifestation of SCD, we investigated the frequency of chromosome damage in peripheral lymphocytes from sickle cells patients and healthy controls. We found an increased frequency of chromosome damage and percentage of aberrant metaphases in these patients when compared with control subjects, even at basal values (p < 0.05). In the cytogenetic sensitivity assay, the results showed that these patients presented a marked decrease in the mitotic index values compared with healthy controls. Cisplatin-induced chromosomal damage in lymphocytes from these patients was significantly higher than the frequency measured in healthy controls. The results obtained in the present study showed that more investigations are needed in order to elucidate the susceptibility to genomic instability of SCD patients.

Bacterial expression of two human aryl sulfotransferases

The effect of replacing a single codon in the N-terminal of human aryl sulfotransferase (HAST) 1 and 3 with one that is more commonly found in E. coli genes was assessed. The pKK233-2 E. coli expression vector was employed and the polymerase chain reaction (PCR) was used to introduce the 5' nucleotide substitution, at the same time maintaining the fidelity of the amino acid sequence. The data indicates that this change had a minimal effect on protein production, subcellular localization or, in the case of HAST3, catalytic activity. In general, the pKK233-2 E. coli vector has been less than optimal for expressing human sulfotransferase cDNAs. (C) 1998 Elsevier Science Ireland Ltd. All rights reserved.

Dynamics of Hepatitis D (delta) virus genotype 3 in the Amazon region of South America

Hepatitis delta virus (HDV) is widely distributed and associated with fulminant hepatitis epidemics in areas with high prevalence of HBV. Several studies performed in the 1980s showed data on HDV infection in South America, but there are no studies on the viral dynamics of this virus. The aim of this study was to conduct an evolutionary analysis of hepatitis delta genotype 3 (HDV/3) prevalent in South America: estimate its nucleotide substitution rate, determine the time of most recent ancestor (TMRCA) and characterize the epidemic history and evolutionary dynamics. Furthermore, we characterized the presence of HBV/HDV infection in seven samples collected from patients who died due to fulminant hepatitis from Amazon region in Colombia and included them in the evolutionary analysis. This is the first study reporting HBV and HDV sequences from the Amazon region of Colombia. Of the seven Colombian patients, five were positive for HBV-DNA and HDV-RNA. Of them, two samples were successfully sequenced for HBV (subgenotypes F3 and Fib) and the five samples HDV positive were classified as HDV/3. By using all HDV/3 available reference sequences with sampling dates (n = 36), we estimated the HDV/3 substitution rate in 1.07 x 10(-3) substitutions per site per year (s/s/y), which resulted in a time to the most recent common ancestor (TMRCA) of 85 years. Also, it was determined that HDV/3 spread exponentially from early 1950s to the 1970s in South America. This work discusses for the first time the viral dynamics for the HDV/3 circulating in South America. We suggest that the measures implemented to control HBV transmission resulted in the control of HDV/3 spreading in South America, especially after the important raise in this infection associated with a huge mortality during the 1950s up to the 1970s. The differences found among HDV/3 and the other HDV genotypes concerning its diversity raises the hypothesis of a different origin and/or a different transmission route. (C) 2011 Elsevier B.V. All rights reserved.

THE TEMPO AND MODE OF EVOLUTION OF TRANSPOSABLE ELEMENTS AS REVEALED BY MOLECULAR PHYLOGENIES RECONSTRUCTED FROM MOSQUITO GENOMES

Although many mathematical models exist predicting the dynamics of transposable elements (TEs), there is a lack of available empirical data to validate these models and inherent assumptions. Genomes can provide a snapshot of several TE families in a single organism, and these could have their demographics inferred by coalescent analysis, allowing for the testing of theories on TE amplification dynamics. Using the available genomes of the mosquitoes Aedes aegypti and Anopheles gambiae, we indicate that such an approach is feasible. Our analysis follows four steps: (1) mining the two mosquito genomes currently available in search of TE families; (2) fitting, to selected families found in (1), a phylogeny tree under the general time-reversible (GTR) nucleotide substitution model with an uncorrelated lognormal (UCLN) relaxed clock and a nonparametric demographic model; (3) fitting a nonparametric coalescent model to the tree generated in (2); and (4) fitting parametric models motivated by ecological theories to the curve generated in (3).

High rates of molecular evolution in hantaviruses

Hantaviruses are rodent-borne Bunyaviruses that infect the Arvicolinae, Murinae, and Sigmodontinae subfamilies of Muridae. The rate of molecular evolution in the hantaviruses has been previously estimated at approximately 10(-7) nucleotide substitutions per site, per year (substitutions/site/year), based on the assumption of codivergence and hence shared divergence times with their rodent hosts. If substantiated, this would make the hantaviruses among the slowest evolving of all RNA viruses. However, as hantaviruses replicate with an RNA-dependent RNA polymerase, with error rates in the region of one mutation per genome replication, this low rate of nucleotide substitution is anomalous. Here, we use a Bayesian coalescent approach to estimate the rate of nucleotide substitution from serially sampled gene sequence data for hantaviruses known to infect each of the 3 rodent subfamilies: Araraquara virus ( Sigmodontinae), Dobrava virus ( Murinae), Puumala virus ( Arvicolinae), and Tula virus ( Arvicolinae). Our results reveal that hantaviruses exhibit shortterm substitution rates of 10(-2) to 10(-4) substitutions/site/year and so are within the range exhibited by other RNA viruses. The disparity between this substitution rate and that estimated assuming rodent-hantavirus codivergence suggests that the codivergence hypothesis may need to be reevaluated.

No significant association between progesterone receptor exon 4 Val660Leu G/T polymorphism and risk of ovarian cancer

Epidemiological studies suggest that ovarian cancer is an endocrine-related tumour, and progesterone exposure specifically may decrease the risk of ovarian cancer. To assess whether the progesterone receptor (PR) exon 4 valine to leucine amino acid variant is associated with specific tumour characteristics or with overall risk of ovarian cancer, we examined 551 cases of epithelial ovarian cancer and 298 unaffected controls for the underlying G-->T nucleotide substitution polymorphism. Stratification of the ovarian cancer cases according to tumour behaviour (low malignant potential or invasive), histology, grade or stage failed to reveal any heterogeneity with respect to the genotype defined by the PR exon 4 polymorphism. Furthermore, the genotype distribution did not differ significantly between ovarian cancer cases and unaffected controls. Compared with the GG genotype, the age-adjusted odds ratio (95% confidence interval) for risk of ovarian cancer was 0.78 (0.57-1.08) for the GT genotype, and 1.39 (0.47-4.14) for the TT genotype. In conclusion, the PR exon 4 codon 660 leucine variant encoded by the T allele does not appear to be associated with ovarian tumour behaviour, histology, stage or grade. This variant is also not associated with an increased risk of ovarian cancer, and is unlikely to be associated with a large decrease in ovarian cancer risk, although we cannot rule out a moderate inverse association between the GT genotype and ovarian cancer.