Mating system shifts and transposable element evolution in the plant genus Capsella

Background

Despite having predominately deleterious fitness effects, transposable elements (TEs) are major constituents of eukaryote genomes in general and of plant genomes in particular. Although the proportion of the genome made up of TEs varies at least four-fold across plants, the relative importance of the evolutionary forces shaping variation in TE abundance and distributions across taxa remains unclear. Under several theoretical models, mating system plays an important role in governing the evolutionary dynamics of TEs. Here, we use the recently sequenced Capsella rubella reference genome and short-read whole genome sequencing of multiple individuals to quantify abundance, genome distributions, and population frequencies of TEs in three recently diverged species of differing mating system, two self-compatible species (C. rubella and C. orientalis) and their self-incompatible outcrossing relative, C. grandiflora.

Results

We detect different dynamics of TE evolution in our two self-compatible species; C. rubella shows a small increase in transposon copy number, while C. orientalis shows a substantial decrease relative to C. grandiflora. The direction of this change in copy number is genome wide and consistent across transposon classes. For insertions near genes, however, we detect the highest abundances in C. grandiflora. Finally, we also find differences in the population frequency distributions across the three species.

Conclusion

Overall, our results suggest that the evolution of selfing may have different effects on TE evolution on a short and on a long timescale. Moreover, cross-species comparisons of transposon abundance are sensitive to reference genome bias, and efforts to control for this bias are key when making comparisons across species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-602) contains supplementary material, which is available to authorized users.     

Genomic Identification of Founding Haplotypes Reveals the History of the Selfing Species Capsella rubella

The shift from outcrossing to self-fertilization is among the most common evolutionary transitions in flowering plants. Until recently, however, a genome-wide view of this transition has been obscured by both a dearth of appropriate data and the lack of appropriate population genomic methods to interpret such data. Here, we present a novel population genomic analysis detailing the origin of the selfing species, Capsella rubella, which recently split from its outcrossing sister, Capsella grandiflora. Due to the recency of the split, much of the variation within C. rubella is also found within C. grandiflora. We can therefore identify genomic regions where two C. rubella individuals have inherited the same or different segments of ancestral diversity (i.e. founding haplotypes) present in C. rubella''s founder(s). Based on this analysis, we show that C. rubella was founded by multiple individuals drawn from a diverse ancestral population closely related to extant C. grandiflora, that drift and selection have rapidly homogenized most of this ancestral variation since C. rubella''s founding, and that little novel variation has accumulated within this time. Despite the extensive loss of ancestral variation, the approximately 25% of the genome for which two C. rubella individuals have inherited different founding haplotypes makes up roughly 90% of the genetic variation between them. To extend these findings, we develop a coalescent model that utilizes the inferred frequency of founding haplotypes and variation within founding haplotypes to estimate that C. rubella was founded by a potentially large number of individuals between 50 and 100 kya, and has subsequently experienced a twenty-fold reduction in its effective population size. As population genomic data from an increasing number of outcrossing/selfing pairs are generated, analyses like the one developed here will facilitate a fine-scaled view of the evolutionary and demographic impact of the transition to self-fertilization.     

Bypassing reproductive barriers in hybrid seeds using chemically induced epimutagenesis

The triploid block, which prevents interploidy hybridizations in flowering plants, is characterized by a failure in endosperm development, arrest in embryogenesis, and seed collapse. Many genetic components of triploid seed lethality have been successfully identified in the model plant Arabidopsis thaliana, most notably the paternally expressed genes (PEGs), which are upregulated in tetraploid endosperm with paternal excess. Previous studies have shown that the paternal epigenome is a key determinant of the triploid block response, as the loss of DNA methylation in diploid pollen suppresses the triploid block almost completely. Here, we demonstrate that triploid seed collapse is bypassed in Arabidopsis plants treated with the DNA methyltransferase inhibitor 5-Azacytidine during seed germination and early growth. We identified strong suppressor lines showing stable transgenerational inheritance of hypomethylation in the CG context, as well as normalized expression of PEGs in triploid seeds. Importantly, differentially methylated loci segregate in the progeny of “epimutagenized” plants, which may allow epialleles involved in the triploid block response to be identified in future studies. Finally, we demonstrate that chemically induced epimutagenesis facilitates hybridization between different Capsella species, thus potentially emerging as a strategy for producing triploids and interspecific hybrids with high agronomic interest.

