Using nuclear loci and allelic variation to disentangle the phylogeny of Phyllostachys (Poaceae, Bambusoideae)

Testing complete plastomes and nuclear ribosomal DNA sequences for species identification in a taxonomically difficult bamboo genus Fargesia

Fargesia, the largest genus within the temperate bamboo tribe Arundinarieae, has more than 90 species mainly distributed in the mountains of Southwest China. The Fargesia bamboos are important components of the subalpine forest ecosystems that provide food and habitat for many endangered animals, including the giant panda. However, species-level identification of Fargesia is difficult. Moreover, the rapid radiation and slow molecular evolutionary rate of Fargesia pose a significant challenge to using DNA barcoding with standard plant barcodes (rbcL, matK, and ITS) in bamboos. With progress in the sequencing technologies, complete plastid genomes (plastomes) and nuclear ribosomal DNA (nrDNA) sequences have been proposed as organelle barcodes for species identification; however, these have not been tested in bamboos. We collected 196 individuals representing 62 species of Fargesia to comprehensively evaluate the discriminatory power of plastomes and nrDNA sequences compared to standard barcodes. Our analysis indicates that complete plastomes have substantially higher discriminatory power (28.6%) than standard barcodes (5.7%), whereas nrDNA sequences show a moderate improvement (65.4%) compared to ITS (47.2%). We also found that nuclear markers performed better than plastid markers, and ITS alone had higher discriminatory power than complete plastomes. The study also demonstrated that plastomes and nrDNA sequences can contribute to intrageneric phylogenetic resolution in Fargesia. However, neither of these sequences were able to discriminate all the sampled species, and therefore, more nuclear markers need to be identified.

Phylogenomics and the flowering plant tree of life

The advances accelerated by next-generation sequencing and long-read sequencing technologies continue to provide an impetus for plant phylogenetic study. In the past decade, a large number of phylogenetic studies adopting hundreds to thousands of genes across a wealth of clades have emerged and ushered plant phylogenetics and evolution into a new era. In the meantime, a roadmap for researchers when making decisions across different approaches for their phylogenomic research design is imminent. This review focuses on the utility of genomic data (from organelle genomes, to both reduced representation sequencing and whole-genome sequencing) in phylogenetic and evolutionary investigations, describes the baseline methodology of experimental and analytical procedures, and summarizes recent progress in flowering plant phylogenomics at the ordinal, familial, tribal, and lower levels. We also discuss the challenges, such as the adverse impact on orthology inference and phylogenetic reconstruction raised from systematic errors, and underlying biological factors, such as whole-genome duplication, hybridization/introgression, and incomplete lineage sorting, together suggesting that a bifurcating tree may not be the best model for the tree of life. Finally, we discuss promising avenues for future plant phylogenomic studies.

Andean uplift, drainage basin formation, and the evolution of plants living in fast-flowing aquatic ecosystems in northern South America

Northern South America is a geologically dynamic and species-rich region. Fossil and stratigraphic data show that mountain uplift in the tropical Andes reconfigured river drainages. These landscape changes shaped the evolution of the flora in the region, yet the impacts on aquatic taxa have been overlooked. We explore the role of landscape change on the evolution of plants living strictly in rivers across drainage basins in northern South America by conducting population structure, phylogenetic inference, and divergence-dating analyses for two species in the genus Marathrum (Podostemaceae). Mountain uplift and drainage basin formation isolated populations of M. utile and M. foeniculaceum in northern South America and created barriers to gene flow across river drainages. Sympatric species hybridize and the hybrids show the phenotype of one parental line. We propose that the pattern of divergence of populations reflects the formation of river drainages, which was not complete until < 4.1 million yr ago (Ma). Our study provides a clear picture of the role of landscape change on the evolution of plants living strictly in rivers in northern South America. By shifting the focus to aquatic taxa, we provide a novel perspective on the processes shaping the evolution of the Neotropical flora.

Analysis of 427 genomes reveals moso bamboo population structure and genetic basis of property traits

