ddRAD analyses reveal a credible phylogenetic relationship of the four main genera of Bambusa-Dendrocalamus-Gigantochloa complex (Poaceae: Bambusoideae)

The origin and morphological character evolution of the paleotropical woody bamboos

The woody bamboos (Bambusoideae) exhibit distinctive biological traits within Poaceae, such as highly lignified culms, rapid shoot growth, monocarpic mass flowering and nutlike or fleshy caryopses. Much of the remarkable morphological diversity across the subfamily exists within a single hexaploid clade, the paleotropical woody bamboos (PWB), making it ideal to investigate the factors underlying morphological evolution in woody bamboos. However, the origin and biogeographical history of PWB remain elusive, as does the effect of environmental factors on the evolution of their morphological characters. We generated a robust and time-calibrated phylogeny of PWB using single nucleotide polymorphisms retrieved from optimized double digest restriction site associated DNA sequencing, and explored the evolutionary trends of habit, inflorescence, and caryopsis type in relation to environmental factors including climate, soil, and topography. We inferred that the PWB started to diversify across the Oligocene-Miocene boundary and formed four major clades, that is, Melocanninae, Racemobambosinae s.l. (comprising Dinochloinae, Greslanlinae, Racemobambosinae s.str. and Temburongiinae), Hickeliinae and Bambusinae s.l. (comprising Bambusinae s.str. plus Holttumochloinae). The ancestor of PWB was reconstructed as having erect habit, indeterminate inflorescence and basic caryopsis. The characters including climbing/scrambling habit, determinate inflorescence, and nucoid/bacoid caryopsis have since undergone multiple changes and reversals during the diversification of PWB. The evolution of all three traits was correlated with, and hence likely influenced by, aspects of climate, topography, and soil, with climate factors most strongly correlated with morphological traits, and soil factors least so. However, topography had more influence than climate or soil on the evolution of erect habit, whereas both factors had greater effect on the evolution of bacoid caryopsis than did soil. Our results provide novel insights into morphological diversity and adaptive evolution in bamboos for future ecological and evolutionary research.

Genetic Diversity Analysis and Prediction of Potential Suitable Areas for the Rare and Endangered Wild Plant Henckelia longisepala

Henckelia longisepala (H. W. Li) D. J. Middleton & Mich. Möller is a rare and endangered plant species found only in Southeastern Yunnan, China, and Northern Vietnam. It is listed as a threatened species in China and recognized as a plant species with extremely small populations (PSESP), while also having high ornamental value and utilization potential. This study used ddRAD-seq technology to quantify genetic diversity and structure for 32 samples from three extant populations of H. longisepala. The H. longisepala populations were found to have low levels of genetic diversity (Ho = 0.1216, He = 0.1302, Pi = 0.1731, FIS = 0.1456), with greater genetic differentiation observed among populations (FST = 0.3225). As indicated by genetic structure and phylogenetic analyses, samples clustered into three distinct genetic groups that corresponded to geographically separate populations. MaxEnt modeling was used to identify suitable areas for H. longisepala across three time periods and two climate scenarios (SSP1-2.6, SSP5-8.5). High-suitability areas were identified in Southeastern Yunnan Province, Northern Vietnam, and Eastern Laos. Future H. longisepala distribution was predicted to remain centered in these areas, but with a decrease in the total amount of suitable habitat. The present study provides key data on H. longisepala genetic diversity, as well as a theoretical basis for the conservation, development, and utilization of its germplasm resources.

Deep genome skimming reveals the hybrid origin of Pseudosasa gracilis (Poaceae: Bambusoideae)