Genome-wide loss of DNA methylation induced by 5-Azacytidine allows interploidy and interspecific hybridization barriers to be bypassed in Arabidopsis and Capsella.     

Independent FLC Mutations as Causes of Flowering-Time Variation in Arabidopsis thaliana and Capsella rubella

Capsella rubella is an inbreeding annual forb closely related to Arabidopsis thaliana, a model species widely used for studying natural variation in adaptive traits such as flowering time. Although mutations in dozens of genes can affect flowering of A. thaliana in the laboratory, only a handful of such genes vary in natural populations. Chief among these are FRIGIDA (FRI) and FLOWERING LOCUS C (FLC). Common and rare FRI mutations along with rare FLC mutations explain a large fraction of flowering-time variation in A. thaliana. Here we document flowering time under different conditions in 20 C. rubella accessions from across the species’ range. Similar to A. thaliana, vernalization, long photoperiods and elevated ambient temperature generally promote flowering. In this collection of C. rubella accessions, we did not find any obvious loss-of-function FRI alleles. Using mapping-by-sequencing with two strains that have contrasting flowering behaviors, we identified a splice-site mutation in FLC as the likely cause of early flowering in accession 1408. However, other similarly early C. rubella accessions did not share this mutation. We conclude that the genetic basis of flowering-time variation in C. rubella is complex, despite this very young species having undergone an extreme genetic bottleneck when it split from C. grandiflora a few tens of thousands of years ago.     

The Genomic Aftermath of Hybridization in the Opportunistic Pathogen Candida metapsilosis

Candida metapsilosis is a rarely-isolated, opportunistic pathogen that belongs to a clade of pathogenic yeasts known as the C. parapsilosis sensu lato species complex. To gain insight into the recent evolution of C. metapsilosis and the genetic basis of its virulence, we sequenced the genome of 11 clinical isolates from various locations, which we compared to each other and to the available genomes of the two remaining members of the complex: C. orthopsilosis and C. parapsilosis. Unexpectedly, we found compelling genomic evidence that C. metapsilosis is a highly heterozygous hybrid species, with all sequenced clinical strains resulting from the same past hybridization event involving two parental lineages that were approximately 4.5% divergent in sequence. This result indicates that the parental species are non-pathogenic, but that hybridization between them formed a new opportunistic pathogen, C. metapsilosis, that has achieved a worldwide distribution. We show that these hybrids are diploid and we identified strains carrying loci for both alternative mating types, which supports mating as the initial mechanism for hybrid formation. We trace the aftermath of this hybridization at the genomic level, and reconstruct the evolutionary relationships among the different strains. Recombination and introgression -resulting in loss of heterozygosis- between the two subgenomes have been rampant, and includes the partial overwriting of the MTLa mating locus in all strains. Collectively, our results shed light on the recent genomic evolution within the C. parapsilosis sensu lato complex, and argue for a re-definition of species within this clade, with at least five distinct homozygous lineages, some of which having the ability to form hybrids.     

Pleistocene dynamics of the Eurasian steppe as a driving force of evolution: Phylogenetic history of the genus Capsella (Brassicaceae)

