Reticulate evolution on a global scale: a nuclear phylogeny for New World Dryopteris (Dryopteridaceae)

Niche and phenotypic differentiation in fern hybrid speciation, a case study of Pteris fauriei (Pteridaceae)

BACKGROUND AND AIMS

Niche differentiation is a crucial issue in speciation. Although it has a well-known role in adaptive processes of hybrid angiosperms, it is less understood in hybrid ferns. Here, we investigate whether an intermediate ecological niche of a fern hybrid is a novel adaptation that provides insights into fern hybrid speciation.

METHODS

Pteris fauriei (Pteridaceae) is a natural hybrid fern, occurring in environments between its parent species. The maternal Pteris minor is found in sunny areas, but the habitat of the paternal Pteris latipinna is shady. We combined data from morphology, leaf anatomy and photosynthetic traits to explore adaptation and differentiation, along with measuring the environmental features of their niches. We also performed experiments in a common garden to understand ecological plasticity.

KEY RESULTS

The hybrid P. fauriei was intermediate between the parent species in stomatal density, leaf anatomical features and photosynthetic characteristics in both natural habitats and a common garden. Interestingly, the maternal P. minor showed significant environmental plasticity and was more similar to the hybrid P. fauriei in the common garden, suggesting that the maternal species experiences stress in its natural habitats but thrives in environments similar to those of the hybrid.

CONCLUSIONS

Based on the similar niche preferences of the hybrid and parents, we propose hybrid superiority. Our results indicate that the hybrid P. fauriei exhibits greater fitness and can compete with and occupy the initial niches of the maternal P. minor. Consequently, we suggest that the maternal P. minor has experienced a niche shift, elucidating the pattern of niche differentiation in this hybrid group. These findings offer a potential explanation for the frequent occurrence of hybridization in ferns and provide new insights into fern hybrid speciation, enhancing our understanding of fern diversity.

Compositae-ParaLoss-1272: A complementary sunflower-specific probe set reduces paralogs in phylogenomic analyses of complex systems

Premise

A family-specific probe set for sunflowers, Compositae-1061, enables family-wide phylogenomic studies and investigations at lower taxonomic levels, but may lack resolution at genus to species levels, especially in groups complicated by polyploidy and hybridization.

Methods

We developed a Hyb-Seq probe set, Compositae-ParaLoss-1272, that targets orthologous loci in Asteraceae. We tested its efficiency across the family by simulating target enrichment sequencing in silico. Additionally, we tested its effectiveness at lower taxonomic levels in the historically complex genus Packera. We performed Hyb-Seq with Compositae-ParaLoss-1272 for 19 Packera taxa that were previously studied using Compositae-1061. The resulting sequences from each probe set, plus a combination of both, were used to generate phylogenies, compare topologies, and assess node support.

Results

We report that Compositae-ParaLoss-1272 captured loci across all tested Asteraceae members, had less gene tree discordance, and retained longer loci than Compositae-1061. Most notably, Compositae-ParaLoss-1272 recovered substantially fewer paralogous sequences than Compositae-1061, with only ~5% of the recovered loci reporting as paralogous, compared to ~59% with Compositae-1061.

Discussion

Given the complexity of plant evolutionary histories, assigning orthology for phylogenomic analyses will continue to be challenging. However, we anticipate Compositae-ParaLoss-1272 will provide improved resolution and utility for studies of complex groups and lower taxonomic levels in the sunflower family.

Interploidy gene flow via a 'pentaploid bridge' and ploidy reduction in Cystopteris fragilis fern complex (Cystopteridaceae: Polypodiales)

KEY MESSAGE

Our results indicate the existence of interploidy gene flow in Cystopteris fragilis, resulting in sexual triploid and diploid gametophytes from pentaploid parents. Similar evolutionary dynamics might operate in other fern complexes and need further investigation. Polyploidization and hybridization are a key evolutionary processes in ferns. Here, we outline an interploidy gene flow pathway operating in the polyploid Cystopteris fragilis complex. The conditions necessary for the existence of this pathway were tested. A total of 365 C. fragilis individuals were collected covering representatives of all three predominant ploidy levels (tetraploid, pentaploid, and hexaploid), cultivated, had their ploidy level estimated by flow cytometry, and their spores collected. The spores, as well as gametophytes and sporophytes established from them, were analysed by flow cytometry. Spore abortion rate was also estimated. In tetraploids, we observed the formation of unreduced (tetraploid) spores (ca 2%). Collected pentaploid individuals indicate ongoing hybridization between ploidy levels. Pentaploids formed up to 52% viable spores, ca 79% of them reduced, i.e. diploid and triploid. Reduced spores formed viable gametophytes, and, in the case of triploids, filial hexaploid sporophytes, showing evidence of sexual reproduction. Some tetraploid sporophytes reproduce apomictically (based on uniform ploidy of their metagenesis up to filial sporophytes). Triploid and diploid gametophytes from pentaploid parents are able to mate among themselves, or with "normal" reduced gametophytes from the sexual tetraploid sporophytes (the dominant ploidy level in the sporophytes in this populations), to produce tetraploid, pentaploid, and hexaploid sporophytes, allowing for geneflow from the pentaploids to both the tetraploid and hexaploid populations. Similar evolutionary dynamics might operate in other fern complexes and need further investigation.

