Assessing the utility of SoilGrids250 for biogeographic inference of plant populations.

Inclusion of edaphic conditions in biogeographical studies typically provides a better fit and deeper understanding of plant distributions. Increased reliance on soil data calls for easily accessible data layers providing continuous soil predictions worldwide. Although SoilGrids provides a potentially useful source of predicted soil data for biogeographic applications, its accuracy for estimating the soil characteristics experienced by individuals in small-scale populations is unclear. We used a biogeographic sampling approach to obtain soil samples from 212 sites across the midwestern and eastern United States, sampling only at sites where there was a population of one of the 22 species in Lobelia sect. Lobelia. We analyzed six physical and chemical characteristics in our samples and compared them with predicted values from SoilGrids. Across all sites and species, soil texture variables (clay, silt, sand) were better predicted by SoilGrids (R 2: .25-.46) than were soil chemistry variables (carbon and nitrogen, R 2 ≤ .01; pH, R 2: .19). While SoilGrids predictions rarely matched actual field values for any variable, we were able to recover qualitative patterns relating species means and population-level plant characteristics to soil texture and pH. Rank order of species mean values from SoilGrids and direct measures were much more consistent for soil texture (Spearman r S = .74-.84; all p < .0001) and pH (r S = .61, p = .002) than for carbon and nitrogen (p > .35). Within the species L. siphilitica, a significant association, known from field measurements, between soil texture and population sex ratios could be detected using SoilGrids data, but only with large numbers of sites. Our results suggest that modeled soil texture values can be used with caution in biogeographic applications, such as species distribution modeling, but that soil carbon and nitrogen contents are currently unreliable, at least in the region studied here.

Do Genetic Drift and Gene Flow Affect the Geographic Distribution of Female Plants in Gynodioecious Lobelia siphilitica?

Variation in population sex ratio is particularly pronounced in gynodioecious angiosperms. Extremely high female frequencies in gynodioecious populations cannot be readily explained by selective forces alone. To assess the contributions of drift and gene flow to extreme sex-ratio variation, we documented sex ratio and population size in 92 populations of Lobelia siphilitica across its range and genotyped plants using plastid and nuclear genetic markers. Similarity in spatial patterns of genetic and demographic variables may suggest that drift and/or gene flow have contributed to population sex-ratio variation in L. siphilitica. We found strong spatial structuring of extremely high female frequencies: populations with >50% female plants are restricted to the south−central portion of the range. However, we did not detect any spatial structuring in population size nor metrics of genetic diversity, suggesting that extreme variation in female frequency is not strongly affected by drift or gene flow. Extreme sex-ratio variation is frequently observed in gynodioecious plants, but its causes are difficult to identify. Further investigation into mechanisms that create or maintain the spatial structure of sex ratios in gynodioecious species will provide much needed insight.

Beyond balancing selection: frequent mitochondrial recombination contributes to high-female frequencies in gynodioecious Lobelia siphilitica (Campanulaceae).

Gynodioecy is a sexual system in which females and hermaphrodites co-occur. In most gynodioecious angiosperms, sex is determined by an interaction between mitochondrial male-sterility genes (CMS) that arise via recombination and nuclear restorer alleles that evolve to suppress them. In theory, gynodioecy occurs when multiple CMS types are maintained at equilibrium frequencies by balancing selection. However, some gynodioecious populations contain very high frequencies of females. High female frequencies are not expected under balancing selection, but could be explained by the repeated introduction of novel CMS types. To test for balancing selection and/or the repeated introduction of novel CMS, we characterised cytoplasmic haplotypes from 61 populations of Lobelia siphilitica that vary widely in female frequency. We confirmed that mitotype diversity and female frequency were positively correlated across populations, consistent with balancing selection. However, while low-female populations hosted mostly common mitotypes, high-female populations and female plants hosted mostly rare, recombinant mitotypes likely to carry novel CMS types. Our results suggest that balancing selection maintains established CMS types across this species, but extreme female frequencies result from frequent invasion by novel CMS types. We conclude that balancing selection alone cannot account for extreme population sex-ratio variation within a gynodioecious species.

Physiological effects of temperature do not explain prevalence of females in populations of gynodioecious Lobelia siphilitica growing in warmer climates.

