Individualized mating system estimation using genomic data

Using Runs of Homozygosity and Machine Learning to Disentangle Sources of Inbreeding and Infer Self-Fertilization Rates

Runs of homozygosity (ROHs) are indicative of elevated homozygosity and inbreeding due to mating of closely related individuals. Self-fertilization can be a major source of inbreeding which elevates genome-wide homozygosity and thus should also create long ROHs. While ROHs are frequently used to understand inbreeding in the context of conservation and selective breeding, as well as for consanguinity of populations and their demographic history, it remains unclear how ROH characteristics are altered by selfing and if this confounds expected signatures of inbreeding due to demographic change. Using simulations, we study the impact of the mode of reproduction and demographic history on ROHs. We apply random forests to identify unique characteristics of ROHs, indicative of different sources of inbreeding. We pinpoint distinct features of ROHs that can be used to better characterize the type of inbreeding the population was subjected to and to predict outcrossing rates and complex demographic histories. Using additional simulations and four empirical datasets, two from highly selfing species and two from mixed-maters, we predict the selfing rate and validate our estimations. We find that self-fertilization rates are successfully identified even with complex demography. Population genetic summary statistics improve algorithm accuracy particularly in the presence of additional inbreeding, e.g. from population bottlenecks. Our findings highlight the importance of ROHs in disentangling confounding factors related to various sources of inbreeding and demonstrate situations where such sources cannot be differentiated. Additionally, our random forest models provide a novel tool to the community for inferring selfing rates using genomic data.

Supergene potential of a selfish centromere

Selfishly evolving centromeres bias their transmission by exploiting the asymmetry of female meiosis and preferentially segregating to the egg. Such female meiotic drive systems have the potential to be supergenes, with multiple linked loci contributing to drive costs or enhancement. Here, we explore the supergene potential of a selfish centromere (D) in Mimulus guttatus, which was discovered in the Iron Mountain (IM) Oregon population. In the nearby Cone Peak population, D is still a large, non-recombining and costly haplotype that recently swept, but shorter haplotypes and mutational variation suggest a distinct population history. We detected D in five additional populations spanning more than 200 km; together, these findings suggest that selfish centromere dynamics are widespread in M. guttatus. Transcriptome comparisons reveal elevated differences in expression between driving and non-driving haplotypes within, but not outside, the drive region, suggesting large-scale cis effects of D's spread on gene expression. We use the expression data to refine linked candidates that may interact with drive, including Nuclear Autoantigenic Sperm Protein (NASP^SIM3), which chaperones the centromere-defining histone CenH3 known to modify Mimulus drive. Together, our results show that selfishly evolving centromeres may exhibit supergene behaviour and lay the foundation for future genetic dissection of drive and its costs. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

The impact of urbanization on outcrossing rate and population genetic variation in the native wildflower, Impatiens capensis

Cities are one of the fastest growing ecosystems on the planet, and conserving urban biodiversity is of primary importance. Urbanization increases habitat fragmentation and may be particularly problematic for native plant species which often exist in small, remnant populations in cities. We studied the effects of urbanization on Impatiens capensis, a self-compatible native wildflower, which is an important nectar and pollen source for native bees and hummingbirds. We sampled I. capensis from six populations located in urban and rural habitats in Toronto, Ontario, Canada. We sequenced the DNA of 43 families (N = 86 individuals) using genotype-by-sequencing to obtain 5627 single nucleotide polymorphisms. From each parent and offspring, we estimated individual outcrossing rates, population-level genetic diversity and genetic structure among populations. We found that 95% of plants were outcrossed, and populations were genetically differentiated, where urban populations contained a subset of the genetic variation found in rural populations. Urban populations exhibited lower genetic diversity than rural populations, and we detected a relationship between population census size and habitat on genetic diversity. Despite high outcrossing rates, our results suggest that urbanization reduces the genetic diversity of I. capensis populations, potentially increasing the vulnerability of these populations to long-term population declines and extirpation in response to urbanization.

Population genomic consequences of life-history and mating system adaptation to a geothermal soil mosaic in yellow monkeyflowers