Where to draw the boundaries? Using landscape genomics to disentangle the scribbly gum species complex

PREMISE

Species delimitation is an integral part of evolution and ecology and is vital in conservation science. However, in some groups, species delimitation is difficult, especially where ancestral relationships inferred from morphological or genetic characters are discordant, possibly due to a complicated demographic history (e.g., recent divergences between lineages). Modern genetic techniques can take into account complex histories to distinguish species at a reasonable cost and are increasingly used in numerous applications. We focus on the scribbly gums, a group of up to five closely related and morphologically similar "species" within the eucalypts.

METHODS

Multiple populations of each recognized scribbly gum species were sampled over a wide region across climates, and genomewide scans were used to resolve species boundaries.

RESULTS

None of the taxa were completely divergent, and there were two genetically distinct entities: the inland distributed Eucalyptus rossii and a coastal conglomerate consisting of four species forming three discernible, but highly admixed groups. Divergence among taxa was likely driven by temporal vicariant processes resulting in partial separation across biogeographic barriers. High interspecific gene flow indicated separated taxa reconnected at different points in time, blurring species boundaries.

CONCLUSIONS

Our results highlight the need for genetic screening when dealing with closely related taxonomic entities, particularly those with modest morphological differences. We show that high-throughput sequencing can be effective at identifying species groupings and processes driving divergence, even in the most taxonomically complex groups, and be used as a standard practice for disentangling species complexes.

Dynamic hybridization between two spleenworts, Asplenium incisum and Asplenium ruprechtii in Korea

Natural hybridization between Asplenium incisum and A. ruprechtii has been observed in Northeast Asia and its allotetraploid species, A. castaneoviride, was reported. However, the hybridization process between the parental species and the origin of the allotetraploid taxon remains obscure. Additionally, the systematic affinities of the recently described hybrid A. bimixtum, considered to have originated from the hybridization of A. ruprechtii, A. trichomanes, and A. incisum, is unresolved owing to its similarity to A. castaneoviride. The goals of this study were to (1) investigate the hybridization between A. ruprechtii and A. incisum; (2) verify the origin of A. castaneoviride occurring in Korea, whether it independently arose from 2x sterile hybrids; and (3) elucidate the reliability of identifying A. bimixtum. Three genotypes, A. incisum, A. ruprechtii, and their hybrid, were identified based on the nuclear gene pgiC sequence and finally divided them into six types by ploidy levels: diploid A. incisum, A. ruprechtii, and four hybrid types (diploid A. × castaneoviride, triploid A. × castaneoviride, allotetraploid A. castaneoviride, and A. bimixtum). In the analyses of plastid DNA, all hybrids had an A. ruprechtii-type rbcL gene. In addition, the four plastomes of A. ruprechtii and the hybrids had high pairwise sequence identities greater than 98.48%. They increased up to 99.88% when a large deletion of A. x castaneoriviride (2x) collected from Buramsan populations was ignored. Notably, this large deletion was also found in triploid A. × castaneoviride and allotetraploid A. castaneoviride in the same populations. Sequence data of the nuclear and plastid genes showed that hybridization is unidirectional, and A. ruprechtii is the maternal parent. The large deletion of rpoC2-rps2 commonly found in the different ploidy hybrids of the Buramsan population suggests that the allotetraploid A. castaneoviride can be created independently from sterile hybrids. We assume that both polyploidization driving allopolyploidy and minority cytotype exclusion took place independently in the population, since A castaenoviride co-occurs with A. ruprechtii in small populations. Furthermore, it was also observed that an enlarged noncoding region in fern organelle (ENRIFO) of the plastome was found in the genus Asplenium.

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.

Phylogenomic resolution of order- and family-level monocot relationships using 602 single-copy nuclear genes and 1375 BUSCO genes