Capsella is a model plant genus of the Brassicaceae closely related to Arabidopsis. To disentangle its biogeographical history and intrageneric phylogenetic relationships, 282 individuals of all five currently recognized Capsella species were genotyped using a restriction digest‐based next‐generation sequencing method. Our analysis retrieved two main lineages within Capsella that split c. one million years ago, with western C. grandiflora and C. rubella forming a sister lineage to the eastern lineage consisting of C. orientalis. The split was attributed to continuous latitudinal displacements of the Eurasian steppe belt to the south during Early Pleistocene glacial cycles. During the interglacial cycles of the Late Pleistocene, hybridization of the two lineages took place in the southwestern East European Plain, leading to the allotetraploid C. bursa‐pastoris. Extant genetic variation within C. orientalis postdated any extensive glacial events. Ecological niche modeling showed that suitable habitat for C. orientalis existed during the Last Glacial Maximum around the north coast of the Black Sea and in southern Kazakhstan. Such a scenario is also supported by population genomic data that uncovered the highest genetic diversity in the south Kazakhstan cluster, suggesting that C. orientalis originated in continental Asia and migrated north‐ and possibly eastwards after the last ice age. Post‐glacial hybridization events between C. bursa‐pastoris and C. grandiflora/rubella in the southwestern East European Plain and the Mediterranean gave rise to C. thracica. Introgression of C. grandiflora/rubella into C. bursa‐pastoris resulted in a new Mediterranean cluster within the already existing Eurasian C. bursa‐pastoris cluster. This study shows that the continuous displacement and disruption of the Eurasian steppe belt during the Pleistocene was the driving force in the evolution of Capsella.     

Increased Maternal Genome Dosage Bypasses the Requirement of the FIS Polycomb Repressive Complex 2 in Arabidopsis Seed Development

Seed development in flowering plants is initiated after a double fertilization event with two sperm cells fertilizing two female gametes, the egg cell and the central cell, leading to the formation of embryo and endosperm, respectively. In most species the endosperm is a polyploid tissue inheriting two maternal genomes and one paternal genome. As a consequence of this particular genomic configuration the endosperm is a dosage sensitive tissue, and changes in the ratio of maternal to paternal contributions strongly impact on endosperm development. The FERTILIZATION INDEPENDENT SEED (FIS) Polycomb Repressive Complex 2 (PRC2) is essential for endosperm development; however, the underlying forces that led to the evolution of the FIS-PRC2 remained unknown. Here, we show that the functional requirement of the FIS-PRC2 can be bypassed by increasing the ratio of maternal to paternal genomes in the endosperm, suggesting that the main functional requirement of the FIS-PRC2 is to balance parental genome contributions and to reduce genetic conflict. We furthermore reveal that the AGAMOUS LIKE (AGL) gene AGL62 acts as a dosage-sensitive seed size regulator and that reduced expression of AGL62 might be responsible for reduced size of seeds with increased maternal genome dosage.     

Flower morphology and pollen germination in the genus Capsella (Brassicaceae)

Evolutive changes in mating systems are often accompanied by changes in flower morphology, such as the reduction in size or even loss of petals, changes in production of volatiles, pollen/ovule ratio, the position between anthers and stigma and the germination time of pollen after pollination. These changes have been merged under the term “selfing syndrome” and often result in new taxonomic species. The evolutionary shift frequently happens parallel within many families and genera, for example within the Brassicaceae family. Within the genus Capsella, which is closely related to the molecular model species pair Arabidopsis lyrata (SI)/A. thaliana (SC), we studied self-incompatible and self-compatible species. SC species C. rubella and C. bursa-pastoris produce in comparison with the SI species C. grandiflora (i) smaller petals as the result of decreased cell division and only less of decreasing cell volume, (ii) less production of pollen in one flower, (iii) show a lesser incision between the two valves of the fruits, in combination with a shorter style, and (iv) have a much quicker fertilization of SC pollen after pollination. Crossing success between the diploid species, between different provenances of the tetraploid C. bursa-pastoris, and between the two diploid species and particular individuals of the self-incompatible C. grandiflora has been proven.     

Aspartate aminotransferase isozymes in the genusCapsella (Brassicaceae): Subcellular location,gene duplication,and polymorphism

