Effective population size and tests of neutrality at cytoplasmic genes in Arabidopsis

A two-step self-pollination mechanism maximizes fertility in Brassicaceae

Self-pollination in self-compatible plant species often occurs prior to flower opening. By tracking the temporal progress of pollination in Arabidopsis, we observed that pollen predominantly targets the lateral region of the stigma in unopened flowers. Notably, approximately 7 h after flower opening, flowers close, thereby pressing anthers toward the central region of the stigma for a second self-pollination. This two-step self-pollination results in a doubling of pollen deposition, which significantly increases the ovule-targeting ratio and improves fertility under pollen-limiting conditions, as evident in the anther-dehiscence-defective mutant myb108 and under environmental stress conditions. Analysis using gamete-interaction-defective mutants hap2/gcs1 and dmp8 dmp9 revealed that the timely separation of both pollination events promotes fertilization recovery efficiency. A similar two-step pollination was observed in two other self-pollinating but not in outcrossing Brassicaceae species. This mechanism represents a reproductive assurance strategy in predominantly self-pollinating annuals to maximize fertility under unfavorable conditions.

Divergence in the Aquilegia ecalcarata complex is correlated with geography and climate oscillations: Evidence from plastid genome data

Quaternary climate oscillations and geographical heterogeneity play important roles in determining species and genetic diversity distribution patterns, but how these factors affect the migration and differentiation of East Asian plants species at the population level remains poorly understood. The Aquilegia ecalcarata complex, a group that originated in the Late Tertiary and is widely distributed throughout East Asia, displays high genetic variation that is suitable for studying elaborate phylogeographic patterns and demographic history related to the impact of Quaternary climate and geography. We used plastid genome data from 322 individuals in 60 populations of the A. ecalcarata complex to thoroughly explore the impact of Quaternary climate oscillations and geography on the phylogeographic patterns and demographic history of the A. ecalcarata complex through a series of phylogenetic, divergence time estimation, and demographic history analyses. The dry, cold climate and frequent climate oscillations that occurred during the early Pleistocene and the Mid-Pleistocene transition led to the differentiation of the A. ecalcarata complex, which was isolated in various areas. Geographically, the A. ecalcarata complex can be divided into Eastern and Western Clades and five subclades, which conform to the divergence of the East Asian flora. Our results clearly show the impact of Quaternary climate and geography on evolutionary history at the population level. These findings promote the understanding of the relationship between plant genetic differentiation and climate and geographical factors of East Asia at the population level.

Retrieving the deleterious mutations before extinction: genome-wide comparison of shared derived mutations in liver cancer and normal population

STUDY PURPOSE

Deleterious mutations would be rapidly purged from natural populations along with the extinction of their carriers. The currently observed mutations in existing species are mostly neutral. The inaccessibility of deleterious mutations impedes the functional studies on how these mutations affect the fitness at individual level.

STUDY DESIGN

The connection between the deleterious genotype and the non-adaptive phenotype could be bridged by sequencing the genome before extinction. Although this approach is no longer feasible for evolutionary biologists, it is feasible for cancer biologists by profiling the mutations in tumour samples which are so deleterious that the carriers hardly live.

RESULTS

By comparing the derived mutation profile between normal populations and patients with liver cancer, we found that the shared mutations, which are highly deleterious, are suppressed to low allele frequencies in normal populations and tissues, but show remarkably high frequency in tumours. The density of shared mutations is negatively correlated with gene conservation and expression levels.

CONCLUSIONS

Deleterious mutations are suppressed in functionally important genes as well as in normal populations. This work deepened our understanding on how natural selection act on deleterious mutations by analogising the cancer evolution to species evolution, which are essentially the same molecular process but at different time scales.

An evaluation of alternative explanations for widespread cytonuclear discordance in annual sunflowers (Helianthus)

