Male-driven evolution of mitochondrial and chloroplastidial DNA sequences in plants

When and why are mitochondria paternally inherited?

In contrast with nuclear genes that are passed on through both parents, mitochondrial genes are maternally inherited in most species, most of the time. The genetic conflict stemming from this transmission asymmetry is well-documented, and there is an abundance of population-genetic theory associated with it. While occasional or aberrant paternal inheritance occurs, there are only a few cases where exclusive paternal inheritance of mitochondrial genomes is the evolved state. Why this is remains poorly understood. By examining commonalities between species with exclusive paternal inheritance, we discuss what they may tell us about the evolutionary forces influencing mitochondrial inheritance patterns. We end by discussing recent technological advances that make exploring the causes and consequences of paternal inheritance feasible.

Phylogenetic origin of two Japanese Torreya taxa found in two regions with strongly contrasting snow depth

The Japanese archipelago exhibits a notable difference in snow depth in winter, deep snow on the Sea of Japan side and low snow cover on the Pacific Ocean side. This contrasting pattern has shaped the distribution of infraspecific taxon pairs in a range of woody plants, with taxa found on the Sea of Japan side typically exhibiting a stunted shrub form with multiple decumbent stems. The phylogenetic origin of these taxon pairs is unknown, i.e., whether the two taxa diverged from the same species or if they have different origins. This study aimed to reveal the phylogenetic origin of two varieties of Torreya nucifera (Taxaceae); var. nucifera is a tree found on the Pacific Ocean side, whereas var. radicans is a shrub found on the Sea of Japan side. We examined the phylogenetic relationships of the two varieties and worldwide Torreya taxa using whole chloroplast genomes, chloroplast DNA fragments, and the nuclear ribosomal internal transcribed spacer (ITS). The whole chloroplast genome phylogeny indicated that T. nucifera var. radicans was a sister taxon to Chinese T. grandis, rather than to var. nucifera. In contrast, the nuclear ITS phylogeny indicated that while several haplotypes of T. nucifera var. radicans were closely related to T. grandis, most haplotypes of T. nucifera var. radicans formed a single clade with those of var. nucifera. This implies that the homogenization of the ITS has occurred between the two taxa, while taxon-specific chloroplast DNA haplotypes were retained. These discordant phylogenies suggested that the two taxa have different phylogenetic origins, but have an intricate evolutionary history, involving inter-taxa hybridization and gene flow, possibly when their distributions were confined to sympatric refugia. Given the genetic evidence and distinct difference in growth form, we propose that T. nucifera var. radicans should be taxonomically treated as a distinct species, T. fruticosa.

Molecular evolution and SSRs analysis based on the chloroplast genome of Callitropsis funebris

Chloroplast genome sequences have been used to understand evolutionary events and to infer efficiently phylogenetic relationships. Callitropsis funebris (Cupressaceae) is an endemic species in China. Its phylogenetic position is controversial due to morphological characters similar to those of Cupressus, Callitropsis, and Chamaecyparis. This study used next-generation sequencing technology to sequence the complete chloroplast genome of Ca. funebris and then constructed the phylogenetic relationship between Ca. funebris and its related species based on a variety of data sets and methods. Simple sequence repeats (SSRs) and adaptive evolution analysis were also conducted. Our results showed that the monophyletic branch consisting of Ca. funebris and Cupressus tonkinensis is a sister to Cupressus, while Callitropsis is not monophyletic; Ca. nootkatensis and Ca. vietnamensis are nested in turn at the base of the monophyletic group Hesperocyparis. The statistical results of SSRs supported the closest relationship between Ca. funebris and Cupressus. By performing adaptive evolution analysis under the phylogenetic background of Cupressales, the Branch model detected three genes and the Site model detected 10 genes under positive selection; and the Branch-Site model uncovered that rpoA has experienced positive selection in the Ca. funebries branch. Molecular analysis from the chloroplast genome highly supported that Ca. funebris is at the base of Cupressus. Of note, SSR features were found to be able to shed some light on phylogenetic relationships. In short, this chloroplast genomic study has provided new insights into the phylogeny of Ca. funebris and revealed multiple chloroplast genes possibly undergoing adaptive evolution.

Population Genomic Analyses Confirm Male-Biased Mutation Rates in Snakes

Male-biased mutation rates occur in a diverse array of organisms. The ratio of male-to-female mutation rate may have major ramifications for evolution across the genome, and for sex-linked genes in particular. In ZW species, the Z chromosome is carried by males two-thirds of the time, leading to the prediction that male-biased mutation rates will have a disproportionate effect on the evolution of Z-linked genes relative to autosomes and the W chromosome. Colubroid snakes (including colubrids, elapids, and viperids) have ZW sex determination, yet male-biased mutation rates have not been well studied in this group. Here we analyze a population genomic dataset from rattlesnakes to quantify genetic variation within and genetic divergence between species. We use a new method for unbiased estimation of population genetic summary statistics to compare variation between the Z chromosome and autosomes and to calculate net nucleotide differentiation between species. We find evidence for a 2.03-fold greater mutation rate in male rattlesnakes relative to females, corresponding to an average μZ/μA ratio of 1.1. Our results from snakes are quantitatively similar to birds, suggesting that male-biased mutation rates may be a common feature across vertebrate lineages with ZW sex determination.

Inheritance of chloroplast and mitochondrial genomes in cucumber revealed by four reciprocal F1 hybrid combinations

Both genomes in chloroplasts and mitochondria of plant cell are usually inherited from maternal parent, with rare exceptions. To characterize the inheritance patterns of the organelle genomes in cucumber (Cucumis sativus var. sativus), two inbred lines and their reciprocal F₁ hybrids were analyzed using an next generation whole genome sequencing data. Their complete chloroplast genome sequences were de novo assembled, and a single SNP was identified between the parental lines. Two reciprocal F₁ hybrids have the same chloroplast genomes with their maternal parents. Meanwhile, 292 polymorphic sites were identified between mitochondrial genomes of the two parental lines, which showed the same genotypes with their paternal parents in the two reciprocal F₁ hybrids, without any recombination. The inheritance patterns of the chloroplast and mitochondria genomes were also confirmed in four additional cucumber accessions and their six reciprocal F₁ hybrids using molecular markers derived from the identified polymorphic sites. Taken together, our results indicate that the cucumber chloroplast genome is maternally inherited, as is typically observed in other plant species, whereas the large cucumber mitochondrial genome is paternally inherited. The combination of DNA markers derived from the chloroplast and mitochondrial genomes will provide a convenient system for purity test of F₁ hybrid seeds in cucumber breeding.

Reassessing Banana Phylogeny and Organelle Inheritance Modes Using Genome Skimming Data

