Genetic factors predict hybrid formation in the British flora.

Natural hybridization can have a profound evolutionary impact, with consequences ranging from the extinction of rare taxa to the origin of new species. Natural hybridization is particularly common in plants; however, our understanding of the general factors that promote or prevent hybridization is hampered by the highly variable outcomes in different lineages. Here, we quantify the influence of different predictors on hybrid formation across species from an entire flora. We combine estimates of hybridization with ecological attributes and a new species-level phylogeny for over 1,100 UK flowering plant species. Our results show that genetic factors, particularly parental genetic distance, as well as phylogenetic position and ploidy, are key determinants of hybrid formation, whereas many other factors such as range overlap and genus size explain much less variation in hybrid formation. Overall, intrinsic genetic factors shape the evolutionary and ecological consequences of natural hybridization across species in a flora.

The emerging importance of cross-ploidy hybridisation and introgression.

Natural hybridisation is now recognised as pervasive in its occurrence across the Tree of Life. Resurgent interest in natural hybridisation fuelled by developments in genomics has led to an improved understanding of the genetic factors that promote or prevent species cross-mating. Despite this body of work overturning many widely held assumptions about the genetic barriers to hybridisation, it is still widely thought that ploidy differences between species will be an absolute barrier to hybridisation and introgression. Here, we revisit this assumption, reviewing findings from surveys of polyploidy and hybridisation in the wild. In a case study in the British flora, 203 hybrids representing 35% of hybrids with suitable data have formed via cross-ploidy matings, while a wider literature search revealed 59 studies (56 in plants and 3 in animals) in which cross-ploidy hybridisation has been confirmed with genetic data. These results show cross-ploidy hybridisation is readily overlooked, and potentially common in some groups. General findings from these studies include strong directionality of hybridisation, with introgression usually towards the higher ploidy parent, and cross-ploidy hybridisation being more likely to involve allopolyploids than autopolyploids. Evidence for adaptive introgression across a ploidy barrier and cases of cross-ploidy hybrid speciation shows the potential for important evolutionary outcomes.

Hybridization and divergent climatic preferences drive divergence of two allopatric Gentiana species on the Qinghai-Tibet Plateau.

BACKGROUND AND AIMS: Exploring how species diverge is vital for understanding the drivers of speciation. Factors such as geographical separation and ecological selection, hybridization, polyploidization and shifts in mating system are all major mechanisms of plant speciation, but their contributions to divergence are rarely well understood. Here we test these mechanisms in two plant species, Gentiana lhassica and G. hoae, with the goal of understanding recent allopatric species divergence on the Qinghai-Tibet Plateau (QTP). METHODS: We performed Bayesian clustering, phylogenetic analysis and estimates of hybridization using 561 302 nuclear genomic single nucleotide polymorphisms (SNPs). We performed redundancy analysis, and identified and annotated species-specific SNPs (ssSNPs) to explore the association between climatic preference and genetic divergence. We also estimated genome sizes using flow cytometry to test for overlooked polyploidy. KEY RESULTS: Genomic evidence confirms that G. lhassica and G. hoae are closely related but distinct species, while genome size estimates show divergence occurred without polyploidy. Gentiana hoae has significantly higher average FIS values than G. lhassica. Population clustering based on genomic SNPs shows no signature of recent hybridization, but each species is characterized by a distinct history of hybridization with congeners that has shaped genome-wide variation. Gentiana lhassica has captured the chloroplast and experienced introgression with a divergent gentian species, while G. hoae has experienced recurrent hybridization with related taxa. Species distribution modelling suggested range overlap in the Last Interglacial Period, while redundancy analysis showed that precipitation and temperature are the major climatic differences explaining the separation of the species. The species differ by 2993 ssSNPs, with genome annotation showing missense variants in genes involved in stress resistance. CONCLUSIONS: This study suggests that the distinctiveness of these species on the QTP is driven by a combination of hybridization, geographical isolation, mating system differences and evolution of divergent climatic preferences.

Is there hybridization between diploid and tetraploid Euphrasia in a secondary contact zone?