The subcellular location of aspartate aminotransferase isozymes (EC 2.6.1.1) in the genusCapsella(Brassicaceae) was studied. The diploid speciesC. grandiflora andC. rubella have three AAT isozymes, including one located in the plastids. Each locus is duplicated in the tetraploidCapsella bursa-pastoris. Variation at the plastid-coding locus exceeded that at the other loci.C. bursa-pastoris had some unique alleles not detected in the diploid species. Segregation in open-pollinated families revealed thatCapsella grandiflora was outcrossing, whereasC. rubella was highly inbred, with most populations homozygous or uniform at all three loci. Inheritance in the tetraploid colonizerC. bursa-pastoris is disomic. This species was also predominantly selfing with outcrossing rates between 2% and 10%.Financial support by the German Research Foundation DFG is gratefully acknowledged.     

Multiple Origins of the Pathogenic Yeast Candida orthopsilosis by Separate Hybridizations between Two Parental Species

Mating between different species produces hybrids that are usually asexual and stuck as diploids, but can also lead to the formation of new species. Here, we report the genome sequences of 27 isolates of the pathogenic yeast Candida orthopsilosis. We find that most isolates are diploid hybrids, products of mating between two unknown parental species (A and B) that are 5% divergent in sequence. Isolates vary greatly in the extent of homogenization between A and B, making their genomes a mosaic of highly heterozygous regions interspersed with homozygous regions. Separate phylogenetic analyses of SNPs in the A- and B-derived portions of the genome produces almost identical trees of the isolates with four major clades. However, the presence of two mutually exclusive genotype combinations at the mating type locus, and recombinant mitochondrial genomes diagnostic of inter-clade mating, shows that the species C. orthopsilosis does not have a single evolutionary origin but was created at least four times by separate interspecies hybridizations between parents A and B. Older hybrids have lost more heterozygosity. We also identify two isolates with homozygous genomes derived exclusively from parent A, which are pure non-hybrid strains. The parallel emergence of the same hybrid species from multiple independent hybridization events is common in plant evolution, but is much less documented in pathogenic fungi.     

Introgression Threatens the Survival of the Critically Endangered Freshwater Crayfish Cherax tenuimanus (Decapoda: Parastacidae) in the Wild

Hybridization and genetic introgression following the introduction of exotic species can pose a significant threat to the survival of geographically restricted species. A remnant population of the critically endangered freshwater crayfish Cherax tenuimanus in the upper reaches of the Margaret River in southwestern Australia is under threat following the introduction and spread of its congener Cherax cainii. Here, we examine the extent of hybridization and introgression between the two species using twelve polymorphic microsatellite loci. Our study reveals there are three times more C. cainii than C. tenuimanus at our study site in the upper reaches of the Margaret River. There is also evidence of hybridization and introgression between C. tenuimanus and C. cainii at this site, with F₁, F₂ and backcrossed individuals identified. While interbreeding was confirmed in this study, our simulations suggest that the levels of introgression are much lower than would be expected under random mating, indicating partial reproductive barriers exist. Nevertheless, it is apparent that hybridization and introgression with C. cainii pose a serious threat to C. tenuimanus and their survival in the wild will require active adaptive management and continued genetic monitoring to evaluate management effectiveness.     

The Control Region of Maternally and Paternally Inherited Mitochondrial Genomes of Three Species of the Sea Mussel Genus Mytilus

Species of the mussel genus Mytilus possess maternally and paternally transmitted mitochondrial genomes. In the interbreeding taxa Mytilus edulis and M. galloprovincialis, several genomes of both types have been fully sequenced. The genome consists of the coding part (which, in addition to protein and RNA genes, contains several small noncoding sequences) and the main control region (CR), which in turn consists of three distinct parts: the first variable (VD1), the conserved (CD), and the second variable (VD2) domain. The maternal and paternal genomes are very similar in gene content and organization, even though they differ by >20% in primary sequence. They differ even more at VD1 and VD2, yet they are remarkably similar at CD. The complete sequence of a genome from the closely related species M. trossulus was previously reported and found to consist of a maternal-like coding part and a paternal-like and a maternal-like CR. From this and from the fact that it was extracted from a male individual, it was inferred that this is a genome that switched from maternal to paternal transmission. Here we provide clear evidence that this genome is the maternal genome of M. trossulus. We have found that in this genome the tRNA^Gln in the coding region is apparently defective and that an intact copy of this tRNA occurs in the CR, that one of the two conserved domains is missing essential motifs, and that one of the two first variable domains has a high rate of divergence. These features may explain the large size and mosaic structure of the CR of the maternal genome of M. trossulus. We have also obtained CR sequences of the maternal and paternal genomes of M. californianus, a more distantly related species. We compare the control regions from all three species, focusing on the divergence among genomes of different species origin and among genomes of different transmission routes.     