Pseudosasa gracilis (Poaceae: Bambusoideae) is a temperate woody bamboo species endemic to South-central China with a narrow distribution. Previous phylogenetic studies revealed an unexpected, isolated phylogenetic position of Ps. gracilis. Here we conducted phylogenomic analysis by sampling populations of Ps. gracilis and its sympatric species Ps. nanunica and Sinosasa polytricha reflecting different genomic signals, by deep genome skimming. Integrating molecular evidence from chloroplast genes and genome-wide SNPs, we deciphered the phylogenetic relationships of Ps. gracilis. Both plastid and nuclear data indicate that Ps. gracilis is more closely related to Sinosasa, which is discordant with the taxonomic treatment. To further explore this molecular-morphological conflict, we screened 411 "perfect-copy" syntenic genes to reconstruct phylogenies using both the concatenation and coalescent methods. We observed extensive discordance between gene trees and the putative species tree. A significant hybridization event was detected based on 411 genes from the D subgenome, showing Ps. gracilis was a hybrid descendant between Sinosasa longiligulata and Ps. nanunica, with 63.56% and 36.44% inheritance probabilities of each parent. Moreover, introgression events were detected in the C subgenome between Ps. gracilis and S. polytricha in the same distribution region. Our findings suggest that sympatric hybridization and introgression play a crucial role in the origin of Ps. gracilis. By providing an empirical example of bamboo of hybrid origin using comprehensive analyses based on genomic data from different inheritance systems and morphological characters, our study represents a step forward in understanding of reticulate evolution of bamboos.

The complete chloroplast genome of Dendrocalamus liboensis Hsueh & D. Z. Li 1985 and its phylogenetic analysis

Dendrocalamus liboensis Hsueh & D. Z. Li 1985 is a unique member of the Bambusoideae subfamily found in Guizhou, China. The species has both economic importance and ornamental value. This study represents the first report of the sequencing and assembly of the complete chloroplast genome of D. liboensis. The total length of the genome was 139,483 bp, with a conventional quadripartite framework consisting of a large single-copy (LSC) region (83,001 bp in length), a small single-copy (SSC) region (12,896 bp in length), and two inverted repeats (IR) regions (both 21,793 bp in length). Overall, the D. liboensis chloroplast genome contained 128 functional genes, including 83 protein-coding genes, 37 tRNAs, and 8 rRNAs. Phylogenetic analysis showed that D. liboensis closely resembled D. sapidus, with both found on a strongly supported branch of the phylogenetic tree.

Phylogenomic analyses reveal reticulate evolution between Neomicrocalamus and Temochloa (Poaceae: Bambusoideae)

Neomicrocalamus and Temochloa are closely related to bamboo genera. However, when considered with newly discovered and morphologically similar material from China and Vietnam, the phylogenetic relationship among these three groups was ambiguous in the analyses based on DNA regions. Here, as a means of investigating the relationships among the three bamboo groups and exploring potential sources of genomic conflicts, we present a phylogenomic examination based on the whole plastome, single-nucleotide polymorphism (SNP), and single-copy nuclear (SCN) gene datasets. Three different phylogenetic hypotheses were found. The inconsistency is attributed to the combination of incomplete lineage sorting and introgression. The origin of newly discovered bamboos is from introgressive hybridization between Temochloa liliana (which contributed 80.7% of the genome) and Neomicrocalamus prainii (19.3%), indicating that the newly discovered bamboos are closer to T. liliana in genetics. The more similar morphology and closer distribution elevation also imply a closer relationship between Temochloa and newly discovered bamboos.

The identity of Dinochloa species and enumeration of Melocalamus (Poaceae: Bambusoideae) in China

Three woody bamboo species collected in Hainan, China in 1940 have been described as Dinochloa based on vegetative specimens. However, the identity of these species has long been in doubt, largely because the vegetative phase in species of Dinochloa is morphologically similar to that in species of Melocalamus, a climbing or scrambling bamboo genus of the paleotropical woody bamboos (Poaceae: Bambusoideae) that consists of about 15 species and one variety. To determine the phylogenetic affinity of the three Dinochloa species from Hainan, we sampled almost all recognized Chinese species of Melocalamus and representative species of Dinochloa as well as other closely related genera, performed molecular phylogenetic analysis, and compared their morphology based on herbarium and fieldwork investigation. Our ddRAD data indicate that the three species from Hainan are closely related to Melocalamus, not Dinochloa. Morphological analysis showed that these three species have a climbing habit but do not grow spirally, their culm leaves have smooth bases, and there is a ring of powder and/or tomenta above and below the nodes. Taken together our findings indicate that the three species from Hainan originally published in Dinochloa should be transferred to Melocalamus, i.e., Melocalamus orenudus (McClure) D.Z. Li & J.X. Liu, Melocalamus puberulus (McClure) D.Z. Li & J.X. Liu, and Melocalamus utilis (McClure) D.Z. Li & J.X. Liu, respectively. This study concludes with an enumeration of Chinese species of Melocalamus, with a key to nine recognized species and one variety, and a lectotypification for M. compatiflorus.

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.