PREMISE: Strong postzygotic reproductive isolating barriers are usually expected to limit the extent of natural hybridization between species with contrasting ploidy. However, genomic sequencing has revealed previously overlooked examples of natural cross-ploidy hybridization in some flowering plant genera, suggesting that the phenomenon may be more common than once thought. We investigated potential cross-ploidy hybridization in British eyebrights (Euphrasia, Orobanchaceae), a group from which 13 putative cross-ploidy hybrid combinations have been reported based on morphology. METHODS: We analyzed a contact zone between diploid Euphrasia rostkoviana and tetraploid E. arctica in Wales. We sequenced part of the internal transcribed spacer (ITS) of nuclear ribosomal DNA and used genotyping by sequencing (GBS) to look for evidence of cross-ploidy hybridization and introgression. RESULTS: Common variant sites in the ITS region were fixed between diploids and tetraploids, indicating a strong barrier to hybridization. Clustering analyses of 356 single-nucleotide polymorphisms (SNPs) generated using GBS clearly separated samples by ploidy and revealed strong genetic structure (FST = 0.44). However, the FST distribution across all SNPs was bimodal, indicating potential differential selection on loci between diploids and tetraploids. Demographic inference suggested potential gene flow, limited to around one or fewer migrants per generation. CONCLUSIONS: Our results suggest that recent cross-ploidy hybridization is rare or absent in a site of secondary contact in Euphrasia. While a strong ploidy barrier prevents hybridization over ecological timescales, such hybrids may form in stable populations over evolutionary timescales, potentially allowing cross-ploidy introgression to take place.

The symmetry spectrum in a hybridising, tropical group of rhododendrons.

Many diverse plant clades possess bilaterally symmetrical flowers and specialised pollination syndromes, suggesting that these traits may promote diversification. We examined the evolution of diverse floral morphologies in a species-rich tropical radiation of Rhododendron. We used restriction-site associated DNA sequencing on 114 taxa from Rhododendron sect. Schistanthe to reconstruct phylogenetic relationships and examine hybridisation. We then captured and quantified floral variation using geometric morphometric analyses, which we interpreted in a phylogenetic context. We uncovered phylogenetic conflict and uncertainty caused by introgression within and between clades. Morphometric analyses revealed flower symmetry to be a morphological continuum without clear transitions between radial and bilateral symmetry. Tropical Rhododendron species that began diversifying into New Guinea c. 6 million years ago expanded into novel floral morphological space. Our results showed that the evolution of tropical Rhododendron is characterised by recent speciation, recurrent hybridisation and the origin of floral novelty. Floral variation evolved via changes to multiple components of the corolla that are only recognised in geometric morphometrics with both front and side views of flowers.

A phylogeny of Antirrhinum reveals parallel evolution of alpine morphology.

Parallel evolution of similar morphologies in closely related lineages provides insight into the repeatability and predictability of evolution. In the genus Antirrhinum (snapdragons), as in other plants, a suite of morphological characters are associated with adaptation to alpine environments. We tested for parallel trait evolution in Antirrhinum by investigating phylogenetic relationships using restriction-site associated DNA (RAD) sequencing. We then associated phenotypic information to our phylogeny to reconstruct the patterns of morphological evolution and related this to evidence for hybridisation between emergent lineages. Phylogenetic analyses showed that the alpine character syndrome is present in multiple groups, suggesting that Antirrhinum has repeatedly colonised alpine habitats. Dispersal to novel environments happened in the presence of intraspecific and interspecific gene flow. We found support for a model of parallel evolution in Antirrhinum. Hybridisation in natural populations, and a complex genetic architecture underlying the alpine morphology syndrome, support an important role of natural selection in maintaining species divergence in the face of gene flow.

Performance of generalist hemiparasitic Euphrasia across a phylogenetically diverse host spectrum.

Generalist hemiparasites may attach to many different host species and experience complex parasite-host interactions. How these parasite-host interactions impact on the fitness of hemiparasitic plants remain largely unknown. We used experimentally tractable eyebrights (Euphrasia, Orobanchaceae) to understand parasite-host species interactions affecting the performance of a generalist hemiparasitic plant. Common garden experiments were carried out measuring Euphrasia performance across 45 diverse hosts and in different parasite-host combinations. We show that variation in hemiparasite performance can be attributed mainly to host species and host phylogenetic relationships (λ = 0.82; 0.17-1.00 CI). When variation in performance is considered temporally, annual host species cause earlier flowering, and lead to poorer performance late in the season. While Euphrasia species typically perform similarly on a given host species, some eyebrights show more specialized parasite-host species interactions. Our results show that generalist hemiparasites only benefit from attaching to a limited, but phylogenetically divergent, subset of hosts. The conserved responses of divergent Euphrasia species suggest hemiparasite performance is affected by common host attributes. However, evidence for more complex parasite-host species interactions show that a generalist hemiparasite can potentially respond to individual host selection pressures and may adapt to local host communities.