Interaction of endosperm size factors in maize

Crosses involving certain B-A translocations produce a reduced size of endosperm when those regions of the A chromosomes are missing in the sperm that fertilizes the polar nuclei. Previous studies involving the long arm of chromosome 10 showed that additional copies of this segment introduced through the maternal side could not rescue the reduced size phenotype conditioned by the corresponding deficiency in the paternal gamete. In this paper, experiments are described showing that other segments introduced maternally produce an even smaller kernel when fertilized by a sperm missing the same A chromosome segment or other ones that carry factors affecting endosperm size.—The example analyzed in detail involves reciprocal crosses between TB-4Sa and TB-10L19. Extra doses of 4S enhance the small kernel effect normally produced by TB-10L19. The additional copies of 4S have no effect on kernel mass when the 10L segment is present in the paternal contribution to the endosperm. The maternal enhancement by 4S is also effective with crosses by TB-1La but not by TB-1Sb. A survey of inter se crosses of B-A translocations shows that, when the maternal enhancement occurs, it is confined to those regions that themselves give a small kernel effect when used as a male. This correlation is strengthened by the observations that TB-10L19 enhances the small kernel effect of TB-1Sb, but TB-10L32 will not. Since these two B-10L translocations span the best localized small kernel effect region, this result supports the correlation of maternal enhancement regions with the paternal small kernel effect ones.—Because the enhancement can be attributed to a dosage effect and because the enhancement regions are coincident with the small kernel segments, it is postulated that this interacting system is analogous to aneuploid effects in diploid tissues but exhibits unique properties because of the evolutionary history and triploid condition of the endosperm.     

Expression of maternal and paternal histone genes during early cleavage stages of the echinoderm hybrid Strongylocentrotus purpuratus × Lytechinus pictus

Interspecific hybrids of the sea urchins Strongylocentrotus purpuratus (♀) and Lytechinus pictus (♂) were used to estimate the contributions of the maternal and paternal genomes to histone mRNA synthesis during early development. Radiolabeled histone mRNAs from the two sea urchin species were identified by hybridization to cloned histone genes from both S. purpuratus and L. pictus and shown to be electrophoretically distinguishable. The synthesis of maternal and paternal histone mRNA in these hybrid embryos is evident as early as the two-cell stage. By at least the 16-cell stage, both maternal and paternal histone mRNAs are associated with polysomes. The relative amounts of the maternal and paternal histone mRNAs synthesized by the zygote appear to be similar.     

Multiple lines of evidence supports the two varieties of Halenia elliptica (Gentianaceae) as two species

Delimiting species requires multiple sources of evidence. Here, we delimited two varieties of Halenia elliptica (Gentianaceae) using several lines of evidence, including morphological traits and mating system in a sympatric population, phylogenetic relationships based on nrITS and cpDNA (rpl16) data, and complete chloroplast genome sequences. Comparative analysis of 21 morphological traits clearly separates the two varieties of H. elliptica. Examination of the flowering process and pollination treatments indicate that H. elliptica var. grandiflora produces seeds via outcrossing, whereas H. elliptica var. elliptica produces seeds via mixed mating. Furthermore, hand-pollinated hybridization of the two varieties produced no seeds. Observations of pollinators showed that when bees began a pollination bout on H. elliptica var. grandiflora they preferred to continue pollinating this variety; however, when they began a pollination bout on H. elliptica var. elliptica, they showed no preference for either variety. Phylogenetic analysis confirmed the monophyly of H. elliptica, which was further divided into two monophyletic clades corresponding to the two varieties. A large number of variants from the chloroplast genomes reflected remarkable genetic dissimilarities between the two varieties of H. elliptica. We recommend that the two varieties of H. elliptica should be revised as two species (H. elliptica and H. grandiflora). Our findings indicate that H. elliptica varieties may have split into two separate species due to a shift in mating system, changes in flowering phenology and/or post-pollination reproductive isolation.     

