Genome Sequences Provide Insights into the Reticulate Origin and Unique Traits of Woody Bamboos

Thinking globally, acting locally in the 21st century: Bamboo to bioproducts and cleaned mine sites

Current solutions to global challenges place tension between global benefits and local impacts. The result is increasing opposition to implementation of beneficial climate policies. Prioritizing investment in projects with tangible local benefits that also contribute to global climate change can resolve this tension and make local communities' partners instead of antagonists to change; the approach advocated is a new take on "thinking globally, acting locally". This approach is a departure from the usual strategy of focusing resources on solutions perceived to have the largest potential global impact, without regards to local concerns. Reclamation of polluted mine sites by using fast growing bamboo to remove heavy metals provides a case study to show what is possible. Effective implementation of thinking globally while acting locally will require increased coordination between different types of researchers, new educational models, and greater stakeholder participation in problem identification and solution development.

Conservation and Divergence of PEPC Gene Family in Different Ploidy Bamboos

Phosphoenolpyruvate carboxylase (PEPC), as a necessary enzyme for higher plants to participate in photosynthesis, plays a key role in photosynthetic carbon metabolism and the stress response. However, the molecular biology of the PEPC family of Bambusoideae has been poorly studied, and the function of its members in the growth and development of Bambusoideae is still unclear. Here, we identified a total of 62 PEPC family members in bamboo. All the PEPC genes in the bamboo subfamily were divided into twelve groups, each group typically containing significantly fewer PEPC members in Olyra latifolia than in Phyllostachys edulis, Dendrocalamus latiflorus and Dendrocalamus brandisii. The results of an intraspecific and interspecies collinearity analysis showed that fragment replication and whole genome replication were the main driving forces of bamboo PEPC members. Furthermore, the Ka/Ks values of collinear genes showed that bamboo PEPC experienced purification selection. In addition, the promoter region of PEPC genes contains cis-acting elements related to light response, plant hormone response and response to stress. An analysis of the expression levels of the PEPC family in different developmental stages and tissues of bamboo shoots has shown that PhePEPC7, PhePEPC9 and PhePEPC10 were highly expressed in the leaves of non-flowering plants and culms. Furthermore, PhePEPC6 was significantly upregulated in leaves after GA treatment. Further research has shown that PhePEPC6 was mainly localized in the cell membrane. This provides a solid bioinformatics foundation for further understanding the biological functions of the bamboo PEPC family.

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.

The CONSTANS-LIKE gene PeCOL13 regulates flowering through intron-retained alternative splicing in Phyllostachys edulis

Woody bamboo exhibits a unique flowering characteristic with a lengthy flowering cycle, often followed by death. In many plant species, alternative splicing (AS) is a common phenomenon involved in controlling flowering. In this study, a PeCOL13 gene in moso bamboo (Phyllostachys edulis) was characterized. It produced two isoforms: PeCOL13α and PeCOL13β, due to an intron-retained AS. The PeCOL13α expressed in the vegetative phase and the reproductive phase, but the PeCOL13β didn't express during the vegetative phase and showed only a weak expression from F1 to F3 during the reproductive phase. Overexpression of PeCOL13α in rice (Oryza sativa) resulted in a delayed heading time through inhibiting the expressions of Hd3a, OsFTL1, and Ehd1 and activating the expressions of Ghd7 and RCN1. However, the PeCOL13β-overexpressed rice didn't show any significant differences in flowering compared with wild-type (WT), and the expressions of downstream flowering genes had no notable changes. Further analysis revealed that both PeCOL13α and PeCOL13β can bind to the PeFT promoter. Meanwhile, PeCOL13α can inhibit the transcription of PeFT, but PeCOL13β showed no effect. When PeCOL13α and PeCOL13β coexist, the inhibitory effect of PeCOL13α on PeFT transcription was weakened by PeCOL13β. This study provides new insights into the mechanism of bamboo flowering research.

Chromosomal-level genome and metabolome analyses of highly heterozygous allohexaploid Dendrocalamus brandisii elucidate shoot quality and developmental characteristics

Dendrocalamus brandisii (Munro) Kurz is a sympodial bamboo species with inimitable taste and flavorful shoots. Its rapid growth and use as high-quality material make this bamboo species highly valued for both food processing and wood applications. However, genome information for D. brandisii is lacking, primarily due to its polyploidy and large genome size. Here, we assembled a high-quality genome for hexaploid D. brandisii, which comprises 70 chromosomes with a total size of 2,756 Mb, using long-read HiFi sequencing. Furthermore, we accurately separated the genome into its three constituent subgenomes. We used Oxford Nanopore Technologies long reads to construct a transcriptomic dataset covering 15 tissues for gene annotation to complement our genome assembly, revealing differential gene expression and post-transcriptional regulation. By integrating metabolome analysis, we unveiled that well-balanced lignin formation, as well as abundant flavonoid and fructose contents, contribute to the superior quality of D. brandisii shoots. Integrating genomic, transcriptomic, and metabolomic datasets provided a solid foundation for enhancing bamboo shoot quality and developing efficient gene-editing techniques. This study should facilitate research on D. brandisii and enhance its use as a food source and wood material by providing crucial genomic resources.

Variations and trade-offs in leaf and culm functional traits among 77 woody bamboo species

BACKGROUND

Woody bamboos are the only diverse large perennial grasses in mesic-wet forests and are widely distributed in the understory and canopy. The functional trait variations and trade-offs in this taxon remain unclear due to woody bamboo syndromes (represented by lignified culm of composed internodes and nodes). Here, we examined the effects of heritable legacy and occurrence site climates on functional trait variations in leaf and culm across 77 woody bamboo species in a common garden. We explored the trade-offs among leaf functional traits, the connection between leaf nitrogen (N), phosphorus (P) concentrations and functional niche traits, and the correlation of functional traits between leaves and culms.

RESULTS

The Bayesian mixed models reveal that the combined effects of heritable legacy (phylogenetic distances and other evolutionary processes) and occurrence site climates accounted for 55.10-90.89% of the total variation among species for each studied trait. The standardized major axis analysis identified trade-offs among leaf functional traits in woody bamboo consistent with the global leaf economics spectrum; however, compared to non-bamboo species, the woody bamboo exhibited lower leaf mass per area but higher N, P concentrations and assimilation, dark respiration rates. The canonical correlation analysis demonstrated a positive correlation (ρ = 0.57, P-value < 0.001) between leaf N, P concentrations and morphophysiology traits. The phylogenetic principal components and trait network analyses indicated that leaf and culm traits were clustered separately, with leaf assimilation and respiration rates associated with culm ground diameter.

CONCLUSION

Our study confirms the applicability of the leaf economics spectrum and the biogeochemical niche in woody bamboo taxa, improves the understanding of woody bamboo leaf and culm functional trait variations and trade-offs, and broadens the taxonomic units considered in plant functional trait studies, which contributes to our comprehensive understanding of terrestrial forest ecosystems.

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.

Phylogenomic profiles of whole-genome duplications in Poaceae and landscape of differential duplicate retention and losses among major Poaceae lineages

Poaceae members shared a whole-genome duplication called rho. However, little is known about the evolutionary pattern of the rho-derived duplicates among Poaceae lineages and implications in adaptive evolution. Here we present phylogenomic/phylotranscriptomic analyses of 363 grasses covering all 12 subfamilies and report nine previously unknown whole-genome duplications. Furthermore, duplications from a single whole-genome duplication were mapped to multiple nodes on the species phylogeny; a whole-genome duplication was likely shared by woody bamboos with possible gene flow from herbaceous bamboos; and recent paralogues of a tetraploid Oryza are implicated in tolerance of seawater submergence. Moreover, rho duplicates showing differential retention among subfamilies include those with functions in environmental adaptations or morphogenesis, including ACOT for aquatic environments (Oryzoideae), CK2β for cold responses (Pooideae), SPIRAL1 for rapid cell elongation (Bambusoideae), and PAI1 for drought/cold responses (Panicoideae). This study presents a Poaceae whole-genome duplication profile with evidence for multiple evolutionary mechanisms that contribute to gene retention and losses.

