Pattern of synonymous and nonsynonymous substitutions: An indicator of mechanisms of molecular evolution

Pheromone Response and Mating Behavior in Fission Yeast

Most ascomycete fungi, including the fission yeast Schizosaccharomyces pombe, secrete two peptidyl mating pheromones: C-terminally modified and unmodified peptides. S. pombe has two mating types, plus and minus, which secrete two different pheromones, P-factor (unmodified) and M-factor (modified), respectively. These pheromones are specifically recognized by receptors on the cell surface of cells of opposite mating types, which trigger a pheromone response. Recognition between pheromones and their corresponding receptors is important for mate discrimination; therefore, genetic changes in pheromone or receptor genes affect mate recognition and cause reproductive isolation that limits gene flow between populations. Such genetic variation in recognition via the pheromone/receptor system may drive speciation. Our recent studies reported that two pheromone receptors in S. pombe might have different stringencies in pheromone recognition. In this review, we focus on the molecular mechanism of pheromone response and mating behavior, emphasizing pheromone diversification and its impact on reproductive isolation in S. pombe and closely related fission yeast species. We speculate that the "asymmetric" system might allow flexible adaptation to pheromone mutational changes while maintaining stringent recognition of mating partners. The loss of pheromone activity results in the extinction of an organism's lineage. Therefore, genetic changes in pheromones and their receptors may occur gradually and/or coincidently before speciation. Our findings suggest that the M-factor plays an important role in partner discrimination, whereas P-factor communication allows flexible adaptation to create variations in S. pombe. Our inferences provide new insights into the evolutionary mechanisms underlying pheromone diversification.

Mitogenome selection in the evolution of key ecological strategies in the ancient hexapod class Collembola

A longstanding question in evolutionary biology is how natural selection and environmental pressures shape the mitochondrial genomic architectures of organisms. Mitochondria play a pivotal role in cellular respiration and aerobic metabolism, making their genomes functionally highly constrained. Evaluating selective pressures on mitochondrial genes can provide functional and ecological insights into the evolution of organisms. Collembola (springtails) are an ancient hexapod group that includes the oldest terrestrial arthropods in the fossil record, and that are closely associated with soil environments. Of interest is the diversity of habitat stratification preferences (life forms) exhibited by different species within the group. To understand whether signals of positive selection are linked to the evolution of life forms, we analysed 32 published Collembola mitogenomes in a phylomitogenomic framework. We found no evidence that signatures of selection are correlated with the evolution of novel life forms, but rather that mutations have accumulated as a function of time. Our results highlight the importance of nuclear-mitochondrial interactions in the evolution of collembolan life forms and that mitochondrial genomic data should be interpreted with caution, as complex selection signals may complicate evolutionary inferences.

Combined Methylome, Transcriptome and Proteome Analyses Document Rapid Acclimatization of a Bacterium to Environmental Changes

Polynucleobacter asymbioticus strain QLW-P1DMWA-1^T represents a group of highly successful heterotrophic ultramicrobacteria that is frequently very abundant (up to 70% of total bacterioplankton) in freshwater habitats across all seven continents. This strain was originally isolated from a shallow Alpine pond characterized by rapid changes in water temperature and elevated UV radiation due to its location at an altitude of 1300 m. To elucidate the strain’s adjustment to fluctuating environmental conditions, we recorded changes occurring in its transcriptomic and proteomic profiles under contrasting experimental conditions by simulating thermal conditions in winter and summer as well as high UV irradiation. To analyze the potential connection between gene expression and regulation via methyl group modification of the genome, we also analyzed its methylome. The methylation pattern differed between the three treatments, pointing to its potential role in differential gene expression. An adaptive process due to evolutionary pressure in the genus was deduced by calculating the ratios of non-synonymous to synonymous substitution rates for 20 Polynucleobacter spp. genomes obtained from geographically diverse isolates. The results indicate purifying selection.

Comprehensive study of the genes involved in chlorophyll synthesis and degradation pathways in some monocot and dicot plant species

