Evolutionary and demographic history among Maghrebian Medicago species (Fabaceae) based on the nucleotide sequences of the chloroplast DNA barcode trnH-psbA

Data from trnH-psbA intergenic spacer (cpDNA) were analyzed to elucidate molecular evolution within and among Maghrebian species of Medicago. The spacer highlighted a high interspecific variation and a low intraspecific diversity among species. Haplotype and nucleotide diversities revealed high level of variation. Parsimony and median-joining Network methods revealed (1) the segregation into 17 haplotypes; (2) the ancestral behaviour of the annual Medicago minima and (3) the clusters are independent of the geographic origin. The neutral evolution of Wright and Fisher is rejected since the Tajima's D values deviated from 0. Besides, the statistical analyses are in agreement with an evolution into stable populations' size.     

Mutational process in protein-coding genes of human mitochondrial genome in context of evolution of Homo genus

The human mitochondrial genome, although small in size, shows a high level of variation that differs across nucleotide groups. In this work, mutation rates in mtDNA were compared in species of the Homo genus, including humans, Neanderthals, Denisova hominins, and other primate species. It was found that more than half (56.5%) of the polymorphisms in protein-coding genes of human mtDNA are actually reverse mutations to the pre-H. sapiens state of the mitochondrial genome. Among hypervariable nucleotide positions, only a small portion of mutations are specific to H. sapiens, while the majority of mutations (both nucleotide and amino acid substitutions) result in a loss of Homo-specific variants of polymorphisms. Most commonly, polymorphism variants specific to H. sapiens arise as a result of unique forward mutations and disappear mainly due to multiple reverse mutations, including those in mutational hot spots.     

Gene diversity patterns at 10 X-chromosomal loci in humans and chimpanzees

We have investigated the pattern and extent of nucleotide diversityin 10 X-chromosomal genes where mutations are known to causemental retardation in humans. For each gene, we sequenced theentire coding region from cDNA in humans, chimpanzees, and orangutans,as well as about 3 kb of genomic DNA in 20 humans sampled worldwideand in 10 chimpanzees representing two "subspecies." Overallnucleotide diversity in these genes is about twofold lower inhumans than in chimpanzees, and nucleotide diversity withinand between species is low, suggesting that a high level offunctional constraint acts on these genes. Strikingly, we findthat a summary of the allele frequency spectrum is significantlycorrelated in humans and chimpanzees, perhaps reflecting verysimilar levels of constraint at these genes in the two species.A possible exception is FMR2, which shows a higher number ofnonsynonymous than synonymous substitutions on the human lineage,suggesting the action of positive selection.     

Synonymous and Nonsynonymous Substitutions in Genes from Gramineae: Intragenic Correlations

In this work, we have investigated the relationships between synonymous and nonsynonymous rates and base composition in coding sequences from Gramineae to analyze the factors underlying the variation in substitutional rates. We have shown that in these genes the rates of nucleotide divergence, both synonymous and nonsynonymous, are, to some extent, dependent on each other and on the base composition. In the first place, the variation in nonsynonymous rate is related to the GC level at the second codon position (the higher the GC₂ level, the higher the amino acid replacement rate). The correlation is especially strong with T₂, the coefficients being significant in the three data sets analyzed. This correlation between nonsynonymous rate and base composition at the second codon position is also detectable at the intragenic level, which implies that the factors that tend to increase the intergenic variance in nonsynonymous rates also affect the intragenic variance. On the other hand, we have shown that the synonymous rate is strongly correlated with the GC₃ level. This correlation is observed both across genes and at the intragenic level. Similarly, the nonsynonymous rate is also affected at the intragenic level by GC₃ level, like the silent rate. In fact, synonymous and nonsynonymous rates exhibit a parallel behavior in relation to GC₃ level, indicating that the intragenic patterns of both silent and amino acid divergence rates are influenced in a similar way by the intragenic variation of GC₃. This result, taken together with the fact that the number of genes displaying intragenic correlation coefficients between synonymous and nonsynonymous rates is not very high, but higher than random expectation (in the three data sets analyzed), strongly suggests that the processes of silent and amino acid replacement divergence are, at least in part, driven by common evolutionary forces in genes from Gramineae. Received: 2 July 1998 / Accepted: 18 April 1999     