Bananas (Musa spp.) are some of the most important fruit crops in the world, contributing up to US$10 billion in export values annually. In this study, we use high-throughput sequencing to obtain genomic resources of high-copy DNA molecules in bananas. We sampled 13 wild species and eight cultivars that represent the three genera (Ensete, Musa, and Musella) of the banana family (Musaceae). Their plastomic, 45S rDNA, and mitochondrial scaffolds were recovered from genome skimming data. Two major clades (Clades I & II) within Musa are strongly supported by the three genomic compartment data. We document, for the first time, that the plastomes of Musaceae have expanded inverted repeats (IR) after they diverged from their two close relatives, Heliconiaceae (the lobster-claws) and Strelitziaceae (the traveler's bananas). The presence/absence of rps19 within IR regions reinforces the two intra-generic clades within Musa. Our comparisons of the bananas' plastomic and mitochondrial DNA sequence trees aid in identifying hybrid bananas' parentage. As the mitochondrial genes of Musa have elevated substitution rates, paternal inheritance likely plays an influential role on the Musa mitogenome evolution. We propose genome skimming as a useful method for reliable genealogy tracing and phylogenetics in bananas.

Siberian larch (Larix sibirica Ledeb.) mitochondrial genome assembled using both short and long nucleotide sequence reads is currently the largest known mitogenome

Background

Plant mitochondrial genomes (mitogenomes) can be structurally complex while their size can vary from ~ 222 Kbp in Brassica napus to 11.3 Mbp in Silene conica. To date, in comparison with the number of plant species, only a few plant mitogenomes have been sequenced and released, particularly for conifers (the Pinaceae family). Conifers cover an ancient group of land plants that includes about 600 species, and which are of great ecological and economical value. Among them, Siberian larch (Larix sibirica Ledeb.) represents one of the keystone species in Siberian boreal forests. Yet, despite its importance for evolutionary and population studies, the mitogenome of Siberian larch has not yet been assembled and studied.

Results

Two sources of DNA sequences were used to search for mitochondrial DNA (mtDNA) sequences: mtDNA enriched samples and nucleotide reads generated in the de novo whole genome sequencing project, respectively. The assembly of the Siberian larch mitogenome contained nine contigs, with the shortest and the largest contigs being 24,767 bp and 4,008,762 bp, respectively. The total size of the genome was estimated at 11.7 Mbp. In total, 40 protein-coding, 34 tRNA, and 3 rRNA genes and numerous repetitive elements (REs) were annotated in this mitogenome. In total, 864 C-to-U RNA editing sites were found for 38 out of 40 protein-coding genes. The immense size of this genome, currently the largest reported, can be partly explained by variable numbers of mobile genetic elements, and introns, but unlikely by plasmid-related sequences. We found few plasmid-like insertions representing only 0.11% of the entire Siberian larch mitogenome.

Conclusions

Our study showed that the size of the Siberian larch mitogenome is much larger than in other so far studied Gymnosperms, and in the same range as for the annual flowering plant Silene conica (11.3 Mbp). Similar to other species, the Siberian larch mitogenome contains relatively few genes, and despite its huge size, the repeated and low complexity regions cover only 14.46% of the mitogenome sequence.

How to identify sex chromosomes and their turnover

Although sex is a fundamental component of eukaryotic reproduction, the genetic systems that control sex determination are highly variable. In many organisms the presence of sex chromosomes is associated with female or male development. Although certain groups possess stable and conserved sex chromosomes, others exhibit rapid sex chromosome evolution, including transitions between male and female heterogamety, and turnover in the chromosome pair recruited to determine sex. These turnover events have important consequences for multiple facets of evolution, as sex chromosomes are predicted to play a central role in adaptation, sexual dimorphism, and speciation. However, our understanding of the processes driving the formation and turnover of sex chromosome systems is limited, in part because we lack a complete understanding of interspecific variation in the mechanisms by which sex is determined. New bioinformatic methods are making it possible to identify and characterize sex chromosomes in a diverse array of non-model species, rapidly filling in the numerous gaps in our knowledge of sex chromosome systems across the tree of life. In turn, this growing data set is facilitating and fueling efforts to address many of the unanswered questions in sex chromosome evolution. Here, we synthesize the available bioinformatic approaches to produce a guide for characterizing sex chromosome system and identity simultaneously across clades of organisms. Furthermore, we survey our current understanding of the processes driving sex chromosome turnover, and highlight important avenues for future research.

Rapid functional divergence after small-scale gene duplication in grasses

BACKGROUND

Gene duplication has played an important role in the evolution and domestication of flowering plants. Yet little is known about how plant duplicate genes evolve and are retained over long timescales, particularly those arising from small-scale duplication (SSD) rather than whole-genome duplication (WGD) events.

RESULTS

We address this question in the Poaceae (grass) family by analyzing gene expression data from nine tissues of Brachypodium distachyon, Oryza sativa japonica (rice), and Sorghum bicolor (sorghum). Consistent with theoretical predictions, expression profiles of most grass genes are conserved after SSD, suggesting that functional conservation is the primary outcome of SSD in grasses. However, we also uncover support for widespread functional divergence, much of which occurs asymmetrically via the process of neofunctionalization. Moreover, neofunctionalization preferentially targets younger (child) duplicate gene copies, is associated with RNA-mediated duplication, and occurs quickly after duplication. Further analysis reveals that functional divergence of SSD-derived genes is positively correlated with both sequence divergence and tissue specificity in all three grass species, and particularly with anther expression in B. distachyon.

CONCLUSIONS

Our results suggest that SSD-derived grass genes often undergo rapid functional divergence that may be driven by natural selection on male-specific phenotypes. These observations are consistent with those in several animal species, suggesting that duplicate genes take similar evolutionary trajectories in plants and animals.

Complete plastome sequences of Picea asperata and P. crassifolia and comparative analyses with P. abies and P. morrisonicola

Picea asperata and P. crassifolia have sympatric ranges and are closely related, but the differences between these species at the plastome level are unknown. To better understand the patterns of variation among Picea plastomes, the complete plastomes of P. asperata and P. crassifolia were sequenced. Then, the plastomes were compared with the complete plastomes of P. abies and P. morrisonicola, which are closely and distantly related to the focal species, respectively. We also used these sequences to construct phylogenetic trees to determine the relationships among and between the four species as well as additional taxa from Pinaceae and other gymnosperms. Analysis of our sequencing data allowed us to identify 438 single nucleotide polymorphism (SNPs) point mutation events, 95 indel events, four inversion events, and seven highly variable regions, including six gene spacer regions (psbJ-petA, trnT-psaM, trnS-trnD, trnL-rps4, psaC-ccsA, and rps7-trnL) and one gene (ycf1). The highly variable regions are appropriate targets for future use in the phylogenetic reconstructions of closely related, sympatric species of Picea as well as Pinaceae in general.

Lineage-Specific Expression Divergence in Grasses Is Associated with Male Reproduction, Host-Pathogen Defense, and Domestication

Poaceae (grasses) is an agriculturally important and widely distributed family of plants with extraordinary phenotypic diversity, much of which was generated under recent lineage-specific evolution. Yet, little is known about the genes and functional modules involved in the lineage-specific divergence of grasses. Here, I address this question on a genome-wide scale by applying a novel branch-based statistic of lineage-specific expression divergence, LED, to RNA-seq data from nine tissues of the wild grass Brachypodium distachyon and its domesticated relatives Oryza sativa japonica (rice) and Sorghum bicolor (sorghum). I find that LED is generally smallest in B. distachyon and largest in O. sativa japonica, which underwent domestication earlier than S. bicolor, supporting the hypothesis that domestication may increase the rate of lineage-specific expression divergence in grasses. Moreover, in all three species, LED is positively correlated with protein-coding sequence divergence and tissue specificity, and negatively correlated with network connectivity. Further analysis reveals that genes with large LED are often primarily expressed in anther, implicating lineage-specific expression divergence in the evolution of male reproductive phenotypes. Gene ontology enrichment analysis also identifies an overrepresentation of terms related to male reproduction in the two domesticated grasses, as well as to those involved in host-pathogen defense in all three species. Last, examinations of genes with the largest LED reveal that their lineage-specific expression divergence may have contributed to antimicrobial functions in B. distachyon, to enhanced adaptation and yield during domestication in O. sativa japonica, and to defense against a widespread and devastating fungal pathogen in S. bicolor. Together, these findings suggest that lineage-specific expression divergence in grasses may increase under domestication and preferentially target rapidly evolving genes involved in male reproduction, host-pathogen defense, and the origin of domesticated phenotypes.

Phylogenomic evidence for ancient recombination between plastid genomes of the Cupressus-Juniperus-Xanthocyparis complex (Cupressaceae)

BACKGROUND

Phylogenetic relationships among Eastern Hemisphere cypresses, Western Hemisphere cypresses, junipers, and their closest relatives are controversial, and generic delimitations have been in flux for the past decade. To address relationships and attempt to produce a more robust classification, we sequenced 11 new plastid genomes (plastomes) from the five variously described genera in this complex (Callitropsis, Cupressus, Hesperocyparis, Juniperus, and Xanthocyparis) and compared them with additional plastomes from diverse members of Cupressaceae.

RESULTS

Phylogenetic analysis of protein-coding genes recovered a topology in which Juniperus is sister to Cupressus, whereas a tree based on whole plastomes indicated that the Callitropsis-Hesperocyparis-Xanthocyparis (CaHX) clade is sister to Cupressus. A sliding window analysis of site-specific phylogenetic support identified a ~ 15 kb region, spanning the genes ycf1 and ycf2, which harbored an anomalous signal relative to the rest of the genome. After excluding these genes, trees based on the remainder of the genes and genome consistently recovered a topology grouping the CaHX clade and Cupressus with strong bootstrap support. In contrast, trees based on the ycf1 and ycf2 region strongly supported a sister relationship between Cupressus and Juniperus.

CONCLUSIONS

These results demonstrate that standard phylogenomic analyses can result in strongly supported but conflicting trees. We suggest that the conflicting plastomic signals result from an ancient introgression event involving ycf1 and ycf2 that occurred in an ancestor of this species complex. The introgression event was facilitated by plastomic recombination in an ancestral heteroplasmic individual carrying distinct plastid haplotypes, offering further evidence that recombination occurs between plastomes. Finally, we provide strong support for previous proposals to recognize five genera in this species complex: Callitropsis, Cupressus, Hesperocyparis, Juniperus, and Xanthocyparis.

Organellar Genomes of White Spruce (Picea glauca): Assembly and Annotation

The genome sequences of the plastid and mitochondrion of white spruce (Picea glauca) were assembled from whole-genome shotgun sequencing data using ABySS. The sequencing data contained reads from both the nuclear and organellar genomes, and reads of the organellar genomes were abundant in the data as each cell harbors hundreds of mitochondria and plastids. Hence, assembly of the 123-kb plastid and 5.9-Mb mitochondrial genomes were accomplished by analyzing data sets primarily representing low coverage of the nuclear genome. The assembled organellar genomes were annotated for their coding genes, ribosomal RNA, and transfer RNA. Transcript abundances of the mitochondrial genes were quantified in three developmental tissues and five mature tissues using data from RNA-seq experiments. C-to-U RNA editing was observed in the majority of mitochondrial genes, and in four genes, editing events were noted to modify ACG codons to create cryptic AUG start codons. The informatics methodology presented in this study should prove useful to assemble organellar genomes of other plant species using whole-genome shotgun sequencing data.

Longevity Is Linked to Mitochondrial Mutation Rates in Rockfish: A Test Using Poisson Regression

The mitochondrial theory of ageing proposes that the cumulative effect of biochemical damage in mitochondria causes mitochondrial mutations and plays a key role in ageing. Numerous studies have applied comparative approaches to test one of the predictions of the theory: That the rate of mitochondrial mutations is negatively correlated with longevity. Comparative studies face three challenges in detecting correlates of mutation rate: Covariation of mutation rates between species due to ancestry, covariation between life-history traits, and difficulty obtaining accurate estimates of mutation rate. We address these challenges using a novel Poisson regression method to examine the link between mutation rate and lifespan in rockfish (Sebastes). This method has better performance than traditional sister-species comparisons when sister species are too recently diverged to give reliable estimates of mutation rate. Rockfish are an ideal model system: They have long life spans with indeterminate growth and little evidence of senescence, which minimizes the confounding tradeoffs between lifespan and fecundity. We show that lifespan in rockfish is negatively correlated to rate of mitochondrial mutation, but not the rate of nuclear mutation. The life history of rockfish allows us to conclude that this relationship is unlikely to be driven by the tradeoffs between longevity and fecundity, or by the frequency of DNA replications in the germline. Instead, the relationship is compatible with the hypothesis that mutation rates are reduced by selection in long-lived taxa to reduce the chance of mitochondrial damage over its lifespan, consistent with the mitochondrial theory of ageing.

Exploring the Relationships between Mutation Rates, Life History, Genome Size, Environment, and Species Richness in Flowering Plants

A new view is emerging of the interplay between mutation at the genomic level, substitution at the population level, and diversification at the lineage level. Many studies have suggested that rate of molecular evolution is linked to rate of diversification, but few have evaluated competing hypotheses. By analyzing sequences from 130 families of angiosperms, we show that variation in the synonymous substitution rate is correlated among genes from the mitochondrial, chloroplast, and nuclear genomes and linked to differences in traits among families (average height and genome size). Within each genome, synonymous rates are correlated to nonsynonymous substitution rates, suggesting that increasing the mutation rate results in a faster rate of genome evolution. Substitution rates are correlated with species richness in protein-coding sequences from the chloroplast and nuclear genomes. These data suggest that species traits contribute to lineage-specific differences in the mutation rate that drive both synonymous and nonsynonymous rates of change across all three genomes, which in turn contribute to greater rates of divergence between populations, generating higher rates of diversification. These observations link mutation in individuals to population-level processes and to patterns of lineage divergence.

Why are most organelle genomes transmitted maternally?

Why the DNA-containing organelles, chloroplasts, and mitochondria, are inherited maternally is a long standing and unsolved question. However, recent years have seen a paradigm shift, in that the absoluteness of uniparental inheritance is increasingly questioned. Here, we review the field and propose a unifying model for organelle inheritance. We argue that the predominance of the maternal mode is a result of higher mutational load in the paternal gamete. Uniparental inheritance evolved from relaxed organelle inheritance patterns because it avoids the spread of selfish cytoplasmic elements. However, on evolutionary timescales, uniparentally inherited organelles are susceptible to mutational meltdown (Muller's ratchet). To prevent this, fall-back to relaxed inheritance patterns occurs, allowing low levels of sexual organelle recombination. Since sexual organelle recombination is insufficient to mitigate the effects of selfish cytoplasmic elements, various mechanisms for uniparental inheritance then evolve again independently. Organelle inheritance must therefore be seen as an evolutionary unstable trait, with a strong general bias to the uniparental, maternal, mode.

Outcrossing rates and organelle inheritance estimated from two natural populations of the Japanese endemic conifer Sciadopitys verticillata

The Japanese endemic conifer Sciadopitys verticillata is one of the most phylogenetically isolated species of all plants. Occurring in small and scattered populations, the species is currently classified as Near Threatened by the International Union for Conservation of Nature and Natural Resources (IUCN) and as Vulnerable in three prefectures of Japan. This study investigated two major factors that should impact the genetic structure of the species at both the nuclear and organelle DNA level, the mating system and the inheritance of both the chloroplast and mitochondrial genomes. The mating system is crucial to determining the degree of outcrossing of plant species and thus should have a key role in shaping the species' population level genetic diversity and gene flow between populations but as yet has not been studied in S. verticillata. Nine mother trees and their seedling progeny from two natural populations were genotyped using genetic markers from three plant genomes (eight nuclear microsatellites and DNA sequence for the chloroplast and mitochondria). Using a maximum likelihood method implemented in the software MLTR, the study found an outcrossing rate in the seedling stage of 0.49 and 0.79 for Aburazaka and Mount Shirotori populations, respectively, and an average of 0.66 at the species level. These outcrossing rates were low for conifers and therefore may have potential deleterious implications for the conservation of the species. The test of organelle inheritance supported paternal transmission of both the chloroplast and mitochondria consistent with previous microscopic evidence.

Selection-driven evolution of sex-biased genes is consistent with sexual selection in Arabidopsis thaliana

Sex-biased genes are genes with a preferential or specific expression in one sex and tend to show an accelerated rate of evolution in animals. Various hypotheses--which are not mutually exclusive--have been put forth to explain observed patterns of rapid evolution. One possible explanation is positive selection, but this has been shown only in few animal species and mostly for male-specific genes. Here, we present a large-scale study that investigates evolutionary patterns of sex-biased genes in the predominantly self-fertilizing plant Arabidopsis thaliana. Unlike most animal species, A. thaliana does not possess sex chromosomes, its flowers develop both male and female sexual organs, and it is characterized by low outcrossing rates. Using cell-specific gene expression data, we identified genes whose expression is enriched in comparison with all other tissues in the male and female gametes (sperm, egg, and central cell), as well as in synergids, pollen, and pollen tubes, which also play an important role in reproduction. Genes specifically expressed in gametes and synergids show higher rates of protein evolution compared with the genome-wide average and no evidence for positive selection. In contrast, pollen- and pollen tube-specific genes not only have lower rates of protein evolution but also exhibit a higher proportion of adaptive amino acid substitutions. We show that this is the result of increased levels of purifying and positive selection among genes with pollen- and pollen tube-specific expression. The increased proportion of adaptive substitutions cannot be explained by the fact that pollen- and pollen tube-expressed genes are enriched in segmental duplications, are on average older, or have a larger effective population size. Our observations are consistent with prezygotic sexual selection as a result of interactions during pollination and pollen tube growth such as pollen tube competition.

Does the mode of plastid inheritance influence plastid genome architecture?

Plastid genomes show an impressive array of sizes and compactnesses, but the forces responsible for this variation are unknown. It has been argued that species with small effective genetic population sizes are less efficient at purging excess DNA from their genomes than those with large effective population sizes. If true, one may expect the primary mode of plastid inheritance to influence plastid DNA (ptDNA) architecture. All else being equal, biparentally inherited ptDNAs should have a two-fold greater effective population size than those that are uniparentally inherited, and thus should also be more compact. Here, we explore the relationship between plastid inheritance pattern and ptDNA architecture, and consider the role of phylogeny in shaping our observations. Contrary to our expectations, we found no significant difference in plastid genome size or compactness between ptDNAs that are biparentally inherited relative to those that are uniparentally inherited. However, we also found that there was significant phylogenetic signal for the trait of mode of plastid inheritance. We also found that paternally inherited ptDNAs are significantly smaller (n = 19, p = 0.000001) than those that are maternally, uniparentally (when isogamous), or biparentally inherited. Potential explanations for this observation are discussed.

Female parthenogenetic apomixis and androsporogenetic parthenogenesis in embryonal cells of Araucaria angustifolia: interpolation of progenesis and asexual heterospory in an artificial sporangium

Cell fate, development timing and occurrence of reproductive versus apomictic development in gymnosperms are shown to be influenced by culture conditions in vitro. In this study, female parthenogenetic apomixis (fPA), androsporogenetic parthenogenesis (mAP) and progenesis were demonstrated using embryonal initials of Araucaria angustifolia in scaled-up cell suspensions passing through a single-cell bottleneck in darkness and in an artificial sporangium (AS). Expression was based on defined nutrition, hormones and feedforward-adaptive feedback process controls at 23-25 °C and in darkness. In fPA, the nucleus of an embryonal initial undergoes endomitosis and amitosis, forming a diploid egg-equivalent and an apoptotic ventral canal nucleus in a transdifferentiated archegonial tube. Discharge of egg-equivalent cells as parthenospores and their dispersal into the aqueous culture medium were followed by free-nuclear conifer-type proembryogenesis. This replaced the plesiomorphic and central features of proembryogenesis in Araucariaceae. Protoplasmic fusions of embryonal initials were used to reconstruct heterokaryotic expressions of fPA in multiwell plates. In mAP, restitutional meiosis (automixis) was responsible for androsporogenesis and the discharge of monads, dyads, tetrads and polyads. In a display of progenesis, reproductive development was brought to an earlier ontogenetic stage and expressed by embryonal initials. Colchicine increased polyploidy, but androspore formation became aberrant and fragmented. Aberrant automixis led to the formation of chromosomal bouquets, which contributed to genomic silencing in embryonal initials, cytomixis and the formation of pycnotic micronucleated cells. Dispersal of female and male parthenospores displayed heteromorphic asexual heterospory in an aqueous environment.

Is rarity of pinedrops (Pterospora andromedea) in eastern North America linked to rarity of its unique fungal symbiont?

Like other myco-heterotrophic plants, Pterospora andromedea (pinedrops) is dependent upon its specific fungal symbionts for survival. The rarity of pinedrops fungal symbiont was investigated in the eastern United States where pinedrops are rare. Wild populations of eastern pinedrops were sampled, and the plant haplotypes and fungal symbionts were characterized with molecular techniques; these data were compared to those from the West with phylogenetic analyses. The frequency of the fungal symbiont in eastern white pine forests was assessed using a laboratory soil bioassay and in situ pinedrops seed baiting. Only one plant haplotype and fungal symbiont was detected. The plant haplotype was not unique to the East. The fungal symbiont appears to be a new species within the genus Rhizopogon, closely related to the western symbionts. This fungal species was not frequent in soils with or without pinedrops, but was less frequent in the latter and in comparison to the fungal symbionts in western forests. Seed baiting resulted in few germinants, suggesting that mycelial networks produced by the eastern fungal symbiont were rare. Results suggest that eastern pinedrops rarity is influenced by the distribution and rarity of its fungal symbiont.

Plant DNA barcodes and the influence of gene flow

Success of species assignment using DNA barcodes has been shown to vary among plant lineages because of a wide range of different factors. In this study, we confirm the theoretical prediction that gene flow influences species assignment with simulations and a literature survey. We show that the genome experiencing the highest gene flow is, in the majority of the cases, the best suited for species delimitation. Our results clearly suggest that, for most angiosperm groups, plastid markers will not be the most appropriate for use as DNA barcodes. We therefore advocate shifting the focus from plastid to nuclear markers to achieve an overall higher success using DNA barcodes.

Recent acceleration of plastid sequence and structural evolution coincides with extreme mitochondrial divergence in the angiosperm genus Silene

The angiosperm genus Silene exhibits some of the most extreme and rapid divergence ever identified in mitochondrial genome architecture and nucleotide substitution rates. These patterns have been considered mitochondrial specific based on the absence of correlated changes in the small number of available nuclear and plastid gene sequences. To better assess the relationship between mitochondrial and plastid evolution, we sequenced the plastid genomes from four Silene species with fully sequenced mitochondrial genomes. We found that two species with fast-evolving mitochondrial genomes, S. noctiflora and S. conica, also exhibit accelerated rates of sequence and structural evolution in their plastid genomes. The nature of these changes, however, is markedly different from those in the mitochondrial genome. For example, in contrast to the mitochondrial pattern, which appears to be genome wide and mutationally driven, the plastid substitution rate accelerations are restricted to a subset of genes and preferentially affect nonsynonymous sites, indicating that altered selection pressures are acting on specific plastid-encoded functions in these species. Indeed, some plastid genes in S. noctiflora and S. conica show strong evidence of positive selection. In contrast, two species with more slowly evolving mitochondrial genomes, S. latifolia and S. vulgaris, have correspondingly low rates of nucleotide substitution in plastid genes as well as a plastid genome structure that has remained essentially unchanged since the origin of angiosperms. These results raise the possibility that common evolutionary forces could be shaping the extreme but distinct patterns of divergence in both organelle genomes within this genus.

Female parthenogenetic apomixis and androsporogenesis in Douglas-fir embryonal initials in an artificial sporangium

Control of female parthenogenetic apomixis and androsporogenesis of Douglas-fir embryonal initials was studied using an experimental culture system in which changes in growth condition can mediate changes in cell identity and outcomes. This culture system constitutes an artificial sporangium in which myriad culture conditions can be simulated and should be applicable for research on a variety of gymnosperms. In this study, embryonal initials from developing seeds from two Douglas-fir trees were rescued and became reprogrammed for female parthenogenetic apomixis (fPA) and parthenogenetic androsporogenesis (mPA). Female PA was initiated by endomitosis forming a binucleate cell with a diploid egg-equivalent and an apoptotic ventral canal nucleus in an archegonial tube. Egg-equivalent nuclei formed cells (parthenotes) that were discharged into an aqueous culture medium. Parthenotes developed axial tiers atypical of early embryogenesis in seeds. Earlier in the year, androsporangial tubes were parthenogenetically differentiated and released monads, dyads, triads, and tetrads into the culture medium. Spores showed chromosomal aberrations. PA demonstrated a temporal separation in gender expression (dichogamy). Embryonal initials brought forward and by-passed the long juvenile phases normally needed for cells to develop into trees and express reproductive maturity. Expressions of fPA and mPA indicated that the specialized culture flasks served as an artificial sporangium (AS). Awareness is raised for the value of an AS for research in gymnosperm life cycles and as a teaching and research laboratory.

Male mutation bias and possible long-term effects of human activities

The ability of a population to adapt to changing environments depends critically on the amount and kind of genetic variability it possesses. Mutations are an important source of new genetic variability and may lead to new adaptations, especially if the population size is large. Mutation rates are extremely variable between and within species, and males usually have higher mutation rates as a result of elevated rates of male germ cell division. This male bias affects the overall mutation rate. We examined the factors that influence male mutation bias, and focused on the effects of classical life-history parameters, such as the average age at reproduction and elevated rates of sperm production in response to sexual selection and sperm competition. We argue that human-induced changes in age at reproduction or in sexual selection will affect male mutation biases and hence overall mutation rates. Depending on the effective population size, these changes are likely to influence the long-term persistence of a population.

Evolutionary rates in Veronica L. (Plantaginaceae): disentangling the influence of life history and breeding system

The evolutionary rate at which DNA sequences evolve is known to differ between different groups of organisms. However, the reasons for these different rates are seldom known. Among plants, the generation-time hypothesis, which states that organisms that reproduce faster also have more DNA substitutions per time, has gained most popularity. We evaluate the generation-time hypothesis using 131 DNA sequences from the plastid trnLF region and the nuclear ribosomal ITS region of the genus Veronica (Plantaginaceae). We also examine the alternative hypothesis that a higher substitution rate is correlated with selfing breeding system. Selfing is associated with annual life history in many organisms and may thus often be the underlying reason for observed correlations of annual life history with other characters. We provide evidence that annual life history is more likely to be the responsible factor for higher substitution rates in Veronica than a selfing breeding system. Nevertheless, the way in which annual life history may influence substitution rate in detail remains unknown, and some possibilities are discussed.

Relative rates of synonymous substitutions in the mitochondrial, chloroplast and nuclear genomes of seed plants

Previous studies have estimated that, in angiosperms, the synonymous substitution rate of chloroplast genes is three times higher than that of mitochondrial genes and that of nuclear genes is twelve times higher than that of mitochondrial genes. Here we used 12 genes in 27 seed plant species to investigate whether these relative rates of substitutions are common to diverse seed plant groups. We find that the overall relative rate of synonymous substitutions of mitochondrial, chloroplast and nuclear genes of all seed plants is 1:3:10, that these ratios are 1:2:4 in gymnosperms but 1:3:16 in angiosperms and that they go up to 1:3:20 in basal angiosperms. Our results show that the mitochondrial, chloroplast and nuclear genomes of seed plant groups have different synonymous substitutions rates, that these rates are different in different seed plant groups and that gymnosperms have smaller ratios than angiosperms.

Evolutionary rate variation at multiple levels of biological organization in plant mitochondrial DNA

We examined patterns of mitochondrial polymorphism and divergence in the angiosperm genus Silene and found substantial variation in evolutionary rates among species and among lineages within species. Moreover, we found corresponding differences in the amount of polymorphism within species. We argue that, along with our earlier findings of rate variation among genes, these patterns of rate heterogeneity at multiple phylogenetic scales are most likely explained by differences in underlying mutation rates. In contrast, no rate variation was detected in nuclear or chloroplast loci. We conclude that mutation rate heterogeneity is a characteristic of plant mitochondrial sequence evolution at multiple biological scales and may be a crucial determinant of how much polymorphism is maintained within species. These dramatic patterns of variation raise intriguing questions about the mechanisms driving and maintaining mutation rate heterogeneity in plant mitochondrial genomes. Additionally, they should alter our interpretation of many common phylogenetic and population genetic analyses.

Sex: differences in mutation, recombination, selection, gene flow, and genetic drift

In many instances, there are large sex differences in mutation rates, recombination rates, selection, rates of gene flow, and genetic drift. Mutation rates are often higher in males, a difference that has been estimated both directly and indirectly. The higher male mutation rate appears related to the larger number of cell divisions in male lineages but mutation rates also appear gene- and organism-specific. When there is recombination in only one sex, it is always the homogametic sex. When there is recombination in both sexes, females often have higher recombination but there are many exceptions. There are a number of hypotheses to explain the sex differences in recombination. Sex-specific differences in selection may result in stable polymorphisms or for sex chromosomes, faster evolutionary change. In addition, sex-dependent selection may result in antagonistic pleiotropy or sexually antagonistic genes. There are many examples of sex-specific differences in gene flow (dispersal) and a number of adaptive explanations for these differences. The overall effective population size (genetic drift) is dominated by the lower sex-specific effective population size. The mean of the mutation, recombination, and gene flow rates over the two sexes can be used in a population genetics context unless there are sex-specific differences in selection or genetic drift. Sex-specific differences in these evolutionary factors appear to be unrelated to each other. The evolutionary explanations for sex-specific differences for each factor are multifaceted and, in addition, explanations may include chance, nonadaptive differences, or mechanistic, nonevolutionary factors.

Sperm competition can drive a male-biased mutation rate

A pattern of male-biased mutation has been found in a wide range of species. The standard explanation for this bias is that there are greater numbers of mitotic cell divisions in the history of the average sperm, compared to the average egg, and that mutations typically result from errors made during replication. However, this fails to provide an ultimate evolutionary explanation for why the male germline would tolerate more mutations that are typically deleterious. One possibility is that if there is a tradeoff between producing large numbers of sperm and expending energetic resources in maintaining a lower mutation rate, sperm competition would select for males that produce larger numbers of sperm despite a higher resulting mutation rate. Here I describe a model that jointly considers the fitness consequences of deleterious mutation and mating success in the face of sperm competition. I show that a moderate level of sperm competition can account for the observation that the male germline tolerates a higher mutation rate than the female germline.

An unusual case of gender-associated mitochondrial DNA heteroplasmy: the mytilid Musculista senhousia (Mollusca Bivalvia)

BACKGROUND

Doubly Uniparental Inheritance (DUI) represents the most outstanding exception to matrilinear inheritance of mitochondrial DNA (mtDNA), typical of Metazoa. In a few bivalve mollusks, two sex-linked mtDNAs (the so-called M and F) are inherited in a peculiar way: both daughters and sons receive their F from the mother, whereas sons inherit M from the father (males do not transmit F to their progeny). This realizes a double mechanism of transmission, in which M and F mtDNAs are inherited uniparentally. DUI systems represent a unique experimental model for testing the evolutionary mechanisms that apply to mitochondrial genomes and their transmission patterns as well as to mtDNA recombination.

RESULTS

A new case of DUI is described in Musculista senhousia (Mollusca: Bivalvia: Mytilidae). Its heteroplasmy pattern is in line with standard DUI. Sequence variability analysis evidenced two main results: F haplotypes sequence variability is higher than that of M haplotypes, and F mitochondrial haplotypes experience a higher mutation rate in males' somatic tissues than in females' ones. Phylogenetic analysis revealed also that M. senhousia M and F haplotypes cluster separately from that of the other mytilids.

CONCLUSION

Sequence variability analysis evidenced some unexpected traits. The inverted variability pattern (the F being more variable than M) was new and it challenges most of the rationales proposed to account for sex-linked mtDNA evolution. We tentatively related this to the history of the Northern Adriatic populations analyzed. Moreover, F sequences evidenced a higher mutation level in male's soma, this variability being produced de novo each generation. This suggests that mechanisms evolved to protect mtDNA in females (f.i. antioxidant gene complexes) might be under relaxed selection in males. Phylogenetic analysis of sex-linked haplotypes confirmed that they have switched their roles during the evolutionary history of mytilids, at variance to what has been observed in unionids. Consequently, reciprocal monophyly of M and F lineages got easily lost because of role-reversals and consequent losses of M lineages, as already observed in Mytilus.

Extensive variation in synonymous substitution rates in mitochondrial genes of seed plants

Background

It has long been known that rates of synonymous substitutions are unusually low in mitochondrial genes of flowering and other land plants. Although two dramatic exceptions to this pattern have recently been reported, it is unclear how often major increases in substitution rates occur during plant mitochondrial evolution and what the overall magnitude of substitution rate variation is across plants.

Results

A broad survey was undertaken to evaluate synonymous substitution rates in mitochondrial genes of angiosperms and gymnosperms. Although most taxa conform to the generality that plant mitochondrial sequences evolve slowly, additional cases of highly accelerated rates were found. We explore in detail one of these new cases, within the genus Silene. A roughly 100-fold increase in synonymous substitution rate is estimated to have taken place within the last 5 million years and involves only one of ten species of Silene sampled in this study. Examples of unusually slow sequence evolution were also identified. Comparison of the fastest and slowest lineages shows that synonymous substitution rates vary by four orders of magnitude across seed plants. In other words, some plant mitochondrial lineages accumulate more synonymous change in 10,000 years than do others in 100 million years. Several perplexing cases of gene-to-gene variation in sequence divergence within a plant were uncovered. Some of these probably reflect interesting biological phenomena, such as horizontal gene transfer, mitochondrial-to-nucleus transfer, and intragenomic variation in mitochondrial substitution rates, whereas others are likely the result of various kinds of errors.

Conclusion

The extremes of synonymous substitution rates measured here constitute by far the largest known range of rate variation for any group of organisms. These results highlight the utility of examining absolute substitution rates in a phylogenetic context rather than by traditional pairwise methods. Why substitution rates are generally so low in plant mitochondrial genomes yet occasionally increase dramatically remains mysterious.

Gender-specific selection on codon usage in plant genomes

Background

Currently, there is little data available regarding the role of gender-specific gene expression on synonymous codon usage (translational selection) in most organisms, and particularly plants. Using gender-specific EST libraries (with > 4000 ESTs) from Zea mays and Triticum aestivum, we assessed whether gender-specific gene expression per se and gender-specific gene expression level are associated with selection on codon usage.

Results

We found clear evidence of a greater bias in codon usage for genes expressed in female than in male organs and gametes, based on the variation in GC content at third codon positions and the frequency of species-preferred codons. This finding holds true for both highly and for lowly expressed genes. In addition, we found that highly expressed genes have greater codon bias than lowly expressed genes for both female- and male-specific genes. Moreover, in both species, genes with female-specific expression show a greater usage of species-specific preferred codons for each of the 18 amino acids having synonymous codons. A supplemental analysis of Brassica napus suggests that bias in codon usage could also be higher in genes expressed in male gametophytic tissues than in heterogeneous (flower) tissues.

Conclusion

This study reports gender-specific bias in codon usage in plants. The findings reported here, based on the analysis of 1 497 876 codons, are not caused either by differences in the biological functions of the genes or by differences in protein lengths, nor are they likely attributable to mutational bias. The data are best explained by gender-specific translational selection. Plausible explanations for these findings and the relevance to these and other organisms are discussed.

Asymmetric postmating isolation: Darwin's corollary to Haldane's rule

Asymmetric postmating isolation, where reciprocal interspecific crosses produce different levels of fertilization success or hybrid sterility/inviability, is very common. Darwin emphasized its pervasiveness in plants, but it occurs in all taxa assayed. This asymmetry often results from Dobzhansky-Muller incompatibilities (DMIs) involving uniparentally inherited genetic factors (e.g., gametophyte-sporophyte interactions in plants or cytoplasmic-nuclear interactions). Typically, unidirectional (U) DMIs act simultaneously with bidirectional (B) DMIs between autosomal loci that affect reciprocal crosses equally. We model both classes of two-locus DMIs to make quantitative and qualitative predictions concerning patterns of isolation asymmetry in parental species crosses and in the hybrid F(1) generation. First, we find conditions that produce expected differences. Second, we present a stochastic analysis of DMI accumulation to predict probable levels of asymmetry as divergence time increases. We find that systematic interspecific differences in relative rates of evolution for autosomal vs. nonautosomal loci can lead to different expected F(1) fitnesses from reciprocal crosses, but asymmetries are more simply explained by stochastic differences in the accumulation of U DMIs. The magnitude of asymmetry depends primarily on the cumulative effects of U vs. B DMIs (which depend on heterozygous effects of DMIs), the average number of DMIs required to produce complete reproductive isolation (more asymmetry occurs when fewer DMIs are required), and the shape of the function describing how fitness declines as DMIs accumulate. Comparing our predictions to data from diverse taxa indicates that unidirectional DMIs, specifically involving sex chromosomes, cytoplasmic elements, and maternal effects, are likely to play an important role in postmating isolation.

Characteristics, causes and evolutionary consequences of male-biased mutation

Mutation has traditionally been considered a random process, but this paradigm is challenged by recent evidence of divergence rate heterogeneity in different genomic regions. One facet of mutation rate variation is the propensity for genetic change to correlate with the number of germ cell divisions, reflecting the replication-dependent origin of many mutations. Haldane was the first to connect this association of replication and mutation to the difference in the number of cell divisions in oogenesis (low) and spermatogenesis (usually high), and the resulting sex difference in the rate of mutation. The concept of male-biased mutation has been thoroughly analysed in recent years using an evolutionary approach, in which sequence divergence of autosomes and/or sex chromosomes are compared to allow inference about the relative contribution of mothers and fathers in the accumulation of mutations. For instance, assuming that a neutral sequence is analysed, that rate heterogeneity owing to other factors is cancelled out by the investigation of many loci and that the effect of ancestral polymorphism is properly taken into account, the male-to-female mutation rate ratio, alpham, can be solved from the observed difference in rate of X and Y chromosome divergence. The male mutation bias is positively correlated with the relative excess of cell divisions in the male compared to the female germ line, as evidenced by a generation time effect: in mammals, alpham is estimated at approximately 4-6 in primates, approximately 3 in carnivores and approximately 2 in small rodents. Another life-history correlate is sexual selection: when there is intense sperm competition among males, increased sperm production will be associated with a larger number of mitotic cell divisions in spermatogenesis and hence an increase in alpham. Male-biased mutation has implications for important aspects of evolutionary biology such as mate choice in relation to mutation load, sexual selection and the maintenance of genetic diversity despite strong directional selection, the tendency for a disproportionate large role of the X (Z) chromosome in post-zygotic isolation, and the evolution of sex.

The ovotestis: an underdeveloped organ of evolution

In animals that have separate sexes (gonochorists), many sperm are produced to fertilise a few eggs. As the male germline undergoes more mitoses, so the accumulated mutation frequency is elevated in sperm compared with ova, and evolution is 'male-driven'. In contrast, in many hermaphroditic animals, a single organ--the ovotestis--produces both ova and sperm. Since self-renewing cells in the ovotestis may give rise to both cell types throughout life, ova in hermaphrodites could in theory have undergone as many cell divisions as sperm. Here, I consider some possible effects of the ovotestis on evolution. In particular, I hypothesise that the accumulated mutation frequency of nuclear genes in hermaphrodites (including species that change sex) may reach twice that compared with gonochorists. There may be an even greater increase in the mitochondrial mutation frequency. Further developmental studies and the accumulation of comparative data should allow hypothesis testing. If the prediction is correct, then it may provide the most-straightforward explanation for the extraordinary diversity of mitochondrial DNA in some hermaphrodites, especially molluscs.

Mammalian Male Mutation Bias: Impacts of Generation Time and Regional Variation in Substitution Rates

In mammals, males undergo a greater number of germline cell divisions compared with females. Thus, the male germline accumulates more DNA replication errors, which result in male mutation bias-a higher mutation rate for males than for females. The phenomenon of male mutation bias has been investigated mostly for rodents and primates, however, it has not been studied in detail for other mammalian orders. Here we sequenced and analyzed five introns of three genes (DBX/DBY, UTX/UTY, and ZFX/ZFY) homologous between X and Y chromosomes in several species of perissodactyls (horses and rhinos) and of primates. Male mutation bias was evident: substitution rate was higher for a Y chromosome intron than for its X chromosome homologue for all five intron pairs studied. Substitution rates varied regionally among introns sequenced on the same chromosome and this variation influenced male mutation bias inferred from each intron pair. Interestingly, we observed a positive correlation in substitution rates between homologous X and homologous Y introns as well as between orthologous primate and perissodactyl introns. The male-to-female mutation rate ratio estimated from concatenated sequences of five perissodactyl introns was 3.88 (95% CI = 2.90-6.07). Using the data generated here and estimates available in the literature, we compared male mutation bias among several mammalian orders. We conclude that male mutation bias is significantly higher for organisms with long generation times (primates, perissodactyls, and felids) than for organisms with short generation times (e.g., rodents) since the former undergo a greater number of male germline cell divisions.

Substitution Rate Heterogeneity and the Male Mutation Bias

Germline mutation rates have been found to be higher in males than in females in many organisms, a likely consequence of cell division being more frequent in spermatogenesis than in oogenesis. If the majority of mutations are due to DNA replication error, the male-to-female mutation rate ratio (alpha(m)) is expected to be similar to the ratio of the number of germ line cell divisions in males and females (c), an assumption that can be tested with proper estimates of alpha(m) and c. Alpha(m) is usually estimated by comparing substitution rates in putatively neutral sequences on the sex chromosomes. However, substantial regional variation in substitution rates across chromosomes may bias estimates of alpha(m) based on the substitution rates of short sequences. To investigate regional substitution rate variation, we estimated sequence divergence in 16 gametologous introns located on the Z and W chromosomes of five bird species of the order Galliformes. Intron ends and potentially conserved blocks were excluded to reduce the effect of using sequences subject to negative selection. We found significant substitution rate variation within Z chromosome (G15 = 37.6, p = 0.0010) as well as within W chromosome introns (G15 = 44.0, p = 0.0001). This heterogeneity also affected the estimates of alpha(m), which varied significantly, from 1.53 to 3.51, among the introns (ANOVA: F(13,14) = 2.68, p = 0.04). Our results suggest the importance of using extensive data sets from several genomic regions to avoid the effects of regional mutation rate variation and to ensure accurate estimates of alpha(m).

The influence of environmental factors, the pollen : ovule ratio and seed bank persistence on molecular evolutionary rates in plants

One of the main goals of molecular evolutionary biology is to determine the factors that influence the evolutionary rate of selectively neutral DNA, but much remains unknown, especially for plants. Key factors that could alter the mutation rate include environmental tolerances (because they reflect a plants vulnerability to changes in habitat), the pollen:ovule ratio (as it is associated with the number of mitotic divisions) and seed longevity (because this influences the number of generations per unit time in plants). This is the first study to demonstrate that seed bank persistence and drought tolerance are positively associated with molecular evolutionary rates in plants and that pollen:ovule ratio, shade tolerance and salinity tolerance have no detectable relationship. The implications of the findings to our understanding of the impact of environmental agents, the number of cell divisions and cell aging on neutral DNA sequence evolution are discussed.

Multiple major increases and decreases in mitochondrial substitution rates in the plant family Geraniaceae

BACKGROUND

Rates of synonymous nucleotide substitutions are, in general, exceptionally low in plant mitochondrial genomes, several times lower than in chloroplast genomes, 10-20 times lower than in plant nuclear genomes, and 50-100 times lower than in many animal mitochondrial genomes. Several cases of moderate variation in mitochondrial substitution rates have been reported in plants, but these mostly involve correlated changes in chloroplast and/or nuclear substitution rates and are therefore thought to reflect whole-organism forces rather than ones impinging directly on the mitochondrial mutation rate. Only a single case of extensive, mitochondrial-specific rate changes has been described, in the angiosperm genus Plantago.

RESULTS

We explored a second potential case of highly accelerated mitochondrial sequence evolution in plants. This case was first suggested by relatively poor hybridization of mitochondrial gene probes to DNA of Pelargonium hortorum (the common geranium). We found that all eight mitochondrial genes sequenced from P. hortorum are exceptionally divergent, whereas chloroplast and nuclear divergence is unexceptional in P. hortorum. Two mitochondrial genes were sequenced from a broad range of taxa of variable relatedness to P. hortorum, and absolute rates of mitochondrial synonymous substitutions were calculated on each branch of a phylogenetic tree of these taxa. We infer one major, approximately 10-fold increase in the mitochondrial synonymous substitution rate at the base of the Pelargonium family Geraniaceae, and a subsequent approximately 10-fold rate increase early in the evolution of Pelargonium. We also infer several moderate to major rate decreases following these initial rate increases, such that the mitochondrial substitution rate has returned to normally low levels in many members of the Geraniaceae. Finally, we find unusually little RNA editing of Geraniaceae mitochondrial genes, suggesting high levels of retroprocessing in their history.

CONCLUSION

The existence of major, mitochondrial-specific changes in rates of synonymous substitutions in the Geraniaceae implies major and reversible underlying changes in the mitochondrial mutation rate in this family. Together with the recent report of a similar pattern of rate heterogeneity in Plantago, these findings indicate that the mitochondrial mutation rate is a more plastic character in plants than previously realized. Many molecular factors could be responsible for these dramatic changes in the mitochondrial mutation rate, including nuclear gene mutations affecting the fidelity and efficacy of mitochondrial DNA replication and/or repair and--consistent with the lack of RNA editing--exceptionally high levels of "mutagenic" retroprocessing. That the mitochondrial mutation rate has returned to normally low levels in many Geraniaceae raises the possibility that, akin to the ephemerality of mutator strains in bacteria, selection favors a low mutation rate in plant mitochondria.

Inheritance and recombination of mitochondrial genomes in plants, fungi and animals

It is generally assumed that mitochondrial genomes are uniparentally transmitted, homoplasmic and nonrecombining. However, these assumptions draw largely from early studies on animal mitochondrial DNA (mtDNA). In this review, we show that plants, animals and fungi are all characterized by episodes of biparental inheritance, recombination among genetically distinct partners, and selfish elements within the mitochondrial genome, but that the extent of these phenomena may vary substantially across taxa. We argue that occasional biparental mitochondrial transmission may allow organisms to achieve the best of both worlds by facilitating mutational clearance but continuing to restrict the spread of selfish genetic elements. We also show that methodological biases and disproportionately allocated study effort are likely to have influenced current estimates of the extent of biparental inheritance, heteroplasmy and recombination in mitochondrial genomes from different taxa. Despite these complications, there do seem to be discernible similarities and differences in transmission dynamics and likelihood of recombination of mtDNA in plant, animal and fungal taxa that should provide an excellent opportunity for comparative investigation of the evolution of mitochondrial genome dynamics.

Male-driven evolution in closely related species of the mouse genus Mus

Recently, other researchers have found that closely related primate species had a lower male-to-female mutation rate ratio (alpha) than distantly related species. To determine if this is a general phenomenon affecting other mammalian orders, eleven species or subspecies of the rodent genus Mus and two outgroup species were compared. Intron sequences from a gene in the nonrecombining region of the Y chromosome Jarid1d (Smcy) and its X chromosomal gametolog, Jarid1c (Smcx), were analyzed in a phylogenetic context. The male-to-female mutation rate ratio for all thirteen taxa is approximately 2.5, which is similar to previous estimates in more distantly related rodents. However, when branches with lengths of more than 2.5% were removed from the analysis, the male-to-female mutation rate ratio dropped to 0.9. Thus, in closely related rodents, as in closely related primates, the male-to-female mutation rate ratio is lower than expected.

Standardizing for microsatellite length in comparisons of genetic diversity

Mutation rates at microsatellites tend to increase with the number of repeats of the motif, leading to higher levels of polymorphism at long microsatellites. To standardize levels of diversity when microsatellites differ in size, we investigate the relationship between tract length and variation and provide a formula to adjust allelic richness to a fixed mean number of repeats in the specific case of chloroplast microsatellites. A comparison between 39 loci from eight species of conifers (where chloroplast DNA is paternally inherited) and 64 loci from 12 species of angiosperms (where chloroplast DNA is generally predominantly maternally inherited) indicates that the greater allelic richness found in conifers remains significant after controlling for number of repeats. The approach stresses the advantage of reporting variation in number of repeats instead of relative fragment sizes.

Rate variation in parasitic plants: correlated and uncorrelated patterns among plastid genes of different function

BACKGROUND

The analysis of synonymous and nonsynonymous rates of DNA change can help in the choice among competing explanations for rate variation, such as differences in constraint, mutation rate, or the strength of genetic drift. Nonphotosynthetic plants of the Orobanchaceae have increased rates of DNA change. In this study 38 taxa of Orobanchaceae and relatives were used and 3 plastid genes were sequenced for each taxon.

RESULTS

Phylogenetic reconstructions of relative rates of sequence evolution for three plastid genes (rbcL, matK and rps2) show significant rate heterogeneity among lineages and among genes. Many of the non-photosynthetic plants have increases in both synonymous and nonsynonymous rates, indicating that both (1) selection is relaxed, and (2) there has been a change in the rate at which mutations are entering the population in these species. However, rate increases are not always immediate upon loss of photosynthesis. Overall there is a poor correlation of synonymous and nonsynonymous rates. There is, however, a strong correlation of synonymous rates across the 3 genes studied and the lineage-speccific pattern for each gene is strikingly similar. This indicates that the causes of synonymous rate variation are affecting the whole plastid genome in a similar way. There is a weaker correlation across genes for nonsynonymous rates. Here the picture is more complex, as could be expected if there are many causes of variation, differing from taxon to taxon and gene to gene.

CONCLUSIONS

The distinctive pattern of rate increases in Orobanchaceae has at least two causes. It is clear that there is a relaxation of constraint in many (though not all) non-photosynthetic lineages. However, there is also some force affecting synonymous sites as well. At this point, it is not possible to tell whether it is generation time, speciation rate, mutation rate, DNA repair efficiency or some combination of these factors.

Mitochondrial substitution rates are extraordinarily elevated and variable in a genus of flowering plants

Plant mitochondrial (mt) genomes have long been known to evolve slowly in sequence. Here we show remarkable departure from this pattern of conservative evolution in a genus of flowering plants. Substitution rates at synonymous sites vary substantially among lineages within Plantago. At the extreme, rates in Plantago exceed those in exceptionally slow plant lineages by approximately 4,000-fold. The fastest Plantago lineages set a new benchmark for rapid evolution in a DNA genome, exceeding even the fastest animal mt genome by an order of magnitude. All six mt genes examined show similarly elevated divergence in Plantago, implying that substitution rates are highly accelerated throughout the genome. In contrast, substitution rates show little or no elevation in Plantago for each of four chloroplast and three nuclear genes examined. These results, combined with relatively modest elevations in rates of nonsynonymous substitutions in Plantago mt genes, indicate that major, reversible changes in the mt mutation rate probably underlie the extensive variation in synonymous substitution rates. These rate changes could be caused by major changes in any number of factors that control the mt mutation rate, from the production and detoxification of oxygen free radicals in the mitochondrion to the efficacy of mt DNA replication and/or repair.