Clear species boundaries are crucial for plans and actions on biodiversity conservation. However, morphological similarities among allied species can result in taxonomic difficulties, thus impeding conservation efforts. In China, Cinnamomum japonicum Siebold is a well-known endangered plant, yet suffers from longstanding taxonomic issues. Here, we explicitly evaluate whether C. japonicum, C. chenii, and C. chekiangense are the same phylogenetic species on the basis of a multi-individual sampling strategy. We identified three sets of low-copy orthologous genes from 19 Lauraceae taxa for phylogenetic inferences. Both the concatenation and coalescent-based phylogenies supported that C. chenii individuals were embedded in the C. japonicum clade, indicating these two taxa are conspecific. Meanwhile, C. chekiangense accessions formed a monophyly which was not sister to C. japonicum. This result, together with the morphological differences that the leaves of C. japonicum are glabrous with a faveolate pattern of venation while those of C. chekiangense have trichomes and inevident lateral veins, led us to consider both as two distinct species. Based on 17 728 neutral single nucleotide polymorphisms (SNPs), the ADMIXTURE analysis suggested that the Chinese C. japonicum populations in Zhoushan Archipelago (=C. chenii) were genetically differentiated from the Japanese and Korean ones. Furthermore, ecological niche modeling predicted that the present distribution area of Chinese C. japonicum is likely to be unsuitable under global warming scenarios. Together with its limited distribution and genetic uniqueness, we recommend that Chinese C. japonicum deserves conservation priorities.

Hybrid genomes usually harbor asymmetrical parental contributions. However, it is challenging to infer the functional significance of asymmetrical retention of parental alleles in hybrid populations of conifer trees. Here we investigated the diversity in the glutathione S-transferase (GST) gene family in a hybrid pine Pinus densata and its parents (Pinus tabuliformis and Pinus yunnanensis). Plant GSTs play major roles in protecting plants against biotic and abiotic stresses. In this study, 19 orthologous groups of GST genes were identified and cloned from these three species. We examined their expression in different tissues, and then purified the corresponding proteins to characterize their enzymatic activities and specificities toward different substrates. We found that among the 19 GST orthologous groups, divergence in gene expression and in enzymatic activities toward different substrates was prevalent. P. densata preferentially retained P. yunnanensis-like GSTs for 17 out of the 19 gene loci. We determined the first GST crystal structure from conifer species at a resolution of 2.19??. Based on this structure, we performed site-directed mutagenesis to replace amino acid residuals in different wild-types of GSTs to understand their functional impacts. Reciprocal replacement of amino acid residuals in native GSTs of P. densata and P. tabuliformis demonstrated significant changes in enzyme functions and identified key sites controlling GSTs activities. This study illustrates an approach to evaluating the functional significance of sequence variations in conifer genomes. Our study also sheds light on plausible mechanisms for controlling the selective retention of parental alleles in the P. densata genome.

Chinese Tajiks are an Indo-Iranian-speaking population in Xinjiang, northwest China. Although the complex demographic history has been characterized, the ancestral sources and genetic admixture of Indo-Iranian-speaking groups in this region remain poorly understood. We here provide the genome-wide genotyping data for over 700?000 single-nucleotide polymorphisms (SNPs) and mtDNA multiplex sequencing data in 64 Chinese male Tajik individuals from two dialect groups, Wakhi and Selekur. We applied principal component analysis (PCA), ADMIXTURE, f-statistics, treemix, qpWave/qpAdm, Admixture-induced Linkage Disequilibrium for Evolutionary Relationships (ALDER), and Fst analyses to infer a fine-scale population genetic structure and admixture history. Our results reveal that Chinese Tajiks showed the closest affinity and similar genetic admixture pattern with ancient Xinjiang populations, especially Xinjiang samples in the historical era. Chinese Tajiks also have gene flow from European and Neolithic Iran farmers-related populations. We observed a genetic substructure in the two Tajik dialect groups. The Selekur-speaking group who lived in the county had more gene flow from East Asians than Wakhi-speaking people who inhabited the village. These results document the population movements contributed to the influx of diverse ancestries in the Xinjiang region.