Analysis of the evolutionary forces shaping mitochondrial genomes of a Neotropical malaria vector complex

Many vectors of human malaria belong to complexes of morphologically indistinguishable cryptic species. Here we report the analysis of the newly sequenced complete mitochondrial DNA molecules from six recognized or putative species of one such group, the Neotropical Anopheles albitarsis complex. The molecular evolution of these genomes had been driven by purifying selection, particularly strongly acting on the RNA genes. Directional mutation pressure associated with the strand-asynchronous asymmetric mtDNA replication mechanism may have shaped a pronounced DNA strand asymmetry in the nucleotide composition in these and other Anopheles species. The distribution of sequence polymorphism, coupled with the conflicting phylogenetic trees inferred from the mitochondrial DNA and from the published white gene fragment sequences, indicates that the evolution of the complex may have involved ancient mtDNA introgression. Six protein coding genes (nad5, nad4, cox3, atp6, cox1 and nad2) have high levels of sequence divergence and are likely informative for population genetics studies. Finally, the extent of the mitochondrial DNA variation within the complex supports the notion that the complex consists of a larger number of species than until recently believed.     

Unusual signatures of highly adaptable R-loci in closely-related Arabidopsis species

Plant resistance genes (R-genes) evolve rapidly in response to changing environments. What are the most remarkable signatures of fast adaptive genes, besides the commonly revealed rapid divergence and high non-synonymous substitution rate? Here we investigated these changes in five R-loci following recent differentiation between Arabidopsis thaliana and Arabidopsis lyrata. Extreme differences in evolutionary rates were observed: e.g., an overall 5.46-9.83-fold different nucleotide diversity at two R-loci between species, ten-fold higher non-synonymous substitution rates within one species versus the other, significantly different Ka/Ks ratios between species for the same R-gene, and high interspecific divergence at one R-locus. Particularly, we observed an elevated level of trans-specific polymorphism at one R-locus and a differentially maintained presence/absence polymorphism at another. The high frequency of ancestral polymorphisms amongst R-genes suggests that the persistence of some functional variation is an important evolutionary mechanism shaping genetic variation in R-genes, while the variation of presence/absence polymorphisms provides a potential mechanism for malleable activation of adaptive resistance response pathways. The distinct patterns among R-genes suggest that the same R-gene ortholog can be quickly shaped by different evolutionary processes, e.g., purifying selection in one species but positive selection in a closely-related species.     

Mammalian gene evolution: Nucleotide sequence divergence between mouse and rat

As a paradigm of mammalian gene evolution, the nature and extent of DNA sequence divergence between homologous protein-coding genes from mouse and rat have been investigated. The data set examined includes 363 genes totalling 411 kilobases, making this by far the largest comparison conducted between a single pair of species. Mouse and rat genes are on average 93.4% identical in nucleotide sequence and 93.9% identical in amino acid sequence. Individual genes vary substantially in the extent of nonsynonymous nucleotide substitution, as expected from protein evolution studies; here the variation is characterized. The extent of synonymous (or silent) substitution also varies considerably among genes, though the coefficient of variation is about four times smaller than for nonsynonymous substitutions. A small number of genes mapped to the X-chromosome have a slower rate of molecular evolution than average, as predicted if molecular evolution is male-driven. Base composition at silent sites varies from 33% to 95% G + C in different genes; mouse and rat homologues differ on average by only 1.7% in silent-site G + C, but it is shown that this is not necessarily due to any selective constraint on their base composition. Synonymous substitution rates and silent site base composition appear to be related (genes at intermediate G + C have on average higher rates), but the relationship is not as strong as in our earlier analyses. Rates of synonymous and nonsynonymous substitution are correlated, apparently because of an excess of substitutions involving adjacent pairs of nucleotides. Several factors suggest that synonymous codon usage in rodent genes is not subject to selection.     

