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

Specimen Identification Through Multilocus Species Tree Constructed From Single-Copy Orthologs (SCOs): A Case Study in Cymbidium Subgenus Jensoa

Standard barcodes and ultra-barcode encounter significant challenges when delimiting and discriminating closely related species characterized by deep coalescence, hybrid speciation, gene flow, or low sequence variation. Single-copy orthologs (SCOs) have been widely recognized as standardized nuclear markers in metazoan DNA taxonomy, yet their application in plant taxonomy remains unexplored. This study evaluates the efficacy of SCOs for identifying recently diverged species within the Cymbidium subgenus Jensoa, where ultra-barcodes have previously shown limited resolution. Remarkably, over 90% of the 9094 targeted reference SCOs, inferred from three Cymbidium genomes, were successfully retrieved for all 11 representative species in subg. Jensoa using ALiBaSeq at a minimal 5× depth from whole genome shotgun sequences. The species tree, reconstructed from multiple refined SCO matrices under the coalescent model, effectively distinguished all species and identified mislabeled or misidentified specimens. The comprehensive and refined SCO matrices produced by our pipeline not only enhance phylogenetic analysis but also improve the precision of species diagnosis. Additionally, biparentally inherited SCOs, serving as multi-locus markers, not only augment the effectiveness of DNA barcoding but also support a transition to multi-locus, species-tree-based specimen assignment strategies.

Phylogenomics and biogeography of Guadua: Insights into a neotropical woody bamboo genus

PREMISE

The genus Guadua includes some of the largest bamboo species in the neotropics, with certain species of significant economic importance and used since pre-Hispanic times to build houses using the traditional bajareque technique. Guadua species are distributed from Mexico to South America. The potential monophyly of this genus has been suggested based on plastid markers and limited sampling. Here we included more species and nuclear data to study the morphological classification of Guadua and to reconstruct its ancestral area distribution.

METHODS

Samples were collected for 16 Guadua taxa (13 species) to use nuclear single nucleotide polymorphisms (SNPs) derived from restriction-site associated DNA sequencing (RADSeq) data to construct a phylogenetic hypothesis using maximum likelihood and coalescent methods. We estimated divergence times using the RelTime method and reconstructed ancestral geographic areas using S-DEC analysis.

RESULTS

The SNP data supported the monophyly of the genus Guadua and the existence of two distinct clades, Amplexifolia and Angustifolia. The origin of Guadua was estimated as approximately 10.47 million years ago in Mesoamerica. The biogeographic distribution of Guadua can likely be explained by a combination of dispersal and vicariant events.

CONCLUSIONS

Our study sheds new light on the phylogenetic relationships within Guadua and on its evolutionary history and biogeography, enhancing our understanding of its diversification and distribution across various biogeographic regions.

The wheat NLR pair RXL/Pm5e confers resistance to powdery mildew

Powdery mildew poses a significant threat to global wheat production and most cloned and deployed resistance genes for wheat breeding encode nucleotide-binding and leucine-rich repeat (NLR) immune receptors. Although two genetically linked NLRs function together as an NLR pair have been reported in other species, this phenomenon has been relatively less studied in wheat. Here, we demonstrate that two tightly linked NLR genes, RXL and Pm5e, arranged in a head-to-head orientation, function together as an NLR pair to mediate powdery mildew resistance in wheat. The resistance function of the RXL/Pm5e pair is validated by mutagenesis, gene silencing, and gene-editing assays. Interestingly, both RXL and Pm5e encode atypical NLRs, with RXL possessing a truncated NB-ARC (nucleotide binding adaptor shared by APAF-1, plant R proteins and CED-4) domain and Pm5e featuring an atypical coiled-coil (CC) domain. Notably, RXL and Pm5e lack an integrated domain associated with effector recognition found in all previously reported NLR pairs. Additionally, RXL and Pm5e exhibit a preference for forming hetero-complexes rather than homo-complexes, highlighting their cooperative role in disease resistance. We further show that the CC domain of Pm5e specifically suppresses the hypersensitive response induced by the CC domain of RXL through competitive interaction, revealing regulatory mechanisms within this NLR pair. Our study sheds light on the molecular mechanism underlying RXL/Pm5e-mediated powdery mildew resistance and provides a new example of an NLR pair in wheat disease resistance.

A 23-million-year record of morphological evolution within Neotropical grass pollen

Grass-dominated biomes in South America comprise c. 20 million years of history, yet their evolution and underlying drivers remain poorly understood. Here we apply a novel approach that combines scanning electron microscopy imaging with computational analysis to quantify the morphometrics of grass (Poaceae) pollen micro-ornamentation from the Neotropics since the Early Miocene (23 million years ago). Three spatial-temporal pollen sets were assembled to further elucidate the variation and evolutionary traits of grasses through space and time. Our results reveals that three spatial-temporal pollen groups occupy unique, partially overlapping regions of their exine morphospace. The direction of this shift is consistent over time, progressing towards less dense ornamentation. Interestingly, the extent of the occupied morphospace did not vary significantly. This is the first time that the true morphological variation in Poaceae pollen micro-ornamentation becomes apparent through time. We hypothesize that changes in grass pollen exine since the Early Miocene were driven by evolutionary processes (evolutionary drift and/or directional selection), and potentially migration at the continental scale. The high diversity in pollen micro-ornamentation is likely related to their evolutionary success in the Neogene.

A wheat tandem kinase and NLR pair confers resistance to multiple fungal pathogens

Tandem kinase proteins underlie the innate immune systems of cereal plants, but how they initiate plant immune responses remains unclear. This report identifies wheat protein wheat tandem NBD 1 (WTN1), a noncanonical nucleotide-binding leucine-rich repeat (NLR) receptor featuring tandem nucleotide binding adaptor shared by APAF-1, plant R proteins, and CED-4 (NB-ARC) domains, required for WTK3-mediated disease resistance. Both WTK3 and its allelic variant Rwt4-known for conferring resistance to wheat powdery mildew and blast, respectively-are capable of recognizing the blast effector PWT4. They activate WTN1 to form calcium-permeable channels, akin to ZAR1 and Sr35. Thus, tandem kinase proteins and their associated NLRs operate as "sensor-executor" pairs against fungal pathogens. Additionally, evolutionary analyses reveal a coevolutionary trajectory of the tandem kinase-NLR module, highlighting their cooperative role in triggering plant immunity.

Constraint of accessible chromatins maps regulatory loci involved in maize speciation and domestication

Comparative genomic studies can identify genes under evolutionary constraint or specialized for trait innovation. Growing evidence suggests that evolutionary constraint also acts on non-coding regulatory sequences, exerting significant impacts on fitness-related traits, although it has yet to be thoroughly explored in plants. Using the assay for transposase-accessible chromatin by sequencing (ATAC-seq), we profile over 80,000 maize accessible chromatin regions (ACRs), revealing that ACRs evolve faster than coding genes, with about one-third being maize-specific and regulating genes associated with speciation. We highlight the role of transposable elements (TEs) in driving intraspecific innovation of ACRs and identify hundreds of candidate ACRs potentially involved in transcriptional rewiring during maize domestication. Additionally, we demonstrate the importance of accessible chromatin in maintaining subgenome dominance and controlling complex trait variations. This study establishes a framework for analyzing the evolutionary trajectory of plant regulatory sequences and offers candidate loci for downstream exploration and application in maize breeding.

Genome assembly of Stewartia sinensis reveals origin and evolution of orphan genes in Theaceae

Orphan genes play crucial roles in diverse biological processes, but the evolutionary trajectories and functional divergence remain largely unexplored. The Theaceae family, including the economically and culturally important tea plant, offers a distinctive model to examine these aspects. Here, we integrated Nanopore long-read sequencing, Illumina short-read sequencing, and Hi-C methods to decode a pseudo-chromosomal genome assembly of Stewartia sinensis, from the earliest-diverging tribe of Theaceae, spanning 2.95 Gb. Comparative genomic analysis revealed the absence of recent whole-genome duplication events in the Theaceae ancestor, highlighting tandem duplications as the predominant mechanism of gene expansion. We identified 31,331 orphan genes, some of which appear to have ancient origins, suggesting early emergence with frequent gains and losses, while others seem more specific and recent. Notably, orphan genes are distinguished by shorter lengths, fewer exons and functional domains compared to genes that originate much earlier, like transcription factors. Moreover, tandem duplication contributes significantly to the adaptive evolution and characteristic diversity of Theaceae, and it is also a major mechanism driving the origination of orphan genes. This study illuminates the evolutionary dynamics of orphan genes, providing a valuable resource for understanding the origin and evolution of tea plant flavor and enhancing genetic breeding efforts.

Loss of Pleiotropic Regulatory Functions in Tannin1, the Sorghum Ortholog of Arabidopsis Master Regulator TTG1