Transposable elements (TEs) are prevalent components of diverse genomes, and play an important role on the genomic stability and expression regulation of their adjacent genes. It is interesting to know the variation of TE expression and the effects of the presence/absence of TEs on gene expression after hybridization. Here we assessed the expression variation of TEs and the impacts of TEs on expression of nearby genes after hybridization based on comparisons of three pairs of reciprocal F₁ hybrids and four parents in Capsella rubella. Of the 480 TE families expressed in all the four parents and six F₁ hybrids, 7–23 (1.5%–4.2%) TE families were significantly differentially expressed between in silico and real F₁ hybrids, indicating the expression levels of these TE families were affected during hybridization. In particular, there was a Copia TE superfamily and a non-long terminal repeat (non-LTR) TE differentially expressed between the reciprocal F₁ hybrids of 879 and 86IT1, indicating maternal effects may have impacts on expression of TEs in these F₁ hybrids. Besides the impacts on the expression of TE families of the hybridization, genes adjacent to polymorphic TEs tended to show a higher proportion (24.83%) of allele-specific expression (ASE) in F₁ hybrids. Overall, our results highlight the impacts of hybridization on the expression level variation of TEs, and the effects of TEs on ASE after hybridization.

Rhododendron is the largest genus in Ericaceae and is well known for its diversity and beauty of flowers present in different species, making it a much-revered lineage of ornamental plants. Many species of Rhododendron are intolerant of high temperatures, which are becoming more common and intense in urban areas under global climate change. Therefore, the discovery and description of genes from heat-tolerant Rhododendron lineages are essential in the development of new climate-resilient cultivars. One such species known to be heat tolerant is Rhododendron × pulchrum Sweet. To better understand the genomics of heat tolerance in this species, we assembled a haplotype-resolved and chromosome-scale genome for R. × pulchrum, which had a genome size of 509 Mb; a scaffold N50 of 37 251 370 bp; and contained 35 610 genes. In addition, based on the same reannotation pipeline, we conducted pan-genomic analyses for all seven available chromosome-scale Rhododendron genomes and found 14 415 gene groups shared across all species and 18 018 gene groups distributed in the other species, including 1879 gene groups found in only a single species. Finally, we analyzed the transcriptomic data from heat-treated and non-heat-treated R. × pulchrum plants to quantify the genes that are most important during heat stress in an effort to inform the development of climate-resilient cultivars. This study provides insight into the genome diversity in Rhododendron and targets several genes related to agronomic traits that may help in further analysis.

A fundamental question in speciation genomics is how evolutionary processes shape the genomic landscape of differentiation between species. Regions of elevated differentiation, referred to as genomic islands, could be shared among closely related species (shared islands) or specific to a lineage (lineage-specific islands). Shared islands are typically assumed to result from background selection. However, simulations and empirical studies have suggested that positive selection contributes to both shared and lineage-specific islands. Here, we utilized comparative population genomics to examine the contributions of different evolutionary processes to patterns of genetic differentiation when gene flow and incomplete lineage sorting are minimal. We used whole-genome resequencing data for 135 individuals from four oak species, including two independent species pairs, Quercus variabilis Blume and Quercus acutissima Carruth. in the subgenus Cerris, and Quercus dentata Thunb. and Quercus griffithii Hook.f. & Thomson ex Miq. in the subgenus Quercus. We found that both shared and subgenus-specific islands were caused by positive selection, including selective sweeps in current populations and in their most recent common ancestors. Moreover, the recombination rate was a better predictor of genomic differentiation than gene density. Overall, our results reveal that recombination and positive selection impacted genomic differentiation considerably and provide a more precise grasp of how genomic islands formed in Quercus.

The West Liao River (WLR) and Yellow River (YR) basins are two major centers of millet farming in northern China. The result from flotation analyses and the spatial distribution of archeological sites indicate that two distinct survival strategies—agriculture and pastoralism were adopted in the southern and western regions of the WLR. Previous studies of ancient populations from the western area of the WLR suggested a correlation between a pastoral economy in the Bronze Age Upper Xiajiadian culture with a decreased genetic affinity with YR farmers. However, the population history of the southern WLR is unknown mainly due to the lack of ancient genetic data. Here we report the genomic data of an ancient individual from the Majiazishan site from the Late Bronze Age southern WLR region associated with Upper Xiajiadian culture. Unlike individuals from western WLR, this individual derived ancestry entirely from Late Neolithic YR farmers. We found a genetic substructure of the ancient human population of Upper Xiajiadian culture, which is consistent with the differences in the subsistence strategies of western and southern WLR. Climate deterioration led to different populations occupying the west and the south, respectively, in the WLR: the nomadic population from the Amur River (AR) in the west and the agricultural population from the YR in the south.