PREMISE OF THE STUDY: Gynodioecy is a sexual polymorphism whereby female and hermaphroditic plants co-occur within populations. In many gynodioecious species, stressful abiotic environments are associated with higher frequencies of females. This association suggests that abiotic stress affects the relative fitness of females and hermaphrodites and, thus, the maintenance of gynodioecy. METHODS: To test whether abiotic stress affects the fitness of females and hermaphrodites, we grew open-pollinated Lobelia siphilitica families in temperature regimes characteristic of the southern portion of the species' range (where females are common) and the northern portion of the range (where females are rare). We measured physiological and phenological traits that are indicative of heat stress, and fitness components of females and hermaphrodites that could affect the maintenance of gynodioecy. KEY RESULTS: Contrary to expectations if growth at high temperatures is stressful, we found that the hot treatment increased leaf chlorophyll content, decreased the percentage of plants that delayed flowering initiation, and did not affect the quantum efficiency of photosystem II. Growth at high temperatures did not affect the magnitude of the difference in rosette size (a correlate of flower number) between females and hermaphrodites, or the variance in pollen viability among hermaphrodites. CONCLUSIONS: We found that growing-season temperatures typical of high female L. siphilitica populations were not stressful and did not affect either the fitness of females compared to hermaphrodites or variation in fitness among hermaphrodites. Consequently, further research is necessary to explain correlations between abiotic environmental factors and the frequency of females in this and other gynodioecious species.

Selfish evolution of cytonuclear hybrid incompatibility in Mimulus.

Intraspecific coevolution between selfish elements and suppressors may promote interspecific hybrid incompatibility, but evidence of this process is rare. Here, we use genomic data to test alternative models for the evolution of cytonuclear hybrid male sterility in Mimulus In hybrids between Iron Mountain (IM) Mimulus guttatus × Mimulus nasutus, two tightly linked M. guttatus alleles (Rf1/Rf2) each restore male fertility by suppressing a local mitochondrial male-sterility gene (IM-CMS). Unlike neutral models for the evolution of hybrid incompatibility loci, selfish evolution predicts that the Rf alleles experienced strong selection in the presence of IM-CMS. Using whole-genome sequences, we compared patterns of population-genetic variation in Rf at IM to a neighbouring population that lacks IM-CMS. Consistent with local selection in the presence of IM-CMS, the Rf region shows elevated FST, high local linkage disequilibrium and a distinct haplotype structure at IM, but not at Cone Peak (CP), suggesting a recent sweep in the presence of IM-CMS. In both populations, Rf2 exhibited lower polymorphism than other regions, but the low-diversity outliers were different between CP and IM. Our results confirm theoretical predictions of ubiquitous cytonuclear conflict in plants and provide a population-genetic mechanism for the evolution of a common form of hybrid incompatibility.

Why is gynodioecy a rare but widely distributed sexual system? Lessons from the Lamiaceae.

Gynodioecy, a sexual system where females and hermaphrodites co-occur, is found in << 1% of angiosperm species. To understand why gynodioecy is rare, we need to understand why females are maintained in some lineages, but not in others. We modelled the evolution of gynodioecy in the Lamiaceae, and investigated whether transition rates between gynodioecious and nongynodioecious states varied across the family. We also investigated whether the evolution of gynodioecy was correlated with the evolution of a herbaceous growth form and temperate distribution. Transition rates differed between Lamiaceae subfamilies. In the Nepetoideae, there were many transitions towards gynodioecy (n = 11), but also many reversions to nongynodioecy (n = 29). In addition, a herbaceous growth form, but not a temperate distribution, affected the rate of transitions both towards and away from gynodioecy; transitions towards gynodioecy occurred ˜16 times more frequently and transitions away from gynodioecy occurred ˜11 times more frequently in herbaceous lineages than in woody lineages. Within the Lamiaceae, lineages in which gynodioecy has frequently evolved also have a high rate of reversions to the nongynodioecious state. Consequently, to understand why gynodioecy is rare, we need to understand why sexual systems are more evolutionarily labile in some lineages than in others.

Frequency-dependent fitness in gynodioecious Lobelia siphilitica.

