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

A roadmap of phylogenomic methods for studying polyploid plant genera

Phylogenetic inference of polyploid species is the first step towards understanding their patterns of diversification. In this paper, we review the challenges and limitations of inferring species relationships of polyploid plants using traditional phylogenetic sequencing approaches, as well as the mischaracterization of the species tree from single or multiple gene trees. We provide a roadmap to infer interspecific relationships among polyploid lineages by comparing and evaluating the application of current phylogenetic, phylogenomic, transcriptomic, and whole-genome approaches using different sequencing platforms. For polyploid species tree reconstruction, we assess the following criteria: (1) the amount of prior information or tools required to capture the genetic region(s) of interest; (2) the probability of recovering homeologs for polyploid species; and (3) the time efficiency of downstream data analysis. Moreover, we discuss bioinformatic pipelines that can reconstruct networks of polyploid species relationships. In summary, although current phylogenomic approaches have improved our understanding of reticulate species relationships in polyploid-rich genera, the difficulties of recovering reliable orthologous genes and sorting all homeologous copies for allopolyploids remain a challenge. In the future, assembled long-read sequencing data will assist the recovery and identification of multiple gene copies, which can be particularly useful for reconstructing the multiple independent origins of polyploids.

PacBio high-throughput multi-locus sequencing reveals high genetic diversity in mushroom-forming fungi

Multi-locus sequence data are widely used in fungal systematic and taxonomic studies to delimit species and infer evolutionary relationships. We developed and assessed the efficacy of a multi-locus pooled sequencing method using PacBio long-read high-throughput sequencing. Samples included fresh and dried voucher specimens, cultures and archival DNA extracts of Agaricomycetes with an emphasis on the order Cantharellales. Of the 283 specimens sequenced, 93.6% successfully amplified at one or more loci with a mean of 3.3 loci amplified. Our method recovered multiple sequence variants representing alleles of rDNA loci and single copy protein-coding genes rpb1, rpb2 and tef1. Within-sample genetic variation differed by locus and taxonomic group, with the greatest genetic divergence observed among sequence variants of rpb2 and tef1 from corticioid Cantharellales. Our method is a cost-effective approach for generating accurate multi-locus sequence data coupled with recovery of alleles from polymorphic samples and multi-organism specimens. These results have important implications for understanding intra-individual genomic variation among genetic loci commonly used in species delimitation of fungi.

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.

Phylogenetic Analysis of Allotetraploid Species Using Polarized Genomic Sequences

Phylogenetic analysis of polyploid hybrid species has long posed a formidable challenge as it requires the ability to distinguish between alleles of different ancestral origins in order to disentangle their individual evolutionary history. This problem has been previously addressed by conceiving phylogenies as reticulate networks, using a two-step phasing strategy that first identifies and segregates homoeologous loci and then, during a second phasing step, assigns each gene copy to one of the subgenomes of an allopolyploid species. Here, we propose an alternative approach, one that preserves the core idea behind phasing-to produce separate nucleotide sequences that capture the reticulate evolutionary history of a polyploid-while vastly simplifying its implementation by reducing a complex multistage procedure to a single phasing step. While most current methods used for phylogenetic reconstruction of polyploid species require sequencing reads to be pre-phased using experimental or computational methods-usually an expensive, complex, and/or time-consuming endeavor-phasing executed using our algorithm is performed directly on the multiple-sequence alignment (MSA), a key change that allows for the simultaneous segregation and sorting of gene copies. We introduce the concept of genomic polarization that, when applied to an allopolyploid species, produces nucleotide sequences that capture the fraction of a polyploid genome that deviates from that of a reference sequence, usually one of the other species present in the MSA. We show that if the reference sequence is one of the parental species, the polarized polyploid sequence has a close resemblance (high pairwise sequence identity) to the second parental species. This knowledge is harnessed to build a new heuristic algorithm where, by replacing the allopolyploid genomic sequence in the MSA by its polarized version, it is possible to identify the phylogenetic position of the polyploid's ancestral parents in an iterative process. The proposed methodology can be used with long-read and short-read high-throughput sequencing data and requires only one representative individual for each species to be included in the phylogenetic analysis. In its current form, it can be used in the analysis of phylogenies containing tetraploid and diploid species. We test the newly developed method extensively using simulated data in order to evaluate its accuracy. We show empirically that the use of polarized genomic sequences allows for the correct identification of both parental species of an allotetraploid with up to 97% certainty in phylogenies with moderate levels of incomplete lineage sorting (ILS) and 87% in phylogenies containing high levels of ILS. We then apply the polarization protocol to reconstruct the reticulate histories of Arabidopsis kamchatica and Arabidopsis suecica, two allopolyploids whose ancestry has been well documented. [Allopolyploidy; Arabidopsis; genomic polarization; homoeologs; incomplete lineage sorting; phasing; polyploid phylogenetics; reticulate evolution.].

Phasing Gene Copies into Polyploid Subgenomes Using a Bayesian Phylogenetic Approach

This chapter describes the usage of homologizer to phase gene copies into polyploid subgenomes. Allopolyploids contain multiple copies of each genetic locus, where each copy potentially belongs to a different subgenome with its own distinct evolutionary history. If gene copies across different loci are incorrectly phased (i.e., assigned to the wrong subgenome), then the bifurcating tree assumption underlying multilocus phylogenetic inference and related analyses will be violated, leading to unsound results. homologizer is a highly flexible Bayesian method that uses a phylogenetic framework to infer the posterior probabilities of the phasing of gene copies into subgenomes. We describe how to prepare input data and other considerations needed to perform homologizer analyses and demonstrate how to visualize and interpret the results. We first walk through a basic example using homologizer to phase gene copies into polyploid subgenomes and then demonstrate how homologizer can be used as a hypothesis-testing tool to detect non-homeologous sequences such as hidden paralogs or allelic variation through the tools of Bayesian model comparison.

PURC Provides Improved Sequence Inference for Polyploid Phylogenetics and Other Manifestations of the Multiple-Copy Problem

Inferring the true biological sequences from amplicon mixtures remains a difficult bioinformatics problem. The traditional approach is to cluster sequencing reads by similarity thresholds and treat the consensus sequence of each cluster as an "operational taxonomic unit" (OTU). Recently, this approach has been improved by model-based methods that correct PCR and sequencing errors in order to infer "amplicon sequence variants" (ASVs). To date, ASV approaches have been used primarily in metagenomics, but they are also useful for determining homeologs in polyploid organisms. To facilitate the usage of ASV methods among polyploidy researchers, we incorporated ASV inference alongside OTU clustering in PURC v2.0, a major update to PURC (Pipeline for Untangling Reticulate Complexes). In addition, PURC v2.0 features faster demultiplexing than the original version and has been updated to be compatible with Python 3. In this chapter we present results indicating that using the ASV approach is more likely to infer the correct biological sequences in comparison to the earlier OTU-based PURC and describe how to prepare sequencing data, run PURC v2.0 under several different modes, and interpret the output.

