Transmission ratio distortion in Arabidopsis lyrata: effects of population divergence and the S-locus

Genomic mechanisms and consequences of diverse postzygotic barriers between monkeyflower species

The evolution of genomic incompatibilities causing postzygotic barriers to hybridization is a key step in species divergence. Incompatibilities take 2 general forms-structural divergence between chromosomes leading to severe hybrid sterility in F1 hybrids and epistatic interactions between genes causing reduced fitness of hybrid gametes or zygotes (Dobzhansky-Muller incompatibilities). Despite substantial recent progress in understanding the molecular mechanisms and evolutionary origins of both types of incompatibility, how each behaves across multiple generations of hybridization remains relatively unexplored. Here, we use genetic mapping in F2 and recombinant inbred line (RIL) hybrid populations between the phenotypically divergent but naturally hybridizing monkeyflowers Mimulus cardinalis and M. parishii to characterize the genetic basis of hybrid incompatibility and examine its changing effects over multiple generations of experimental hybridization. In F2s, we found severe hybrid pollen inviability (<50% reduction vs parental genotypes) and pseudolinkage caused by a reciprocal translocation between Chromosomes 6 and 7 in the parental species. RILs retained excess heterozygosity around the translocation breakpoints, which caused substantial pollen inviability when interstitial crossovers had not created compatible heterokaryotypic configurations. Strong transmission ratio distortion and interchromosomal linkage disequilibrium in both F2s and RILs identified a novel 2-locus genic incompatibility causing sex-independent gametophytic (haploid) lethality. The latter interaction eliminated 3 of the expected 9 F2 genotypic classes via F1 gamete loss without detectable effects on the pollen number or viability of F2 double heterozygotes. Along with the mapping of numerous milder incompatibilities, these key findings illuminate the complex genetics of plant hybrid breakdown and are an important step toward understanding the genomic consequences of natural hybridization in this model system.

Haploid selection, sex ratio bias, and transitions between sex-determining systems

Sex determination is remarkably dynamic; many taxa display shifts in the location of sex-determining loci or the evolution of entirely new sex-determining systems. Predominant theories for why we observe such transitions generally conclude that novel sex-determining systems are favoured by selection if they equalise the sex ratio or increase linkage with a locus that experiences different selection in males versus females. We use population genetic models to extend these theories in two ways: (1) We consider the dynamics of loci very tightly linked to the ancestral sex-determining loci, e.g., within the nonrecombining region of the ancestral sex chromosomes. Variation at such loci can favour the spread of new sex-determining systems in which the heterogametic sex changes (XY to ZW or ZW to XY) and the new sex-determining region is less closely linked (or even unlinked) to the locus under selection. (2) We consider selection upon haploid genotypes either during gametic competition (e.g., pollen competition) or meiosis (i.e., nonmendelian segregation), which can cause the zygotic sex ratio to become biased. Haploid selection can drive transitions between sex-determining systems without requiring selection to act differently in diploid males versus females. With haploid selection, we find that transitions between male and female heterogamety can evolve so that linkage with the sex-determining locus is either strengthened or weakened. Furthermore, we find that sex ratio biases may increase or decrease with the spread of new sex chromosomes, which implies that transitions between sex-determining systems cannot be simply predicted by selection to equalise the sex ratio. In fact, under many conditions, we find that transitions in sex determination are favoured equally strongly in cases in which the sex ratio bias increases or decreases. Overall, our models predict that sex determination systems should be highly dynamic, particularly when haploid selection is present, consistent with the evolutionary lability of this trait in many taxa.

Role of seed germination in adaptation and reproductive isolation in Arabidopsis lyrata

Seed germination is an important developmental and life history stage. Yet, the evolutionary impact of germination has mainly been studied in the context of dormancy, or for its role in reproductive isolation between species. Here, we aim to examine multiple consequences of genetic divergence on germination traits between two Arabidopsis lyrata subspecies: ssp. petraea (Eurasia) and ssp. lyrata (North America). Postdormancy germination time, a potentially adaptive trait, showed differentiation between the populations, and quantitative trait loci (QTL) mapping revealed that the trait variation is mainly controlled by two antagonistic loci. These QTL areas contain several candidate genes with known function in postdormancy germination in A. thaliana. The sequence variation of three genes was consistent with differential selection, and they also included fixed nonsynonymous substitutions with potential to account for the phenotypic differentiation. We further show that the divergence between the subspecies has led to a slight but significant reduction in hybrid germination proportions, indicating incipient reproductive isolation. Comparison of reciprocal F₁ and F₂ progenies suggests that Bateson-Dobzhansky-Muller incompatibilities likely act through uniparentally inherited factors. Examination of genomewide transmission ratio distortion further revealed that cytonuclear interactions cause substantial pregermination inviability in the hybrids. These results confirm that seed germination has adaptive potential beyond the dormancy stage and that hybrid seed inviability can be one of the first reproductive barriers to arise during divergence.