Moso bamboo (Phyllostachys edulis) is an economically and ecologically important nontimber forestry species. Further development of this species as a sustainable bamboo resource has been hindered by a lack of population genome information. Here, we report a moso bamboo genomic variation atlas of 5.45 million single-nucleotide polymorphisms (SNPs) from whole-genome resequencing of 427 individuals covering 15 representative geographic areas. We uncover low genetic diversity, high genotype heterozygosity, and genes under balancing selection underlying moso bamboo population adaptation. We infer its demographic history with one bottleneck and its recently small population without a rebound. We define five phylogenetic groups and infer that one group probably originated by a single-origin event from East China. Finally, we conduct genome-wide association analysis of nine important property-related traits to identify candidate genes, many of which are involved in cell wall, carbohydrate metabolism, and environmental adaptation. These results provide a foundation and resources for understanding moso bamboo evolution and the genetic mechanisms of agriculturally important traits.

The complete chloroplast genome of a solid type of Phyllostachys nidularia (Bambusoideae: Poaceae), a species endemic to China

Phyllostachys nidularia (Bambusoideae: Poaceae), widely distributed in the Yangtze River Basin and various provinces (regions) in southern China, is one of the most important small and medium-sized bamboo species used in both bamboo shoots and timber. In the present study, we assembled a complete chloroplast genome of the economically important bamboo form Phyllostachys nidularia f. farcta H.R.Zhao & A.T.Liu using whole genome sequencing data previously reported. The complete chloroplast (cp) genome is 139,706 bp in length. A total of 129 unique genes were annotated, including 82 protein-coding, 39 tRNA, and eight rRNA genes. Phylogenetic analysis results supported that P. nidularia f. farcta was closely related to Phyllostachys reticulata. This work would help us better understand the evolution of the Phyllostachys cp genome.

Targeted enrichment of novel chloroplast-based probes reveals a large-scale phylogeny of 412 bamboos

BACKGROUND

The subfamily Bambusoideae belongs to the grass family Poaceae and has significant roles in culture, economy, and ecology. However, the phylogenetic relationships based on large-scale chloroplast genomes (CpGenomes) were elusive. Moreover, most of the chloroplast DNA sequencing methods cannot meet the requirements of large-scale CpGenome sequencing, which greatly limits and impedes the in-depth research of plant genetics and evolution.

RESULTS