We assess relationships among 192 species in all 12 monocot orders and 72 of 77 families, using 602 conserved single-copy (CSC) genes and 1375 benchmarking single-copy ortholog (BUSCO) genes extracted from genomic and transcriptomic datasets. Phylogenomic inferences based on these data, using both coalescent-based and supermatrix analyses, are largely congruent with the most comprehensive plastome-based analysis, and nuclear-gene phylogenomic analyses with less comprehensive taxon sampling. The strongest discordance between the plastome and nuclear gene analyses is the monophyly of a clade comprising Asparagales and Liliales in our nuclear gene analyses, versus the placement of Asparagales and Liliales as successive sister clades to the commelinids in the plastome tree. Within orders, around six of 72 families shifted positions relative to the recent plastome analysis, but four of these involve poorly supported inferred relationships in the plastome-based tree. In Poales, the nuclear data place a clade comprising Ecdeiocoleaceae+Joinvilleaceae as sister to the grasses (Poaceae); Typhaceae, (rather than Bromeliaceae) are resolved as sister to all other Poales. In Commelinales, nuclear data place Philydraceae sister to all other families rather than to a clade comprising Haemodoraceae+Pontederiaceae as seen in the plastome tree. In Liliales, nuclear data place Liliaceae sister to Smilacaceae, and Melanthiaceae are placed sister to all other Liliales except Campynemataceae. Finally, in Alismatales, nuclear data strongly place Tofieldiaceae, rather than Araceae, as sister to all the other families, providing an alternative resolution of what has been the most problematic node to resolve using plastid data, outside of those involving achlorophyllous mycoheterotrophs. As seen in numerous prior studies, the placement of orders Acorales and Alismatales as successive sister lineages to all other extant monocots. Only 21.2% of BUSCO genes were demonstrably single-copy, yet phylogenomic inferences based on BUSCO and CSC genes did not differ, and overall functional annotations of the two sets were very similar. Our analyses also reveal significant gene tree-species tree discordance despite high support values, as expected given incomplete lineage sorting (ILS) related to rapid diversification. Our study advances understanding of monocot relationships and the robustness of phylogenetic inferences based on large numbers of nuclear single-copy genes that can be obtained from transcriptomes and genomes.

Deep vicariance and frequent transoceanic dispersal shape the evolutionary history of a globally distributed fern family

PREMISE

The historical biogeography of ferns is typically expected to be dominated by long-distance dispersal due to their minuscule spores. However, few studies have inferred the historical biogeography of a large and widely distributed group of ferns to test this hypothesis. Our aims were to determine the extent to which long-distance dispersal vs. vicariance have shaped the history of the fern family Blechnaceae, to explore ecological correlates of dispersal and diversification, and to determine whether these patterns differ between the northern and southern hemispheres.

METHODS

We used sequence data for three chloroplast loci to infer a time-calibrated phylogeny for 154 of 265 species of Blechnaceae, including representatives of all genera in the family. This tree was used to conduct ancestral range reconstruction and stochastic character mapping, estimate diversification rates, and identify ecological correlates of diversification.

RESULTS

Blechnaceae originated in Eurasia and began diversifying in the late Cretaceous. A lineage comprising most extant diversity diversified principally in the austral Pacific region around the Paleocene-Eocene Thermal Maximum. Land connections that existed near the poles during periods of warm climates likely facilitated migration of several lineages, with subsequent climate-mediated vicariance shaping current distributions. Long-distance dispersal is frequent and asymmetrical, with New Zealand/Pacific Islands, Australia, and tropical America being major source areas.

CONCLUSIONS

Ancient vicariance and extensive long-distance dispersal have shaped the history of Blechnaceae in both the northern and southern hemispheres. The exceptional diversity in austral regions appears to reflect rapid speciation in these areas; mechanisms underlying this evolutionary success remain uncertain.

An open and continuously updated fern tree of life

