Transcriptome-mining for single-copy nuclear markers in ferns

Insights into the evolutionary history and taxonomic status of Sinopteris (Pteridaceae)

As an endemic Chinese genus, Sinopteris C. Chr. & Ching was once considered an early diverged taxon of cheilanthoid ferns, and its taxonomic status has long been controversial. In this study, eight datasets spanning the complete chloroplast genomes and three nuclear genes were used to reconstruct the phylogeny of Sinopteris and its relatives. In addition, combining morphological analyses, divergence time estimation, and ancestral trait reconstruction, the origin and evolutionary history of Sinopteris were comprehensively discussed. Based on the complete chloroplast genome dataset, our analyses yielded a phylogram with all clades strongly supported (ML-BS = 100, BI-PP = 1.0), and the topology was almost identical to that based on the concatenated sequences of nrDNA, CRY2, and IBR3. Two species of Sinopteris were united and sister to Aleuritopteris niphobola (C. Chr.) Ching. They constituted a stable monophyletic group embedded in Aleuritopteris Fée. This was also consistent with the results of morphological analyses. Divergence time estimation indicated that the clade of Aleuritopteris and Sinopteris originated in the early Miocene (ca. 16.80 Ma) and experienced two rapid diversifications, which could coincide with environmental heterogeneity caused by the progressive uplift of the Himalayas and the intense uplift of the Hengduan Mountains. Sinopteris originated in the late Miocene (ca. 6.96 Ma), accompanied by the sharp intensifications of Asian Monsoon, and began to diversify at 2.34 Ma, following the intense uplift of the Hengduan Mountains. Ancestral character reconstruction showed that monangial sori and subsessile sporangia were clearly late derived states rather than early diverged states. Both the molecular phylogenetic and morphological analyses support the inclusion of Sinopteris in Aleuritopteris.

Establishment of an allotetraploid fern species, Lepisorus yamaokae Seriz., between two highly niche-differentiated parental species

PREMISE

The successful establishment of polyploid species is hypothesized to be promoted by niche differentiation from the parental species or by range shifts during climate oscillations. However, few studies have considered both of these factors simultaneously. We resolved the origin of a tetraploid fern, Lepisorus yamaokae, and explored a pattern of niche differentiation among the allotetraploid and parental species in past and current climates.

METHODS

We reconstructed phylogenetic trees based on plastid marker and single-copy nuclear genes to resolve the allopolyploid origin of L. yamaokae. We also evaluated climatic niche differentiation among L. yamaokae and its two parental species using species distribution models in geographic space and principal component analysis.

RESULTS

We infer that L. yamaokae had a single allotetraploid origin from L. annuifrons and L. uchiyamae. Climatic niche analyses show that the parental species currently occupy different niche spaces. The predicted distribution of the parental species at the Last Glacial Maximum (LGM) suggests more opportunities for hybridization during the LGM or during other recent temporary range shifts. Lepisorus yamaokae has a narrower niche than the additive niche of the parental species. We also observed niche conservatism in L. yamaokae.

CONCLUSIONS

Range shifts of the parental species during climatic oscillations in the Quaternary likely facilitated the formation and establishment of L. yamaokae. Further, the genetic intermediacy of L. yamaokae may have enabled a niche shift in its microenvironment, resulting in its successful establishment without a macroclimatic niche shift in L. yamaokae.

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.

PlantTribes2: Tools for comparative gene family analysis in plant genomics