RNA-binding proteins (RBPs) are components of the post-transcriptional regulatory system, but their regulatory effects on complex traits remain unknown. Using an integrated strategy involving map-based cloning, functional characterizations, and transcriptomic and population genomic analyses, we revealed that RBP-K (LOC_Os08g23120), RBP-A (LOC_Os11g41890), and RBP-J (LOC_Os10g33230) encode proteins that form an RBP-A-J-K complex that negatively regulates rice yield-related traits. Examinations of the RBP-A-J-K complex indicated RBP-K functions as a relatively non-specific RBP chaperone that enables RBP-A and RBP-J to function normally. Additionally, RBP-J most likely affects GA pathways, resulting in considerable increases in grain and panicle lengths, but decreases in grain width and thickness. In contrast, RBP-A negatively regulates the expression of genes most likely involved in auxin-regulated pathways controlling cell wall elongation and carbohydrate transport, with substantial effects on the rice grain filling process as well as grain length and weight. Evolutionarily, RBP-K is relatively ancient and highly conserved, whereas RBP-J and RBP-A are more diverse. Thus, the RBP-A-J-K complex may represent a typical functional model for many RBPs and protein complexes that function at transcriptional and post-transcriptional levels in plants and animals for increased functional consistency, efficiency, and versatility, as well as increased evolutionary potential. Our results clearly demonstrate the importance of RBP-mediated post-transcriptional regulation for the diversity of complex traits. Furthermore, rice grain yield and quality may be enhanced by introducing various complete or partial loss-of-function mutations to specific RBP genes using clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein 9 technology and by exploiting desirable natural tri-genic allelic combinations at the loci encoding the components of the RBP-A-J-K complex through marker-assisted selection.

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

Many examples of phenotypic modifications resulting from high-elevation adaptation have been documented, however, the underlying processes responsible for these modifications and whether the continuity of the adaptation process remain elusive, particularly in plants. The alpine plants distributed along wide elevational gradients provide an ideal system to address this question. Here, we collected transcriptomes from multiple tissues of three species with different elevations (~1500, ~2500, and ~3600 m in the Hengduan Mountains, Southwest China) in two genera Fargesia and Yushania of alpine bamboos, respectively, and conducted evolutionary and expressional analyses. Results showed that high-elevation adaptation emerged earlier in the evolutionary history of both genera and evolved continuously as the elevation increased. Moreover, convergence of genetic changes was observed in the two genera, with amounts of candidate genes responsible for high-elevation adaptation identified under positive selection. Overall, our study provides an empirical example and valuable genetic resource for further investigation of high-elevation adaptation in plants and sheds new light on how plants adapting to high-elevation environments in a biodiversity hotspot.

Plants are highly susceptible to abiotic stresses, particularly heat stress during the reproductive stage. However, the specific molecular mechanisms underlying this sensitivity remain largely unknown. In the current study, we demonstrate that the Nuclear Transcription Factor, X-box Binding Protein 1-Like 1 (NFXL1), directly regulates the expression of DEHYDRATION-RESPONSIVE ELEMENT-BINDING PROTEIN 2A (DREB2A), which is crucial for reproductive thermotolerance in Arabidopsis. NFXL1 is upregulated by heat stress, and its mutation leads to a reduction in silique length (seed number) under heat stress conditions. RNA-Seq analysis reveals that NFXL1 has a global impact on the expression of heat stress responsive genes, including DREB2A, Heat Shock Factor A3 (HSFA3) and Heat Shock Protein 17.6 (HSP17.6) in flower buds. Interestingly, NFXL1 is enriched in the promoter region of DREB2A, but not of either HSFA3 or HSP17.6. Further experiments using electrophoretic mobility shift assay have confirmed that NFXL1 directly binds to the DNA fragment derived from the DREB2A promoter. Moreover, effector–reporter assays have shown that NFXL1 activates the DREB2A promoter. The DREB2A mutants are also heat stress sensitive at the reproductive stage, and DEREB2A is epistatic to NFXL1 in regulating thermotolerance in flower buds. It is known that HSFA3, a direct target of DREB2A, regulates the expression of heat shock proteins genes under heat stress conditions. Thus, our findings establish NFXL1 as a critical upstream regulator of DREB2A in the transcriptional cassette responsible for heat stress responses required for reproductive thermotolerance in Arabidopsis.

Transcriptional regulation is essential for balancing multiple metabolic pathways that influence oil accumulation in seeds. Thus far, the transcriptional regulatory mechanisms that govern seed oil accumulation remain largely unknown. Here, we identified the transcriptional regulatory network composed of MADS-box transcription factors SEEDSTICK (STK) and SEPALLATA3 (SEP3), which bridges several key genes to regulate oil accumulation in seeds. We found that STK, highly expressed in the developing embryo, positively regulates seed oil accumulation in Arabidopsis (Arabidopsis thaliana). Furthermore, we discovered that SEP3 physically interacts with STK in vivo and in vitro. Seed oil content is increased by the SEP3 mutation, while it is decreased by SEP3 overexpression. The chromatin immunoprecipitation, electrophoretic mobility shift assay, and transient dual-luciferase reporter assays showed that STK positively regulates seed oil accumulation by directly repressing the expression of MYB5, SEP3, and SEED FATTY ACID REDUCER 4 (SFAR4). Moreover, genetic and molecular analyses demonstrated that STK and SEP3 antagonistically regulate seed oil production and that SEP3 weakens the binding ability of STK to MYB5, SEP3, and SFAR4. Additionally, we demonstrated that TRANSPARENT TESTA 8 (TT8) and ACYL-ACYL CARRIER PROTEIN DESATURASE 3 (AAD3) are direct targets of MYB5 during seed oil accumulation in Arabidopsis. Together, our findings provide the transcriptional regulatory network antagonistically orchestrated by STK and SEP3, which fine tunes oil accumulation in seeds.