Transcriptional master regulators are often targeted to improve plant traits, but antagonistic pleiotropic effects of these regulators can hamper this approach. The Myb-bHLH-WDR (MBW) complex is a broadly conserved transcriptional regulator affecting pigmentation, biotic stress resistance, and abiotic stress tolerance. We investigated the function of sorghum grain pigmentation regulator Tannin1, the ortholog of Arabidopsis pleiotropic WD40 regulator TTG1, to test for conserved pleiotropic regulatory effects and to better understand the evolution of the MBW complex in Poaceae. We characterized genome-wide differential expression of leaf tissue using RNA sequencing in near-isogenic lines (NILs) that contrasted wildtype Tan1 and loss-of-function tan1-b alleles, under optimal temperature and chilling stress. Notably, Gene Ontology analyses revealed no pathways with differential expression between Tan1 and tan1-b NILs, suggesting that, in contrast to Arabidopsis TTG1, Tannin1 has no pleiotropic regulatory role in leaves. Further, NILs had no visible difference in anthocyanin pigmentation, and no genes with known or expected function in flavonoid synthesis were differentially expressed. Genome-wide, only 18 total genes were differentially expressed between NILs, with six of these genes located inside the NIL introgression region, an observation most parsimoniously explained by cis-regulatory effects unrelated to Tannin1 regulation. Comparing our findings with known function of TTG1 orthologs in maize, rice, and Arabidopsis, we conclude that pleiotropic regulatory function in leaf tissue was likely lost in panicoid grass evolution before the sorghum-maize split. These findings inform future molecular breeding of MBW regulated traits and highlight the benefit of subfunctionalization to relieve pleiotropic constraints.

Historic rewiring of grass flowering time pathways and implications for crop improvement under climate change

Grasses are fundamental to human survival, providing a large percentage of our calories, fuel, and fodder for livestock, and an enormous global carbon sink. A particularly important part of the grass plant is the grain-producing inflorescence that develops in response to both internal and external signals that converge at the shoot tip to influence meristem behavior. Abiotic signals that trigger reproductive development vary across the grass family, mostly due to the unique ecological and phylogenetic histories of each clade. The time it takes a grass to flower has implications for its ability to escape harsh environments, while also indirectly affecting abiotic stress tolerance, inflorescence architecture, and grain yield. Here, we synthesize recent insights into the evolution of grass flowering time in response to past climate change, particularly focusing on genetic convergence in underlying traits. We then discuss how and why the rewiring of a shared ancestral flowering pathway affects grass yields, and outline ways in which researchers are using this and other information to breed higher yielding, climate-proof cereal crops.

Phylogenomic analyses shed new light on the spatiotemporal evolution of global larches: Implications for the dynamics of boreal forests

As the Earth warms, understanding the long-term dynamics of forest ecosystems is essential for guiding forest management and biodiversity conservation. Insights from past dynamics may provide valuable lessons for managing today's forests. Here, we investigated the spatiotemporal evolution of global larches to gain further insights into how boreal forests change over time. We first reconstructed a highly resolved and robust phylogeny of Larix covering all widely recognized species, using both transcriptome-based 1,301 orthologous genes (OGs) and plastid genomes. In sharp contrast to previous studies, an unexpected deep split between the circumboreal and Qinghai-Tibetan Plateau (QTP) larches was revealed in our study. Within each lineage, two geographically distinct clades were further resolved. Biogeographical analyses suggest that Larix might have an origin of Eocene in high-latitude uplands, and during the Miocene, all extant species have appeared. Cenozoic climate- and orogeny-triggered vicariance likely played a major role in the divergence of global larches. Our results also demonstrate that the proto-boreal forest biome may have a relatively old origin back to the early Miocene, and significant winnowing and species alteration would have occurred as the climate shifted to much colder and drier conditions during the Neogene. Ecological niche analyses show various responses of the circumboreal and QTP larches under different climate scenarios, but both lineages are negatively impacted by warming climates. These findings have important conservation implications given the sensitivity of boreal forests in the face of global warming. Our work further emphasizes the importance of a solid phylogenetic framework for evolutionary and biogeographical inferences.

Grass Phylogeny Working Group III. A nuclear phylogenomic tree of grasses (Poaceae) recovers current classification despite gene tree incongruence

Grasses (Poaceae) comprise c. 11 800 species and are central to human livelihoods and terrestrial ecosystems. Knowing their relationships and evolutionary history is key to comparative research and crop breeding. Advances in genome-scale sequencing allow for increased breadth and depth of phylogenomic analyses, making it possible to infer a new reference species tree of the family. We inferred a comprehensive species tree of grasses by combining new and published sequences for 331 nuclear genes from genome, transcriptome, target enrichment and shotgun data. Our 1153-tip tree covers 79% of grass genera (including 21 genera sequenced for the first time) and all but two small tribes. We compared it to a newly inferred 910-tip plastome tree. We recovered most of the tribes and subfamilies previously established, despite pervasive incongruence among nuclear gene trees. The early diversification of the PACMAD clade could represent a hard polytomy. Gene tree-species tree reconciliation suggests that reticulation events occurred repeatedly. Nuclear-plastome incongruence is rare, with very few cases of supported conflict. We provide a robust framework for the grass tree of life to support research on grass evolution, including modes of reticulation, and genetic diversity for sustainable agriculture.

Exaptation of ancestral cell-identity networks enables C4 photosynthesis

C4 photosynthesis is used by the most productive plants on the planet, and compared with the ancestral C3 pathway, it confers a 50% increase in efficiency1. In more than 60 C4 lineages, CO2 fixation is compartmentalized between tissues, and bundle-sheath cells become photosynthetically activated2. How the bundle sheath acquires this alternate identity that allows efficient photosynthesis is unclear. Here we show that changes to bundle-sheath gene expression in C4 leaves are associated with the gain of a pre-existing cis-code found in the C3 leaf. From single-nucleus gene-expression and chromatin-accessibility atlases, we uncover DNA binding with one finger (DOF) motifs that define bundle-sheath identity in the major crops C3 rice and C4 sorghum. Photosynthesis genes that are rewired to be strongly expressed in the bundle-sheath cells of C4 sorghum acquire cis-elements that are recognized by DOFs. Our findings are consistent with a simple model in which C4 photosynthesis is based on the recruitment of an ancestral cis-code associated with bundle-sheath identity. Gain of such elements harnessed a stable patterning of transcription factors between cell types that are found in both C3 and C4 leaves to activate photosynthesis in the bundle sheath. Our findings provide molecular insights into the evolution of the complex C4 pathway, and might also guide the rational engineering of C4 photosynthesis in C3 crops to improve crop productivity and resilience3,4.

Comparative analysis of chloroplast genomes and phylogenetic relationships in the endemic Chinese bamboo Gelidocalamus (Bambusoideae)

Introduction

Gelidocalamus Wen is a small yet taxonomically challenging genus within the Arundinarieae tribe. Recent molecular studies have suggested it may not be monophyletic. However, limited species sampling and insufficient molecular marker information have resulted in poorly resolved phylogenetic relationships within this genus.

Methods

The complete chloroplast genomes covering all 16 species and one variant of Gelidocalamus were sequenced, and comparative analyses were conducted. Phylogenetic analyses were performed using different molecular markers, including chloroplast data, the nuclear ribosomal DNA (nrDNA) repeats region, and 29 mitochondrial protein-coding genes. Additionally, the divergence times of Gelidocalamus were estimated to reveal their evolutionary history.

Results

The plastomes of Gelidocalamus ranged in size from 139,500 bp to 139,801 bp, with a total of 137 identified genes, including 90 protein-coding genes, 39 tRNA genes, and 8 rRNA genes. The size of the nrDNA repeats ranged from 5,802 bp to 5,804 bp. Phylogenetic analysis based on chloroplast data revealed that Gelidocalamus is polyphyletic, with different subclades distributed within the IV and V clades. However, phylogenetic analysis based on nrDNA and mitochondrial genes did not effectively resolve the relationships within the genus.

Discussion

Comparative analysis of chloroplast genomes indicated that Gelidocalamus shares a high degree of similarity with closely related genera in terms of chloroplast genome collinearity, codon usage bias, and repetitive sequences. Divergence time estimation suggests that it is a relatively young group, with all members appearing successively over the past four million years. The complex phylogenetic patterns may arise from the rapid radiation of Arundinarieae. This study provides a preliminary foundation for further in-depth research on the phylogeny, genomic structural features, and divergence times of this genus.

Peering through the hedge: Multiple datasets yield insights into the phylogenetic relationships and incongruences in the tribe Lilieae (Liliaceae)