Plant genome-scale resources are being generated at an increasing rate as sequencing technologies continue to improve and raw data costs continue to fall; however, the cost of downstream analyses remains large. This has resulted in a considerable range of genome assembly and annotation qualities across plant genomes due to their varying sizes, complexity, and the technology used for the assembly and annotation. To effectively work across genomes, researchers increasingly rely on comparative genomic approaches that integrate across plant community resources and data types. Such efforts have aided the genome annotation process and yielded novel insights into the evolutionary history of genomes and gene families, including complex non-model organisms. The essential tools to achieve these insights rely on gene family analysis at a genome-scale, but they are not well integrated for rapid analysis of new data, and the learning curve can be steep. Here we present PlantTribes2, a scalable, easily accessible, highly customizable, and broadly applicable gene family analysis framework with multiple entry points including user provided data. It uses objective classifications of annotated protein sequences from existing, high-quality plant genomes for comparative and evolutionary studies. PlantTribes2 can improve transcript models and then sort them, either genome-scale annotations or individual gene coding sequences, into pre-computed orthologous gene family clusters with rich functional annotation information. Then, for gene families of interest, PlantTribes2 performs downstream analyses and customizable visualizations including, (1) multiple sequence alignment, (2) gene family phylogeny, (3) estimation of synonymous and non-synonymous substitution rates among homologous sequences, and (4) inference of large-scale duplication events. We give examples of PlantTribes2 applications in functional genomic studies of economically important plant families, namely transcriptomics in the weedy Orobanchaceae and a core orthogroup analysis (CROG) in Rosaceae. PlantTribes2 is freely available for use within the main public Galaxy instance and can be downloaded from GitHub or Bioconda. Importantly, PlantTribes2 can be readily adapted for use with genomic and transcriptomic data from any kind of organism.

Transcriptome-wide SNPs for Botrychium lunaria ferns enable fine-grained analysis of ploidy and population structure

Ferns are the second most diverse group of land plants after angiosperms. Extant species occupy a wide range of habitats and contribute significantly to ecosystem functioning. Despite the importance of ferns, most taxa are poorly covered by genomic resources and within‐species studies based on high‐resolution markers are entirely lacking. The genus Botrychium belongs to the family Ophioglossaceae, which includes species with very large genomes and chromosome numbers (e.g., Ophioglossum reticulatum 2n = 1520). The genus has a cosmopolitan distribution with 35 species, half of which are polyploids. Here, we establish a transcriptome for Botrychium lunaria (L.) Sw., a diploid species with an extremely large genome of about ~19.0–23.7 Gb. We assembled 25,677 high‐quality transcripts with an average length of 1,333 bp based on deep RNA‐sequencing of a single individual. We sequenced 11 additional transcriptomes of individuals from two populations in Switzerland, including the population of the reference individual. Based on read mapping to reference transcript sequences, we identified 374,463 single nucleotide polymorphisms (SNPs) segregating among individuals for an average density of 14 SNPs per kilobase. We found that all 12 transcriptomes were most likely from diploid individuals. The transcriptome‐wide markers provided unprecedented resolution of the population genetic structure, revealing substantial variation in heterozygosity among individuals. We also constructed a phylogenomic tree of 92 taxa representing all fern orders to ascertain the placement of the genus Botrychium. High‐quality transcriptomic resources and SNP sets constitute powerful population genomic resources to investigate the ecology, and evolution of fern populations.

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.

Target Nuclear and Off-Target Plastid Hybrid Enrichment Data Inform a Range of Evolutionary Depths in the Orchid Genus Epidendrum