Paternal-specific S-allele transmission in sweet cherry (Prunus avium L.): the potential for sexual selection

Homomorphic self-incompatibility is a well-studied example of a physiological process that is thought to increase population diversity and reduce the expression of inbreeding depression. Whereas theoretical models predict the presence of a large number of S-haplotypes with equal frequencies at equilibrium, unequal allele frequencies have been repeatedly reported and attributed to sampling effects, population structure, demographic perturbation, sheltered deleterious mutations or selection pressure on linked genes. However, it is unclear to what extent unequal segregations are the results of gametophytic or sexual selection. Although these two forces are difficult to disentangle, testing S-alleles in the offspring of controlled crosses provides an opportunity to separate these two phenomena. In this work, segregation and transmission of S-alleles have been characterized in progenies of mixed donors and fully compatible pollinations under field conditions in Prunus avium. Seed set patterns and pollen performance have also been characterized. The results reveal paternal-specific distorted transmission of S-alleles in most of the crosses. Interestingly, S-allele segregation within any given paternal or maternal S-locus was random. Observations on pollen germination, pollen tube growth rate, pollen tube cohort size, seed set dynamics and transmission patterns strongly suggest post-pollination, prezygotic sexual selection, with male-male competition as the most likely mechanism. According to these results, post-pollination sexual selection takes precedence over frequency-dependent selection in explaining unequal S-haplotype frequencies.

Genetics of cryptic speciation within an Arctic mustard, Draba nivalis

Crossing experiments indicate that hybrid sterility barriers frequently have developed within diploid, circumpolar plant species of the genus Draba. To gain insight into the rapid evolution of postzygotic reproductive isolation in this system, we augmented the linkage map of one of these species, D. nivalis, and searched for quantitative trait loci (QTLs) associated with reproductive isolation. The map adds 63 new dominant markers to a previously published dataset of 31 co-dominant microsatellites. These markers include 52 amplified fragment length polymorphisms (AFLPs) and 11 sequence-specific amplified polymorphisms (SSAPs) based on retrotransposon sequence. 22 markers displaying transmission ratio distortion were further included in the map. We resolved eight linkage groups with a total map length of 894 cM. Significant genotype-trait associations, or quantitative trait loci (QTL), were detected for reproductive phenotypes including pollen fertility (4 QTLs), seed set (3 QTLs), flowering time (3 QTLs) and number of flowers (4 QTLs). Observed patterns of inheritance were consistent with the influence of both nuclear-nuclear interactions and chromosomal changes on these traits. All seed set QTLs and one pollen fertility QTL displayed underdominant effects suggestive of the involvement of chromosomal rearrangements in hybrid sterility. Interestingly, D. nivalis is predominantly self-fertilizing, which may facilitate the establishment of underdominant loci and contribute to reproductive isolation.

Investigating incipient speciation in Arabidopsis lyrata from patterns of transmission ratio distortion

Our understanding of the development of intrinsic reproductive isolation is still largely based on theoretical models and thorough empirical studies on a small number of species. Theory suggests that reproductive isolation develops through accumulation of epistatic genic incompatibilities, also known as Bateson-Dobzhansky-Muller (BDM) incompatibilities. We can detect these from marker transmission ratio distortion (TRD) in hybrid progenies of crosses between species or populations, where TRD is expected to result from selection against heterospecific allele combinations in hybrids. TRD may also manifest itself because of intragenomic conflicts or competition between gametes or zygotes. We studied early stage speciation in Arabidopsis lyrata by investigating patterns of TRD across the genome in F2 progenies of three reciprocal crosses between four natural populations. We found that the degree of TRD increases with genetic distance between crossed populations, but also that reciprocal progenies may differ substantially in their degree of TRD. Chromosomes AL6 and especially AL1 appear to be involved in many single- and two-locus distortions, but the location and source of TRD vary between crosses and between reciprocal progenies. We also found that the majority of single- and two-locus TRD appears to have a gametic, as opposed to zygotic, origin. Thus, while theory on BDM incompatibilities is typically illustrated with derived nuclear alleles proving incompatible in hybrid zygotes, our results suggest a prominent role for distortions emerging before zygote formation.

Time scales of divergence and speciation among natural populations and subspecies of Arabidopsis lyrata (Brassicaceae)

PREMISE OF THE STUDY