Lack of ITS sequence homogenization in Erysimum species (Brassicaceae) with different ploidy levels

The internal transcribed spacers (ITS) exhibit concerted evolution by the fast homogenization of these sequences at the intragenomic level. However, the rate and extension of this process are unclear and might be conditioned by the number and divergence of the different ITS copies. In some cases, such as hybrid species and polyploids, ITS sequence homogenization appears incomplete, resulting in multiple haplotypes within the same organism. Here, we studied the dynamics of concerted evolution in 85 individuals of seven plant species of the genus Erysimum (Brassicaceae) with multiple ploidy levels. We estimated the rate of concerted evolution and the degree of sequence homogenization separately for ITS1 and ITS2 and whether these varied with ploidy. Our results showed incomplete sequence homogenization, especially for polyploid samples, indicating a lack of concerted evolution in these taxa. Homogenization was usually higher in ITS2 than in ITS1, suggesting that concerted evolution operates more efficiently on the former. Furthermore, the hybrid origin of several species appears to contribute to the maintenance of high haplotype diversity, regardless of the level of ploidy. These findings indicate that sequence homogenization of ITS is a dynamic and complex process that might result in varying intra- and inter-genomic diversity levels.

Adaptive responses drive the success of polyploid yellowcresses (Rorippa, Brassicaceae) in the Hengduan Mountains, a temperate biodiversity hotspot

Polyploids contribute substantially to plant evolution and biodiversity; however, the mechanisms by which they succeed are still unclear. According to the polyploid adaptation hypothesis, successful polyploids spread by repeated adaptive responses to new environments. Here, we tested this hypothesis using two tetraploid yellowcresses (Rorippa), the endemic Rorippa elata and the widespread Rorippa palustris, in the temperate biodiversity hotspot of the Hengduan Mountains. Speciation modes were resolved by phylogenetic modeling using 12 low-copy nuclear loci. Phylogeographical patterns were then examined using haplotypes phased from four plastid and ITS markers, coupled with historical niche reconstruction by ecological niche modeling. We inferred the time of hybrid origins for both species as the mid-Pleistocene, with shared glacial refugia within the southern Hengduan Mountains. Phylogeographic and ecological niche reconstruction indicated recurrent northward colonization by both species after speciation, possibly tracking denuded habitats created by glacial retreat during interglacial periods. Common garden experiment involving perennial R. elata conducted over two years revealed significant changes in fitness-related traits across source latitudes or altitudes, including latitudinal increases in survival rate and compactness of plant architecture, suggesting gradual adaptation during range expansion. These findings support the polyploid adaptation hypothesis and suggest that the spread of polyploids was aided by adaptive responses to environmental changes during the Pleistocene. Our results thus provide insight into the evolutionary success of polyploids in high-altitude environments.

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.

Recurrent allopolyploidization, Y-chromosome introgression and the evolution of sexual systems in the plant genus Mercurialis

The plant genus Mercurialis includes dioecious, monoecious and androdioecious species (where males coexist with hermaphrodites). Its diversification involved reticulate evolution via hybridization and polyploidization. The Y chromosome of the diploid species Mercurialis annua shows only mild signs of degeneration. We used sequence variation at a Y-linked locus in several species and at multiple autosomal and pseudoautosomal loci to investigate the origin and evolution of the Y chromosome across the genus. Our study provides evidence for further cases of allopolyploid speciation. It also reveals that all lineages with separate sexes (with one possible exception) share the same ancestral Y chromosome. Surprisingly, males in androdioecious populations of hexaploid M. annua carry a Y chromosome that is not derived from either of its two putative progenitor lineages but from a more distantly related perennial dioecious lineage via introgression. These results throw new light on the evolution of sexual systems and polyploidy in Mercurialis and secure it as a promising model for further study of plant sex chromosomes. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Identifying cryptic fern gametophytes using DNA barcoding: A review

Ferns and lycophytes are unique among land plants in having sporophyte (diploid) and gametophyte (haploid) generations that can grow independently of each other. While most studies of fern ecology focus on the more visible sporophytic stage, the gametophyte is critically important, as it is the sexual phase of the life cycle. Yet, fern gametophytes have long been neglected in field studies due to their small size and cryptic morphology. DNA barcoding is a powerful method that can be used to identify field-collected gametophytes to species and allow for detailed study of their ecology. Here, we review the state of DNA barcoding as applied to fern gametophytes. First, we trace the history of DNA barcoding and how it has come to be applied to fern gametophytes. Next, we summarize case studies that show how DNA barcoding has been used to better understand fern species distributions, gametophyte ecology, and community ecology. Finally, we propose avenues for future research using this powerful tool, including next-generation DNA sequencing for in-field identification of cryptic gametophytes.

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.

Phylogeny, classification, and character evolution of tribe Citharexyleae (Verbenaceae)

PREMISE

As a family of Neotropical origin and primarily Neotropical distribution, the Verbenaceae are a good but understudied system with which to understand Neotropical evolution. Tribe Citharexyleae comprises three genera: Baillonia, Citharexylum-one of the largest genera in Verbenaceae-and Rehdera. A molecular phylogenetic approach was taken to resolve intergeneric relationships in Citharexyleae and infrageneric relationships in Citharexylum. The phylogeny is used to elucidate character evolution in a widespread, morphologically diverse Neotropical genus.

METHODS

Seven plastid regions, two nuclear ribosomal spacers, and six low-copy nuclear loci were analyzed for 64 species of Citharexyleae. Phylogenetic analyses were conducted using maximum likelihood, Bayesian inference, and multispecies coalescent approaches. Habit, presence or absence of thorns, inflorescence architecture, flower color, fruit color, and geography were examined to identify diagnostic character states for clades within Citharexylum.

RESULTS

Rehdera is resolved as sister to Citharexylum, and Baillonia nested within Citharexylum. Two species, C. oleinum and C. tetramerum, are not closely related to tribe Citharexyleae, but may be related to members of tribe Duranteae instead. Seven clades within Citharexylum are inferred, each characterized by a combination of geography, fruit color and/or maturation, and inflorescence architecture. There is evidence of correlated evolution between habit, axillary inflorescences, and flower number per inflorescence. Shrubs with reduced inflorescences have evolved repeatedly.

CONCLUSIONS

