Nuclear loci developed from multiple transcriptomes yield high resolution in phylogeny of scaly tree ferns (Cyatheaceae) from China and Vietnam

Comparative and population genomics of buckwheat species reveal key determinants of flavor and fertility

Common buckwheat (Fagopyrum esculentum) is an ancient crop with a world-wide distribution. Due to its excellent nutritional quality and high economic and ecological value, common buckwheat is becoming increasingly important throughout the world. The availability of a high-quality reference genome sequence and population genomic data will accelerate the breeding of common buckwheat, but the high heterozygosity due to the outcrossing nature has greatly hindered the genome assembly. Here we report the assembly of a chromosome-scale high-quality reference genome of F. esculentum var. homotropicum, a homozygous self-pollinating variant of common buckwheat. Comparative genomics revealed that two cultivated buckwheat species, common buckwheat (F. esculentum) and Tartary buckwheat (F. tataricum), underwent metabolomic divergence and ecotype differentiation. The expansion of several gene families in common buckwheat, including FhFAR genes, is associated with its wider distribution than Tartary buckwheat. Copy number variation of genes involved in the metabolism of flavonoids is associated with the difference of rutin content between common and Tartary buckwheat. Furthermore, we present a comprehensive atlas of genomic variation based on whole-genome resequencing of 572 accessions of common buckwheat. Population and evolutionary genomics reveal genetic variation associated with environmental adaptability and floral development between Chinese and non-Chinese cultivated groups. Genome-wide association analyses of multi-year agronomic traits with the content of flavonoids revealed that Fh05G014970 is a potential major regulator of flowering period, a key agronomic trait controlling the yield of outcrossing crops, and that Fh06G015130 is a crucial gene underlying flavor-associated flavonoids. Intriguingly, we found that the gene translocation and sequence variation of FhS-ELF3 contribute to the homomorphic self-compatibility of common buckwheat. Collectively, our results elucidate the genetic basis of speciation, ecological adaptation, fertility, and unique flavor of common buckwheat, and provide new resources for future genomics-assisted breeding of this economically important crop.

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

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

Complete chloroplast genome molecular structure, comparative and phylogenetic analyses of Sphaeropteris lepifera of Cyatheaceae family: a tree fern from China

Sphaeropteris lepifera is a tree fern in the Cyatheaceae, a family that has played an important role in the evolution of plant systems. This study aimed to analyze the complete chloroplast genome of S. lepifera and compared it with previously published chloroplast genomes Cyatheaceae family. The chloroplast genome of S. lepifera comprised 162,114 bp, consisting of a large single copy (LSC) region of 86,327 bp, a small single copy (SSC) region of 27,731 bp and a pair of inverted repeats (IRa and IRb) of 24,028 bp each. The chloroplast genome encoded 129 genes, comprising 32 transfer RNAs, 8 ribosomal RNAs, and 89 protein-coding genes. Comparison of the genomes of 7 Cyatheaceae plants showed that the chloroplast genome of S. lepifera was missing the gene trnV-UAC. Expansion of the SSC region led to the difference in the chloroplast genome size of S. lepifera. Eight genes, atpI, ccsA, petA, psaB, rpl16, rpoA, rpoC1, and ycf2 have high nucleic acid diversity and can be regarded as potential molecular markers. The genes trnG-trnR and atpB were suitable for DNA barcodes between different communities of S. lepifera. The S. lepifera groups in Zhejiang Province probably diffused from Pingtan and Ningde, Fujian. The results will provide a basis for species identification, biological studies, and endangerment mechanism of S. lepifera.

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.

The flying spider-monkey tree fern genome provides insights into fern evolution and arborescence

To date, little is known about the evolution of fern genomes, with only two small genomes published from the heterosporous Salviniales. Here we assembled the genome of Alsophila spinulosa, known as the flying spider-monkey tree fern, onto 69 pseudochromosomes. The remarkable preservation of synteny, despite resulting from an ancient whole-genome duplication over 100 million years ago, is unprecedented in plants and probably speaks to the uniqueness of tree ferns. Our detailed investigations into stem anatomy and lignin biosynthesis shed new light on the evolution of stem formation in tree ferns. We identified a phenolic compound, alsophilin, that is abundant in xylem, and we provided the molecular basis for its biosynthesis. Finally, analysis of demographic history revealed two genetic bottlenecks, resulting in rapid demographic declines of A. spinulosa. The A. spinulosa genome fills a crucial gap in the plant genomic landscape and helps elucidate many unique aspects of tree fern biology.

First Multi-Organ Full-Length Transcriptome of Tree Fern Alsophila spinulosa Highlights the Stress-Resistant and Light-Adapted Genes

Alsophila spinulosa, a relict tree fern, is a valuable plant for investigating environmental adaptations. Its genetic resources, however, are scarce. We used the PacBio and Illumina platforms to sequence the polyadenylated RNA of A. spinulosa root, rachis, and pinna, yielding 125,758, 89,107, and 89,332 unigenes, respectively. Combining the unigenes from three organs yielded a non-redundant reference transcriptome with 278,357 unigenes and N50 of 4141 bp, which were further reconstructed into 38,470 UniTransModels. According to functional annotation, pentatricopeptide repeat genes and retrotransposon-encoded polyprotein genes are the most abundant unigenes. Clean reads mapping to the full-length transcriptome is used to assess the expression of unigenes. The stress-induced ASR genes are highly expressed in all three organs, which is validated by qRT-PCR. The organ-specific upregulated genes are enriched for pathways involved in stress response, secondary metabolites, and photosynthesis. Genes for five types of photoreceptors, CRY signaling pathway, ABA biosynthesis and transduction pathway, and stomatal movement-related ion channel/transporter are profiled using the high-quality unigenes. The gene expression pattern coincides with the previously identified stomatal characteristics of fern. This study is the first multi-organ full-length transcriptome report of a tree fern species, the abundant genetic resources and comprehensive analysis of A. spinulosa, which provides the groundwork for future tree fern research.

Allopolyploid Speciation Accompanied by Gene Flow in a Tree Fern

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

Inferring the Phylogenetic Positions of Two Fig Wasp Subfamilies of Epichrysomallinae and Sycophaginae Using Transcriptomes and Mitochondrial Data

Fig wasps are a group of insects (Hymenoptera: Chalcidoidea) that live in the compact syconia of fig trees (Moraceae: Ficus). Accurate classification and phylogenetic results are very important for studies of fig wasps, but the taxonomic statuses of some fig wasps, especially the non-pollinating subfamilies are difficult to determine, such as Epichrysomallinae and Sycophaginae. To resolve the taxonomic statuses of Epichrysomallinae and Sycophaginae, we obtained transcriptomes and mitochondrial genome (mitogenome) data for four species of fig wasps. These newly added data were combined with the data of 13 wasps (data on 11 fig wasp species were from our laboratory and two wasp species were download from NCBI). Based on the transcriptome and genome data, we obtained 145 single-copy orthologous (SCO) genes in 17 wasp species, and based on mitogenome data, we obtained 13 mitochondrial protein-coding genes (PCGs) for each of the 17 wasp species. Ultimately, we used 145 SCO genes, 13 mitochondrial PCGs and combined SCO genes and mitochondrial genes data to reconstruct the phylogenies of fig wasps using both maximum likelihood (ML) and Bayesian inference (BI) analyses. Our results suggest that both Epichrysomallinae and Sycophaginae are more closely related to Agaonidae with a high statistical support.