Nucleotide variation in the p53 tumor-suppressor gene of voles from Chernobyl, Ukraine

The 1986 Chernobyl disaster contaminated vast regions of Ukraine and Belarus with a variety of radioactive isotopes and heavy metals. While over 90% of the radioactive isotopes have decayed into stable compounds, radiation levels in contaminated areas are still extraordinarily high. In fact, some rodents living near the reactor have internal 134,137Cs concentrations approaching 80 000 Bq/g. Several recent genetic analyses of vertebrates have illustrated that mutation rates of organisms exposed to radiation from Chernobyl are higher than in control groups, but none have studied DNA sequences. Nucleotide sequences of rodent mitochondrial genes were originally reported to have been hypervariable, but those results were subsequently retracted. Herein, I report the results of a pilot study to determine the extent of nucleotide variation at the p53 gene in four species of rodents (voles) from Chernobyl and from control sites. I sequenced a 788 bp region (coding and non-coding) of p53 in 30 different mice comprising four different species of Microtus. Nucleotide variation at the population level was due to deletions and substitutions; both were limited to introns. There were no significant differences between the number of haplotypes in radioactive and control populations (p=0.60). Rare or private alleles might have arisen due to unique mutational pressures at Chernobyl. Alternatively, natural selection might have favored one allele over others (i.e., a selective sweep). Neither scenario is strongly supported by these data. Thus, no apparent genetic effects of the Chernobyl disaster on the p53 gene of resident voles were revealed; more extensive surveys will be necessary to determine if mutation rates are indeed elevated in mice from Chernobyl. However, two salient points emerge; the first involves the utility of introns as markers for mutations in coding regions and the second considers the relative merits of cloning in mutation detection studies.     

Synonymous and Nonsynonymous Substitutions in Mammalian Genes: Intragenic Correlations

Previous investigations indicated that synonymous and nonsynonymous substitution rates are correlated in mammalian genes. In the present work, this correlation has been studied at the intragenic level using a dataset of 48 orthologous genes from species belonging to at least four different mammalian orders. The results obtained show that the intragenic variability in synonymous rates is correlated with that of nonsynonymous rates. Moreover, the variation in GC level (and especially of C level) of silent positions along each gene is correlated with the variation in synonymous rate. These results reinforce the previous conclusions that synonymous and nonsynonymous rates as well as GC levels of silent positions are to some extent under common selective constraints. Received: 10 July 1997 / Accepted: 13 August 1997     

Comparison of DNA and protein polymorphisms between humans and chimpanzees

To examine the nucleotide diversity at silent (synonymous + intron + untranslated) and non-silent (nonsynonymous) sites in chimpanzees and humans, genes at six nuclear loci from two chimpanzees were sequenced. The average silent diversity was 0.19%, which was significantly higher than that in humans (0.05%). This observation suggests a significantly larger effective population size and a higher extent of neutral polymorphism in chimpanzees than in humans. On the other hand, the non-silent nucleotide diversity is similar in both species, resulting in a larger fraction of neutral mutations at non-silent sites in humans than in chimpanzees. Other types of polymorphism data were collected from the literature or databases to examine whether or not they are consistent with the nuclear DNA sequence polymorphism observed here. The nucleotide diversity at both silent and non-silent sites in mitochondrial (mt) DNA genes was compatible with that of the nuclear genes. Microsatellite loci showed a similar high extent of heterozygosity in both species, perhaps due to the combined effect of a high mutation rate and a recent population expansion in humans. At protein loci, humans are more heterozygous than chimpanzees, and the estimated fraction of neutral alleles in humans (0.84) is much larger than that in chimpanzees (0.26). These data show that the neutral fraction in non-silent changes is relatively large in the human population. This difference may be due to a relaxation of the functional constraint against proteins in the human lineage. To evaluate this possibility, it will be necessary to examine nucleotide sequences in relation to the physiological or biochemical properties of proteins.     

