Incomplete lineage sorting and hybridization underlie tree discordance in Petunia and related genera (Petunieae, Solanaceae)

Despite the overarching history of species divergence, phylogenetic studies often reveal distinct topologies across regions of the genome. The sources of these gene tree discordances are variable, but incomplete lineage sorting (ILS) and hybridization are among those with the most biological importance. Petunia serves as a classic system for studying hybridization in the wild. While field studies suggest that hybridization is frequent, the extent of reticulation within Petunia and its closely related genera has never been examined from a phylogenetic perspective. In this study, we used transcriptomic data from 11 Petunia, 16 Calibrachoa, and 10 Fabiana species to illuminate the relationships between these species and investigate whether hybridization played a significant role in the diversification of the clade. We inferred that gene tree discordance within genera is linked to hybridization events along with high levels of ILS due to their rapid diversification. Moreover, network analyses estimated deeper hybridization events between Petunia and Calibrachoa, genera that have different chromosome numbers. Although these genera cannot hybridize at the present time, ancestral hybridization could have played a role in their parallel radiations, as they share the same habitat and life history.

Genome-wide data reveal bi-direction and asymmetrical hybridization origin of a fern species Microlepia matthewii

Introduction

Natural hybridization is common and plays a crucial role in driving biodiversity in nature. Despite its significance, the understanding of hybridization in ferns remains inadequate. Therefore, it is imperative to study fern hybridization to gain a more comprehensive understanding of fern biodiversity. Our study delves into the role of hybridization in shaping fern species, employing Microlepia matthewii as a case study to investigate its origins of hybridization.

Methods

We performed double digest Genotyping-by-sequencing (dd-GBS) on M. matthewii and its potential parent species, identifying nuclear and chloroplast SNPs. Initially, nuclear SNPs were employed to construct the three cluster analysis: phylogenetic tree, principal component analysis, and population structure analysis. Subsequently, to confirm whether the observed genetic mixture pattern resulted from hybridization, we utilized two methods: ABBA-BABA statistical values in the D-suite program and gene frequency covariance in the Treemix software to detect gene flow. Finally, we employed chloroplast SNPs to construct a phylogenetic tree, tracing the maternal origin.

Results and discussion

The analysis of the nuclear SNP cluster revealed that M. matthewii possesses a genetic composition that is a combination of M. hancei and M. calvescens. Furthermore, the analysis provided strong evidence of significant gene flow signatures from the parental species to the hybrid, as indicated by the two gene flow analyses. The samples of M. matthewii cluster separately with M. hancei or M. calvescens on the chloroplast systematic tree. However, the parentage ratio significantly differs from 1:1, suggesting that M. matthewii is a bidirectional and asymmetrical hybrid offspring of M. hancei and M. calvescens.

Domestication is associated with increased interspecific hybrid compatibility in landfowl (order: Galliformes)