Local selection can promote phenotypic divergence despite gene flow across habitat mosaics, but adaptation itself may generate substantial barriers to genetic exchange. In plants, life-history, phenology, and mating system divergence have been proposed to promote genetic differentiation in sympatry. In this study, we investigate phenotypic and genetic variation in Mimulus guttatus (yellow monkeyflowers) across a geothermal soil mosaic in Yellowstone National Park (YNP). Plants from thermal annual and nonthermal perennial habitats were heritably differentiated for life-history and mating system traits, consistent with local adaptation to the ephemeral thermal-soil growing season. However, genome-wide genetic variation primarily clustered plants by geographic region, with little variation sorting by habitat. The one exception was an extreme thermal population also isolated by a 200 m geographical gap of no intermediate habitat. Individual inbreeding coefficients (F_IS ) were higher (and predicted by trait variation) in annual plants and annual pairs showed greater isolation by distance at local (<1 km) scales. Finally, YNP adaptation does not reuse a widespread inversion that underlies M. guttatus life-history ecotypes range-wide, suggesting a novel genetic mechanism. Overall, this work suggests that life-history and mating system adaptation strong enough to shape individual mating patterns does not necessarily generate incipient speciation without geographical barriers.

Adaptive and nonadaptive causes of heterogeneity in genetic differentiation across the Mimulus guttatus genome

Genetic diversity becomes structured among populations over time due to genetic drift and divergent selection. Although population structure is often treated as a uniform underlying factor, recent resequencing studies of wild populations have demonstrated that diversity in many regions of the genome may be structured quite dissimilar to the genome-wide pattern. Here, we explored the adaptive and nonadaptive causes of such genomic heterogeneity using population-level, whole genome resequencing data obtained from annual Mimulus guttatus individuals collected across a rugged environment landscape. We found substantial variation in how genetic differentiation is structured both within and between chromosomes, although, in contrast to other studies, known inversion polymorphisms appear to serve only minor roles in this heterogeneity. In addition, much of the genome can be clustered into eight among-population genetic differentiation patterns, but only two of these clusters are particularly consistent with patterns of isolation by distance. By performing genotype-environment association analysis, we also identified genomic intervals where local adaptation to specific climate factors has accentuated genetic differentiation among populations, and candidate genes in these windows indicate climate adaptation may proceed through changes affecting specialized metabolism, drought resistance, and development. Finally, by integrating our findings with previous studies, we show that multiple aspects of plant reproductive biology may be common targets of balancing selection and that variants historically involved in climate adaptation among populations have probably also fuelled rapid adaptation to microgeographic environmental variation within sites.

From pollen dispersal to plant diversification: genetic consequences of pollination mode

Pollinators influence patterns of plant speciation, and one intuitive hypothesis is that pollinators affect rates of plant diversification through their effects on pollen dispersal. By specifying mating events and pollen flow across the landscape, distinct types of pollinators may cause different opportunities for allopatric speciation. This pollen dispersal-dependent speciation hypothesis predicts that pollination mode has effects on the spatial context of mating events that scale up to impact population structure and rates of species formation. Here I consider recent comparative studies, including genetic analyses of plant mating events, population structure and comparative phylogenetic analyses, to examine evidence for this model. These studies suggest that highly mobile pollinators conduct greater gene flow within and among populations, compared to less mobile pollinators. These differences influence patterns of population structure across the landscape. However, the effects of pollination mode on speciation rates is less predictable. In some contexts, the predicted effects of pollen dispersal are outweighed by other factors that govern speciation rates. A multiscale approach to examine effects of pollination mode on plant mating system, population structure and rates of diversification is key to determining the role of pollen dispersal on plant speciation for model clades.

How and When Does Outcrossing Occur in the Predominantly Selfing Species Medicago truncatula?

Empirical studies on natural populations of Medicago truncatula revealed selfing rates higher than 80%, but never up to 100%. Similarly, several studies of predominantly selfing species show variability in the level of residual outcrossing between populations and also between temporal samples of the same population. However, these studies measure global selfing rates at the scale of the population and we do not know whether there is intra-population variation and how outcrossing events are distributed, between genotypes, plants, flowers, or seeds. Theoretical studies predict the maintenance of residual outcrossing in highly selfing species due to environmental (e.g., pollen biology) and/or genetic determinants and decompositions of the variation in outcrossing rate using experimental data can be very informative to test these hypotheses. Here, we focus on one natural population of M. truncatula in order to describe precisely its mating system. In particular, we investigated the determinants of the selfing rate by testing for seasonal variations (environmental determinism) and variations between genotypes (genetic determinism). We measured selfing rates in maternal progenies from plants collected widely across a natural population. For each plant, we collected pods from flowers produced at the beginning and at the end of the flowering season to test for a seasonal variation in the outcrossing rate. For each collected offspring, we also estimated the likelihood that it was issued from a self-fertilization event and assessed the genetic component of variation of this mating system measure. We found a significant, albeit small, increase in outcrossing rate in progenies collected at the end [t _m = 0.137 (SD = 0.025)] compared to those collected at the beginning [t _m = 0.083 (0.016)] of the flowering season. A significant between genotypes variation in selfing rate was also detected, resulting in a heritability of 9% for the rate of residual outcrossing. Altogether, our work shows that despite a predominantly selfing reproductive mode, M. truncatula displays variation in residual outcrossing rate, and that this trait is likely under a complex determinism combining environmental and genetic factors. We discuss the evolutionary implications of our results for the population.