A subgeneric classification for Citharexylum is proposed. Although suites of associated traits are found, character morphology has been labile throughout Citharexylum's evolutionary history. Morphological diversity may be related to adaptation to differing mesic and xeric habitats.

The diversity and community structure of symbiotic cyanobacteria in hornworts inferred from long-read amplicon sequencing

PREMISE

Nitrogen-fixing endosymbioses with cyanobacteria have evolved independently in five very different plant lineages. Expanding knowledge of these symbioses promises to improve the understanding of symbiosis evolution and broaden the toolkit for agricultural engineering to reduce artificial fertilizer use. Here we focused on hornworts, a bryophyte lineage in which all members host cyanobacteria, and investigated factors shaping the diversity of their cyanobiont communities.

METHODS

We sampled hornworts and adjacent soils in upstate New York throughout the hornwort growing season. We included all three sympatric hornwort species in the area, allowing us to directly compare partner selectivity. To profile cyanobacteria communities, we established a metabarcoding protocol targeting rbcL-X with PacBio long reads.

RESULTS

The hornwort cyanobionts detected were phylogenetically diverse, including clades that do not contain other known plant symbionts. We found significant overlap between hornwort cyanobionts and soil cyanobacteria, a pattern not previously reported in other plant-cyanobacteria symbioses. Cyanobiont communities differed between host plants only centimeters apart, but we did not detect an effect of sampling time or host species on the cyanobacterial community structure.

CONCLUSIONS

This study expands the phylogenetic diversity of known symbiotic cyanobacteria. Our analyses suggest that hornwort cyanobionts have a tight connection to the soil background, and we found no evidence that time within growing season, host species, or distance at the scale of meters strongly govern cyanobacteria community assembly. This study provides a critical foundation for further study of the ecology, evolution, and interaction dynamics of plant-cyanobacteria symbiosis.

Perspectives and Benefits of High-Throughput Long-Read Sequencing in Microbial Ecology

Short-read, high-throughput sequencing (HTS) methods have yielded numerous important insights into microbial ecology and function. Yet, in many instances short-read HTS techniques are suboptimal, for example, by providing insufficient phylogenetic resolution or low integrity of assembled genomes. Single-molecule and synthetic long-read (SLR) HTS methods have successfully ameliorated these limitations. In addition, nanopore sequencing has generated a number of unique analysis opportunities, such as rapid molecular diagnostics and direct RNA sequencing, and both Pacific Biosciences (PacBio) and nanopore sequencing support detection of epigenetic modifications. Although initially suffering from relatively low sequence quality, recent advances have greatly improved the accuracy of long-read sequencing technologies. In spite of great technological progress in recent years, the long-read HTS methods (PacBio and nanopore sequencing) are still relatively costly, require large amounts of high-quality starting material, and commonly need specific solutions in various analysis steps. Despite these challenges, long-read sequencing technologies offer high-quality, cutting-edge alternatives for testing hypotheses about microbiome structure and functioning as well as assembly of eukaryote genomes from complex environmental DNA samples.

Unveiling the Patterns of Reticulated Evolutionary Processes with Phylogenomics: Hybridization and Polyploidy in the genus Rosa

Reticulation, caused by hybridization and allopolyploidization, is considered an important and frequent phenomenon in the evolution of numerous plant lineages. Although both processes represent important driving forces of evolution, they are mostly ignored in phylogenetic studies involving a large number of species. Indeed only a scattering of methods exists to recover a comprehensive reticulated evolutionary history for a broad taxon sampling. Among these methods, comparisons of topologies obtained from plastid markers with those from a few nuclear sequences are favored, even though they restrict in-depth studies of hybridization and polyploidization. The genus Rosa encompasses c. 150 species widely distributed throughout the northern hemisphere and represents a challenging taxonomic group in which hybridization and polyploidization are prominent. Our main objective was to develop a general framework that would take patterns of reticulation into account in the study of the phylogenetic relationships among Rosa species. Using amplicon sequencing we targeted allele variation in the nuclear genome as well as haploid sequences in the chloroplast genome. We successfully recovered robust plastid and nuclear phylogenies and performed in-depth tests for several scenarios of hybridization using a maximum pseudo-likelihood approach on taxon subsets. Our diploid-first approach followed by hybrid and polyploid grafting resolved most of the evolutionary relationships among Rosa subgenera, sections, and selected species. Based on these results, we provide new directions for a future revision of the infrageneric classification in Rosa. The stepwise strategy proposed here can be used to reconstruct the phylogenetic relationships of other challenging taxonomic groups with large numbers of hybrid and polyploid taxa.

Allele Sorting as a Novel Approach to Resolving the Origin of Allotetraploids Using Hyb-Seq Data: A Case Study of the Balkan Mountain Endemic Cardamine barbaraeoides

Mountains of the Balkan Peninsula are significant biodiversity hotspots with great species richness and a large proportion of narrow endemics. Processes that have driven the evolution of the rich Balkan mountain flora, however, are still insufficiently explored and understood. Here we focus on a group of Cardamine (Brassicaceae) perennials growing in wet, mainly mountainous habitats. It comprises several Mediterranean endemics, including those restricted to the Balkan Peninsula. We used target enrichment with genome skimming (Hyb-Seq) to infer their phylogenetic relationships, and, along with genomic in situ hybridization (GISH), to resolve the origin of tetraploid Cardamine barbaraeoides endemic to the Southern Pindos Mts. (Greece). We also explored the challenges of phylogenomic analyses of polyploid species and developed a new approach of allele sorting into homeologs that allows identifying subgenomes inherited from different progenitors. We obtained a robust phylogenetic reconstruction for diploids based on 1,168 low-copy nuclear genes, which suggested both allopatric and ecological speciation events. In addition, cases of plastid-nuclear discordance, in agreement with divergent nuclear ribosomal DNA (nrDNA) copy variants in some species, indicated traces of interspecific gene flow. Our results also support biogeographic links between the Balkan and Anatolian-Caucasus regions and illustrate the contribution of the latter region to high Balkan biodiversity. An allopolyploid origin was inferred for C. barbaraeoides, which highlights the role of mountains in the Balkan Peninsula both as refugia and melting pots favoring species contacts and polyploid evolution in response to Pleistocene climate-induced range dynamics. Overall, our study demonstrates the importance of a thorough phylogenomic approach when studying the evolution of recently diverged species complexes affected by reticulation events at both diploid and polyploid levels. We emphasize the significance of retrieving allelic and homeologous variation from nuclear genes, as well as multiple nrDNA copy variants from genome skim data.

Polyploid phylogenetics