Some species are able to hybridize despite being exceptionally diverged. The causes of this variation in accumulation of reproductive isolation remain poorly understood, and domestication as an impetus or hindrance to reproductive isolation remains to be characterized. In this study, we investigated the role of divergence time, domestication, and mismatches in morphology, habitat, and clutch size among hybridizing species on reproductive isolation in the bird order Galliformes. We compiled and analyzed hybridization occurrences from literature and recorded measures of postzygotic reproductive isolation. We used a text-mining approach leveraging a historical aviculture magazine to quantify the degree of domestication across species. We obtained divergence time, morphology, habitat, and clutch size data from open sources. We found 123 species pairs (involving 77 species) with known offspring fertility (sterile, only males fertile, or both sexes fertile). We found that divergence time and clutch size were significant predictors of reproductive isolation (McFadden's Pseudo-R2 = 0.59), but not habitat or morphological mismatch. Perhaps most interesting, we found a significant relationship between domestication and reproductive compatibility after correcting for phylogeny, removing extreme values, and addressing potential biases (F1,74 = 5.43, R2 = 0.06, P-value = 0.02). We speculate that the genetic architecture and disruption in selective reproductive regimes associated with domestication may impact reproductive isolation, causing domesticated species to be more reproductively labile.

Subgenome evolutionary dynamics in allotetraploid ferns: insights from the gene expression patterns in the allotetraploid species Phegopteris decursivepinnata (Thelypteridacea, Polypodiales)

Allopolyploidization often leads to disruptive conflicts among more than two sets of subgenomes, leading to genomic modifications and changes in gene expression. Although the evolutionary trajectories of subgenomes in allopolyploids have been studied intensely in angiosperms, the dynamics of subgenome evolution remain poorly understood in ferns, despite the prevalence of allopolyploidization. In this study, we have focused on an allotetraploid fern-Phegopteris decursivepinnata-and its diploid parental species, P. koreana (K) and P. taiwaniana (T). Using RNA-seq analyses, we have compared the gene expression profiles for 9,540 genes among parental species, synthetic F1 hybrids, and natural allotetraploids. The changes in gene expression patterns were traced from the F1 hybrids to the natural allopolyploids. This study has revealed that the expression patterns observed in most genes in the F1 hybrids are largely conserved in the allopolyploids; however, there were substantial differences in certain genes between these groups. In the allopolyploids compared with the F1 hybrids, the number of genes showing a transgressive pattern in total expression levels was increased. There was a slight reduction in T-dominance and a slight increase in K-dominance, in terms of expression level dominance. Interestingly, there is no obvious bias toward the T- or K-subgenomes in the number and expression levels overall, showing the absence of subgenome dominance. These findings demonstrated the impacts of the substantial transcriptome change after hybridization and the moderate modification during allopolyploid establishment on gene expression in ferns and provided important insights into subgenome evolution in polyploid ferns.

Frequent allopolyploidy with distant progenitors in the moss genera Physcomitrium and Entosthodon (Funariaceae) identified via subgenome phasing of targeted nuclear genes

Allopolyploids represent a new frontier in species discovery among embryophytes. Within mosses, allopolyploid discovery is challenged by low morphological complexity. The rapid expansion of sequencing approaches in addition to computational developments to identifying genome merger and whole-genome duplication using variation among nuclear loci representing homeologs has allowed for increased allopolyploid discovery among mosses. Here, we test a novel approach to phasing homeologs within loci and phasing loci across subgenomes, or subgenome assignment, called Homologizer, in the family Funariaceae. We confirm the intergeneric hybrid nature of Entosthodon hungaricus, and the allopolyploid origin of Physcomitrium eurystomum and one population of Physcomitrium collenchymatum. We also reveal that hybridization gave rise to Physcomitrium immersum, as well as to yet unrecognized lineages sharing the phenotype of Physcomitrium pyriforme and Physcomitrium sphaericum. Our findings demonstrate the utility of our approach when working with polyploid genomes, and its value in identifying progenitor species using target capture data.

Transcriptome-Based Study on the Phylogeny and Hybridization of Marattialean Ferns (Marattiaceae)

Marattiaceae is a phylogenetically isolated family of tropical eusporangiate ferns including six genera with more than one-hundred species. In Marattiaceae, monophyly of genera has been well-supported phylogenetically. However, the phylogenetic relationships among them were elusive and controversial. Here, a dataset of 26 transcriptomes (including 11 newly generated) were used to assess single-copy nuclear genes and to obtain the organelle gene sequences. Through phylotranscriptomic analysis, the phylogeny and hybridization events of Marattiaceae were explored and a robust phylogenomic framework for the evolution of Marattiaceae was provided. Using both concatenation- and coalescent-based phylogenies, the gene-tree discordance, incomplete lineage sorting (ILS) simulations, and network inference were examined. Except the low support with mitochondrial genes of Marattiaceae, nuclear genes and chloroplast genes strongly supported a sister relationship between Marattiaceae and leptosporangiate ferns. At the genus level, all phylogenetic analysis based on nuclear genes datasets recovered five genera in Marattiaceae as monophyletic with strong support. Danaea and Ptisana were the first two diverged clades in turn. Christensenia was a sister clade to the clade Marattia + Angiopteris s.l. In Angiopteris s.l., three clades (Angiopteris s.s., the Archangiopteris group, and An. sparsisora) were well identified with maximum support. The Archangiopteris group was derived from Angiopteris s.s. at ca. 18 Ma. The putative hybrid species An. sparsisora between Angiopteris s.s. and the Archangiopteris group was verified by the species network analyses and the maternal plastid genes. This study will improve our understanding for using the phylotranscriptomic method to explore phylogeny and investigate hybridization events for difficult taxa in ferns.

Deep reticulation: the long legacy of hybridization in vascular plant evolution

Hybridization has long been recognized as a fundamental evolutionary process in plants but, until recently, our understanding of its phylogenetic distribution and biological significance across deep evolutionary scales has been largely obscure. Over the past decade, genomic and phylogenomic datasets have revealed, perhaps not surprisingly, that hybridization, often associated with polyploidy, has been common throughout the evolutionary history of plants, particularly in various lineages of flowering plants. However, phylogenomic studies have also highlighted the challenges of disentangling signals of ancient hybridization from other sources of genomic conflict (in particular, incomplete lineage sorting). Here, we provide a critical review of ancient hybridization in vascular plants, outlining well-documented cases of ancient hybridization across plant phylogeny, as well as the challenges unique to documenting ancient versus recent hybridization. We provide a definition for ancient hybridization, which, to our knowledge, has not been explicitly attempted before. Further documenting the extent of deep reticulation in plants should remain an important research focus, especially because published examples likely represent the tip of the iceberg in terms of the total extent of ancient hybridization. However, future research should increasingly explore the macroevolutionary significance of this process, in terms of its impact on evolutionary trajectories (e.g. how does hybridization influence trait evolution or the generation of biodiversity over long time scales?), as well as how life history and ecological factors shape, or have shaped, the frequency of hybridization across geologic time and plant phylogeny. Finally, we consider the implications of ubiquitous ancient hybridization for how we conceptualize, analyze, and classify plant phylogeny. Networks, as opposed to bifurcating trees, represent more accurate representations of evolutionary history in many cases, although our ability to infer, visualize, and use networks for comparative analyses is highly limited. Developing improved methods for the generation, visualization, and use of networks represents a critical future direction for plant biology. Current classification systems also do not generally allow for the recognition of reticulate lineages, and our classifications themselves are largely based on evidence from the chloroplast genome. Updating plant classification to better reflect nuclear phylogenies, as well as considering whether and how to recognize hybridization in classification systems, will represent an important challenge for the plant systematics community.

Genome-wide data reveal cryptic diversity and hybridization in a group of tree ferns

Discovery of cryptic diversity is essential to understanding both the process of speciation and the conservation of species. Determining species boundaries in fern lineages represents a major challenge due to lack of morphologically diagnostic characters and frequent hybridization. Genomic data has substantially enhanced our understanding of the speciation process, increased the resolution of species delimitation studies, and led to the discovery of cryptic diversity. Here, we employed restriction-site-associated DNA sequencing (RAD-seq) and integrated phylogenomic and population genomic analyses to investigate phylogenetic relationships and evolutionary history of 16 tree ferns with marginate scales (Cyatheaceae) from China and Vietnam. We conducted multiple species delimitation analyses using the multispecies coalescent (MSC) model and novel approaches based on genealogical divergence index (gdi) and isolation by distance (IBD). In addition, we inferred species trees using concatenation and several coalescent-based methods, and assessed hybridization patterns and rate of gene flow across the phylogeny. We obtained highly supported and generally congruent phylogenies inferred from concatenated and summary-coalescent methods, and the monophyly of all currently recognized species were strongly supported. Our results revealed substantial evidence of cryptic diversity in three widely distributed Gymnosphaera species, each of which was composite of two highly structure lineages that may correspond to cryptic species. We found that hybridization was fairly common between not only closely related species, but also distantly related species. Collectively, it appears that scaly tree ferns may contain cryptic diversity and hybridization has played an important role throughout the evolutionary history of this group.

A need for standardized reporting of introgression: Insights from studies across eukaryotes

With the rise of affordable next-generation sequencing technology, introgression-or the exchange of genetic materials between taxa-has become widely perceived to be a ubiquitous phenomenon in nature. Although this claim is supported by several keystone studies, no thorough assessment of the frequency of introgression across eukaryotes in nature has been performed to date. In this manuscript, we aim to address this knowledge gap by examining patterns of introgression across eukaryotes. We collated a single statistic, Patterson's D, which can be used as a test for introgression across 123 studies to further assess how taxonomic group, divergence time, and sequencing technology influence reports of introgression. Overall, introgression has mostly been measured in plants and vertebrates, with less attention given to the rest of the Eukaryotes. We find that the most frequently used metrics to detect introgression are difficult to compare across studies and even more so across biological systems due to differences in study effort, reporting standards, and methodology. Nonetheless, our analyses reveal several intriguing patterns, including the observation that differences in sequencing technologies may bias values of Patterson's D and that introgression may differ throughout the course of the speciation process. Together, these results suggest the need for a unified approach to quantifying introgression in natural communities and highlight important areas of future research that can be better assessed once this unified approach is met.

Genetic confirmation of a hybrid between two highly divergent cardinalid species: A rose-breasted grosbeak (Pheucticus ludovicianus) and a scarlet tanager (Piranga olivacea)

Using low-coverage whole-genome sequencing, analysis of vocalizations, and inferences from natural history, we document a first-generation hybrid between a rose-breasted grosbeak (Pheucticus ludovicianus) and a scarlet tanager (Piranga olivacea). These two species occur sympatrically throughout much of eastern North America, although were not previously known to interbreed. Following the field identification of a putative hybrid, we use genetic and bioacoustic data to show that a rose-breasted grosbeak was the maternal parent and a scarlet tanager was the paternal parent of the hybrid, whose song was similar to the latter species. These two species diverged >10 million years ago, and thus it is surprising to find a hybrid formed under natural conditions in the wild. Notably, the hybrid has an exceptionally heterozygous genome, with a conservative estimate of a heterozygous base every 100 bp. The observation that this hybrid of such highly divergent parental taxa has survived until adulthood serves as another example of the capacity for hybrid birds to survive with an exceptionally divergent genomic composition.

The Syngameon Enigma

Despite their evolutionary relevance, multispecies networks or syngameons are rarely reported in the literature. Discovering how syngameons form and how they are maintained can give insight into processes such as adaptive radiations, island colonizations, and the creation of new hybrid lineages. Understanding these complex hybridization networks is even more pressing with anthropogenic climate change, as syngameons may have unique synergistic properties that will allow participating species to persist. The formation of a syngameon is not insurmountable, as several ways for a syngameon to form have been proposed, depending mostly on the magnitude and frequency of gene flow events, as well as the relatedness of its participants. Episodic hybridization with small amounts of introgression may keep syngameons stable and protect their participants from any detrimental effects of gene flow. As genomic sequencing becomes cheaper and more species are included in studies, the number of known syngameons is expected to increase. Syngameons must be considered in conservation efforts as the extinction of one participating species may have detrimental effects on the survival of all other species in the network.

Gene-rich germline-restricted chromosomes in black-winged fungus gnats evolved through hybridization

Germline-restricted DNA has evolved in diverse animal taxa and is found in several vertebrate clades, nematodes, and flies. In these lineages, either portions of chromosomes or entire chromosomes are eliminated from somatic cells early in development, restricting portions of the genome to the germline. Little is known about why germline-restricted DNA has evolved, especially in flies, in which 3 diverse families, Chironomidae, Cecidomyiidae, and Sciaridae, carry germline-restricted chromosomes (GRCs). We conducted a genomic analysis of GRCs in the fungus gnat Bradysia (Sciara) coprophila (Diptera: Sciaridae), which has 2 large germline-restricted "L" chromosomes. We sequenced and assembled the genome of B. coprophila and used differences in sequence coverage and k-mer frequency between somatic and germline tissues to identify GRC sequence and compare it to the other chromosomes in the genome. We found that the GRCs in B. coprophila are large, gene rich, and have many genes with divergent homologs on other chromosomes in the genome. We also found that 2 divergent GRCs exist in the population we sequenced. GRC genes are more similar in sequence to genes from another Dipteran family (Cecidomyiidae) than to homologous genes from Sciaridae. This unexpected finding suggests that these chromosomes likely arose in Sciaridae through hybridization with a related lineage. These results provide a foundation from which to answer many questions about the evolution of GRCs in Sciaridae, such as how this hybridization event resulted in GRCs and what features on these chromosomes cause them to be restricted to the germline.

Inferring the Total-Evidence Timescale of Marattialean Fern Evolution in the Face of Model Sensitivity

Phylogenetic divergence-time estimation has been revolutionized by two recent developments: 1) total-evidence dating (or "tip-dating") approaches that allow for the incorporation of fossils as tips in the analysis, with their phylogenetic and temporal relationships to the extant taxa inferred from the data and 2) the fossilized birth-death (FBD) class of tree models that capture the processes that produce the tree (speciation, extinction, and fossilization) and thus provide a coherent and biologically interpretable tree prior. To explore the behavior of these methods, we apply them to marattialean ferns, a group that was dominant in Carboniferous landscapes prior to declining to its modest extant diversity of slightly over 100 species. We show that tree models have a dramatic influence on estimates of both divergence times and topological relationships. This influence is driven by the strong, counter-intuitive informativeness of the uniform tree prior, and the inherent nonidentifiability of divergence-time models. In contrast to the strong influence of the tree models, we find minor effects of differing the morphological transition model or the morphological clock model. We compare the performance of a large pool of candidate models using a combination of posterior-predictive simulation and Bayes factors. Notably, an FBD model with epoch-specific speciation and extinction rates was strongly favored by Bayes factors. Our best-fitting model infers stem and crown divergences for the Marattiales in the mid-Devonian and Late Cretaceous, respectively, with elevated speciation rates in the Mississippian and elevated extinction rates in the Cisuralian leading to a peak diversity of ${\sim}$2800 species at the end of the Carboniferous, representing the heyday of the Psaroniaceae. This peak is followed by the rapid decline and ultimate extinction of the Psaroniaceae, with their descendants, the Marattiaceae, persisting at approximately stable levels of diversity until the present. This general diversification pattern appears to be insensitive to potential biases in the fossil record; despite the preponderance of available fossils being from Pennsylvanian coal balls, incorporating fossilization-rate variation does not improve model fit. In addition, by incorporating temporal data directly within the model and allowing for the inference of the phylogenetic position of the fossils, our study makes the surprising inference that the clade of extant Marattiales is relatively young, younger than any of the fossils historically thought to be congeneric with extant species. This result is a dramatic demonstration of the dangers of node-based approaches to divergence-time estimation, where the assignment of fossils to particular clades is made a priori (earlier node-based studies that constrained the minimum ages of extant genera based on these fossils resulted in much older age estimates than in our study) and of the utility of explicit models of morphological evolution and lineage diversification. [Bayesian model comparison; Carboniferous; divergence-time estimation; fossil record; fossilized birth-death; lineage diversification; Marattiales; models of morphological evolution; Psaronius; RevBayes.].