Universal angiosperm enrichment probe sets designed to enrich hundreds of putatively orthologous nuclear single-copy loci are increasingly being applied to infer phylogenetic relationships of different lineages of angiosperms at a range of evolutionary depths. Studies applying such probe sets have focused on testing the universality and performance of the target nuclear loci, but they have not taken advantage of off-target data from other genome compartments generated alongside the nuclear loci. Here we do so to infer phylogenetic relationships in the orchid genus Epidendrum and closely related genera of subtribe Laeliinae. Our aims are to: 1) test the technical viability of applying the plant anchored hybrid enrichment (AHE) method (Angiosperm v.1 probe kit) to our focal group, 2) mine plastid protein coding genes from off-target reads; and 3) evaluate the performance of the target nuclear and off-target plastid loci in resolving and supporting phylogenetic relationships along a range of taxonomical depths. Phylogenetic relationships were inferred from the nuclear data set through coalescent summary and site-based methods, whereas plastid loci were analyzed in a concatenated partitioned matrix under maximum likelihood. The usefulness of target and flanking non-target nuclear regions and plastid loci was assessed through the estimation of their phylogenetic informativeness. Our study successfully applied the plant AHE probe kit to Epidendrum, supporting the universality of this kit in angiosperms. Moreover, it demonstrated the feasibility of mining plastome loci from off-target reads generated with the Angiosperm v.1 probe kit to obtain additional, uniparentally inherited sequence data at no extra sequencing cost. Our analyses detected some strongly supported incongruences between nuclear and plastid data sets at shallow divergences, an indication of potential lineage sorting, hybridization, or introgression events in the group. Lastly, we found that the per site phylogenetic informativeness of the ycf1 plastid gene surpasses that of all other plastid genes and several nuclear loci, making it an excellent candidate for assessing phylogenetic relationships at medium to low taxonomic levels in orchids.

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.

Seed size evolution and biogeography of Plukenetia (Euphorbiaceae), a pantropical genus with traditionally cultivated oilseed species

BACKGROUND

Plukenetia is a small pantropical genus of lianas and vines with variably sized edible oil-rich seeds that presents an ideal system to investigate neotropical and pantropical diversification patterns and seed size evolution. We assessed the biogeography and seed evolution of Plukenetia through phylogenetic analyses of a 5069 character molecular dataset comprising five nuclear and two plastid markers for 86 terminals in subtribe Plukenetiinae (representing 20 of ~ 23 Plukenetia species). Two nuclear genes, KEA1 and TEB, were used for phylogenetic reconstruction for the first time. Our goals were: (1) produce a robust, time-dependent evolutionary framework for Plukenetia using BEAST; (2) reconstruct its biogeographical history with ancestral range estimation in BIOGEOBEARS; (3) define seed size categories; (4) identify patterns of seed size evolution using ancestral state estimation; and (5) conduct regression analyses with putative drivers of seed size using the threshold model.

RESULTS

Plukenetia was resolved into two major groups, which we refer to as the pinnately- and palmately-veined clades. Our analyses suggest Plukenetia originated in the Amazon or Atlantic Forest of Brazil during the Oligocene (28.7 Mya) and migrated/dispersed between those regions and Central America/Mexico throughout the Miocene. Trans-oceanic dispersals explain the pantropical distribution of Plukenetia, including from the Amazon to Africa in the Early Miocene (17.4 Mya), followed by Africa to Madagascar and Africa to Southeast Asia in the Late Miocene (9.4 Mya) and Pliocene (4.5 Mya), respectively. We infer a single origin of large seeds in the ancestor of Plukenetia. Seed size fits a Brownian motion model of trait evolution and is moderately to strongly associated with plant size, fruit type/dispersal syndrome, and seedling ecology. Biome shifts were not drivers of seed size, although there was a weak association with a transition to fire prone semi-arid savannas.

CONCLUSIONS

The major relationships among the species of Plukenetia are now well-resolved. Our biogeographical analyses support growing evidence that many pantropical distributions developed by periodic trans-oceanic dispersals throughout the Miocene and Pliocene. Selection on a combination of traits contributed to seed size variation, while movement between forest edge/light gap and canopy niches likely contributed to the seed size extremes in Plukenetia.

Highly Resolved Phylogenetic Relationships within Order Acipenseriformes According to Novel Nuclear Markers