Ferns, with about 12,000 species, are the second most diverse lineage of vascular plants after angiosperms. They have been the subject of numerous molecular phylogenetic studies, resulting in the publication of trees for every major clade and DNA sequences from nearly half of all species. Global fern phylogenies have been published periodically, but as molecular systematics research continues at a rapid pace, these become quickly outdated. Here, we develop a mostly automated, reproducible, open pipeline to generate a continuously updated fern tree of life (FTOL) from DNA sequence data available in GenBank. Our tailored sampling strategy combines whole plastomes (few taxa, many loci) with commonly sequenced plastid regions (many taxa, few loci) to obtain a global, species-level fern phylogeny with high resolution along the backbone and maximal sampling across the tips. We use a curated reference taxonomy to resolve synonyms in general compliance with the community-driven Pteridophyte Phylogeny Group I classification. The current FTOL includes 5,582 species, an increase of ca. 40% relative to the most recently published global fern phylogeny. Using an updated and expanded list of 51 fern fossil constraints, we find estimated ages for most families and deeper clades to be considerably older than earlier studies. FTOL and its accompanying datasets, including the fossil list and taxonomic database, will be updated on a regular basis and are available via a web portal (https://fernphy.github.io) and R packages, enabling immediate access to the most up-to-date, comprehensively sampled fern phylogeny. FTOL will be useful for anyone studying this important group of plants over a wide range of taxonomic scales, from smaller clades to the entire tree. We anticipate FTOL will be particularly relevant for macroecological studies at regional to global scales and will inform future taxonomic systems with the most recent hypothesis of fern phylogeny.

Phylotranscriptomics Illuminates the Placement of Whole Genome Duplications and Gene Retention in Ferns

Ferns are the second largest clade of vascular plants with over 10,000 species, yet the generation of genomic resources for the group has lagged behind other major clades of plants. Transcriptomic data have proven to be a powerful tool to assess phylogenetic relationships, using thousands of markers that are largely conserved across the genome, and without the need to sequence entire genomes. We assembled the largest nuclear phylogenetic dataset for ferns to date, including 2884 single-copy nuclear loci from 247 transcriptomes (242 ferns, five outgroups), and investigated phylogenetic relationships across the fern tree, the placement of whole genome duplications (WGDs), and gene retention patterns following WGDs. We generated a well-supported phylogeny of ferns and identified several regions of the fern phylogeny that demonstrate high levels of gene tree-species tree conflict, which largely correspond to areas of the phylogeny that have been difficult to resolve. Using a combination of approaches, we identified 27 WGDs across the phylogeny, including 18 large-scale events (involving more than one sampled taxon) and nine small-scale events (involving only one sampled taxon). Most inferred WGDs occur within single lineages (e.g., orders, families) rather than on the backbone of the phylogeny, although two inferred events are shared by leptosporangiate ferns (excluding Osmundales) and Polypodiales (excluding Lindsaeineae and Saccolomatineae), clades which correspond to the majority of fern diversity. We further examined how retained duplicates following WGDs compared across independent events and found that functions of retained genes were largely convergent, with processes involved in binding, responses to stimuli, and certain organelles over-represented in paralogs while processes involved in transport, organelles derived from endosymbiotic events, and signaling were under-represented. To date, our study is the most comprehensive investigation of the nuclear fern phylogeny, though several avenues for future research remain unexplored.

Reticulate evolution in the Pteris fauriei group (Pteridaceae)

The Pteris fauriei group (Pteridaceae) has a wide distribution in Eastern Asia and includes 18 species with similar but varied morphology. We collected more than 300 specimens of the P. fauriei group and determined ploidy by flow cytometry and inferred phylogenies by molecular analyses of chloroplast and nuclear DNA markers. Our results reveal a complicated reticulate evolution, consisting of seven parental taxa and 58 hybrids. The large number of hybrid taxa have added significant morphological complexity to the group leading to difficult taxonomic issues. The hybrids generally had broader ranges and more populations than their parental taxa. Genetic combination of different pairs of parental species created divergent phenotypes of hybrids, exhibited by both morphological characteristics and ecological fidelities. Niche novelty could facilitate hybrid speciation. Apogamy is common in this group and potentially contributes to the sustainability of the whole group. We propose that frequent hybridizations among members of the P. fauriei group generate and maintain genetic diversity, via novel genetic combinations, niche differentiation, and apogamy.

Homoploid hybrids, allopolyploids, and high ploidy levels characterize the evolutionary history of a western North American quillwort (Isoëtes) complex

Polyploidy and hybridization are important processes in seed-free plant evolution. However, a major gap lies in our understanding of how these processes affect the evolutionary history of high-ploidy systems. The heterosporous lycophyte genus Isoëtes is a lineage with many putative hybrids and high-level polyploid taxa (ranging from tetraploid to dodecaploid). Here, we use a complex of western North American Isoëtes, to understand the role of hybridization and high-level polyploidy in generating and maintaining novel diversity. To uncover these processes, we use restriction-site associated DNA sequencing (RADseq), multiple alleles of a single low-copy nuclear marker, whole plastomes, cytology (genome size estimates and chromosome counts), and reproductive status (fertile or sterile). With this dataset, we show that hybridization occurs easily between species in this complex and is bidirectional between identical, but not different, cytotypes. Furthermore, we show that fertile allopolyploids appear to have formed repeatedly from sterile homoploid and interploid hybrids. We propose that low prezygotic reproductive barriers and a high frequency of whole-genome duplication allow for high-level polyploid systems to generate novel lineages, and that these mechanisms may be important in shaping extant Isoëtes diversity.

Allopolyploid Speciation Accompanied by Gene Flow in a Tree Fern

Hybridization in plants may result in hybrid speciation or introgression and, thus, is now widely understood to be an important mechanism of species diversity on an evolutionary timescale. Hybridization is particularly common in ferns, as is polyploidy, which often results from hybrid crosses. Nevertheless, hybrid speciation as an evolutionary process in fern lineages remains poorly understood. Here, we employ flow cytometry, phylogeny, genomewide single nucleotide polymorphism data sets, and admixture and coalescent modeling to show that the scaly tree fern, Gymnosphaera metteniana is a naturally occurring allotetraploid species derived from hybridization between the diploids, G. denticulata and G. gigantea. Moreover, we detected ongoing gene flow between the hybrid species and its progenitors, and we found that G. gigantea and G. metteniana inhabit distinct niches, whereas climatic niches of G. denticulata and G. metteniana largely overlap. Taken together, these results suggest that either some degree of intrinsic genetic isolation between the hybrid species and its parental progenitors or ecological isolation over short distances may be playing an important role in the evolution of reproductive barriers. Historical climate change may have facilitated the origin of G. metteniana, with the timing of hybridization coinciding with a period of intensification of the East Asian monsoon during the Pliocene and Pleistocene periods in southern China. Our study of allotetraploid G. metteniana represents the first genomic-level documentation of hybrid speciation in scaly tree ferns and, thus, provides a new perspective on evolution in the lineage.

Non-adaptive evolutionary processes governed the diversification of a temperate conifer lineage after its migration into the tropics

Constructing phylogenetic relationships among closely related species is a recurrent challenge in evolutionary biology, particularly for long-lived taxa with large effective population sizes and uncomplete reproductive isolation, like conifers. Conifers further have slow evolutionary rates, which raises the question of whether adaptive or non/adaptive processes were predominantly involved when they rapidly diversified after migrating from temperate regions into the tropical mountains. Indeed, fine-scale phylogenetic relationships within several conifer genus remain under debate. Here, we studied the phylogenetic relationships of endemic firs (Abies, Pinaceae) discontinuously distributed in the montane forests from the Southwestern United States to Guatemala, and addressed several hypotheses related to adaptive and non-adaptive radiations. We derived over 80 K SNPs from genotyping by sequencing (GBS) for 45 individuals of nine Mesoamerican species to perform phylogenetic analyses. Both Maximum Likelihood and quartets-inference phylogenies resulted in a well-resolved topology, showing a single fir lineage divided in four subgroups that coincided with the main mountain ranges of Mesoamerica; thus having important taxonomic implications. Such subdivision fitted a North-South isolation by distance framework, in which non-adaptive allopatric processes seemed the rule. Interestingly, several reticulations were observed within subgroups, especially in the central-south region, which may explain past difficulties for generating infrageneric phylogenies. Further evidence for non-adaptive processes was obtained from analyses of 21 candidate-gene regions, which exhibited diminishing values of π_a/π_s and K_a/K_s with latitude, thus indicating reduced efficiency of purifying selection towards the Equator. Our study indicates that non-adaptive allopatric processes may be key generators of species diversity and endemism in the tropics.

Palynological analysis of the genus Dryopteris Adans. (Dryopteridaceae) in Argentina

The spore morphology and wall ultrastructure of Dryopteris filix-mas, D. patula and D. wallichiana from Argentina were studied using light microscope, scanning and transmission electron microscope. The study was carried out with herbarium material from Argentine institutions. Equatorial diameters, polar diameters and laesura length were measured. The spores are monoletes with a rugate ornamentation. The folds are short to long, inflate, irregular in shape and size, and varying from subglobose to elongate. The perispore surface is rugulate. The exospore of all the species analyzed is two-layered in section. Simple and branched channels are also present. The perispore is composed of two layers, the inner one forms the ornamentation and the outer covers all the outer and inner surfaces. Some abnormalities, such as globose, triangular or twisted spores were observed. The morphology and ultrastructure of the species are very similar. The differences observed are related to the length and thickness of the perispore folds. The characteristics of these spores would not provide relevant information to differentiate species or sections within the genus, but can provide information for phylogenetic studies as well as for alterations in the biological cycles.

Cryptic and extensive hybridization between ancient lineages of American crows

Most species and therefore most hybrid zones have historically been defined using phenotypic characters. However, both speciation and hybridization can occur with negligible morphological differentiation. Recently developed genomic tools provide the means to better understand cryptic speciation and hybridization. The Northwestern Crow (Corvus caurinus) and American Crow (Corvus brachyrhynchos) are continuously distributed sister taxa that lack reliable traditional characters for identification. In this first population genomic study of Northwestern and American crows, we use genomic SNPs (nuDNA) and mtDNA to investigate the degree of genetic differentiation between these crows and the extent to which they may hybridize. Our results indicate that American and Northwestern crows have distinct evolutionary histories, supported by two nuDNA ancestry clusters and two 1.1%-divergent mtDNA clades dating to the late Pleistocene, when glacial advances may have isolated crow populations in separate refugia. We document extensive hybridization, with geographic overlap of mtDNA clades and admixture of nuDNA across >900 km of western Washington and western British Columbia. This broad hybrid zone consists of late-generation hybrids and backcrosses, but not recent (e.g., F1) hybrids. Nuclear DNA and mtDNA clines had concordant widths and were both centred in southwestern British Columbia, farther north than previously postulated. Overall, our results suggest a history of reticulate evolution in American and Northwestern crows, perhaps due to recurring neutral expansion(s) from Pleistocene glacial refugia followed by lineage fusion(s). However, we do not rule out a contributing role for more recent potential drivers of hybridization, such as expansion into human-modified habitats.

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.

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.

The systematic position of Dryopteris blanfordii subsp. nigrosquamosa (Ching) Fraser-Jenkins within the genus Dryopteris Adans

Dryopteris blanfordii (C.Hope) C.Chr. is a member of the Dryopteridaceae, growing in high altitude Picea or Abies forests (2900-3500 m) in China and India. Phylogenetic relationships between D. blanfordii subsp. nigrosquamosa and closely related species of Dryopteris were investigated using a combined analysis of multiple molecular data sets (the protein-coding region of rbcL and matK genes and intergenic spacers psbA-trnH, trnP-petG, rps4-trnS, trnL-trnF and rbcL-accD). An assumption about the position of D. blanfordii subsp. nigrosquamosa within Dryopteris was made by using the Maximum Likelihood and Bayesian Inference approach and chloroplast marker sequences of Dryopteris species from GenBank. The results demonstrated that Asian taxa D. blanfordii subsp. nigrosquamosa and D. laeta as well as two American species D. arguta and D. marginalis belong to the same clade, all four of them being part of Dryopteris section Dryopteris.

Conflicting phylogenomic signals reveal a pattern of reticulate evolution in a recent high-Andean diversification (Asteraceae: Astereae: Diplostephium)

High-throughput sequencing is helping biologists to overcome the difficulties of inferring the phylogenies of recently diverged taxa. The present study analyzes the phylogenetic signal of genomic regions with different inheritance patterns using genome skimming and ddRAD-seq in a species-rich Andean genus (Diplostephium) and its allies. We analyzed the complete nuclear ribosomal cistron, the complete chloroplast genome, a partial mitochondrial genome, and a nuclear-ddRAD matrix separately with phylogenetic methods. We applied several approaches to understand the causes of incongruence among datasets, including simulations and the detection of introgression using the D-statistic (ABBA-BABA test). We found significant incongruence among the nuclear, chloroplast, and mitochondrial phylogenies. The strong signal of hybridization found by simulations and the D-statistic among genera and inside the main clades of Diplostephium indicate reticulate evolution as a main cause of phylogenetic incongruence. Our results add evidence for a major role of reticulate evolution in events of rapid diversification. Hybridization and introgression confound chloroplast and mitochondrial phylogenies in relation to the species tree as a result of the uniparental inheritance of these genomic regions. Practical implications regarding the prevalence of hybridization are discussed in relation to the phylogenetic method.

Reticulate evolution and sea-level fluctuations together drove species diversification of slipper orchids (Paphiopedilum) in South-East Asia

South-East Asia covers four of the world's biodiversity hotspots, showing high species diversity and endemism. Owing to the successive expansion and contraction of distribution and the fragmentation by geographical barriers, the tropical flora greatly diversified in this region during the Tertiary, but the evolutionary tempo and mode of species diversity remain poorly investigated. Paphiopedilum, the largest genus of slipper orchids comprising nearly 100 species, is mainly distributed in South-East Asia, providing an ideal system for exploring how plant species diversity was shaped in this region. Here, we investigated the evolutionary history of this genus with eight cpDNA regions and four low-copy nuclear genes. Discordance between gene trees and network analysis indicates that reticulate evolution occurred in the genus. Ancestral area reconstruction suggests that vicariance and long-distance dispersal together led to its current distribution. Diversification rate variation was detected and strongly correlated with the species diversity in subg. Paphiopedilum (~80 species). The shift of speciation rate in subg. Paphiopedilum was coincident with sea-level fluctuations in the late Cenozoic, which could have provided ecological opportunities for speciation and created bridges or barriers for gene flow. Moreover, some other factors (e.g. sympatric distribution, incomplete reproductive barriers and clonal propagation) might also be advantageous for the formation and reproduction of hybrid species. In conclusion, our study suggests that the interplay of reticulate evolution and sea-level fluctuations has promoted the diversification of the genus Paphiopedilum and sheds light into the evolution of Orchidaceae and the historical processes of plant species diversification in South-East Asia.

Reticulate evolution in the apogamous Dryopteris varia complex (Dryopteridaceae, subg. Erythrovariae, sect. Variae) and its related sexual species in Japan

Apogamous fern species are often difficult to distinguish from related species because of their continuous morphological variations. To clarify the genetic relationships among the members of the Dryopteris varia complex, we analyzed the nucleotide sequences of the plastid gene rbcL and the nuclear gene PgiC. We also analyzed the diploid sexual species D. caudipinna and D. chinensis, which have not been included in the complex, but were recently shown to be closely related to the complex in a molecular phylogenetic study. The PgiC sequences of the diploid sexual species, D. varia, D. saxifraga, D. sp. 'protobissetiana' (undescribed diploid sexual species), D. caudipinna, and D. chinensis, were well differentiated and hence designated A, B, C, D, and E, respectively. Thus, the PgiC constitution of apogamous species in the complex was as follows: D. bissetiana, B + C; D. kobayashii, B + C + E); D. pacifica, A + C, A + B + C, or A + C + D; D. sacrosancta, A + C + E; and D. saxifragivaria, B + C. These results suggest that these apogamous species are formed by hybridizations of species including not only the three diploid sexual species of the D. varia complex (A, B, and C) but also the two diploid sexual species D. caudipinna (D) and D. chinensis (E), which do not belong to the complex.

