Domestication has shaped the population structure and agronomic traits of tea plants, yet the complexity of tea population structure and genetic variation that determines these traits remains unclear. We here investigated the resequencing data of 363 diverse tea accessions collected extensively from almost all tea distributions and found that the population structure of tea plants was divided into eight subgroups, which were basically consistent with their geographical distributions. The genetic diversity of tea plants in China decreased from southwest to east as latitude increased. Results also indicated that Camellia sinensis var. assamica (CSA) illustrated divergent selection signatures with Camellia sinensis var. sinensis (CSS). The domesticated genes of CSA were mainly involved in leaf development, flavonoid and alkaloid biosynthesis, while the domesticated genes in CSS mainly participated in amino acid metabolism, aroma compounds biosynthesis, and cold stress. Comparative population genomics further identified ~730 Mb novel sequences, generating 6,058 full-length protein-encoding genes, significantly expanding the gene pool of tea plants. We also discovered 217,376 large-scale structural variations and 56,583 presence and absence variations (PAVs) across diverse tea accessions, some of which were associated with tea quality and stress resistance. Functional experiments demonstrated that two PAV genes (CSS0049975 and CSS0006599) were likely to drive trait diversification in cold tolerance between CSA and CSS tea plants. The overall findings not only revealed the genetic diversity and domestication of tea plants, but also underscored the vital role of structural variations in the diversification of tea plant traits.

Gene innovation plays an essential role in trait evolution. Rhizobial symbioses, the most important N₂-fixing agent in agricultural systems that exists mainly in Leguminosae, is one of the most attractive evolution events. However, the gene innovations underlying Leguminosae root nodule symbiosis (RNS) remain largely unknown. Here, we investigated the gene gain event in Leguminosae RNS evolution through comprehensive phylogenomic analyses. We revealed that Leguminosae-gain genes were acquired by gene duplication and underwent a strong purifying selection. Kyoto Encyclopedia of Genes and Genomes analyses showed that the innovated genes were enriched in flavonoid biosynthesis pathways, particular downstream of chalcone synthase (CHS). Among them, Leguminosae-gain type Ⅱ chalcone isomerase (CHI) could be further divided into CHI1A and CHI1B clades, which resulted from the products of tandem duplication. Furthermore, the duplicated CHI genes exhibited exon–intron structural divergences evolved through exon/intron gain/loss and insertion/deletion. Knocking down CHI1B significantly reduced nodulation in Glycine max (soybean) and Medicago truncatula; whereas, knocking down its duplication gene CHI1A had no effect on nodulation. Therefore, Leguminosae-gain type Ⅱ CHI participated in RNS and the duplicated CHI1A and CHI1B genes exhibited RNS functional divergence. This study provides functional insights into Leguminosae-gain genetic innovation and sub-functionalization after gene duplication that contribute to the evolution and adaptation of RNS in Leguminosae.

Pineapple is the third most crucial tropical fruit worldwide and available in five varieties. Genomes of different pineapple varieties have been released to date; however, none of them are complete, with all exhibiting substantial gaps and representing only two of the five pineapple varieties. This significantly hinders the advancement of pineapple breeding efforts. In this study, we sequenced the genomes of three varieties: a wild pineapple variety, a fiber pineapple variety, and a globally cultivated edible pineapple variety. We constructed the first gap-free reference genome (Ref) for pineapple. By consolidating multiple sources of evidence and manually revising each gene structure annotation, we identified 26,656 protein-coding genes. The BUSCO evaluation indicated a completeness of 99.2%, demonstrating the high quality of the gene structure annotations in this genome. Utilizing these resources, we identified 7,209 structural variations across the three varieties. Approximately 30.8% of pineapple genes were located within ±5 kb of structural variations, including 30 genes associated with anthocyanin synthesis. Further analysis and functional experiments demonstrated that the high expression of AcMYB528 aligns with the accumulation of anthocyanins in the leaves, both of which may be affected by a 1.9-kb insertion fragment. In addition, we developed the Ananas Genome Database, which offers data browsing, retrieval, analysis, and download functions. The construction of this database addresses the lack of pineapple genome resource databases. In summary, we acquired a seamless pineapple reference genome with high-quality gene structure annotations, providing a solid foundation for pineapple genomics and a valuable reference for pineapple breeding.

The structural and functional diversity of plant metabolites is largely created via chemical modification of a basic backbone. However, metabolite modifications in plants have still not been thoroughly investigated by metabolomics approaches. In this study, a widely targeted metabolite modificomics (WTMM) strategy was developed based on ultra-high performance liquid chromatography-quadrupole-linear ion trap (UHPLC-Q-Trap) and UHPLC-Q-Exactive-Orbitrap (UHPLC-QE-Orbitrap), which greatly improved the detection sensitivity and the efficiency of identification of modified metabolites. A metabolite modificomics study was carried out using tomato as a model, and over 34,000 signals with MS2 information were obtained from approximately 232 neutral loss transitions. Unbiased metabolite profiling was also performed by utilizing high-resolution mass spectrometry data to annotate a total of 2,118 metabolites with 125 modification types; of these, 165 modified metabolites were identified in this study. Next, the WTMM database was used to assess diseased tomato tissues and 29 biomarkers were analyzed. In summary, the WTMM strategy is not only capable of large-scale detection and quantitative analysis of plant-modified metabolites in plants, but also can be used for plant biomarker development.