Predicting evolutionary change at the DNA level in a natural Mimulus population

Evolution by natural selection occurs when the frequencies of genetic variants change because individuals differ in Darwinian fitness components such as survival or reproductive success. Differential fitness has been demonstrated in field studies of many organisms, but it remains unclear how well we can quantitatively predict allele frequency changes from fitness measurements. Here, we characterize natural selection on millions of Single Nucleotide Polymorphisms (SNPs) across the genome of the annual plant Mimulus guttatus. We use fitness estimates to calibrate population genetic models that effectively predict allele frequency changes into the next generation. Hundreds of SNPs experienced "male selection" in 2013 with one allele at each SNP elevated in frequency among successful male gametes relative to the entire population of adults. In the following generation, allele frequencies at these SNPs consistently shifted in the predicted direction. A second year of study revealed that SNPs had effects on both viability and reproductive success with pervasive trade-offs between fitness components. SNPs favored by male selection were, on average, detrimental to survival. These trade-offs (antagonistic pleiotropy and temporal fluctuations in fitness) may be essential to the long-term maintenance of alleles. Despite the challenges of measuring selection in the wild, the strong correlation between predicted and observed allele frequency changes suggests that population genetic models have a much greater role to play in forward-time prediction of evolutionary change.

How rapidly do self-compatible populations evolve selfing? Mating system estimation within recently evolved self-compatible populations of Azorean Tolpis succulenta (Asteraceae)

Genome-wide genotyping and Bayesian inference method (BORICE) were employed to estimate outcrossing rates and paternity in two small plant populations of Tolpis succulenta (Asteraceae) on Graciosa island in the Azores. These two known extant populations of T. succulenta on Graciosa have recently evolved self-compatibility. Despite the expectation that selfing would occur at an appreciable rate (self-incompatible populations of the same species show low but nonzero selfing), high outcrossing was found in progeny arrays from maternal plants in both populations. This is inconsistent with an immediate transition to high selfing following the breakdown of a genetic incompatibility system. This finding is surprising given the small population sizes and the recent colonization of an island from self-incompatible colonists of T. succulenta from another island in the Azores, and a potential paucity of pollinators, all factors selecting for selfing through reproductive assurance. The self-compatible lineage(s) likely have high inbreeding depression (ID) that effectively halts the evolution of increased selfing, but this remains to be determined. Like their progeny, all maternal plants in both populations are fully outbred, which is consistent with but not proof of high ID. High multiple paternity was found in both populations, which may be due in part to the abundant pollinators observed during the flowering season.

Genome-wide genotyping estimates mating system parameters and paternity in the island species Tolpis succulenta

PREMISE

The mating system has profound consequences, not only for ecology and evolution, but also for the conservation of threatened or endangered species. Unfortunately, small populations are difficult to study owing to limits on sample size and genetic marker diversity. Here, we estimated mating system parameters in three small populations of an island plant using genomic genotyping. Although self-incompatible (SI) species are known to often set some self-seed, little is known about how "leaky SI" affects selfing rates in nature or the role that multiple paternity plays in small populations.

METHODS

We generalized the BORICE mating system program to determine the siring pattern within maternal families. We applied this algorithm to maternal families from three populations of Tolpis succulenta from Madeira Island and genotyped the progeny using RADseq. We applied BORICE to estimate each individual offspring as outcrossed or selfed, the paternity of each outcrossed offspring, and the level of inbreeding of each maternal plant.

RESULTS

Despite a functional self-incompatibility system, these data establish T. succulenta as a pseudo-self-compatible (PSC) species. Two of 75 offspring were strongly indicated as products of self-fertilization. Despite selfing, all adult maternal plants were fully outbred. There was high differentiation among and low variation within populations, consistent with a history of genetic isolation of these small populations. There were generally multiple sires per maternal family. Twenty-two percent of sib contrasts (between outcrossed offspring within maternal families) shared the same sire.

CONCLUSIONS

Genome-wide genotyping, combined with appropriate analytical methods, enables estimation of mating system and multiple paternity in small populations. These data address questions about the evolution of reproductive traits and the conservation of threatened populations.