The increasing use of genome-scale data has significantly facilitated phylogenetic analyses, contributing to the dissection of the underlying evolutionary mechanisms that shape phylogenetic incongruences, such as incomplete lineage sorting (ILS) and hybridization. Lilieae, a prominent member of the Liliaceae family, comprises four genera and approximately 260 species, representing 43% of all species within Liliaceae. They possess high ornamental, medicinal and edible values. Yet, no study has explored the validity of various genome-scale data in phylogenetic analyses within this tribe, nor have potential evolutionary mechanisms underlying its phylogenetic incongruences been investigated. Here, transcriptome, Angiosperms353, plastid and mitochondrial data, were collected from 50 to 93 samples of Lilieae, covering all four recognized genera. Multiple datasets were created and used for phylogenetic analyses based on concatenated and coalescent-based methods. Evolutionary rates of different datasets were calculated, and divergence times were estimated. Various approaches, including coalescence simulation, Quartet Sampling (QS), calculation of concordance factors (gCF and sCF), as well as MSCquartets and reticulate network inference, were carried out to infer the phylogenetic discordances and analyze their underlying mechanisms using a reduced 33-taxon dataset. Despite extensive phylogenetic discordances among gene trees, robust phylogenies were inferred from nuclear and plastid data compared to mitochondrial data, with lower synonymous substitution detected in mitochondrial genes than in nuclear and plastid genes. Significant ILS was detected across the phylogeny of Lilieae, with clear evidence of reticulate evolution identified. Divergence time estimation indicated that most of lineages in Lilieae diverged during a narrow time frame (ranging from 5.0 Ma to 10.0 Ma), consistent with the notion of rapid radiation evolution. Our results suggest that integrating transcriptomic and plastid data can serve as cost-effective and efficient tools for phylogenetic inference and evolutionary analysis within Lilieae, and Angiosperms353 data is also a favorable choice. Mitochondrial data are more suitable for phylogenetic analyses at higher taxonomic levels due to their stronger conservation and lower synonymous substitution rates. Significant phylogenetic incongruences detected in Lilieae were caused by both incomplete lineage sorting (ILS) and reticulate evolution, with hybridization and "ghost introgression" likely prevalent in the evolution of Lilieae species. Our findings provide new insights into the phylogeny of Lilieae, enhancing our understanding of the evolution of species in this tribe.

Phylogenomics resolves the backbone of Poales and identifies signals of hybridization and polyploidy

Poales, as one of the largest orders of angiosperm, holds crucial economic and ecological importance. Nevertheless, achieving a consensus topology has been challenging in previous studies due to limited molecular data and sparse taxon sampling. The uneven distribution of species diversity among families and the factors leading to elevated species richness in certain lineages have also been subjects of ongoing discussion and investigation. In this study, we conducted a comprehensive sampling, including representatives from all 14 families and 85 taxa of Poales, along with five additional outgroups. To reconstruct the phylogeny of Poales, we employed a combination of coalescent and concatenation methods on three nuclear gene sets (1093, 491, 143) and one plastid gene set (53), which were inferenced from genomic data. We also conducted phylogenetic hypothesis analyses to evaluate two major conflicting nodes detected in phylogenetic analyses. As a result, we successfully resolved the backbone of Poales and provided a timeline for its evolutionary history. We recovered the sister relationship between Typhaceae and Bromeliaceae as the earliest diverging families within Poales. The clade consisting of Ecdeiocoleaceae and Joinvilleaceae was recovered as the sister group of Poaceae. Within the xyrid clade, Mayacaceae and Erioaculaceae + Xyridaceae successively diverged along the backbone of Poales. The topology of [Aristidoideae, ((Micrairoideae, Panicoideae), (Arundinoideae, (Chloridoideae, Danthonioideae)))] within the PACMAD clade has received strong support from multiple findings. We also delved into the underlying biological factors that contributed to the conflicting nodes observed in the phylogenetic analysis. Apart from the uncertainty regarding the sister group of Poaceae caused by cytonuclear discordance, frequent hybridization and polyploidy may have contributed to other conflicting nodes. We identified 26 putative whole-genome duplication (WGD) events within Poales. However, apart from the σ-WGD and the ρ-WGD, we did not observe any potential polyploid events that could be directly linked to the species diversification in specific lineages. Furthermore, there was a significant increase in the net diversification rate of Poales following the K-Pg boundary.

Three-fold expansion of the genetic diversity of blunerviruses through plant (meta)transcriptome data-mining

Kitaviruses are plant-infecting, non-systemic disease-causing viruses with multipartite single-stranded RNA genomes. Despite their importance, knowledge on kitaviruses is limited in comparison with other plant virus groups, mainly because of the lesser number of identified and characterized kitaviruses and their isolates. In the present study, we explored plant (meta)transcriptome data available in public domain and identified genome sequences of eighteen putative novel blunerviruses in eighteen plant species, including four gymnosperm and four monocot species. Four RNA segments (RNAs 1-4) of eleven identified viruses were recovered, whilst at least two RNA segments were recovered for the remaining viruses. Phylogenetic analysis grouped the identified viruses with known blunerviruses. Based on genome organization, sequence identities of encoded proteins with known blunerviruses and phylogeny, the identified viruses are regarded as new members of the genus Blunervirus. The study paves way for initiating further studies on understanding biological properties, economic importance and geographical distribution of identified blunerviruses.

Comparative plastome analysis of Arundinelleae (Poaceae, Panicoideae), with implications for phylogenetic relationships and plastome evolution

BACKGROUND

Arundinelleae is a small tribe within the Poaceae (grass family) possessing a widespread distribution that includes Asia, the Americas, and Africa. Several species of Arundinelleae are used as natural forage, feed, and raw materials for paper. The tribe is taxonomically cumbersome due to a paucity of clear diagnostic morphological characters. There has been scant genetic and genomic research conducted for this group, and as a result the phylogenetic relationships and species boundaries within Arundinelleae are poorly understood.

RESULTS

We compared and analyzed 11 plastomes of Arundinelleae, of which seven plastomes were newly sequenced. The plastomes range from 139,629 base pairs (bp) (Garnotia tenella) to 140,943 bp (Arundinella barbinodis), with a standard four-part structure. The average GC content was 38.39%, but varied in different regions of the plastome. In all, 110 genes were annotated, comprising 76 protein-coding genes, 30 tRNA genes, and four rRNA genes. Furthermore, 539 simple sequence repeats, 519 long repeats, and 10 hyper-variable regions were identified from the 11 plastomes of Arundinelleae. A phylogenetic reconstruction of Panicoideae based on 98 plastomes demonstrated the monophyly of Arundinella and Garnotia, but the circumscription of Arundinelleae remains unresolved.

CONCLUSION

Complete chloroplast genome sequences can improve phylogenetic resolution relative to single marker approaches, particularly within taxonomically challenging groups. All phylogenetic analyses strongly support the monophyly of Arundinella and Garnotia, respectively, but the monophylly of Arundinelleae was not well supported. The intergeneric phylogenetic relationships within Arundinelleae require clarification, indicating that more data is necessary to resolve generic boundaries and evaluate the monophyly of Arundinelleae. A comprehensive taxonomic revision for the tribe is necessary. In addition, the identified hyper-variable regions could function as molecular markers for clarifying phylogenetic relationships and potentially as barcoding markers for species identification in the future.

Chromosome-level genome assemblies reveal genome evolution of an invasive plant Phragmites australis

Biological invasions pose a significant threat to ecosystems, disrupting local biodiversity and ecosystem functions. The genomic underpinnings of invasiveness, however, are still largely unknown, making it difficult to predict and manage invasive species effectively. The common reed (Phragmites australis) is a dominant grass species in wetland ecosystems and has become particularly invasive when transferred from Europe to North America. Here, we present a high-quality gap-free, telomere-to-telomere genome assembly of Phragmites australis consisting of 24 pseudochromosomes and a B chromosome. Fully phased subgenomes demonstrated considerable subgenome dominance and revealed the divergence of diploid progenitors approximately 30.9 million years ago. Comparative genomics using chromosome-level scaffolds for three other lineages and a previously published draft genome assembly of an invasive lineage revealed that gene family expansions in the form of tandem duplications may have contributed to the invasiveness of the lineage. This study sheds light on the genome evolution of Arundinoideae grasses and suggests that genetic drivers, such as gene family expansions and tandem duplications, may underly the processes of biological invasion in plants. These findings provide a crucial step toward understanding and managing the genetic basis of invasiveness in plant species.

EPAD1 Orthologs Play a Conserved Role in Pollen Exine Patterning