Rice (Oryza sativa) is a significant crop worldwide with a genome shaped by various evolutionary factors. Rice centromeres are crucial for chromosome segregation, and contain some unreported genes. Due to the diverse and complex centromere region, a comprehensive understanding of rice centromere structure and function at the population level is needed. We constructed a high-quality centromere map based on the rice super pan-genome consisting of a 251-accession panel comprising both cultivated and wild species of Asian and African rice. We showed that rice centromeres have diverse satellite repeat CentO, which vary across chromosomes and subpopulations, reflecting their distinct evolutionary patterns. We also revealed that long terminal repeats (LTRs), especially young Gypsy-type LTRs, are abundant in the peripheral CentO-enriched regions and drive rice centromere expansion and evolution. Furthermore, high-quality genome assembly and complete telomere-to-telomere (T2T) reference genome enable us to obtain more centromeric genome information despite mapping and cloning of centromere genes being challenging. We investigated the association between structural variations and gene expression in the rice centromere. A centromere gene, OsMAB, which positively regulates rice tiller number, was further confirmed by expression quantitative trait loci, haplotype analysis and clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein 9 methods. By revealing the new insights into the evolutionary patterns and biological roles of rice centromeres, our finding will facilitate future research on centromere biology and crop improvement.

Anthocyanins are secondary metabolites induced by environmental stimuli and developmental signals. The positive regulators of anthocyanin biosynthesis have been reported, whereas the anthocyanin repressors have been neglected. Although the signal transduction pathways of gibberellin (GA) and jasmonic acid (JA) and their regulation of anthocyanin biosynthesis have been investigated, the cross-talk between GA and JA and the antagonistic mechanism of regulating anthocyanin biosynthesis remain to be investigated. In this study, we identified the anthocyanin repressor MdbHLH162 in apple and revealed its molecular mechanism of regulating anthocyanin biosynthesis by integrating the GA and JA signals. MdbHLH162 exerted passive repression by interacting with MdbHLH3 and MdbHLH33, which are two recognized positive regulators of anthocyanin biosynthesis. MdbHLH162 negatively regulated anthocyanin biosynthesis by disrupting the formation of the anthocyanin-activated MdMYB1-MdbHLH3/33 complexes and weakening transcriptional activation of the anthocyanin biosynthetic genes MdDFR and MdUF3GT by MdbHLH3 and MdbHLH33. The GA repressor MdRGL2a antagonized MdbHLH162-mediated inhibition of anthocyanins by sequestering MdbHLH162 from the MdbHLH162-MdbHLH3/33 complex. The JA repressors MdJAZ1 and MdJAZ2 interfered with the antagonistic regulation of MdbHLH162 by MdRGL2a by titrating the formation of the MdRGL2a-MdbHLH162 complex. Our findings reveal that MdbHLH162 integrates the GA and JA signals to negatively regulate anthocyanin biosynthesis. This study provides new information for discovering more anthocyanin biosynthesis repressors and explores the cross-talk between hormone signals.

As an important halophyte in the Yellow River Delta, the Amaranthaceae C₃ Suaeda salsa (L.) Pall. has attracted much attention for the “red carpet” landscape, and could be simply divided into red and green phenotypes according to the betacyanin content in the fleshy leaves. However, S. salsa has not been sequenced yet, which limited people's understanding of this species at the molecular level. We constructed a high-quality chromosome-scale reference genome by combining high-throughput sequencing, PacBio single molecule real-time sequencing, and Hi-C sequencing techniques with a genome size of 445.10 Mb and contigs N50 of 2.94 Mb. Through the annotation of the S. salsa genome, 298.76 Mb of the repetitive sequences and 23 965 protein-coding genes were identified, of which the proportion of long terminal repeats type in the repetitive sequences was the most abundant, about 50.74% of the S. salsa genome. Comparative genomics indicated that S. salsa underwent a whole-genome duplication event about 146.15 million years ago (Ma), and the estimated divergence time between S. salsa and Suaeda aralocaspica was about 16.9 Ma. A total of four betacyanins including betanidin, celosianin II, amaranthin and 6′-O-malonyl-celosianin II were identified and purified in both phenotypes, while two significantly up-regulated betacyanins (celosianin II and amaranthin) may be the main reason for the red color in red phenotype. In addition, we also performed transcriptomics and metabolomics in both phenotypes to explore the molecular mechanisms of pigment synthesis, and a series of structural genes and transcription factors concerning with betacyanin production were selected in S. salsa.

