Molecular population genetics of accessory gland protein genes and testis-expressed genes in Drosophila mojavensis and D. arizonae

Molecular population genetic investigation of Drosophila male reproductive genes has focused primarily on melanogaster subgroup accessory gland protein genes (Acp's). Consistent with observations from male reproductive genes of numerous taxa, Acp's evolve more rapidly than nonreproductive genes. However, within the Drosophila genus, large data sets from additional types of male reproductive genes and from different species groups are lacking. Here we report findings from a molecular population genetics analysis of male reproductive genes of the repleta group species, Drosophila arizonae and D. mojavensis. We find that Acp's have dramatically higher average pairwise Ka/Ks (0.93) than testis-enriched genes (0.19) and previously reported melanogaster subgroup Acp's (0.42). Overall, 10 of 19 Acp's have Ka/Ks > 1 either in nonpolarized analyses or in at least one lineage of polarized analyses. Of the nine Acp's for which outgroup data were available, average Ka/Ks was considerably higher in D. mojavensis (2.08) than in D. arizonae (0.87). Contrasts of polymorphism and divergence suggest that adaptive protein evolution at Acp's is more common in D. mojavensis than in D. arizonae.     

SRY基因在人猿超科和旧大陆猴中具有不同的进化规律

通过ＰＣＲ扩增、测序,得到了白臀叶猴和红面猴的ＳＲＹ基因全序列。结合现有的灵长类其他物种序列进行分析,验证了ＨＭＧ盒的保守性。通过构建系统发育树,比较旧大陆猴和人猿超科两个类群内和类群间ＨＭＧ盒侧翼序列Ｋａ／Ｋｓ的比率。有趣的是,人猿超科两物种比较呈现较高的Ｋａ／Ｋｓ比值,但在旧大陆猴中及旧大陆猴与狨猴间的Ｋａ／Ｋｓ比值显著低于人猿超科的,呈现很不同的格局。同时,对于ＨＭＧ盒序列,Ｋａ／Ｋｓ比值在     

A statistical analysis of nucleotide substitutions in the Drosophila Adh region reflects irregularities in molecular clocks

Substitutions rates are expected to be rather constant when a gene is compared between species. To analyze this feature, Ka/Ks ratios have been studied for Alcohol dehydrogenase (Adh) and Alcohol dehydrogenase duplication (Adh-dup) genes in Drosophila species. Adh Ka/Ks values are lower in intrasubgenus comparisons involving species of the Sophophora group than when these species are compared to the D. immigrans and S. lebanonensis, and this difference does not occur in the Adh-dup comparisons.     

Reductive genome evolution from the mother of Rickettsia

The Rickettsia genus is a group of obligate intracellular α-proteobacteria representing a paradigm of reductive evolution. Here, we investigate the evolutionary processes that shaped the genomes of the genus. The reconstruction of ancestral genomes indicates that their last common ancestor contained more genes, but already possessed most traits associated with cellular parasitism. The differences in gene repertoires across modern Rickettsia are mainly the result of differential gene losses from the ancestor. We demonstrate using computer simulation that the propensity of loss was variable across genes during this process. We also analyzed the ratio of nonsynonymous to synonymous changes (Ka/Ks) calculated as an average over large sets of genes to assay the strength of selection acting on the genomes of Rickettsia, Anaplasmataceae, and free-living γ-proteobacteria. As a general trend, Ka/Ks were found to decrease with increasing divergence between genomes. The high Ka/Ks for closely related genomes are probably due to a lag in the removal of slightly deleterious nonsynonymous mutations by natural selection. Interestingly, we also observed a decrease of the rate of gene loss with increasing divergence, suggesting a similar lag in the removal of slightly deleterious pseudogene alleles. For larger divergence (Ks > 0.2), Ka/Ks converge toward similar values indicating that the levels of selection are roughly equivalent between intracellular α-proteobacteria and their free-living relatives. This contrasts with the view that obligate endocellular microorganisms tend to evolve faster as a consequence of reduced effectiveness of selection, and suggests a major role of enhanced background mutation rates on the fast protein divergence in the obligate intracellular α-proteobacteria.     

Average allozyme heterozygosity in vertebrates correlates with Ka/Ks measured in the human-mouse lineage

It is well established that different allozyme proteins vary in heterozygosity in averages made over large numbers of species. For example, the enzyme 6-phosphogluconate dehydrogenase has a much higher average heterozygosity than glutamate dehydrogenase. Allozyme data alone provide insufficient power to determine the evolutionary cause of such a difference. Many studies have now been carried out on the DNA sequences coding for allozymes. These have identified diverse selective and nonselective causes of polymorphisms at individual loci. However the studies are mainly in a small number of model species; thus, it is difficult to identify from these DNA studies specific causes of global average heterozygosity differences among allozyme proteins. Here we demonstrate that estimates of average heterozygosity for 37 allozyme proteins in vertebrates correlate positively with Ka and Ka/Ks but not with Ks, measured in the human-mouse lineage. The values of Ka/Ks are less than 0.25, and Ka/Ks is negatively correlated with subunit number (quaternary structure), a measure of structural constraint. Proteins with lower levels of constraint have higher values of both Ka/Ks and heterozygosity. These results better support the hypothesis that differences in average allozyme diversity between proteins are more closely related to differences in the level of purifying selection than to differences in the underlying mutation rate or level of positive selection.     

Evaluation of whether accelerated protein evolution in chordates has occurred before, after, or simultaneously with gene duplication

Gene duplication and loss are predicted to be at least of the order of the substitution rate and are key contributors to the development of novel gene function and overall genome evolution. Although it has been established that proteins evolve more rapidly after gene duplication, we were interested in testing to what extent this reflects causation or association. Therefore, we investigated the rate of evolution prior to gene duplication in chordates. Two patterns emerged; firstly, branches, which are both preceded by a duplication and followed by a duplication, display an elevated rate of amino acid replacement. This is reflected in the ratio of nonsynonymous to synonymous substitution (mean nonsynonymous to synonymous nucleotide substitution rate ratio [Ka:Ks]) of 0.44 compared with branches preceded by and followed by a speciation (mean Ka:Ks of 0.23). The observed patterns suggest that there can be simultaneous alteration in the selection pressures on both gene duplication and amino acid replacement, which may be consistent with co-occurring increases in positive selection, or alternatively with concurrent relaxation of purifying selection. The pattern is largely, but perhaps not completely, explained by the existence of certain families that have elevated rates of both gene duplication and amino acid replacement. Secondly, we observed accelerated amino acid replacement prior to duplication (mean Ka:Ks for postspeciation preduplication branches was 0.27). In some cases, this could reflect adaptive changes in protein function precipitating a gene duplication event. In conclusion, the circumstances surrounding the birth of new proteins may frequently involve a simultaneous change in selection pressures on both gene-copy number and amino acid replacement. More precise modeling of the relative importance of preduplication, postduplication, and simultaneous amino acid replacement will require larger and denser genomic data sets from multiple species, allowing simultaneous estimation of lineage-specific fluctuations in mutation rates and adaptive constraints.     

Covariation of GC content and the silent site substitution rate in rodents: implications for methodology and for the evolution of isochores

Many attempts to test selectionist and neutralist models employ estimates of synonymous (Ks) and non-synonymous (Ka) substitution rates of orthologous genes. For example, a stronger Ka-Ks correlation than expected under neutrality has been argued to indicate a role for selection and the absence of a Ks-GC4 correlation has been argued to be inconsistent with neutral models for isochore evolution. However, both of these results, we have shown previously, are sensitive to the method by which Ka and Ks are estimated. Using a maximum likelihood (ML) estimator (GY94) we found a positive correlation between Ks and GC4 and only a weak correlation between Ka and Ks, lower than expected under neutral expectations. This ML method is computationally slow. Recently, a new ad hoc approximation of this ML method has been provided (YN00). This is effectively an extension of Li's protocol but that also allows for codon usage bias. This method is computationally near-instantaneous and therefore potentially of great utility for analysis of large datasets. Here we ask whether this method might have such applicability. To this end we ask whether it too recovers the two unusual results. We report that when the ML and earlier ad hoc methods disagree, YN00 recovers the results described by the ML methods, i.e. a positive correlation between GC4 and Ks and only a weak correlation between Ks and Ka. If the ML method can be trusted, then YN00 can also be considered an adequately reliable method for analysis of large datasets. Assuming this to be so we also analyze further the patterns. We show, for example, that the positive correlation between GC4 and Ks is probably in part a mutational bias, there being more methyl induced CpG-->TpG mutations in GC rich regions. As regards the evolution of isochores, it seems inappropriate to use the claimed lack of a correlation between GC and Ks as definitive evidence either against or for any model. If the positive correlation is real then, we argue, this is hard to reconcile with the biased gene conversion model for isochore formation as this predicts a negative correlation.     

Patterns of divergence among conifer ESTs and polymorphism in Pinus sylvestris identify putative selective sweeps

Finding genes that are under positive selection is a difficult task, especially in non-model organisms. Here, we have analyzed expressed sequence tag (EST) data from 4 species (Pinus pinaster, Pinus taeda, Picea glauca, and Pseudotsuga menziesii) to investigate selection patterns during their evolution and to identify genes likely to be under positive selection. To confirm selection, population samples of these genes have been sequenced in Pinus sylvestris, a species that was not included in the EST data set. The estimates of branch-specific Ka/Ks (nonsynonymous/synonymous substitution rates) across all genes in the EST data set were similar or smaller than estimates from other higher plant species. There was no evidence for the traditional indication of positive selection, Ka/Ks above 1. However, several lines of evidence based on polymorphism patterns suggest that genes with high Ka/Ks (0.20-0.52) in the EST data set are in fact more affected by positive selection in P. sylvestris than genes with low Ka/Ks (0.01-0.04). The high Ka/Ks genes have a lower level of polymorphism and more negative Tajima's D than the low Ka/Ks genes. Further, in the high Ka/Ks group, the Hudson-Kreitman-Aguade test is significant. This suggests that the EST data set is a good starting point for finding genes under positive selection in conifers and that even moderate Ka/Ks values could be indicative of selection. A group of 5 genes with high Ka/Ks collectively show evidence for positive selection within P. sylvestris.     

Molecular evolution of an imprinted gene: repeatability of patterns of evolution within the mammalian insulin-like growth factor type II receptor.

The repeatability of patterns of variation in Ka/Ks and Ks is expected if such patterns are the result of deterministic forces. We have contrasted the molecular evolution of the mammalian insulin-like growth factor type II receptor (Igf2r) in the mouse-rat comparison with that in the human-cow comparison. In so doing, we investigate explanations for both the evolution of genomic imprinting and for Ks variation (and hence putatively for mutation rate evolution). Previous analysis of Igf2r, in the mouse-rat comparison, found Ka/Ks patterns that were suggested to be contrary to those expected under the conflict theory of imprinting. We find that Ka/Ks variation is repeatable and hence confirm these patterns. However, we also find that the molecular evolution of Igf2r signal sequences suggests that positive selection, and hence conflict, may be affecting this region. The variation in Ks across Igf2r is also repeatable. To the best of our knowledge this is the first demonstration of such repeatability. We consider three explanations for the variation in Ks across the gene: (1) that it is the result of mutational biases, (2) that it is the result of selection on the mutation rate, and (3) that it is the product of selection on codon usage. Explanations 2 and 3 predict a Ka-Ks correlation, which is not found. Explanation 3 also predicts a negative correlation between codon bias and Ks, which is also not found. However, in support of explanation 1 we do find that in rodents the rate of silent C --> T mutations at CpG sites does covary with Ks, suggesting that methylation-induced mutational patterns can explain some of the variation in Ks. We find evidence to suggest that this CpG effect is due to both variation in CpG density, and to variation in the frequency with which CpGs mutate. Interestingly, however, a GC4 analysis shows no covariance with Ks, suggesting that to eliminate methyl-associated effects CpG rates themselves must be analyzed. These results suggest that, in contrast to previous studies of intragenic variation, Ks patterns are not simply caused by the same forces responsible for Ka/Ks correlations.     

Drosophila auraria复合种Period基因Thr-Gly区段的初步研究

测定了金色果蝇复合种（Drosophila auraria species complex）5个姐妹种（sibling species）和D.rufa的period(per)基因的Thr-Gly区段序列。该区段序列分析表明：DNA序列的碱其组成拥有果蝇其他基因的共同特点;颠换数多于转换数,两两种种间的颠换率与转换率的总比值为2～5,密码子第3位的颠换与转换的比值为2.5～5;同义替换/异义替换（Ks/Ka）值远大于10,且有的物种间根本不存在非同义突变,低的Ka值说明该复合种的per基因Thr-Gly区段在进化过程中可能承受着较强的选择压力。运用所得的核苷酸序列构建Drosophila auraria复合种的系统发育树,为澄清该类群的系统演化关系提供了新的线索。     

Analysis of VH gene diversity in rainbow trout (Oncorhynchus mykiss): both nonsynonymous and synonymous nucleotide changes are more frequent in CDRs than in FRs

 Forty-six immunoglobulin VH gene sequences of rainbow trout were compiled to analyze the extent of variations and the frequency of nucleotide changes in CDRs and FRs. The results show that the frequency of nonsynonymous (amino acid replacing) changes (Ka) are on average 4.9 times higher in complementarity determining regions (CDRs) than in FRs, thus contributing more diversity in CDRs. Unexpectedly, however, the frequency of synonymous (silent) changes (Ks) show the same tendency: it was 5.3 times higher in CDRs than in framework regions (FRs). The distribution of Ks/Ka ratios of each comparison shows no segregation between CDRs and FRs. The same analysis applied to five germline VH genes of Heterodontus francisci shows the same result as was found with the rainbow trout. In contrast, the results from mouse data show that, while the CDR/FR ratio for Ka is much higher (7.4), the CDR/FR ratio for Ks is only slightly higher (1.8). The distribution of Ks/Ka ratios in mouse indicates clear segregation between CDRs and FRs. This suggests that CDR germline diversity is largely generated by gene conversion in VHs of lower vertebrates such as rainbow trout or shark. This mechanism might be advantageous to lower vertebrates in generating V gene diversity faster than other mechanisms such as point mutation and selection. Received: 10 July 1996 / Revised: 9 September 1996     

A transcriptomic scan for positively selected genes in two closely related marine fishes: Sebastes caurinus and S. rastrelliger

Comparative genomic analyses can provide valuable insight into functional evolutionary divergence among closely related species. Here we employ a comparative evolutionary analysis of expressed sequence tags (ESTs) from two closely related species of marine fishes (genus Sebastes--rockfish). Sebastes is a highly diverse group of marine fishes that inhabit a wide array of marine habitats and the study of this group can provide insights into speciation in the marine environment. ESTs were developed for S. caurinus (23,668 from brain, kidney, and spleen tissues) and S. rastrelliger (11,207 from brain and pituitary tissues). Following assembly we were able to identify, with high confidence, 257 orthologous sequence pairs between the two species through a reciprocal best hit blast search. An analysis of functional divergence between orthologs revealed that 19.46% had Ka/Ks values greater than 0.5 and 8.17% had Ka/Ks values greater than one, identifying a large pool of candidate genes to further study adaptive divergence in the group. Genes with elevated Ka/Ks values belonged to the following functional categories: immune function, metabolism, longevity, and reproductive behavior, indicating that adaptive divergence in these functional groups may be important in the diversification of this group of fishes. This study provides the ground work to better understand the molecular evolution of genes involved in a radiation of marine fishes.     

Expansion Mechanisms and Functional Divergence of the Glutathione S-Transferase Family in Sorghum and Other Higher Plants

Glutathione S-transferases (GSTs) exist in various eukaryotes and function in detoxification of xenobiotics and in response to abiotic and biotic stresses. We have carried out a genome-wide survey of this gene family in 10 plant genomes. Our data show that tandem duplication has been regarded as the major expansion mechanism and both monocot and dicot plants may have practiced different expansion and evolutionary history. Non-synonymous substitutions per site (Ka) and synonymous substitutions per site (Ks) analyses showed that N- and C-terminal functional domains of GSTs (GST_N and GST_C) seem to have evolved under a strong purifying selection (Ka/Ks < 1) under different selective pressures. Differential evolutionary rates between GST_N and GST_C and high degree of expression divergence have been regarded as the major drivers for the retention of duplicated genes and the adaptability to various stresses. Expression profiling also indicated that the gene family plays a role not only in stress-related biological processes but also in the sugar-signalling pathway. Our survey provides additional annotation of the plant GST gene family and advance the understanding of plant GSTs in lineage-specific expansion and species diversification.     

Pericentromeric effects shape the patterns of divergence, retention, and expression of duplicated genes in the paleopolyploid soybean

The evolutionary forces that govern the divergence and retention of duplicated genes in polyploids are poorly understood. In this study, we first investigated the rates of nonsynonymous substitution (Ka) and the rates of synonymous substitution (Ks) for a nearly complete set of genes in the paleopolyploid soybean (Glycine max) by comparing the orthologs between soybean and its progenitor species Glycine soja and then compared the patterns of gene divergence and expression between pericentromeric regions and chromosomal arms in different gene categories. Our results reveal strong associations between duplication status and Ka and gene expression levels and overall low Ks and low levels of gene expression in pericentromeric regions. It is theorized that deleterious mutations can easily accumulate in recombination-suppressed regions, because of Hill-Robertson effects. Intriguingly, the genes in pericentromeric regions-the cold spots for meiotic recombination in soybean-showed significantly lower Ka and higher levels of expression than their homoeologs in chromosomal arms. This asymmetric evolution of two members of individual whole genome duplication (WGD)-derived gene pairs, echoing the biased accumulation of singletons in pericentromeric regions, suggests that distinct genomic features between the two distinct chromatin types are important determinants shaping the patterns of divergence and retention of WGD-derived genes.     

Strong Purifying Selection on the Odysseus Gene in Two Clades of Sibling Species of the Drosophila montium Species Subgroup

The Odysseus (OdsH) gene was duplicated from its ancestral neuron-expressed gene, unc-4, and then evolved very rapidly under strong positive Darwinian selection as a speciation gene causing hybrid-male sterility between closely related species of the Drosophila simulans clade. Has OdsH also experienced similar positive selection between Drosophila sibling species other than those of the simulans clade? We cloned and sequenced OdsH and unc-4 from two clades of the Drosophila montium species subgroup, the Drosophila lini and the Drosophila kikkawai clades. The ratios of Ka/Ks for OdsH were remarkably low between sibling species of these two clades, suggesting that OdsH has been subjected to strong purifying selection in these two clades.     

Comparison of three methods for estimating rates of synonymous and nonsynonymous nucleotide substitutions

Three frequently used methods for estimating the synonymous and nonsynonymous substitution rates (Ks and Ka) were evaluated and compared for their accuracies; these methods are denoted by LWL85, LPB93, and GY94, respectively. For this purpose, we used a codon-evolution model to obtain the expected Ka and Ks values for the above three methods and compared the values with those obtained by the three methods. We also proposed some modifications of LWL85 and LPB93 to increase their accuracies. Our computer simulations under the codon-evolution model showed that for sequences < or =300 codons, the performance of GY94 may not be reliable. For longer sequences, GY94 is more accurate for estimating the Ka/Ks ratio than the modified LPB93 and LWL85 in the majority of the cases studied. This is particularly so when k > or = 3, which is the transition/transversion (mutation) rate ratio. However, when k is approximately 2 and when the sequence divergence is relatively large, the modified LWL85 performed better than GY94 and the modified LPB93. The inferiority of LPB93 to LWL85 is surprising because LPB93 was intended to improve LWL85. Also, it has been thought that the codon-based method of GY94 is better than the heuristic method of LWL85, but our simulation results showed that in many cases, the opposite was true, even though our simulation was based on the codon-evolution model.     

Selective constraints in cold‐region wild boars may defuse the effects of small effective population size on molecular evolution of mitogenomes

Spatial range expansion during population colonization is characterized by demographic events that may have significant effects on the efficiency of natural selection. Population genetics suggests that genetic drift brought by small effective population size (N_e) may undermine the efficiency of selection, leading to a faster accumulation of nonsynonymous mutations. However, it is still unknown whether this effect might be balanced or even reversed by strong selective constraints. Here, we used wild boars and local domestic pigs from tropical (Vietnam) and subarctic region (Siberia) as animal model to evaluate the effects of functional constraints and genetic drift on shaping molecular evolution. The likelihood‐ratio test revealed that Siberian clade evolved significantly different from Vietnamese clades. Different datasets consistently showed that Siberian wild boars had lower Ka/Ks ratios than Vietnamese samples. The potential role of positive selection for branches with higher Ka/Ks was evaluated using branch‐site model comparison. No signal of positive selection was found for the higher Ka/Ks in Vietnamese clades, suggesting the interclade difference was mainly due to the reduction in Ka/Ks for Siberian samples. This conclusion was further confirmed by the result from a larger sample size, among which wild boars from northern Asia (subarctic and nearby region) had lower Ka/Ks than those from southern Asia (temperate and tropical region). The lower Ka/Ks might be due to either stronger functional constraints, which prevent nonsynonymous mutations from accumulating in subarctic wild boars, or larger N_e in Siberian wild boars, which can boost the efficacy of purifying selection to remove functional mutations. The latter possibility was further ruled out by the Bayesian skyline plot analysis, which revealed that historical N_e of Siberian wild boars was smaller than that of Vietnamese wild boars. Altogether, these results suggest stronger functional constraints acting on mitogenomes of subarctic wild boars, which may provide new insights into their local adaptation of cold resistance.     

The complete mitochondrial genomes of two species from <Emphasis Type="Italic">Sinocyclocheilus</Emphasis> (Cypriniformes: Cyprinidae) and a phylogenetic analysis within Cyprininae

We determined the complete mitochondrial DNA sequences for two species of surface- and cave-dwelling-cyprinid fishes, Sinocyclocheilus grahami and S. altishoulderus. Sequence comparison of 13 protein-coding genes shows that the mutation pattern of each single gene is quite similar to those of other vertebrate animal species. Analysis of the ratios of Ka/Ks at these loci between Sinocyclocheilus and two other cyprinid species (Cyprinus carpio and Procypris rabaudi) show that Ka/Ks ratios are differed, consistent with purifying selection and variation in functional constraint among genes. Bayesian analysis and maximum likelihood analysis of the concatenated mitochondrial protein sequences for 14 cyprinid taxa support the monophyly of the family Cyprininae, and further confirm the monophyly of the genus Sinocyclocheilus. The two Sinocyclocheilus species fall within the Cyprinion-Onychostoma lineage, including Cyprinus, Carassius, and Procypris, rather than among the Barbinae, as previously suggested on morphological grounds.