Order Acipenseriformes contains 27 extant species distributed across the northern hemisphere, including so-called "living fossil" species of garfish and sturgeons. Previous studies have focused on their mitochondrial genetics and have rarely used nuclear genetic data, leaving questions as to their phylogenetic relationships. This study aimed to utilize a bioinformatics approach to screen for candidate single-copy nuclear genes, using transcriptomic data from sturgeon species and genomic data from the spotted gar, Lepisosteus oculatus. We utilized nested polymerase chain reaction (PCR) and degenerate primers to identify nuclear protein-coding (NPC) gene markers to determine phylogenetic relationships among the Acipenseriformes. We identified 193 nuclear single-copy genes, selected from 1850 candidate genes with at least one exon larger than 700 bp. Forty-three of these genes were used for primer design and development of 30 NPC markers, which were sequenced for at least 14 Acipenseriformes species. Twenty-seven NPC markers were found completely in 16 species. Gene trees according to Bayesian inference (BI) and maximum likelihood (ML) were calculated based on the 30 NPC markers (20,946 bp total). Both gene and species trees produced very similar topologies. A molecular clock model estimated the divergence time between sturgeon and paddlefish at 204.1 Mya, approximately 10% later than previous estimates based on cytochrome b data (184.4 Mya). The successful development and application of NPC markers provides a new perspective and insight for the phylogenetic relationships of Acipenseriformes. Furthermore, the newly developed nuclear markers may be useful in further studies on the conservation, evolution, and genomic biology of this group.

Phylogenetic analysis reveals the origins of tetraploid and hexaploid species in the Japanese Lepisorus thunbergianus (Polypodiaceae) complex

The Japanese Lepisorus thunbergianus complex contains diploid and tetraploid races of L. thunbergianus and a hexaploid species, L. mikawanus. Here, we performed molecular phylogenetic analysis on this complex to delimit species and to elucidate the evolutionary origins of tetraploid and hexaploid species. Chloroplast DNA (cpDNA) phylogeny supported the monophyly of the complex. Based on a single-copy nuclear gene (PgiC) tree, the tetraploid L. thunbergianus samples could be classified into two variants: an allotetraploid of hybrid origin between diploid L. thunbergianus and Japanese L. angustus and another allotetraploid of hybrid origin between diploid L. thunbergianus and an unknown diploid race of L. tosaensis. These variants can be recognized morphologically and distinguished from their parent species. Hence, here we described these allopolyploids as new species, L. nigripes and L. kuratae, respectively. The hexaploid species L. mikawanus has three types of PgiC alleles, each of which was derived from diploid L. thunbergianus, L. tosaensis, and Japanese L. angustus, while cpDNA shows that it is included in Japanese L. thunbergianus clade. Based on the cpDNA phylogeny and PgiC nucleotide sequences, we therefore concluded that L. mikawanus is an allohexaploid that originated through hybridization between tetraploid species, L. nigripes and an unknown ancestral diploid race of L. tosaensis.

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.

Low-copy nuclear markers in Isoëtes (Isoëtaceae) identified with transcriptomes

PREMISE OF THE STUDY

Few genetic markers provide phylogenetic information in closely related species of Isoëtes (Isoëtaceae). We describe the development of primers for several putative low-copy nuclear markers to resolve the phylogeny of Isoëtes, particularly in the southeastern United States.

METHODS AND RESULTS

We identified regions of interest in Isoëtes transcriptomes based on low-copy genes in other plants. Primers were designed for these regions and tested with 16 taxa of Isoëtes and one species of Lycopodium. Parts of the pgiC, gapC, and IBR3 gene regions show phylogenetic signal within the North American and Mediterranean clades of Isoëtes.

CONCLUSIONS

Transcriptome data prove useful for identification and primer design of low-copy genes. Three new markers show potential for inferring phylogenies in regional clades of Isoëtes, and possibly across the entire genus.

Organelle Genome Inheritance in Deparia Ferns (Athyriaceae, Aspleniineae, Polypodiales)