The pollen wall protects pollen during dispersal and is critical for pollination recognition. In the Poaceae family, the pollen exine stereostructure exhibits a high degree of conservation with similar patterns across species. However, there remains controversy regarding the conservation of key factors involved in its formation among various Poaceae species. EPAD1, as a gene specific to the Poaceae family, and its orthologous genes play a conserved role in pollen wall formation in wheat and rice. However, they do not appear to have significant functions in maize. To further confirm the conserved function of EPAD1 in Poaceae, we performed an analysis on four EPAD1 orthologs from two distinct sub-clades within the Poaceae family. The two functional redundant barley EPAD1 genes (HvEPAD1 and HvEPAD2) from the BOP clade, along with the single copy of sorghum (SbEPAD1) and millet (SiEPAD1) from the PACMAD clade were examined. The CRISPR-Cas9-generated mutants all exhibited defects in pollen wall formation, consistent with previous findings on EPAD1 in rice and wheat. Interestingly, in barley, hvepad2 single mutant also showed apical spikelets abortion, aligning with a decreased expression level of HvEPAD1 and HvEPAD2 from the apical to the bottom of the spike. Our finding provides evidence that EPAD1 orthologs contribute to Poaceae specific pollen exine pattern formation via maintaining primexine integrity despite potential variations in copy numbers across different species.

Complete mitochondrial genome assembly of Zizania latifolia and comparative genome analysis

Zizania latifolia (Griseb.) Turcz. ex Stapf has been cultivated as a popular aquatic vegetable in China due to its important nutritional, medicinal, ecological, and economic values. The complete mitochondrial genome (mitogenome) of Z. latifolia has not been previously studied and reported, which has hindered its molecular systematics and understanding of evolutionary processes. Here, we assembled the complete mitogenome of Z. latifolia and performed a comprehensive analysis including genome organization, repetitive sequences, RNA editing event, intercellular gene transfer, phylogenetic analysis, and comparative mitogenome analysis. The mitogenome of Z. latifolia was estimated to have a circular molecule of 392,219 bp and 58 genes consisting of three rRNA genes, 20 tRNA genes, and 35 protein-coding genes (PCGs). There were 46 and 20 simple sequence repeats (SSRs) with different motifs identified from the mitogenome and chloroplast genome of Z. latifolia, respectively. Furthermore, 49 homologous fragments were observed to transfer from the chloroplast genome to the mitogenome of Z. latifolia, accounting for 47,500 bp, presenting 12.1% of the whole mitogenome. In addition, there were 11 gene-containing homologous regions between the mitogenome and chloroplast genome of Z. latifolia. Also, approximately 85% of fragments from the mitogenome were duplicated in the Z. latifolia nuclear genome. Selection pressure analysis revealed that most of the mitochondrial genes were highly conserved except for ccmFc, ccmFn, matR, rps1, and rps3. A total of 93 RNA editing sites were found in the PCGs of the mitogenome. Z. latifolia and Oryza minuta are the most closely related, as shown by collinear analysis and the phylogenetic analysis. We found that repeat sequences and foreign sequences in the mitogenomes of Oryzoideae plants were associated with genome rearrangements. In general, the availability of the Z. latifolia mitogenome will contribute valuable information to our understanding of the molecular and genomic aspects of Zizania.

Highly active repeat-mediated recombination in the mitogenome of the aquatic grass Hygroryza aristata

BACKGROUND

Floating bamboo (Hygroryza aristata) is an endangered species with a narrow native distribution and is renowned for its unique aesthetic qualities, which holds significant ecological and ornamental value. However, the lack of genetic information research, with only one complete plastome available, significantly hampers conservation efforts and further research for this species.

RESULTS

In this research, we sequenced and assembled the organelle genomes of floating bamboo, including the mitogenome (587,847 bp) and plastome (135,675 bp). The mitogenome can recombine into various configurations, which are mediated by 25 repeat pairs (13 SRs, 6 MRs, 1 LR, and 5 CRs). LR1 and SR5 are particularly notable as they have the ability to combine with other contigs, forming complex repeat units that facilitate further homologous recombination. The rate of homologous recombination varies significantly among species, yet there is still a pronounced positive correlation observed between the length of these repeat pairs and the rate of recombination they mediate. The mitogenome integrates seven intact protein-coding genes from the chloroplast. The codon usage patterns in both organelles are similar, with a noticeable bias towards C and T on the third codon. The gene map of Poales shows the entire loss of rpl6, succinate dehydrogenase subunits (sdh3 and sdh4). Additionally, the BOP clade retained more variable genes compared to the PACMAD clade.

CONCLUSIONS

We provided a high-quality and well-annotated mitogenome for floating bamboo and demonstrated the presence of diverse configurations. Our study has revealed the correlation between repeat length and their corresponding recombination rate despite variations among species. Although the mitogenome can potentially exist in the form of a unicircular in vivo, this occurrence is rare and may not be stable.

A highly resolved nuclear phylogeny uncovers strong phylogenetic conservatism and correlated evolution of fruit color and size in Solanum L

Mutualisms between plants and fruit-eating animals were key to the radiation of angiosperms. Still, phylogenetic uncertainties limit our understanding of fleshy-fruit evolution, as in the case of Solanum, a genus with remarkable fleshy-fruit diversity, but with unresolved phylogenetic relationships. We used 1786 nuclear genes from 247 species, including 122 newly generated transcriptomes/genomes, to reconstruct the Solanum phylogeny and examine the tempo and mode of the evolution of fruit color and size. Our analysis resolved the backbone phylogeny of Solanum, providing high support for its clades. Our results pushed back the origin of Solanum to 53.1 million years ago (Ma), with most major clades diverging between 35 and 27 Ma. Evolution of Solanum fruit color and size revealed high levels of trait conservatism, where medium-sized berries that remain green when ripe are the likely ancestral form. Our analyses revealed that fruit size and color are evolutionary correlated, where dull-colored fruits are two times larger than black/purple and red fruits. We conclude that the strong phylogenetic conservatism shown in the color and size of Solanum fruits could limit the influences of fruit-eating animals on fleshy-fruit evolution. Our findings highlight the importance of phylogenetic constraints on the diversification of fleshy-fruit functional traits.

Diversity and distribution of bamboo-feeding true bugs in China

The Bambusoideae subfamily, originating in the late Cretaceous, has evolved to include over 1500 species globally. Notably, China hosts the richest diversity of Bambusoideae, with 728 species documented. After a long period of coevolution, plenty of animals could feed on these plants rich in cellulose and lignin. As an important group of pests and participants in the ecosystem, bamboo-feeding true bugs (BFTBs, or bamboo-feeding Heteropteran insects) have attracted the attention of researchers. However, the diversity and distribution of BFTBs still lack systematic and generalized research. In this study, we reviewed the BFTBs in China and simulated the diversity pattern and the driving forces of this pattern. A list of 36 genera with 69 species of BFTBs in China was obtained through paper review and field surveys. And their bamboo-feeding habit had multiple independent origins. The spatial diversity pattern showed that the biodiversity hotspots of BFTBs are located in and around the tropics of southern China. Environmental driving force analysis showed that the minimum temperature of coldest month and annual precipitation were the dominant environmental factors shaping the spatial diversity of BFTBs. Our work quantified the diversity and distribution of BFTBs in China, providing fundamental data support for pest control and evolutionary research.

Progress in systematics and biogeography of Orchidaceae

Orchidaceae are one of the largest families of angiosperms in terms of species richness. In the last decade, numerous studies have delved into reconstructing the phylogenetic framework of Orchidaceae, leveraging data from plastid, mitochondrial and nuclear sources. These studies have provided new insights into the systematics, diversification and biogeography of Orchidaceae, establishing a robust foundation for future research. Nevertheless, pronounced controversies persist regarding the precise placement of certain lineages within these phylogenetic frameworks. To address these discrepancies and deepen our understanding of the phylogenetic structure of Orchidaceae, we provide a comprehensive overview and analysis of phylogenetic studies focusing on contentious groups within Orchidaceae since 2015, delving into discussions on the underlying reasons for observed topological conflicts. We also provide a novel phylogenetic framework at the subtribal level. Furthermore, we examine the tempo and mode underlying orchid species diversity from the perspective of historical biogeography, highlighting factors contributing to extensive speciation. Ultimately, we delineate avenues for future research aimed at enhancing our understanding of Orchidaceae phylogeny and diversity.

Probing of plant transcriptomes reveals the hidden genetic diversity of the family Secoviridae

Secoviruses are single-stranded RNA viruses that infect plants. In the present study, we identified 61 putative novel secoviral genomes in various plant species by mining publicly available plant transcriptome data. These viral sequences represent the genomes of 13 monopartite and 48 bipartite secovirids. The genome sequences of 52 secovirids were coding-complete, and nine were partial. Except for small open reading frames (ORFs) determined in waikaviral genomes and RNA2 of torradoviruses, all of the recovered genomes/genome segments contained a large ORF encoding a polyprotein. Based on genome organization and phylogeny, all but three of the novel secoviruses were assigned to different genera. The genome organization of two identified waika-like viruses resembled that of the recently identified waika-like virus Triticum aestivum secovirus. Phylogenetic analysis revealed a pattern of host-virus co-evolution in a few waika- and waika-like viruses and increased phylogenetic diversity of nepoviruses. The study provides a basis for further investigation of the biological properties of these novel secoviruses.