Nucleotide Variation at the Gpdh Locus in the Genus Drosophila

The Gpdh locus was sequenced in a broad range of Drosophila species. In contrast to the extreme evolutionary constraint seen at the amino acid level, the synonymous sites evolve at rates comparable to those of other genes. Gpdh nucleotide sequences were used to infer a phylogenetic tree, and the relationships among the species of the obscura group were examined in detail. A survey of nucleotide polymorphism within D. pseudoobscura revealed no amino acid variation in this species. Applying a modified McDonald-Kreitman test, the amino acid divergence between species in the obscura group does not appear to be excessive, implying that drift is adequate to explain the patterns of amino acid change at this locus. In addition, the level of polymorphism at the Gpdh locus in D. pseudoobscura is comparable to that found at other loci, as determined by a Hudson-Kreitman-Aguade test. Thus, the pattern of nucleotide variation within and between species at the Gpdh locus is consistent with a neutral model.     

Trimorphic DNA variation in the receptor-like protein kinase gene in the F18L15-130 region of the wild plant Arabidopsis thaliana

DNA variation in the F18L15-130 region, which contains a receptor-like protein kinase gene, was analyzed for the wild plants Arabidopsis thaliana and Arabis gemmifera. In A. thaliana, at least three divergent sequence types were observed (trimorphism), instead of two distinct sequence types (dimorphism) detected in most of other nuclear gene regions studied for this plant species. Intragenic recombinations among the divergent sequence types generated four recombinant sequence types, although pattern of intragenic recombinations was complex. The estimated nucleotide variation in A. thaliana was the highest (pi = 0.0226 and theta = 0.0228) among genes investigated in this plant so far. The high level of nucleotide variation is due to divergence among the three distinct sequence types, each of which has a low level of nucleotide variation. Tests of Tajima, and Fu and Li did not detected deviation from neutrality, except for replacement sites, for which significantly negative Fu and Li's test value was detected. These results suggest that the receptor-like protein kinase gene in this region is functional. Possible causes for the trimorphic pattern of DNA polymorphism was discussed.     

Nucleotide sequence variation at two genes of the phenylpropanoid pathway, the FAH1 and F3H genes, in Arabidopsis thaliana

The FAH1 and F3H genes encode ferulate-5-hydroxylase and flavanone-3-hydroxylase, which are enzymes in the pathways leading to the synthesis of sinapic acid esters and flavonoids, respectively. Nucleotide variation at these genes was surveyed by sequencing a sample of 20 worldwide Arabidopsis thaliana ecotypes and one Arabidopsis lyrata spp. petraea stock. In contrast with most previously studied genes, the percentage of singletons was rather low in both the FAH1 and the F3H gene regions. There was, therefore, no footprint of a recent species expansion in the pattern of nucleotide variation in these regions. In both FAH1 and F3H, nucleotide variation was structured into two major highly differentiated haplotypes. In both genes, there was a peak of silent polymorphism in the 5' part of the coding region without a parallel increase in silent divergence. In FAH1, the peak was centered at the beginning of the second exon. In F3H, nucleotide diversity was highest at the beginning of the gene. The observed pattern of variation in both FAH1 and F3H, although suggestive of balancing selection, was compatible with a neutral model with no recombination.     

Genetic Diversity and Population Structure for the Conservation of Giant Spiny Frog(Quasipaa spinosa) Using Microsatellite Loci and Mitochondrial DNA