The nature of intraspecific and interspecific genome size variation in taxonomically complex eyebrights.

BACKGROUND AND AIMS: Genome size varies considerably across the diversity of plant life. Although genome size is, by definition, affected by genetic presence/absence variants, which are ubiquitous in population sequencing studies, genome size is often treated as an intrinsic property of a species. Here, we studied intra- and interspecific genome size variation in taxonomically complex British eyebrights (Euphrasia, Orobanchaceae). Our aim is to document genome size diversity and investigate underlying evolutionary processes shaping variation between individuals, populations and species. METHODS: We generated genome size data for 192 individuals of diploid and tetraploid Euphrasia and analysed genome size variation in relation to ploidy, taxonomy, population affiliation and geography. We further compared the genomic repeat content of 30 samples. KEY RESULTS: We found considerable intraspecific genome size variation, and observed isolation-by-distance for genome size in outcrossing diploids. Tetraploid Euphrasia showed contrasting patterns, with genome size increasing with latitude in outcrossing Euphrasia arctica, but with little genome size variation in the highly selfing Euphrasia micrantha. Interspecific differences in genome size and the genomic proportions of repeat sequences were small. CONCLUSIONS: We show the utility of treating genome size as the outcome of polygenic variation. Like other types of genetic variation, such as single nucleotide polymorphisms, genome size variation may be affected by ongoing hybridization and the extent of population subdivision. In addition to selection on associated traits, genome size is predicted to be affected indirectly by selection due to pleiotropy of the underlying presence/absence variants.

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.

Life history evolution, species differences, and phenotypic plasticity in hemiparasitic eyebrights (Euphrasia).

PREMISE: Species delimitation in parasitic organisms is challenging because traits used to identify species are often plastic and vary depending on the host. Here, we use species from a recent radiation of generalist hemiparasitic Euphrasia to investigate trait variation and trait plasticity. We tested whether Euphrasia species show reliable trait differences, investigated whether these differences correspond to life history trade-offs between growth and reproduction, and quantified plasticity in response to host species. METHODS: Common garden experiments were used to evaluate trait differences between 11 Euphrasia taxa grown on a common host, document phenotypic plasticity when a single Euphrasia species is grown on eight different hosts, and relate observations to trait differences recorded in the wild. RESULTS: Euphrasia exhibited variation in life history strategies; some individuals transitioned rapidly to flowering at the expense of early season growth, while others invested in vegetative growth and delayed flowering. Life history differences were present between some species, though many related taxa lacked clear trait differences. Species differences were further blurred by phenotypic plasticity-many traits were plastic and changed with host type or between environments. CONCLUSIONS: Phenotypic plasticity in response to host and environment confounds species delimitation in Euphrasia. When grown in a common garden environment, some morphologically distinct taxa can be identified, though others represent morphologically similar shallow segregates. Trait differences present between some species and populations demonstrate the rapid evolution of distinct life history strategies in response to local ecological conditions.

Dispersal assembly of rain forest tree communities across the Amazon basin.

We investigate patterns of historical assembly of tree communities across Amazonia using a newly developed phylogeny for the species-rich neotropical tree genus Inga We compare our results with those for three other ecologically important, diverse, and abundant Amazonian tree lineages, Swartzia, Protieae, and Guatteria Our analyses using phylogenetic diversity metrics demonstrate a clear lack of geographic phylogenetic structure, and show that local communities of Inga and regional communities of all four lineages are assembled by dispersal across Amazonia. The importance of dispersal in the biogeography of Inga and other tree genera in Amazonian and Guianan rain forests suggests that speciation is not driven by vicariance, and that allopatric isolation following dispersal may be involved in the speciation process. A clear implication of these results is that over evolutionary timescales, the metacommunity for any local or regional tree community in the Amazon is the entire Amazon basin.

Loss of Color Pigmentation Is Maintained at High Frequency in a Monkey Flower Population.