Guard cell and subsidiary cell sizes are key determinants for stomatal kinetics and drought adaptation in cereal crops

Climate change-induced drought is a major threat to agriculture. C4 crops have a higher water use efficiency (WUE) and better adaptability to drought than C3 crops due to their smaller stomatal morphology and faster response. However, our understanding of stomatal behaviours in both C3 and C4 Poaceae crops is limited by knowledge gaps in physical traits of guard cell (GC) and subsidiary cell (SC). We employed infrared gas exchange analysis and a stomatal assay to explore the relationship between GC/SC sizes and stomatal kinetics across diverse drought conditions in two C3 (wheat and barley) and three C4 (maize, sorghum and foxtail millet) upland Poaceae crops. Through statistical analyses, we proposed a GCSC-τ model to demonstrate how morphological differences affect stomatal kinetics in C4 Poaceae crops. Our findings reveal that morphological variations specifically correlate with stomatal kinetics in C4 Poaceae crops, but not in C3 ones. Subsequent modelling and experimental validation provide further evidence that GC/SC sizes significantly impact stomatal kinetics, which affects stomatal responses to different drought conditions and thereby WUE in C4 Poaceae crops. These findings emphasize the crucial advantage of GC/SC morphological characteristics and stomatal kinetics for the drought adaptability of C4 Poaceae crops, highlighting their potential as future climate-resilient crops.

Habitat differentiation and environmental adaptability contribute to leaf size variations globally in C3 and C4 grasses

The grass family (Poaceae) dominates ~43 % of Earth's land area and contributes 33 % of terrestrial primary productivity that is critical to naturally regulating atmosphere CO2 concentration and global climate change. Currently grasses comprise ~11,780 species and ~50 % of them (~6000 species) utilize C4 photosynthetic pathway. Generally, grass species have smaller leaves under colder and drier environments, but it is unclear whether the primary drivers of leaf size differ between C3 and C4 grasses on a global scale. Here, we analyzed 34 environmental variables, such as latitude, elevation, mean annual temperature, mean annual precipitation, and solar radiation etc., through a comparatively comprehensive database of ~3.0 million occurrence records from 1380 C3 and 978 C4 grass species (2358 species in total). Results from this study confirm that C4 grasses have occupied habitats with lower latitudes and elevations, characterized by warmer, sunnier, drier and less fertile environmental conditions. Grass leaf size correlates positively with mean annual temperature and precipitation as expected. Our results also demonstrate that the mean temperature of the wettest quarter of the year is the primary control for C3 leaf size, whereas C4 leaf size is negatively correlated with the difference between summer and winter temperatures. For C4 grasses, phylogeny exerts a significant effect on leaf size but is less important than environmental factors. Our findings highlight the importance of evolutionarily contrasting variations in leaf size between C3 and C4 grasses for shaping their geographical distribution and habitat suitability at the global scale.

Functional Divergence in Orthologous Transcription Factors: Insights from AtCBF2/3/1 and OsDREB1C

Despite traditional beliefs of orthologous genes maintaining similar functions across species, growing evidence points to their potential for functional divergence. C-repeat binding factors/dehydration-responsive element binding protein 1s (CBFs/DREB1s) are critical in cold acclimation, with their overexpression enhancing stress tolerance but often constraining plant growth. In contrast, a recent study unveiled a distinctive role of rice OsDREB1C in elevating nitrogen use efficiency (NUE), photosynthesis, and grain yield, implying functional divergence within the CBF/DREB1 orthologs across species. Here, we delve into divergent molecular mechanisms of OsDREB1C and AtCBF2/3/1 by exploring their evolutionary trajectories across rice and Arabidopsis genomes, regulatomes, and transcriptomes. Evolutionary scrutiny shows discrete clades for OsDREB1C and AtCBF2/3/1, with the Poaceae-specific DREB1C clade mediated by a transposon event. Genome-wide binding profiles highlight OsDREB1C's preference for GCCGAC compared to AtCBF2/3/1's preference for A/GCCGAC, a distinction determined by R12 in the OsDREB1C AP2/ERF domain. Cross-species multiomic analyses reveal shared gene orthogroups (OGs) and underscore numerous specific OGs uniquely bound and regulated by OsDREB1C, implicated in NUE, photosynthesis, and early flowering, or by AtCBF2/3/1, engaged in hormone and stress responses. This divergence arises from gene gains/losses (∼16.7% to 25.6%) and expression reprogramming (∼62.3% to 66.2%) of OsDREB1C- and AtCBF2/3/1-regulated OGs during the extensive evolution following the rice-Arabidopsis split. Our findings illustrate the regulatory evolution of OsDREB1C and AtCBF2/3/1 at a genomic scale, providing insights on the functional divergence of orthologous transcription factors following gene duplications across species.

Molecular Evolution of RAMOSA1 (RA1) in Land Plants

RAMOSA1 (RA1) is a Cys2-His2-type (C2H2) zinc finger transcription factor that controls plant meristem fate and identity and has played an important role in maize domestication. Despite its importance, the origin of RA1 is unknown, and the evolution in plants is only partially understood. In this paper, we present a well-resolved phylogeny based on 73 amino acid sequences from 48 embryophyte species. The recovered tree topology indicates that, during grass evolution, RA1 arose from two consecutive SUPERMAN duplications, resulting in three distinct grass sequence lineages: RA1-like A, RA1-like B, and RA1; however, most of these copies have unknown functions. Our findings indicate that RA1 and RA1-like play roles in the nucleus despite lacking a traditional nuclear localization signal. Here, we report that copies diversified their coding region and, with it, their protein structure, suggesting different patterns of DNA binding and protein-protein interaction. In addition, each of the retained copies diversified regulatory elements along their promoter regions, indicating differences in their upstream regulation. Taken together, the evidence indicates that the RA1 and RA1-like gene families in grasses underwent subfunctionalization and neofunctionalization enabled by gene duplication.

Are cereal grasses a single genetic system?

In 1993, a passionate and provocative call to arms urged cereal researchers to consider the taxon they study as a single genetic system and collaborate with each other. Since then, that group of scientists has seen their discipline blossom. In an attempt to understand what unity of genetic systems means and how the notion was borne out by later research, we survey the progress and prospects of cereal genomics: sequence assemblies, population-scale sequencing, resistance gene cloning and domestication genetics. Gene order may not be as extraordinarily well conserved in the grasses as once thought. Still, several recurring themes have emerged. The same ancestral molecular pathways defining plant architecture have been co-opted in the evolution of different cereal crops. Such genetic convergence as much as cross-fertilization of ideas between cereal geneticists has led to a rich harvest of genes that, it is hoped, will lead to improved varieties.

Unprecedented variation pattern of plastid genomes and the potential role in adaptive evolution in Poales

BACKGROUND

The plastid is the photosynthetic organelle in plant cell, and the plastid genomes (plastomes) are generally conserved in evolution. As one of the most economically and ecologically important order of angiosperms, Poales was previously documented to exhibit great plastomic variation as an order of photoautotrophic plants.

RESULTS

We acquired 93 plastomes, representing all the 16 families and 5 major clades of Poales to reveal the extent of their variation and evolutionary pattern. Extensive variation including the largest one in monocots with 225,293 bp in size, heterogeneous GC content, and a wide variety of gene duplication and loss were revealed. Moreover, rare occurrences of three inverted repeat (IR) copies in angiosperms and one IR loss were observed, accompanied by short IR (sIR) and small direct repeat (DR). Widespread structural heteroplasmy, diversified inversions, and unusual genomic rearrangements all appeared in Poales, occasionally within a single species. Extensive repeats in the plastomes were found to be positively correlated with the observed inversions and rearrangements. The variation all showed a "small-large-moderate" trend along the evolution of Poales, as well as for the sequence substitution rate. Finally, we found some positively selected genes, mainly in C4 lineages, while the closely related lineages of those experiencing gene loss tended to have undergone more relaxed purifying selection.

CONCLUSIONS

The variation of plastomes in Poales may be related to its successful diversification into diverse habitats and multiple photosynthetic pathway transitions. Our order-scale analyses revealed unusual evolutionary scenarios for plastomes in the photoautotrophic order of Poales and provided new insights into the plastome evolution in angiosperms as a whole.

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.

Eleven new species of Trichoderma (Hypocreaceae, Hypocreales) from China