Plant populations that face new environments adapt and diverge simultaneously, and both processes leave footprints in their genetic diversity. Arabidopsis lyrata is an excellent species for studying these processes. Pairs of populations and subspecies of A. lyrata represent different stages of divergence. These populations are also known to be locally adapted and display various stages of emerging reproductive isolation.

METHODS

We used nucleotide diversity data from 19 loci to estimate divergence times and levels of diversity among nine A. lyrata populations. Traditional distance-based methods and model-based clustering analysis were used to supplement pairwise coalescence-based analysis of divergence.

KEY RESULTS

Estimated divergence times varied from 130,000 generations between North American and European subspecies to 39,000 generations between central European and Scandinavian populations. In concordance with previous studies, the highest level of diversity was found in Central Europe and the lowest in North America and a diverged Russian Karhumäki population. Local adaptation among Northern and central European populations has emerged during the last 39,000 generations. Populations that are reproductively isolated by prezygotic mechanisms have been separated for a longer time period of ∼70,000 generations but still have shared nucleotide polymorphism.

CONCLUSIONS

In A. lyrata, reproductively isolated populations started to diverge ∼70,000 generations ago and more closely related, locally adapted populations have been separate lineages for ∼39,000 generations. However, based on the posterior distribution of divergence times, the processes leading to reproductive isolation and local adaptation are likely to temporally coincide.

An EST-derived SNP and SSR genetic linkage map of cassava (Manihot esculenta Crantz)

Cassava (Manihot esculenta Crantz) is one of the most important food security crops in the tropics and increasingly being adopted for agro-industrial processing. Genetic improvement of cassava can be enhanced through marker-assisted breeding. For this, appropriate genomic tools are required to dissect the genetic architecture of economically important traits. Here, a genome-wide SNP-based genetic map of cassava anchored in SSRs is presented. An outbreeder full-sib (F1) family was genotyped on two independent SNP assay platforms: an array of 1,536 SNPs on Illumina's GoldenGate platform was used to genotype a first batch of 60 F1. Of the 1,358 successfully converted SNPs, 600 which were polymorphic in at least one of the parents and was subsequently converted to KBiosciences' KASPar assay platform for genotyping 70 additional F1. High-precision genotyping of 163 informative SSRs using capillary electrophoresis was also carried out. Linkage analysis resulted in a final linkage map of 1,837 centi-Morgans (cM) containing 568 markers (434 SNPs and 134 SSRs) distributed across 19 linkage groups. The average distance between adjacent markers was 3.4 cM. About 94.2% of the mapped SNPs and SSRs have also been localized on scaffolds of version 4.1 assembly of the cassava draft genome sequence. This more saturated genetic linkage map of cassava that combines SSR and SNP markers should find several applications in the improvement of cassava including aligning scaffolds of the cassava genome sequence, genetic analyses of important agro-morphological traits, studying the linkage disequilibrium landscape and comparative genomics.

The recombination landscape in Arabidopsis thaliana F2 populations

Recombination during meiosis shapes the complement of alleles segregating in the progeny of hybrids, and has important consequences for phenotypic variation. We examined allele frequencies, as well as crossover (XO) locations and frequencies in over 7000 plants from 17 F(2) populations derived from crosses between 18 Arabidopsis thaliana accessions. We observed segregation distortion between parental alleles in over half of our populations. The potential causes of distortion include variation in seed dormancy and lethal epistatic interactions. Such a high occurrence of distortion was only detected here because of the large sample size of each population and the number of populations characterized. Most plants carry only one or two XOs per chromosome pair, and therefore inherit very large, non-recombined genomic fragments from each parent. Recombination frequencies vary between populations but consistently increase adjacent to the centromeres. Importantly, recombination rates do not correlate with whole-genome sequence differences between parental accessions, suggesting that sequence diversity within A. thaliana does not normally reach levels that are high enough to exert a major influence on the formation of XOs. A global knowledge of the patterns of recombination in F(2) populations is crucial to better understand the segregation of phenotypic traits in hybrids, in the laboratory or in the wild.

Nuclear-cytoplasmic interactions reduce male fertility in hybrids of Arabidopsis lyrata subspecies

We examined the level of postzygotic reproductive isolation in F(1) and F(2) hybrids of reciprocal crosses between the Arabidopsis lyrata subspecies lyrata (North American) and petraea (European). Our main results are: first, the percentage of fertile pollen was significantly reduced in the F(1) and F(2) compared to the parental populations. Second, mean pollen fertility differed markedly between reciprocal crosses: 84% in the F(2) with ssp. lyrata cytoplasm and 61% in the F(2) with ssp. petraea cytoplasm. Third, 17% of the F(2) with ssp. petraea cytoplasm showed male sterility (produced less than 30 pollen grains in our subsample). The hybrids were female fertile. We used QTL mapping to find the genomic regions that determine pollen fertility and that restore cytoplasmic male sterility (CMS). In the F(2) with ssp. lyrata cytoplasm, an epistatic pair of QTLs was detected. In the reciprocal F(2) progeny, four QTLs demonstrated within-population polymorphism for hybrid male sterility. In addition, in the F(2) with ssp. petraea cytoplasm, there was a strong male fertility restorer locus on chromosome 2 where a cluster of CMS restorer gene-related PPR genes have been found in A. lyrata. Our results underline the importance of cytonuclear interactions in understanding genetics of the early stages of speciation.