Phylogeny and evolutionary histories of Pyrus L. revealed by phylogenetic trees and networks based on data from multiple DNA sequences

Reconstructing the phylogeny of Pyrus has been difficult due to the wide distribution of the genus and lack of informative data. In this study, we collected 110 accessions representing 25 Pyrus species and constructed both phylogenetic trees and phylogenetic networks based on multiple DNA sequence datasets. Phylogenetic trees based on both cpDNA and nuclear LFY2int2-N (LN) data resulted in poor resolution, especially, only five primary species were monophyletic in the LN tree. A phylogenetic network of LN suggested that reticulation caused by hybridization is one of the major evolutionary processes for Pyrus species. Polytomies of the gene trees and star-like structure of cpDNA networks suggested rapid radiation is another major evolutionary process, especially for the occidental species. Pyrus calleryana and P. regelii were the earliest diverged Pyrus species. Two North African species, P. cordata, P. spinosa and P. betulaefolia were descendent of primitive stock Pyrus species and still share some common molecular characters. Southwestern China, where a large number of P. pashia populations are found, is probably the most important diversification center of Pyrus. More accessions and nuclear genes are needed for further understanding the evolutionary histories of Pyrus.

Phylogeny of a genomically diverse group of elymus (poaceae) allopolyploids reveals multiple levels of reticulation