Plastid phylogenomics and morphological character evolution of Chloridoideae (Poaceae)

Chloridoideae is one of the largest subfamilies of Poaceae, containing many species of great economic and ecological value; however, phylogenetic relationships among the subtribes and genera of Cynodonteae are controversial. In the present study, we combined 111 plastomes representing all five tribes, including 25 newly sequenced plastomes that are mostly from Cynodonteae. Phylogenetic analyses supported the five monophyletic tribes of Chloridoideae, including Centropodieae, Triraphideae, Eragrostideae, Zoysieae and Cynodonteae. Simultaneously, nine monophyletic lineages were revealed in Cynodonteae: supersubtribe Boutelouodinae, subtribes Tripogoninae, Aeluropodinae, Eleusininae, Dactylocteniinae, supersubtribe Gouiniodinae, Cleistogenes and Orinus, and subtribe Triodiinae. Within the tribe of Cynodonteae, the basal lineage is supersubtribe Boutelouodinae and Tripogoninae is sister to the remaining lineages. The clade formed of Aeluropodinae and Eleusininae is sister to the clade composed of Dactylocteniinae, supersubtribe Gouiniodinae, Cleistogenes and Orinus, and subtribe Triodiinae. The clade comprising Dactylocteniinae and supersubtribe Gouiniodinae is sister to the clade comprising Cleistogenes, Orinus, and Triodiinae. Acrachne is a genus within Eleusininae but not within Dactylocteniinae. Molecular evidence determined that Diplachne is not clustered with Leptochloa, which indicated that Diplachne should not be combined into Leptochloa. Cleistogenes is sister to a clade composed of Orinus and Triodia, whereas the recently proposed subtribe Orininae was not supported. Cynodonteae was estimated to have experienced rapid divergence within a short period, which could be a major obstacle in resolving its phylogenetic relationships. Ancestral state reconstructions of morphological characters showed that the most recent common ancestor (MRCA) of Chloridoideae has a panicle, multiple florets in each spikelet, the peaked type of stomatal subsidiary cells, and a saddle-shaped phytoliths, while the ancestral morphological characters of Cynodonteae are the panicle, peaked type of stomatal subsidiary cells, sharp-cap cell typed and equal-base-cell microhair, and square-shaped phytoliths. Overall, plastome phylogenomics provides new insights into the phylogenetic relationships and morphological character evolution of Chloridoideae.

Complete Chloroplast Genome Features of Dendrocalamusfarinosus and Its Comparison and Evolutionary Analysis with Other Bambusoideae Species

Dendrocalamus farinosus is one of the essential bamboo species mainly used for food and timber in the southwestern region of China. In this study, the complete chloroplast (cp) genome of D. farinosus is sequenced, assembled, and the phylogenetic relationship analyzed. The cp genome has a circular and quadripartite structure, has a total length of 139,499 bp and contains 132 genes: 89 protein-coding genes, eight rRNAs and 35 tRNAs. The repeat analyses showed that three types of repeats (palindromic, forward and reverse) are present in the genome. A total of 51 simple sequence repeats are identified in the cp genome. The comparative analysis between different species belonging to Dendrocalamus revealed that although the cp genomes are conserved, many differences exist between the genomes. The analysis shows that the non-coding regions were more divergent than the coding regions, and the inverted repeat regions are more conserved than the single-copy regions. Moreover, these results also indicate that rpoC2 may be used to distinguish between different bamboo species. Phylogenetic analysis results supported that D. farinosus was closely related to D. latiflorus. Furthermore, these bamboo species’ geographical distribution and rhizome types indicate two evolutionary pathways: one is from the tropics to the alpine zone, and the other is from the tropics to the warm temperate zone. Our study will be helpful in the determination of the cp genome sequences of D. farinosus, and provides new molecular data to understand the Bambusoideae evolution.

