Genomic studies on the nature of species: adaptation and speciation inMimulus

Suillus: an emerging model for the study of ectomycorrhizal ecology and evolution

Research on mycorrhizal symbiosis has been slowed by a lack of established study systems. To address this challenge, we have been developing Suillus, a widespread ecologically and economically relevant fungal genus primarily associated with the plant family Pinaceae, into a model system for studying ectomycorrhizal (ECM) associations. Over the last decade, we have compiled extensive genomic resources, culture libraries, a phenotype database, and protocols for manipulating Suillus fungi with and without their tree partners. Our efforts have already resulted in a large number of publicly available genomes, transcriptomes, and respective annotations, as well as advances in our understanding of mycorrhizal partner specificity and host communication, fungal and plant nutrition, environmental adaptation, soil nutrient cycling, interspecific competition, and biological invasions. Here, we highlight the most significant recent findings enabled by Suillus, present a suite of protocols for working with the genus, and discuss how Suillus is emerging as an important model to elucidate the ecology and evolution of ECM interactions.

Genomic Approaches Are Improving Taxonomic Representation in Genetic Studies of Speciation

Until recently, our understanding of the genetics of speciation was limited to a narrow group of model species with a specific set of characteristics that made genetic analysis feasible. Rapidly advancing genomic technologies are eliminating many of the distinctions between laboratory and natural systems. In light of these genomic developments, we review the history of speciation genetics, advances that have been gleaned from model and non-model organisms, the current state of the field, and prospects for broadening the diversity of taxa included in future studies. Responses to a survey of speciation scientists across the world reveal the ongoing division between the types of questions that are addressed in model and non-model organisms. To bridge this gap, we suggest integrating genetic studies from model systems that can be reared in the laboratory or greenhouse with genomic studies in related non-models where extensive ecological knowledge exists.

Genomic hotspots of chromosome rearrangements explain conserved synteny despite high rates of chromosome evolution in a holocentric lineage

Holocentric organisms, unlike typical monocentric organisms, have kinetochore activity distributed along almost the whole length of the chromosome. Because of this, chromosome rearrangements through fission and fusion are more likely to become fixed in holocentric species, which may account for the extraordinary rates of chromosome evolution that many holocentric lineages exhibit. Long blocks of genome synteny have been reported in animals with holocentric chromosomes despite high rates of chromosome rearrangements. Nothing is known from plants, however, despite the fact that holocentricity appears to have played a key role in the diversification of one of the largest angiosperm genera, Carex (Cyperaceae). In the current study, we compared genomes of Carex species and a distantly related Cyperaceae species to characterize conserved and rearranged genome regions. Our analyses span divergence times ranging between 2 and 50 million years. We also compared a C. scoparia chromosome-level genome assembly with a linkage map of the same species to study rearrangements at a population level and suppression of recombination patterns. We found longer genome synteny blocks than expected under a null model of random rearrangement breakpoints, even between very distantly related species. We also found repetitive DNA to be non-randomly associated with holocentromeres and rearranged regions of the genome. The evidence of conserved synteny in sedges despite high rates of chromosome fission and fusion suggests that conserved genomic hotspots of chromosome evolution related to repetitive DNA shape the evolution of recombination, gene order and crossability in sedges. This finding may help explain why sedges are able to maintain species cohesion even in the face of high interspecific chromosome rearrangements.

Speciation and development

Understanding general principles about the origin of species remains one of the foundational challenges in evolutionary biology. The genomic divergence between groups of individuals can spawn hybrid inviability and hybrid sterility, which presents a tantalizing developmental problem. Divergent developmental programs may yield either conserved or divergent phenotypes relative to ancestral traits, both of which can be responsible for reproductive isolation during the speciation process. The genetic mechanisms of developmental evolution involve cis- and trans-acting gene regulatory change, protein-protein interactions, genetic network structures, dosage, and epigenetic regulation, all of which also have roots in population genetic and molecular evolutionary processes. Toward the goal of demystifying Darwin's "mystery of mysteries," this review integrates microevolutionary concepts of genetic change with principles of organismal development, establishing explicit links between population genetic process and developmental mechanisms in the production of macroevolutionary pattern. This integration aims to establish a more unified view of speciation that binds process and mechanism.