Chlorophyll (Chl) biosynthesis is one of the most important cellular processes essential for plant photosynthesis. Chl degradation pathway is also important catabolic process occurs during leaf senescence, fruit ripening and under biotic or abiotic stress conditions. Here we have systematically investigated the molecular evolution, gene structure, compositional analysis along with ENc plot, correspondence analysis and codon usage bias of the proteins and encoded genes involved in Chl metabolism from monocots and dicots. The gene and species specific phylogenetic trees using amino acid sequences showed clear clustering formation of the selected species based on monocots and dicots but not supported by 18S rRNA. Nucleotide composition of the encoding genes showed that average GC%, GC1%, GC2% and GC3% were higher in monocots. RSCU analysis depicts that genes from monocots for both pathways and genes for synthesis pathway from dicots only biased to G/C-ending synonymous codons but in degradation pathway most optimal codons (except UUG) in dicots biased to A/U-ending synonymous codons. We found strong evidence of episodic diversifying selection at several amino acid sites in all genes investigated. Conserved domain and gene structures were observed for the genes with varying lengths of introns and exons, involved in Chl metabolism along with some intronless genes within synthesis pathway. ENc and correspondence analyses suggested the mutational or selection constraint on the genes to shape the codon usage. These comprehensive studies may be helpful in further research in molecular phylogenetics and genomics and to better understand the evolutionary dynamics of Chl metabolic pathway.Communicated by Ramaswamy H. Sarma.

Genome Sequencing and Analysis of the Peanut B-Genome Progenitor (Arachis ipaensis)

Peanut (Arachis hypogaea L.), an important leguminous crop, is widely cultivated in tropical and subtropical regions. Peanut is an allotetraploid, having A and B subgenomes that maybe have originated in its diploid progenitors Arachis duranensis (A-genome) and Arachis ipaensis (B-genome), respectively. We previously sequenced the former and here present the draft genome of the latter, expanding our knowledge of the unique biology of Arachis. The assembled genome of A. ipaensis is ~1.39 Gb with 39,704 predicted protein-encoding genes. A gene family analysis revealed that the FAR1 family may be involved in regulating peanut special fruit development. Genomic evolutionary analyses estimated that the two progenitors diverged ~3.3 million years ago and suggested that A. ipaensis experienced a whole-genome duplication event after the divergence of Glycine max. We identified a set of disease resistance-related genes and candidate genes for biological nitrogen fixation. In particular, two and four homologous genes that may be involved in the regulation of nodule development were obtained from A. ipaensis and A. duranensis, respectively. We outline a comprehensive network involved in drought adaptation. Additionally, we analyzed the metabolic pathways involved in oil biosynthesis and found genes related to fatty acid and triacylglycerol synthesis. Importantly, three new FAD2 homologous genes were identified from A. ipaensis and one was completely homologous at the amino acid level with FAD2 from A. hypogaea. The availability of the A. ipaensis and A. duranensis genomic assemblies will advance our knowledge of the peanut genome.

New insights into the evolutionary history of plant sorbitol dehydrogenase

BACKGROUND

Sorbitol dehydrogenase (SDH, EC 1.1.1.14) is the key enzyme involved in sorbitol metabolism in higher plants. SDH genes in some Rosaceae species could be divided into two groups. L-idonate-5-dehydrogenase (LIDH, EC 1.1.1.264) is involved in tartaric acid (TA) synthesis in Vitis vinifera and is highly homologous to plant SDHs. Despite efforts to understand the biological functions of plant SDH, the evolutionary history of plant SDH genes and their phylogenetic relationship with the V. vinifera LIDH gene have not been characterized.

RESULTS

A total of 92 SDH genes were identified from 42 angiosperm species. SDH genes have been highly duplicated within the Rosaceae family while monocot, Brassicaceae and most Asterid species exhibit singleton SDH genes. Core Eudicot SDHs have diverged into two phylogenetic lineages, now classified as SDH Class I and SDH Class II. V. vinifera LIDH was identified as a Class II SDH. Tandem duplication played a dominant role in the expansion of plant SDH family and Class II SDH genes were positioned in tandem with Class I SDH genes in several plant genomes. Protein modelling analyses of V. vinifera SDHs revealed 19 putative active site residues, three of which exhibited amino acid substitutions between Class I and Class II SDHs and were influenced by positive natural selection in the SDH Class II lineage. Gene expression analyses also demonstrated a clear transcriptional divergence between Class I and Class II SDH genes in V. vinifera and Citrus sinensis (orange).