Selection is frequency dependent when an individual's fitness depends on the frequency of its phenotype. Frequency-dependent selection should be common in gynodioecious plants, where individuals are female or hermaphroditic; if the fitness of females is limited by the availability of pollen to fertilize their ovules, then they should have higher fitness when rare than when common. To test whether the fitness of females is frequency dependent, we manipulated the sex ratio in arrays of gynodioecious Lobelia siphilitica. To test whether fitness was frequency dependent because of variation in pollen availability, we compared open-pollinated and supplemental hand-pollinated plants. Open-pollinated females produced more seeds when they were rare than when they were common, as expected if fitness is negatively frequency dependent. However, hand-pollinated females also produced more seeds when they were rare, indicating that variation in pollen availability was not the cause of frequency-dependent fitness. Instead, fitness was frequency dependent because both hand- and open-pollinated females opened more flowers when they were rare than when they were common. This plasticity in the rate of anthesis could cause fitness to be frequency dependent even when reproduction is not pollen limited, and thus expand the conditions under which frequency-dependent selection operates in gynodioecious species.

Conspecific plant-soil feedback scales with population size in Lobelia siphilitica (Lobeliaceae).

Plant-soil interactions directly affect plant success in terms of establishment, survival, growth and reproduction. Negative plant-soil feedback on such traits may therefore reduce the density and abundance of plants of a given species at a given site. Furthermore, if conspecific feedback varies among population sites, it could help explain geographic variation in plant population size. We tested for among-site variation in conspecific plant-soil feedback in a greenhouse experiment using seeds and soils from 8 natural populations of Lobelia siphilitica hosting 30-330 plants. The first cohort of seeds was grown on soil collected from each native site, while the second cohort was grown on the soil conditioned by the first. Our goal was to distinguish site-specific effects mediated by biotic and/or abiotic soil properties from those inherent in seed sources. Cohort 1 plants grown from seeds produced in small populations performed better in terms of germination, growth, and survival compared to plants produced in large populations. Plant performance decreased substantially between cohorts, indicating strong negative feedback. Most importantly, the strength of negative feedback scaled linearly (i.e., was less negative) with increasing size of the native plant population, particularly for germination and survival, and was better explained by soil- rather than seed-source effects. Even with a small number of sites, our results suggest that the potential for negative plant-soil feedback varies among populations of L. siphilitica, and that small populations were more susceptible to negative feedback. Conspecific plant-soil feedback may contribute to plant population size variation within a species' native range.

Testing models of sex ratio evolution in a gynodioecious plant: female frequency covaries with the cost of male fertility restoration.

In many gynodioecious species, cytoplasmic male sterility genes (CMS) and nuclear male fertility restorers (Rf) jointly determine whether a plant is female or hermaphrodite. Equilibrium models of cytonuclear gynodioecy, which describe the effect of natural selection within populations on the sex ratio, predict that the frequency of females in a population will primarily depend on the cost of male fertility restoration, a negative pleiotropic effect of Rf alleles on hermaphrodite fitness. Specifically, when the cost of restoration is higher, the frequency of females at equilibrium is predicted to be higher. To test this prediction, we estimated variation in the cost of restoration across 26 populations of Lobelia siphilitica, a species in which Rf alleles can have negative pleiotropic effects on pollen viability. We found that L. siphilitica populations with many females were more likely to contain hermaphrodites with low pollen viability. This is consistent with the prediction that the cost of restoration is a key determinant of variation in female frequency. Our results suggest that equilibrium models can explain variation in sex ratio among natural populations of gynodioecious species.

The evolutionary ecology of cytonuclear interactions in angiosperms.

Interactions between cytoplasmic and nuclear genomes have significant evolutionary consequences. In angiosperms, the most common cytonuclear interaction is between mitochondrial genes that disrupt pollen production (cytoplasmic male sterility, CMS) and nuclear genes that restore it (nuclear male fertility restorers, Rf). The outcome of CMS/Rf interactions can depend on whether Rf alleles have negative pleiotropic effects on fitness. Although these fitness costs are often considered to be independent of the ecological context, we argue that the effects of Rf alleles on fitness should be context dependent. Thus, measuring the cost of restoration across a range of environments could help explain geographic and phylogenetic variation in the distribution of Rf alleles and the outcome of CMS/Rf interactions.

Evidence against equimolarity of large repeat arrangements and a predominant master circle structure of the mitochondrial genome from a monkeyflower (Mimulus guttatus) lineage with cryptic CMS.