Phylogenomic analyses in Phrymaceae reveal extensive gene tree discordance in relationships among major clades

PREMISE

Phylogenomic datasets using genomes and transcriptomes provide rich opportunities beyond resolving bifurcating phylogenetic relationships. Monkeyflower (Phrymaceae) is a model system for evolutionary ecology. However, it lacks a well-supported phylogeny as a basis for a stable taxonomy and for macroevolutionary comparisons.

METHODS

We sampled 24 genomes and transcriptomes in Phrymaceae and closely related families, including eight newly sequenced transcriptomes. We reconstructed the phylogeny using IQ-TREE and ASTRAL, evaluated gene tree discordance using PhyParts, Quartet Sampling, and a cloudogram, and carried out reticulation analyses using PhyloNet and HyDe. We searched for whole genome duplication (WGD) events using chromosome numbers, synonymous distances, and gene duplication events as evidence.

RESULTS

Most gene trees support the monophyly of Phrymaceae and each of its tribes. Most gene trees also support tribe Mimuleae being sister to Phrymeae + Diplaceae + Leucocarpeae, with extensive gene tree discordance among the latter three. Despite the discordance, the monophyly of Mimulus s.l. is rejected, and no individual reticulation event among the Phrymaceae tribes is well-supported. Reticulation likely occurred among Erythranthe bicolor and closely related species. No ancient WGD was detected in Phrymaceae. Instead, small-scale duplications are among potential drivers of macroevolutionary diversification of Phrymaceae.

CONCLUSIONS

We show that analysis of reticulate evolution is sensitive to taxon sampling and methods used. We also demonstrate that phylogenomic datasets using genomes and transcriptomes present rich opportunities to investigate gene family evolution and genome duplication events involved in lineage diversification and adaptation.

The regulatory network for petal anthocyanin pigmentation is shaped by the MYB5a/NEGAN transcription factor in Mimulus

Much of the visual diversity of angiosperms is due to the frequent evolution of novel pigmentation patterns in flowers. The gene network responsible for anthocyanin pigmentation, in particular, has become a model for investigating how genetic changes give rise to phenotypic innovation. In the monkeyflower genus Mimulus, an evolutionarily recent gain of petal lobe anthocyanin pigmentation in M. luteus var. variegatus was previously mapped to genomic region pla2. Here, we use sequence and expression analysis, followed by transgenic manipulation of gene expression, to identify MYB5a—orthologous to the NEGAN transcriptional activator from M. lewisii—as the gene responsible for the transition to anthocyanin-pigmented petals in M. l. variegatus. In other monkeyflower taxa, MYB5a/NEGAN is part of a reaction-diffusion network that produces semi-repeating spotting patterns, such as the array of spots in the nectar guides of both M. lewisii and M. guttatus. Its co-option for the evolution of an apparently non-patterned trait—the solid petal lobe pigmentation of M. l. variegatus—illustrates how reaction-diffusion can contribute to evolutionary novelty in non-obvious ways. Transcriptome sequencing of a MYB5a RNAi line of M. l. variegatus reveals that this genetically simple change, which we hypothesize to be a regulatory mutation in cis to MYB5a, has cascading effects on gene expression, not only on the enzyme-encoding genes traditionally thought of as the targets of MYB5a but also on all of its known partners in the anthocyanin regulatory network.

Digest: The role of postzygotic isolation in Mimulus speciation

Recently diverged species often show incomplete reproductive isolation and subsequently experience hybridization and introgression. The plant genus Mimulus includes many such examples of closely related taxa that prove useful for studying incipient speciation. However, Sandstedt et al. show that in contrast to other Mimulus species, species of the M. tilingii complex are characterized by strong postzygotic isolation mediated by multiple barriers. This finding highlights the diverse speciation histories of related plant groups.

Multi-level patterns of genetic structure and isolation by distance in the widespread plant Mimulus guttatus