CONCLUSIONS

Phylogenetic, natural selection and synteny analyses provided strong support for the emergence of SDH Class II by positive natural selection after tandem duplication in the common ancestor of core Eudicot plants. The substitutions of three putative active site residues might be responsible for the unique enzyme activity of V. vinifera LIDH, which belongs to SDH Class II and represents a novel function of SDH in V. vinifera that may be true also of other Class II SDHs. Gene expression analyses also supported the divergence of SDH Class II at the expression level. This study will facilitate future research into understanding the biological functions of plant SDHs.

Causes and consequences of variability in peptide mating pheromones of ascomycete fungi

The reproductive genes of fungi, like those of many other organisms, are thought to diversify rapidly. This phenomenon could be associated with the formation of reproductive barriers and speciation. Ascomycetes produce two classes of mating type-specific peptide pheromones. These are required for recognition between the mating types of heterothallic species. Little is known regarding the diversity or the extent of species specificity in pheromone peptides among these fungi. We compared the putative protein-coding DNA sequences of the 2 pheromone classes from 70 species of Ascomycetes. The data set included previously described pheromones and putative pheromones identified from genomic sequences. In addition, pheromone genes from 12 Fusarium species in the Gibberella fujikuroi complex were amplified and sequenced. Pheromones were largely conserved among species in this complex and, therefore, cannot alone account for the reproductive barriers observed between these species. In contrast, pheromone peptides were highly diverse among many other Ascomycetes, with evidence for both positive diversifying selection and relaxed selective constraint. Repeats of the α-factor-like pheromone, which occur in tandem arrays of variable copy number, were found to be conserved through purifying selection and not concerted evolution. This implies that sequence specificity may be important for pheromone reception and that interspecific differences may indeed be associated with functional divergence. Our findings also suggest that frequent duplication and loss causes the tandem repeats to experience "birth-and-death" evolution, which could in fact facilitate interspecific divergence of pheromone peptide sequences.

Phylogeny of the nematode genus Pristionchus and implications for biodiversity, biogeography and the evolution of hermaphroditism

BACKGROUND

The nematode Pristionchus pacificus has originally been developed as a satellite organism for comparison to Caenorhabditis elegans. A 10X coverage of the whole genome of P. pacificus is available, making P. pacificus the first non-Caenorhabditis nematode with a fully sequenced genome. The macroevolutionary comparison between P. pacificus and C. elegans has been complemented by microevolutionary studies of closely related strains and species within the genus Pristionchus. In addition, new understanding of the biology of Pristionchus from field studies, demonstrating a close association with various scarab beetles and the Colorado potato beetle, supports consideration of this nematode in studies of ecosystems. In the course of field studies on four continents more than 1,200 isolates were established from 15,000 beetle specimens representing 18 Pristionchus species. Two remarkable features of the Pristionchus-beetle association are the high species specificity of the interaction and the interception of the beetle's sex communication system for host recognition by the nematodes, as suggested by chemotaxis studies. Evolutionary interpretations of differences in developmental, behavioral and ecological patterns require a phylogenetic framework of the genus Pristionchus.

RESULTS

Here, we provide a robust phylogeny of all 18 available Pristionchus species based on a set of 27 ribosomal protein genes encompassing a total of 10,971 bp. The phylogenetic tree provides evidence for North American and European clades, which are embedded in a deeper clade that includes Asian species. It also indicates putative invasion events. Of the 18 Pristionchus species, 13 are gonochoristic and five are hermaphroditic. The phylogeny indicates that all hermaphroditic species have arisen independently within the genus Pristionchus.

CONCLUSION

Combined ribosomal protein cDNA data can provide the basis for reconstruction of a robust phylogenetic framework for microevolutionary and biogeographic analyses. An additional major implication of our studies is the use of Pristionchus for nematode biodiversity assessments. While some species are represented by more than 100 isolates, others were found less than four times. Such patterns were observed on all continents and in all phylogenetic clades indicating that species asymmetry is a widespread phenomenon, which can now be further investigated by molecular tools.