Cytonuclear discordance is commonly observed in phylogenetic studies, yet few studies have tested whether these patterns reflect incomplete lineage sorting or organellar introgression. Here, we used whole-chloroplast sequence data in combination with over 1000 nuclear single-nucleotide polymorphisms to clarify the extent of cytonuclear discordance in wild annual sunflowers (Helianthus), and to test alternative explanations for such discordance. Our phylogenetic analyses indicate that cytonuclear discordance is widespread within this group, both in terms of the relationships among species and among individuals within species. Simulations of chloroplast evolution show that incomplete lineage sorting cannot explain these patterns in most cases. Instead, most of the observed discordance is better explained by cytoplasmic introgression. Molecular tests of evolution further indicate that selection may have played a role in driving patterns of plastid variation - although additional experimental work is needed to fully evaluate the importance of selection on organellar variants in different parts of the geographic range. Overall, this study represents one of the most comprehensive tests of the drivers of cytonuclear discordance and highlights the potential for gene flow to lead to extensive organellar introgression in hybridizing taxa.

Patterns of Genome-Wide Nucleotide Diversity in the Gynodioecious Plant Thymus vulgaris Are Compatible with Recent Sweeps of Cytoplasmic Genes

Gynodioecy is a sexual dimorphism where females coexist with hermaphrodite individuals. In most cases, this dimorphism involves the interaction of cytoplasmic male sterility (CMS) genes and nuclear restorer genes. Two scenarios can account for how these interactions maintain gynodioecy. Either CMS genes recurrently enter populations at low frequency via mutation or migration and go to fixation unimpeded (successive sweeps), or CMS genes maintain polymorphism over evolutionary time through interactions with a nuclear restorer allele (balanced polymorphism). To distinguish between these scenarios, we used transcriptome sequencing in gynodioecious Thymus vulgaris and surveyed genome-wide diversity in 18 naturally occurring individuals sampled from populations at a local geographic scale. We contrast the amount and patterns of nucleotide diversity in the nuclear and cytoplasmic genome, and find ample diversity at the nuclear level (π = 0.019 at synonymous sites) but reduced genetic diversity and an excess of rare polymorphisms in the cytoplasmic genome relative to the nuclear genome. Our finding is incompatible with the maintenance of gynodioecy via scenarios invoking long-term balancing selection, and instead suggests the recent fixation of CMS lineages in the populations studied.

Radical amino acid mutations persist longer in the absence of sex

Harmful mutations are ubiquitous and inevitable, and the rate at which these mutations are removed from populations is a critical determinant of evolutionary fate. Closely related sexual and asexual taxa provide a particularly powerful setting to study deleterious mutation elimination because sexual reproduction should facilitate mutational clearance by reducing selective interference between sites and by allowing the production of offspring with different mutational complements than their parents. Here, we compared the rate of removal of conservative (i.e., similar biochemical properties) and radical (i.e., distinct biochemical properties) nonsynonymous mutations from mitochondrial genomes of sexual versus asexual Potamopyrgus antipodarum, a New Zealand freshwater snail characterized by coexisting and ecologically similar sexual and asexual lineages. Our analyses revealed that radical nonsynonymous mutations are cleared at higher rates than conservative changes and that sexual lineages eliminate radical changes more rapidly than asexual counterparts. These results are consistent with reduced efficacy of purifying selection in asexual lineages allowing harmful mutations to remain polymorphic longer than in sexual lineages. Together, these data illuminate some of the population-level processes contributing to mitochondrial mutation accumulation and suggest that mutation accumulation could influence the outcome of competition between sexual and asexual lineages.

Rapid genome shrinkage in a self-fertile nematode reveals sperm competition proteins

To reveal impacts of sexual mode on genome content, we compared chromosome-scale assemblies of the outcrossing nematode Caenorhabditis nigoni to its self-fertile sibling species, C. briggsaeC. nigoni's genome resembles that of outcrossing relatives but encodes 31% more protein-coding genes than C. briggsaeC. nigoni genes lacking C. briggsae orthologs were disproportionately small and male-biased in expression. These include the male secreted short (mss) gene family, which encodes sperm surface glycoproteins conserved only in outcrossing species. Sperm from mss-null males of outcrossing C. remanei failed to compete with wild-type sperm, despite normal fertility in noncompetitive mating. Restoring mss to C. briggsae males was sufficient to enhance sperm competitiveness. Thus, sex has a pervasive influence on genome content that can be used to identify sperm competition factors.

Recombination-dependent replication and gene conversion homogenize repeat sequences and diversify plastid genome structure