An understanding of genetic structure is essential for answering many questions in population genetics. However, complex population dynamics and scale-dependent processes can make it difficult to detect if there are distinct genetic clusters present in natural populations. Inferring discrete population structure is particularly challenging in the presence of continuous genetic variation such as isolation by distance. Here, we use the plant species Mimulus guttatus as a case study for understanding genetic structure at three spatial scales. We use reduced-representation sequencing and marker-based genotyping to understand dispersal dynamics and to characterise genetic structure. Our results provide insight into the spatial scale of genetic structure in a widespread plant species, and demonstrate how dispersal affects spatial genetic variation at the local, regional, and range-wide scale. At a fine-spatial scale, we show dispersal is rampant with little evidence of spatial genetic structure within populations. At a regional-scale, we show continuous differentiation driven by isolation by distance over hundreds of kilometres, with broad geographic genetic clusters that span major barriers to dispersal. Across Western North America, we observe geographic genetic structure and the genetic signature of multiple postglacial recolonisation events, with historical gene flow linking isolated populations. Our genetic analyses show M. guttatus is highly dispersive and maintains large metapopulations with high intrapopulation variation. This high diversity and dispersal confounds the inference of genetic structure, with multi-level sampling and spatially-explicit analyses required to understand population history.

Afro-alpine flagships revisited: Parallel adaptation, intermountain admixture and shallow genetic structuring in the giant senecios (Dendrosenecio)

Distantly related lineages of the enigmatic giant rosette plants of tropical alpine environments provide classical examples of convergent adaptation. For the giant senecios (Dendrosenecio), the endemic landmarks of the East African sky islands, it has also been suggested that parallel adaptation has been important for within-lineage differentiation. To test this hypothesis and to address potential gene flow and hybridization among the isolated sky islands, we organized field expeditions to all major mountains. We sampled all currently accepted species and all but one subspecies and genotyped 460 plants representing 109 populations. We tested whether genetic structuring corresponds to geography, as predicted by a parallel adaptation hypothesis, or to altitudinal belt and habitat rather than mountains, as predicted by a hypothesis of a single origin of adaptations. Bayesian and Neighbor-Net analyses showed that the main genetic structure is shallow and largely corresponds to geography, supporting a hypothesis of recent, rapid radiation via parallel altitude/habitat adaptation on different mountains. We also found evidence for intermountain admixture, suggesting several long-distance dispersals by wind across vast areas of unsuitable habitat. The combination of parallel adaptation, secondary contact, and hybridization may explain the complex patterns of morphological variation and the contradicting taxonomic treatments of these rare enigmatic giants, supporting the use of wide taxonomic concepts. Notably, the within-population genetic diversity was very low and calls for increased conservation efforts.

Patterns of Hybrid Seed Inviability in the Mimulus guttatus sp. Complex Reveal a Potential Role of Parental Conflict in Reproductive Isolation

Genomic conflicts may play a central role in the evolution of reproductive barriers. Theory predicts that early-onset hybrid inviability may stem from conflict between parents for resource allocation to offspring. Here, we describe M. decorus: a group of cryptic species within the M. guttatus species complex that are largely reproductively isolated by hybrid seed inviability (HSI). HSI between M. guttatus and M. decorus is common and strong, but populations of M. decorus vary in the magnitude and directionality of HSI with M. guttatus. Patterns of HSI between M. guttatus and M. decorus, as well as within M. decorus, conform to the predictions of parental conflict: first, reciprocal F1s exhibit size differences and parent-of-origin-specific endosperm defects; second, the extent of asymmetry between reciprocal F1 seed size is correlated with asymmetry in HSI; and third, inferred differences in the extent of conflict predict the extent of HSI between populations. We also find that HSI is rapidly evolving, as populations that exhibit the most HSI are each others' closest relative. Lastly, although all populations appear largely outcrossing, we find that the differences in the inferred strength of conflict scale positively with π, suggesting that demographic or life history factors other than transitions to self-fertilization may influence the rate of parental-conflict-driven evolution. Overall, these patterns suggest the rapid evolution of parent-of-origin-specific resource allocation alleles coincident with HSI within and between M. guttatus and M. decorus. Parental conflict may therefore be an important evolutionary driver of reproductive isolation.