Differentiation of the two rice subspecies indica and japonica: a Gene Ontology perspective

Cultivated rice (Oryza sativa) is comprised of two subspecies: japonica and indica. Polymorphism levels between putative homologues were determined for genes whose japonica homologue had been classified into functional categories using the Gene Ontology (GO) system. Genes were partitioned into below-average and above-average polymorphism groups, and then the set of genes having each GO term was checked for the randomness of its distribution into these polymorphism groups using a series of False Discovery Rate (FDR) tests. The robustness of the conclusions was enhanced by employing different cutoff values and sequence samplings in the FDR tests. Significant nonrandom polymorphism distributions were found for protein-coding sequences in many GO categories. In contrast, a random distribution for nearly all GO terms was seen with intron sequences. These results were extended by measuring the nonsynonymous to synonymous codon usage ratio (dN/dS) using a permutation test, which showed that some above-average polymorphism GO categories also had a high proportion of genes with a dN/dS ratio greater than one, suggesting positive selection on these GO categories during indica-japonica differentiation. An analysis of predominant gene names in the significant GO categories divided them into four functional classes: production of defense-related compounds, cell wall, cell signaling, and transcription factors.

Nucleotide and predicted amino acid sequence-based analysis of the avian metapneumovirus type C cell attachment glycoprotein gene: phylogenetic analysis and molecular epidemiology of U.S. pneumoviruses

A serologically distinct avian metapneumovirus (aMPV) was isolated in the United States after an outbreak of turkey rhinotracheitis (TRT) in February 1997. The newly recognized U.S. virus was subsequently demonstrated to be genetically distinct from European subtypes and was designated aMPV serotype C (aMPV/C). We have determined the nucleotide sequence of the gene encoding the cell attachment glycoprotein (G) of aMPV/C (Colorado strain and three Minnesota isolates) and predicted amino acid sequence by sequencing cloned cDNAs synthesized from intracellular RNA of aMPV/C-infected cells. The nucleotide sequence comprised 1,321 nucleotides with only one predicted open reading frame encoding a protein of 435 amino acids, with a predicted M(r) of 48,840. The structural characteristics of the predicted G protein of aMPV/C were similar to those of the human respiratory syncytial virus (hRSV) attachment G protein, including two mucin-like regions (heparin-binding domains) flanking both sides of a CX3C chemokine motif present in a conserved hydrophobic pocket. Comparison of the deduced G-protein amino acid sequence of aMPV/C with those of aMPV serotypes A, B, and D, as well as hRSV revealed overall predicted amino acid sequence identities ranging from 4 to 16.5%, suggesting a distant relationship. However, G-protein sequence identities ranged from 72 to 97% when aMPV/C was compared to other members within the aMPV/C subtype or 21% for the recently identified human MPV (hMPV) G protein. Ratios of nonsynonymous to synonymous nucleotide changes were greater than one in the G gene when comparing the more recent Minnesota isolates to the original Colorado isolate. Epidemiologically, this indicates positive selection among U.S. isolates since the first outbreak of TRT in the United States.

Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models

Approximate methods for estimating the numbers of synonymous and nonsynonymous substitutions between two DNA sequences involve three steps: counting of synonymous and nonsynonymous sites in the two sequences, counting of synonymous and nonsynonymous differences between the two sequences, and correcting for multiple substitutions at the same site. We examine complexities involved in those steps and propose a new approximate method that takes into account two major features of DNA sequence evolution: transition/transversion rate bias and base/codon frequency bias. We compare the new method with maximum likelihood, as well as several other approximate methods, by examining infinitely long sequences, performing computer simulations, and analyzing a real data set. The results suggest that when there are transition/transversion rate biases and base/codon frequency biases, previously described approximate methods for estimating the nonsynonymous/synonymous rate ratio may involve serious biases, and the bias can be both positive and negative. The new method is, in general, superior to earlier approximate methods and may be useful for analyzing large data sets, although maximum likelihood appears to always be the method of choice.