Organelle genomes of land plants are predominately inherited maternally but in some cases can also be transmitted paternally or biparentally. Compared to seed plants (>83% genera of angiosperms and >12% genera of gymnosperms), plastid genome (plastome) inheritance has only been investigated in fewer than 2% of fern genera, and mitochondrial genome (mitogenome) from only one fern genus. We developed a new and efficient method to examine plastome and mitogenome inheritance in a fern species-Deparia lancea (Athyriaceae, Aspleniineae, Polypodiales), and found that plastid and mitochondrial DNAs were transmitted from only the maternal parentage to a next generation. To further examine whether both organelle genomes have the same manner of inheritance in other Deparia ferns, we sequenced both plastid and mitochondrial DNA regions of inter-species hybrids, and performed phylogenetic analyses to identify the origins of organellar DNA. Evidence from our experiments and phylogenetic analyses support that both organelle genomes in Deparia are uniparentally and maternally inherited. Most importantly, our study provides the first report of mitogenome inheritance in eupolypod ferns, and the second one among all ferns.

A novel chloroplast gene reported for flagellate plants

PREMISE OF THE STUDY

Gene space in plant plastid genomes is well characterized and annotated, yet we discovered an unrecognized open reading frame (ORF) in the fern lineage that is conserved across flagellate plants.

METHODS

We initially detected a putative uncharacterized ORF by the existence of a highly conserved region between rps16 and matK in a series of matK alignments of leptosporangiate ferns. We mined available plastid genomes for this ORF, which we now refer to as ycf94, to infer evolutionary selection pressures and assist in functional prediction. To further examine the transcription of ycf94, we assembled the plastid genome and sequenced the transcriptome of the leptosporangiate fern Adiantum shastense Huiet & A.R. Sm.

KEY RESULTS

The ycf94 predicted protein has a distinct transmembrane domain but with no sequence homology to other proteins with known function. The nonsynonymous/synonymous substitution rate ratio of ycf94 is on par with other fern plastid protein-encoding genes, and additional homologs can be found in a few lycophyte, moss, hornwort, and liverwort plastid genomes. Homologs of ycf94 were not found in seed plants. In addition, we report a high level of RNA editing for ycf94 transcripts-a hallmark of protein-coding genes in fern plastomes.

CONCLUSIONS

The degree of sequence conservation, together with the presence of a distinct transmembrane domain and RNA-editing sites, suggests that ycf94 is a protein-coding gene of functional significance in ferns and, potentially, bryophytes and lycophytes. However, the origin and exact function of this gene require further investigation.

Insights into evolution in Andean Polystichum (Dryopteridaceae) from expanded understanding of the cytosolic phosphoglucose isomerase gene

Cytosolic phosphoglucose isomerase (pgiC) is an enzyme essential to glycolysis found universally in eukaryotes, but broad understanding of variation in the gene coding for pgiC is lacking for ferns. We used a substantially expanded representation of the gene for Andean species of the fern genus Polystichum to characterize pgiC in ferns relative to angiosperms, insects, and an amoebozoan; assess the impact of selection versus neutral evolutionary processes on pgiC; and explore evolutionary relationships of selected Andean species. The dataset of complete sequences comprised nine accessions representing seven species and one hybrid from the Andes and Serra do Mar. The aligned sequences of the full data set comprised 3376 base pairs (70% of the entire gene) including 17 exons and 15 introns from two central areas of the gene. The exons are highly conserved relative to angiosperms and retain substantial homology to insect pgiC, but intron length and structure are unique to the ferns. Average intron size is similar to angiosperms; intron number and location in insects are unlike those of the plants we considered. The introns included an array of indels and, in intron 7, an extensive microsatellite array with potential utility in analyzing population-level histories. Bayesian and maximum-parsimony analysis of 129 variable nucleotides in the Andean polystichums revealed that 59 (1.7% of the 3376 total) were phylogenetically informative; most of these united sister accessions. The phylogenetic trees for the Andean polystichums were incongruent with previously published cpDNA trees for the same taxa, likely the result of rapid evolutionary change in the introns and contrasting stability in the exons. The exons code a total of seven amino-acid substitutions. Comparison of non-synonymous to synonymous substitutions did not suggest that the pgiC gene is under selection in the Andes. Variation in pgiC including two additional accessions represented by incomplete sequences provided new insights into reticulate relationships among Andean taxa.

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

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