Color polymorphisms have long been of evolutionary interest for their diverse roles, including mate choice, predator avoidance, and pollinator attraction. While color variation is often under strong selection, some taxa demonstrate unexpectedly high frequencies of presumed deleterious color forms. Here we show that a genetic variant underlying complete loss of anthocyanin pigmentation has risen to an unexpectedly high frequency of >0.2 in a natural population of the plant Mimulus guttatus. Decreased expression of MYB5 transcription factor is associated with unpigmented morphs. While the allele was found only in heterozygote adults in the wild, suggesting negative selection, experiments were unable to demonstrate a fitness cost for unpigmented plants, suggesting a cryptic selection pressure in the wild. However, life-history differences among morphs suggests that unpigmented individuals benefit from later flowering and clonal growth. Overall, our study highlights the complex interplay of factors maintaining variation in nature, even for genes of major effect.

New insights into the population biology of endoparasitic Rafflesiaceae.

Parasitic plants demonstrate a diversity of growth strategies, life histories, and developmental and physiological characteristics. Most research to date has focused on a narrow range of parasitic taxa, particularly in the Orobanchaceae, while the other independent origins of parasitism have largely gone unstudied. One type of parasite that has received relatively little attention are the endophytic parasites, which have a fascinating growth strategy where the parasite is embedded within the host tissue, with the flower the only externally visibly plant part. Endophytic growth makes it challenging to understand basic aspects of species biology, such as the size of a given parasite, the number of parasites per host, and the genetic diversity of populations. Recent studies by Barkman et al. (2017) and Pelser et al. (2017) have used microsatellite genotyping to investigate the population biology of endoparasitic Rafflesiaceae species in Asia. They show the potential for extensive parasite spread within a host vine and the strong partitioning of genetic diversity by host. These species are also shown to have an outcrossing mating system. However, these studies suggest different reproductive strategies, one supporting monoecy and one suggesting dioecy. Overall, these studies partly "lift the lid" on the cryptic biology of Rafflesia and the Rafflesiaceae and open the door for future comparative studies between endophytic and free-living parasitic plants.

Maintenance of species boundaries in a Neotropical radiation of Begonia.

A major goal of evolutionary biology is to determine the mechanisms generating biodiversity. In Begonia, one of the largest plant genera (1900+ species), it has been postulated that the high number of endemic species is a by-product of low gene flow among populations, which predisposes the group to speciation. However, this model of divergence requires that reproductive barriers accumulate rapidly among diverging species that overlap in their geographic ranges, otherwise speciation will be opposed by homogenizing gene flow in zones of secondary contact. Here, we test the outcomes of secondary contact in Begonia by genotyping multiple sympatric sites with 12 nuclear and seven plastid loci. We show that three sites of secondary contact between B. heracleifolia and B. nelumbiifolia are highly structured, mostly containing parental genotypes, with few F1 hybrids. A sympatric site between B. heracleifolia and B. sericoneura contains a higher proportion of F1s, but little evidence of introgression. The lack of later-generation hybrids contrasts with that documented in many other plant taxa, where introgression is extensive. Our results, in conjunction with previous genetic work, show that Begonia demonstrate properties making them exceptionally prone to speciation, at multiple stages along the divergence continuum. Not only are populations weakly connected by gene flow, promoting allopatric speciation, but species often show strong reproductive barriers in secondary contact. Whether similar mechanisms contribute to diversification in other large genera remains to be tested.

The extent and genetic basis of phenotypic divergence in life history traits in Mimulus guttatus.

Differential natural selection acting on populations in contrasting environments often results in adaptive divergence in multivariate phenotypes. Multivariate trait divergence across populations could be caused by selection on pleiotropic alleles or through many independent loci with trait-specific effects. Here, we assess patterns of association between a suite of traits contributing to life history divergence in the common monkey flower, Mimulus guttatus, and examine the genetic architecture underlying these correlations. A common garden survey of 74 populations representing annual and perennial strategies from across the native range revealed strong correlations between vegetative and reproductive traits. To determine whether these multitrait patterns arise from pleiotropic or independent loci, we mapped QTLs using an approach combining high-throughput sequencing with bulk segregant analysis on a cross between populations with divergent life histories. We find extensive pleiotropy for QTLs related to flowering time and stolon production, a key feature of the perennial strategy. Candidate genes related to axillary meristem development colocalize with the QTLs in a manner consistent with either pleiotropic or independent QTL effects. Further, these results are analogous to previous work showing pleiotropy-mediated genetic correlations within a single population of M. guttatus experiencing heterogeneous selection. Our findings of strong multivariate trait associations and pleiotropic QTLs suggest that patterns of genetic variation may determine the trajectory of adaptive divergence.