The genus Corydalis, with ca. 530 species, has long been considered taxonomically challenging because of its great variability. Previous molecular analyses, based on a few molecular markers and incomplete taxonomic sampling, were clearly inadequate to delimit sections and subgenera. We have performed phylogenetic analyses of Corydalis and related taxa, using 65 shared protein-coding plastid genes from 313 accessions (including 280 samples of ca. 226 species of Corydalis) and 152 universal low-copy nuclear genes from 296 accessions (including 271 samples of Corydalis) covering all 42 previously recognized sections and five independent “series”. Phylogenetic trees were inferred using Bayesian Inference and Maximum Likelihood. Eight selected morphological characters were estimated using ancestral state reconstructions. Results include: (i) of the three subgenera of Corydalis, two are fully supported by both the plastid and nuclear data; the third, subg. Cremnocapnos, is weakly supported by plastid DNA only, whereas in the nuclear data the two included sections form successive outgroups to the rest of the genus; (ii) among all 42 sections and five “series”, 25 sections and one “series” are resolved as monophyletic in both data sets; (iii) the common ancestor of Corydalis is likely to be a perennial plant with a taproot, yellow flowers with a short saccate spur, linear fruits with recurved fruiting pedicels, and seeds with elaiosomes; (iv) we provide a new classification of Corydalis with four subgenera (of which subg. Bipapillatae is here newly described), 39 sections, 16 of which are consistent with the previous classification, 16 sections have been recircumscribed, one section has been reinstated and six new sections are established. Characters associated with lifespan, underground structures, floral spur, fruit and elaiosomes are important for the recognition of subgenera and sections. These new phylogenetic analyses combined with ancestral character reconstructions uncovered previously unrecognized relationships, and greatly improved our understanding of the evolution of the genus.

Evolutionary convergence is one of the most striking examples of adaptation driven by natural selection. However, genomic evidence for convergent adaptation to extreme environments remains scarce. Here, we assembled reference genomes of two alpine plants, Saussurea obvallata (Asteraceae) and Rheum alexandrae (Polygonaceae), with 37,938 and 61,463 annotated protein-coding genes. By integrating an additional five alpine genomes, we elucidated genomic convergence underlying high-altitude adaptation in alpine plants. Our results detected convergent contractions of disease-resistance genes in alpine genomes, which might be an energy-saving strategy for surviving in hostile environments with only a few pathogens present. We identified signatures of positive selection on a set of genes involved in reproduction and respiration (e.g., MMD1, NBS1, and HPR), and revealed signatures of molecular convergence on genes involved in self-incompatibility, cell wall modification, DNA repair and stress resistance, which may underlie adaptation to extreme cold, high ultraviolet radiation and hypoxia environments. Incorporating transcriptomic data, we further demonstrated that genes associated with cuticular wax and flavonoid biosynthetic pathways exhibit higher expression levels in leafy bracts, shedding light on the genetic mechanisms of the adaptive “greenhouse” morphology. Our integrative data provide novel insights into convergent evolution at a high-taxonomic level, aiding in a deep understanding of genetic adaptation to complex environments.

Angiosperms dominate the Earth's ecosystems and provide most of the basic necessities for human life. The major angiosperm clades comprise 64 orders, as recognized by the APG IV classification. However, the phylogenetic relationships of angiosperms remain unclear, as phylogenetic trees with different topologies have been reconstructed depending on the sequence datasets utilized, from targeted genes to transcriptomes. Here, we used currently available de novo genome data to reconstruct the phylogenies of 366 angiosperm species from 241 genera belonging to 97 families across 43 of the 64 orders based on orthologous genes from the nuclear, plastid, and mitochondrial genomes of the same species with compatible datasets. The phylogenetic relationships were largely consistent with previously constructed phylogenies based on sequence variations in each genome type. However, there were major inconsistencies in the phylogenetic relationships of the five Mesangiospermae lineages when different genomes were examined. We discuss ways to address these inconsistencies, which could ultimately lead to the reconstruction of a comprehensive angiosperm tree of life. The angiosperm phylogenies presented here provide a basic framework for further updates and comparisons. These phylogenies can also be used as guides to examine the evolutionary trajectories among the three genome types during lineage radiation.