Polyploidy is a dominant feature of extant plant diversity. However, major research questions, including whether polyploidy is important to long-term evolution or is just 'evolutionary noise', remain unresolved due to difficulties associated with the generation and analysis of data from polyploid lineages. Many of these difficulties have been recently overcome, such that it is now often relatively straightforward to infer the full and often reticulate phylogenetic history of groups with recently formed polyploids. This nascent field of 'polyploid phylogenetics' allows researchers to tackle long-standing questions of polyploid macroevolution, supplies the foundation for mechanistic models of ploidy change, and provides the opportunity to include a more complete and representative sample of plant taxa in our analyses in general.

Resolving phylogeny and polyploid parentage using genus-wide genome-wide sequence data from birch trees

Numerous plant genera have a history including frequent hybridisation and polyploidisation (allopolyploidisation), which means that their phylogeny is a network of reticulate evolution that cannot be accurately depicted as a bifurcating tree with a single tip per species. The genus Betula, which contains many ecologically important tree species, is a case in point. We generated genome-wide sequence reads for 27 diploid and 36 polyploid Betula species or subspecies using restriction site associated DNA (RAD) sequences. These reads were assembled into contigs with a mean length of 675 bp. We reconstructed the evolutionary relationships among diploid Betula species using both supermatrix (concatenation) and species tree methods. We identified the closest diploid relatives of the polyploids according to the relative rates at which reads from polyploids mapped to contigs from different diploid species within a concatenated reference sequence. By mapping reads from allopolyploids to their different putative diploid relatives we assembled contigs from the putative sub-genomes of allopolyploid taxa. We used these to build new phylogenies that included allopolyploid sub-genomes as separate tips. This approach yielded a highly evidenced phylogenetic hypothesis for the genus Betula, including the complex reticulate origins of the majority of its polyploid taxa. Our phylogeny divides the genus into two well supported clades, which, interestingly, differ in their seed-wing morphology. We therefore propose to split Betula into two subgenera.

A phylogenomic approach resolves the backbone of Prunus (Rosaceae) and identifies signals of hybridization and allopolyploidy

The genus Prunus, which contains 250-400 species, has ample genomic resources for the economically important taxa in the group including cherries, peaches, and almonds. However, the backbone of Prunus, specifically the position of the racemose group relative to the solitary and corymbose groups, remains phylogenetically uncertain. Surprisingly, phylogenomic analyses to resolve relationships in the genus are lacking. Here, we assemble transcriptomes from 17 Prunus species representing four subgenera, and use existing transcriptome assemblies, to resolve key relationships in the genus using a phylogenomic approach. From the transcriptomes, we constructed 21-taxon datasets of putatively single-copy nuclear genes with 591 and 379 genes, depending on taxon-occupancy filtering. Plastome sequences were obtained or assembled for all species present in the nuclear data set. The backbone of Prunus was resolved consistently in the nuclear and chloroplast phylogenies, but we found substantial cytonuclear discord within subgenera. Our nuclear phylogeny recovered a monophyletic racemose group, contrasting with previous studies finding paraphyly that suggests repeated allopolyploidy early in the evolutionary history of the genus. However, we detected multiple species with histories consistent with hybridization and allopolyploidy, including a deep hybridization event involving subgenus Amygdalus and the Armeniaca clade in subgenus Prunus. Analyses of gene tree conflict revealed substantial discord at several nodes, including the crown node of the racemose group. Alternative gene tree topologies that conflicted with the species tree were consistent with a paraphyletic racemose group, highlighting the complex reticulated evolutionary history of this group.

TaCKX gene family, at large, is associated with thousand-grain weight and plant height in common wheat

KEY MESSAGE

We used SMRT sequencing and explored the haplotypes of TaCKX genes, linked with thousand-grain weight and plant height, and developed the functionally validated markers, which can be used in the marker-assisted breeding program. Cytokinin oxidase/dehydrogenase (CKX) enzymes catalyze the permanent degradation of cytokinins. Identification of the TaCKX alleles associated with yield traits and the development of functional markers is the first step in using these alleles in marker-assisted breeding program. To identify the alleles, we sequenced the genome fragments, containing TaCKX genes from 48 wheat genotypes, by PacBio® sequencing. Six out of 22 TaCKX genes were found polymorphic, forming 14 distinct haplotypes. Functional markers were developed and validated for all the polymorphic TaCKX genes. Four specific haplotypes, i.e., TaCKX2A_2, TaCKX4A_2, TaCKX5A_3, and TaCKX9A_2, were found significantly associated with high thousand-grain weight (TGW) and short plant height (PH) in Chinese wheat micro-core collection (MCC) and GWAS open population (GWAS-OP), whereas TaCKX1B_2 in GWAS-OP and TaCKX11A_3 in MCC were significantly associated with high TGW and short PH. The mean values of TGW and PH for cumulative favorable haplotypes from chromosome 3A, i.e., TaCKX2A_2, TaCKX4A_2, and TaCKX5A_3, were significantly higher as compared to the cumulative unfavored haplotypes, and the change was additive in manner. Frequency distribution analysis revealed that since the 1960s, the frequency of the favorable haplotypes and TGW has gradually increased in Chinese wheat cultivars. Expression profiling in the seed tissue excised at 2, 4, 6, and 8 days after anthesis depicted that the favorable haplotypes are significantly less expressive as compared to the unfavored haplotypes. We conclude that the functional markers developed in this study can be used to select the favorable haplotypes of TaCKX genes in wheat marker-assisted breeding programs.

Diet of a rare herbivore based on DNA metabarcoding of feces: Selection, seasonality, and survival

In herbivores, survival and reproduction are influenced by quality and quantity of forage, and hence, diet and foraging behavior are the foundation of an herbivore's life history strategy. Given the importance of diet to most herbivores, it is imperative that we know the species of plants they prefer, especially for herbivorous species that are at risk for extinction. However, it is often difficult to identify the diet of small herbivores because: (a) They are difficult to observe, (b) collecting stomach contents requires sacrificing animals, and (c) microhistology requires accurately identifying taxa from partially digested plant fragments and likely overemphasizes less-digestible taxa. The northern Idaho ground squirrel (Urocitellus brunneus) is federally threatened in the United States under the Endangered Species Act. We used DNA metabarcoding techniques to identify the diet of 188 squirrels at 11 study sites from fecal samples. We identified 42 families, 126 genera, and 120 species of plants in the squirrel's diet. Our use of three gene regions was beneficial because reliance on only one gene region (e.g., only trnL) would have caused us to miss >30% of the taxa in their diet. Northern Idaho ground squirrel diet differed between spring and summer, frequency of many plants in the diet differed from their frequency within their foraging areas (evidence of selective foraging), and several plant genera in their diet were associated with survival. Our results suggest that while these squirrels are generalists (they consume a wide variety of plant species), they are also selective and do not eat plants relative to availability. Consumption of particular genera such as Perideridia may be associated with higher overwinter survival.