In plants, the genome structure of hybrids changes compared with their parents, but the effects of these changes in hybrids remain elusive. Comparing reciprocal crosses between Col×C24 and C24×Col in Arabidopsis using high-throughput chromosome conformation capture assay (Hi-C) analysis, we found that hybrid three-dimensional (3D) chromatin organization had more long-distance interactions relative to parents, and this was mainly located in promoter regions and enriched in genes with heterosis-related pathways. The interactions between euchromatin and heterochromatin were increased, and the compartment strength decreased in hybrids. In compartment domain (CD) boundaries, the distal interactions were more in hybrids than their parents. In the hybrids of CURLY LEAF (clf) mutants clf_Col×clf_C24 and clf_C24×clf_Col, the heterosis phenotype was damaged, and the long-distance interactions in hybrids were fewer than in their parents with lower H3K27me3. ChIP-seq data revealed higher levels of H3K27me3 in the region adjacent to the CD boundary and the same interactional homo-trans sites in the wild-type (WT) hybrids, which may have led to more long-distance interactions. In addition, the differentially expressed genes (DEGs) located in the boundaries of CDs and loop regions changed obviously in WT, and the functional enrichment for DEGs was different between WT and clf in the long-distance interactions and loop regions. Our findings may therefore propose a new epigenetic explanation of heterosis in the Arabidopsis hybrids and provide new insights into crop breeding and yield increase.

Drought is a major threat to alfalfa (Medicago sativa L.) production. The discovery of important alfalfa genes regulating drought response will facilitate breeding for drought-resistant alfalfa cultivars. Here, we report a genome-wide association study of drought resistance in alfalfa. We identified and functionally characterized an MYB-like transcription factor gene (MsMYBH), which increases the drought resistance in alfalfa. Compared with the wild-types, the biomass and forage quality were enhanced in MsMYBH overexpressed plants. Combined RNA-seq, proteomics and chromatin immunoprecipitation analysis showed that MsMYBH can directly bind to the promoters of MsMCP1, MsMCP2, MsPRX1A and MsCARCAB to improve their expression. The outcomes of such interactions include better water balance, high photosynthetic efficiency and scavenge excess H₂O₂ in response to drought. Furthermore, an E3 ubiquitin ligase (MsWAV3) was found to induce MsMYBH degradation under long-term drought, via the 26S proteasome pathway. Furthermore, variable-number tandem repeats in MsMYBH promoter were characterized among a collection of germplasms, and the variation is associated with promoter activity. Collectively, our findings shed light on the functions of MsMYBH and provide a pivotal gene that could be leveraged for breeding drought-resistant alfalfa. This discovery also offers new insights into the mechanisms of drought resistance in alfalfa.

Nypa fruticans (Wurmb), a mangrove palm species with origins dating back to the Late Cretaceous period, is a unique species for investigating long-term adaptation strategies to intertidal environments and the early evolution of palms. Here, we present a chromosome-level genome sequence and assembly for N. fruticans. We integrated the genomes of N. fruticans and other palm family members for a comparative genomic analysis, which confirmed that the common ancestor of all palms experienced a whole-genome duplication event around 89 million years ago, shaping the distinctive characteristics observed in this clade. We also inferred a low mutation rate for the N. fruticans genome, which underwent strong purifying selection and evolved slowly, thus contributing to its stability over a long evolutionary period. Moreover, ancient duplicates were preferentially retained, with critical genes having experienced positive selection, enhancing waterlogging tolerance in N. fruticans. Furthermore, we discovered that the pseudogenization of Early Methionine-labelled 1 (EM1) and EM6 in N. fruticans underly its crypto-vivipary characteristics, reflecting its intertidal adaptation. Our study provides valuable genomic insights into the evolutionary history, genome stability, and adaptive evolution of the mangrove palm. Our results also shed light on the long-term adaptation of this species and contribute to our understanding of the evolutionary dynamics in the palm family.

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 Sui people living in Guizhou province have a unique ethnic culture and population history due to their long-time isolation from other populations. To investigate the genetic structure of Sui populations in different regions of Guizhou, we genotyped 89 individuals from four Sui populations using genome-wide single nucleotide polymorphisms arrays. We analyzed the data using principal component analysis, ADMIXTURE analysis, f-statistics, qpWave/qpAdm, TreeMix analysis, fineSTRUCTURE, and GLOBETROTTER. We found that Sui populations in Guizhou were genetically homogeneous and had a close genetic affinity with Tai-Kadai-speaking populations, Hmong-Mien-speaking Hmong, and some ancient populations from southern China. The Sui populations could be modeled as an admixture of 33.5%–37.9% of Yellow River Basin farmer-related ancestry and 62.1%–66.5% of Southeast Asian-related ancestry, indicating that the southward expansion of northern East Asian-related ancestry influenced the formation of the Tai-Kadai-speaking Sui people. Future publications of more ancient genomics in southern China could effectively provide further insight into the demographic history and population structure of the Sui people.