Despite intense investigation for over 25 years, the in vivo structure of plant mitochondrial genomes remains uncertain. Mapping studies and genome sequencing generally produce large circular chromosomes, whereas electrophoretic and microscopic studies typically reveal linear and multibranched molecules. To more fully assess the structure of plant mitochondrial genomes, the complete sequence of the monkeyflower (Mimulus guttatus DC. line IM62) mitochondrial DNA was constructed from a large (35 kb) paired-end shotgun sequencing library to a high depth of coverage (~30×). The complete genome maps as a 525,671 bp circular molecule and exhibits a fairly conventional set of features including 62 genes (encoding 35 proteins, 24 transfer RNAs, and 3 ribosomal RNAs), 22 introns, 3 large repeats (2.7, 9.6, and 29 kb), and 96 small repeats (40-293 bp). Most paired-end reads (71%) mapped to the consensus sequence at the expected distance and orientation across the entire genome, validating the accuracy of assembly. Another 10% of reads provided clear evidence of alternative genomic conformations due to apparent rearrangements across large repeats. Quantitative assessment of these repeat-spanning read pairs revealed that all large repeat arrangements are present at appreciable frequencies in vivo, although not always in equimolar amounts. The observed stoichiometric differences for some arrangements are inconsistent with a predominant master circular structure for the mitochondrial genome of M. guttatus IM62. Finally, because IM62 contains a cryptic cytoplasmic male sterility (CMS) system, an in silico search for potential CMS genes was undertaken. The three chimeric open reading frames (ORFs) identified in this study, in addition to the previously identified ORFs upstream of the nad6 gene, are the most likely CMS candidate genes in this line.

New perspectives on the evolution of plant mating systems.

BACKGROUND: The remarkable diversity of mating patterns and sexual systems in flowering plants has fascinated evolutionary biologists for more than a century. Enduring questions about this topic include why sexual polymorphisms have evolved independently in over 100 plant families, and why proportions of self- and cross-fertilization often vary dramatically within and among populations. Important new insights concerning the evolutionary dynamics of plant mating systems have built upon a strong foundation of theoretical models and innovative field and laboratory experiments. However, as the pace of advancement in this field has accelerated, it has become increasingly difficult for researchers to follow developments outside their primary area of research expertise. SCOPE: In this Viewpoint paper we highlight three important themes that span and integrate different subdisciplines: the changes in morphology, phenology, and physiology that accompany the transition to selfing; the evolutionary consequences of pollen pool diversity in flowering plants; and the evolutionary dynamics of sexual polymorphisms. We also highlight recent developments in molecular techniques that will facilitate more efficient and cost-effective study of mating patterns in large natural populations, research on the dynamics of pollen transport, and investigations on the genetic basis of sexual polymorphisms. This Viewpoint also serves as the introduction to a Special Issue on the Evolution of Plant Mating Systems. The 15 papers in this special issue provide inspiring examples of recent discoveries, and glimpses of exciting developments yet to come.

Metacommunity organization of soil microorganisms depends on habitat defined by presence of Lobelia siphilitica plants.

We tested regional-scale spatial patterns in soil microbial community composition for agreement with species sorting and dispersal limitation, two alternative mechanisms behind different models of metacommunity organization. Furthermore, we tested whether regional metacommunity organization depends on local habitat type. We sampled from sites across Ohio and West Virginia hosting populations of Lobelia siphilitica, and compared the metacommunity organization of soil microbial communities under L. siphilitica to those in adjacent areas at each site. In the absence of L. siphilitica, bacterial community composition across the region was consistent with species sorting. However, under L. siphilitica, bacterial community composition was consistent with dispersal limitation. Fungal community composition remained largely unexplained, although fungal communities under L. siphilitica were both significantly different in composition and less variable in composition than in adjacent areas. Our results show that communities in different local habitat types (e.g., in the presence or absence of a particular plant) may be structured on a regional scale by different processes, despite being separated by only centimeters at the local scale.

A novel approach to estimating the cost of male fertility restoration in gynodioecious plants.