RICE INDETERMINATE 1 (RID1) plays a critical role in controlling floral transition in rice (Oryza sativa). However, the molecular basis for this effect, particularly the target genes and regulatory specificity, remains largely unclear. Here, we performed chromatin immunoprecipitation followed by sequencing (ChIP-seq) in young leaves at the pre-floral-transition stage to identify the target genes of RID1, identifying 2,680 genes associated with RID1 binding sites genome-wide. RID1 binding peaks were highly enriched for TTTGTC, the direct binding motif of the INDETERMINATE DOMAIN protein family that includes RID1. Interestingly, CACGTG and GTGGGCCC, two previously uncharacterized indirect binding motifs, were enriched through the interactions of RID1 with the novel flowering-promoting proteins OsPIL12 and OsTCP11, respectively. Moreover, the ChIP-seq data demonstrated that RID1 bound to numerous rice heading-date genes, such as HEADING DATE 1 (HD1) and FLAVIN-BINDING, KELCH REPEAT, F-BOX 1 (OsFKF1). Notably, transcriptome sequencing (RNA-seq) analysis revealed roles of RID1 in diverse developmental pathways. Genetic analysis combined with genome-wide ChIP-seq and RNA-seq results showed that RID1 directly binds to the promoter of OsERF#136 (a repressor of rice flowering) and negatively regulates its expression. Overall, our findings provide new insights into the molecular and genetic mechanisms underlying rice floral transition and characterize OsERF#136 as a previously unrecognized direct target of RID1.

Phylogenomic evidence from an increasing number of studies has demonstrated that different data sets and analytical approaches often reconstruct strongly supported but conflicting relationships. In this study, 785 single-copy nuclear genes and 75 complete plastomes were used to infer the phylogenetic relationships and estimate the historical biogeography of the apple genus Malus sensu lato, an economically important lineage disjunctly distributed in the Northern Hemisphere and involved in known and suspected hybridization and allopolyploidy events. The nuclear phylogeny recovered the monophyly of Malus s.l. (including Docynia); however, the genus was supported to be biphyletic in the plastid phylogeny. An ancient chloroplast capture event in the Eocene in western North America best explains the cytonuclear discordance. Our conflict analysis demonstrated that ILS, hybridization, and allopolyploidy could explain the widespread nuclear gene tree discordance. One deep hybridization event (Malus doumeri) and one recent event (Malus coronaria) were detected in Malus s.l. Furthermore, our historical biogeographic analysis integrating living and fossil data supported a widespread East Asian-western North American origin of Malus s.l. in the Eocene, followed by several extinction and dispersal events in the Northern Hemisphere. We also propose a general workflow for assessing phylogenomic discordance and biogeographic analysis using deep genome skimming data sets.

Trichoderma harzianum is a plant‐beneficial fungus that secretes small cysteine‐rich proteins that induce plant defense responses; however, the molecular mechanism involved in this induction is largely unknown. Here, we report that the class II hydrophobin ThHyd1 acts as an elicitor of induced systemic resistance (ISR) in plants. Immunogold labeling and immunofluorescence revealed ThHyd1 localized on maize (Zea mays) root cell plasma membranes. To identify host plant protein interactors of Hyd1, we screened a maize B73 root cDNA library. ThHyd1 interacted directly with ubiquilin 1‐like (UBL). Furthermore, the N‐terminal fragment of UBL was primarily responsible for binding with Hyd1 and the eight‐cysteine amino acid of Hyd1 participated in the protein‐protein interactions. Hyd1 from T. harzianum (Thhyd1) and ubl from maize were co‐expressed in Arabidopsis thaliana, they synergistically promoted plant resistance against Botrytis cinerea. RNA‐sequencing analysis of global gene expression in maize leaves 24 h after spraying with Curvularia lunata spore suspension showed that Thhyd1‐induced systemic resistance was primarily associated with brassinosteroid signaling, likely mediated through BAK1. Jasmonate/ethylene (JA/ET) signaling was also involved to some extent in this response. Our results suggest that the Hyd1‐UBL axis might play a key role in inducing systemic resistance as a result of Trichoderma‐plant interactions.

Commercial varieties of upland cotton (Gossypium hirsutum) have undergone extensive breeding for agronomic traits, such as fiber quality, disease resistance, and yield. Cotton breeding programs have widely used Chinese upland cotton source germplasm (CUCSG) with excellent agronomic traits. A better understanding of the genetic diversity and genomic characteristics of these accessions could accelerate the identification of desirable alleles. Here, we analyzed 10,522 high‐quality single‐nucleotide polymorphisms (SNP) with the CottonSNP63K microarray in 137 cotton accessions (including 12 hybrids of upland cotton). These data were used to investigate the genetic diversity, population structure, and genomic characteristics of each population and the contribution of these loci to heterosis. Three subgroups were identified, in agreement with their known pedigrees, geographical distributions, and times since introduction. For each group, we identified lineage‐specific genomic divergence regions, which potentially harbor key alleles that determine the characteristics of each group, such as early maturity‐related loci. Investigation of the distribution of heterozygous loci, among 12 commercial cotton hybrids, revealed a potential role for these regions in heterosis. Our study provides insight into the population structure of upland cotton germplasm. Furthermore, the overlap between lineage‐specific regions and heterozygous loci, in the high‐yield hybrids, suggests a role for these regions in cotton heterosis.