Applying DNA Barcodes to Identify Closely Related Species of Ferns: A Case Study of the Chinese Adiantum (Pteridaceae)

DNA barcoding is a fast-developing technique to identify species by using short and standard DNA sequences. Universal selection of DNA barcodes in ferns remains unresolved. In this study, five plastid regions (rbcL, matK, trnH-psbA, trnL-F and rps4-trnS) and eight nuclear regions (ITS, pgiC, gapC, LEAFY, ITS2, IBR3_2, DET1, and SQD1_1) were screened and evaluated in the fern genus Adiantum from China and neighboring areas. Due to low primer universality (matK) and/or the existence of multiple copies (ITS), the commonly used barcodes matK and ITS were not appropriate for Adiantum. The PCR amplification rate was extremely low in all nuclear genes except for IBR3_2. rbcL had the highest PCR amplification rate (94.33%) and sequencing success rate (90.78%), while trnH-psbA had the highest species identification rate (75%). With the consideration of discriminatory power, cost-efficiency and effort, the two-barcode combination of rbcL+ trnH-psbA seems to be the best choice for barcoding Adiantum, and perhaps basal polypod ferns in general. The nuclear IBR3_2 showed 100% PCR amplification success rate in Adiantum, however, it seemed that only diploid species could acquire clean sequences without cloning. With cloning, IBR3_2 can successfully distinguish cryptic species and hybrid species from their related species. Because hybridization and allopolyploidy are common in ferns, we argue for including a selected group of nuclear loci as barcodes, especially via the next-generation sequencing, as it is much more efficient to obtain single-copy nuclear loci without the cloning procedure.

Transcriptome sequencing reveals genome-wide variation in molecular evolutionary rate among ferns

BACKGROUND

Transcriptomics in non-model plant systems has recently reached a point where the examination of nuclear genome-wide patterns in understudied groups is an achievable reality. This progress is especially notable in evolutionary studies of ferns, for which molecular resources to date have been derived primarily from the plastid genome. Here, we utilize transcriptome data in the first genome-wide comparative study of molecular evolutionary rate in ferns. We focus on the ecologically diverse family Pteridaceae, which comprises about 10 % of fern diversity and includes the enigmatic vittarioid ferns-an epiphytic, tropical lineage known for dramatically reduced morphologies and radically elongated phylogenetic branch lengths. Using expressed sequence data for 2091 loci, we perform pairwise comparisons of molecular evolutionary rate among 12 species spanning the three largest clades in the family and ask whether previously documented heterogeneity in plastid substitution rates is reflected in their nuclear genomes. We then inquire whether variation in evolutionary rate is being shaped by genes belonging to specific functional categories and test for differential patterns of selection.

RESULTS

We find significant, genome-wide differences in evolutionary rate for vittarioid ferns relative to all other lineages within the Pteridaceae, but we recover few significant correlations between faster/slower vittarioid loci and known functional gene categories. We demonstrate that the faster rates characteristic of the vittarioid ferns are likely not driven by positive selection, nor are they unique to any particular type of nucleotide substitution.

CONCLUSIONS

Our results reinforce recently reviewed mechanisms hypothesized to shape molecular evolutionary rates in vittarioid ferns and provide novel insight into substitution rate variation both within and among fern nuclear genomes.

Unraveling the origin of the Appalachian gametophyte, Vittaria appalachiana

PREMISE OF STUDY

Ferns and lycophytes are distinct among plants in producing two free-living life stages: a long-lived sporophyte phase and a (usually) short-lived gametophyte phase. Notably, however, some species have perennial, vegetatively reproducing gametophytes. Vittaria appalachiana is one of just three species in which mature sporophytes are unknown. It has a wide range throughout the Appalachian Mountains and Plateau, where it reproduces asexually via gemmae. The origin of V. appalachiana, however, has long been a mystery, with most previous studies suggesting it may have resulted from hybridization of two closely related Vittaria species (V. graminifolia and V. lineata).