Incompatibility Phylogenetic Signals between Double-Digest Restriction Site-Associated DNA Sequencing and Plastid Genomes in Chinese Curcuma (Zingiberaceae)—A Recent Qinghai–Tibetan Plateau Diversification Genera

Curcuma is of high economic value, credited to its medicinal, edible, and ornamental properties, which possess all signatures of adaptability, and rapid radiation, especially species of Curcuma (Chinese Curcuma, a recent Qinghai–Tibetan Plateau diversification genera) scattered in China. However, little is known about the incongruent phylogenetic signals within this genera from different inheritance patterns that will militate against the further development of this genera. In this research, we applied complete chloroplast genome data together with double-digest restriction site-associated DNA sequencing data (ddRAD-seq) strategy to investigate phylogenetic signals of Chinese Curcuma species, clustering using two RAD analysis pipelines (STACKS and pyRAD). Phylogenetic trees were obtained from each locus based on the maximum likelihood (ML) and multispecies coalescent (BEAST) methods. For visual comparison, multi-method and different datasets were used to infer the phylogeny. We discovered inconsistent relationships for the Chinese Curcuma with varying degrees of support using different methods and datasets.

Integrated Genomic Analyses From Low-Depth Sequencing Help Resolve Phylogenetic Incongruence in the Bamboos (Poaceae: Bambusoideae)

The bamboos (Bambusoideae, Poaceae) comprise a major grass lineage with a complex evolutionary history involving ancient hybridization and allopolyploidy. About 1700 described species are classified into three tribes, Olyreae (herbaceous bamboos), Bambuseae (tropical woody bamboos), and Arundinarieae (temperate woody bamboos). Nuclear analyses strongly support monophyly of the woody tribes, whereas plastome analyses strongly support paraphyly, with Bambuseae sister to Olyreae. Our objectives were to clarify the origin(s) of the woody bamboo tribes and resolve the nuclear vs. plastid conflict using genomic tools. For the first time, plastid and nuclear genomic information from the same bamboo species were combined in a single study. We sampled 51 species of bamboos representing the three tribes, estimated their genome sizes and generated low-depth sample sequence data, from which plastomes were assembled and nuclear repeats were analyzed. The distribution of repeat families was found to agree with nuclear gene phylogenies, but also provides novel insights into nuclear evolutionary history. We infer two early, independent hybridization events, one between an Olyreae ancestor and a woody ancestor giving rise to the two Bambuseae lineages, and another between two woody ancestors giving rise to the Arundinarieae. Retention of the Olyreae plastome associated with differential dominance of nuclear genomes and subsequent diploidization in some lineages explains the paraphyly observed in plastome phylogenetic estimations. We confirm ancient hybridization and allopolyploidy in the origins of the extant woody bamboo lineages and propose biased fractionation and diploidization as important factors in their evolution.

The complete plastid genome of Thyrsostachys siamensis (Poaceae, Bambusoideae)

Thyrsostachys is oligotypic genus of Bambusinae, while its phylogenetic position had been unclear. Here, the complete plastid genome of the type species, T. siamensis, was sequenced and analyzed in this work. The complete genome is a typical quadripartite structure with 139,522 bp in length, comprising of a large single-copy region (LSC, 83,032 bp), a small single-copy region (SSC, 12,892 bp), and a pair of invert repeats regions (IR, 21,799 bp). The genome contains 138 genes, 89 protein-coding genes, 41 tRNA genes, and 8 rRNA genes. The GC content of genome was 38.9%. Phylogenetic analysis indicated T. siamensis was sister to Dendrocalamus birmanicus within Bambusinae.

Phylogenomic resolution of the Ceratitis FARQ complex (Diptera: Tephritidae)