PySmooth: a Python tool for the removal and correction of genotyping errors

In genetic mapping studies involving many individuals, genome-wide markers such as single nucleotide polymorphisms (SNPs) can be detected using different methods. However, it comes with some errors. Some SNPs associated with diseases can be in regions encoding long noncoding RNAs (lncRNAs). Therefore, identifying the errors in genotype file and correcting them is crucial for accurate genetic mapping studies. We develop a Python tool called PySmooth, that offers an easy-to-use command line interface for the removal and correction of genotyping errors. PySmooth uses the approach of a previous tool called SMOOTH with some modifications. It inputs a genotype file, detects errors and corrects them. PySmooth provides additional features such as imputing missing data, better user-friendly usage, generates summary and visualization files, has flexible parameters, and handles more genotype codes. AVAILABILITY AND IMPLEMENTATION: PySmooth is available at https://github.com/lncRNAAddict/PySmooth .

Genome assemblies of 11 bamboo species highlight diversification induced by dynamic subgenome dominance

Polyploidy (genome duplication) is a pivotal force in evolution. However, the interactions between parental genomes in a polyploid nucleus, frequently involving subgenome dominance, are poorly understood. Here we showcase analyses of a bamboo system (Poaceae: Bambusoideae) comprising a series of lineages from diploid (herbaceous) to tetraploid and hexaploid (woody), with 11 chromosome-level de novo genome assemblies and 476 transcriptome samples. We find that woody bamboo subgenomes exhibit stunning karyotype stability, with parallel subgenome dominance in the two tetraploid clades and a gradual shift of dominance in the hexaploid clade. Allopolyploidization and subgenome dominance have shaped the evolution of tree-like lignified culms, rapid growth and synchronous flowering characteristic of woody bamboos as large grasses. Our work provides insights into genome dominance in a remarkable polyploid system, including its dependence on genomic context and its ability to switch which subgenomes are dominant over evolutionary time.

Genome architecture and genetic diversity of allopolyploid okra (Abelmoschus esculentus)

The allopolyploid okra (Abelmoschus esculentus) unveiled telomeric repeats flanking distal gene-rich regions and short interstitial TTTAGGG telomeric repeats, possibly representing hallmarks of chromosomal speciation. Ribosomal RNA (rRNA) genes organize into 5S clusters, distinct from the 18S-5.8S-28S units, indicating an S-type rRNA gene arrangement. The assembly, in line with cytogenetic and cytometry observations, identifies 65 chromosomes and a 1.45 Gb genome size estimate in a haploid sibling. The lack of aberrant meiotic configurations implies limited to no recombination among sub-genomes. k-mer distribution analysis reveals 75% has a diploid nature and 15% heterozygosity. The configurations of Benchmarking Universal Single-Copy Ortholog (BUSCO), k-mer, and repeat clustering point to the presence of at least two sub-genomes one with 30 and the other with 35 chromosomes, indicating the allopolyploid nature of the okra genome. Over 130 000 putative genes, derived from mapped IsoSeq data and transcriptome data from public okra accessions, exhibit a low genetic diversity of one single nucleotide polymorphisms per 2.1 kbp. The genes are predominantly located at the distal chromosome ends, declining toward central scaffold domains. Long terminal repeat retrotransposons prevail in central domains, consistent with the observed pericentromeric heterochromatin and distal euchromatin. Disparities in paralogous gene counts suggest potential sub-genome differentiation implying possible sub-genome dominance. Amino acid query sequences of putative genes facilitated phenol biosynthesis pathway annotation. Comparison with manually curated reference KEGG pathways from related Malvaceae species reveals the genetic basis for putative enzyme coding genes that likely enable metabolic reactions involved in the biosynthesis of dietary and therapeutic compounds in okra.

Widespread incomplete lineage sorting and introgression shaped adaptive radiation in the Gossypium genus

Cotton (Gossypium) stands as a crucial economic crop, serving as the primary source of natural fiber for the textile sector. However, the evolutionary mechanisms driving speciation within the Gossypium genus remain unresolved. In this investigation, we leveraged 25 Gossypium genomes and introduced four novel assemblies-G. harknessii, G. gossypioides, G. trilobum, and G. klotzschianum (Gklo)-to delve into the speciation history of this genus. Notably, we encountered intricate phylogenies potentially stemming from introgression. These complexities are further compounded by incomplete lineage sorting (ILS), a factor likely to have been instrumental in shaping the swift diversification of cotton. Our focus subsequently shifted to the rapid radiation episode during a concise period in Gossypium evolution. For a recently diverged lineage comprising G. davidsonii, Gklo, and G. raimondii, we constructed a finely detailed ILS map. Intriguingly, this analysis revealed the non-random distribution of ILS regions across the reference Gklo genome. Moreover, we identified signs of robust natural selection influencing specific ILS regions. Noteworthy variations pertaining to speciation emerged between the closely related sister species Gklo and G. davidsonii. Approximately 15.74% of speciation structural variation genes and 12.04% of speciation-associated genes were estimated to intersect with ILS signatures. These findings enrich our understanding of the role of ILS in adaptive radiation, shedding fresh light on the intricate speciation history of the Gossypium genus.

Mysterious Bamboo flowering phenomenon: A literature review and new perspectives

Bamboo, a globally distributed non-timber forest resource, plays a critical role in local ecosystems and economies. Despite its significance, the understanding of bamboo's long and unpredictable flowering cycles remains limited. Our bibliometric analysis of bamboo flowering-related literature from the Web of Science database reveals an initial focus on regeneration studies, with a recent trend shifting towards microscopic and molecular perspectives. Furthermore, our narrative review emphasizes the importance of considering factors such as the proportion of flowering culms and the duration of flowering in classifying bamboo flowering phenomena. While numerous studies have endorsed the predator saturation hypothesis as a suitable explanation for the synchronicity of bamboo flowering, no existing theory explains bamboo's prolonged flowering cycles. We propose a new natural selection hypothesis as a potential explanation for these extraordinary cycles, underscoring the need for further research in this area. Despite the substantial volume of data accumulated on bamboo flowering, these resources have not been fully exploited in recent research. Future studies would benefit from more comprehensive data collection methods, encompassing field observations, satellite remote sensing data, and omics data. The convergence of traditional ecological studies with molecular techniques may pave the way for significant advancements in bamboo flowering research.

Haplotype-based phylogenetic analysis and population genomics uncover the origin and domestication of sweetpotato

The hexaploid sweetpotato (Ipomoea batatas) is one of the most important root crops worldwide. However, its genetic origin remains controversial, and its domestication history remains unknown. In this study, we used a range of genetic evidence and a newly developed haplotype-based phylogenetic analysis to identify two probable progenitors of sweetpotato. The diploid progenitor was likely closely related to Ipomoea aequatoriensis and contributed the B1 subgenome, IbT-DNA2, and the lineage 1 type of chloroplast genome to sweetpotato. The tetraploid progenitor of sweetpotato was most likely I. batatas 4x, which donated the B2 subgenome, IbT-DNA1, and the lineage 2 type of chloroplast genome. Sweetpotato most likely originated from reciprocal crosses between the diploid and tetraploid progenitors, followed by a subsequent whole-genome duplication. In addition, we detected biased gene exchanges between the subgenomes; the rate of B1 to B2 subgenome conversions was nearly three times higher than that of B2 to B1 subgenome conversions. Our analyses revealed that genes involved in storage root formation, maintenance of genome stability, biotic resistance, sugar transport, and potassium uptake were selected during the speciation and domestication of sweetpotato. This study sheds light on the evolution of sweetpotato and paves the way for improvement of this crop.

Characteristics of Amorphophallus konjac as indicated by its genome