Deep ancestral introgression shapes evolutionary history of dragonflies and damselflies

Introgression is an important biological process affecting at least 10% of the extant species in the animal kingdom. Introgression significantly impacts inference of phylogenetic species relationships where a strictly binary tree model cannot adequately explain reticulate net-like species relationships. Here we use phylogenomic approaches to understand patterns of introgression along the evolutionary history of a unique, non-model insect system: dragonflies and damselflies (Odonata). We demonstrate that introgression is a pervasive evolutionary force across various taxonomic levels within Odonata. In particular, we show that the morphologically "intermediate" species of Anisozygoptera (one of the three primary suborders within Odonata besides Zygoptera and Anisoptera), which retain phenotypic characteristics of the other two suborders, experienced high levels of introgression likely coming from zygopteran genomes. Additionally, we find evidence for multiple cases of deep inter-superfamilial ancestral introgression.

Cryptic diversity in the model fern genus Ceratopteris (Pteridaceae)

Cryptic species are present throughout the tree of life. They are especially prevalent in ferns, because of processes such hybridization, polyploidy, and reticulate evolution. In addition, the simple morphology of ferns limits phenotypic variation and makes it difficult to detect cryptic species. The model fern genus Ceratopteris has long been suspected to harbor cryptic diversity, in particular within the highly polymorphic C. thalictroides. Yet no studies have included samples from throughout its pan-tropical range or utilized genomic sequencing, making it difficult to assess the full extent of cryptic variation within this genus. Here, we present the first multilocus phylogeny of the genus using reduced representation genomic sequencing (RADseq) and examine population structure, phylogenetic relationships, and ploidy level variation. We recover similar species relationships found in previous studies, find support for the cryptic species C. gaudichaudii as genetically distinct, and identify novel genomic variation within two of the mostly broadly distributed species in the genus, C. thalictroides and C. cornuta. Finally, we detail the utility of our approach for working on cryptic, reticulate groups of ferns. Specifically, it does not require a reference genome, of which there are very few available for ferns. RADseq is a cost-effective way to work with study groups lacking genomic resources, and to obtain the thousands of nuclear markers needed to untangle species complexes.