Ancestral Admixture Is the Main Determinant of Global Biodiversity in Fission Yeast

Mutation and recombination are key evolutionary processes governing phenotypic variation and reproductive isolation. We here demonstrate that biodiversity within all globally known strains of Schizosaccharomyces pombe arose through admixture between two divergent ancestral lineages. Initial hybridization was inferred to have occurred ∼20-60 sexual outcrossing generations ago consistent with recent, human-induced migration at the onset of intensified transcontinental trade. Species-wide heritable phenotypic variation was explained near-exclusively by strain-specific arrangements of alternating ancestry components with evidence for transgressive segregation. Reproductive compatibility between strains was likewise predicted by the degree of shared ancestry. To assess the genetic determinants of ancestry block distribution across the genome, we characterized the type, frequency, and position of structural genomic variation using nanopore and single-molecule real-time sequencing. Despite being associated with double-strand break initiation points, over 800 segregating structural variants exerted overall little influence on the introgression landscape or on reproductive compatibility between strains. In contrast, we found strong ancestry disequilibrium consistent with negative epistatic selection shaping genomic ancestry combinations during the course of hybridization. This study provides a detailed, experimentally tractable example that genomes of natural populations are mosaics reflecting different evolutionary histories. Exploiting genome-wide heterogeneity in the history of ancestral recombination and lineage-specific mutations sheds new light on the population history of S. pombe and highlights the importance of hybridization as a creative force in generating biodiversity.

Gene regulatory divergence between locally adapted ecotypes in their native habitats

Local adaptation is a key driver of ecological specialization and the formation of new species. Despite its importance, the evolution of gene regulatory divergence among locally adapted populations is poorly understood, especially how that divergence manifests in nature. Here, we evaluate gene expression divergence and allele-specific gene expression responses for locally adapted coastal perennial and inland annual accessions of the yellow monkeyflower, Mimulus guttatus, in a field reciprocal transplant experiment. Overall, 6765 (73%) of surveyed genes were differentially expressed between coastal and inland habitats, while 7213 (77%) were differentially expressed between the coastal perennial and inland annual accessions. Cis-regulatory variation was pervasive, affecting 79% (5532) of differentially expressed genes. We detected trans effects for 52% (3611) of differentially expressed genes. Expression plasticity of alleles across habitats (G × E interactions) appears to be relatively common (affecting 18% of transcripts) and could minimize fitness trade-offs at loci that contribute to local adaptation. We also found evidence that at least one chromosomal inversion may act as supergene by holding together haplotypes of differentially expressed genes, but this pattern depends on habitat context. Our results highlight multiple key patterns regarding the relationship between gene expression and the evolution of locally adapted populations.

Notes from a snail island: Littorinid evolution and adaptation

The most successful study systems are built on a foundation of decades of research on the basic biology, ecology and life history of the organisms in question. Combined with new technologies, this can provide a formidable means to address important issues in evolutionary biology and molecular ecology. Littorinid marine snails are a good example of this, with a rich literature on their taxonomy, speciation, thermal tolerance and behavioural adaptations. In August 2017, an international meeting on Littorinid evolution was held at the Tjärnö Marine Research Laboratory in Western Sweden. In this meeting review, I provide a summary of some of the exciting work on parallel evolution, sexual selection and adaptation to environmental stress presented there. I argue that newly available genomic resources present an opportunity for integrating the traditionally divergent fields of speciation and environmental adaptation in Littorinid research.

A Whole-Transcriptome Approach to Evaluating Reference Genes for Quantitative Gene Expression Studies: A Case Study in Mimulus