The grass tribe Triticeae (=Hordeeae) comprises only about 300 species, but it is well known for the economically important crop plants wheat, barley, and rye. The group is also recognized as a fascinating example of evolutionary complexity, with a history shaped by numerous events of auto- and allopolyploidy and apparent introgression involving diploids and polyploids. The genus Elymus comprises a heterogeneous collection of allopolyploid genome combinations, all of which include at least one set of homoeologs, designated St, derived from Pseudoroegneria. The current analysis includes a geographically and genomically diverse collection of 21 tetraploid Elymus species, and a single hexaploid species. Diploid and polyploid relationships were estimated using four molecular data sets, including one that combines two regions of the chloroplast genome, and three from unlinked nuclear genes: phosphoenolpyruvate carboxylase, β-amylase, and granule-bound starch synthase I. Four gene trees were generated using maximum likelihood, and the phylogenetic placement of the polyploid sequences reveals extensive reticulation beyond allopolyploidy alone. The trees were interpreted with reference to numerous phenomena known to complicate allopolyploid phylogenies, and introgression was identified as a major factor in their history. The work illustrates the interpretation of complicated phylogenetic results through the sequential consideration of numerous possible explanations, and the results highlight the value of careful inspection of multiple independent molecular phylogenetic estimates, with particular focus on the differences among them.