Auxin regulates flower and fruit abscission, but how developmental signals mediate auxin transport in abscission remains unclear. Here, we reveal the role of the transcription factor BEL1-LIKE HOMEODOMAIN11 (SlBEL11) in regulating auxin transport during abscission in tomato (Solanum lycopersicum). SlBEL11 is highly expressed in the fruit abscission zone, and its expression increases during fruit development. Knockdown of SlBEL11 expression by RNA interference (RNAi) caused premature fruit drop at the breaker (Br) and 3d post-breaker (Br+3) stages of fruit development. Transcriptome and metabolome analysis of SlBEL11-RNAi lines revealed impaired flavonoid biosynthesis and decreased levels of most flavonoids, especially quercetin, which functions as an auxin transport inhibitor. This suggested that SlBEL11 prevents premature fruit abscission by modulating auxin efflux from fruits, which is crucial for the formation of an auxin response gradient. Indeed, quercetin treatment suppressed premature fruit drop in SlBEL11-RNAi plants. DNA affinity purification sequencing (DAP-seq) analysis indicated that SlBEL11 induced expression of the transcription factor gene SlMYB111 by directly binding to its promoter. Chromatin immunoprecipitation-quantitative polymerase chain reaction and electrophoretic mobility shift assay showed that S. lycopersicum MYELOBLASTOSIS VIRAL ONCOGENE HOMOLOG111 (SlMYB111) induces the expression of the core flavonoid biosynthesis genes SlCHS1, SlCHI, SlF3H, and SlFLS by directly binding to their promoters. Our findings suggest that the SlBEL11-SlMYB111 module modulates flavonoid biosynthesis to fine-tune auxin efflux from fruits and thus maintain an auxin response gradient in the pedicel, thereby preventing premature fruit drop.

Achieving high yield and good quality in crops is essential for human food security and health. However, there is usually disharmony between yield and quality. Seed storage protein (SSP) and starch, the predominant components in cereal grains, determine yield and quality, and their coupled synthesis causes a yield–quality trade-off. Therefore, dissection of the underlying regulatory mechanism facilitates simultaneous improvement of yield and quality. Here, we summarize current findings about the synergistic molecular machinery underpinning SSP and starch synthesis in the leading staple cereal crops, including maize, rice and wheat. We further evaluate the functional conservation and differentiation of key regulators and specify feasible research approaches to identify additional regulators and expand insights. We also present major strategies to leverage resultant information for simultaneous improvement of yield and quality by molecular breeding. Finally, future perspectives on major challenges are proposed.

The basis of modern pharmacology is the human ability to exploit the production of specialized metabolites from medical plants, for example, terpenoids, alkaloids, and phenolic acids. However, in most cases, the availability of these valuable compounds is limited by cellular or organelle barriers or spatio-temporal accumulation patterns within different plant tissues. Transcription factors (TFs) regulate biosynthesis of these specialized metabolites by tightly controlling the expression of biosynthetic genes. Cutting-edge technologies and/or combining multiple strategies and approaches have been applied to elucidate the role of TFs. In this review, we focus on recent progress in the transcription regulation mechanism of representative high-value products and describe the transcriptional regulatory network, and future perspectives are discussed, which will help develop high-yield plant resources.

ChinaMu is the largest sequence-indexed Mutator (Mu) transposon insertional library in maize (Zea mays). In this study, we made significant improvements to the size and quality of the ChinaMu library. We developed a new Mu-tag isolation method Mu-Tn5-seq (MuT-seq). Compared to the previous method used by ChinaMu, MuT-seq recovered 1/3 more germinal insertions, while requiring only about 1/14 of the sequencing volume and 1/5 of the experimental time. Using MuT-seq, we identified 113,879 germinal insertions from 3,168 Mu-active F₁ families. We also assembled a high-quality genome for the Mu-active line Mu-starter, which harbors the initial active MuDR element and was used as the pollen donor for the mutation population. Using the Mu-starter genome, we recovered 33,662 (15.6%) additional germinal insertions in 3,244 (7.4%) genes in the Mu-starter line. The Mu-starter genome also improved the assignment of 117,689 (54.5%) germinal insertions. The newly upgraded ChinaMu dataset currently contains 215,889 high-quality germinal insertions. These insertions cover 32,224 pan-genes in the Mu-starter and B73Ref5 genomes, including 23,006 (80.4%) core genes shared by the two genomes. As a test model, we investigated Mu insertions in the pentatricopeptide repeat (PPR) superfamily, discovering insertions for 92% (449/487) of PPR genes in ChinaMu, demonstrating the usefulness of ChinaMu as a functional genomics resource for maize.