A consensus phylogenomic approach highlights paleopolyploid and rapid radiation in the history of Ericales

PREMISE

Large genomic data sets offer the promise of resolving historically recalcitrant species relationships. However, different methodologies can yield conflicting results, especially when clades have experienced ancient, rapid diversification. Here, we analyzed the ancient radiation of Ericales and explored sources of uncertainty related to species tree inference, conflicting gene tree signal, and the inferred placement of gene and genome duplications.

METHODS

We used a hierarchical clustering approach, with tree-based homology and orthology detection, to generate six filtered phylogenomic matrices consisting of data from 97 transcriptomes and genomes. Support for species relationships was inferred from multiple lines of evidence including shared gene duplications, gene tree conflict, gene-wise edge-based analyses, concatenation, and coalescent-based methods, and is summarized in a consensus framework.

RESULTS

Our consensus approach supported a topology largely concordant with previous studies, but suggests that the data are not capable of resolving several ancient relationships because of lack of informative characters, sensitivity to methodology, and extensive gene tree conflict correlated with paleopolyploidy. We found evidence of a whole-genome duplication before the radiation of all or most ericalean families, and demonstrate that tree topology and heterogeneous evolutionary rates affect the inferred placement of genome duplications.

CONCLUSIONS

We provide several hypotheses regarding the history of Ericales, and confidently resolve most nodes, but demonstrate that a series of ancient divergences are unresolvable with these data. Whether paleopolyploidy is a major source of the observed phylogenetic conflict warrants further investigation.