The effect of tandem substitutions on the correlation between synonymous and nonsynonymous rates in rodents

Nonsynonymous substitutions in DNA cause amino acid substitutions while synonymous substitutions in DNA leave amino acids unchanged. The cause of the correlation between the substitution rates at nonsynonymous (K(A)) and synonymous (K(S)) sites in mammals is a contentious issue, and one that impacts on many aspects of molecular evolution. Here we use a large set of orthologous mammalian genes to investigate the causes of the K(A)-K(S) correlation in rodents. The strength of the K(A)-K(S) correlation exceeds the neutral theory expectation when substitution rates are estimated using algorithmic methods, but not when substitution rates are estimated by maximum likelihood. Irrespective of this methodological uncertainty the strength of the K(A)-K(S) correlation appears mostly due to tandem substitutions, an excess of which is generated by substitutional nonindependence. Doublet mutations cannot explain the excess of tandem synonymous-nonsynonymous substitutions, and substitution patterns indicate that selection on silent sites is the likely cause. We find no evidence for selection on codon usage. The nature of the relationship between synonymous divergence and base composition is unclear because we find a significant correlation if we use maximum-likelihood methods but not if we use algorithmic methods. Finally, we find that K(S) is reduced at the start of genes, which suggests that selection for RNA structure may affect silent sites in mammalian protein-coding genes.

Inferring the fitness effects of DNA mutations from polymorphism and divergence data: statistical power to detect directional selection under stationarity and free recombination

The fitness effects of classes of DNA mutations can be inferred from patterns of nucleotide variation. A number of studies have attributed differences in levels of polymorphism and divergence between silent and replacement mutations to the action of natural selection. Here, I investigate the statistical power to detect directional selection through contrasts of DNA variation among functional categories of mutations. A variety of statistical approaches are applied to DNA data simulated under Sawyer and Hartl's Poisson random field model. Under assumptions of free recombination and stationarity, comparisons that include both the frequency distributions of mutations segregating within populations and the numbers of mutations fixed between populations have substantial power to detect even very weak selection. Frequency distribution and divergence tests are applied to silent and replacement mutations among five alleles of each of eight Drosophila simulans genes. Putatively "preferred" silent mutations segregate at higher frequencies and are more often fixed between species than "unpreferred" silent changes, suggesting fitness differences among synonymous codons. Amino acid changes tend to be either rare polymorphisms or fixed differences, consistent with a combination of deleterious and adaptive protein evolution. In these data, a substantial fraction of both silent and replacement DNA mutations appear to affect fitness.

Antarctic fish hemoglobins: evidence for adaptive evolution at subzero temperature

Notothenioids represent a large group of marine teleosts that are mostly endemic to the Antarctic Ocean. In this environment, the low metabolic demand and the high oxygen concentration reduce the need for hemoglobin(s) [Hb(s)]. The extreme condition is represented by the icefish (Channichthyidae, Notothenioidei), the only vertebrates that lack Hb. We obtained the nucleotide sequence coding for the beta-globin chain of the single major Hb form in six red-blooded notothenioids. These included Gymnodraco acuticeps, one of the closest species to the Hb-less icefish, which is also the only known fish having a single Hb without Bohr effect. This species shows a higher rate of nonsynonymous substitutions (KA), in contrast with the homogeneity of synonymous substitution (KS) rates, and KA/KS ratios significantly greater than one in the majority of comparisons. These results are suggestive of positive selection, diversifying the single major Hb toward specialized functions. A single Hb that is free to diversify means that its role in routine oxygen transport can be reduced in the presence of a combination of physiological, ecological, and environmental factors. Although a reduced "routine" function for Hb, as is apparent in G. acuticeps, might, indeed, evoke the lack of Hb in icefish, evidence of diversifying selection reported here is at variance with the hypothesis of a simple trend from a single Hb toward the Hb-less condition.

Synonymous and nonsynonymous rate variation in nuclear genes of mammals