In plants, thousands of nucleus-encoded proteins translated in the cytosol are sorted to chloroplasts and mitochondria by binding to specific receptors of the TOC (translocon on the outer chloroplast membrane) and the TOM (translocon on the outer mitochondrial membrane) complexes for import into those organelles. The degradation pathways for these receptors are unclear. Here, we discovered a converged ubiquitin-proteasome pathway for the degradation of Arabidopsis thaliana TOC and TOM tail-anchored receptors. The receptors are ubiquitinated by E3 ligase(s) and pulled from the outer membranes by the AAA⁺ adenosine triphosphatase CDC48, after which a previously uncharacterized cytosolic protein, transmembrane domain (TMD)-binding protein for tail-anchored outer membrane proteins (TTOP), binds to the exposed TMDs at the C termini of the receptors and CDC48, and delivers these complexes to the 26S proteasome.

Camptothecin is a complex monoterpenoid indole alkaloid with remarkable antitumor activity. Given that two C-10 modified camptothecin derivatives, topotecan and irinotecan, have been approved as potent anticancer agents, there is a critical need for methods to access other aromatic ring-functionalized congeners (e.g., C-9, C-10, etc.). However, contemporary methods for chemical oxidation are generally harsh and low-yielding when applied to the camptothecin scaffold, thereby limiting the development of modified derivatives. Reported herein, we have identified four tailoring enzymes responsible for C-9 modifications of camptothecin from Nothapodytes tomentosa, via metabolomic and transcriptomic analysis. These consist of a cytochrome P450 (NtCPT9H) which catalyzes the regioselective oxidation of camptothecin to 9-hydroxycamptothecin, as well as two methyltransferases (NtOMT1/2, converting 9-hydroxycamptothecin to 9-methoxycamptothecin), and a uridine diphosphate-glycosyltransferase (NtUGT5, decorating 9-hydroxycamptothecin to 9-β-D-glucosyloxycamptothecin). Importantly, the critical residues that contribute to the specific catalytic activity of NtCPT9H have been elucidated through molecular docking and mutagenesis experiments. This work provides a genetic basis for producing camptothecin derivatives through metabolic engineering. This will hasten the discovery of novel C-9 modified camptothecin derivatives, with profound implications for pharmaceutical manufacture.

Reynoutria multiflora is a widely used medicinal plant in China. Its medicinal compounds are mainly stilbenes and anthraquinones which possess important pharmacological activities in anti‐aging, anti‐inflammatory and anti‐oxidation, but their biosynthetic pathways are still largely unresolved. Here, we reported a near‐complete genome assembly of R. multiflora consisting of 1.39 Gb with a contig N50 of 122.91 Mb and only one gap left. Genome evolution analysis revealed that two recent bursts of long terminal repeats (LTRs) contributed significantly to the increased genome size of R. multiflora, and numerous large chromosome rearrangements were observed between R. multiflora and Fagopyrum tataricum genomes. Comparative genomics analysis revealed that a recent whole‐genome duplication specific to Polygonaceae led to a significant expansion of gene families associated with disease tolerance and the biosynthesis of stilbenes and anthraquinones in R. multiflora. Combining transcriptomic and metabolomic analyses, we elucidated the molecular mechanisms underlying the dynamic changes in content of medicinal ingredients in R. multiflora roots across different growth years. Additionally, we identified several putative key genes responsible for anthraquinone and stilbene biosynthesis. We identified a stilbene synthase gene PM0G05131 highly expressed in roost, which may exhibit an important role in the accumulation of stilbenes in R. multiflora. These genomic data will expedite the discovery of anthraquinone and stilbenes biosynthesis pathways in medicinal plants.