Trichoderma spp. are globally distributed and are considered significant fungal resources. This study presents the discovery of 11 new species of Trichoderma: T. caeruleum, T. gongcheniae, T. graminicola, T. graminis, T. hongkuii, T. parapeberdyi, T. neoguizhouense, T. neohongkuii, T. parahamatum, T. parahongkuii, and T. shaanxiensis. All of these new species were isolated from soils, except for T. caeruleum, T. graminicola, T. graminis, and T. neohongkuii, which were found as endophytes in Poaceae plants. The phylogenetic position of these novel species was determined by analysing the concatenated sequences of the second largest nuclear RNA polymerase subunit encoding gene (rpb2) and the translation elongation factor 1-alpha encoding gene (tef1). The results of the phylogenetic analysis revealed that each new species formed a distinct lineage: T. gongcheniae, T. graminicola, T. graminis, T. neoguizhouense, T. parapeberdyi, and T. shaanxiensis belong to the Harzianum Clade, T. hongkuii, T. parahongkuii, and T. neohongkuii are new members of the Koningii Clade, T. parahamatum belongs to the Hamatum Clade, and T. caeruleum does not fall within any of the named clades. The study also provided a detailed description of the morphology and cultural characteristics of the newly discovered Trichoderma species. The discovery contributes to the advancement of knowledge about Trichoderma species resources in China.

The newly assembled chloroplast genome of Aeluropus littoralis: molecular feature characterization and phylogenetic analysis with related species

Aeluropus littoralis, a halophyte grass, is widely distributed from the Mediterranean to the Indian subcontinent through the Mongolian Gobi. This model halophyte has garnered increasing attention owing to its use as forage and its high tolerance to environmental stressors. The chloroplast genomes of many plants have been extensively examined for molecular, phylogenetic and transplastomic applications. However, no published research on the A. littoralis chloroplast (cp) genome was discovered. Here, the entire chloroplast genome of A. littoralis was assembled implementing accurate long-read sequences. The entire chloroplast genome, with an estimated length of 135,532 bp (GC content: 38.2%), has a quadripartite architecture and includes a pair of inverted repeat (IR) regions, IRa and IRb (21,012 bp each), separated by a large and a small single-copy regions (80,823 and 12,685 bp, respectively). The features of A. littoralis consist of 133 genes that synthesize 87 peptides, 38 transfer RNAs, and 8 ribosomal RNAs. Of these genes, 86 were unique, whereas 19 were duplicated in IR regions. Additionally, a total of forty-six simple sequence repeats, categorized into 32-mono, four-di, two-tri, and eight-tetranucleotides, were discovered. Furthermore, ten sets of repeats greater than 20 bp were located primarily in the LSC region. Evolutionary analysis based on chloroplast sequence data revealed that A. littoralis with A. lagopoides and A. sinensis belong to the Aeluropodinae subtribe, which is a sister to the Eleusininae in the tribe Cynodonteae and the subfamily Chloridoideae. This subfamily belongs to the PACMAD clade, which contains the majority of the C4 photosynthetic plants in the Poaceae. The newly constructed A. littoralis cp genome offers valuable knowledge for DNA barcoding, phylogenetic, transplastomic research, and other biological studies.

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

Angiosperms (flowering plants) are by far the most diverse land plant group with over 300,000 species. The sudden appearance of diverse angiosperms in the fossil record was referred to by Darwin as the "abominable mystery," hence contributing to the heightened interest in angiosperm evolution. Angiosperms display wide ranges of morphological, physiological, and ecological characters, some of which have probably influenced their species richness. The evolutionary analyses of these characteristics help to address questions of angiosperm diversification and require well resolved phylogeny. Following the great successes of phylogenetic analyses using plastid sequences, dozens to thousands of nuclear genes from next-generation sequencing have been used in angiosperm phylogenomic analyses, providing well resolved phylogenies and new insights into the evolution of angiosperms. In this review we focus on recent nuclear phylogenomic analyses of large angiosperm clades, orders, families, and subdivisions of some families and provide a summarized Nuclear Phylogenetic Tree of Angiosperm Families. The newly established nuclear phylogenetic relationships are highlighted and compared with previous phylogenetic results. The sequenced genomes of Amborella, Nymphaea, Chloranthus, Ceratophyllum, and species of monocots, Magnoliids, and basal eudicots, have facilitated the phylogenomics of relationships among five major angiosperms clades. All but one of the 64 angiosperm orders were included in nuclear phylogenomics with well resolved relationships except the placements of several orders. Most families have been included with robust and highly supported placements, especially for relationships within several large and important orders and families. Additionally, we examine the divergence time estimation and biogeographic analyses of angiosperm on the basis of the nuclear phylogenomic frameworks and discuss the differences compared with previous analyses. Furthermore, we discuss the implications of nuclear phylogenomic analyses on ancestral reconstruction of morphological, physiological, and ecological characters of angiosperm groups, limitations of current nuclear phylogenomic studies, and the taxa that require future attention.

Identification and Analysis of PEPC Gene Family Reveals Functional Diversification in Orchidaceae and the Regulation of Bacterial-Type PEPC

Phosphoenolpyruvate carboxylase (PEPC) gene family plays a crucial role in both plant growth and response to abiotic stress. Approximately half of the Orchidaceae species are estimated to perform CAM pathway, and the availability of sequenced orchid genomes makes them ideal subjects for investigating the PEPC gene family in CAM plants. In this study, a total of 33 PEPC genes were identified across 15 orchids. Specifically, one PEPC gene was found in Cymbidium goeringii and Platanthera guangdongensis; two in Apostasia shenzhenica, Dendrobium chrysotoxum, D. huoshanense, Gastrodia elata, G. menghaiensis, Phalaenopsis aphrodite, Ph. equestris, and Pl. zijinensis; three in C. ensifolium, C. sinense, D. catenatum, D. nobile, and Vanilla planifolia. These PEPC genes were categorized into four subgroups, namely PEPC-i, PEPC-ii, and PEPC-iii (PTPC), and PEPC-iv (BTPC), supported by the comprehensive analyses of their physicochemical properties, motif, and gene structures. Remarkably, PEPC-iv contained a heretofore unreported orchid PEPC gene, identified as VpPEPC4. Differences in the number of PEPC homolog genes among these species were attributed to segmental duplication, whole-genome duplication (WGD), or gene loss events. Cis-elements identified in promoter regions were predominantly associated with light responsiveness, and circadian-related elements were observed in each PEPC-i and PEPC-ii gene. The expression levels of recruited BTPC, VpPEPC4, exhibited a lower expression level than other VpPEPCs in the tested tissues. The expression analyses and RT-qPCR results revealed diverse expression patterns in orchid PEPC genes. Duplicated genes exhibited distinct expression patterns, suggesting functional divergence. This study offered a comprehensive analysis to unveil the evolution and function of PEPC genes in Orchidaceae.

Angiosperm-wide analysis of fruit and ovary evolution aided by a new nuclear phylogeny supports association of the same ovary type with both dry and fleshy fruits

Fruit functions in seed protection and dispersal and belongs to many dry and fleshy types, yet their evolutionary pattern remains unclear in part due to uncertainties in the phylogenetic relationships among several orders and families. Thus we used nuclear genes of 502 angiosperm species representing 231 families to reconstruct a well supported phylogeny, with resolved relationships for orders and families with previously uncertain placements. Using this phylogeny as a framework, molecular dating supports a Triassic origin of the crown angiosperms, followed by the emergence of most orders in the Jurassic and Cretaceous and their rise to ecological dominance during the Cretaceous Terrestrial Revolution. The robust phylogeny allowed an examination of the evolutionary pattern of fruit and ovary types, revealing a trend of parallel carpel fusions during early diversifications in eudicots, monocots, and magnoliids. Moreover, taxa in the same order or family with the same ovary type can develop either dry or fleshy fruits with strong correlations between specific types of dry and fleshy fruits; such associations of ovary, dry and fleshy fruits define several ovary-fruit "modules" each found in multiple families. One of the frequent modules has an ovary containing multiple ovules, capsules and berries, and another with an ovary having one or two ovules, achenes (or other single-seeded dry fruits) and drupes. This new perspective of relationships among fruit types highlights the closeness of specific dry and fleshy fruit types, such as capsule and berry, that develop from the same ovary type and belong to the same module relative to dry and fleshy fruits of other modules (such as achenes and drupes). Further analyses of gene families containing known genes for ovary and fruit development identified phylogenetic nodes with multiple gene duplications, supporting a possible role of whole-genome duplications, in combination with climate changes and animal behaviors, in angiosperm fruit and ovary diversification.

Invasive alien plants are phylogenetically distinct from other alien species across spatial and taxonomic scales in China

Introduction

Phylogenetic relatedness is one of the important factors in the community assembly process. Here, we aimed to understand the large-scale phylogenetic relationship between alien plant species at different stages of the invasion process and how these relationships change in response to the environmental filtering process at multiple spatial scales and different phylogenetic extents.

Methods