*In many gynodioecious plants, sex is determined by cytoplasmic male sterility genes (CMS) and nuclear male fertility restorers (Rf). Models predict that the costs of restoration are important determinants of population sex ratios. However, current approaches to the estimation of these costs require prior identification of CMS genotypes, information that is available for few species. *We tested a novel approach to estimating the cost of restoration in natural populations without determining CMS or Rf genotypes. We used estimates of pollen viability and offspring sex ratios from open- and hand-pollinated families of Lobelia siphilitica to test whether the cost of restoration, expressed as low pollen viability, is higher in populations with more females. *Among populations with CMS, we found that variation in pollen viability was higher in small populations with more females, as expected if the proportion of females within populations increases with the maximum cost of restoration. In controlled crosses, families with low pollen viability also produced fewer females, suggesting that variation in viability is primarily determined by the number and frequency of Rf alleles carried. *This approach to estimating the cost of restoration can be applied to other cytonuclear gynodioecious species, offering new opportunities for testing gynodioecy models in the wild.

Hybrid male sterility in Mimulus (Phrymaceae) is associated with a geographically restricted mitochondrial rearrangement.

Cytoplasmic male sterility (CMS) and nuclear fertility restoration (Rf) involves intergenomic coevolution. Although male-sterile phenotypes are rarely expressed in natural populations of angiosperms, CMS genes are thought to be common. The evolutionary dynamics of CMS/Rf systems are poorly understood, leaving gaps in our understanding of mechanisms and consequences of cytonuclear interactions. We characterized the molecular basis and geographic distribution of a CMS gene in Mimulus guttatus. We used outcrossing M. guttatus (with CMS and Rf) to self-fertilizing M. nasutus (lacking CMS and Rf) to generate hybrids segregating for CMS. Mitochondrial transcripts containing an essential gene (nad6) were perfectly associated with male sterility. The CMS mitotype was completely absent in M. nasutus, present in all genotypes collected from the original collection site, but in only two individuals from 34 other M. guttatus populations. This pattern suggests that the CMS likely originated at a single locality, spread to fixation within the population, but has not spread to other populations, indicating possible ecological or genetic constraints on dispersal of this CMS mitotype between populations. Extreme localization may be characteristic of CMS in hermaphroditic species, in contrast to geographically widespread mitotypes commonly found in gynodioecious species, and could directly contribute to hybrid incompatibilities in nature.

An experimental test of the effects of resources and sex ratio on maternal fitness and phenotypic selection in gynodioecious Fragaria virginiana.

Resources, sex ratio, and seed production by hermaphrodites covary among natural populations of many gynodioecious plant species, such that they are functionally "more dioecious" as resources become more limiting. Strong correlations among these three factors confound our understanding of their relative roles in maintaining polymorphic sexual systems. We manipulated resource availability and sex ratio and measured their effects on relative fertility and phenotypic selection through the maternal fitness of females and hermaphrodites of Fragaria virginiana. Two results were particularly surprising. First, hermaphrodites showed little variability in fecundity across resource treatments and showed strong positive and context-dependent selection for fruit set. This suggests that variation in hermaphrodite seed production along resource gradients in nature may result from adaptation rather than plasticity. Second, although females increased their fecundity with higher resources, their fertility was unaffected by sex ratio, which is predicted to mediate pollen limitation of females in natural populations where they are common. Selection on petal size of females was also weak, indicating a minimal effect of pollinator attraction on variation in the fertility of female plants. Hence, we found no mechanistic explanation for the complete absence of high-resource high female populations in nature. Despite strong selection for increased fruit set of hermaphrodites, both the strength of selection and its contribution to the maintenance of gynodioecy are severely reduced under conditions where females have high relative fecundity (i.e., low resources and high-female sex ratios). High relative fertility plus high female frequency means that the evolution of phenotypic traits in hermaphrodites (i.e., response to selection via seed function) should be manifested through females because most hermaphrodites will have female mothers. Fruit set was never under strong selection in females; hence, selection to increase fruit set hermaphrodites will be less effective in maintaining their fruiting ability in natural populations with low resources and high female frequency. In sum, both sex ratio and resource availability influence trait evolution indirectly-through their effects on relative fertility of the sexes and patterns of selection. Sex ratio did not impose strong pollen limitation on females but did directly moderate the outcome of natural selection by biasing the maternal sex of the next generation. This direct effect of sex ratio on the manifestation of natural selection is expected to have far greater impact on the evolution of traits, such as seed-producing ability in hermaphrodites and the maintenance of sexual polymorphisms in nature, compared to indirect effects of sex ratio on relative fertility of the sexes.