Infraspecific diversification of the star cloak fern (Notholaena standleyi) in the deserts of the United States and Mexico

PREMISE

Not all ferns grow in moist and shaded habitats. One well-known example is Notholaena standleyi, a species that thrives in deserts of the southwestern United States and Mexico. This species exhibits several "chemotypes" that differ in farina (flavonoid exudates) color and chemistry. By integrating data from molecular phylogenetics, cytology, biochemistry, and biogeography, we circumscribed the major evolutionary lineages within N. standleyi and reconstructed their diversification histories.

METHODS

Forty-eight samples were selected from across the geographic distribution of N. standleyi. Phylogenetic relationships were inferred using four plastid and five nuclear markers. Ploidy levels were inferred using spore sizes calibrated by chromosome counts, and farina chemistry was compared using thin-layer chromatography.

RESULTS

Four clades are recognized, three of which roughly correspond to previously recognized chemotypes. The diploid clades G and Y are found in the Sonoran and Chihuahuan deserts, respectively; they are estimated to have diverged in the Pleistocene, congruent with the postulated timing of climatological events separating these two deserts. Clade P/YG is tetraploid and partially overlaps the distribution of clade Y in the eastern Chihuahuan Desert. It is apparently confined to limestone, a geologic substrate rarely occupied by members of the other clades. The cryptic (C) clade, a diploid group known only from southern Mexico and highly disjunct from the other three clades, is newly recognized here.

CONCLUSIONS

Our results reveal a complex intraspecific diversification history of N. standleyi, traceable to a variety of evolutionary drivers including classic allopatry, parapatry with or without changes in geologic substrate, and sympatric divergence through polyploidization.

A Guide to Carrying Out a Phylogenomic Target Sequence Capture Project

High-throughput DNA sequencing techniques enable time- and cost-effective sequencing of large portions of the genome. Instead of sequencing and annotating whole genomes, many phylogenetic studies focus sequencing effort on large sets of pre-selected loci, which further reduces costs and bioinformatic challenges while increasing coverage. One common approach that enriches loci before sequencing is often referred to as target sequence capture. This technique has been shown to be applicable to phylogenetic studies of greatly varying evolutionary depth. Moreover, it has proven to produce powerful, large multi-locus DNA sequence datasets suitable for phylogenetic analyses. However, target capture requires careful considerations, which may greatly affect the success of experiments. Here we provide a simple flowchart for designing phylogenomic target capture experiments. We discuss necessary decisions from the identification of target loci to the final bioinformatic processing of sequence data. We outline challenges and solutions related to the taxonomic scope, sample quality, and available genomic resources of target capture projects. We hope this review will serve as a useful roadmap for designing and carrying out successful phylogenetic target capture studies.

Phylogenomic Relationships of Diploids and the Origins of Allotetraploids in Dactylorhiza (Orchidaceae)

Disentangling phylogenetic relationships proves challenging for groups that have evolved recently, especially if there is ongoing reticulation. Although they are in most cases immediately isolated from diploid relatives, sets of sibling allopolyploids often hybridize with each other, thereby increasing the complexity of an already challenging situation. Dactylorhiza (Orchidaceae: Orchidinae) is a genus much affected by allopolyploid speciation and reticulate phylogenetic relationships. Here, we use genetic variation at tens of thousands of genomic positions to unravel the convoluted evolutionary history of Dactylorhiza. We first investigate circumscription and relationships of diploid species in the genus using coalescent and maximum likelihood methods, and then group 16 allotetraploids by maximum affiliation to their putative parental diploids, implementing a method based on genotype likelihoods. The direction of hybrid crosses is inferred for each allotetraploid using information from maternally inherited plastid RADseq loci. Starting from age estimates of parental taxa, the relative ages of these allotetraploid entities are inferred by quantifying their genetic similarity to the diploids and numbers of private alleles compared with sibling allotetraploids. Whereas northwestern Europe is dominated by young allotetraploids of postglacial origins, comparatively older allotetraploids are distributed further south, where climatic conditions remained relatively stable during the Pleistocene glaciations. Our bioinformatics approach should prove effective for the study of other naturally occurring, nonmodel, polyploid plant complexes.

TRANSNAP: a web database providing comprehensive information on Japanese pear transcriptome