Amorphophallus konjac, belonging to the genus Amorphophallus of the Araceae family, is an economically important crop widely used in health products and biomaterials. In the present work, we performed the whole-genome assembly of A. konjac based on the NovaSeq platform sequence data. The final genome assembly was 4.58 Gb with a scaffold N50 of 3212 bp. The genome includes 39,421 protein-coding genes, and 71.75% of the assemblies were repetitive sequences. Comparative genomic analysis showed 1647 gene families have expanded and 2685 contracted in the A. konjac genome. Likewise, genome evolution analysis indicated that A. konjac underwent whole-genome duplication, possibly contributing to the expansion of certain gene families. Furthermore, we identified many candidate genes involved in the tuber formation and development, cellulose and lignification synthesis. The genome of A. konjac obtained in this work provides a valuable resource for the further study of the genetics, genomics, and breeding of this economically important crop, as well as for evolutionary studies of Araceae family.

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.

A comparative full-length transcriptomic resource provides insight into the perennial monocarpic mass flowering

Perennial monocarpic mass flowering represents as a key developmental innovation in flowering time diversity in several biological and economical essential families, such as the woody bamboos and the shrubby Strobilanthes. However, molecular and genetic mechanisms underlying this important biodiversity remain poorly investigated. Here, we generated a full-length transcriptome resource incorporated into the BlueOmics database (http://blueomics.iflora.cn) for two Strobilanthes species, which feature contrasting flowering time behaviors. Using about 112 and 104 Gb Iso-seq reads together with ~185 and ~75 Gb strand-specific RNA seq data, we annotated 80 971 and 79 985 non-redundant full-length transcripts for the perennial polycarpic Strobilanthes tetrasperma and the perennial monocarpic Strobilanthes biocullata, respectively. In S. tetrasperma, we identified 8794 transcripts showing spatiotemporal expression in nine tissues. In leaves and shoot apical meristems at two developmental stages, 977 and 1121 transcripts were differentially accumulated in S. tetrasperma and S. biocullata, respectively. Interestingly, among the 33 transcription factors showing differential expression in S. tetrasperma but without differential expression in S. biocullata, three were involved potentially in the photoperiod and circadian-clock pathway of flowering time regulation (FAR1 RELATED SEQUENCE 12, FRS12; NUCLEAR FACTOR Y A1, NFYA1; PSEUDO-RESPONSE REGULATOR 5, PRR5), hence provides an important clue in deciphering the flowering diversity mechanisms. Our data serve as a key resource for further dissection of molecular and genetic mechanisms underpinning key biological innovations, here, the perennial monocarpic mass flowering.

Phylotranscriptomic analyses reveal deep gene tree discordance in Camellia (Theaceae)

Gene tree discordance is a significant legacy of biological evolution. Multiple factors can result in incongruence among genes, such as introgression, incomplete lineage sorting (ILS), gene duplication or loss. Resolving the background of gene tree discordance is a critical way to uncover the process of species diversification. Camellia, the largest genus in Theaceae, has controversial taxonomy and systematics due in part to a complex evolutionary history. We used 60 transcriptomes of 55 species, which represented 15 sections of Camellia to investigate its phylogeny and the possible causes of gene tree discordance. We conducted gene tree discordance analysis based on 1,617 orthologous low-copy nuclear genes, primarily using coalescent species trees and polytomy tests to distinguish hard and soft conflict. A selective pressure analysis was also performed to assess the impact of selection on phylogenetic topology reconstruction. Our results detected different levels of gene tree discordance in the backbone of Camellia, and recovered rapid diversification as one of the possible causes of gene tree discordance. Furthermore, we confirmed that none of the currently proposed sections of Camellia was monophyletic. Comparisons among datasets partitioned under different selective pressure regimes showed that integrating all orthologous genes provided the best phylogenetic resolution of the species tree of Camellia. The findings of this study reveal rapid diversification as a major source of gene tree discordance in Camellia and will facilitate future investigation of reticulate relationships at the species level in this important plant genus.

Phylogenomic Analyses Reveal an Allopolyploid Origin of Core Didymocarpinae (Gesneriaceae) Followed by Rapid Radiation

Allopolyploid plants have long been regarded as possessing genetic advantages under certain circumstances due to the combined effects of their hybrid origins and duplicated genomes. However, the evolutionary consequences of allopolyploidy in lineage diversification remain to be fully understood. Here, we investigate the evolutionary consequences of allopolyploidy using 138 transcriptomic sequences of Gesneriaceae, including 124 newly sequenced, focusing particularly on the largest subtribe Didymocarpinae. We estimated the phylogeny of Gesneriaceae using concatenated and coalescent-based methods based on five different nuclear matrices and 27 plastid genes, focusing on relationships among major clades. To better understand the evolutionary affinities in this family, we applied a range of approaches to characterize the extent and cause of phylogenetic incongruence. We found that extensive conflicts between nuclear and chloroplast genomes and among nuclear genes were caused by both incomplete lineage sorting (ILS) and reticulation, and we found evidence of widespread ancient hybridization and introgression. Using the most highly supported phylogenomic framework, we revealed multiple bursts of gene duplication throughout the evolutionary history of Gesneriaceae. By incorporating molecular dating and analyses of diversification dynamics, our study shows that an ancient allopolyploidization event occurred around the Oligocene-Miocene boundary, which may have driven the rapid radiation of core Didymocarpinae.

Identification and characterization of the cupin_1 domain-containing proteins in ma bamboo (Dendrocalamus latiflorus) and their potential role in rhizome sprouting

Cupin_1 domain-containing protein (CDP) family, which is a member of the cupin superfamily with the most diverse functions in plants, has been found to be involved in hormone pathways that are closely related to rhizome sprouting (RS), a vital form of asexual reproduction in plants. Ma bamboo is a typical clumping bamboo, which mainly reproduces by RS. In this study, we identified and characterized 53 Dendrocalamus latiflorus CDP genes and divided them into seven subfamilies. Comparing the genetic structures among subfamilies showed a relatively conserved gene structure within each subfamily, and the number of cupin_1 domains affected the conservation among D. latiflorus CDP genes. Gene collinearity results showed that segmental duplication and tandem duplication both contributed to the expansion of D. latiflorus CDP genes, and lineage-specific gene duplication was an important factor influencing the evolution of CDP genes. Expression patterns showed that CDP genes generally had higher expression levels in germinating underground buds, indicating that they might play important roles in promoting shoot sprouting. Transcription factor binding site prediction and co-expression network analysis indicated that D. latiflorus CDPs were regulated by a large number of transcription factors, and collectively participated in rhizome buds and shoot development. This study significantly provided new insights into the evolutionary patterns and molecular functions of CDP genes, and laid a foundation for further studying the regulatory mechanisms of plant rhizome sprouting.

Advancements and challenges in bamboo breeding for sustainable development

Bamboo is a highly renewable biomass resource with outstanding ecological, economic and social benefits. However, its lengthy vegetative growth stage and uncertain flowering period have hindered the application of traditional breeding methods. In recent years, significant progress has been made in bamboo breeding. While technical advances in bamboo breeding have been impressive, it is essential to also consider the broader implications we can learn from bamboo's extraordinary features for sustainable development. This review provides an overview of the current status of bamboo breeding technology, including a detailed history of bamboo breeding divided into four eras, a comprehensive map of bamboo germplasm gardens worldwide, with a focus on China, and a summary of available transgenic technologies for gene function verification and genetic improvement. As the demand for bamboo as a sustainable and renewable resource increases continuously, breeding objectives should be focused on enhancing yield, wood properties and adaptability to diverse environmental conditions. In particular, priority should be given to improving fiber length, internode length and wall thickness, as well as regulating lignin and cellulose content for papermaking, substitute for plastic and other applications. Furthermore, we highlight the challenges and opportunities for future research and development in bamboo breeding, including the application of omics technologies, artificial intelligence and the development of new breeding methods. Finally, by integrating the technical advances in bamboo breeding with a discussion of its broader implications for sustainable development, this review provides a comprehensive framework for the development of bamboo industry.

Deciphering complex reticulate evolution of Asian Buddleja (Scrophulariaceae): insights into the taxonomy and speciation of polyploid taxa in the Sino-Himalayan region

BACKGROUND AND AIMS

Species of the genus Buddleja in Asia are mainly distributed in the Sino-Himalayan region and form a challenging taxonomic group, with extensive hybridization and polyploidization. A phylogenetic approach to unravelling the history of reticulation in this lineage will deepen our understanding of the speciation in biodiversity hotspots.

METHODS

For this study, we obtained 80 accessions representing all the species in the Asian Buddleja clade, and the ploidy level of each taxon was determined by flow cytometry analyses. Whole plastid genomes, nuclear ribosomal DNA, single nucleotide polymorphisms and a large number of low-copy nuclear genes assembled from genome skimming data were used to investigate the reticulate evolutionary history of Asian Buddleja. Complex cytonuclear conflicts were detected through a comparison of plastid and species trees. Gene tree incongruence was also analysed to detect any reticulate events in the history of this lineage.

KEY RESULTS

Six hybridization events were detected, which are able to explain the cytonuclear conflict in Asian Buddleja. Furthermore, PhyloNet analysis combining species ploidy data indicated several allopolyploid speciation events. A strongly supported species tree inferred from a large number of low-copy nuclear genes not only corrected some earlier misinterpretations, but also indicated that there are many Asian Buddleja species that have been lumped mistakenly. Divergent time estimation shows two periods of rapid diversification (8-10 and 0-3 Mya) in the Asian Buddleja clade, which might coincide with the final uplift of the Hengduan Mountains and Quaternary climate fluctuations, respectively.

CONCLUSIONS

This study presents a well-supported phylogenetic backbone for the Asian Buddleja species, elucidates their complex and reticulate evolutionary history and suggests that tectonic activity, climate fluctuations, polyploidization and hybridization together promoted the diversification of this lineage.

Molecular identification of Bambusa changningensis is the natural bamboo hybrid of B. rigida × Dendrocalamus farinosus

Bamboo is one of the fastest-growing plants commonly used in food, fibre, paper, biofuel, ornamental and medicinal industries. Natural hybridization in bamboo is rare due to its long vegetative period followed by gregarious flowering and death of the entire population. In the current study, a new bamboo species, Bambusa changningensis, shows intermediate characteristics of Dendrocalamus farinosus and B. rigida morphologically, but it is unknown whether B. changningensis is a natural hybrid. Moreover, B. changningensis has been identified as a superior variety of Sichuan Province with high pulping yield, fibre length and width. Therefore, we analyzed the morphological characteristics, DNA markers, DNA barcoding and chloroplast genomes to identify the hybrid origin of B. changningensis and possible maternal parent. We have developed the transcriptomic data for B. changningensis and mined the SSR loci. The putative parental lines and hybrid were screened for 64 SSR makers and identified that SSR14, SSR28, SSR31 and SSR34 markers showed both alleles of the parental species in B. changningensis, proving heterozygosity. Sequencing nuclear gene GBSSI partial regions and phylogenetic analysis also confirm the hybrid nature of B. changningensis. Further, we have generated the complete chloroplast genome sequence (139505 bp) of B. changningensis. By analyzing the cp genomes of both parents and B. changningensis, we identified that B. rigida might be the female parent. In conclusion, our study identified that B. changningensis is a natural hybrid, providing evidence for bamboo's natural hybridization. This is the first report on confirming a natural bamboo hybrid and its parents through SSR and chloroplast genome sequence.

Whole-genome Duplication Reshaped Adaptive Evolution in A Relict Plant Species, Cyclocarya paliurus

Cyclocarya paliurus is a relict plant species that survived the last glacial period and shows a population expansion recently. Its leaves have been traditionally used to treat obesity and diabetes with the well-known active ingredient cyclocaric acid B. Here, we presented three C. paliurus genomes from two diploids with different flower morphs and one haplotype-resolved tetraploid assembly. Comparative genomic analysis revealed two rounds of recent whole-genome duplication events and identified 691 genes with dosage effects that likely contribute to adaptive evolution through enhanced photosynthesis and increased accumulation of triterpenoids. Resequencing analysis of 45 C. paliurus individuals uncovered two bottlenecks, consistent with the known events of environmental changes, and many selectively swept genes involved in critical biological functions, including plant defense and secondary metabolite biosynthesis. We also proposed the biosynthesis pathway of cyclocaric acid B based on multi-omics data and identified key genes, in particular gibberellin-related genes, associated with the heterodichogamy in C. paliurus species. Our study sheds light on evolutionary history of C. paliurus and provides genomic resources to study the medicinal herbs.

Ecological niche modelling and population genetic analysis of Indian temperate bamboo Drepanostachyum falcatum in the western Himalayas

The present study was conducted to understand the key ecological and biological questions of conservation importance in Drepanostachyum falcatum which aimed to map potential distribution in the western Himalayas and decipher spatial genetic structure. Eco-distribution maps were generated through ecological niche modelling using the Maximum Entropy (MaxEnt) algorithm implemented with 228 geocoordinates of species presence and 12 bioclimatic variables. Concomitantly, 26 natural populations in the western Himalayas were genetically analysed using ten genomic sequence-tagged microsatellite (STMS) markers. Model-derived distribution was adequately supported with appropriate statistical measures, such as area under the 'receiver operating characteristics (ROC)' curve (AUC; 0.917 ± 0.034)", Kappa (K; 0.418), normalized mutual information (NMI; 0.673) and true skill statistic (TSS; 0.715). Further, Jackknife test and response curves showed that the precipitation (pre- and post-monsoon) and temperature (average throughout the year and pre-monsoon) maximize the probabilistic distribution of D. falcatum. We recorded a wide and abundant (4096.86 km²) distribution of D. falcatum in the western Himalayas with maximum occurrence at 1500 to 2500 m asl. Furthermore, marker analysis exemplified high gene diversity with low genetic differentiation in D. falcatum. Relatively, the populations of Uttarakhand are more genetically diverse than Himachal Pradesh, whereas within the Uttarakhand, the Garhwal region captured a higher allelic diversity than Kumaon. Clustering and structure analysis indicated two major gene pools, where genetic admixing appeared to be controlled by long-distance gene flow, horizontal geographical distance, aspect, and precipitation. Both the species distribution map and population genetic structure derived herein may serve as valuable resources for conservation and management of Himalayan hill bamboos.

Deep reticulation: the long legacy of hybridization in vascular plant evolution

Hybridization has long been recognized as a fundamental evolutionary process in plants but, until recently, our understanding of its phylogenetic distribution and biological significance across deep evolutionary scales has been largely obscure. Over the past decade, genomic and phylogenomic datasets have revealed, perhaps not surprisingly, that hybridization, often associated with polyploidy, has been common throughout the evolutionary history of plants, particularly in various lineages of flowering plants. However, phylogenomic studies have also highlighted the challenges of disentangling signals of ancient hybridization from other sources of genomic conflict (in particular, incomplete lineage sorting). Here, we provide a critical review of ancient hybridization in vascular plants, outlining well-documented cases of ancient hybridization across plant phylogeny, as well as the challenges unique to documenting ancient versus recent hybridization. We provide a definition for ancient hybridization, which, to our knowledge, has not been explicitly attempted before. Further documenting the extent of deep reticulation in plants should remain an important research focus, especially because published examples likely represent the tip of the iceberg in terms of the total extent of ancient hybridization. However, future research should increasingly explore the macroevolutionary significance of this process, in terms of its impact on evolutionary trajectories (e.g. how does hybridization influence trait evolution or the generation of biodiversity over long time scales?), as well as how life history and ecological factors shape, or have shaped, the frequency of hybridization across geologic time and plant phylogeny. Finally, we consider the implications of ubiquitous ancient hybridization for how we conceptualize, analyze, and classify plant phylogeny. Networks, as opposed to bifurcating trees, represent more accurate representations of evolutionary history in many cases, although our ability to infer, visualize, and use networks for comparative analyses is highly limited. Developing improved methods for the generation, visualization, and use of networks represents a critical future direction for plant biology. Current classification systems also do not generally allow for the recognition of reticulate lineages, and our classifications themselves are largely based on evidence from the chloroplast genome. Updating plant classification to better reflect nuclear phylogenies, as well as considering whether and how to recognize hybridization in classification systems, will represent an important challenge for the plant systematics community.

Morphological and Anatomical Analysis of the Internodes of a New Dwarf Variant of Moso Bamboo, Phyllostachys edulis f. exaurita

The lack of mutants due to the long periods between flowering of bamboo plants is one of the limiting factors inhibiting research progress in the culm development of bamboo plants. In this study, a stable new dwarf variant of Phyllostachys edulis (Moso bamboo), Phyllostachys edulis f. exaurita T. G. Chen, was discovered and was characterized morphologically, anatomically, and physiologically. The height, diameter at breast height, number of internodes, length and wall thickness of internodes, length, width and number of parenchyma cells of internodes, and morphology of the wide-type (WT) and dwarf variant vascular bundles were compared. The height of the variant was only 49% that of the WT Moso bamboo. It was concluded that the decrease in internode number and length was the cause of dwarfism in P. edulis f. exaurita. The decreased internode length was the result of a decrease in cell number and cell length in the internode. In addition, the laws of change of internode length, internode thickness, cell length, and cell number differed between the WT Moso bamboo and the variant. Furthermore, lower IAA and zeatin concentrations were detected in the buds of the variant. These results suggest that P. edulis f. exaurita is a variant with inhibited primary thickening growth, which is valuable for interpretating the molecular mechanisms underlying the primary thickening growth of bamboo that are still largely unknown.

Phylogeny, transcriptional profile, and auxin-induced phosphorylation modification characteristics of conserved PIN proteins in Moso bamboo (Phyllostachys edulis)

Auxin polar transport is an important way for auxin to exercise its function, and auxin plays an irreplaceable role in the rapid growth of Moso bamboo. We identified and performed the structural analysis of PIN-FORMED auxin efflux carriers in Moso bamboo and obtained a total of 23 PhePIN genes from five gene subfamilies. We also performed chromosome localization and intra- and inter-species synthesis analysis. Phylogenetic analyses of 216 PIN genes showed that PIN genes are relatively conserved in the evolution of the Bambusoideae and have undergone intra-family segment replication in Moso bamboo. The PIN genes' transcriptional patterns showed that the PIN1 subfamily plays a major regulatory role. PIN genes and auxin biosynthesis maintain a high degree of consistency in spatial and temporal distribution. Phosphoproteomics analysis identified many phosphorylated protein kinases that respond to auxin regulation through autophosphorylation and phosphorylation of PIN proteins. The protein interaction network showed that there is a plant hormone interaction regulatory network with PIN protein as the core. We provide a comprehensive PIN protein analysis that complements the auxin regulatory pathway in Moso bamboo and paves the way for further auxin regulatory studies in bamboo.

Genome survey sequencing-based SSR marker development and their validation in Dendrocalamus longispathus

Bamboo is an important genetic resource in India, supporting rural livelihood and industries. Unfortunately, most Indian bamboo taxa are devoid of basic genomic or marker information required to comprehend the genetic processes for further conservation and management. In this study, we perform genome survey sequencing for development of de novo genomic SSRs in Dendrocalamus longispathus, a socioeconomically important bamboo species of northeast India. Using Illumina platform, 69.49 million raw reads were generated and assembled into 1,145,321 contig with GC content 43% and N50 1228 bp. In total, 46,984 microsatellite repeats were mined-out wherein di-nucleotide repeats were most abundant (54.71%) followed by mono- (31.91%) and tri-repeats (9.85%). Overall, AT-rich repeats were predominant in the genome, but GC-rich motifs were more frequent in tri-repeats. Afterwards, 21,596 SSR loci were successfully tagged with the primer pairs, and a subset of 50 were validated through polymerase chain reaction amplification. Of these, 36 SSR loci were successfully amplified, and 16 demonstrated polymorphism. Using 13 polymorphic SSRs, a moderate level of gene diversity (He = 0.480; Ar = 3.52) was recorded in the analysed populations of D. longispathus. Despite the high gene flow (Nm = 4.928) and low genetic differentiation (F_ST = 0.119), severe inbreeding (F_IS = 0.407) was detected. Further, genetic clustering and STRUCTURE analysis revealed that the entire genetic variability is captured under two major gene pools. Conclusively, we present a comprehensive set of novel SSR markers in D. longispathus as well as other taxa of tropical woody bamboos.

New insights into intergeneric relationships of Hickeliinae (Poaceae: Bambusoideae) revealed by complete plastid genomes

The Hickeliinae (Poaceae: Bambusoideae) is an ecologically and economically significant subtribe of tropical bamboos restricted to Madagascar, Comoros, Reunion Island, and a small part of continental Africa (Tanzania). Because these bamboos rarely flower, field identification is challenging, and inferring the evolutionary history of Hickeliinae from herbarium specimens is even more so. Molecular phylogenetic work is critical to understanding this group of bamboos. Here, comparative analysis of 22 newly sequenced plastid genomes showed that members of all genera of Hickeliinae share evolutionarily conserved plastome structures. We also determined that Hickeliinae plastome sequences are informative for phylogenetic reconstructions. Phylogenetic analysis showed that all genera of Hickeliinae are monophyletic, except for Nastus, which is paraphyletic and forms two distant clades. The type species of Nastus (Clade II) is endemic to Reunion Island and is not closely related to other sampled species of Nastus endemic to Madagascar (Clade VI). Clade VI (Malagasy Nastus) is sister to the Sokinochloa + Hitchcockella clade (Clade V), and both clades have a clumping habit with short-necked pachymorph rhizomes. The monotypic Decaryochloa is remarkable in having the longest floret in Bambuseae and forms a distinct Clade IV. Clade III, which has the highest generic diversity, consists of Cathariostachys, Perrierbambus, Sirochloa, and Valiha, which are also morphologically diverse. This work provides significant resources for further genetic and phylogenomic studies of Hickeliinae, an understudied subtribe of bamboo.

Comparative proteomic analysis and antioxidant enzyme activity provide new insights into the embryogenic competence of Guadua chacoensis (Bambusoideae, Poaceae)

Somatic embryogenesis (SE) involves modifications of cellular, biochemical, genetic, and epigenetic patterns. Our work investigated proteins as markers of embryogenic response and characterized the redox state of embryogenic cultures (EC) of Guadua chacoensis. We identified a total of 855 proteins; 129 were up- and 136 down-accumulated in EC as compared with non-embryogenic culture (NEC). Additionally, 37 and 22 proteins were identified as unique in EC and NEC, respectively. Heat-shock proteins as unique proteins and increased activity in Superoxide Dismutase and Guaiacol Peroxidase in EC suggest that the embryogenic response requires activation of the stress response mechanism. Ribosomal, translational, and glycolytic proteins in EC seem to be associated with protein synthesis and energy sources for embryo development, respectively. Accumulation of cell wall-related proteins, such as Arabinogalactan and Polygalacturonase inhibitors, and signaling transduction proteins, including Chitinase, Phospholipase, and Guanine nucleotide-binding proteins in EC seems to be associated with embryogenic response. Enhancement of H₂O₂ content in EC compared to NEC suggests a possible role as a secondary messenger in SE. Altogether, the present study identified marker proteins of embryogenic response in G. chacoensis and revealed the activation of ROS scavenging enzymes to assure cell redox homeostasis and SE responses. SIGNIFICANCE: Somatic embryogenesis is a promising technique for the propagation and conservation of bamboo species; however, this route has been the least understood and studied until now. This study corresponds to the first work approaching proteomics complemented with biochemical analyses in the somatic embryogenesis of bamboo, bringing robust and precise information that can improve our understanding of this complex morphogenetic route.

Development of an efficient expression system with large cargo capacity for interrogation of gene function in bamboo based on bamboo mosaic virus

Bamboo is one of the fastest growing plants among monocotyledonous species and is grown extensively in subtropical regions. Although bamboo has high economic value and produces much biomass quickly, gene functional research is hindered by the low efficiency of genetic transformation in this species. We therefore explored the potential of a bamboo mosaic virus (BaMV)-mediated expression system to investigate genotype-phenotype associations. We determined that the sites between the triple gene block proteins (TGBps) and the coat protein (CP) of BaMV are the most efficient insertion sites for the expression of exogenous genes in both monopodial and sympodial bamboo species. Moreover, we validated this system by individually overexpressing the two endogenous genes ACE1 and DEC1, which resulted in the promotion and suppression of internode elongation, respectively. In particular, this system was able to drive the expression of three 2A-linked betalain biosynthesis genes (more than 4 kb in length) to produce betalain, indicating that it has high cargo capacity and may provide the prerequisite basis for the development of a DNA-free bamboo genome editing platform in the future. Since BaMV can infect multiple bamboo species, we anticipate that the system described in this study will greatly contribute to gene function research and further promote the molecular breeding of bamboo.

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.

Population diversity analyses provide insights into key horticultural traits of Chinese native thymes

Chinese native thymes (CNTs) in the genus Thymus (family Lamiaceae) are rich in bioactive terpenes, which exert antiviral, anti-inflammatory, antioxidation, immunological, and antimicrobial effects. Plants exhibit morphological variation, including erect-type and creeping-type growth forms; however, the molecular mechanisms underlying important horticultural traits have not been determined. Here, we collected 39 CNTs providing strategic plant resources for studies of lignin, terpenoids, and glandular trichomes of thymes. Using resequencing data as well as phenotypic, metabonomic, phylogenetic, population genetic, and transcriptomic analyses, we identified and characterized key genes involved in lignin biosynthesis, terpenoid biosynthesis, and glandular trichome formation. We found many regulatory genes or transcription factors related to these three important horticultural traits, including genes encoding caffeic acid O-methyltransferase (COMT), terpene synthase (TPS), v-myb avian myeloblastosis viral oncogene homolog (MYB), and homeodomain-leucine zipper (HD-ZIP). Population diversity analyses provided insights into growth form, terpenoid, and glandular trichome evolution in CNTs. Furthermore, our results revealed that T. mongolicus accessions might be wild ancestors, and T. quinquecostatus, T. quinquecostatus var. asiaticus, and T. quinquecostatus var. przewalskii might be transitional accessions that derived from T. mongolicus accessions. Finally, T. nervulosus, T. inaequalis, T. mandschuricus, T. curtus, T. amurensis, T. proximus, T. altaicus, T. roseus, and T. marschallianus showed high divergence. We found evidence for introgression between erect-type European cultivated thymes and CNTs. These findings improve our understanding of the determinants of variation in horticultural traits and provide candidate loci for research and breeding.

Overexpression of <italic>PvSVP1</italic>, an <italic>SVP</italic>-like gene of bamboo, causes early flowering and abnormal floral organs in <italic>Arabidopsis</italic> and rice

<p indent="0mm">Bamboo is a nontimber woody plant featuring a long vegetative stage and uncertain flowering time. Therefore, the genes belonging to flowering repressors might be essential in regulating the transition from the vegetative to reproductive stage in bamboo. The <italic>Short Vegetative Phase</italic> ( <italic>SVP</italic>) gene plays a pivotal role in floral transition and development. However, little is known about the bamboo <italic>SVP</italic> homologues. In this study, <italic>Phyllostachys violascens</italic> <italic>PvSVP1</italic> is isolated by analysis of the <italic>P</italic>. <italic>edulis</italic> transcriptome database. Phylogenetic analysis shows that <italic>PvSVP1</italic> is closely related to <italic>OsMADS55</italic> (rice <italic>SVP</italic> homolog). <italic>PvSVP1</italic> is ubiquitously expressed in various tissues, predominantly in vegetative tissues. To investigate the function of <italic>PvSVP1</italic>, <italic>PvSVP1</italic> is overexpressed in <italic>Arabidopsis</italic> and rice under the influence of the 35S promoter. Overexpression of <italic>PvSVP1</italic> in <italic>Arabidopsis</italic> causes early flowering and produces abnormal petals and sepals. Quantitative real-time PCR reveals that overexpression in <italic>Arabidopsis</italic> produces an early flowering phenotype by downregulating <italic>FLC</italic> and upregulating <italic>FT</italic> and produces abnormal floral organs by upregulating <italic>AP1</italic>, <italic>AP3</italic> and <italic>PI</italic> expressions. Simultaneously, overexpression of <italic>PvSVP1</italic> in rice alters the expressions of flowering-related genes such as <italic>Hd3a</italic>, <italic>RFT1</italic>, <italic>OsMADS56</italic> and <italic>Ghd7</italic> and promotes flowering under field conditions. In addition, PvSVP1 may be a nuclear protein which interacts with PvVRN1 and PvMADS56 on the yeast two-hybrid and BiFC systems. Our study suggests that <italic>PvSVP1</italic> may play a vital role in flowering time and development by interacting with PvVRN1 and PvMADS56 in the nucleus. Furthermore, this study paves the way toward understanding the complex flowering process of bamboo. </p>.

Sequencing and Assembly of Polyploid Genomes

Polyploidy has been observed throughout major eukaryotic clades and has played a vital role in the evolution of angiosperms. Recent polyploidizations often result in highly complex genome structures, posing challenges to genome assembly and phasing. Recent advances in sequencing technologies and genome assembly algorithms have enabled high-quality, near-complete chromosome-level assemblies of polyploid genomes. Advances in novel sequencing technologies include highly accurate single-molecule sequencing with HiFi reads, chromosome conformation capture with Hi-C technique, and linked reads sequencing. Additionally, new computational approaches have also significantly improved the precision and reliability of polyploid genome assembly and phasing, such as HiCanu, hifiasm, ALLHiC, and PolyGembler. Herein, we review recently published polyploid genomes and compare the various sequencing, assembly, and phasing approaches that are utilized in these genome studies. Finally, we anticipate that accurate and telomere-to-telomere chromosome-level assembly of polyploid genomes could ultimately become a routine procedure in the near future.

Genome-wide identification and functional analysis of silicon transporter family genes in moso bamboo (Phyllostachys edulis)

Silicon (Si) has crucial effects on plant development and stress resistance. Silicon transporters regulate Si absorption, transport, and distribution in plants. In this study, we identified and characterized the Si transporter gene family of moso bamboo (Phyllostachys edulis) and cloned seven putative Si transporter genes. Moso bamboo Si transporters contain conserved functional domains that mediate the accumulation of considerable amounts of Si. The analysis of gene duplication patterns and divergence times suggested that the expansion of the moso bamboo Si transporter family was mainly due to segmental duplications. The expression of moso bamboo Si transporter genes, which varied among organs, was significantly modulated by Si treatments. The subcellular localization analysis showed that Si transporters are plasma membrane proteins. The Si content increased in transgenic Arabidopsis overexpressing PeLsi1-1 or PeLsi1-2, which affected vegetative and reproductive growth. Our single-particle tracking analysis revealed the four diffusion modes of PeLsi1-1 on the plasma membrane. Moreover, the particle velocity, dwell time, and motion range of PeLsi1-1 decreased in response to Si treatments. The results of this study will further clarify the molecular mechanisms underlying Si absorption and accumulation in bamboo plants.

Transcriptome-wide analysis reveals core transcriptional regulators associated with culm development and variation in Dendrocalamus sinicus, the strongest woody bamboo in the world

Transcription factors (TFs) play indispensable roles in plant development and stress responses. As the largest woody bamboo species in the world, Dendrocalamus sinicus is endemic to Yunnan Province, China, and possesses two natural variants characterized by culm shape, namely straight or bent culms. Understanding the transcriptional regulation network of D. sinicus provides a unique opportunity to clarify the growth and development characteristics of woody bamboos. In this study, 10,236 TF transcripts belonging to 57 families were identified from transcriptome data of two variants at different developmental stages, from which we constructed a transcriptional regulatory network and unigene-coding protein-TFs interactive network of culm development for this attractive species. Gene function enrichment analysis revealed that hormone signaling and MAPK signaling pathways were two most enriched pathways in TF-regulated network. Based on PPI analysis, 50 genes interacting with nine TFs were screened as the core regulation components related to culm development. Of them, 18 synergistic genes of seven TFs, including nuclear cap-binding protein subunit 1, transcription factor GTE9-like, and ATP-dependent DNA helicase DDX11 isoform X1, involved in culm-shape variation. Most of these genes would interact with MYB, C3H, and ARF transcription factors. Six members with two each from ARF, C3H, and MYB transcription factor families and six key interacting genes (IAA3, IAA19, leucine-tRNA ligase, nuclear cap-binding protein subunit 1, elongation factor 2, and coiled-coil domain-containing protein 94) cooperate with these transcription factors were differentially expressed at development stage of young culms, and were validated by quantitative PCR. Our results represent a crucial step towards understanding the regulatory mechanisms of TFs involved in culm development and variation of D. sinicus.

Rapid growth of Moso bamboo (Phyllostachys edulis): Cellular roadmaps, transcriptome dynamics, and environmental factors

Moso bamboo (Phyllostachys edulis) shows remarkably rapid growth (114.5 cm/day), but the underlying biological mechanisms remain unclear. After examining more than 12,750 internodes from more than 510 culms from 17 Moso populations, we identified internode 18 as a representative internode for rapid growth. This internode includes a 2-cm cell division zone (DZ), a cell elongation zone up to 12 cm, and a secondary cell wall (SCW) thickening zone. These zones elongated 11.8 cm, produced approximately 570,000,000 cells, and deposited ∼28 mg g-1 dry weight (DW) lignin and ∼44 mg g-1 DW cellulose daily, far exceeding vegetative growth observed in other plants. We used anatomical, mathematical, physiological, and genomic data to characterize development and transcriptional networks during rapid growth in internode 18. Our results suggest that (1) gibberellin may directly trigger the rapid growth of Moso shoots, (2) decreased cytokinin and increased auxin accumulation may trigger cell DZ elongation, and (3) abscisic acid and mechanical pressure may stimulate rapid SCW thickening via MYB83L. We conclude that internode length involves a possible tradeoff mediated by mechanical pressure caused by rapid growth, possibly influenced by environmental temperature and regulated by genes related to cell division and elongation. Our results provide insight into the rapid growth of Moso bamboo.

Identification of the 14-3-3 Gene Family in Bamboo and Characterization of Pe14-3-3b Reveals Its Potential Role in Promoting Growth

The 14-3-3 protein family plays an important role in regulating plant growth and development. The genes of the 14-3-3 family have been reported in multiple species. However, little is known about the 14-3-3 gene family in bamboo. In this study, a total of 58 genes belonging to the 14-3-3 family were identified in three representative bamboo species, i.e., Olyra latifolia, Phyllostachys edulis, and Bonia amplexicaulis, whose encoding proteins were grouped into ε and non-ε groups by phylogeny analysis with 14-3-3 proteins from Arabidopsis thaliana and Oryza sativa. The 14-3-3s had diverse gene structures and motif characteristics among the three bamboo species. Collinearity analysis suggested that the genes of the 14-3-3 family in bamboo had undergone a strong purification selection during evolution. Tissue-specific expression analysis showed the expression of Pe14-3-3s varied in different tissues of P. edulis, suggesting that they had functional diversity during growth and development. Co-expression analysis showed that four Pe14-3-3s co-expressed positively with eight ribosomal genes. Yeast two-hybrid (Y2H) assays showed that Pe14-3-3b/d could interact with Pe_ribosome-1/5/6, and qPCR results demonstrated that Pe14-3-3b/d and Pe_ribosome-1/5/6 had similar expression trends with the increase in shoot height, which further confirmed that they would work together to participate in the shoot growth and development of bamboo. Additionally, the transgenic Arabidopsis plants overexpressing Pe14-3-3b had longer roots, a larger stem diameter, an earlier bolting time and a faster growth rate than wild-type Arabidopsis, indicating that Pe14-3-3b acted as a growth promoter. Our results provide comprehensive information on 14-3-3 genes in bamboo and highlight Pe14-3-3b as a potential target for bamboo improvement.

A comprehensive analysis of the floral transition in ma bamboo (Dendrocalamus latiflorus) reveals the roles of DlFTs involved in flowering

Bamboo has a unique flowering characteristics of long and unpredictable vegetative period, which differs from annual herbs and perennial woody plants. In order to understand the molecular regulatory mechanism of bamboo flowering, a comprehensive study was conducted in ma bamboo (Dendrocalamus latiflorus Munro), including morphological, physiological and transcriptiome analyses. Differentially expressed genes related to the flowering pathway were identified by comparative transcriptome analysis. DlFT1, a homologous gene of FT/Hd3a, was significantly upregulated in flowering bamboo. Direct differentiation of spikelets from calli occurred and the downstream gene AP1 was upregulated in the transgenic bamboo overexpressing DlFT1. Transgenic rice overexpressing DlFT1 showed a strong early flowering phenotype. DlFT1 and DlTFL1 could interact with DlFD, and DlTFL1 delayed flowering. It is presumed that DlTFL1 plays an antagonistic role with DlFT1 in ma bamboo. In addition, the expression of DlFT1 was regulated by DlCO1, indicating that a CO-FT regulatory module might exist in ma bamboo. These results suggest that DlFT1 is a florigen candidate gene with conservative function in promoting flowering. Interestingly, the results have shown for the first time that DlFT2 can specifically interact with E3 ubiquitin ligase WAV3, while DlFT3 transcripts are mainly nonsense splicing. These findings provide better understanding of the roles of the florigen gene in bamboo and lay a theoretical basis for regulating bamboo flowering in the future.

Comprehensive Analyses of Simple Sequence Repeat (SSR) in Bamboo Genomes and Development of SSR Markers with Peroxidase Genes

Simple sequence repeats (SSRs) are one of the most important molecular markers, which are widespread in plants. Bamboos are important forest resources worldwide. Here, the comprehensive identification and comparative analysis of SSRs were performed in three woody and two herbaceous bamboo species. Altogether 567,175 perfect SSRs and 71,141 compound SSRs were identified from 5737.8 Mb genome sequences of five bamboo species. Di-nucleotide SSRs were the most predominant type, with an average of ~50,152.2 per species. Most SSRs were located in intergenic regions, while those located in genic regions were relatively less. Moreover, the results of annotation distribution indicated that terms with P450, peroxidase and ATP-binding cassette transporter related to lignin biosynthesis might play important roles in woody and herbaceous bamboos under the mediation of SSRs. Furthermore, the peroxidase gene family consisted of a large number of genes containing SSRs was selected for the evolutionary relationship analysis and SSR markers development. Fifteen SSR markers derived from peroxidase family genes of Phyllostachys edulis were identified as polymorphic in 34 accessions belonging to seven genera in Bambusoideae. These results provided a comprehensive insight of SSR markers into bamboo genomes, which would facilitate bamboo research related to comparative genomics, evolution and marker-assisted selection.

PeVDE, a violaxanthin de-epoxidase gene from moso bamboo, confers photoprotection ability in transgenic Arabidopsis under high light

Plants employ an array of photoprotection mechanisms to alleviate the harmful effects of high light intensity. The violaxanthin cycle, which is associated with non-photochemical quenching (NPQ), involves violaxanthin de-epoxidase (VDE), and zeaxanthin epoxidase (ZEP) and is one of the most rapid and efficient mechanisms protecting plants under high light intensity. Woody bamboo is a class of economically and ecologically important evergreen grass species widely distributed in tropical and subtropical areas. However, the function of VDE in bamboo has not yet been elucidated. In this study, we found that high light intensity increased NPQ and stimulated the de-epoxidation of violaxanthin cycle components in moso bamboo (Phyllostachys edulis), whereas, samples treated with the VDE inhibitor (dithiothreitol) exhibited lower NPQ capacity, suggesting that violaxanthin cycle plays an important role in the photoprotection of bamboo. Further analysis showed that not only high light intensity but also extreme temperatures (4 and 42°C) and drought stress upregulated the expression of PeVDE in bamboo leaves, indicating that PeVDE is induced by multiple abiotic stresses. Overexpression of PeVDE under the control of the CaMV 35S promoter in Arabidopsis mutant npq1 mutant could rescue its NPQ, indicating that PeVDE functions in dissipating the excess absorbed light energy as thermal energy in bamboo. Moreover, compared with wild-type (Col-0) plants, the transgenic plants overexpressing PeVDE displayed enhanced photoprotection ability, higher NPQ capacity, slower decline in the maximum quantum yield of photosystem II (F _v /F _m ) under high light intensity, and faster recovery under optimal conditions. These results suggest that PeVDE positively regulates the response to high light intensity in bamboo plants growing in the natural environment, which could improve their photoprotection ability through the violaxanthin cycle and NPQ.

Single-base methylome analysis reveals dynamic changes of genome-wide DNA methylation associated with rapid stem growth of woody bamboos

CG and CHG methylation levels in the rapid shoot growth stages (ST2-ST4) of woody bamboos were obviously decreased, which might regulate the internode elongation during rapid shoot growth, while CHH methylation was strongly associated with shoot developmental time or age. DNA methylation plays a critical role in the regulation of plant growth and development. Woody bamboos have a unique trait of rapid stem growth resulted from internode elongation at the shooting period. However, it is still unclear whether DNA methylation significantly controls the bamboo rapid stem growth. Here we present whole-genome DNA methylation profiles of the paleotropical woody bamboo Bonia amplexicaulis at five newly defined stages of shoot growth, named ST1-ST5. We found that CG and CHG methylation levels in the rapid shoot growth stages (ST2-ST4) were significantly lower than in the incubation (ST1) and plateau stages (ST5). The changes in methylation levels mainly occurred in flanking regions of genes and gene body regions, and 23647 differentially methylated regions (DMRs) were identified between ST1 and rapid shoot growth stages (ST2-ST4). Combined with transcriptome analysis, we found that DMR-related genes enriched in the auxin and jasmonic acid (JA) signal transduction, and other pathways closely related to plant growth. Intriguingly, CHH methylation was not involved in the rapid shoot growth, but strongly associated with shoot developmental time by gradually accumulating in transposable elements (TEs) regions. Overall, our results reveal the importance of DNA methylation in regulating the bamboo rapid shoot growth and suggest a role of DNA methylation associated with development time or age in woody bamboos.

The numbers of fungi: contributions from traditional taxonomic studies and challenges of metabarcoding

The global diversity of fungi has been estimated using several different approaches. There is somewhere between 2–11 million estimated species, but the number of formally described taxa is around 150,000, a tiny fraction of the total. In this paper, we examine 12 ascomycete genera as case studies to establish trends in fungal species descriptions, and introduce new species in each genus. To highlight the importance of traditional morpho-molecular methods in publishing new species, we introduce novel taxa in 12 genera that are considered to have low species discovery. We discuss whether the species are likely to be rare or due to a lack of extensive sampling and classification. The genera are Apiospora, Bambusicola, Beltrania, Capronia, Distoseptispora, Endocalyx, Neocatenulostroma, Neodeightonia, Paraconiothyrium, Peroneutypa, Phaeoacremonium and Vanakripa. We discuss host-specificity in selected genera and compare the number of species epithets in each genus with the number of ITS (barcode) sequences deposited in GenBank and UNITE. We furthermore discuss the relationship between the divergence times of these genera with those of their hosts. We hypothesize whether there might be more species in these genera and discuss hosts and habitats that should be investigated for novel species discovery.

Genome-Wide Identification, Expansion, and Evolution Analysis of Homeobox Gene Family Reveals TALE Genes Important for Secondary Cell Wall Biosynthesis in Moso Bamboo (Phyllostachys edulis)

The TALE gene family is a subfamily of the homeobox gene family and has been implicated in regulating plant secondary growth. However, reports about the evolutionary history and function of the TALE gene family in bamboo are limited. Here, the homeobox gene families of moso bamboo Olyra latifolia and Bonia amplexicaulis were identified and compared. Many duplication events and obvious expansions were found in the TALE family of woody bamboo. PhTALEs were found to have high syntenies with TALE genes in rice. Through gene co-expression analysis and quantitative real-time PCR analysis, the candidate PhTALEs were thought to be involved in regulating secondary cell wall development of moso bamboo during the fast-growing stage. Among these candidate PhTALEs, orthologs of OsKNAT7, OSH15, and SH5 in moso bamboo may regulate xylan synthesis by regulating the expression of IRX-like genes. These results suggested that PhTALEs may participate in the secondary cell wall deposition in internodes during the fast-growing stage of moso bamboo. The expansion of the TALE gene family may be implicated in the increased lignification of woody bamboo when divergent from herbaceous bamboos.

A well-supported nuclear phylogeny of Poaceae and implications for the evolution of C4 photosynthesis

Poaceae (the grasses) includes rice, maize, wheat, and other crops, and is the most economically important angiosperm family. Poaceae is also one of the largest plant families, consisting of over 11 000 species with a global distribution that contributes to diverse ecosystems. Poaceae species are classified into 12 subfamilies, with generally strong phylogenetic support for their monophyly. However, many relationships within subfamilies, among tribes and/or subtribes, remain uncertain. To better resolve the Poaceae phylogeny, we generated 342 transcriptomic and seven genomic datasets; these were combined with other genomic and transcriptomic datasets to provide sequences for 357 Poaceae species in 231 genera, representing 45 tribes and all 12 subfamilies. Over 1200 low-copy nuclear genes were retrieved from these datasets, with several subsets obtained using additional criteria, and used for coalescent analyses to reconstruct a Poaceae phylogeny. Our results strongly support the monophyly of 11 subfamilies; however, the subfamily Puelioideae was separated into two non-sister clades, one for each of the two previously defined tribes, supporting a hypothesis that places each tribe in a separate subfamily. Molecular clock analyses estimated the crown age of Poaceae to be ∼101 million years old. Ancestral character reconstruction of C₃/C₄ photosynthesis supports the hypothesis of multiple independent origins of C₄ photosynthesis. These origins are further supported by phylogenetic analysis of the ppc gene family that encodes the phosphoenolpyruvate carboxylase, which suggests that members of three paralogous subclades (ppc-aL1a, ppc-aL1b, and ppc-B2) were recruited as functional C₄ppc genes. This study provides valuable resources and a robust phylogenetic framework for evolutionary analyses of the grass family.

Allele-aware chromosome-scale assembly of the allopolyploid genome of hexaploid Ma bamboo (Dendrocalamus latiflorus Munro)

Dendrocalamus latiflorus Munro is a woody clumping bamboo with rapid shoot growth. Both genetic transformation and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) gene editing techniques are available for D. latiflorus, enabling reverse genetic approaches. Thus, D. latiflorus has the potential to be a model bamboo species. However, the genome sequence of D. latiflorus has remained unreported due to its polyploidy and large genome size. Here, we sequenced the D. latiflorus genome and assembled it into three allele-aware subgenomes (AABBCC), representing the largest genome of a major bamboo species. We assembled 70 allelic chromosomes (2, 737 Mb) for hexaploid D. latiflorus using both single-molecule sequencing from the Pacific Biosciences (PacBio) Sequel platform and chromosome conformation capture sequencing (Hi-C). Repetitive sequences comprised 52.65% of the D. latiflorus genome. We annotated 135 231 protein-coding genes in the genome based on transcriptomes from eight different tissues. Transcriptome sequencing using RNA-Seq and PacBio single-molecule real-time long-read isoform sequencing revealed highly differential alternative splicing (AS) between non-abortive and abortive shoots, suggesting that AS regulates the abortion rate of bamboo shoots. This high-quality hexaploid genome and comprehensive strand-specific transcriptome datasets for this Poaceae family member will pave the way for bamboo research using D. latiflorus as a model species.