While quantitative PCR (qPCR) is widely recognized as being among the most accurate methods for quantifying gene expression, it is highly dependent on the use of reliable, stably expressed reference genes. With the increased availability of high-throughput methods for measuring gene expression, whole-transcriptome approaches may be increasingly utilized for reference gene selection and validation. In this study, RNA-seq was used to identify a set of novel qPCR reference genes and evaluate a panel of traditional “housekeeping” reference genes in two species of the evolutionary model plant genus Mimulus. More broadly, the methods proposed in this study can be used to harness the power of transcriptomes to identify appropriate reference genes for qPCR in any study organism, including emerging and nonmodel systems. We find that RNA-seq accurately estimates gene expression means in comparison to qPCR, and that expression means are robust to moderate environmental and genetic variation. However, measures of expression variability were only in agreement with qPCR for samples obtained from a shared environment. This result, along with transcriptome-wide comparisons, suggests that environmental changes have greater impacts on expression variability than on expression means. We discuss how this issue can be addressed through experimental design, and suggest that the ever-expanding pool of published transcriptomes represents a rich and low-cost resource for developing better reference genes for qPCR.

The genetics of extreme microgeographic adaptation: an integrated approach identifies a major gene underlying leaf trichome divergence in Yellowstone Mimulus guttatus

Microgeographic adaptation provides a particularly interesting context for understanding the genetic basis of phenotypic divergence and may also present unique empirical challenges. In particular, plant adaptation to extreme soil mosaics may generate barriers to gene flow or shifts in mating system that confound simple genomic scans for adaptive loci. Here, we combine three approaches - quantitative trait locus (QTL) mapping of candidate intervals in controlled crosses, population resequencing (PoolSeq) and analyses of wild recombinant individuals - to investigate one trait associated with Mimulus guttatus (yellow monkeyflower) adaptation to geothermal soils in Yellowstone National Park. We mapped a major QTL causing dense leaf trichomes in thermally adapted plants to a <50-kb region of linkage Group 14 (Tr14) previously implicated in trichome divergence between independent M. guttatus populations. A PoolSeq scan of Tr14 region revealed a cluster of six genes, coincident with the inferred QTL peak, with high allele frequency differences sufficient to explain observed phenotypic differentiation. One of these, the R2R3 MYB transcription factor Migut.N02661, is a plausible functional candidate and was also strongly associated (r²  = 0.27) with trichome phenotype in analyses of wild-collected admixed individuals. Although functional analyses will be necessary to definitively link molecular variants in Tr14 with trichome divergence, our analyses are a major step in that direction. They point to a simple, and parallel, genetic basis for one axis of Mimulus guttatus adaptation to an extreme habitat, suggest a broadly conserved genetic basis for trichome variation across flowering plants and pave the way for further investigations of this challenging case of microgeographic incipient speciation.

Transcriptome sequencing and simple sequence repeat marker development for three Macaronesian endemic plant species

PREMISE OF THE STUDY

Oceanic islands offer unparalleled opportunities to investigate evolutionary processes such as adaptation and speciation. However, few genomic resources are available for oceanic island endemics. In this study, we publish transcriptome sequences from three Macaronesian endemic plant species (Argyranthemum broussonetii [Asteraceae], Descurainia bourgaeana [Brassicaceae], and Echium wildpretii [Boraginaceae]) that are representative of lineages that have radiated in the region. In addition, the utility of transcriptome data for marker development is demonstrated.

METHODS AND RESULTS

Transcriptomes from the three plant species were sequenced, assembled, and annotated. Between 1972 and 2282 simple sequence repeats (SSRs) were identified for each taxon. Primers were designed and tested for 30 of the candidate SSRs identified in Argyranthemum, of which 12 amplified well across three species and eight were polymorphic.

CONCLUSIONS

We demonstrate here that a single transcriptome sequence is sufficient to identify hundreds of polymorphic SSR markers. The SSRs are applicable to a wide range of questions relating to the evolution of island lineages.

Strongly asymmetric hybridization barriers shape the origin of a new polyploid species and its hybrid ancestor

PREMISE OF THE STUDY

Hybridization between diploids and tetraploids can lead to new allopolyploid species, often via a triploid intermediate. Viable triploids are often produced asymmetrically, with greater success observed for "maternal-excess" crosses where the mother has a higher ploidy than the father. Here we investigated the evolutionary origins of Mimulus peregrinus, an allohexaploid recently derived from the triploid M. ×robertsii, to determine whether reproductive asymmetry has shaped the formation of this new species.

METHODS