PREMISE OF THE STUDY

There is a misinterpretation in the literature regarding the variable orientation of the small single copy region of plastid genomes (plastomes). The common phenomenon of small and large single copy inversion, hypothesized to occur through intramolecular recombination between inverted repeats (IR) in a circular, single unit-genome, in fact, more likely occurs through recombination-dependent replication (RDR) of linear plastome templates. If RDR can be primed through both intra- and intermolecular recombination, then this mechanism could not only create inversion isomers of so-called single copy regions, but also an array of alternative sequence arrangements.

METHODS

We used Illumina paired-end and PacBio single-molecule real-time (SMRT) sequences to characterize repeat structure in the plastome of Monsonia emarginata (Geraniaceae). We used OrgConv and inspected nucleotide alignments to infer ancestral nucleotides and identify gene conversion among repeats and mapped long (>1 kb) SMRT reads against the unit-genome assembly to identify alternative sequence arrangements.

RESULTS

Although M. emarginata lacks the canonical IR, we found that large repeats (>1 kilobase; kb) represent ∼22% of the plastome nucleotide content. Among the largest repeats (>2 kb), we identified GC-biased gene conversion and mapping filtered, long SMRT reads to the M. emarginata unit-genome assembly revealed alternative, substoichiometric sequence arrangements.

CONCLUSION

We offer a model based on RDR and gene conversion between long repeated sequences in the M. emarginata plastome and provide support that both intra-and intermolecular recombination between large repeats, particularly in repeat-rich plastomes, varies unit-genome structure while homogenizing the nucleotide sequence of repeats.

On the adaptive value of cytoplasmic genomes in plants

Is DNA variation maintained in organelle genomes selectively neutral? The answer to this question has important implications for many aspects of ecology and evolution. While traditionally the answer has been 'yes', recent studies in animals have shown that, on the contrary, mitochondrial DNA polymorphism is frequently adaptive. In plants, however, the neutrality assumption has not been strongly challenged. Here, we begin with a critical evaluation of arguments in favour of this long-held view. We then discuss the latest empirical evidence for the opposing prediction that sequence variation in plant cytoplasmic genomes is frequently adaptive. While outstanding research progress is being made towards understanding this fundamental topic, we highlight the need for studies that combine information ranging from field experiments to physiology to molecular evolutionary biology. Such an interdisciplinary approach provides a means for determining the frequency, drivers and evolutionary significance of adaptive organelle DNA variation.

Mating systems and selection efficacy: a test using chloroplastic sequence data in Angiosperms

Selfing is assumed to reduce selection efficacy, especially purifying selection. This can be tested using molecular data, for example by comparing the Dn/Ds ratio between selfing and outcrossing lineages. So far, little evidence of relaxed selection against weakly deleterious mutations (as inferred by a higher Dn/Ds ratio) in selfers as compared to outcrossers has been found, contrary to the pattern often observed between asexual and sexual lineages. However, few groups have been studied to date. To further test this hypothesis, we compiled and analysed chloroplastic sequence data sets in several plant groups. We found a general trend towards relaxed selection in selfers in our data sets but with weak statistical support. Simulations suggested that the results were compatible with weak-to-moderate Dn/Ds ratio differences in selfing lineages. Simple theoretical predictions also showed that the ability to detect relaxed selection in selfers could strongly depend on the distribution of the effects of deleterious mutations on fitness. Our results are compatible with a recent origin of selfing lineages whereby deleterious mutations potentially have a strong impact on population extinction or with a more ancient origin but without a marked effect of deleterious mutations on the extinction dynamics.

A high-resolution genetic map of yellow monkeyflower identifies chemical defense QTLs and recombination rate variation

Genotyping-by-sequencing methods have vastly improved the resolution and accuracy of genetic linkage maps by increasing both the number of marker loci as well as the number of individuals genotyped at these loci. Using restriction-associated DNA sequencing, we construct a dense linkage map for a panel of recombinant inbred lines derived from a cross between divergent ecotypes of Mimulus guttatus. We used this map to estimate recombination rate across the genome and to identify quantitative trait loci for the production of several secondary compounds (PPGs) of the phenylpropanoid pathway implicated in defense against herbivores. Levels of different PPGs are correlated across recombinant inbred lines suggesting joint regulation of the phenylpropanoid pathway. However, the three quantitative trait loci identified in this study each act on a distinct PPG. Finally, we map three putative genomic inversions differentiating the two parental populations, including a previously characterized inversion that contributes to life-history differences between the annual/perennial ecotypes.