METHODS

A four-gene plastid data set including 32 samples of six Vittaria species, plus samples of five outgroup species, was analyzed to uncover phylogenetic relationships. Additional analyses of nuclear DET1 gene sequences allowed for the examination of hypotheses involving a hybrid origin for V. appalachiana.

KEY RESULTS

In the plastid phylogeny, V. appalachiana is well supported as monophyletic, but is embedded within V. graminifolia. With the exception of a single aberrant allele, this result is mirrored in the nuclear tree.

CONCLUSIONS

Through analyses of plastid and nuclear data sets, this study demonstrates that a hybrid origin for V. appalachiana is unlikely. Instead, it appears that this species emerged from within the V. graminifolia lineage. Further work is needed to fully elucidate the genetic structure within this group.

Congruent Deep Relationships in the Grape Family (Vitaceae) Based on Sequences of Chloroplast Genomes and Mitochondrial Genes via Genome Skimming

Vitaceae is well-known for having one of the most economically important fruits, i.e., the grape (Vitis vinifera). The deep phylogeny of the grape family was not resolved until a recent phylogenomic analysis of 417 nuclear genes from transcriptome data. However, it has been reported extensively that topologies based on nuclear and organellar genes may be incongruent due to differences in their evolutionary histories. Therefore, it is important to reconstruct a backbone phylogeny of the grape family using plastomes and mitochondrial genes. In this study,next-generation sequencing data sets of 27 species were obtained using genome skimming with total DNAs from silica-gel preserved tissue samples on an Illumina NextSeq 500 instrument [corrected]. Plastomes were assembled using the combination of de novo and reference genome (of V. vinifera) methods. Sixteen mitochondrial genes were also obtained via genome skimming using the reference genome of V. vinifera. Extensive phylogenetic analyses were performed using maximum likelihood and Bayesian methods. The topology based on either plastome data or mitochondrial genes is congruent with the one using hundreds of nuclear genes, indicating that the grape family did not exhibit significant reticulation at the deep level. The results showcase the power of genome skimming in capturing extensive phylogenetic data: especially from chloroplast and mitochondrial DNAs.

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

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

A genome-scale mining strategy for recovering novel rapidly-evolving nuclear single-copy genes for addressing shallow-scale phylogenetics in Hydrangea

Background

Identifying orthologous molecular markers that potentially resolve relationships at and below species level has been a major challenge in molecular phylogenetics over the past decade. Non-coding regions of nuclear low- or single-copy markers are a vast and promising source of data providing information for shallow-scale phylogenetics. Taking advantage of public transcriptome data from the One Thousand Plant Project (1KP), we developed a genome-scale mining strategy for recovering potentially orthologous single-copy markers to address low-scale phylogenetics. Our marker design targeted the amplification of intron-rich nuclear single-copy regions from genomic DNA. As a case study we used Hydrangea section Cornidia, one of the most recently diverged lineages within Hydrangeaceae (Cornales), for comparing the performance of three of these nuclear markers to other “fast” evolving plastid markers.

Results

Our data mining and filtering process retrieved 73 putative nuclear single-copy genes which are potentially useful for resolving phylogenetic relationships at a range of divergence depths within Cornales. The three assessed nuclear markers showed considerably more phylogenetic signal for shallow evolutionary depths than conventional plastid markers. Phylogenetic signal in plastid markers increased less markedly towards deeper evolutionary divergences. Potential phylogenetic noise introduced by nuclear markers was lower than their respective phylogenetic signal across all evolutionary depths. In contrast, plastid markers showed higher probabilities for introducing phylogenetic noise than signal at the deepest evolutionary divergences within the tribe Hydrangeeae (Hydrangeaceae).