Origins, genetic structure, and systematics of the narrow endemic peatmosses (Sphagnum): S. guwassanense and S. triseriporum (Sphagnaceae)

PREMISE OF THE STUDY

Sphagnum dominates vast expanses of wetland habitats throughout the northern hemisphere and species delimitation within the genus is important because floristic changes associated with a warming global climate may have measureable impacts on large-scale ecological processes. Most northern hemisphere peatmoss species (Sphagnum) have circumboreal ranges, but the Japanese species generally known as S. calymmatophyllum is endemic to Honshu Island. This prompted a population genetic and phylogenetic analysis to resolve the origin(s), population structure, and phylogenetic relationships of this morphologically variable species. •

METHODS

Sixty plants collected from Mt. Gassan and Mt. Hakkoda were genotyped for 12 microsatellite loci. Two plastid loci and three anonymous nuclear loci were sequenced in a subset of the plants, plus representatives from 10 closely related species. •

KEY RESULTS

Gametophytes exhibited fixed or nearly fixed heterozygosity at 9-10 of the 12 microsatellite loci. Two genetic groups were resolved by the microsatellite data, individuals showed no evidence of admixture, and the two groups of plants differ in morphology. They are heterozygous for different sets of alleles. The two taxa share plastid DNA sequences with two species that are common in Alaska. •

CONCLUSIONS

Two taxa were distinguished: S. guwassanense and S. triseriporum. Both are allopolyploids; they originated independently from different but closely related progenitors. The maternal progenitor was likely either S. orientale or S. inexspectatum. The two allopolyploid taxa are heterozygous for (different) private microsatellite alleles, and one progenitor could be extinct.

Possible natural hybridization of two morphologically distinct species of Acropora (Cnidaria, Scleractinia) in the Pacific: fertilization and larval survival rates

Natural hybridization of corals in the Indo-Pacific has been considered rather rare. However, field studies have observed many corals with intermediate interspecific or unusual morphologies. Given that the existence of F1 hybrids with intermediate interspecific morphologies has been proven in the Caribbean, hybrids may also inhabit the Indo-Pacific and occur more frequently than expected. In this study, we focused on two morphologically different species, Acropora florida and A. intermedia, and performed crossing experiments at Akajima Island, Japan. Results showed that these species could hybridize in both directions via eggs and sperm, but that fertilization rates significantly differed according to which species provided eggs. These results are similar to those reported from the Caribbean. Although all embryos developed normally to the planular larval stage, the developmental processes of some hybrid embryos were delayed by approximately 1 h compared with conspecific embryos, suggesting that fertilization occurred 1 h later in interspecific crosses than in intraspecific crosses. More successful hybridization could occur under conditions with low numbers of conspecific colonies. Additionally, a comparison of survival rates between hybrid and intraspecific larvae revealed that intra- and interspecific larvae produced from eggs of A. florida survived for significantly longer than those produced from eggs of A. intermedia. Considering these data, under specific conditions, hybrids can be expected to be produced and survive in nature in the Pacific. Furthermore, we identified one colony with intermediate morphology between A. florida and A. intermedia in the field. This colony was fertilized only by eggs of A. florida, with high fertilization rates, suggesting that this colony would be a hybrid of these two species and might be backcrossed.