The giant spiny frog(Quasipaa spinosa) is an endangered species with a relatively small distribution limited to southern China and Northern Vietnam. This species is becoming increasingly threatened because of over-exploitation and habitat degradation. This study provides data on the genetic diversity and population genetic structure of the giant spiny frog to facilitate the further development of effective conservation recommendations for this economically important but threatened species. We examined 10 species-specific microsatellite loci and Cyt b genes(562 bp) collected from 13 wild populations across the entire range of this species. Results of 10 microsatellite loci analysis showed a generally high level of genetic diversity. Moreover, the genetic differentiation among all 12 populations was moderate to large(overall F_(ST) = 0.1057). A total of 51 haplotypes were identified for Cyt b, which suggests high haplotype nucleotide diversities. Phylogeographic and population structure analyses using both DNA markers suggested that the wild giant spiny frog can be divided into four distinct major clades, i.e., Northern Vietnam, Western China, Central China, and Eastern China. The clades with significant genetic divergence are reproductively isolated, as evidenced by a high number of private alleles and strong incidence of failed amplification in microsatellite loci. Our research, coupled with other studies, suggests that Q. spinosa might be a species complex within which no detectable morphological variation has been revealed. The four phylogenetic clades and some subclades with distinct geographical distribution should be regarded as independent management units for conservation purposes.     

DNA序列进化过程中核苷酸替代的非独立性研究

本文评述了DNA序列间核苷酸替代数的估计方法,并通过对七个物种中组蛋白基因的比较对DNA进化的模型进行了考察。发现H2A基因第三位点上的碱基组成在物种间变异很大,并且跟H2A基因第一位点、H4基因第一、三位点及H2A上游,下游序列中的碱基组成有强正相关,提示DNA序列进化过程中存在着物种特异的区域性约束力。可能的原因是高等真核生物中GC含量升高,或者是染色体重组使这些同源序列位于不同的等质区段,从而受到不同的选择突变压。密码内各位点上核苷酸替代的相关性分析表明不同位点的替代是非独立的,其原因可能是一次替代事件引起多个位点的变化。文中讨论了这些结果对进化树推断的意义。     

Characterization of histone H3/H4 gene region and phylogenetic affinity of Ichthyophthirius multifiliis based on H4 DNA sequence variation

We sequenced the amino-terminal third of the histone H3 and H4 genes and the intergenic region from Ichthyophthirius multifiliis. Fourteen recombinant clones of 646 bp were sequenced and the level of sequence variation detected among these clones was similar to that reported among closely related species of Tetrahymena and to levels of sequence variation detected within other ciliates. The intergenic region is 417 bp and approximately 92% AT rich, making it the longest and most AT-rich ciliate H3/H4 intergenic region yet identified. Similar to Tetrahymena, the intergenic region of Ichthyophthirius contains two CCAAT regions arranged in a complementary orientation. A neighbor-joining tree was constructed based on nucleotide sequence variation among H4 genes to evaluate evolutionary relationships within and among six classes of Ciliophora. The single shortest neighbor-joining tree depicted a sister-group relationship of Ichthyophthirius with taxa of Tetrahymenina, thereby supporting monophyly of Oligohymenophorea.     

Low codon bias and high rates of synonymous substitution in Drosophila hydei and D. melanogaster histone genes

We have evaluated codon usage bias in Drosophila histone genes and have obtained the nucleotide sequence of a 5,161-bp D. hydei histone gene repeat unit. This repeat contains genes for all five histone proteins (H1, H2a, H2b, H3, and H4) and differs from the previously reported one by a second EcoRI site. These D. hydei repeats have been aligned to each other and to the 5.0-kb (i.e., long) and 4.8-kb (i.e., short) histone repeat types from D. melanogaster. In each species, base composition at synonymous sites is similar to the average genomic composition and approaches that in the small intergenic spacers of the histone gene repeats. Accumulation of synonymous changes at synonymous sites after the species diverged is quite high. Both of these features are consistent with the relatively low codon usage bias observed in these genes when compared with other Drosophila genes. Thus, the generalization that abundantly expressed genes in Drosophila have high codon bias and low rates of silent substitution does not hold for the histone genes.      