Seed-mediated connectivity among fragmented populations of Quercus castanea (Fagaceae) in a Mexican landscape

PREMISE OF STUDY

Anthropogenic fragmentation is an ongoing process in many forested areas that may create loss of connectivity among tree populations and constitutes a serious threat to ecological and genetic processes. We tested the central hypothesis that seed dispersal mitigates the impact of fragmentation by comparing connectivity and genetic diversity of adult vs. seedling populations in recently fragmented populations of the Mexican red oak Quercus castanea.

METHODS

Adult individuals, established before fragmentation, and seedlings, established after fragmentation, were sampled at 33 forest fragments of variable size (0.2 to 294 ha) within the Cuitzeo basin, Michoacán state, and genotyped using seven highly polymorphic chloroplast microsatellite markers (cpSSRs). To test whether seed dispersal retains connectivity among fragmented populations, we compared genetic diversity and connectivity networks between adults and progeny and determined the effect of fragment size on these values.

KEY RESULTS

Seventy haplotypes were identified, 63 in the adults and 60 in the seedlings, with average within-population diversity (hS) values of 0.624 in the adults and 0.630 in the seedlings. A positive correlation of genetic diversity values with fragment size was found in the seedling populations but not in the adult populations. The network connectivity analysis revealed lower connectivity among seedling populations than among adults. The number of connections (edges) as well as other network properties, such as betweenness centrality, node degree and closeness, were significantly lower in the seedlings network.

CONCLUSIONS

Habitat fragmentation in this landscape is disrupting seed-dispersal-mediated genetic connectivity among extant populations.

How closely does genetic diversity in finite populations conform to predictions of neutral theory? Large deficits in regions of low recombination

Levels of genetic diversity in finite populations are crucial in conservation and evolutionary biology. Genetic diversity is required for populations to evolve and its loss is related to inbreeding in random mating populations, and thus to reduced population fitness and increased extinction risk. Neutral theory is widely used to predict levels of genetic diversity. I review levels of genetic diversity in finite populations in relation to predictions of neutral theory. Positive associations between genetic diversity and population size, as predicted by neutral theory, are observed for microsatellites, allozymes, quantitative genetic variation and usually for mitochondrial DNA (mtDNA). However, there are frequently significant deviations from neutral theory owing to indirect selection at linked loci caused by balancing selection, selective sweeps and background selection. Substantially lower genetic diversity than predicted under neutrality was found for chromosomes with low recombination rates and high linkage disequilibrium (compared with 'normally' recombining chromosomes within species and adjusted for different copy numbers and mutation rates), including W (median 100% lower) and Y (89% lower) chromosomes, dot fourth chromosomes in Drosophila (94% lower) and mtDNA (67% lower). Further, microsatellite genetic and allelic diversity were lost at 12 and 33% faster rates than expected in populations adapting to captivity, owing to widespread selective sweeps. Overall, neither neutral theory nor most versions of the genetic draft hypothesis are compatible with all empirical results.

Twenty-fold difference in evolutionary rates between the mitochondrial and plastid genomes of species with secondary red plastids

Within plastid-bearing species, the relative rates of evolution between mitochondrial and plastid genomes are poorly studied, but for the few lineages in which they have been explored, including land plants and green algae, the mitochondrial DNA mutation rate is nearly always estimated to be lower than or equal to that of the plastid DNA. Here, we show that in protists from three distinct lineages with secondary, red algal-derived plastids, the opposite is true: their mitochondrial genomes are evolving 5-30 times faster than their plastid genomes, even when the plastid is nonphotosynthetic. These findings have implications for understanding the origins and evolution of organelle genome architecture and the genes they encode.

Genetic consequences of seed banks in the perennial herb Arabidopsis lyrata subsp. petraea (Brassicaceae)

PREMISE OF THE STUDY