Geography is essential for reproductive isolation between florally diversified morning glory species from Amazon canga savannahs

The variety, relative importance and eco-evolutionary stability of reproductive barriers are critical to understanding the processes of speciation and species persistence. Here we evaluated the strength of the biotic prezygotic and postzygotic isolation barriers between closely related morning glory species from Amazon canga savannahs. The flower geometry and flower visitor assemblage analyses supported pollination by the bees in lavender-flowered Ipomoea marabaensis and recruitment of hummingbirds as pollinators in red-flowered Ipomoea cavalcantei. Nevertheless, native bee species and alien honeybees foraged on flowers of both species. Real-time interspecific hybridization underscored functionality of the overlap in flower visitor assemblages, questioning the strength of prezygotic isolation underpinned by diversification in flower colour and geometry. Interspecific hybrids were fertile and produced offspring in nature. No significant asymmetry in interspecific hybridization and hybrid incompatibilities among offspring were found, indicating weak postmating and postzygotic isolation. The results suggested that despite floral diversification, the insular-type geographic isolation remains a major barrier to gene flow. Findings set a framework for the future analysis of contemporary evolution of plant-pollinator networks at the population, community, and ecosystem levels in tropical ecosystems that are known to be distinct from the more familiar temperate climate models.

Correlates of hybridization in plants

Hybridization is a biological phenomenon increasingly recognized as an important evolutionary process in both plants and animals, as it is linked to speciation, radiation, extinction, range expansion and invasion, and allows for increased trait diversity in agricultural and horticultural systems. Estimates of hybridization frequency vary across taxonomic groups, but causes of this variation are unknown. Here, we ask on a global scale whether hybridization is linked to any of 11 traits related to plant life history, reproduction, genetic predisposition, and environment or opportunity. Given that hybridization is not evenly distributed across the plant tree of life, we use phylogenetic generalized least squares regression models and phylogenetic path analysis to detect statistical associations between hybridization and plant traits at both the family and genus levels. We find that perenniality and woodiness are each weakly associated with an increased frequency of hybridization in univariate analyses, but path analysis suggests that the direct linkage is between perenniality and increased hybridization (with woodiness having only an indirect relationship with hybridization via perenniality). Weak associations between higher rates of hybridization and higher outcrossing rates, abiotic pollination syndromes, vegetative reproductive modes, larger genomes, and less variable genome sizes are detectable in some cases but not others. We argue that correlational evidence at the global scale, such as that presented here, provides a robust framework for forming hypotheses to examine and test drivers of hybridization at a more mechanistic level.

Asymmetric hybridization in Central European populations of the Dryopteris carthusiana group

PREMISE

Hybridization is a key process in plant speciation. Despite its importance, there is no detailed study of hybridization rates in fern populations. A proper estimate of hybridization rates is needed to understand factors regulating hybridization.

METHODS

We studied hybridization in the European Dryopteris carthusiana group, represented by one diploid and two tetraploid species and their hybrids. We sampled ~100 individuals per population in 40 mixed populations of the D. carthusiana group across Europe. All plants were identified by measuring genome size (DAPI staining) using flow cytometry. To determine the maternal parentage of hybrids, we sequenced the chloroplast region trnL-trnF of all taxa involved.