We used reciprocal crosses between the diploid (M. guttatus) and tetraploid (M. luteus) progenitors to determine the viability of triploid M. ×robertsii hybrids resulting from paternal- vs. maternal-excess crosses. To investigate whether experimental results predict patterns seen in the field, we performed parentage analyses comparing natural populations of M. peregrinus to its diploid, tetraploid, and triploid progenitors. Organellar sequences obtained from pre-existing genomic data, supplemented with additional genotyping was used to establish the maternal ancestry of multiple M. peregrinus and M. ×robertsii populations.

KEY RESULTS

We found strong evidence for asymmetric origins of M. peregrinus, but opposite to the common pattern, with paternal-excess crosses significantly more successful than maternal-excess crosses. These results successfully predicted hybrid formation in nature: 111 of 114 M. ×robertsii individuals, and 27 of 27 M. peregrinus, had an M. guttatus maternal haplotype.

CONCLUSION

This study, which includes the first Mimulus chloroplast genome assembly, demonstrates the utility of parentage analysis through genome skimming. We highlight the benefits of complementing genomic analyses with experimental approaches to understand asymmetry in allopolyploid speciation.

Absence of postmating barriers between a selfing vs. outcrossing Chilean Mimulus species pair

PREMISE OF THE STUDY

Reproductive isolation between sympatric species pairs may be maintained by both pre- and postmating barriers. Here we evaluate potential barriers to mating between the outcrossing Mimulus luteus and its more highly selfing sympatric congener, M. cupreus, two members of the South American luteus complex of Mimulus.

METHODS

Seed set was compared following autonomous self-pollination, manual pollination, conspecific outcrossing, and sympatric and allopatric hybridization, for laboratory-maintained inbred lines and wild-collected accessions. Survival and reproductive fitness of hybrids relative to parental species were examined across environments that differed with respect to temperature and soil nutrients, two factors that vary across the ranges of M. luteus and M. cupreus.

KEY RESULTS

Mimulus luteus was minimally capable of autonomous self-fertilization, consistent with reliance on an animal pollinator, whereas M. cupreus was a successful selfer across all tested accessions. Postmating barriers to hybridization are negligible, in both low- and high-stress environments, across multiple sympatric and allopatric populations.

CONCLUSION

As in the North American M. guttatus-M. nasutus species pair, postmating barriers contribute little to isolation between M. luteus and M. cupreus. This result reinforces the importance of premating barriers, specifically species differences in reliance on, and accessibility to, animal pollinators. A unique aspect of the M. luteus-M. cupreus pair is the recent gain of red floral anthocyanin pigmentation in M. cupreus. On the basis of species differences in vegetative anthocyanin production, a facultative stress-protective response, we propose a potential stress-protective role for the constitutive floral anthocyanins of M. cupreus.

Comparative Transcriptomics Indicates a Role for SHORT VEGETATIVE PHASE (SVP) Genes in Mimulus guttatus Vernalization Response

The timing of reproduction in response to variable environmental conditions is critical to plant fitness, and is a major driver of taxon differentiation. In the yellow monkey flower, Mimulus guttatus, geographically distinct North American populations vary in their photoperiod and chilling (vernalization) requirements for flowering, suggesting strong local adaptation to their surroundings. Previous analyses revealed quantitative trait loci (QTL) underlying short-day mediated vernalization responsiveness using two annual M. guttatus populations that differed in their vernalization response. To narrow down candidate genes responsible for this variation, and to reveal potential downstream genes, we conducted comparative transcriptomics and quantitative PCR (qPCR) in shoot apices of parental vernalization responsive IM62, and unresponsive LMC24 inbred lines grown under different photoperiods and temperatures. Our study identified several metabolic, hormone signaling, photosynthetic, stress response, and flowering time genes that are differentially expressed between treatments, suggesting a role for their protein products in short-day-mediated vernalization responsiveness. Only a small subset of these genes intersected with candidate genes from the previous QTL study, and, of the main candidates tested with qPCR under nonpermissive conditions, only SHORT VEGETATIVE PHASE (SVP) gene expression met predictions for a population-specific short-day-repressor of flowering that is repressed by cold.