Seed banks may increase the effective population size (N(e)) of plants as a result of elevated coalescence times for alleles residing in the populations. This has been empirically demonstrated in populations of the annual Arabidopsis thaliana, whereas comparable data for perennial species are currently lacking. We studied the contribution of seed banks to effective sizes of natural populations of the self-incompatible, perennial Arabidopsis lyrata subsp. petraea, a close relative of A. thaliana. •

METHODS

Fourteen populations of A. lyrata collected throughout the Norwegian distribution range were analyzed using microsatellite markers. •

KEY RESULTS

The genetic composition of seed-bank and aboveground cohorts was found to be highly similar, with little genetic differentiation between cohorts in most populations. However, the proportion of private alleles was higher in aboveground than in seed-bank cohorts. The presence of seed banks significantly increased total N(e), but the contribution from seed banks to overall N(e) were lower than the contribution from aboveground cohorts in most populations. Estimated historical N(e) values, reflecting the effective sizes of populations throughout the history of the species, were considerably higher than estimates of contemporary N(e), reflecting number of reproducing individuals within the past few generations. •

CONCLUSIONS

The results show that the seed bank contributes to total N(e) in the perennial herb A. lyrata. However, the contribution is similar to or lower than that of the above-ground fraction of the population and markedly weaker than that previously documented in the annual A. thaliana.

Correlation between nuclear plastid DNA abundance and plastid number supports the limited transfer window hypothesis

The abundance of nuclear plastid DNA-like sequences (NUPTs) in nuclear genomes can vary immensely; however, the forces responsible for this variation are poorly understood. “The limited transfer window hypothesis” predicts that species with only one plastid per cell will have fewer NUPTs than those with many plastids per cell, but a lack of genome sequence data from monoplastidic species has made this hypothesis difficult to test. Here, by analyzing newly available genome sequences from diverse mono- and polyplastidic taxa, we show that the hypothesis holds. On average, the polyplastidic species we studied had 80 times more NUPTs than those that were monoplastidic. Moreover, NUPT content was positively related to nuclear genome size, indicating that in addition to plastid number, NUPTs are influenced by the forces controlling the expansion and contraction of noncoding nuclear DNA. These findings are consistent with data on nuclear DNAs of mitochondrial origin (NUMTs), suggesting that similar processes govern the abundance of both NUPTs and NUMTs.

Sex-linked inheritance in macaque monkeys: implications for effective population size and dispersal to Sulawesi

Sex-specific differences in dispersal, survival, reproductive success, and natural selection differentially affect the effective population size (N(e)) of genomic regions with different modes of inheritance such as sex chromosomes and mitochondrial DNA. In papionin monkeys (macaques, baboons, geladas, mandrills, drills, and mangabeys), for example, these factors are expected to reduce N(e) of paternally inherited portions of the genome compared to maternally inherited portions. To explore this further, we quantified relative N(e) of autosomal DNA, X and Y chromosomes, and mitochondrial DNA using molecular polymorphism and divergence information from pigtail macaque monkeys (Macaca nemestrina). Consistent with demographic expectations, we found that N(e) of the Y is lower than expected from a Wright-Fisher idealized population with an equal proportion of males and females, whereas N(e) of mitochondrial DNA is higher. However, N(e) of 11 loci on the X chromosome was lower than expected, a finding that could be explained by pervasive hitchhiking effects on this chromosome. We evaluated the fit of these data to various models involving natural selection or sex-biased demography. Significant support was recovered for natural selection acting on the Y chromosome. A demographic model with a skewed sex ratio was more likely than one with sex-biased migration and explained the data about as well as an ideal model without sex-biased demography. We then incorporated these results into an evaluation of macaque divergence and migration on Borneo and Sulawesi islands. One X-linked locus was not monophyletic on Sulawesi, but multilocus data analyzed in a coalescent framework failed to reject a model without migration between these islands after both were colonized.

Low nucleotide diversity for the expanded organelle and nuclear genomes of Volvox carteri supports the mutational-hazard hypothesis