Conclusions

While nuclear single-copy markers are highly informative for shallow evolutionary depths without introducing phylogenetic noise, plastid markers might be more appropriate for resolving deeper-level divergences such as the backbone relationships of the Hydrangeaceae family and deeper, at which non-coding parts of nuclear markers could potentially introduce noise due to elevated rates of evolution. The herein developed and demonstrated transcriptome based mining strategy has a great potential for the design of novel and highly informative nuclear markers for a range of plant groups and evolutionary scales.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0416-z) contains supplementary material, which is available to authorized users.

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

•

PREMISE OF THE STUDY

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

METHODS

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

KEY RESULTS

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

CONCLUSIONS

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

MarkerMiner 1.0: A new application for phylogenetic marker development using angiosperm transcriptomes

PREMISE OF THE STUDY

Targeted sequencing using next-generation sequencing (NGS) platforms offers enormous potential for plant systematics by enabling economical acquisition of multilocus data sets that can resolve difficult phylogenetic problems. However, because discovery of single-copy nuclear (SCN) loci from NGS data requires both bioinformatics skills and access to high-performance computing resources, the application of NGS data has been limited.

METHODS AND RESULTS

We developed MarkerMiner 1.0, a fully automated, open-access bioinformatic workflow and application for discovery of SCN loci in angiosperms. Our new tool identified as many as 1993 SCN loci from transcriptomic data sampled as part of four independent test cases representing marker development projects at different phylogenetic scales.

CONCLUSIONS

MarkerMiner is an easy-to-use and effective tool for discovery of putative SCN loci. It can be run locally or via the Web, and its tabular and alignment outputs facilitate efficient downstream assessments of phylogenetic utility, locus selection, intron-exon boundary prediction, and primer or probe development.

AliView: a fast and lightweight alignment viewer and editor for large datasets

SUMMARY

AliView is an alignment viewer and editor designed to meet the requirements of next-generation sequencing era phylogenetic datasets. AliView handles alignments of unlimited size in the formats most commonly used, i.e. FASTA, Phylip, Nexus, Clustal and MSF. The intuitive graphical interface makes it easy to inspect, sort, delete, merge and realign sequences as part of the manual filtering process of large datasets. AliView also works as an easy-to-use alignment editor for small as well as large datasets.

AVAILABILITY AND IMPLEMENTATION

AliView is released as open-source software under the GNU General Public License, version 3.0 (GPLv3), and is available at GitHub (www.github.com/AliView). The program is cross-platform and extensively tested on Linux, Mac OS X and Windows systems. Downloads and help are available at http://ormbunkar.se/aliview

CONTACT

anders.larsson@ebc.uu.se

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Horizontal transfer of an adaptive chimeric photoreceptor from bryophytes to ferns

Ferns are well known for their shade-dwelling habits. Their ability to thrive under low-light conditions has been linked to the evolution of a novel chimeric photoreceptor--neochrome--that fuses red-sensing phytochrome and blue-sensing phototropin modules into a single gene, thereby optimizing phototropic responses. Despite being implicated in facilitating the diversification of modern ferns, the origin of neochrome has remained a mystery. We present evidence for neochrome in hornworts (a bryophyte lineage) and demonstrate that ferns acquired neochrome from hornworts via horizontal gene transfer (HGT). Fern neochromes are nested within hornwort neochromes in our large-scale phylogenetic reconstructions of phototropin and phytochrome gene families. Divergence date estimates further support the HGT hypothesis, with fern and hornwort neochromes diverging 179 Mya, long after the split between the two plant lineages (at least 400 Mya). By analyzing the draft genome of the hornwort Anthoceros punctatus, we also discovered a previously unidentified phototropin gene that likely represents the ancestral lineage of the neochrome phototropin module. Thus, a neochrome originating in hornworts was transferred horizontally to ferns, where it may have played a significant role in the diversification of modern ferns.