RESULTS

We found hybrids in 85% of populations. Triploid D. ×ambroseae occurred in every population that included both parent species and is most abundant when the parent species are equally abundant. By contrast, tetraploid D. ×deweveri was rare (15 individuals total) and triploid D. ×sarvelae was absent. The parentage of hybrid taxa is asymmetric. Despite expectations from previous studies, tetraploid D. dilatata is the predominant male parent of its triploid hybrid.

CONCLUSIONS

This is a thorough investigation of hybridization rates in natural populations of ferns. Hybridization rates differ greatly even among closely related fern taxa. In contrast to angiosperms, our data suggest that hybridization rates are highest in balanced parent populations and support the notion that some ferns possess very weak barriers to hybridization. Our results from sequencing cpDNA challenge established notions about the correlation of ploidy level and mating tendencies.

Genetic and morphological identification of a recurrent Dicksonia tree fern hybrid in New Zealand

Hybridization is common in many ferns and has been a significant factor in fern evolution and speciation. However, hybrids are rare between the approximately 30 species of Dicksonia tree ferns world-wide, and none are well documented. In this study we examine the relationship of a newly-discovered Dicksonia tree fern from Whirinaki, New Zealand, which does not fit the current taxonomy of the three species currently recognized in New Zealand. Our microsatellite genotyping and ddRAD-seq data indicate these plants are F1 hybrids that have formed multiple times between D. fibrosa and D. lanata subsp. lanata. The Whirinaki plants have intermediate morphology between D. fibrosa and D. lanata subsp. lanata and their malformed spores are consistent with a hybrid origin. The Whirinaki plants–Dicksonia fibrosa × D. lanata subsp. lanata–are an example of hybridization between distantly related fern lineages, with the two parent species estimated to have diverged 55–25 mya. Our chloroplast sequencing indicates asymmetric chloroplast inheritance in the Whirinaki morphology with D. lanata subsp. lanata always contributing the chloroplast genome.

The rise of the Andes promoted rapid diversification in Neotropical Phlegmariurus (Lycopodiaceae)

Tropical mountains are disproportionately biodiverse relative to their surface area, but the processes underlying their exceptional diversity require further study. Here, we use comparative phylogenetic methods to examine the impact of the Andean orogeny on the diversification of Neotropical Phlegmariurus, a species-rich lycophyte clade. We generated a time-calibrated phylogeny of 105 species of Neotropical Phlegmariurus and estimated lineage diversification rates. We tested for correlations between lineage diversification rates and species range size, niche breadth, elevational range amplitude, and mean elevation of occurrence. A recently developed macroevolutionary model was used to incorporate geological data and test for an association between diversification rates and the Andean uplift. Diversification rates of Neotropical Phlegmariurus are negatively correlated with species range size and positively correlated with mean elevation of species occurrence. The rise of the Andes is strongly associated with increased rates of diversification in Neotropical Phlegmariurus during the last 10 Myr. Our study demonstrates the importance of mountain-building events and geographical isolation of alpine populations as drivers of rapid diversification, even in spore-dispersed plants. This work also highlights the usefulness of combined phylogenetic, geological and ecological datasets, and the promise of comparative environment-dependent diversification models in better understanding the evolutionary origins of biodiversity.

Whole-Genome Duplication and Plant Macroevolution

Whole-genome duplication (WGD) is characteristic of almost all fundamental lineages of land plants. Unfortunately, the timings of WGD events are loosely constrained and hypotheses of evolutionary consequence are poorly formulated, making them difficult to test. Using examples from across the plant kingdom, we show that estimates of timing can be improved through the application of molecular clock methodology to multigene datasets. Further, we show that phenotypic change can be quantified in morphospaces and that relative phenotypic disparity can be compared in the light of WGD. Together, these approaches facilitate tests of hypotheses on the role of WGD in plant evolution, underscoring the potential of plants as a model system for investigating the role WGD in macroevolution.

Dissection by genomic and plumage variation of a geographically complex hybrid zone between two Australian non-sister parrot species, Platycercus adscitus and Platycercus eximius

The study of hybrid zones advances understanding of the speciation process, and approaches incorporating genomic data are increasingly used to draw significant conclusions about the impact of hybridisation. Despite the progress made, the complex interplay of factors that can lead to substantially variable hybridisation outcomes are still not well understood, and many systems and/or groups remain comparatively poorly studied. Our study aims to broaden the literature on avian hybrid zones, investigating a potentially geographically and temporally complex putative hybrid zone between two native Australian non-sister parrot species, the pale-headed and eastern rosellas (Platycercus adscitus and Platycercus eximius, respectively). We analysed six plumage traits and >1400 RADseq loci and detected hybrid individuals and an unexpectedly complex geographic structure. The hybrid zone is larger than previously described due to either observer bias or its movement over recent decades. It comprises different subregions where genetic and plumage signals of admixture vary markedly in their concordance. Evidence of contemporary hybridisation (later generation and backcrossed individuals) both within and beyond the previously defined zone, when coupled with a lack of F1 hybrids and differential patterns of introgression among potentially diagnostic loci, indicates a lack of post-zygotic barriers to gene flow between species. Despite ongoing gene flow, species boundaries are likely maintained largely by strong pre-mating barriers. These findings are discussed in detail and future avenues for research into this system are proposed, which would be of benefit to the speciation and hybrid zone literature.