The noncoding-DNA content of organelle and nuclear genomes can vary immensely. Both adaptive and nonadaptive explanations for this variation have been proposed. This study addresses a nonadaptive explanation called the mutational-hazard hypothesis and applies it to the mitochondrial, plastid, and nuclear genomes of the multicellular green alga Volvox carteri. Given the expanded architecture of the V. carteri organelle and nuclear genomes (60-85% noncoding DNA), the mutational-hazard hypothesis would predict them to have less silent-site nucleotide diversity (π(silent)) than their more compact counterparts from other eukaryotes-ultimately reflecting differences in 2N(g)μ (twice the effective number of genes per locus in the population times the mutation rate). The data presented here support this prediction: Analyses of mitochondrial, plastid, and nuclear DNAs from seven V. carteri forma nagariensis geographical isolates reveal low values of π(silent) (0.00038, 0.00065, and 0.00528, respectively), much lower values than those previously observed for the more compact organelle and nuclear DNAs of Chlamydomonas reinhardtii (a close relative of V. carteri). We conclude that the large noncoding-DNA content of the V. carteri genomes is best explained by the mutational-hazard hypothesis and speculate that the shift from unicellular to multicellular life in the ancestor that gave rise to V. carteri contributed to a low V. carteri population size and thus a reduced 2N(g)μ. Complete mitochondrial and plastid genome maps for V. carteri are also presented and compared with those of C. reinhardtii.

Signature of diversifying selection on members of the pentatricopeptide repeat protein family in Arabidopsis lyrata

Pentatricopeptide repeat (PPR) proteins compose a family of nuclear-encoded transcriptional regulators of cytoplasmic genes. They have shown dramatic expansion in copy number in plants, and although the functional importance of many remains unclear, a subset has been repeatedly implicated as nuclear restorers for cytoplasmic male sterility. Here we investigate the molecular population genetics and molecular evolution of seven single-copy PPR genes in the outcrossing model plant Arabidopsis lyrata. In comparison with neutral reference loci, we find, on average, elevated levels of polymorphism and an excess of high-frequency variants at these PPR genes, suggesting that natural selection is maintaining polymorphism at some of these loci. This elevation in diversity persists when we control for divergence and generally decreases in the flanking regions, suggesting that these genes are themselves the targets of selection. Some of the PPR genes also demonstrate elevated population differentiation, which is consistent with spatially varying selection. In contrast, no comparable patterns are observed at these loci in A. thaliana, providing no evidence for the action of balancing selection in this selfing species. Taken together, these results suggest that a subset of PPR genes may be subject to balancing selection associated with ongoing cytonuclear coevolution in the outcrossing A. lyrata, which is possibly mediated either by intergenomic conflict or by compensatory evolution.

Nucleotide diversity of the Chlamydomonas reinhardtii plastid genome: addressing the mutational-hazard hypothesis

BACKGROUND

The mutational-hazard hypothesis argues that the noncoding-DNA content of a genome is a consequence of the mutation rate (mu) and the effective number of genes per locus in the population (N(g)). The hypothesis predicts that genomes with a high N(g)mu will be more compact than those with a small N(g)mu. Approximations of N(g)mu can be gained by measuring the nucleotide diversity at silent sites (pi(silent)). We addressed the mutation-hazard hypothesis apropos plastid-genome evolution by measuring pi(silent) of the Chlamydomonas reinhardtii plastid DNA (ptDNA), the most noncoding-DNA-dense plastid genome observed to date. The data presented here in conjunction with previously published values of pi(silent) for the C. reinhardtii mitochondrial and nuclear genomes, which are respectively compact and bloated, allow for a complete analysis of nucleotide diversity and genome compactness in all three genetic compartments of this model organism.

RESULTS

In C. reinhardtii, the mean estimate of pi(silent) for the ptDNA (14.5 x 10(-3)) is less than that of the nuclear DNA (32 x 10(-3)) and greater than that of the mitochondrial DNA (8.5 x 10(-3)). On average, C. reinhardtii has approximately 4 times more silent-site ptDNA diversity than the mean value reported for land plants, which have more compact plastid genomes. The silent-site nucleotide diversity of the different ptDNA loci that were studied varied significantly: from 0 to 71 x 10(-3) for synonymous sites and from 0 to 42 x 10(-3) for intergenic regions.

CONCLUSION

