A study of the phylogeny of Brassica rapa,B. nigra,Raphanus sativus,and their related genera using noncoding regions of chloroplast DNA

There are two evolutionary lineages in the genus Brassica: the rapa/oleracea lineage and the nigra lineage. Using nuclear DNA sequences such as the intergenic spacer between 5S rRNA genes and the internal transcribed spacer between 18S and 25S rRNA genes, we and others had previously demonstrated that Raphanus sativus is closely related to the nigra lineage. In the present study, we sequenced the chloroplast noncoding region between trnT and trnF and that between trnD and trnT in seven species and showed that R. sativus is more closely related to the rapa/oleracea lineage than to the nigra lineage. The conflicting results from nuclear DNA and chloroplast DNA support the hypothesis that Raphanus was derived from a hybridization between the rapa/oleracea and the nigra lineages. We estimated the date of this hybridization event to be 60% of the divergence time between the two Brassica lineages. In addition, the pattern and rate of nucleotide substitution were studied. There were more transversions than transitions in these noncoding regions, which have a high AT content. Furthermore, the proportion of transversions among the substitutions at a site increases with increasing A + T content of its two adjacent nucleotides. An influence of immediate 5(') pyrimidine on substitution pattern is also observed when both adjacent bases in the two DNA strands are A or T. The rate of nucleotide substitution in the trnL group I intron is only about one third of the rate in the nearby intergenic spacers in the trnT-trnF fragment. The rate of nucleotide substitution in the rapa/oleracea lineage is at least 1.5 times that in the nigra lineage.     

Assessing Substitution Variation Across Sites in Grass Chloroplast DNA

We assess the similarity of base substitution processes, described by empirically derived 4 × 4 matrices, using chi-square homogeneity tests. Such significance analyses allow us to assess variation in sequence evolution across sites and we apply them to matrices derived from noncoding sites in different contexts in grass chloroplast DNA. We show that there is statistically significant variation in rates and patterns of mutation among noncoding sites in different contexts and then demonstrate a similar and significant influence of context on substitutions at fourfold degenerate sites of coding regions from grass chloroplast DNA. These results show that context has the same general effect on substitution bias in coding and noncoding DNA: the A+T content of flanking bases is correlated with rate of substitution, transition bias, and GC → AT pressure, while the number of flanking pyrimidines on a single strand is correlated with a mutational bias, or skew, toward pyrimidines. Despite the similarity in general trends, however, when we compare coding and noncoding matrices we find that there is a statistically significant difference between them even when we control for context. Most noticeably, fourfold degenerate sites in coding sequences are undergoing substitution at a higher rate and there are also significant differences in the relationship between pyrimidines skew and the number of flanking pyrimidines. Possible reasons for the differences between coding and noncoding sites are discussed. Furthermore, our analysis illustrates a simple statistical way for comparing substitution processes across sites allowing us to better study variation in evolutionary processes across a genome. [Reviewing Editor: Dr. Martin Kreitman]     

Parallel rate heterogeneity in chloroplast and mitochondrial genomes of Brazil nut trees (Lecythidaceae) is consistent with lineage effects

We investigated whether relative rates of divergence were correlated between the mitochondrial and chloroplast genomes as expected under lineage effects or were genome specific as expected with locus-specific effects. Five mitochondrial noncoding regions (nad1B_C, nad4exon1_2, nad7exon2_3, nad7exon3_4, and rps14-cob) for 21 samples from Lecythidaceae were sequenced. Three chloroplast regions (rpl20-5'rps12, trnS-trnG, and psbA-trnH) were sequenced to expand the taxa in an existing data set. Absolute rates of nucleotide and insertion and deletion (indel) changes were 13 times faster in the chloroplast genome than in the mitochondrial genome. Similar indel length frequency distributions for both organelles suggested that common mechanisms were responsible for generating indels. Molecular clock tests applied to phylogenetic trees estimated from mitochondrial and chloroplast sequences revealed global rate heterogeneity of nucleotide substitution. Maximum likelihood and Tajima's 1D relative rate tests show that Lecythis zabucajo exhibited a rate acceleration for both the mitochondrial and chloroplast sequences. Whereas Eschweilera romeu-cardosoi showed a significant rate slowdown for chloroplast sequences, the mitochondrial sequences for 3 Eschweilera taxa showed evidence for a rate slowdown only when compared with L. zabucajo. Significant rate heterogeneity was also observed for indel changes in the mitochondrial genome but not for the chloroplast. The lack of mitochondrial nucleotide changes for some taxa as well as chloroplast indel homoplasy may have limited the power of relative rate tests to detect rate variation. Relative ratio tests consistently indicated rate proportionality among branch lengths between the mitochondrial and chloroplast phylogenetic trees. The relative ratio tests showed that taxa possessing rate heterogeneity had parallel relative divergence rates in both mitochondrial and chloroplast sequences as expected under lineage effects. A neutral replication-dependent model of rate heterogeneity for both nucleotide and indel changes provides a simple explanation for common patterns of rate heterogeneity across the 2 organelle genomes in Lecythidaceae. The lineage effects observed here were uncoupled from annual/perennial habit because all the species from this study are perennial.     

The Influence of Specific Neighboring Bases on Substitution Bias in Noncoding Regions of the Plant Chloroplast Genome

Substitutions occurring in noncoding sequences of the plant chloroplast genome violate the independence of sites that is assumed by substitution models in molecular evolution. The probability that a substitution at a site is a transversion, as opposed to a transition, increases significantly with increasing A + T content of the two adjacent nucleotides. In the present study, this dependency of substitutions on local context is examined further in a number of noncoding regions from the chloroplast genome of members of the grass family (Poaceae). Two features were examined; the influence of specific neighboring bases, as opposed to the general A + T content, on transversion proportion and an influence on substitutions by nucleotides other than the two immediately adjacent to the site of substitution. In both cases, a significant effect was found. In the case of specific nucleotides, transversion proportion is significantly higher at sites with a pyrimidine immediately 5′ on either strand. Substitutions at sites of the type YNR, where N is the site of substitution, have the highest rate of transversion. This specific effect is secondary to the A + T content effect such that, in terms of proportion of substitutions that are transversions, the nucleotides are ranked T > A > C > G as to their effect when they are immediately 5′ to the site of substitution. In the case of nucleotides other than the immediate neighbors, a significant influence on substitution dynamics is observed in the case where the two neighboring bases are both A and/or T. Thus, substitutions are primarily, but not exclusively, influenced by the composition of the two nucleotides that are immediately adjacent. These results indicate that the pattern of molecular evolution of the plant chloroplast genome is extremely complex as a result of a variety of inter-site dependencies. Received: 18 October 1996 / Accepted: 12 April 1997     

Evolution of the mitochondrial rps3 intron in perennial and annual angiosperms and homology to nad5 intron 1.

The plant mitochondrial rps3 intron was analyzed for substitution and indel rate variation among 15 monocot and dicot angiosperms from 10 genera, including perennial and annual taxa. Overall, the intron sequence was very conserved among angiosperms. Based on length polymorphism, 10 different alleles were identified among the 10 genera. These allelic differences were mainly attributable to large indels. An insertion of 133 nucleotides, observed in the Alnus intron was partially or completely absent in the other lineages of the family Betulaceae. This insertion was located within domain IV of the secondary-structure model of this group IIA intron. A mobile element of 47 nucleotides that showed homology to sequences located in rice rps3 intron and in intergenic plant mitochondrial genomes was found within this insertion. Both substitution and indel rates were low among the Betulaceae sequences, but substitution rates were increasingly larger than indel rates in comparisons involving more distantly related taxa. From a secondary-structure model, regions involved in helical structures were shown to be well preserved from indels as compared to substitutions, but compensatory changes were not observed among the angiosperm sequences analyzed. Using approximate divergence times based on the fossil record, substitution and indel rate heterogeneity was observed between different pairs of annual and perennial taxa. In particular, the annual petunia and primrose evolved more than 15 and 10 times faster, for substitution and indel rates respectively, than the perennial birch and alder. This is the first demonstration of an evolutionary rate difference between perennial and annual forms in noncoding DNA, lending support to neutral causes such as the generation time, population size, and speciation rate effects to explain such rate heterogeneity. Surprisingly, the sequence from the rps3 intron had a high identity with the sequence of intron 1 from the angiosperm mitochondrial nad5 gene, suggesting a common origin of these two group IIA introns.     

Compositional evolution of noncoding DNA in the human and chimpanzee genomes

We have examined the compositional evolution of noncoding DNA in the primate genome by comparison of lineage-specific substitutions observed in 1.8 Mb of genomic alignments of human, chimpanzee, and baboon with 6542 human single-nucleotide polymorphisms (SNPs) rooted using chimpanzee sequence. The pattern of compositional evolution, measured in terms of the numbers of GC-->AT and AT-->GC changes, differs significantly between fixed and polymorphic sites, and indicates that there is a bias toward fixation of AT-->GC mutations, which could result from weak directional selection or biased gene conversion in favor of high GC content. Comparison of the frequency distributions of a subset of the SNPs revealed no significant difference between GC-->AT and AT-->GC polymorphisms, although AT-->GC polymorphisms in regions of high GC segregate at slightly higher frequencies on average than GC-->AT polymorphisms, which is consistent with a fixation bias favoring high GC in these regions. However, the substitution data suggest that this fixation bias is relatively weak, because the compositional structure of the human and chimpanzee genomes is becoming homogenized, with regions of high GC decreasing in GC content and regions of low GC increasing in GC content. The rate and pattern of nucleotide substitution in 333 Alu repeats within the human-chimpanzee-baboon alignments are not significantly affected by the GC content of the region in which they are inserted, providing further evidence that, since the time of the human-chimpanzee ancestor, there has been little or no regional variation in mutation bias.     

The Role of Context-Dependent Mutations in Generating Compositional and Codon Usage Bias in Grass Chloroplast DNA

Abstract The influence of local base composition on mutations in chloroplast DNA (cpDNA) is studied in detail and the resulting, empirically derived, mutation dynamics are used to analyze both base composition and codon usage bias. A 4 × 4 substitution matrix is generated for each of the 16 possible flanking base combinations (contexts) using 17,253 noncoding sites, 1309 of which are variable, from an alignment of three complete grass chloroplast genome sequences. It is shown that substitution bias at these sites is correlated with flanking base composition and that the A+T content of these flanking sites as well as the number of flanking pyrimidines on the same strand appears to have general influences on substitution properties. The context-dependent equilibrium base frequencies predicted from these matrices are then applied to two analyses. The first examines whether or not context dependency of mutations is sufficient to generate average compositional differences between noncoding cpDNA and silent sites of coding sequences. It is found that these two classes of sites exist, on average, in very different contexts and that the observed mutation dynamics are expected to generate significant differences in overall composition bias that are similar to the differences observed in cpDNA. Context dependency, however, cannot account for all of the observed differences: although silent sites in coding regions appear to be at the equilibrium predicted, noncoding cpDNA has a significantly lower A+T content than expected from its own substitution dynamics, possibly due to the influence of indels. The second study examines the codon usage of low-expression chloroplast genes. When context is accounted for, codon usage is very similar to what is predicted by the substitution dynamics of noncoding cpDNA. However, certain codon groups show significant deviation when followed by a purine in a manner suggesting some form of weak selection other than translation efficiency. Overall, the findings indicate that a full understanding of mutational dynamics is critical to understanding the role selection plays in generating composition bias and sequence structure.     

Silent nucleotide substitutions and G+C content of some mitochondrial and bacterial genes

Summary The G+C content of DNA varies widely in different organisms, especially microorganisms. This variation is accompanied by changes in the nucleotide composition of silent positions in codons. (Silent positions are defined and explained in the text.) These changes are mostly neutral or near neutral, and appear to result from mutation pressure in the direction of increasing either A+T (AT pressure) or G+C(GC pressure) content. Variations in G+C content are also accompanied by substitutions at replacement positions in codons. These substituions produce changes in the amino acid content of homologous proteins. The examples studied were genes for 13 mitochondrial proteins in five species, and A and B genes for bacterial tryptophan synthase in four species.In microorganisms, varying AT and GC mutational pressures, presumably resulting from shifts in the DNA polymerase system, exert strong effects on molecular evolution by changing the G+C content of DNA. These effects may be greater than those of random drift. The effects of GC pressure on silent substitutions in the systems examined are several times as great as the effects on replacement substitutions.GC pressure is exerted on noncoding as well as coding regions in mitochondrial DNA. This is shown by the close correlation (correlation coefficient, 0.99) of the G+C content of the noncoding D loop of mitochondria with the G+C content of silent positions in the corresponding mitochondrial genes.     

The Nature of Nucleotide Sequence Divergence between Barley and Maize Chloroplast DNA

Analysis of a 2175-base pair (bp) SmaI-HindIII fragment of barley chloroplast DNA revealed that rbcL (the gene for the large subunit of ribulose 1,5-bisphosphate carboxylase) and atpB (the gene for the beta subunit of ATPase) are transcribed divergently and are separated by an untranscribed region of 155-166 bp. The rbcL mRNA has a 320-residue untranslated leader region, whereas the atpB mRNA has a 296- to 309-residue leader region. The sequence of these regions, together with the initial 113 bp of the atpB-coding region and the initial 1279 bp of the rbcL-coding region, is compared with the analogous maize chloroplast DNA sequences. Two classes of nucleotide differences are present, substitutions and insertions/deletions. Nucleotide substitutions show a 1.9-fold bias toward transitions in the rbcL-coding region and a 1.5-fold bias toward transitions in the noncoding region. The level of nucleotide substitutions between the barley and maize sequences is about 0.065/bp. Seventy-one percent of the substitutions in the rbcL-coding region are at the third codon position, and 95% of these are synonymous changes. Insertion/deletion events, which are confined to the noncoding regions, are not randomly distributed in these regions and are often associated with short repeated sequences. The extent of change for the noncoding regions (about 0.093 events/bp) is less than the extent of change at the third codon positions in the rbcL-coding region (about 0.135 events/bp), including insertion/delection events. Limited sequence analysis of the analogous DNA from a wild line ( Hordeum spontaneum) and a primitive Iranian barley (H. vulgare) suggested a low rate of chloroplast DNA evolution. Compared to spinach chloroplast DNA, the barley rbcL-atpB untranslated region is extremely diverged, with only the putative rbcL promoters and ribosome-binding site being extensively conserved.     

茄科植物叶绿体基因组插入、缺失和核苷酸替代的发生方式及影响

摘要: 核苷酸替代和indels(插入、缺失统称)发生是进化的重要动力。以茄科植物为研究对象, 探讨茄属中番茄和马铃薯、烟草属中绒毛状烟草和普通烟草分化时叶绿体基因组indels和核苷酸替代的发生方式, 以及这两种突变对基因组造成的影响。结果显示: indels和核苷酸替代的发生都不是随意的。indels发生在A+T丰富的区域, 1 bp indels占据总数的30%以上, 大部分indels都为低于10 bp的较短片段。核苷酸替代表现出Ts(转换)/Tv(颠换)偏差, 但T→G, A→C颠换频率却明显增加。Ts/Tv比值出现种属特异性, 番茄和马铃薯比较时替代的Ts/Tv比值低于绒毛状烟草和普通烟草比较时Ts/Tv比值。不同物种替代的(A+T)/(G+C)比值有一定差异, 从而影响基因组的(G+C)%, 此比值的差异与形成物种的生长习性有一定的关系。     

两种纺织娘线粒体基因组的比较分析

目前关于螽斯科昆虫的线粒体基因组全序列及其分子进化的研究报道很少。本研究利用L-PCR技术结合嵌套步移PCR扩增获得纺织娘Mecopoda elongata和日本纺织娘M. niponensis的线粒体基因组全序列, 同时对二者之间的碱基组成和结构特点进行了比较分析。结果显示： 纺织娘线粒体基因组（GenBank登录号JQ917910）序列全长15 284 bp, A+T含量71.8%; 日本纺织娘线粒体基因组（GenBank登录号 JQ917909）序列全长15 364 bp, A+T含量72.4%; 2种纺织娘序列长度差异主要是控制区长度不同引起（纺织娘控制区长294 bp, 日本纺织娘控制区长393 bp）。2种纺织娘基因组基因含量、 相对位置及转录方向均与其他已报道的螽斯科昆虫一致, 未发现基因重排现象; 基因组中均存在较长的间隔序列, 在trnA/trnR之间的间隔序列长度分别为63 bp与68 bp, 在trnQ/trnM之间的分别为55 bp和26 bp, 在trnS^UCN/nad1之间的均为21 bp。而最长的基因重叠区域在2种纺织娘trnC/trnW之间均为8 bp, 在atp8/atp6和nad4L/nad4L之间均为7 bp。蛋白质编码基因的碱基组成和密码子使用均具有明显的偏倚性; 除nad1和nad2以特殊的TTG作为起始密码子, cox1使用特殊的起始密码子ATGA外, 其余的10种蛋白质编码基因均使用典型的ATN作为起始密码子。在tRNA基因中, 除trnS^AGN外, 均能折叠形成典型的三叶草形二级结构。在这些tRNA基因中均存在一定数目的以G-U错配为主的碱基错配, 类似现象同样存在于其他已测定的六足动物线粒体基因组中, 表明G-U配对在线粒体基因组中很可能是一种完全正常的碱基配对方式。基因组中控制区的A+T含量略低于线粒体基因组的其他区域, 表明高A+T含量并不是该区域的必要特征。本研究结果为螽斯科系统发生关系重建积累了有价值的数据资料。     

Variation patterns of the mitochondrial 16S rRNA gene with secondary structure constraints and their application to phylogeny of cyprinine fishes (Teleostei: Cypriniformes)

The mitochondrial 16S ribosomal RNA (rRNA) gene sequences from 93 cyprinid fishes were examined to reconstruct the phylogenetic relationships within the diverse and economically important subfamily Cyprininae. Within the subfamily a biased nucleotide composition (A>T, C>G) was observed in the loop regions of the gene, and in stem regions apparent selective pressures of base pairing showed a bias in favor of G over C and T over A. The bias may be associated with transition-transversion bias. Rates of nucleotide substitution were lower in stems than in loops. Analysis of compensatory substitutions across these taxa demonstrates 68% covariation in the gene and a logical weighting factor to account for dependence in mutations for phylogenetic inference should be 0.66. Comparisons of varied stem-loop weighting schemes indicate that the down-weightings for stem regions could improve the phylogenetic analysis and the degree of non-independence of stem substitutions was not as important as expected. Bayesian inference under four models of nucleotide substitution indicated that likelihood-based phylogenetic analyses were more effective in improving the phylogenetic performance than was weighted parsimony analysis. In Bayesian analyses, the resolution of phylogenies under the 16-state models for paired regions, incorporating GTR + G + I models for unpaired regions was better than those under other models. The subfamily Cyprininae was resolved as a monophyletic group, as well as tribe Labein and several genera. However, the monophyly of the currently recognized tribes, such as Schizothoracin, Barbin, Cyprinion + Onychostoma lineages, and some genera was rejected. Furthermore, comparisons of the parsimony and Bayesian analyses and results of variable length bootstrap analysis indicates that the mitochondrial 16S rRNA gene should contain important character variation to recover well-supported phylogeny of cyprinid taxa whose divergences occurred within the recent 8 MY, but could not provide resolution power for deep phylogenies spanning 10-19 MYA.     

Neighboring base composition is strongly correlated with base substitution bias in a region of the chloroplast genome

Nucleotide sequence from a region of the chloroplast genome is presented for 12 species spanning four subfamilies of the grass family. The region contains the coding sequence for the rbcL gene and the intergenic spacer between the gene coding the large subunit of ribulose-1,5-bisphosphate carboxylase/oxygenase rbcL and the photosystem I gene psal. This intergenic spacer contains a pseudogene for rpl23 as well as two noncoding segments with different A+T contents. Using the sequence of rbcL a chloroplast phylogeny of this family was constructed by parsimony. Variable sites of the two noncoding segments were traced onto the phylogeny to study the dynamics of base substitution. This was also performed for the fourfold-degenerate sites of the rbcL gene. A wide variation in transversion/transition is observed between the two noncoding segments and between the noncoding DNA and the fourfold-degenerate sites of rbcL This variation is correlated with regional A+T content. As regional A+T content decreases, the ratio of transversions to transitions also decreases. Substitutions were then scored in relation to neighboring base composition. The composition of the two bases immediately flanking each substitution is highly correlated with the transversion/transition bias. When both the 5 and 3 flanking bases are an A or a T, transversions are observed 2.2 times as frequently as transitions. When either or both neighbors are a C or a G, the opposite trend is found; transitions are observed 1.5 times more frequently than transversions. Correspondence to: Brian R. Morton     

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.     

Genomic heterogeneity of background substitutional patterns in Drosophila melanogaster

Mutation is the underlying force that provides the variation upon which evolutionary forces can act. It is important to understand how mutation rates vary within genomes and how the probabilities of fixation of new mutations vary as well. If substitutional processes across the genome are heterogeneous, then examining patterns of coding sequence evolution without taking these underlying variations into account may be misleading. Here we present the first rigorous test of substitution rate heterogeneity in the Drosophila melanogaster genome using almost 1500 nonfunctional fragments of the transposable element DNAREP1_DM. Not only do our analyses suggest that substitutional patterns in heterochromatic and euchromatic sequences are different, but also they provide support in favor of a recombination-associated substitutional bias toward G and C in this species. The magnitude of this bias is entirely sufficient to explain recombination-associated patterns of codon usage on the autosomes of the D. melanogaster genome. We also document a bias toward lower GC content in the pattern of small insertions and deletions (indels). In addition, the GC content of noncoding DNA in Drosophila is higher than would be predicted on the basis of the pattern of nucleotide substitutions and small indels. However, we argue that the fast turnover of noncoding sequences in Drosophila makes it difficult to assess the importance of the GC biases in nucleotide substitutions and small indels in shaping the base composition of noncoding sequences.     

Local recombination and mutation effects on molecular evolution in Drosophila.

I studied the cause of the significant difference in the synonymous-substitution pattern found in the achaete-scute complex genes in two Drosophila lineages, higher codon bias in Drosophila yakuba, and lower bias in D. melanogaster. Besides these genes, the functionally unrelated yellow gene showed the same substitution pattern, suggesting a region-dependent phenomenon in the X-chromosome telomere. Because the numbers of A/T --> G/C substitutions were not significantly different from those of G/C --> A/T in the yellow noncoding regions of these species, a AT/GC mutational bias could not completely account for the synonymous-substitution biases. In contrast, we did find an approximately 14-fold difference in recombination rates in the X-chromosome telomere regions between the two species, suggesting that the reduction of recombination rates in this region resulted in the reduction of the efficacy of selection in D. melanogaster. In addition, the D. orena yellow showed a 5% increase in the G + C content at silent sites in the coding and noncoding regions since the divergence from D. erecta. This pattern was significantly different from those at the orena Adh and Amy loci. These results suggest that local changes in recombination rates and mutational pressures are contributing to the irregular synonymous-substitution patterns in Drosophila.     

The rate and pattern of nucleotide substitution in Drosophila mitochondrial DNA.

The nucleotide sequences of a segment of mitochondrial DNA (mtDNA) have been determined for nine species or subspecies of the subgenus Drosophila of the genus Drosophila. This segment contains two complete protein-coding genes (i.e., NADH dehydrogenase subunit 1 and cytochrome b) and a transfer RNA gene (tRNA(ser)). The G+C content at third-codon positions for the two protein-coding genes was 1.5 times higher than that in the D. melanogaster species group, which belongs to the subgenus Sophophora. However, there was a substantial difference between the nucleotide frequencies of G and C. The number of nucleotide substitutions per silent site was more than three times higher than that for nuclear DNA, although it was only 60% of that for mammalian mtDNA. Both parametric and nonparametric analyses revealed a strong transition-transversion bias in nucleotide substitution, as was observed in mammalian mtDNA. Moreover, the rate of substitution of A and T for G and C is higher than that for the opposite direction. This bias seems to be responsible for the extremely A+T-rich base composition of Drosophila mtDNA. It is also noted that the rate of transitional change between A and G is higher than that between T and C.     

小麦和水稻叶绿体基因组进化中核苷酸替代的种属特异性

以玉米叶绿体基因组为参照序列,采用三序列比较法系统分析了小麦和水稻分化过程中叶绿体基因组核苷酸替代的发生方式.结果表明,小麦中存在(A+T)/(G+C)替代偏差,水稻则无,该差异对小麦和水稻分化后叶绿体基因组G+C含量产生不同的影响,替代使小麦叶绿体基因组G+C含量降低、水稻叶绿体基因组G+C含量表现增加.无论在编码区、非编码区,还是不同功能基因区,小麦叶绿体基因组转换与颠换的比值都显著低于水稻.小麦和水稻叶绿体基因组进化中核苷酸替代呈现种属特异性.     

Phylogeny, rate variation, and genome size evolution of Pelargonium (Geraniaceae)

The phylogeny of 58 Pelargonium species was estimated using five plastid markers (rbcL, matK, ndhF, rpoC1, trnL-F) and one mitochondrial gene (nad5). The results confirmed the monophyly of three major clades and four subclades within Pelargonium but also indicate the need to revise some sectional classifications. This phylogeny was used to examine karyotype evolution in the genus: plotting chromosome sizes, numbers and 2C-values indicates that genome size is significantly correlated with chromosome size but not number. Accelerated rates of nucleotide substitution have been previously detected in both plastid and mitochondrial genes in Pelargonium, but sparse taxon sampling did not enable identification of the phylogenetic distribution of these elevated rates. Using the multigene phylogeny as a constraint, we investigated lineage- and locus-specific heterogeneity of substitution rates in Pelargonium for an expanded number of taxa and demonstrated that both plastid and mitochondrial genes have had accelerated substitution rates but with markedly disparate patterns. In the plastid, the exons of rpoC1 have significantly accelerated substitution rates compared to its intron and the acceleration was mainly due to nonsynonymous substitutions. In contrast, the mitochondrial gene, nad5, experienced substantial acceleration of synonymous substitution rates in three internal branches of Pelargonium, but this acceleration ceased in all terminal branches. Several lineages also have dN/dS ratios significantly greater than one for rpoC1, indicating that positive selection is acting on this gene, whereas the accelerated synonymous substitutions in the mitochondrial gene are the result of elevated mutation rates.     

Nucleotide sequence variation of the chloroplast trnK/matK region in two wild Fagopyrum (Polygonaceae) species, F. leptopodum and F. statice

Nucleotide sequence polymorphisms of the intron of the chloroplast trnK (UUU) gene, including a matK gene, were investigated within two wild Fagopyrum species, F. leptopodum and F. statice, to assess the degree and pattern of the inter- and intraspecific differences in coding and noncoding chloroplast DNA regions in higher plants. Ten and five accessions were used for F. leptopodum and F. statice, respectively. The length of the trnK intron region in these species ranged from 2494 to 2506 bp. In the trnK intron, the net nucleotide substitution number per site (Da) between the two species was 0.00109, lower than the nucleotide diversity (pi), 0.00195 for F. leptopodum and 0.00144 for F. statice, suggesting a low level of interspecific divergence. This result seems to be due to the phylogenetic pattern that both species are interspersed with each other, which was revealed by the phylogenetic analyses based on the nucleotide substitutions and indels. In the matK gene region (1524 bp), seven and two nucleotide substitutions were found within F. leptopodum and F. statice, respectively. All of the nine nucleotide substitutions (eight of which were nonsynonymous) within and between F. leptopodum and F. statice were clustered in the 5' part of the matK gene region, and no variation was found in the 3' part. This suggests that most of the 3' part is occupied by the conserved domains that are important for the binding activity of the gene product to the precursor mRNA, and therefore implies that the 3' part is more functionally constrained than the 5' part.