To develop a set of bamboo probes, we used 99 high-quality CpGenomes with 6 bamboo CpGenomes as representative species for the probe design, and assembled 15 M unique sequences as the final pan-chloroplast genome. A total of 180,519 probes for chloroplast DNA fragments were designed and synthesized by a novel hybridization-based targeted enrichment approach. Another 468 CpGenomes were selected as test data to verify the quality of the newly synthesized probes and the efficiency of the probes for chloroplast capture. We then successfully applied the probes to synthesize, enrich, and assemble 358 non-redundant CpGenomes of woody bamboo in China. Evaluation analysis showed the probes may be applicable to chloroplasts in Magnoliales, Pinales, Poales et al. Moreover, we reconstructed a phylogenetic tree of 412 bamboos (358 in-house and 54 published), supporting a non-monophyletic lineage of the genus Phyllostachys. Additionally, we shared our data by uploading a dataset of bamboo CpGenome into CNGB ( https://db.cngb.org/search/project/CNP0000502/ ) to enrich resources and promote the development of bamboo phylogenetics.

CONCLUSIONS

The development of the CpGenome enrichment pipeline and its performance on bamboos recommended an inexpensive, high-throughput, time-saving and efficient CpGenome sequencing strategy, which can be applied to facilitate the phylogenetics analysis of most green plants.

Parallel ddRAD and Genome Skimming Analyses Reveal a Radiative and Reticulate Evolutionary History of the Temperate Bamboos

Rapid evolutionary radiations are among the most challenging phylogenetic problems, wherein different types of data (e.g., morphology and molecular) or genetic markers (e.g., nuclear and organelle) often yield inconsistent results. The tribe Arundinarieae, that is, the temperate bamboos, is a clade of tetraploid originated 22 Ma and subsequently radiated in East Asia. Previous studies of Arundinarieae have found conflicting relationships and/or low support. Here, we obtain nuclear markers from ddRAD data for 213 Arundinarieae taxa and parallel sampling of chloroplast genomes from genome skimming for 147 taxa. We first assess the feasibility of using ddRAD-seq data for phylogenetic estimates of paleopolyploid and rapidly radiated lineages, optimize clustering thresholds, and analysis workflow for orthology identification. Reference-based ddRAD data assembly approaches perform well and yield strongly supported relationships that are generally concordant with morphology-based taxonomy. We recover five major lineages, two of which are notable (the pachymorph and leptomorph lineages), in that they correspond with distinct rhizome morphologies. By contrast, the phylogeny from chloroplast genomes differed significantly. Based on multiple lines of evidence, the ddRAD tree is favored as the best species tree estimation for temperate bamboos. Using a time-calibrated ddRAD tree, we find that Arundinarieae diversified rapidly around the mid-Miocene corresponding with intensification of the East Asian monsoon and the evolution of key innovations including the leptomorph rhizomes. Our results provide a highly resolved phylogeny of Arundinarieae, shed new light on the radiation and reticulate evolutionary history of this tribe, and provide an empirical example for the study of recalcitrant plant radiations. [Arundinarieae; ddRAD; paleopolyploid; genome skimming; rapid diversification; incongruence.]

Extensive reticulate evolution within Fargesia (s.l.) (Bambusoideae: Poaceae) and its allies: Evidence from multiple nuclear markers

Reticulate evolution resulting from hybridization and introgression has been recognized as a creative source of species and diversification in bamboos. Previous phylogenetic studies revealed that Fargesia (s.l.) (Fargesia and Yushania) was divided into the Fargesia spathe clade and the non-spathe clade. Interestingly, the Fargesia spathe clade may have originated from hybridization among other clades within Fargesia (s.l.). Understanding the hybrid origin of this clade requires a robust phylogenetic framework in which major clades within Fargesia (s.l.) are resolved. Here, we used three nuclear genes to reconstruct the evolutionary history of Fargesia (s.l.) and its allies to identify putative patterns in the origin of the Fargesia spathe clade and to examine the extent to which reticulate evolution has occurred at the interspecific level in bamboos. Bashania species form a clade with Fargesia (s.l.), which is further divided into Group I and Group II. The Fargesia spathe clade, the Alpine Bashania clade, and Fargesia yajiangensis comprise Group I, while the Bashania fargesii clade and the remaining Fargesia (s.l.) species form Group II. Incongruence between the current nuclear-based and previous plastid phylogenies demonstrate several possible hybridization events among Fargesia (s.l.) species and related taxa, which have given rise to the Fargesia spathe clade, the Phyllostachys clade, and the Ampelocalamus clade. We also detected several putative hybrid species of Fargesia (s.l.). Our results show that reticulate evolution has played a prominent role in Fargesia (s.l.) evolution, which could, in part, account for the taxonomic difficulty associated with Fargesia (s.l.) and the alpine bamboos. The study also underscores the importance of hybridization in the evolution of bamboos, at both intergeneric and intrageneric levels.

A partial order and cluster-similarity metric on rooted phylogenetic trees

Metrics on rooted phylogenetic trees are integral to a number of areas of phylogenetic analysis. Cluster-similarity metrics have recently been introduced in order to limit skew in the distribution of distances, and to ensure that trees in the neighbourhood of each other have similar hierarchies. In the present paper we introduce a new cluster-similarity metric on rooted phylogenetic tree space that has an associated local operation, allowing for easy calculation of neighbourhoods, a trait that is desirable for MCMC calculations. The metric is defined by the distance on the Hasse diagram induced by a partial order on the set of rooted phylogenetic trees, itself based on the notion of a hierarchy-preserving map between trees. The partial order we introduce is a refinement of the well-known refinement order on hierarchies. Both the partial order and the hierarchy-preserving maps may also be of independent interest.

Two plastomes of Phyllostachys and reconstruction of phylogenic relationship amongst selected Phyllostachys species using genome skimming

The genus Phyllostachys is economically important; however, only a small amount of complete plastid genomes have been reported to date. Here, we characterized two complete chloroplast genomes of Phyllostachys using genome skimming. The chloroplast genomes of Phyllostachys reticulata and Phyllostachys edulis ‘Pachyloen’ were 136,689 bp and 139,678 bp in length, respectively, and their GC contents were 38.8% and 38.9%, respectively. The sequences of each species contained 132 unique genes, including 39 tRNA, eight rRNA, and 85 protein-coding genes. Phylogenetic analysis shows that all selected Phyllostachys species were grouped into one well-supported clade in the Phyllostachys clade (V) of Arundinarieae. Moreover, in terms of chloroplast genome size, structure, and composition, P. edulis ‘Pachyloen’ is identical to P. edulis, further indicating the affinity between them.