Progress and Promise in using Arabidopsis to Study Adaptation, Divergence, and Speciation

Fundamental questions remain to be answered on how lineages split and new species form. The Arabidopsis genus, with several increasingly well characterized species closely related to the model system A. thaliana, provides a rare opportunity to address key questions in speciation research. Arabidopsis species, and in some cases populations within a species, vary considerably in their habitat preferences, adaptations to local environments, mating system, life history strategy, genome structure and chromosome number. These differences provide numerous open doors for understanding the role these factors play in population divergence and how they may cause barriers to arise among nascent species. Molecular tools available in A. thaliana are widely applicable to its relatives, and together with modern comparative genomic approaches they will provide new and increasingly mechanistic insights into the processes underpinning lineage divergence and speciation. We will discuss recent progress in understanding the molecular basis of local adaptation, reproductive isolation and genetic incompatibility, focusing on work utilizing the Arabidopsis genus, and will highlight several areas in which additional research will provide meaningful insights into adaptation and speciation processes in this genus.

The sheltered genetic load linked to the s locus in plants: new insights from theoretical and empirical approaches in sporophytic self-incompatibility

Inbreeding depression and mating systems evolution are closely linked, because the purging of deleterious mutations and the fitness of individuals may depend on outcrossing vs. selfing rates. Further, the accumulation of deleterious mutations may vary among genomic regions, especially for genes closely linked to loci under balancing selection. Sporophytic self-incompatibility (SSI) is a common genetic mechanism in angiosperm that enables hermaphrodite plants to avoid selfing and promote outcrossing. The SSI phenotype is determined by the S locus and may depend on dominance relationships among alleles. Since most individuals are heterozygous at the S locus and recombination is suppressed in the S-locus region, it has been suggested that deleterious mutations could accumulate at genes linked to the S locus, generating a "sheltered load." In this article, we first theoretically investigate the conditions generating sheltered load in SSI. We show that deleterious mutations can accumulate in linkage with specific S alleles, and particularly if those S alleles are dominant. Second, we looked for the presence of sheltered load in Arabidopsis halleri using CO(2) gas treatment to overcome self-incompatibility. By examining the segregation of S alleles and measuring the relative fitness of progeny, we found significant sheltered load associated with the most dominant S allele (S15) of three S alleles tested. This sheltered load seems to be expressed at several stages of the life cycle and to have a larger effect than genomic inbreeding depression.

Genetic linkage map for Amylostereum areolatum reveals an association between vegetative growth and sexual and self-recognition

Amylostereum areolatum is a filamentous fungus that grows through tip extension, branching and hyphal fusion. In the homokaryotic phase, the hyphae of different individuals are capable of fusing followed by heterokaryon formation, only if they have dissimilar allelic specificities at their mating-type (mat) loci. In turn, hyphal fusion between heterokaryons persists only when they share the same alleles at all of their heterokaryon incompatibility (het) loci. In this study we present the first genetic linkage map for A. areolatum, onto which the mat and het loci, as well as quantitative trait loci (QTLs) for mycelial growth rate are mapped. The recognition loci (mat-A and het-A) are positioned near QTLs associated with mycelial growth, suggesting that the genetic determinants influencing recognition and growth rate in A. areolatum are closely associated. This was confirmed when isolates associated with specific mat and het loci displayed significantly different mycelial growth rates. Although the link between growth and sexual recognition has previously been observed in other fungi, this is the first time that an association between growth and self-recognition has been shown.

Unequal segregation of SRK alleles at the S locus in Brassica cretica

In the Brassicaceae plant family, which includes the Arabidopsis and Brassica genera, self-incompatibility (SI) is controlled by genes at the S locus. Using experimental crosses, we studied the pattern of inheritance of S-locus alleles in the wild species Brassica cretica. Four full-sib families were established and unequal segregation of alleles at the SRK SI gene was found in one family. The segregation distortion acted in favour of a recessive (class II) allele and was best explained by some form of gametic-level selection. Our findings are discussed in the light of theoretical predictions of differential accumulation of deleterious mutations among S-locus alleles.