Target sequence capture of nuclear-encoded genes for phylogenetic analysis in ferns

PREMISE OF THE STUDY

Until recently, most phylogenetic studies of ferns were based on chloroplast genes. Evolutionary inferences based on these data can be incomplete because the characters are from a single linkage group and are uniparentally inherited. These limitations are particularly acute in studies of hybridization, which is prevalent in ferns; fern hybrids are common and ferns are able to hybridize across highly diverged lineages, up to 60 million years since divergence in one documented case. However, it not yet clear what effect such hybridization has on fern evolution, in part due to a paucity of available biparentally inherited (nuclear-encoded) markers.

METHODS

We designed oligonucleotide baits to capture 25 targeted, low-copy nuclear markers from a sample of 24 species spanning extant fern diversity.

RESULTS

Most loci were successfully sequenced from most accessions. Although the baits were designed from exon (transcript) data, we successfully captured intron sequences that should be useful for more focused phylogenetic studies. We present phylogenetic analyses of the new target sequence capture data and integrate these into a previous transcript-based data set.

DISCUSSION

We make our bait sequences available to the community as a resource for further studies of fern phylogeny.

Supervised machine learning reveals introgressed loci in the genomes of Drosophila simulans and D. sechellia

Hybridization and gene flow between species appears to be common. Even though it is clear that hybridization is widespread across all surveyed taxonomic groups, the magnitude and consequences of introgression are still largely unknown. Thus it is crucial to develop the statistical machinery required to uncover which genomic regions have recently acquired haplotypes via introgression from a sister population. We developed a novel machine learning framework, called FILET (Finding Introgressed Loci via Extra-Trees) capable of revealing genomic introgression with far greater power than competing methods. FILET works by combining information from a number of population genetic summary statistics, including several new statistics that we introduce, that capture patterns of variation across two populations. We show that FILET is able to identify loci that have experienced gene flow between related species with high accuracy, and in most situations can correctly infer which population was the donor and which was the recipient. Here we describe a data set of outbred diploid Drosophila sechellia genomes, and combine them with data from D. simulans to examine recent introgression between these species using FILET. Although we find that these populations may have split more recently than previously appreciated, FILET confirms that there has indeed been appreciable recent introgression (some of which might have been adaptive) between these species, and reveals that this gene flow is primarily in the direction of D. simulans to D. sechellia.

Understanding the extent to which species or diverged populations hybridize in nature is crucially important if we are to understand the speciation process. Accordingly numerous research groups have developed methodology for finding the genetic evidence of such introgression. In this report we develop a supervised machine learning approach for uncovering loci which have introgressed across species boundaries. We show that our method, FILET, has greater accuracy and power than competing methods in discovering introgression, and in addition can detect the directionality associated with the gene flow between species. Using whole genome sequences from Drosophila simulans and Drosophila sechellia we show that FILET discovers quite extensive introgression between these species that has occurred mostly from D. simulans to D. sechellia. Our work highlights the complex process of speciation even within a well-studied system and points to the growing importance of supervised machine learning in population genetics.

Next-generation polyploid phylogenetics: rapid resolution of hybrid polyploid complexes using PacBio single-molecule sequencing

Difficulties in generating nuclear data for polyploids have impeded phylogenetic study of these groups. We describe a high-throughput protocol and an associated bioinformatics pipeline (Pipeline for Untangling Reticulate Complexes (Purc)) that is able to generate these data quickly and conveniently, and demonstrate its efficacy on accessions from the fern family Cystopteridaceae. We conclude with a demonstration of the downstream utility of these data by inferring a multi-labeled species tree for a subset of our accessions. We amplified four c. 1-kb-long nuclear loci and sequenced them in a parallel-tagged amplicon sequencing approach using the PacBio platform. Purc infers the final sequences from the raw reads via an iterative approach that corrects PCR and sequencing errors and removes PCR-mediated recombinant sequences (chimeras). We generated data for all gene copies (homeologs, paralogs, and segregating alleles) present in each of three sets of 50 mostly polyploid accessions, for four loci, in three PacBio runs (one run per set). From the raw sequencing reads, Purc was able to accurately infer the underlying sequences. This approach makes it easy and economical to study the phylogenetics of polyploids, and, in conjunction with recent analytical advances, facilitates investigation of broad patterns of polyploid evolution.

Species-rich and polyploid-poor: Insights into the evolutionary role of whole-genome duplication from the Cape flora biodiversity hotspot

PREMISE OF THE STUDY

Whole-genome duplication (WGD) in angiosperms has been hypothesized to be advantageous in unstable environments and/or to increase diversification rates, leading to radiations. Under the first hypothesis, floras in stable environments are predicted to have lower proportions of polyploids than highly, recently disturbed floras, whereas species-rich floras would be expected to have higher than expected proportions of polyploids under the second. The South African Cape flora is used to discriminate between these two hypotheses because it features a hyperdiverse flora predominantly generated by a limited number of radiations (Cape clades), against a backdrop of climatic and geological stability.

METHODS

