Department of Biology, Dalhousie University, Halifax, Nova Scotia, Canada
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
In animals, the germ-lines are maintained separately from the somatic cells; therefore, the mutations in the gametes can arise only from within the germ cells. In plants, the case is quite different. In most plants, the germ-lines develop from apical cells after they have undergone many somatic cell divisions to form the body of a plant. Because of this pattern of growth, it is believed that somatic mutations, arising either from replication errors or from replication-independent causes, can be passed to the gametes; thus, these are the primary source of all new genetic variation (Klekowski 1988
, 1997
). One factor that is overlooked in this perspective, however, is that germ-lines do develop in plants and that the germ-lines could be an important source or a regulator (or both) of gamete mutations (Slatkin 1984
). For plants, the germ-line can be described as the lineage of cells that are committed to give rise to either a sperm or an egg following the transition of a vegetative apex to a reproductive apex (i.e., a separate cell lineage for the sperm and the egg). The cell lineage that ultimately gives rise to the gamete appears well before meiosis and originates at least as early as the formation of the archaesporial cell, the precursor to the archaesporium (Pennell 1988
). The germ-line therefore consists of the cells that arise after the formation of the archaesporium cell, including the pregametic sporogeneous cells, and ends with the formation of the gamete (because the origin of the archaesporium cell is unknown in most taxa, it is possible that the progenitor cells of the archaesporium cell are committed to give rise to the gamete, such that the germ-line begins before the archaesporium cell [Pennell 1988
]). Therefore, any differences between the male and the female germ-lines in plants, either replication dependent or replication independent, could lead to gender-specific mutation rates; thus, the relative evolutionary rates of male- versus female-inherited DNA should be investigated.
The rate of substitution of selectively neutral DNA (introns, intergenic DNA, synonymous sites of coding DNA) should be equal to the mutation rate (Kimura 1983
; Li 1997
). The comparison of evolutionary rates of male- versus female-inherited neutral DNA in plants could therefore help determine the relative mutation rate inherent to sperm versus eggs. Although most plants do not have sex chromosomes (Ginkgo is an exception), differences between the evolutionary rates of male- and female-inherited DNA may be revealed by comparing the genetic distance (i.e., number of nucleotide substitutions per site) at neutral sites of paternally versus maternally inherited organellar DNA, a technique that has been utilized to disclose a male bias in mutation in the mitochondrial DNA of mussels (Rawson and Hilbish 1995
; Liu, Mitton, and Wu 1996
). The evolutionary rates of male- and female-inherited organellar DNA may be compared in gymnosperms because this group of plants has the unusual characteristic that some families inherit their mitochondria or their chloroplasts (or both) along a paternal line of descent, whereas others inherit them maternally (Mogensen 1996
; fig. 1
). The most obvious approach to detect differences between male- and female-inherited organellar DNA would be to determine the genetic distances among three taxa, two that have different organellar inheritance modes (one maternal and one paternal) and an outgroup. Subsequently, the branch length to the taxon with maternal inheritance and the branch length to the taxon with paternal inheritance after their divergence can be compared. This approach is problematic, however, because any gender-specific differences will be confounded with all the other differences between the taxa that could affect the substitution rates (e.g., population size, taxon-specific mutation rates). This problem can be overcome by standardizing the branch length for each taxon (i.e., for the DNA region where inheritance mode differs) relative to the branch length from another part of the genome. For example, for each taxon being compared, the genetic distance for a particular mitochondrial region may be standardized by dividing the branch length for the mitochondrial DNA by the branch length in a chloroplastidial region. When all compared taxa have the same chloroplastidial inheritance mode, the chloroplast will then control for taxon-specific effects for mitochondrial comparisons. Similarly, male and female rates in chloroplastidial DNA can be compared by standardizing the chloroplastidial rates relative to the rates from another DNA region, such as the mitochondrial DNA, where inheritance mode does not differ. Thus, by comparing the standardized distances between two taxa with different inheritance modes, taxon-level effects on substitution rates are removed, and only proportionately higher or lower evolutionary rates associated with maternal or paternal organellar inheritance will be detected.
|
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
DNA Isolation and Amplification
Total genomic DNA was isolated from approximately 1 g of seed or leaf tissue based on the protocol of Doyle and Doyle (1990)
. This DNA was used as a template for PCR amplification of both rbcL and the trnL intron. rbcL was amplified using two pairs of overlapping primers in separate PCR reactions: (1) the forward primer 5'-ATGTCACCACAAACAGAGACTAAAGC-3' combined with reverse primer 5'-CTTCTGCTACAAATAAGAATCGATCTCTCCA-3', and (2) the forward primer 5'-TGAAAACGTGAATTCCCAACCGTTTATGCG-3' and reverse primer 5'-GCAGCAGCTAGTTCCGGGCTCCA-3', as described by Hasebe et al. (1994)
. Two separate reactions were utilized because the PCR reactions using the two end primers did not yield a PCR product for most taxa. The trnL intron was amplified with the forward primer 5'-CGAAATCGGTAGACGCTACG-3' and reverse primer 5'-GGGGATAGAGGGACTTGAAC-3' (Taberlet et al. 1991
). In some instances it was necessary to utilize the reverse primer 5'-ATTTGAACTGGTGACACGAG-3' downstream from the trnL intron. M13 forward and reverse extensions were added to the 5' ends of each primer to allow efficient direct sequencing of PCR products on a LICOR sequencer (M13 forward: 5'-CAGGACGTTGTAAAACGAC-3', M13 reverse: 5'-GGATAACAATTTCACACAGG-3'). Each PCR reaction contained 100 ng genomic DNA, 1 x PCR buffer, 1 x Q solution (Qiagen, Mississauga, Canada), 200 µM each dNTP, 200 nM of each primer, and 2 U Taq DNA polymerase (Qiagen TAQ polymerase #201203). Thermal cycling was conducted with 1 min at 94°C, 1 min at 45°C, and 1 min at 72°C for a total of 40 cycles. The PCR products were then run on 0.7% lowmelting temperature agarose gels at 4°C, the desired fragments cut out and purified using Qiagen Qiaquick Gel Extraction Kit (#28704). PCR products were then quantified on a 1% agarose gel relative to a standard, to ensure that at least 500 ng was available for sequencing.
Alignment of Sequences and the Relative Ratio Test
Amino acid sequences for the translated rbcL and coxI genes and the DNA sequences for trnL intron were each aligned for all 12 taxa using CLUSTALW (Higgins, Thompson, and Gibson 1996
), and the gaps were removed. These sequences were utilized for all further analysis to ensure that the same DNA regions were assessed in all comparisons (coxI = 441 bp, rbcL = 1,293 bp, trnL intron = 230 bp). To evaluate whether synonymous substitution rate differences between taxa remain proportional across loci, a maximum likelihood relative ratio test (Muse and Gaut 1997
) was conducted on branch lengths for coxI and rbcL (across synonymous sites only) based on the codon substitution model, as described by Muse and Gaut (1994)
, and the gymnosperm phylogeny, as described by Chaw et al. (2000
, fig. 1
). The analysis was conducted using the software package HYPHYv0.7 (Muse and Pond 2000
).
Calculation and Statistical Comparison of Standardized Distances
The principal analysis for the detection of differences in the evolutionary rates of male- versus female-inherited DNA consisted of pairwise comparisons between the standardized distance for a taxon with paternal mitochondrial inheritance (Taxon 1 in fig. 2
) versus a taxon with maternal mitochondrial inheritance (Taxon 2), and pairwise comparisons between the standardized distance for a taxon with paternal chloroplastidial inheritance (Taxon 1) versus a taxon with maternal chloroplastidial inheritance (Taxon 2; see fig. 2
). For the mitochondrial comparisons, the chloroplasts were always inherited paternally; for the chloroplastidial comparisons, the mitochondria were always inherited maternally. To calculate the standardized distance of each of the two taxa to be compared, it was first necessary to determine mitochondrial and chloroplastidial branch lengths. The branch length to each of the two taxa after their divergence was determined separately at coxI, rbcL, and the trnL intron as: d01 = (dA1 + d12 - dA2)/2 and d02 = (dA2 + d12 - dA1)/2 (Li 1997
), where d01 = the branch length from the common node to Taxon 1, d02 = the branch length from the common node to Taxon 2, dA1 = the genetic distance from the outgroup to Taxon 1, dA2 = the genetic distance from the outgroup to Taxon 2, and d12 = the genetic distance between Taxon 1 and Taxon 2. These branch lengths were then used to calculate the standardized distance of coxI relative to rbcL: Dp,coxI/rbcL = d01,coxI/d01,rbcL (taxon with paternal inheritance) and Dm,coxI/rbcL = d02,coxI/d02,rbcL (taxon with maternal inheritance) and the standardized distance of coxI relative to the trnL intron: Dp,coxI/trnL = d01,coxI/d01,trnL (taxon with paternal inheritance) and Dm,coxI/trnL = d02,coxI/d02,trnL (taxon with maternal inheritance; fig. 2
), for each of the two compared taxa. The standardized distance was then compared between Taxon 1 (paternal inheritance) and Taxon 2 (maternal inheritance), once for coxI relative to rbcL and separately for coxI relative to the trnL intron. Thus, the mitochondrial rate was standardized by the chloroplastidial rate in both the mitochondrial and chloroplastidial comparisons. Use of the same standardization method in the mitochondrial and chloroplastidial comparisons not only provides consistency but also stabilizes the ratios by having the larger value (i.e., chloroplastidial) in the denominator. Thus, in contrast to the mitochondrial comparisons, in the chloroplastidial comparisons the taxon with the higher evolutionary rate will have a lower standardized distance.
|
To quantify the size of any gender effect on mitochondrial evolutionary rates, the ratio of the mean standardized distances was determined as p/
m for each comparison. For comparisons of maternally and paternally inherited chloroplastidial DNA, in contrast, the evolutionary rate of paternally inherited DNA relative to maternally inherited DNA was quantified as
m/
p because the chloroplastidial DNA is in the denominator of the standardized distance. Thus, for both mitochondrial and chloroplastidial DNA, a ratio of mean standardized distances exceeding unity signifies a higher paternal rate.
Relative Rate Tests
To discern any absolute differences in evolutionary rates among male- and female-inherited DNA, maximum likelihood relative rate tests were conducted for: (1) each pairwise comparison between a taxon with maternal versus a taxon with paternal mitochondrial inheritance at coxI based on the codon substitution model described by Muse and Gaut (1994)
, (2) each pairwise comparison of a taxon with maternal chloroplastidial inheritance versus a taxon with paternal chloroplastidial inheritance at rbcL using the model of Muse and Gaut (1994)
and at the trnL intron based on the model described by Tamura and Nei (1993)
, using the software package HYPHYv0.7. These relative rate tests were conducted to determine whether the standardized distances allowed the detection of rate differences between male- versus female-inherited DNA that would not have been apparent from the comparison of absolute differences in rates between taxa.
Reference Taxa
For relative rate tests and the calculation of branch lengths for standardized distances, Ginkgo biloba was used as the reference taxon for all comparisons not involving itself (i.e., for any comparisons among the conifers: Araucaria spp., Cupressus corneyana, Metasequoia glyptostroboides, Podocarpus spp., Taxodium distichum, Ephedra spp., Gnetum ula, Larix spp., Pinus strobus, or Torreya nucifera) because it is a member of the most closely related family that can be used as an outgroup (fig. 1
). For those cases where evolutionary rate of G. biloba was compared to the conifers, Cycas spp. was used as the outgroup because it is a member of the most closely related family that is an outgroup to both G. biloba and the conifers (fig. 1
).
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
Overall (both mitochondrial and chloroplastidial DNA), 29 of 31 comparisons of standardized distances showed a higher neutral substitution rate in paternally than in maternally inherited DNA. A total of 26 of the 31 comparisons were individually statistically significant, all showing a higher paternal rate.
The relative ratio test between coxI and rbcL across all 12 genera was statistically significant at P < 1.0 x 10-6, indicating that synonymous substitution rates for the two DNA regions are not correlated. Relative rate tests at coxI indicated that taxa with paternal mitochondrial inheritance, namely A. heterophylla, C. corneyana, M. glyptostroboides, and P. macrophyllus, had significantly higher evolutionary rates in this gene than Larix sp., P. strobus, and G. biloba, with maternal mitochondrial inheritance (table 2 ). In addition, Ephedra spp., with maternal mitochondrial and chloroplastidial inheritance, evolved significantly faster at the mitochondrial gene coxI than each taxon with paternal mitochondrial inheritance and also evolved faster at the trnL intron than each taxon with paternal chloroplastidial inheritance (table 2 ).
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
In contrast to the mitochondrial DNA, a higher evolutionary rate in paternally inherited chloroplastidial DNA leads to lower values for the standardized paternal distance than for the standardized maternal distance, resulting from a relatively higher evolutionary rate in the trnL intron and the rbcL gene. The higher standardized distance for Larix spp. and P. strobus, with paternal chloroplastidial inheritance, than the standardized distance for Ephedra spp. and G. biloba, with maternal chloroplastidial inheritance, both for coxI relative to rbcL and for coxI relative to the trnL intron, suggests that paternally inherited chloroplastidial DNA does evolve faster than maternally inherited DNA (table 1 ). A lowered evolutionary rate in the maternally inherited chloroplastidial DNA of Ephedra spp. is remarkable, given that this taxon has a significantly higher absolute substitution rate than Larix spp. and P. strobus for both coxI and the trnL intron (table 2 ). Furthermore, G. biloba showed no significant differences in the relative rate tests when compared with Larix spp. or P. strobus at the trnL intron (table 2 ), suggesting that the use of standardized distances was again critical to the detection of differences in maternally and paternally inherited DNA and that relative rate tests alone would have impeded the detection of such differences. Nevertheless, the relatively higher evolutionary rate in paternally inherited chloroplastidial DNA suggests that sperm tends to have a relatively greater number of mutations in chloroplastidial DNA than eggs.
Magnitude of the Male Bias
The ratio of the standardized distance for taxa with paternal mitochondrial inheritance relative to the standardized distance of taxa with maternal mitochondrial inheritance (p/
m) indicates that the substitution rate of male-inherited DNA may be 6 to 10 times higher than the female rate in these gymnosperms and possibly even higher. This range is derived from the ratio of standardized distances for the comparisons between Pinaceae and the members of the clade containing Podocarpaceae, Cupressaceae, and Araucariaceae. These particular comparisons are the most apt to give accurate estimates of the number of mutations in sperm relative to eggs because these taxonomic groups have probably had contrasting inheritance modes for similar time periods (i.e., because the ancestral form is maternal mitochondrial inheritance, the paternal inheritance mode was probably derived in the lineage that gave rise to Podocarpaceae, Cupressaceae, Araucariaceae, and Taxaceae and occurred before these families diverged from one another, see fig. 1
). Nevertheless, these values represent the lowermost estimates of male-bias because the difference in inheritance modes did not occur immediately following the phylogenetic divergence.
For the comparisons of maternally versus paternally inherited chloroplastidial DNA, the ratio of standardized distances (m/
p) for the significant tests indicates that the mutation rate for the sperm is at least two- to sixfold higher than that for the eggs (table 1
), close to the lowermost estimates for mitochondrial DNA. This suggests that the male effect could potentially be lower in chloroplastidial DNA than mitochondrial DNA, a phenomenon that could arise because of physiological differences between the male and female germ-lines. For example, a higher metabolic rate or a higher number of mitochondrial replications in the male germ-cells (i.e., the product of the number of cell divisions and number of mitochondrial replications per cell division, Birky and Skavaril 1984
) could magnify male effects on the mutation rate in the mitochondrial DNA. Similar to mitochondrial comparisons, these ratios underestimate the relative number of mutations in sperm versus the eggs because the difference in inheritance modes did not occur immediately following the phylogenetic divergence of the compared taxa.
The ratio of standardized distances reflects the overall male bias in mature sperm versus mature eggs. These ratios underestimate the relative number of mutations in the male versus the female germ-line. This is because the number of mutations in the gametes could result from both germ-line mutations and somatic mutations. Therefore, the relative number of mutations in the sperm versus the eggs may differ from the relative number of mutations in the male versus the female germ-line, to the extent that somatic mutations are passed to the gametes. Furthermore, it should be noted that the ratio of standardized distances might be an imprecise measure of the magnitude of the male bias, given the high variance associated with ratios and the fact that these ratios were determined using the means of the standardized distances that themselves have variances. In this regard, the magnitude of the male bias may range considerably beyond that suggested by the ratios of standardized distances.
Possible Causes of the Lack of Significance in Some Comparisons
Even though the totality of the data suggests the presence of a higher substitution rate in paternally inherited DNA, it is noteworthy that the standardized distance for T. nucifera (with maternal mitochondrial inheritance) was not statistically significantly different from the standardized distance for C. corneyana, M. glyptostroboides, Araucaria spp., or T. distichum (with paternal inheritance) in the comparisons of maternally and paternally inherited mitochondrial DNA. It appears that maternal mitochondrial inheritance is very recently derived in T. nucifera of the Taxaceae, which occurs at a branch tip in a clade otherwise exhibiting paternal mitochondrial inheritance (Mogensen 1996
; fig. 1
) The branch length for this taxon in the comparisons to gymnosperm families with paternal mitochondrial inheritance therefore represents both a period of paternal inheritance and only a short subsequent period of maternal inheritance, a situation that may impede the detection of "maternal effects" in this particular taxon. In contrast, the lack of a significantly higher value for the standardized distance for G. biloba (with maternal chloroplastidial inheritance) than for T. nucifera (with paternal inheritance) in the chloroplastidial comparisons is inconsistent with a slow maternal rate, given the long period of divergence of maternal and paternal chloroplastidial inheritance for these two taxa. An important factor that may impede the detection of differences in this case may be "leakiness" of chloroplastidial inheritance (whereby some component of the chloroplasts are transmitted paternally, Yamada, Miyamura, and Hori 1993
), that could potentially dilute the male or female effect on substitution rates. We do not know whether this phenomenon occurs more often in certain taxa. Nevertheless, the fact that the standardized distance for T. nucifera was significantly lower than the standardized distance for Podocarpus spp. in the mitochondrial comparisons and that the standardized distance for G. biloba was considerably higher than the paternal distance for T. nucifera in the chloroplastidial comparisons suggests that the maternally inherited mitochondrial DNA evolves relatively slower and that the paternally inherited chloroplastidial DNA evolves relatively faster in T. nucifera, both consistent with male-driven evolution of DNA sequences.
Germ-Lines are Involved in Determining the Number of Mutations in Their Gametes
The evidence suggesting a gender-specific effect on mutation in these gymnosperms is consistent with the results in mussels, humans, primates, and birds, where the silent substitution rates in male- and female-inherited DNA indicate an excess of male-induced mutations (Shimmin, Chang, and Li 1993
; Chang et al. 1994
; Rawson and Hilbish 1995
; Liu, Mitton, and Wu 1996
; Agulnik et al. 1997
; Ellegren and Fridolfsson 1997
; McVean and Hurst 1997
; Kahn and Quinn 1999
; Bohosslan, Skaletsky, and Page 2000
; Huttley et al. 2000
). In animals, the higher mutation rate in males has been attributed to a higher number of cell divisions, and therefore an increased number of replication errors, in the male germ-line (Li 1997
). In humans, the female germ cells stop dividing before birth and complete with egg maturity, whereas male germ cells continuously divide, resulting in more than 32 male-divisions for each female cell division by age 45 (Vogel and Rathenberg 1975
; Crow 1997
). There is also strong evidence, however, that the higher mutation rate in males results from replication-independent differences between the germ-lines (Shimmin, Chang, and Li 1993
; McVean and Hurst 1997
; Hurst and Ellegren 1998
; Smith and Hurst 1999
; Bohosslan, Skaletsky, and Page 2000
; Huttley et al. 2000
) that may include differences in the expression of DNA repair genes, different levels of cellular mutagens, or a higher level of DNA methylation in the male germ-line (or all) (Bestor 1998
; Blackshear et al. 1998
; Smith and Hurst 1999
; Huttley et al. 2000
). In plants, it is possible that any of these factors, either replication dependent or replication independent, could account for the presence of a higher mutation rate in the sperm than in the eggs; thus, these results cannot contribute to the ongoing debate in animals. The trend toward male-driven evolution in these plants does imply, however, that the germ-lines are somehow involved in determining the number of mutations that exist in their gametes.
Analysis of X-, Y-, and autosomal-linked genes in rodents indicates that the relatively higher evolutionary rate in males is not the result of a higher number of cell divisions in the male germ-line but rather is the result of an adaptive link between the hemizygous expression of X-linked genes and a reduction in the mutation rate. For this reason, Smith and Hurst (1999)
suggested that the higher mutation rate reported in males is really the result of a lowered mutation rate in X-linked DNA than Y-linked DNA and therefore is not caused by a greater number of mutations arising in the male germ-line. For plants, there is evidence from Silene (Caryophyllaceae) indicating that the X-linked DNA has greater neutral genetic variation than Y-linked DNA. In contrast to mammals, however, it was found that this result cannot be explained by different mutation rates in X-linked and Y-linked genes (the purpose of that study was not to compare male and female rates, but it could suggest that the DNA of different chromosomes evolves at different rates, Filatov et al. 2000
, 2001
). These studies using sex chromosomes suggest that the comparison of evolutionary rates in X- and Y-linked DNA may not be the best means to evaluate whether paternally and maternally inherited DNA evolve at different rates because other explanations for rate variation are possible. The trend of a higher mutation rate in the sperm reported in this study, however, cannot be explained by a different mutation rate in different parts of the genome. This is because we compared the evolutionary rate of a single DNA region (i.e., coxI, rbcL, or the trnL intron) when it is carried by the male germ-line with the evolutionary rate when it is carried by the female germ-line. Therefore, it may be inferred that the trend of a higher rate of neutral evolution for paternally inherited DNA is the result of being carried in the male germ-line rather than in the female germ-line. The results presented here therefore suggest, at least for plants, that a fundamental difference between the male and female germ-lines may be responsible for male-driven evolution.
Possible Scenarios that Could Lead to Male-Driven Evolution in Plants
There are many possible scenarios that could explain why male-driven evolution exists in plants, each of which suggests a different role for somatic mutations. One possible scenario is that somatic cell divisions are the primary cause of mutation in plants and that these somatic mutations are passed on more often to the male germ-line than to the female germ-line, thereby causing the higher number of mutations rate in the sperm. This scenario is possible, given the complexity of the pattern of apical cell divisions, such that the male germ-line may develop from a more mutation-prone somatic cell lineage (e.g., more actively dividing) than the female germ-line. A second possible explanation for male-driven evolution in plants is that the somatic mutations are common but that these mutations are rarely transferred to the gametes. Many somatic mutations occur during the formation of internodes after the apical cells have already progressed; therefore, they are not passed to the gametes (Antolin and Strobeck 1985
). In this case, the higher number of mutations in the sperm would necessarily result from differences between the male and female germ-lines, either replication dependent or replication independent. A third possibility is that somatic mutations in plants are rare relative to germ-line mutations; therefore, they are not an important source of mutations in the gametes. According to this scenario, the male bias would have to be caused by differences, either replication dependent or replication independent, in male and female germ-lines. One way that somatic mutations could be rare events compared with germ-line mutations is if the germ-lines arise from a portion of the apical meristem that is largely quiescent during somatic growth, such that the number of somatic mutations arising from replication errors is minimal. Even though quiescent apical cells have been identified in many plant species, whether these regions give rise to the reproductive structures remains uncertain (Klekowski 1988
; Lyndon 1998
). Nevertheless, it is apparent, under each of these possible scenarios, that the germ-lines are somehow involved in determining the number of mutations that are passed to the gametes.
Given that replication-independent factors could at least partially account for male-driven evolution, it is worthwhile to consider the replication-independent differences between the male and the female germ-lines. Upregulation of the putative nucleotide excision DNA repair homologue ERCC1 in the male germ-line cells of Lilium longiflorum (Xu et al. 1998
) indicates that cells of the male germ-line may have increased DNA repair compared with the female germ-line, suggesting that there could be a greater likelihood of DNA repair errors in the sperm because of the higher levels of DNA damage requiring repair than in the eggs. Another possible replication-independent explanation for a higher mutation rate in males is that the sperm is more susceptible to damage-induced mutations because it is released from the parent plant upon maturity and is therefore exposed to solar and UV radiation, chemical mutagens, and fungal and bacterial toxins, and to dehydration (Friedberg 1985
) that could lead to increased DNAdamage-induced mutations in the male germ-line. Nevertheless, because little is currently known regarding the physiology of the male and female germ-lines in plants, it is likely that many other, yet unidentified, factors could account for a higher mutation rate in the male gametes.
Lack of Correlation Among Genes
That the relative ratio test between coxI and rbcL was statistically significant indicates that the synonymous substitution rates for these two DNA regions are not correlated among these gymnosperm taxa. This is in contrast to synonymous substitution rates observed in some studies of angiosperms, where the longer phylogenetic branch length of grasses relative to palms is proportional for three separate DNA regions, the nuclear Adh gene, the chloroplastidial rbcL gene, and the mitochondrial atpA gene (Gaut et al. 1996
; Eyre-Walker and Gaut 1997
). The correlation among different DNA regions in angiosperms suggests that organism-level effects, such as population size or taxon-specific mutation rates, play an important role in determining the evolutionary rates of selectively neutral DNA. The absence of proportionality between coxI and rbcL reported in this study, however, does not necessarily mean that organism-level effects do not exist in gymnosperms. Rather, the lack of correlation among these two genes may indicate that organism-level effects are overridden (or hidden) by gender-specific mutation rates. Organism-level effects on neutral substitution rates may be more readily detected in angiosperms because, unlike gymnosperms, the inheritance modes of organellar DNA do not extensively vary among taxa.
Significant Issues Pertaining to the Study of Gender-Specific Mutation Rates
Even though the results here suggest the presence of male-driven evolution of DNA sequences in plants, it is important to consider some of the limitations inherent to the study. In particular, it should be noted that even though the synonymous sites of coding DNA and of introns are generally believed to be selectively neutral, there is evidence that some level of selection can act on these DNA regions. For example, introns may be conserved at their splicing sites for RNA (Long and Deutsch 1999
), and synonymous sites of coding DNA can undergo some level of selection because of the biochemistry of translation and the patterns of codon usage within a gene (Akashi, Klimin, and Eyre-Walker 1998
). In this regard, it is possible that some of the substitution rates for these taxa were influenced by selective forces that could potentially dilute or magnify the results seen in some of the comparisons. An additional factor to consider is that the mitochondrial replication rate could vary among taxa. Although the mitochondria must replicate before cell division (Brown 1983
), mitochondrial replacement may occur during the lifetime of cells that have a high metabolic rate because of damage to the membrane from metabolic byproducts (Rand 1994
). This would only be an important factor, however, if the number of mitochondrial replications per unit time is associated with the number of mutations and if taxa with paternal inheritance each coincidentally have a higher metabolic rate in the male germ-line. It would be much more likely that such an association, if it existed, would not be coincidental and would suggest that a higher metabolic rate in the male germ-line causes male mutation bias. Another important consideration is that for the comparisons involving Pinaceae and for the comparisons involving Cupressaceae, more than one genus was utilized. Because Pinaceae is such a large taxonomic group and because Cupressaceae s.l. may represent more than one family (i.e., the Taxodiaceae), we concluded that it was worthwhile to include more than one genus in these comparisons. Nevertheless, it should be made clear that the comparisons involving the different genera within the Pinaceae and the different genera within the Cupressaceae are not completely independent of one another because the branch lengths up to the point of the phylogenetic divergence of the genera would have been shared. It should also be realized that although the phylogeny used is not known with certainty, a revised phylogeny would not alter our conclusions as long as Cycas and Ginkgo remained as outgroups. In addition to these three considerations, it is important to recognize that we examined mitochondrial and chloroplastidial DNA in this study because they allow the comparison of the relative evolutionary rates of paternally and maternally inherited DNA by assessing the effect of being carried in the male versus the female germ-line. Equivalent comparisons could not be conducted on nuclear DNA because it is inherited biparentally and therefore we cannot make any inference about the nuclear DNA. Further approaches that may help resolve the relative effect of the male and female germ-lines on mutation rates in plants include the examination of shorter-lived taxa that have differential organellar inheritance modes, such as Cucumis, where the ratio of somatic mutations to gender-specific mutations is apt to be lowered (Havey 1997
; Havey et al. 1998
), and the further examination of gymnosperm taxa across a greater range of mitochondrial and chloroplastidial genes. In addition, increased morphological information regarding the pattern of cell division leading to sperm and egg formation in a range of plant taxa as well as studies about the differential expression of DNA repair genes in sperm and egg will help to further resolve the factors that could cause male-driven evolution of DNA sequences in plants.
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
Keywords: male
mutation rate
plant
germ-line
mitochondria
chloroplasts
Address for correspondence and reprints: C. A. Whittle, Department of Biology, Dalhousie University, 1355 Oxford Street, Halifax, Nova Scotia, Canada B3H 4J1. whitc{at}is2.dal.ca
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Agulnik A. I., C. E. Bishop, J. L. Lerner, S. I. Agulnik, V. V. Solovyev, 1997 Analysis of mutation rates in the Smcy/Smcx genes shows that mammalian evolution is male driven Mamm. Genome 8:134-138[ISI][Medline]
Akashi H., R. M. Klimin, A. Eyre-Walker, 1998 Mutation pressure, natural selection, and the evolution of base composition in Drosophila Genetica 103:49-60[ISI]
Antolin M. F., C. Strobeck, 1985 The population genetics of somatic mutation in plants Am. Nat 126:52-62[ISI]
Besendahl A., Y. L. Qiu, J. Lee, J. D. Palmer, D. Bhattacharya, 2000 The cyanobacterial origin and vertical transmission of the plastid tRNALeu group-I intron Curr. Genet 37:12-23[ISI][Medline]
Bestor T. H., 1998 Cytosine methylation and the unequal developmental potentials of the oocyte and sperm genomes Am. J. Hum. Genet 62:1269-1273[ISI][Medline]
Bhattacharya D., B. Surek, M. Rüsing, S. Damberger, M. Melkonian, 1994 Group-I introns are inherited through common ancestry in the nuclear-encoded rRNA of Zygnematales (Charophyceae) Proc. Natl. Acad. Sci. USA 91:9916-9920
Birky C. W., R. D. Skavaril, 1984 Random partitioning of cytoplasmic organellars at cell division. The effect of organellar and cell volume J. Theor. Biol 106:441-447[ISI][Medline]
Blackshear P. E., S. M. Goldsworthy, J. F. Foley, K. A. McAllister, L. M. Bennet, N. K. Collins, D. O. Bunch, P. Brown, R. W. Wiseman, B. J. Davis, 1998 Brca1 and Brca2 expression patterns in mitotic and meiotic cells of mice Oncogene 16:161-168
Bohosslan H. B., H. Skaletsky, D. C. Page, 2000 Unexpectedly similar rates of nucleotide substitution found in male and female hominids Nature 406:622-625[ISI][Medline]
Bollback J., 2001 URL site: brahms.biology.rochester.edu/software.html
Bousquet J., S. H. Strauss, A. H. Doerksen, R. A. Price, 1992 Extensive variation in evolutionary rate of rbcL gene sequences among seed plants Proc. Natl. Acad. Sci. USA 89:7844-7848[Abstract]
Bowe L. M., C. W. dePamphilis, 1996 Effects of RNA editing and gene processing on phylogenetic reconstruction Mol. Biol. Evol 13:1159-1166[Abstract]
Brown W. M., 1983 Evolution of animal mitochondrial DNA Pp. 6288 in M. Nei and R. K. Koehn, eds. Evolution of genes and proteins. Sinauer, Sunderland, Mass
Chang B. H., L. C. Shimmin, S. K. Shyue, D. Hewett-Emmett, W. H. Li, 1994 Weak male-driven molecular evolution in rodents Proc. Natl. Acad. Sci. USA 91:827-831
Chaw S.-M., C. L. Parkinson, Y. Cheng, T. M. Vincent, J. D. Palmer, 2000 Seed plant phylogeny inferred from all three plant genomes: monophyly of extant gymnosperms and origin of Gnetales from conifers Proc. Natl. Acad. Sci. USA 97:4086-4091
Crow J. F., 1997 The high spontaneous mutation rate: is it a health risk? Proc. Natl. Acad. Sci. USA 94:8380-8386
Doyle J. J., J. L. Doyle, 1990 Isolation of plant DNA from fresh tissue Focus 12:13-15
Ellegren H., A.-K. Fridolfsson, 1997 Male-driven evolution of DNA sequences in birds Nat. Genet 17:182-184[ISI][Medline]
Eyre-Walker A., B. S. Gaut, 1997 Correlated rates of synonymous site evolution across plant genomes Mol. Biol. Evol 14:455-460[Abstract]
Fangan B. M., B. Stedje, O. E. Stabbetorp, E. S. Jensen, K. S. Jakobsen, 1994 A general approach for PCR-amplification and sequencing of chloroplast DNA form crude vascular plant and algal tissue Biotechniques 16:484-494[ISI][Medline]
Felsenstein J., 1993 PHYLIP (phylogeny inference package). Version 3.5c Distributed by the author, Department of Genetics, University of Washington, Seattle
Filatov D. A., V. Laporte, C. Vitte, D. Charlesworth, 2001 DNA diversity in sex-linked and autosomal genes of the plant species Silene latifolia and Silene dioica Mol. Biol. Evol 18:1442-1454
Filatov D. A., F. Moneger, I. Negrutiu, D. Charlesworth, 2000 Low variability in a Y-linked plant gene and its implications for Y-chromosome evolution Nature 404:388-390[ISI][Medline]
Friedberg E. C., 1985 DNA repair Freeman, New York
Gadek P. A., C. J. Quinn, 1993 An analysis of relationships within the Cupressaceae sensu stricto based on rbcL sequences Ann. MO Bot. Gard 80:581-586
Gaut B. S., B. R. Morton, B. C. McCaig, M. T. Clegg, 1996 Substitution rate comparisons between grasses and palms: synonymous rate differences at the nuclear gene Adh parallel rate differences at the plastid gene rbcL Proc. Natl. Acad. Sci. USA 93:10274-10279
Gielly L., Y. Yong-Ming, P. Küpfer, P. Taberlet, 1996 Phylogenetic use of noncoding regions in the genus Gentiana L.: chloroplast trnL (UAA) intron versus nuclear ribosomal internal transcribed spacer sequences Mol. Phylogenet. Evol 5:460-466[ISI][Medline]
Haldane J. B. S., 1947 The mutation rate of the gene for hemophila, and its segregation ratios in males and females Ann. Eugen 13:262-271[ISI]
Hasebe M., T. Omori, M. Nakazawa, T. Sano, M. Kato, K. Iwatsuki, 1994 RbcL gene-sequences provide evidence for the evolutionary lineages of leptosporangiate ferns Proc. Natl. Acad. Sci. USA 91:5730-5734[Abstract]
Havey M. J., 1997 Predominant paternal transmission of the mitochondrial genome in cucumber J. Hered 88:232-235[ISI]
Havey M. J., J. D. McCreight, B. Rhodes, G. Taurick, 1998 Differential transmission of the Cucumis organellar genomes Theor. Appl. Genet 97:122-128[ISI]
Hiesel R., B. Combettes, A. Brennicke, 1994 Evidence for RNA editing in mitochondria of all major groups of land plants except the bryophyta Proc. Natl. Acad. Sci 91:629-663[Abstract]
Higgins D. G., J. D. Thompson, T. J. Gibson, 1996 Using CLUSTAL for multiple sequence alignments Methods Enzymol 266:383-402[ISI][Medline]
Hurst L. D., H. Ellegren, 1998 Sex biases in the mutation rate Trends Genet 14:446-452[ISI][Medline]
Huttley G. A., I. B. Jacobsen, S. R. Wilson, S. Easteal, 2000 How important is DNA replication for mutagenesis? Mol. Biol. Evol 17:929-937
Kahn N. W., T. W. Quinn, 1999 Male-driven evolution among eoaves? A test of the replicative division hypothesis in a heterogametic female (ZW) system J. Mol. Evol 49:750-759[ISI][Medline]
Kimura M., 1983 The neutral theory of molecular evolution Cambridge University Press, Cambridge, UK
Klekowski E. J. Jr., 1988 Developmental selection and plant evolution Columbia University Press, New York
. 1989 Ageing and mutation in plants Nature 340:389-391[ISI]
. 1997 Somatic mutation theory of clonality Pp. 227241 in H. De Kroon and J. van Groenendael, eds. The ecology and evolution of clonal plants. Backhuys Publishers, Leiden
Kron K. A., M. W. Chase, 1993 Systematics of the Ericaceae, Empetraceae, Epacridaceae and related taxa based upon rbcL sequence data Ann. MO Bot. Gard 80:735-741
Kumar S., K. Tamura, M. Nei, 1994 MEGAmolecular evolutionary genetics analysis software for microcomputers Comput. Appl. Biosci 10:189-191[Abstract]
Li W. H., 1997 Molecular evolution Sinauer Associates Inc., Sunderland, Mass
Liu H. P., J. B. Mitton, S. K. Wu, 1996 Paternal mitochondrial DNA differentiation far exceeds maternal mitochondrial DNA and allozyme differentiation in the freshwater mussel, Anodonta grandis grandis Evolution 50:952-957[ISI]
Long M. Y., M. Deutsch, 1999 Association of intron phases with conservation at splice site sequences and the evolution of spliceosomal introns Mol. Biol. Evol 16:1528-1534[Abstract]
Lu M. Z., A. E. Szmidt, X. R. Wang, 1998 RNA editing in gymnosperms and its impact on the evolution of the mitochondrial coxI gene Plant Mol. Biol 37:225-234[ISI][Medline]
Lyndon R. F., 1998 The shoot apical meristem Cambridge University Press, Cambridge
McVean G. T., L. D. Hurst, 1997 Evidence for a selectively favourable reduction in the mutation rate of the X chromosome Nature 386:388-392[ISI][Medline]
Mogensen H. L., 1996 The hows and whys of cytoplasmic inheritance in seed plants Am. J. Bot 83:383-404[ISI]
Morgan D. R., D. E. Soltis, 1993 Phylogenetic relationships among members of the Saxifragaceae sensu lato based on rbcL sequence data Ann. MO Bot. Gard 80:631-660
Muse S. V., B. S. Gaut, 1994 A likelihood approach for comparing synonymous and nonsynonymous nucleotide substitution rates, with application to the chlorophyll genome Mol. Biol. Evol 11:715-724
. 1997 Comparing patterns of nucleotide substitution rates among chloroplast loci using the relative ratio test Genetics 146:393-399
Muse S., S. K. Pond, 2000 HYPHY: hypothesis testing using phylogenetics Available from the authors URL http://peppercat.stat.ncsu.edu/hyphy.
Odintsova M. S., N. P. Yurina, 2000 RNA editing in plant chloroplasts and mitochondria Russ. J. Plant Physiol 47:274-284[ISI]
Pennell R. I., 1988 Sporogenesis in Conifers Adv. Bot. Res 15:179-196[ISI]
Pesole G., L. R. Ceci, C. Gissi, C. Saccone, C. Quagliariello, 1996 Evolution of the nad3-rps12 gene cluster in angiosperm mitochondria: comparison of edited and unedited sequences J. Mol. Evol 43:447-452[ISI][Medline]
Rand D. M., 1994 Thermal habit, metabolic-rate and the evolution of mitochondrial DNA Trends Ecol. Evol 9:125-131[ISI]
Rawson P. D., T. J. Hilbish, 1995 Evolutionary relationships among the male and female mitochondrial DNA lineages in the Mytilus edulis species complex Mol. Biol. Evol 12:893-901[Abstract]
Shimmin L. C., B. H. J. Chang, W. Li, 1993 Male-driven evolution of DNA sequences Nature 362:745-747[ISI][Medline]
Slatkin M., 1984 Somatic mutations as an evolutionary force Pp. 1930 in P. J. Greenwood, P. H. Harvey, and M. Slatkin, eds. Evolution. Cambridge University Press, Cambridge
Smith N. G. C., L. D. Hurst, 1999 The causes of synonymous rate variation in the rodent genome: can substitution rates be used to estimate the sex bias in mutation rate Genetics 152:661-773
Soltis D. E., P. S. Soltis, M. T. Clegg, M. Durbin, 1990 Intraspecific chloroplast DNA variation: systematic and phylogenetic implications Pp. 117150 in P. S. Soltis, D. E. Soltis, and J. J. Doyle, eds. Molecular systematics of plants. Chapman & Hall, New York
Steinhauser S., S. Beckert, I. Capesius, O. Malek, V. Knoop, 1999 Plant mitochondrial RNA editing J. Mol. Evol. 48:303-312[ISI][Medline]
Szmidt A. E., M. Z Lu, X. R. Wang, 2001 Effects of RNA editing on the coxI evolution and reconstruction Euphytica 118:9-18[ISI]
Taberlet P., L. Gielly, G. Pautou, J. Bouvet, 1991 Universal primers for amplification of three non-coding regions of chloroplast DNA Plant Mol. Biol 17:1105-1109[ISI][Medline]
Tamura K., M. Nei, 1993 Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees Mol. Biol. Evol 10:512-526[Abstract]
Vogel E., R. Rathenberg, 1975 Spontaneous mutation in man Pp. 356359 in H. Harris and K. Hirschhorn eds. Advances in human genetics, Vol. 5. Plenum Press, New York
Xu H. L., I. Swoboda, P. L. Bhalla, A. M. Sijbers, C. X. Zhao, E.-K. Ong, J. H. J. Hoeijmakers, M. B. Singh, 1998 Plant homologue of human excision gene ERCC1 points to conservation of DNA repair mechanisms Plant J 13:823-829[ISI][Medline]
Yamada A., S. Miyamura, T. Hori, 1993 Cytological study on plastid inheritance of Cryptomeria japonica D Don. Plant Morphol 5:19-29
Yang Z., 2000 PAML: phylogenetic analysis by maximum likelihood. Version 3.0 University College London, London, UK
Yang Z., R. Neilsen, 2000 Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models Mol. Biol. Evol 17:32-43