Nucleotide diversity and linkage disequilibrium in cold-hardiness- and wood quality-related candidate genes in Douglas fir

Nuclear sequence variation and linkage disequilibrium (LD) were studied in 15 cold-hardiness- and 3 wood quality-related candidate genes in Douglas fir [Pseudotsuga menziesii (Mirb.) Franco]. This set of genes was selected on the basis of its function in other plants and collocation with cold-hardiness-related quantitative trait loci (QTL). The single-nucleotide polymorphism (SNP) discovery panel represented 24 different trees from six regions in Washington and Oregon plus parents of a segregating population used in the QTL study. The frequency of SNPs was one SNP per 46 bp across coding and noncoding regions on average. Haplotype and nucleotide diversities were also moderately high with H(d) = 0.827 +/- 0.043 and pi = 0.00655 +/- 0.00082 on average, respectively. The nonsynonymous (replacement) nucleotide substitutions were almost five times less frequent than synonymous ones and substitutions in noncoding regions. LD decayed relatively slowly but steadily within genes. Haploblock analysis was used to define haplotype tag SNPs (htSNPs). These data will help to select SNPs for association mapping, which is already in progress.     

Variable Rates of Evolution among Drosophila Opsin Genes

DNA sequences and chromosomal locations of four Drosophila pseudoobscura opsin genes were compared with those from Drosophila melanogaster, to determine factors that influence the evolution of multigene families. Although the opsin proteins perform the same primary functions, the comparisons reveal a wide range of evolutionary rates. Amino acid identities for the opsins range from 90% for Rh2 to more than 95% for Rh1 and Rh4. Variation in the rate of synonymous site substitution is especially striking: the major opsin, encoded by the Rh1 locus, differs at only 26.1% of synonymous sites between D. pseudoobscura and D. melanogaster, while the other opsin loci differ by as much as 39.2% at synonymous sites. Rh3 and Rh4 have similar levels of synonymous nucleotide substitution but significantly different amounts of amino acid replacement. This decoupling of nucleotide substitution and amino acid replacement suggests that different selective pressures are acting on these similar genes. There is significant heterogeneity in base composition and codon usage bias among the opsin genes in both species, but there are no consistent relationships between these factors and the rate of evolution of the opsins. In addition to exhibiting variation in evolutionary rates, the opsin loci in these species reveal rearrangements of chromosome elements.     

Unusually limited nucleotide sequence variation of the expressed major histocompatibility complex class I genes of a New World primate species (Saguinus oedipus)

Although major histocompatibility complex (MHC) class I molecules are, as a rule, highly polymorphic in mammalian species, those of the New World primate Saguinus oedipus (cotton-top tamarin) exhibit limited polymorphism. We have cloned and sequenced twelve MHC class I cDNAs from this species. Since cloned cotton-top tamarin cell lines express three to six MHC class I molecules, this species must have at least three functional MHC class I loci. There was, however, no evidence of locus-specific substitutions in the tamarin cDNAs. Unlike all other species studied, tamarin MHC class I cDNAs displayed limited nucleotide sequence variation. The sequence similarity between the two most divergent tamarin cDNAs was 95%. To ensure that the polymerase chain reaction (PCR) primers employed in these studies had amplified all of the tamarins' expressed MHC class I genes, we used another set of primers to amplify only exons 2 and 3 from RNA and DNA. PCR of genomic DNA resulted in the amplification of six distinct clones, of which only three were well expressed. Two of these nonexpressed genes were pseudogenes and the other was a nonclassical gene. Southern blot analysis demonstrated that the tamarin has 8–11 MHC class I genes, suggesting we had indeed cloned the majority of these genes. Cotton-top tamarins are, therefore, unique among mammalian species studied to date in that they express MHC class I molecules with limited nucleotide sequence variation.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers M38403-15.