We compiled all known chromosome counts for species in 21 clades present in the Cape (1653 species, including 24 Cape clades), inferred ploidy levels for these species by inspection or derived from the primary literature, and compared Cape to non-Cape ploidy levels in these clades (17,520 species) using G tests.

KEY RESULTS

The Cape flora has anomalously low proportions of polyploids compared with global levels. This pattern is consistently observed across nearly half the clades and across global latitudinal gradients, although individual lineages seem to be following different paths to low levels of WGD and to differing degrees.

CONCLUSIONS

This pattern shows that the diversity of the Cape flora is the outcome of primarily diploid radiations and supports the hypothesis that WGD may be rare in stable environments.

The challenges to inferring the regulators of biodiversity in deep time

Attempts to infer the ecological drivers of macroevolution in deep time have long drawn inspiration from work on extant systems, but long-term evolutionary and geological changes complicate the simple extrapolation of such theory. Recent efforts to incorporate a more informed ecology into macroevolution have moved beyond the descriptive, seeking to isolate generating mechanisms and produce testable hypotheses of how groups of organisms usurp each other or coexist over vast timespans. This theme issue aims to exemplify this progress, providing a series of case studies of how novel modelling approaches are helping infer the regulators of biodiversity in deep time. In this Introduction, we explore the challenges of these new approaches. First, we discuss how our choices of taxonomic units have implications for the conclusions drawn. Second, we emphasize the need to embrace the interdependence of biotic and abiotic changes, because no living organism ignores its environment. Third, in the light of parts 1 and 2, we discuss the set of dynamic signatures that we might expect to observe in the fossil record. Finally, we ask whether these dynamics represent the most ecologically informative foci for research efforts aimed at inferring the regulators of biodiversity in deep time. The papers in this theme issue contribute in each of these areas.

Origins and diversity of a cosmopolitan fern genus on an island archipelago

Isolated oceanic islands are characterized by patterns of biological diversity different from that on nearby continental mainlands. Isolation can provide the opportunity for evolutionary divergence, but also set the stage for hybridization between related taxa arriving from different sources. Ferns disperse by haploid spores, which are produced in large numbers and can travel long distances in air currents, enabling these plants to become established on most oceanic islands. Here, we examine the origins and patterns of diversity of the cosmopolitan fern genus Pteridium (Dennstaedtiaceae; bracken) on the Galapagos Islands. We use nucleotide sequences from two plastid genes, and two nuclear gene markers, to examine phylogeography of Pteridium on the Galapagos Islands. We incorporate data from a previous study to provide a worldwide context. We also sampled new specimens from South and Central America. We used flow cytometry to estimate genome size of some accessions. We found that both plastid and nuclear haplotypes fall into two distinct clades, consistent with a two-diploid-species taxonomy of P. aquilinum and P. esculentum. As predicted, the allotetraploid P. caudatum possesses nuclear haplotypes from both diploid species. Samples from the Galapagos include P. esculentum subsp. arachnoideum, P. caudatum and possible hybrids between them. Multiple Pteridium taxa were also observed growing together at some sites. We find evidence for multiple origins of Pteridium on the Galapagos Islands and multiple origins of tetraploid P. caudatum throughout its range in Central and South America. We also posit that P. caudatum may include recent diploid hybrids, backcrosses to P. esculentum, as well as allotetraploid plants. The Galapagos Islands are positioned close to the equator where they can receive dispersing propagules from both hemispheres. This may partly explain the high levels of diversity found for this cosmopolitan fern on these islands.

Why are there so few species of ferns?

A recent study has documented a natural hybridization event between two fern lineages that last shared a common ancestor about 60 million years ago. This is one of the deepest hybridization events ever described and has important implications for plant speciation theory.

The evolutionary history of ferns inferred from 25 low-copy nuclear genes

•

PREMISE OF THE STUDY

Understanding fern (monilophyte) phylogeny and its evolutionary timescale is critical for broad investigations of the evolution of land plants, and for providing the point of comparison necessary for studying the evolution of the fern sister group, seed plants. Molecular phylogenetic investigations have revolutionized our understanding of fern phylogeny, however, to date, these studies have relied almost exclusively on plastid data.•

METHODS

Here we take a curated phylogenomics approach to infer the first broad fern phylogeny from multiple nuclear loci, by combining broad taxon sampling (73 ferns and 12 outgroup species) with focused character sampling (25 loci comprising 35877 bp), along with rigorous alignment, orthology inference and model selection.•

KEY RESULTS

Our phylogeny corroborates some earlier inferences and provides novel insights; in particular, we find strong support for Equisetales as sister to the rest of ferns, Marattiales as sister to leptosporangiate ferns, and Dennstaedtiaceae as sister to the eupolypods. Our divergence-time analyses reveal that divergences among the extant fern orders all occurred prior to ∼200 MYA. Finally, our species-tree inferences are congruent with analyses of concatenated data, but generally with lower support. Those cases where species-tree support values are higher than expected involve relationships that have been supported by smaller plastid datasets, suggesting that deep coalescence may be reducing support from the concatenated nuclear data.•

CONCLUSIONS

Our study demonstrates the utility of a curated phylogenomics approach to inferring fern phylogeny, and highlights the need to consider underlying data characteristics, along with data quantity, in phylogenetic studies.