By 2050, the global population is projected to reach 9 billion, underscoring the imperative for innovative solutions to increase grain yield and enhance food security. Nanotechnology has emerged as a powerful tool, providing unique solutions to this challenge. Nanoparticles (NPs) can improve plant growth and nutrition under normal conditions through their high surface-to-volume ratio and unique physical and chemical properties. Moreover, they can be used to monitor crop health status and augment plant resilience against abiotic stresses (such as salinity, drought, heavy metals, and extreme temperatures) that endanger global agriculture. Application of NPs can enhance stress tolerance mechanisms in plants, minimizing potential yield losses and underscoring the potential of NPs to raise crop yield and quality. This review highlights the need for a comprehensive exploration of the environmental implications and safety of nanomaterials and provides valuable guidelines for researchers, policymakers, and agricultural practitioners. With thoughtful stewardship, nanotechnology holds immense promise in shaping environmentally sustainable agriculture amid escalating environmental challenges.

Structural variations (SVs) are a feature of plant genomes that has been largely unexplored despite their significant impact on plant phenotypic traits and local adaptation to abiotic and biotic stress. In this study, we employed woolly grape (Vitis retordii), a species native to the tropical and subtropical regions of East Asia with both coastal and inland habitats, as a valuable model for examining the impact of SVs on local adaptation. We assembled a haplotype-resolved chromosomal reference genome for woolly grape, and conducted population genetic analyses based on whole- genome sequencing (WGS) data from coastal and inland populations. The demographic analyses revealed recent bottlenecks in all populations and asymmetric gene flow from the inland to the coastal population. In total, 1,035 genes associated with plant adaptive regulation for salt stress, radiation, and environmental adaptation were detected underlying local selection by SVs and SNPs in the coastal population, of which 37.29% and 65.26% were detected by SVs and SNPs, respectively. Candidate genes such as FSD2, RGA1, and AAP8 associated with salt tolerance were found to be highly differentiated and selected during the process of local adaptation to coastal habitats in SV regions. Our study highlights the importance of SVs in local adaptation; candidate genes related to salt stress and climatic adaptation to tropical and subtropical environments are important genomic resources for future breeding programs of grapevine and its rootstocks.

Drought stress is a crucial environmental factor that limits plant growth, development, and productivity. Autophagy of misfolded proteins can help alleviate the damage caused in plants experiencing drought. However, the mechanism of autophagy-mediated drought tolerance in plants remains largely unknown. Here, we cloned the gene for a maize (Zea mays) selective autophagy receptor, NEXT TO BRCA1 GENE 1 (ZmNBR1), and identified its role in the response to drought stress. We observed that drought stress increased the accumulation of autophagosomes. RNA sequencing and reverse transcription-quantitative polymerase chain reaction showed that ZmNBR1 is markedly induced by drought stress. ZmNBR1 overexpression enhanced drought tolerance, while its knockdown reduced drought tolerance in maize. Our results established that ZmNBR1 mediates the increase in autophagosomes and autophagic activity under drought stress. ZmNBR1 also affects the expression of genes related to autophagy under drought stress. Moreover, we determined that BRASSINOSTEROID INSENSITIVE 1A (ZmBRI1a), a brassinosteroid receptor of the BRI1-like family, interacts with ZmNBR1. Phenotype analysis showed that ZmBRI1a negatively regulates drought tolerance in maize, and genetic analysis indicated that ZmNBR1 acts upstream of ZmBRI1a in regulating drought tolerance. Furthermore, ZmNBR1 facilitates the autophagic degradation of ZmBRI1a under drought stress. Taken together, our results reveal that ZmNBR1 regulates the expression of autophagy-related genes, thereby increasing autophagic activity and promoting the autophagic degradation of ZmBRI1a under drought stress, thus enhancing drought tolerance in maize. These findings provide new insights into the autophagy degradation of brassinosteroid signaling components by the autophagy receptor NBR1 under drought stress.

The ancient crop broomcorn millet (Panicum miliaceum L.) is an indispensable orphan crop in semi-arid regions due to its short life cycle and excellent abiotic stress tolerance. These advantages make it an important alternative crop to increase food security and achieve the goal of zero hunger, particularly in light of the uncertainty of global climate change. However, functional genomic and biotechnological research in broomcorn millet has been hampered due to a lack of genetic tools such as transformation and genome-editing techniques. Here, we successfully performed genome editing of broomcorn millet. We identified an elite variety, Hongmi, that produces embryogenic callus and has high shoot regeneration ability in in vitro culture. We established an Agrobacterium tumefaciens-mediated genetic transformation protocol and a clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9-mediated genome-editing system for Hongmi. Using these techniques, we produced herbicide-resistant transgenic plants and edited phytoene desaturase (PmPDS), which is involved in chlorophyll biosynthesis. To facilitate the rapid adoption of Hongmi as a model line for broomcorn millet research, we assembled a near-complete genome sequence of Hongmi and comprehensively annotated its genome. Together, our results open the door to improving broomcorn millet using biotechnology.