Japanese pear (Pyrus pyrifolia) is a major fruit tree in the family Rosaceae and is bred for fruit production. To promote the development of breeding strategies and molecular research for Japanese pear, we sequenced the transcripts of Japanese pear variety 'Hosui'. To exhaustively collect information of total gene expression, RNA samples from various organs and stages of Japanese pear were sequenced by three technologies, single-molecule real-time (SMRT) sequencing, 454 pyrosequencing, and Sanger sequencing. Using all those reads, we determined comprehensive reference sequences of Japanese pear. Then, their protein sequences were predicted, and biological functional annotations were assigned. Finally, we developed a web database, TRANSNAP (http://plantomics.mind.meiji.ac.jp/nashi), which is the first web resource of Japanese pear omics information. This database provides highly reliable information via a user-friendly web interface: the reference sequences, gene functional annotations, and gene expression profiles from microarray experiments. In addition, based on sequence comparisons among Japanese, Chinese and European pears, similar protein sequences among the pears and species-specific proteins in Japanese pear can be quickly and efficiently identified. TRANSNAP will aid molecular research and breeding in Japanese pear, and its information is available for comparative analysis among other pear species and families.

An update and reassessment of fern and lycophyte diversity data in the Japanese Archipelago

The fern and lycophyte flora of Japan comprising 721 native taxa (including subspecies and varieties) plus 371 interspecific hybrids was reassessed using a nearly comprehensively sampled distribution map at 10 km resolution vouchered by 216,687 specimens, up-to-date cytotaxonomic information covering 74% of the taxa, and an rbcL sequence dataset covering 97.9% of the taxa. Spatial distribution of species richness and phylogenetic diversity was visualized. Apomixis was observed in 11.0% of the native taxa whose reproductive modes are known. The number of sexually reproducing polyploid taxa (n = 199) is less than sexual diploids (n = 241), and 30 of them are evidently allopolyploid, in contrast with the low number of possible autopolyploids (n = 4). Apomictic taxa were found to have smaller latitudinal ranges than sexual taxa or taxa with multiple reproductive modes. A morphological character dataset in Lucid format is provided for taxonomic identification of the native taxa.

Diversification and independent polyploid origins in the disjunct species Alyssum repens from the Southeastern Alps and the Carpathians

PREMISE

Disjunct distributions have been commonly observed in mountain plant species and have stimulated phylogeographic and phylogenetic research. Here we studied Alyssum repens, a member of the polyploid species complex A. montanum-A. repens, which exhibits SE Alpine-Carpathian disjunctions with a large elevational span and consists of diploid and tetraploid populations. We aimed to investigate the species' genetic and cytotype structure in the context of its distribution patterns, to elucidate the polyploid origins and to propose an appropriate taxonomic treatment.

METHODS

We combined AFLP fingerprinting markers, sequence variation of the highly repetitive ITS region of rDNA and the low-copy DET1 nuclear gene, genome size, and morphometric data.

RESULTS

We identified four geographically structured genetic lineages. One consisted of diploid populations from the foothills of the Southeastern Alps and neighboring regions, and the three others were allopatric montane to alpine groups comprising diploids and tetraploids growing in the Southeastern Carpathians and the Apuseni Mts. in Romania.

CONCLUSIONS

We inferred a vicariance scenario associated with Quaternary climatic oscillations, accompanied by one auto- and two allopolyploidization events most likely involving a northern Balkan relative. Whereas genetic differentiation and allopatric distribution would favor the taxonomic splitting of this species, the genetic lineages largely lack morphological distinguishability, and their ecological, cytotype and genome size divergence is only partial. Even though we probably face here a case of incipient speciation, we propose to maintain the current taxonomic treatment of Alyssum repens as a single, albeit variable, species.

An update and reassessment of fern and lycophyte diversity data in the Japanese Archipelago

The fern and lycophyte flora of Japan comprising 721 native taxa (including subspecies and varieties) plus 371 interspecific hybrids was reassessed using a nearly comprehensively sampled distribution map at 10 km resolution vouchered by 216,687 specimens, up-to-date cytotaxonomic information covering 74% of the taxa, and an rbcL sequence dataset covering 97.9% of the taxa. Spatial distribution of species richness and phylogenetic diversity was visualized. Apomixis was observed in 11.0% of the native taxa whose reproductive modes are known. The number of sexually reproducing polyploid taxa (n = 199) is less than sexual diploids (n = 241), and 30 of them are evidently allopolyploid, in contrast with the low number of possible autopolyploids (n = 4). Apomictic taxa were found to have smaller latitudinal ranges than sexual taxa or taxa with multiple reproductive modes. A morphological character dataset in Lucid format is provided for taxonomic identification of the native taxa.

Worldwide relationships in the fern genus Pteridium (bracken) based on nuclear genome markers

PREMISE

Spore-bearing plants are capable of dispersing very long distances. However, it is not known if gene flow can prevent genetic divergence in widely distributed taxa. Here we address this issue, and examine systematic relationships at a global geographic scale for the fern genus Pteridium.

METHODS

We sampled plants from 100 localities worldwide, and generated nucleotide data from four nuclear genes and two plastid regions. We also examined 2801 single nucleotide polymorphisms detected by a restriction site-associated DNA approach.

RESULTS

We found evidence for two distinct diploid species and two allotetraploids between them. The "northern" species (Pteridium aquilinum) has distinct groups at the continental scale (Europe, Asia, Africa, and North America). The northern European subspecies pinetorum appears to involve admixture among all of these. A sample from the Hawaiian Islands contained elements of both North American and Asian P. aquilinum. The "southern" species, P. esculentum, shows little genetic differentiation between South American and Australian samples. Components of African genotypes are detected on all continents.

CONCLUSIONS

We find evidence of distinct continental-scale genetic differentiation in Pteridium. However, on top of this is a clear signal of recent hybridization. Thus, spore-bearing plants are clearly capable of extensive long-distance gene flow; yet appear to have differentiated genetically at the continental scale. Either gene flow in the past was at a reduced level, or vicariance is possible even in the face of long-distance gene flow.

Phylogeny of Paullinia L. (Paullinieae: Sapindaceae), a diverse genus of lianas with dynamic fruit evolution

Paullinia L. is a genus of c. 220 mostly Neotropical forest-dwelling lianas that display a wide diversity of fruit morphologies. Paullinia resembles other members of the Paullinieae tribe in being a climber with stipulate compound leaves and paired inflorescence tendrils. However, it is distinct in having capsular fruits with woody, coriaceous, or crustaceous pericarps. While consistent in this basic plan, the pericarps of Paullinia fruits are otherwise highly variable-in some species they are winged, whereas in others they are without wings or covered with spines. With the exception of the water-dispersed indehiscent spiny fruits of some members of Paullinia sect. Castanella, all species are dehiscent, opening their capsules while they are still attached to the branch, to reveal arillate animal-dispersed seeds. Here we present a molecular phylogeny of Paullinia derived from 11 molecular markers, including nine newly-developed single-copy nuclear markers amplified by microfluidics PCR. This is the first broadly sampled molecular phylogeny for the genus. Paullinia is supported as monophyletic and is sister to Cardiospermum L., which together are sister to Serjania Mill + Urvillea Kunth. We apply this novel phylogenetic hypothesis to test previous infrageneric classifications and to infer that unwinged fruits represent the ancestral condition, from which there were repeated evolutionary transitions and reversals. However, because the seeds of both winged and unwinged fruits are dispersed by animals, we conclude that the repeated transitions in fruit morphology may relate to visual display strategies to attract animal dispersers, and do not represent transitions to wind dispersal.

Widespread co-occurrence of multiple ploidy levels in fragile ferns (Cystopteris fragilis complex; Cystopteridaceae) probably stems from similar ecology of cytotypes, their efficient dispersal and inter-ploidy hybridization

BACKGROUND AND AIMS

Polyploidy has played an important role in the evolution of ferns. However, the dearth of data on cytotype diversity, cytotype distribution patterns and ecology in ferns is striking in comparison with angiosperms and prevents an assessment of whether cytotype coexistence and its mechanisms show similar patterns in both plant groups. Here, an attempt to fill this gap was made using the ploidy-variable and widely distributed Cystopteris fragilis complex.

METHODS

Flow cytometry was used to assess DNA ploidy level and monoploid genome size (Cx value) of 5518 C. fragilis individuals from 449 populations collected over most of the species' global distributional range, supplemented with data from 405 individuals representing other related species from the complex. Ecological preferences of C. fragilis tetraploids and hexaploids were compared using field-recorded parameters and database-extracted climate data.

KEY RESULTS

Altogether, five different ploidy levels (2x, 4x, 5x, 6x, 8x) were detected and three species exhibited intraspecific ploidy-level variation: C. fragilis, C. alpina and C. diaphana. Two predominant C. fragilis cytotypes, tetraploids and hexaploids, co-occur over most of Europe in a diffuse, mosaic-like pattern. Within this contact zone, 40 % of populations were mixed-ploidy and most also contained pentaploid hybrids. Environmental conditions had only a limited effect on the distribution of cytotypes. Differences were found in the Cx value of tetraploids and hexaploids: between-cytotype divergence was higher in uniform-ploidy than in mixed-ploidy populations.

CONCLUSIONS

High ploidy-level diversity and widespread cytotype coexistence in the C. fragilis complex match the well-documented patterns in some angiosperms. While ploidy coexistence in C. fragilis is not driven by environmental factors, it could be facilitated by the perennial life-form of the species, its reproductive modes and efficient wind dispersal of spores. Independent origins of hexaploids and/or inter-ploidy gene flow may be expected in mixed-ploidy populations according to Cx value comparisons.

Extensive allopolyploidy in the neotropical genus Lachemilla (Rosaceae) revealed by PCR-based target enrichment of the nuclear ribosomal DNA cistron and plastid phylogenomics

PREMISE OF THE STUDY

Polyploidy has been long recognized as an important force in plant evolution. Previous studies had suggested widespread occurrence of polyploidy and the allopolyploid origin of several species in the diverse neotropical genus Lachemilla (Rosaceae). Nonetheless, this evidence has relied mostly on patterns of cytonuclear discordance, and direct evidence from nuclear allelic markers is still needed.

METHODS

Here we used PCR target enrichment in combination with high throughput sequencing to obtain multiple copies of the nuclear ribosomal (nr) DNA cistron and 45 regions of the plastid genome (cpDNA) from 219 accessions representing 48 species of Lachemilla and to explore the allopolyploid origin of species in this group.

KEY RESULTS

We were able to identify multiple nrDNA ribotypes and establish clear evidence of allopolyploidy in 33 species of Lachemilla, showing that this condition is common and widespread in the genus. Additionally, we found evidence for three autopolyploid species. We also established multiple, independent origins of several allopolyploid species. Finally, based solely on the cpDNA phylogeny, we identified that the monotypic genus Farinopsis is the sister group of Lachemilla and allied genera within subtribe Fragariinae.

CONCLUSIONS

Our study demonstrates the utility of the nuclear ribosomal DNA cistron to detect allopolyploidy when concerted evolution of this region is not complete. Additionally, with a robust chloroplast phylogeny in place, the direction of hybridization events can be established, and multiple, independent origins of allopolyploid species can be identified.

Complete sequence of kenaf (Hibiscus cannabinus) mitochondrial genome and comparative analysis with the mitochondrial genomes of other plants

Plant mitochondrial (mt) genomes are species specific due to the vast of foreign DNA migration and frequent recombination of repeated sequences. Sequencing of the mt genome of kenaf (Hibiscus cannabinus) is essential for elucidating its evolutionary characteristics. In the present study, single-molecule real-time sequencing technology (SMRT) was used to sequence the complete mt genome of kenaf. Results showed that the complete kenaf mt genome was 569,915 bp long and consisted of 62 genes, including 36 protein-coding, 3 rRNA and 23 tRNA genes. Twenty-five introns were found among nine of the 36 protein-coding genes, and five introns were trans-spliced. A comparative analysis with other plant mt genomes showed that four syntenic gene clusters were conserved in all plant mtDNAs. Fifteen chloroplast-derived fragments were strongly associated with mt genes, including the intact sequences of the chloroplast genes psaA, ndhB and rps7. According to the plant mt genome evolution analysis, some ribosomal protein genes and succinate dehydrogenase genes were frequently lost during the evolution of angiosperms. Our data suggest that the kenaf mt genome retained evolutionarily conserved characteristics. Overall, the complete sequencing of the kenaf mt genome provides additional information and enhances our better understanding of mt genomic evolution across angiosperms.

Phylogenomic analyses reveal a deep history of hybridization and polyploidy in the Neotropical genus Lachemilla (Rosaceae)

Hybridization, incomplete lineage sorting, and phylogenetic error produce similar incongruence patterns, representing a great challenge for phylogenetic reconstruction. Here, we use sequence capture data and multiple species tree and species network approaches to resolve the backbone phylogeny of the Neotropical genus Lachemilla, while distinguishing among sources of incongruence. We used 396 nuclear loci and nearly complete plastome sequences from 27 species to clarify the relationships among the major groups of Lachemilla, and explored multiple sources of conflict between gene trees and species trees inferred with a plurality of approaches. All phylogenetic methods recovered the four major groups previously proposed for Lachemilla, but species tree methods recovered different topologies for relationships between these four clades. Species network analyses revealed that one major clade, Orbiculate, is likely of ancient hybrid origin, representing one of the main sources of incongruence among the species trees. Additionally, we found evidence for a potential whole genome duplication event shared by Lachemilla and allied genera. Lachemilla shows clear evidence of ancient and recent hybridization throughout the evolutionary history of the group. Also, we show the necessity to use phylogenetic network approaches that can simultaneously accommodate incomplete lineage sorting and gene flow when studying groups that show patterns of reticulation.

Fluidigm2PURC: Automated processing and haplotype inference for double-barcoded PCR amplicons

PREMISE OF THE STUDY

Targeted enrichment strategies for phylogenomic inference are a time- and cost-efficient way to collect DNA sequence data for large numbers of individuals at multiple, independent loci. Automated and reproducible processing of these data is a crucial step for researchers conducting phylogenetic studies.

METHODS AND RESULTS

We present Fluidigm2PURC, an open source Python utility for processing paired-end Illumina data from double-barcoded PCR amplicons. In combination with the program PURC (Pipeline for Untangling Reticulate Complexes), our scripts process raw FASTQ files for analysis with PURC and use its output to infer haplotypes for diploids, polyploids, and samples with unknown ploidy. We demonstrate the use of the pipeline with an example data set from the genus Thalictrum (Ranunculaceae).

CONCLUSIONS

Fluidigm2PURC is freely available for Unix-like operating systems on GitHub (https://github.com/pblischak/fluidigm2purc) and for all operating systems through Docker (https://hub.docker.com/r/pblischak/fluidigm2purc).

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.

Tools for Genetic Studies in Experimental Populations of Polyploids

Polyploid organisms carry more than two copies of each chromosome, a condition rarely tolerated in animals but which occurs relatively frequently in the plant kingdom. One of the principal challenges faced by polyploid organisms is to evolve stable meiotic mechanisms to faithfully transmit genetic information to the next generation upon which the study of inheritance is based. In this review we look at the tools available to the research community to better understand polyploid inheritance, many of which have only recently been developed. Most of these tools are intended for experimental populations (rather than natural populations), facilitating genomics-assisted crop improvement and plant breeding. This is hardly surprising given that a large proportion of domesticated plant species are polyploid. We focus on three main areas: (1) polyploid genotyping; (2) genetic and physical mapping; and (3) quantitative trait analysis and genomic selection. We also briefly review some miscellaneous topics such as the mode of inheritance and the availability of polyploid simulation software. The current polyploid analytic toolbox includes software for assigning marker genotypes (and in particular, estimating the dosage of marker alleles in the heterozygous condition), establishing chromosome-scale linkage phase among marker alleles, constructing (short-range) haplotypes, generating linkage maps, performing genome-wide association studies (GWAS) and quantitative trait locus (QTL) analyses, and simulating polyploid populations. These tools can also help elucidate the mode of inheritance (disomic, polysomic or a mixture of both as in segmental allopolyploids) or reveal whether double reduction and multivalent chromosomal pairing occur. An increasing number of polyploids (or associated diploids) are being sequenced, leading to publicly available reference genome assemblies. Much work remains in order to keep pace with developments in genomic technologies. However, such technologies also offer the promise of understanding polyploid genomes at a level which hitherto has remained elusive.

Practical considerations for plant phylogenomics

The past decade has seen a major breakthrough in our ability to easily and inexpensively sequence genome-scale data from diverse lineages. The development of high-throughput sequencing and long-read technologies has ushered in the era of phylogenomics, where hundreds to thousands of nuclear genes and whole organellar genomes are routinely used to reconstruct evolutionary relationships. As a result, understanding which options are best suited for a particular set of questions can be difficult, especially for those just starting in the field. Here, we review the most recent advances in plant phylogenomic methods and make recommendations for project-dependent best practices and considerations. We focus on the costs and benefits of different approaches in regard to the information they provide researchers and the questions they can address. We also highlight unique challenges and opportunities in plant systems, such as polyploidy, reticulate evolution, and the use of herbarium materials, identifying optimal methodologies for each. Finally, we draw attention to lingering challenges in the field of plant phylogenomics, such as reusability of data sets, and look at some up-and-coming technologies that may help propel the field even further.

Allele phasing has minimal impact on phylogenetic reconstruction from targeted nuclear gene sequences in a case study of Artocarpus

PREMISE OF THE STUDY

Untapped information about allele diversity within populations and individuals (i.e., heterozygosity) could improve phylogenetic resolution and accuracy. Many phylogenetic reconstructions ignore heterozygosity because it is difficult to assemble allele sequences and combine allele data across unlinked loci, and it is unclear how reconstruction methods accommodate variable sequences. We review the common methods of including heterozygosity in phylogenetic studies and present a novel method for assembling allele sequences from target-enriched Illumina sequencing libraries.

METHODS

We performed supermatrix phylogeny reconstruction and species tree estimation of Artocarpus based on three methods of accounting for heterozygous sequences: a consensus method based on de novo sequence assembly, the use of ambiguity characters, and a novel method for incorporating read information to phase alleles. We characterize the extent to which highly heterozygous sequences impeded phylogeny reconstruction and determine whether the use of allele sequences improves phylogenetic resolution or decreases topological uncertainty.

KEY RESULTS

We show here that it is possible to infer phased alleles from target-enriched Illumina libraries. We find that highly heterozygous sequences do not contribute disproportionately to poor phylogenetic resolution and that the use of allele sequences for phylogeny reconstruction does not have a clear effect on phylogenetic resolution or topological consistency.

CONCLUSIONS

We provide a framework for inferring phased alleles from target enrichment data and for assessing the contribution of allelic diversity to phylogenetic reconstruction. In our data set, the impact of allele phasing on phylogeny is minimal compared to the impact of using phylogenetic reconstruction methods that account for gene tree incongruence.

Gene expression atlas of fruit ripening and transcriptome assembly from RNA-seq data in octoploid strawberry (Fragaria × ananassa)

RNA-seq has been used to perform global expression analysis of the achene and the receptacle at four stages of fruit ripening, and of the roots and leaves of strawberry (Fragaria × ananassa). About 967 million reads and 191 Gb of sequence were produced, using Illumina sequencing. Mapping the reads in the related genome of the wild diploid Fragaria vesca revealed differences between the achene and receptacle development program, and reinforced the role played by ethylene in the ripening receptacle. For the strawberry transcriptome assembly, a de novo strategy was followed, generating separate assemblies for each of the ten tissues and stages sampled. The Trinity program was used for these assemblies, resulting in over 1.4 M isoforms. Filtering by a threshold of 0.3 FPKM, and doing Blastx (E-value < 1 e-30) against the UniProt database of plants reduced the number to 472,476 isoforms. Their assembly with the MIRA program (90% homology) resulted in 26,087 contigs. From these, 91.34 percent showed high homology to Fragaria vesca genes and 87.30 percent Fragaria iinumae (BlastN E-value < 1 e-100). Mapping back the reads on the MIRA contigs identified polymorphisms at nucleotide level, using FREEBAYES, as well as estimate their relative abundance in each sample.

Can asexuality confer a short-term advantage? Investigating apparent biogeographic success in the apomictic triploid fern Myriopteris gracilis

PREMISE OF THE STUDY

Although asexual taxa are generally seen as evolutionary dead ends, asexuality appears to provide a short-term benefit in some taxa, including a wider geographic distribution compared to sexual relatives. However, this may be an illusion created by multiple, morphologically cryptic, asexual lineages, each occupying a relatively small area. In this study we investigate the role of multiple lineages in the biogeography of Myriopteris gracilis Fée (Pteridaceae), a North American apomictic triploid fern species with a particularly large range.

METHODS

Range-wide asexuality was assessed by counting spores/sporangium in 606 Myriopteris gracilis specimens from across the species range, and lineage structure was assessed with both plastid DNA sequence and Genotyping By Sequencing (GBS) SNP datasets.

KEY RESULTS

Spore counting of >600 specimens identified no sexual populations, establishing that Myriopteris gracilis is exclusively asexual. The plastid data estimated the crown age of M. gracilis at ca. 2.5 mya and identified two lineages, each largely confined to the eastern or western portions of the range. These groups were further subdivided by the GBS data, revealing at least seven asexual lineages of varying geographic distributions, each occupying a relatively small portion of the total range of M. gracilis.

CONCLUSIONS

Although maintained exclusively through asexual reproduction, the broad distribution of Myriopteris gracilis is a compilation of numerous, independently formed asexual lineages. Since no single asexual lineage occupies the full extent of the species distribution, recurrent lineage formation should be considered when evaluating the short-term benefit of asexuality in this taxon and others.