The Ceratitis FARQ complex (formerly FAR complex) includes four frugivorous tephritids, Ceratitis fasciventris, C. anonae, C. rosa and C. quilicii, the latter two causing important agricultural losses in Africa. Although FARQ species can be identified on the basis of subtle morphological differences, they cannot be resolved as monophyletic when trying phylogenetic tree reconstructions based on mitochondrial or nuclear gene fragments except for microsatellites. In this study, we used mitogenome and genome-wide SNPs to investigate the phylogenetic relationship within the complex as well as between all four Ceratitis subgenera. The analysis of 13 species supported the monophyly of the Ceratitis subgenera Ceratitis, Ceratalaspis, Pardalaspis, and recovered Pterandrus as paraphyletic but could not properly resolve species within the FARQ complex. Conversely, gene and species tree reconstructions based on 785,484 genome-wide SNPs could consistently resolve the FARQ taxa and provide insights into their phylogenetic relationships. Gene flow was detected by TreeMix analysis from C. quilicii to C. fasciventris, suggesting the existence of introgression events in the FARQ complex. Our results suggest that genome-wide SNPs represent a suitable tool for the molecular diagnosis of FARQ species and could possibly be used to develop rapid diagnostic methods or to trace the origins of intercepted samples.

Comparative Plastomes and Phylogenetic Analysis of Cleistogenes and Closely Related Genera (Poaceae)

Cleistogenes (Orininae, Cynodonteae, Chloridoideae, Poaceae) is an ecologically important genus. The phylogenetic placement of Cleistogenes and phylogenetic relationships among Cleistogenes taxa remain controversial for a long time. To resolve the intra- and inter-generic relationships of Cleistogenes, the plastomes of 12 Cleistogenes taxa (including 8 species and 4 varieties), one Orinus species, 15 Triodia species, two Tripogon species, and two Aeluropus species were included in the present study. All the taxa showed a similar pattern in plastome structure, gene order, gene content, and IR boundaries. The number of simple sequence repeats ranged from 145 (O. kokonorica) to 161 (T. plurinervata and T. schinzii). Moreover, 1,687 repeats were identified in these taxa, including 1,012 forward, 650 palindromic, 24 reverse, and one complement. Codon usage analysis revealed that these plastomes contained 16,633 (T. stipoides) to 16,678 (T. tomentosa) codons. Sequence divergence analysis among Cleistogenes and closely related genera identified five non-coding regions (trnS-UGA-psbZ, rpl32-trnL-UAG, trnQ-UUG-psbK, trnD-GUC-psbM, trnT-GGU-trnE-UUC). Phylogenetic analysis of complete plastomes indicated that Cleistogenes is sister to a clade composed of Orinus and Triodia, whereas it did not support the sister relationship between the recently proposed subtribe Orininae (Cleistogenes and Orinus) and Triodia. The subtribe Orininae was not supported by our complete plastome data. The split between Cleistogenes and Orinus-Triodia clade go back to 14.01 Ma. Besides, our findings suggested that C. squarrosa and C. songorica are the successive early diverging groups in the phylogenetic analysis. The other 10 taxa are divided into two groups: a monophyletic group composed of Cleistogenes sp. nov. and C. caespitosa var. ramosa is sister to other eight Cleistogenes taxa. Cleistogenes was estimated to have experienced rapid divergence within a short period, which could be a major obstacle in resolving phylogenetic relationships within Cleistogenes. Collectively, our results provided valuable insights into the phylogenetic study of grass species.

Gigantochloa glabrata (Poaceae, Bambusoideae), a new bamboo species from Yunnan, China

Gigantochloa glabrata N. H. Xia & Y. Zeng ex D. Z. Li & Z. C. Xu, sp. nov., a new species of paleotropical woody bamboo has been described and illustrated from Yunnan, China. The new species is morphologically similar to G. albociliata and G. levis, but differs from them by having erect culm sheath blades; culm sheath ligules 4-6 mm high, truncate, denticulate; and with a ring of white tomentum on the intranode and below the node. The new species was mistakenly identified as Gigantochloa albociliata in the Flora of China and was recognised with description of the vegetative characters in 2014, but it was not effectively published. Here, we designate a complete specimen with inflorescence as the type and describe it in accordance with the Code.