Disruption of Imprinting by Mutator Transposon Insertions in the 5′ Proximal Regions of the Zea mays Mez1 Locus

Imprinting is a form of epigenetic gene regulation in which alleles are differentially regulated according to the parent of origin. The Mez1 gene in maize is imprinted such that the maternal allele is expressed in the endosperm while the paternal allele is not expressed. Three novel Mez1 alleles containing Mutator transposon insertions within the promoter were identified. These mez1-mu alleles do not affect vegetative expression levels or result in morphological phenotypes. However, these alleles can disrupt imprinted expression of Mez1. Maternal inheritance of the mez-m1 or mez1-m4 alleles results in activation of the normally silenced paternal allele of Mez1. Paternal inheritance of the mez1-m2 or mez1-m4 alleles can also result in a loss of silencing of the paternal Mez1 allele. The paternal disruption of imprinting by transposon insertions may reflect a requirement for sequence elements involved in targeting silencing of the paternal allele. The maternal disruption of imprinting by transposon insertions within the Mez1 promoter suggests that maternally produced MEZ1 protein may be involved in silencing of the paternal Mez1 allele. The endosperms with impaired imprinting did not exhibit phenotypic consequences associated with bi-allelic Mez1 expression.     

The Effect of Gametogenesis Regimes on the Chloroplast Genetic System of CHLAMYDOMONAS REINHARDTII

In Chlamydomonas reinhardtii, gamete differentiation is induced by nitrogen deprivation. While cellular nitrogen content and amount of chloroplast DNA in cells of both mating types are reduced during gametogenesis, the spontaneous transmission of paternal (mt^-) chloroplast alleles in crosses is specifically affected by the stringency of the nitrogen starvation regime used for pregrowth and gametogenesis of the mt^- parent. In all cases, reciprocal crosses yielded biparental zygospores whose clones contain predominantly cells expressing only the chloroplast alleles from the maternal (mt⁺) parent. No differences attributable to strain divergence were seen in chloroplast gene inheritance pattern, DNA content, or the relative frequency of transmission of paternal chloroplast alleles to progeny of biparental zygospores.     

Ploidy barrier to endosperm development in maize

Maize kernels inheriting the indeterminate gametophyte mutant (ig) on the female side had endosperms that ranged in ploidy level from diploid (2x) to nonaploid (9x). In crosses with diploid males, only kernels of the triploid endosperm class developed normally. Kernels of the tetraploid endosperm class were half-sized but with well-developed embryos that regularly germinated. Kernels of endosperm composition other than triploid or tetraploid were abortive.-Endosperm ploidy level resulting from mating ig/ig x tetraploid Ig similarly was variable. Most endosperms started to degenerate soon after pollination and remained in an arrested state. Hexaploid endosperm was exceptional; it developed normally during the sequence of stages studied and accounted for plump kernels on mature ears. Since such kernels have diploid maternal tissues (pericarp) but triploid embryos, the present finding favors the view that endosperm failure or success in such circumstances is governed by conditions within the endosperm itself.-Whereas tetraploid endosperm consisting of three maternal genomes and one paternal genome is slightly reduced in size but supports viable seed development, that endosperm having two maternal and two paternal chromosome sets was highly defective and conditioned abortion. Thus, development of maize endosperm evidently is affected by the parental source of its sets of chromosomes.