We identified the alien species in three invasion stages, namely invasive, naturalized, and introduced, in China. The occurrence records of the species were used to quantify two abundance-based phylogenetic metrics [the net relatedness index (NRI) and the nearest taxon index (NTI)] from a highly resolved phylogenetic tree. The metrics were compared between the three categories of alien species. Generalized linear models were used to test the effect of climate on the phylogenetic pattern. All analyses were conducted at four spatial scales and for three major angiosperm families.

Results

We observed significantly higher NRI and NTI values at finer spatial scales, indicating the formation of more clustered assemblages of phylogenetically closely related species in response to the environmental filtering process. Positive NTI values for the invasive and naturalized aliens suggested that the presence of a close relative in the community may help the successful naturalization and invasion of the introduced alien species. In the two-dimensional phylogenetic space, the invasive species communities significantly differed from the naturalized and introduced species, indicating that established alien species need to be phylogenetically different to become invasive. Positive phylogenetic measures for the invasive aliens across the spatial scales suggested that the presence of invasive aliens could facilitate the establishment of other invasive species. Phylogenetic relatedness was more influenced by temperature than precipitation, especially at a finer spatial scale. With decreased temperature, the invasive species showed a more clustered assemblage, indicating conservatism of their phylogenetic niche. The phylogenetic pattern was different at the family level, although there was a consistent tendency across families to form more clustered assemblages.

Discussion

Overall, our study showed that the community assemblage became more clustered with the progression of the invasion process. The phylogenetic measures varied at spatial and taxonomic scales, thereby highlighting the importance of assessing phylogenetic patterns at different gradients of the community assembly process.

A Pilot Study to Identify Grass Species That Mediate Pollen Allergy in Thailand

INTRODUCTION

Some grass species that are either common or widely spread in Thailand have not been used for pollen allergy diagnosis. In order to improve diagnostic accuracy, the aim of this pilot study was to identify the grass species responsible for pollen allergy in Thailand.

METHODS

The skin sensitization of pollen extracts from six different grass species, including rice (Oryza sativa), corn (Zea mays), sorghum (Sorghum bicolor), para grass (Urochloa mutica), ruzi grass (Urochloa eminii), and green panic grass (Megathyrsus maximus), was evaluated by skin prick test (SPT). Serum's IgE specific to each pollen extract was analyzed by Western blot (WB). The ImmunoCAPTM test for Johnson grass was also evaluated.

RESULTS

Of the thirty-six volunteers who participated in this study, eighteen tested positive for at least one of the diagnostic tests, namely SPT, WB analysis, or ImmunoCAPTM. Notably, skin reactivity to para grass, corn, sorghum, and rice was more commonly observed compared to ruzi grass and green panic grass. However, in the WB analysis, individuals with pollen-specific IgE were more frequently detected in sorghum, green panic grass, corn, rice, and ruzi grass than para grass.

CONCLUSION

In this pilot investigation, our findings indicate that the pollen extracts of rice, corn, sorghum, and para grass are associated with pollen allergy in Thailand. These results contribute to the current knowledge on the identification of grass species that are associated with pollen allergy in Thailand and Southeast Asia.

A Machine Learning Framework Identifies Plastid-Encoded Proteins Harboring C3 and C4 Distinguishing Sequence Information

C4 photosynthesis is known to have at least 61 independent origins across plant lineages making it one of the most notable examples of convergent evolution. Of the >60 independent origins, a predicted 22-24 origins, encompassing greater than 50% of all known C4 species, exist within the Panicoideae, Arundinoideae, Chloridoideae, Micrairoideae, Aristidoideae, and Danthonioideae (PACMAD) clade of the Poaceae family. This clade is therefore primed with species ideal for the study of genomic changes associated with the acquisition of the C4 photosynthetic trait. In this study, we take advantage of the growing availability of sequenced plastid genomes and employ a machine learning (ML) approach to screen for plastid genes harboring C3 and C4 distinguishing information in PACMAD species. We demonstrate that certain plastid-encoded protein sequences possess distinguishing and informative sequence information that allows them to train accurate ML C3/C4 classification models. Our RbcL-trained model, for example, informs a C3/C4 classifier with greater than 99% accuracy. Accurate prediction of photosynthetic type from individual sequences suggests biologically relevant, and potentially differing roles of these sequence products in C3 versus C4 metabolism. With this ML framework, we have identified several key sequences and sites that are most predictive of C3/C4 status, including RbcL, subunits of the NAD(P)H dehydrogenase complex, and specific residues within, further highlighting their potential significance in the evolution and/or maintenance of C4 photosynthetic machinery. This general approach can be applied to uncover intricate associations between other similar genotype-phenotype relationships.

Phylogenomic analyses of Camellia support reticulate evolution among major clades

Camellia (Theaceae) is a morphologically highly diverse genus of flowering plants and includes many famous species with high economic value, and the phylogeny of this genus is not fully resolved. We used 95 transcriptomes from 87 Camellia species and identified 1481 low-copy genes to conduct a detailed analysis of the phylogeny of this genus according to various data-screening criteria. The results show that, very different from the two existing classification systems of Camellia, 87 species are grouped into 8 main clades and two independent species, and that all 8 clades except Clade 8 were strongly supported by almost all the coalescent or concatenated trees using different gene subsets. However, the relationships among these clades were weakly supported and different from analyses using different gene subsets; furthermore, they do not agree with the phylogeny from chloroplast genomes of Camellia. Additional analyses support reticulate evolution (probably resulting from introgression or hybridization) among some major Camellia lineages, providing explanation for extensive gene tree conflicts. Furthermore, we inferred that together with the formation of East Asian subtropical evergreen broad-leaved forests, Camellia underwent a radiative divergence of major clades at 23 ∼ 19 Ma in the late Miocene then had a subsequent species burst at 10 ∼ 5 Ma. Principal component and cluster analyses provides new insights into morphological changes underlying the evolution of Camellia and a reference to further clarify subgenus and sections of this genus. The comprehensive study here including a nuclear phylogeny and other analyses reveal the rapid evolutionary history of Camellia.

Comprehensive phylogenetic analyses of Orchidaceae using nuclear genes and evolutionary insights into epiphytism

Orchidaceae (with >28,000 orchid species) are one of the two largest plant families, with economically and ecologically important species, and occupy global and diverse niches with primary distribution in rainforests. Among orchids, 70% grow on other plants as epiphytes; epiphytes contribute up to ~50% of the plant diversity in rainforests and provide food and shelter for diverse animals and microbes, thereby contributing to the health of these ecosystems. Orchids account for over two-thirds of vascular epiphytes and provide an excellent model for studying evolution of epiphytism. Extensive phylogenetic studies of Orchidaceae and subgroups have ;been crucial for understanding relationships among many orchid lineages, although some uncertainties remain. For example, in the largest subfamily Epidendroideae with nearly all epiphytic orchids, relationships among some tribes and many subtribes are still controversial, hampering evolutionary analyses of epiphytism. Here we obtained 1,450 low-copy nuclear genes from 610 orchid species, including 431 with newly generated transcriptomes, and used them for the reconstruction of robust Orchidaceae phylogenetic trees with highly supported placements of tribes and subtribes. We also provide generally well-supported phylogenetic placements of 131 genera and 437 species that were not sampled by previous plastid and nuclear phylogenomic studies. Molecular clock analyses estimated the Orchidaceae origin at ~132 million years ago (Ma) and divergences of most subtribes from 52 to 29 Ma. Character reconstruction supports at least 14 parallel origins of epiphytism; one such origin was placed at the most recent common ancestor of ~95% of epiphytic orchids and linked to modern rainforests. Ten occurrences of rapid increase in the diversification rate were detected within Epidendroideae near and after the K-Pg boundary, contributing to ~80% of the Orchidaceae diversity. This study provides a robust and the largest family-wide Orchidaceae nuclear phylogenetic tree thus far and new insights into the evolution of epiphytism in vascular plants.

The mitochondrial genome of the diploid oat Avena longiglumis

BACKGROUND

Avena longiglumis Durieu (2n = 2x = 14) is a wild relative of cultivated oat (Avena sativa, 2n = 6x = 42) with good agronomic and nutritional traits. The plant mitochondrial genome has a complex organization and carries genetic traits of value in exploiting genetic resources, not least male sterility alleles used to generate F₁ hybrid seeds. Therefore, we aim to complement the chromosomal-level nuclear and chloroplast genome assemblies of A. longiglumis with the complete assembly of the mitochondrial genome (mitogenome) based on Illumina and ONT long reads, comparing its structure with Poaceae species.

RESULTS