Molecular circumscription and major evolutionary lineages of the fern genus Dryopteris (Dryopteridaceae)

BACKGROUND

The fern genus Dryopteris (Dryopteridaceae) is among the most common and species rich fern genera in temperate forests in the northern hemisphere containing 225-300 species worldwide. The circumscription of Dryopteris has been controversial and various related genera have, over the time, been included in and excluded from Dryopteris. The infrageneric phylogeny has largely remained unclear, and the placement of the majority of the supraspecific taxa of Dryopteris has never been tested using molecular data.

RESULTS

In this study, DNA sequences of four plastid loci (rbcL gene, rps4-trnS spacer, trnL intron, trnL-F spacer) were used to reconstruct the phylogeny of Dryopteris. A total of 122 accessions are sampled in our analysis and they represent 100 species of the expanded Dryopteris including Acrophorus, Acrorumohra, Diacalpe, Dryopsis, Nothoperanema, and Peranema. All four subgenera and 19 sections currently recognized in Dryopteris s.s. are included. One species each of Arachniodes, Leptorumohra, and Lithostegia of Dryopteridaceae are used as outgroups. Our study confirms the paraphyly of Dryopteris and provides the first strong molecular evidence on the monophyly of Acrophorus, Diacalpe, Dryopsis, Nothoperanema, and Peranema. However, all these monophyletic groups together with the paraphyletic Acrorumohra are suggested to be merged into Dryopteris based on both molecular and morphological evidence. Our analysis identified 13 well-supported monophyletic groups. Each of the 13 clades is additionally supported by morphological synapomophies and is inferred to represent a major evolutionary lineage in Dryopteris. In contrast, monophyly of the four subgenera and 15 out of 19 sections currently recognized in Dryopteris s.s is not supported by plastid data.

CONCLUSIONS

The genera, Acrophorus, Acrorumohra, Diacalpe, Dryopsis, Nothoperanema, and Peranema, should all be merged into Dryopteris. Most species of these genera share a short rhizome and catadromic arrangement of frond segments, unlike the sister genus of Dryopteris s.l., Arachniodes, which has anadromic arrangement of frond segments. The non-monophyly of the 19 out of the 21 supraspecific taxa (sections, subgenera) in Dryopteris strongly suggests that the current taxonomy of this genus is in need of revision. The disagreement between the previous taxonomy and molecular results in Dryopteris may be due partly to interspecific hybridization and polyplodization. More morphological studies and molecular data, especially from the nuclear genome, are needed to thoroughly elucidate the evolutionary history of Dryopteris. The 13 well-supported clades identified based on our data represent 13 major evolutionary lineages in Dryopteris that are also supported by morphological synapomophies.

Unraveling reticulate evolution in North American Dryopteris (Dryopteridaceae)

BACKGROUND

The thirteen species of Dryopteris in North America have long been suspected of having undergone a complicated history of reticulate evolution via allopolyploid hybridization. Various explanations for the origins of the allopolyploid taxa have been suggested, and though most lines of evidence have supported the so-called "semicristata" hypothesis, contention over the group's history has continued in several recent, conflicting studies.

RESULTS

Sequence data from nine plastid and two nuclear markers were collected from 73 accessions representing 35 species of Dryopteris. Sequences from each of the allopolyploids are most closely related to their progenitor species as predicted by the "semicristata" hypothesis. Allotetraploid D. campyloptera appears to be derived from a hybrid between diploid D. expansa and D. intermedia; D. celsa, from diploid D. ludoviciana x D. goldiana; and D. carthusiana and D. cristata, from diploid "D. semicristata" x D. intermedia and D. ludoviciana, respectively. Allohexaploid D. clintoniana appears to be derived from D. cristata x D.goldiana. The earliest estimated dates of formation of the allopolyploids, based on divergence time analyses, were within the last 6 Ma. We found no evidence for recurrent formation of any of the allopolyploids. The sexual allopolyploid taxa are derived from crosses between parents that show intermediate levels of genetic divergence relative to all pairs of potential progenitors. In addition, the four allotetraploids are transgressive with respect to geographic range relative to one or both of their parents (their ranges extend beyond those of the parents), suggesting that ecological advantages in novel habitats or regions may promote long-term regional coexistence of the hybrid taxa with their progenitors.

CONCLUSIONS

This study provides the first thorough evaluation of the North American complex of woodferns using extensive sampling of taxa and genetic markers. Phylogenies produced from each of three datasets (one plastid and two nuclear) support the "semicristata" hypothesis, including the existence of a missing diploid progenitor, and allow us to reject all competing hypotheses. This study demonstrates the value of using multiple, biparentally inherited markers to evaluate reticulate complexes, assess the frequency of recurrent polyploidization, and determine the relative importance of introgression vs. hybridization in shaping the histories of such groups.