Genomic studies on the nature of species: adaptation and speciation in Mimulus.

Evolutionary biology is in an exciting era, in which powerful genomic tools make the answers accessible to long-standing questions about variation, adaptation and speciation. The availability of a suite of genomic resources, a shared knowledge base and a long history of study have made the phenotypically diverse plant genus Mimulus an important system for understanding ecological and evolutionary processes. An international Mimulus Research Meeting was held at Duke University in June 2014 to discuss developments in ecological and evolutionary genetic studies in Mimulus. Here, we report major recent discoveries presented at the meeting that use genomic approaches to advance our understanding of three major themes: the parallel genetic basis of adaptation; the ecological genomics of speciation; and the evolutionary significance of structural genetic variation. We also suggest future research directions for studies of Mimulus and highlight challenges faced when developing new ecological and evolutionary model systems.

Testing evolutionary hypotheses for DNA barcoding failure in willows.

The goal of DNA barcoding is to enable the rapid identification of taxa from short diagnostic DNA sequence profiles. But how feasible is this objective when many evolutionary processes, such as hybridization and selective sweeps, cause alleles to be shared among related taxa? In this issue of Molecular Ecology, Percy et al. (2014) test the full suite of seven candidate plant barcoding loci in a broad geographic sample of willow species. They show exceptional plastid haplotype sharing between species across continents, with most taxa not possessing a unique barcode sequence. Using population genetic and molecular dating analyses, they implicate hybridization and selective sweeps, but not incomplete lineage sorting, as the historical processes causing widespread haplotype sharing among willow taxa. This study represents an exceptional case of how poorly barcoding can perform, and highlights methodological issues using universal organellar regions for species identification.

The evolution of sex ratio differences and inflorescence architectures in Begonia (Begoniaceae).

PREMISE OF THE STUDY: A major benefit conferred by monoecy is the ability to alter floral sex ratio in response to selection. In monoecious species that produce flowers of a given sex at set positions on the inflorescence, floral sex ratio may be related to inflorescence architecture. We studied the loci underlying differences in inflorescence architecture between two monoecious Begonia species and related this to floral sex ratios. METHODS: We performed trait comparisons and quantitative trait locus (QTL) mapping in a segregating backcross population between Central American Begonia plebeja and B. conchifolia. We focused on traits related to inflorescence architecture, sex ratios, and other reproductive traits. KEY RESULTS: The inflorescence branching pattern of B. conchifolia was more asymmetric than B. plebeja, which in turn affects the floral sex ratio. Colocalizing QTLs of moderate effect influenced both the number of male flowers and the fate decisions of axillary meristems, demonstrating the close link between inflorescence architecture and sex ratio. Additional QTLs were found for stamen number (30% variance explained, VE) and pollen sterility (12.3% VE). CONCLUSIONS: One way in which Begonia species develop different floral sex ratios is through modifications of their inflorescence architecture. The potential pleiotropic action of QTL on inflorescence branching and floral sex ratios may have major implications for trait evolution and responses to selection. The presence of a single QTL of large effect on stamen number may allow rapid divergence for this key floral trait. We propose candidate loci for stamen number and inflorescence branching for future characterization.

The genome sequence of the Annual Mercury, Mercurialis annua L., 1753 (Euphorbiaceae).

We present a genome assembly from a diploid female Mercurialis annua (the Annual Mercury; Tracheophyta; Magnoliopsida; Malpighiales; Euphorbiaceae). The genome sequence is 453.2 megabases in span. Most of the assembly is scaffolded into 8 chromosomal pseudomolecules, including the X chromosome. The organelle genomes have also been assembled, and the mitochondrial genome is 435.28 kilobases in length, while the plastid genome is 169.65 kilobases in length.

The genome sequence of common ivy, Hedera helix L., 1753.

We present a genome assembly from a specimen of Hedera helix (common ivy; Streptophyta; Magnoliopsida; Apiales; Araliaceae). The genome sequence is 1,199.4 megabases in span. Most of the assembly is scaffolded into 24 chromosomal pseudomolecules. The mitochondrial and plastid genomes have also been assembled and are 609.2 and 162.2 kilobases in length respectively.