Our findings on silent-site ptDNA diversity are inconsistent with what would be expected under the mutational-hazard hypothesis and go against the documented trend in other systems of pi(silent) positively correlating with genome compactness. Overall, we highlight the lack of reliable nucleotide-diversity measurements for ptDNA and hope that the values presented here will act as sound data for future research concerning the mutational-hazard hypothesis and plastid evolution in general.

Evolution of the Caenorhabditis elegans genome

A fundamental problem in genome biology is to elucidate the evolutionary forces responsible for generating nonrandom patterns of genome organization. As the first metazoan to benefit from full-genome sequencing, Caenorhabditis elegans has been at the forefront of research in this area. Studies of genomic patterns, and their evolutionary underpinnings, continue to be augmented by the recent push to obtain additional full-genome sequences of related Caenorhabditis taxa. In the near future, we expect to see major advances with the onset of whole-genome resequencing of multiple wild individuals of the same species. In this review, we synthesize many of the important insights to date in our understanding of genome organization and function that derive from the evolutionary principles made explicit by theoretical population genetics and molecular evolution and highlight fertile areas for future research on unanswered questions in C. elegans genome evolution. We call attention to the need for C. elegans researchers to generate and critically assess nonadaptive hypotheses for genomic and developmental patterns, in addition to adaptive scenarios. We also emphasize the potential importance of evolution in the gonochoristic (female and male) ancestors of the androdioecious (hermaphrodite and male) C. elegans as the source for many of its genomic and developmental patterns.

The evolution of plant genomes: scaling up from a population perspective

Plant genomes exhibit tremendous diversity in both their size and structure, with genome sizes across land plants ranging over two to three orders of magnitude and significant variation in structural organization was observed across species (EA Kellogg, JL Bennetzen, The evolution of nuclear genome structure in seed plants, Am J Bot 2004, 91:1709-1725). Five plant genomes are now either completely sequenced or in the draft stage; the grape (O Jaillon et al., The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla, Nature 2007, 449:463-467) and papaya (R Ming et al., The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus), Nature 2008, 452:991-997) whole genome sequences were reported most recently. Moreover, sequencing of 41 additional genomes is in progress. There is now an emerging consensus that understanding genome evolution requires consideration of the population genetics of genome diversification, and that description of evolutionary forces at the level of populations and within species can help identify the features that led to plant genome diversity (M Lynch, JS Conery, The origins of genome complexity, Science 2003, 302:1401-1404). In this review we focus on advances in our understanding of the mechanisms that drive the diversification of genomes. In particular, we look at the extent to which demographic features such as effective population size changes within species can drive genome evolution, discuss population genetic models of genome diversification associated with transposable element (TE) mobilization, and describe recent studies on the evolution of gene families.

On the potential for extinction by Muller's ratchet in Caenorhabditis elegans

Background

The self-fertile hermaphrodite worm C. elegans is an important model organism for biology, yet little is known about the origin and persistence of the self-fertilizing mode of reproduction in this lineage. Recent work has demonstrated an extraordinary degree of selfing combined with a high deleterious mutation rate in contemporary populations. These observations raise the question as to whether the mutation load might rise to such a degree as to eventually threaten the species with extinction. The potential for such a process to occur would inform our understanding of the time since the origin of self-fertilization in C. elegans history.

Results

To address this issue, here we quantify the rate of fitness decline expected to occur via Muller's ratchet for a purely selfing population, using both analytical approximations and globally distributed individual-based simulations from the evolution@home system to compute the rate of deleterious mutation accumulation. Using the best available estimates for parameters of how C. elegans evolves, we conclude that pure selfing can persist for only short evolutionary intervals, and is expected to lead to extinction within thousands of years for a plausible portion of parameter space. Credible lower-bound estimates of nuclear mutation rates do not extend the expected time to extinction much beyond a million years.

Conclusion

Thus we conclude that either the extreme self-fertilization implied by current patterns of genetic variation in C. elegans arose relatively recently or that low levels of outcrossing and other factors are key to the persistence of C. elegans into the present day. We also discuss results for the mitochondrial genome and the implications for C. briggsae, a close relative that made the transition to selfing independently of C. elegans.