A maximum likelihood approach was used to estimate the synonymous and nonsynonymous substitution rates in 48 nuclear genes from primates, artiodactyls, and rodents. A codon-substitution model was assumed, which accounts for the genetic code structure, transition/transversion bias, and base frequency biases at codon positions. Likelihood ratio tests were applied to test the constancy of nonsynonymous to synonymous rate ratios among branches (evolutionary lineages). It is found that at 22 of the 48 nuclear loci examined, the nonsynonymous/synonymous rate ratio varies significantly across branches of the tree. The result provides strong evidence against a strictly neutral model of molecular evolution. Our likelihood estimates of synonymous and nonsynonymous rates differ considerably from previous results obtained from approximate pairwise sequence comparisons. The differences between the methods are explored by detailed analyses of data from several genes. Transition/transversion rate bias and codon frequency biases are found to have significant effects on the estimation of synonymous and nonsynonymous rates, and approximate methods do not adequately account for those factors. The likelihood approach is preferable, even for pairwise sequence comparison, because more realistic models about the mutation and substitution processes can be incorporated in the analysis.

The genome sequence of herpes simplex virus type 2

The genomic DNA sequence of herpes simplex virus type 2 (HSV-2) strain HG52 was determined as 154,746 bp with a G+C content of 70.4%. A total of 74 genes encoding distinct proteins was identified; three of these were each present in two copies, within major repeat elements of the genome. The HSV-2 gene set corresponds closely with that of HSV-1, and the HSV-2 sequence prompted several local revisions to the published HSV-1 sequence (D. J. McGeoch, M. A. Dalrymple, A. J. Davison, A. Dolan, M. C. Frame, D. McNab, L. J. Perry, J. E. Scott, and P. Taylor, J. Gen. Virol. 69:1531-1574, 1988). No compelling evidence for the existence of any additional protein-coding genes in HSV-2 was identified.

Synonymous substitution rates in Drosophila: mitochondrial versus nuclear genes

Synonymous substitution rates in mitochondrial and nuclear genes of Drosophila were compared. To make accurate comparisons, we considered the following: (1) relative synonymous rates, which do not require divergence time estimates, should be used; (2) methods estimating divergence should take into account base composition; (3) only very closely related species should be used to avoid effects of saturation; (4) the heterogeneity of rates should be examined. We modified the methods estimating synonymous substitution numbers to account for base composition bias. By using these methods, we found that mitochondrial genes have 1.7-3.4 times higher synonymous substitution rates than the fastest nuclear genes or 4.5-9.0 times higher rates than the average nuclear genes. The average rate of synonymous transversions was 2.7 (estimated from the melanogaster species subgroup) or 2.9 (estimated from the obscura group) times higher in mitochondrial genes than in nuclear genes. Synonymous transversions in mitochondrial genes occurred at an approximately equivalent rate to those in the fastest nuclear genes. This last result is not consistent with the hypothesis that the difference in turnover rates between mitochondrial and nuclear genomes is the major factor determining higher synonymous substitution rates in mtDNA. We conclude that the difference in synonymous substitution rates is due to a combination of two factors: a higher transitional mutation rate in mtDNA and constraints on nuclear genes due to selection for codon usage.

Molecular evidence from the nuclear genome for the time frame of human evolution

Evolutionary divergence times can be inferred from molecular distances if a molecular clock can be assumed and if the substitution rate can be estimated. We present new evidence from relative rate tests that the rate of substitution at fourfold degenerate sites of nuclear genome-coding DNA is uniform in primate and rodent lineages. We also review recent relative rate test results showing substitution rate uniformity in the nuclear genome of simian primates. DNA distances between a range of mammalian taxa shows that a molecular clock is inconsistent with many assumed divergence times irrespective of the assumed substitution rate. We find that the substitution rate that implies the best compromise fit with divergence times across the range of taxa is 2.0-2.25 x 10(-9). This range of substitution rates implies a divergence time of humans and chimpanzees of 4.0-3.6 million years ago. This postdates the occurrence of Ardipithecus ramidus and the earliest occurrence of Australopithecus afarensis, suggesting that the common ancestor of humans and chimpanzees was bipedal and that the trait has been lost in chimpanzees rather than gained in humans.