The complete mitochondrial genome of A. longiglumis can be represented by one master circular genome being 548,445 bp long with a GC content of 44.05%. It can be represented by linear or circular DNA molecules (isoforms or contigs), with multiple alternative configurations mediated by long (4,100-31,235 bp) and medium (144-792 bp) size repeats. Thirty-five unique protein-coding genes, three unique rRNA genes, and 11 unique tRNA genes are identified. The mitogenome is rich in duplications (up to 233 kb long) and multiple tandem or simple sequence repeats, together accounting for more than 42.5% of the total length. We identify homologous sequences between the mitochondrial, plastid and nuclear genomes, including the exchange of eight plastid-derived tRNA genes, and nuclear-derived retroelement fragments. At least 85% of the mitogenome is duplicated in the A. longiglumis nuclear genome. We identify 269 RNA editing sites in mitochondrial protein-coding genes including stop codons truncating ccmFC transcripts.

CONCLUSIONS

Comparative analysis with Poaceae species reveals the dynamic and ongoing evolutionary changes in mitochondrial genome structure and gene content. The complete mitochondrial genome of A. longiglumis completes the last link of the oat reference genome and lays the foundation for oat breeding and exploiting the biodiversity in the genus.

Nuclear phylogeny and insights into whole-genome duplications and reproductive development of Solanaceae plants

Solanaceae, the nightshade family, have ∼2700 species, including the important crops potato and tomato, ornamentals, and medicinal plants. Several sequenced Solanaceae genomes show evidence for whole-genome duplication (WGD), providing an excellent opportunity to investigate WGD and its impacts. Here, we generated 93 transcriptomes/genomes and combined them with 87 public datasets, for a total of 180 Solanaceae species representing all four subfamilies and 14 of 15 tribes. Nearly 1700 nuclear genes from these transcriptomic/genomic datasets were used to reconstruct a highly resolved Solanaceae phylogenetic tree with six major clades. The Solanaceae tree supports four previously recognized subfamilies (Goetzeioideae, Cestroideae, Nicotianoideae, and Solanoideae) and the designation of three other subfamilies (Schizanthoideae, Schwenckioideae, and Petunioideae), with the placement of several previously unassigned genera. We placed a Solanaceae-specific whole-genome triplication (WGT1) at ∼81 million years ago (mya), before the divergence of Schizanthoideae from other Solanaceae subfamilies at ∼73 mya. In addition, we detected two gene duplication bursts (GDBs) supporting proposed WGD events and four other GDBs. An investigation of the evolutionary histories of homologs of carpel and fruit developmental genes in 14 gene (sub)families revealed that 21 gene clades have retained gene duplicates. These were likely generated by the Solanaceae WGT1 and may have promoted fleshy fruit development. This study presents a well-resolved Solanaceae phylogeny and a new perspective on retained gene duplicates and carpel/fruit development, providing an improved understanding of Solanaceae evolution.

A chromosome-scale genome sequence of sudangrass (Sorghum sudanense) highlights the genome evolution and regulation of dhurrin biosynthesis

Sudangrass is more similar to US commercial sorghums than to cultivated sorghums from Africa sequence-wise and contain significantly lower dhurrin than sorghums. CYP79A1 is linked to dhurrin content in sorghum. Sudangrass [Sorghum sudanense (Piper) Stapf] is a hybrid between grain sorghum and its wild relative S. bicolor ssp. verticilliflorum and is grown as a forage crop due to its high biomass production and low dhurrin content compared to sorghum. In this study, we sequenced the sudangrass genome and showed that the assembled genome was 715.95 Mb with 35,243 protein-coding genes. Phylogenetic analysis with whole genome proteomes demonstrated that the sudangrass genome was more similar to US commercial sorghums than to its wild relatives and cultivated sorghums from Africa. We confirmed that at seedling stage, sudangrass accessions contained significantly lower dhurrin as measured by hydrocyanic acid potential (HCN-p) than cultivated sorghum accessions. Genome-wide association study identified a QTL most tightly associated with HCN-p and the linked SNPs were located in the 3' UTR of Sobic.001G012300 which encodes CYP79A1, the enzyme that catalyzes the first step of dhurrin biosynthesis. As in other grasses such as maize and rice, we also found that copia/gypsy long terminal repeat (LTR) retrotransposons were more abundant in cultivated than in wild sorghums, implying that crop domestication in the grasses was accompanied by increased copia/gypsy LTR retrotransposon insertions in the genomes.

A plastid phylogenomic framework for the palm family (Arecaceae)

BACKGROUND

Over the past decade, phylogenomics has greatly advanced our knowledge of angiosperm evolution. However, phylogenomic studies of large angiosperm families with complete species or genus-level sampling are still lacking. The palms, Arecaceae, are a large family with ca. 181 genera and 2600 species and are important components of tropical rainforests bearing great cultural and economic significance. Taxonomy and phylogeny of the family have been extensively investigated by a series of molecular phylogenetic studies in the last two decades. Nevertheless, some phylogenetic relationships within the family are not yet well-resolved, especially at the tribal and generic levels, with consequent impacts for downstream research.

RESULTS

Plastomes of 182 palm species representing 111 genera were newly sequenced. Combining these with previously published plastid DNA data, we were able to sample 98% of palm genera and conduct a plastid phylogenomic investigation of the family. Maximum likelihood analyses yielded a robustly supported phylogenetic hypothesis. Phylogenetic relationships among all five palm subfamilies and 28 tribes were well-resolved, and most inter-generic phylogenetic relationships were also resolved with strong support.

CONCLUSIONS

The inclusion of nearly complete generic-level sampling coupled with nearly complete plastid genomes strengthened our understanding of plastid-based relationships of the palms. This comprehensive plastid genome dataset complements a growing body of nuclear genomic data. Together, these datasets form a novel phylogenomic baseline for the palms and an increasingly robust framework for future comparative biological studies of this exceptionally important plant family.

Poaceae Chloroplast Genome Sequencing: Great Leap Forward in Recent Ten Years

The first complete chloroplast genome of rice (Oryza sativa) was published in 1989, ushering in a new era of studies of chloroplast genomics in Poaceae. Progresses in Next-Generation Sequencing (NGS) and Third-Generation Sequencing (TGS) technologiesand in the development of genome assembly software, have significantly advanced chloroplast genomics research. Poaceae is one of the most targeted families in chloroplast genome research because of its agricultural, ecological, and economic importance. Over the last 30 years, 2,050 complete chloroplast genome sequences from 40 tribes and 282 genera have been generated, most (97%) of them in the recent ten years. The wealth of data provides the groundwork for studies on species evolution, phylogeny, genetic transformation, and other aspects of Poaceae chloroplast genomes. As a result, we have gained a deeper understanding of the properties of Poaceae chloroplast genomes. Here, we summarize the achievements of the studies of the Poaceae chloroplast genomes and envision the challenges for moving the area ahead.

Improving abiotic stress tolerance of forage grasses - prospects of using genome editing

Due to an increase in the consumption of food, feed, and fuel and to meet global food security needs for the rapidly growing human population, there is a necessity to obtain high-yielding crops that can adapt to future climate changes. Currently, the main feed source used for ruminant livestock production is forage grasses. In temperate climate zones, perennial grasses grown for feed are widely distributed and tend to suffer under unfavorable environmental conditions. Genome editing has been shown to be an effective tool for the development of abiotic stress-resistant plants. The highly versatile CRISPR-Cas system enables increasingly complex modifications in genomes while maintaining precision and low off-target frequency mutations. In this review, we provide an overview of forage grass species that have been subjected to genome editing. We offer a perspective view on the generation of plants resilient to abiotic stresses. Due to the broad factors contributing to these stresses the review focuses on drought, salt, heat, and cold stresses. The application of new genomic techniques (e.g., CRISPR-Cas) allows addressing several challenges caused by climate change and abiotic stresses for developing forage grass cultivars with improved adaptation to the future climatic conditions. Genome editing will contribute towards developing safe and sustainable food systems.

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.

Independent recruitment of FRUITFULL-like transcription factors in the convergent origins of vernalization-responsive grass flowering

Flowering in response to low temperatures (vernalization) has evolved multiple times independently across angiosperms as an adaptation to match reproductive development with the short growing season of temperate habitats. Despite the context of a generally conserved flowering time network, evidence suggests that the genes underlying vernalization responsiveness are distinct across major plant clades. Whether different or similar mechanisms underlie vernalization-induced flowering at narrower (e.g., family-level) phylogenetic scales is not well understood. To test the hypothesis that vernalization responsiveness has evolved convergently in temperate species of the grass family (Poaceae), we carried out flowering time experiments with and without vernalization in several representative species from different subfamilies. We then determined the likelihood that vernalization responsiveness evolved through parallel mechanisms by quantifying the response of Pooideae vernalization pathway FRUITFULL (FUL)-like genes to extended periods of cold. Our results demonstrate that vernalization-induced flowering has evolved multiple times independently in at least five grass subfamilies, and that different combinations of FUL-like genes have been recruited to this pathway on several occasions.