Genome-wide association analyses provide genetic and biochemical insights into natural variation in rice metabolism

Transcripts and genomic intervals associated with variation in metabolite abundance in maize leaves under field conditions

Plants exhibit extensive environment-dependent intraspecific metabolic variation, which likely plays a role in determining variation in whole plant phenotypes. However, much of the work seeking to use natural variation to link genes and transcript's impacts on plant metabolism has employed data from controlled environments. Here, we generated and analyzed data on the variation in the abundance of 26 metabolites across 660 maize inbred lines under field conditions. We employ these data and previously published transcript and whole plant phenotype data reported for the same field experiment to identify both genomic intervals (through genome-wide association studies (GWAS)) and transcripts (using both transcriptome-wide association studies (TWAS) and an explainable artificial intelligence (AI) approach based on random forest (RF)) associated with variation in metabolite abundance. Both genome-wide association and random forest-based methods identified substantial numbers of significant associations including genes with plausible links to the metabolites they are associated with. In contrast, the transcriptome-wide association identified only six significant associations. In three cases, genetic markers associated with metabolic variation in our study colocalized with markers linked to variation in non-metabolic traits scored in the same experiment. We speculate that the poor performance of transcriptome-wide association studies in identifying transcript-metabolite associations may reflect a high prevalence of non-linear interactions between transcripts and metabolites and/or a bias towards rare transcripts playing a large role in determining intraspecific metabolic variation.

Comprehensive Analysis of Ghd7 Variations Using Pan-Genomics and Prime Editing in Rice

The Ghd7 gene in rice plays a crucial role in determining heading date, plant height, and grain yield. However, the variations in Ghd7 and their functional implications across different rice accessions are not fully understood. Based on the release of a large amount of rice genome data in recent years, we investigated Ghd7 through pan-genome analysis of 372 diverse rice varieties and figured out the structural variations (SVs) in the Ghd7 locus. However, due to the high cost of pan-genomes, most genomes are based on next-generation sequencing (NGS) data now. Therefore, we developed a method for identifying SVs using NGS data and Polymerase Chain Reaction (PCR) based on the results of pan-genome analysis and identified 977 accessions carrying such SVs of Ghd7. Furthermore, we identified 46 single-nucleotide polymorphisms (SNPs) and one insertion-deletion (InDel) in the coding region of Ghd7. They are classified into 49 haplotypes. Notably, a splice-site mutation in haplotype H6 causes aberrant mRNA splicing. Using prime editing (PE) technology, we successfully restored the functional of Ghd7 in Yixiang 1B (YX1B), delaying the heading date by approximately 16 days. This modification synchronized the heading date between YX1B and the restorer line Yahui 2115 (YH2115R), enhancing the hybrid rice seed production efficiency. In conclusion, our findings highlight the potential of integrating pan-genomics and precision gene editing to accelerate crop improvement and enhance agronomic traits.

Large-scale metabolomic landscape of edible maize reveals convergent changes in metabolite differentiation and facilitates its breeding improvement

Edible maize is an important food crop that provides energy and nutrients to meet human health and nutritional requirements. However, how environmental pressures and human activity have shaped the metabolome of edible maize remains unclear. In this study, we collected 452 diverse edible maize accessions worldwide, including waxy, sweet, and field maize. A total of 3020 non-redundant metabolites, including 802 annotated metabolites, were identified using a two-step optimized approach, which generated the most comprehensive annotated metabolite dataset in plants to date. Although specific metabolite differentiation was detected between field and sweet maize and between field and waxy maize, convergent metabolite differentiation was the dominant pattern. We identified hub genes in all metabolite classes by hotspot analysis in a metabolite genome-wide association study. Seventeen and 15 hub genes were selected as the key differentiation genes for flavonoids and lipids, respectively. Surprisingly, almost all of these genes were under diversifying selection, suggesting that diversifying selection was the main genetic mechanism of convergent metabolic differentiation. Further genetic and molecular studies revealed the roles and genetic diversifying selection mechanisms of ZmGPAT11 in convergent metabolite differentiation in the lipid pathway. On the basis of our research, we established the first edible maize metabolome database, EMMDB (https://www.maizemdb.site/home/). We successfully used EMMDB for precision improvement of nutritional and flavor traits and bred the elite inbred line 6644_2, with greatly increased contents of flavonoids, lysophosphatidylcholines, lysophosphatidylethanolamines, and vitamins. Collectively, our study sheds light on the genetic mechanisms of metabolite differentiation in edible maize and provides a database for breeding improvement of flavor and nutritional traits in edible maize by metabolome precision design.

Metabolome genome-wide association analyses identify a splice mutation in AADAT affects lysine degradation in duck skeletal muscle

Metabolites in skeletal muscles play an important role in their growth, development, immunity and other physiological activities. However, the genetic basis of metabolites in skeletal muscle remains poorly understood. Here, we identified 247 candidate divergent regions containing 905 protein-coding genes closely related to metabolic pathways, including lysine degradation and fatty acid biosynthesis. We then profiled 3,060 metabolites in 246 skeletal muscle samples from F2 segregating population generated by mallard×Pekin duck crosses using metabolomic approaches. We identified 2,044 significant metabolome-based GWAS signals and 21 candidate genes potentially modulating metabolite contents in skeletal muscle. Among them, the levels of 2-aminoadipic acid in skeletal muscle were significantly correlated with body weight and intramuscular fat content, determined by a 939-bp CR1 LINE insertion in AADAT. We further found that the CR1 LINE insertion most possibly led to a splice mutation in AADAT, resulting in the downregulation of the lysine degradation pathway in skeletal muscle. Moreover, intramuscular fat content and fatty acids biosynthesis pathway was significantly increased in individuals with CR1 LINE insertion. This study enhances our understanding of the genetic basis of skeletal muscle metabolic traits and promotes the efficient utilization of metabolite traits in the genetic improvement of animals.

Advances and Future Prospects of Pigment Deposition in Pigmented Rice

Pigmented rice, particularly the black and red varieties, is popular due to its better nutritional value. Anthocyanins and proanthocyanidins are two major flavonoid subcategories with broad physiological functions and therapeutic significance. However, pigment deposition is a complex process, and the molecular mechanism involved remains unknown. This review explores the metabolites responsible for the pigmentation in various rice tissues. Moreover, the current challenges, feasible strategies, and potential future directions in pigmented rice research are reported.

Transcriptome analysis reveals the key roles of TaSMP1 and ABA signaling pathway in wheat seed dormancy and germination

MAIN CONCLUSION

This study analyzed dynamic transcriptome profiles to reveal differential expression patterns of ABA related and LEA protein family genes and verified that TaSMP1 affects seed germination by interacting with TaABI5. Seed dormancy is a crucial survival strategy for plants to cope with environmental stresses. High levels of seed dormancy result in uneven germination, while low levels of seed dormancy increase the risk of pre-harvest sprouting (PHS), which threatens crop yield and quality. Therefore, achieving the optimal balance between seed dormancy and germination is vital for maximum potential crop yield and quality. This study constructed dynamic transcriptome expression profiles of the germination process for the weakly dormant wheat variety Jing 411 (J411) and the strongly dormant landrace variety Hongsuibai (HSB), revealing the temporal expression of differentially expressed genes. Plant hormone-related genes played a crucial role in the early germination response, particularly the abscisic acid (ABA) signaling gene TaABI5 and the ABA catabolism gene TaCYP707A1. The late embryogenesis abundant (LEA) protein family genes exhibited differential expression patterns during the germination of seeds with varying levels of dormancy. The TaSMP1 gene, a member of the LEA protein family, was identified as a negative regulator of seed dormancy, interacting directly with the key transcription factor TaABI5 in the ABA signaling pathway and influencing the expression of the seed germination gene TaDOG1L1. This study provides essential insights into the molecular mechanisms balancing seed dormancy and germination, offering potential targets for enhancing wheat resistance to PHS.

Genome-wide profiling of polymorphic short tandem repeats and their influence on gene expression and trait variation in diverse rice populations

Short tandem repeats (STRs) modulate gene expression and contribute to trait variation. However, a systematic evaluation of the genomic characteristics of STRs has not been conducted, and their influence on gene expression in rice remains unclear. Here, we construct a map of 137,629 polymorphic STRs in the rice (Oryza sativa L.) genome using a population-scale resequencing dataset. A genome-wide survey encompassing 4726 accessions shows that the occurrence frequency, mutational patterns, chromosomal distribution, and functional properties of STRs are correlated with the sequences and lengths of repeat motifs. Leveraging a transcriptome dataset from 127 rice accessions, we identify 44,672 expression STRs (eSTRs) by modeling gene expression in response to the length variation of STRs. These eSTRs are notably enriched in the regulatory regions of genes with active transcriptional signatures. Population analysis identifies numerous STRs that have undergone genetic divergence among different rice groups and 1726 tagged STRs that may be associated with agronomic traits. By editing the (ACT)7 STR in OsFD1 promoter, we further experimentally validate its role in regulating gene expression and phenotype. Our study highlights the contribution of STRs to transcriptional regulation in plants and establishes the foundation for their potential use as alternative targets for genetic improvement.

Integrative phenomics, metabolomics and genomics analysis provides new insights for deciphering the genetic basis of metabolism in polished rice

BACKGROUND

Metabolomics is one of the most widely used omics tools for deciphering the functional networks of the metabolites for crop improvement. However, it is technically demanding and costly.

RESULTS

We propose a relatively inexpensive approach for metabolomics analysis in which metabolomics is combined with hyperspectral imaging via machine learning. This approach can be used to target important steps in flavonoid and lipid biosynthesis in rice. We extract 1848 hyperspectral indices and 887 metabolites from polished grains of 533 Oryza sativa accessions. Hyperspectral indices are then linked to metabolites through correlation analysis and modelling. Based on this, a total of 554 metabolites and 1313 hyperspectral indices are identified for further genome-wide association study (GWAS). By GWAS, we detect 17,509 significant locus-trait associations with 2882 single nucleotide polymorphisms (SNPs). Colocalization analysis links these SNPs to the corresponding metabolites and hyperspectral indices. We detect 6415 pairs of metabolites and hyperspectral indices within a linkage disequilibrium of 300 kb in the Oryza sativa genome. We then characterize 1761 candidate genes colocalized to these loci by transcriptomic analysis. We further verify novel candidate genes encoding a novel flavonoid (LOC_Os09g18450) and a flavonoid/lipid (LOC_Os07g11020) respectively by gene editing and overexpression in rice.

CONCLUSION

Our findings indicate that hyperspectral imaging combined with machine learning methods could serve as a powerful tool for quickly and inexpensively assessing crop metabolites.

Rationale and design of the Dog Aging Project precision cohort: a multi-omic resource for longitudinal research in geroscience

A significant challenge in multi-omic geroscience research is the collection of high quality, fit-for-purpose biospecimens from a diverse and well-characterized study population with sufficient sample size to detect age-related changes in physiological biomarkers. The Dog Aging Project designed the precision cohort to study the mechanisms underlying age-related change in the metabolome, microbiome, and epigenome in companion dogs, an emerging model system for translational geroscience research. One thousand dog-owner pairs were recruited into cohort strata based on life stage, sex, size, and geography. We designed and built a novel implementation of the REDCap electronic data capture system to manage study participants, logistics, and biospecimen and survey data collection in a secure online platform. In collaboration with primary care veterinarians, we collected and processed blood, urine, fecal, and hair samples from 976 dogs. The resulting data include complete blood count, chemistry profile, immunophenotyping by flow cytometry, metabolite quantification, fecal microbiome characterization, epigenomic profile, urinalysis, and associated metadata characterizing sample conditions at collection and during lab processing. The project, which has already begun collecting second- and third-year samples from precision cohort dogs, demonstrates that scientifically useful biospecimens can be collected from a geographically dispersed population through collaboration with private veterinary clinics and downstream labs. The data collection infrastructure developed for the precision cohort can be leveraged for future studies. Most important, the Dog Aging Project is an open data project. We encourage researchers around the world to apply for data access and utilize this rich, constantly growing dataset in their own work.

Integrating ATAC-seq and RNA-seq to reveal the dynamics of chromatin accessibility and gene expression in regulating aril coloration of Taxus mairei

Fruit coloration results from a complex process. Maire yew (Taxus mairei) is an evergreen tree with red, purple, and yellow fruits (arils). While significant progress has been made in understanding pigment biosynthesis in arils, the role of chromatin accessibility in color development remains less well understood. To gain deeper insights into the genetic and epigenetic factors involved, we employed RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin sequencing (ATAC-seq). By integrating the results, we identified 723 differentially expressed genes (DEGs) with chromatin changes in P vs. R, including 312 up- and 411 down-regulated genes. In Y vs. R, we found 159 DEGs, with 97 up- and 62 down-regulated. KEGG enrichment analysis highlighted the flavonoid and carotenoid pathways as major contributors to color variation. When the arils turned purple, the expression levels of C4H, CHS, C3'H, F3'H, F3H, DFR, PSY, PDS, β-OHase, CYP97A3, and LUT1 were significantly up-regulated, while ZDS was down-regulated. The transition to yellow arils was characterized by the up-regulation of F3H, DFR, ZDS, CYP97A3, β-OHase, and LUT1, accompanied by the down-regulation of C4H, CHS, PSY, and PDS. Additionally, 27 transcription factors (TFs) were identified, including MYB, bHLH, and bZIP. These TFs may potentially influence variation in aril color by regulating downstream genes. In total, eight genes were selected for qRT-PCR validation, indicating the reliability of the transcriptome sequencing data. Our results provide in-depth information regarding the coloration of the arils in Maire yew. The study could provide insights for further genetic improvement in Taxus.

Metabolomic Analysis Reveals the Diversity of Defense Metabolites in Nine Cereal Crops

Cereal crops are important staple foods, and their defense metabolites hold significant research importance. In this study, we employed LC-MS-based untargeted and widely-targeted metabolomics to profile the leaf metabolome of nine cereal species, including rice, wheat, maize, barley, sorghum, common oat, foxtail millet, broomcorn millet, and adlay. A total of 9869 features were detected, among them, 1131 were annotated, encompassing 18 classes such as flavonoids, lipids, and alkaloids. Results revealed that 531 metabolites were detected in all species, while each cereal crop possessed 4 to 12 unique metabolites. Focusing on defense metabolites, we identified eight benzoxazinoids uniquely present in maize, wheat, and adlay. Hierarchical clustering based on metabolite abundance divided all metabolites into nine clusters, and subsequent pathway enrichment analysis revealed that the stress-related flavonoid biosynthesis pathway was enriched in multiple clusters. Further analysis showed that four downstream compounds of HBOA (2-hydroxy-1,4-benzoxazin-3-one) in the benzoxazinoid biosynthesis pathway were enriched in maize. Wheat uniquely accumulated the 4'-methylated product of tricin, trimethoxytricetin, whereas adlay accumulated the tricin precursor tricetin in the flavonoid biosynthesis pathway. In summary, this study elucidates the metabolic diversity in defense metabolites among various cereal crops, providing valuable background information for the improvement of stress resistance in cereal crops.

Unraveling key genes and pathways involved in Verticillium wilt resistance by integrative GWAS and transcriptomic approaches in Upland cotton

Verticillium dahliae Kleb, the cause of Verticillium wilt, is a particularly destructive soil-borne vascular disease that affects cotton, resulting in serious decline in fiber quality and causing significant losses in cotton production worldwide. However, the progress in identification of wilt-resistance loci or genes in cotton has been limited, most probably due to the highly complex genetic nature of the trait. Nevertheless, the molecular mechanism behind the Verticillium wilt resistance remains poorly understood. In the present study, we investigated the phenotypic variations in Verticillium tolerance and conducted a genome wide association study (GWAS) among a natural population containing 383 accessions of upland cotton germplasm and performed transcriptomic analysis of cotton genotypes with differential responses to Verticillium wilt. GWAS detected 70 significant SNPs and 116 genes associated with resistance loci in two peak signals on D02 and D11 in E1. The transcriptome analysis identified a total of 2689 and 13289 differentially expressed genes (DEGs) among the Verticillium wilt-tolerant (J46) and wilt-susceptible (J11) genotypes, respectively. The DEGs were predominantly enriched in metabolism, plant hormone signal transduction, phenylpropanoid pathway, MAPK cascade pathway and plant-pathogen interaction pathway in GO and KEGG analyses. The identified DEGs were found to comprise several transcription factor (TF) gene families, primarily including AP2/ERF, ZF, WRKY, NAC and MYB, in addition to pentatricopeptide repeat (PPR) proteins and Resistance (R) genes. Finally, by integrating the two results, 34 candidate genes were found to overlap between GWAS and RNA-seq analyses, associated with Verticillium-wilt resistance, including WRKY, MYB, CYP and RGA. This work contributes to our knowledge of the molecular processes underlying cotton responses to Verticillium wilt, offering crucial insights for additional research into the genes and pathways implicated in these responses and paving the way for developing Verticillium wilt-resistant cotton varieties through accelerated breeding by providing a plethora of candidate genes.

Identification of UDP-glucosyltransferase involved in the biosynthesis of phloridzin in Gossypium hirsutum

Phloridzin has various functions, including antioxidant properties and the treatment of diabetes, and has long been used in pharmaceutical and physiological research. The glycosylation of phloretin is a key step in the biosynthesis of phloridzin. In this study, a genome-wide association study (GWAS) based on phloridzin content was applied, and the key gene GhUGT88F3 for phloridzin-specific biosynthesis was identified in cotton. A single-base deletion in GhUGT88F3 in haplotype I caused a frameshift mutation, leading to premature translation termination and a significant reduction in phloridzin content. Molecular docking revealed important amino acid residues for GhUGT88F3's UDP-glucose transfer activity. Additionally, the transcription factor GhMYB330 was found to positively regulate GhUGT88F3 expression through population transcriptome analysis and LUC experiment. Moreover, phloridzin content was significantly elevated in both GhUGT88F3 and GhMYB330 overexpression transgenic plants. This study expands the diversity of UDP-glucosyltransferases in plants and offers a potential strategy for the sustainable production of bioactive compounds with therapeutic potential.

Mining genomic regions associated with stomatal traits and their candidate genes in bread wheat through genome-wide association study (GWAS)

KEY MESSAGE

112 candidate quantitative trait loci (QTLs) and 53 key candidate genes have been identified as associated with stomatal traits in wheat. These include bHLH, MADS-box transcription factors, and mitogen-activated protein kinases (MAPKs). Stomata is a common feature of the leaf surface of plants and serve as vital conduits for the exchange of gases (primarily CO₂ and water vapor) between plants and the external environment. In this study, a comprehensive genome analysis was conducted by integrating genome-wide association study (GWAS) and genome prediction to identify the genomic regions and candidate genes of stomatal traits associated with drought resistance and water-saving properties in a panel of 184 diverse bread wheat genotypes. There were significant variations on stomatal traits in the wheat panel across different environmental conditions. GWAS was conducted with the genotypic data from the wheat 660 K single-nucleotide polymorphism (SNP) chip, and the stomatal traits conducted across three environments during two growing seasons. The final GWAS identified 112 candidate QTLs that exhibited at least two significant marker-trait associations. Subsequent analysis identified 53 key candidate genes, including 13 bHLH transcription factor, 2 MADS-box transcription factors, and 4 mitogen-activated protein kinase genes, which may be strongly associated with stomatal traits. The application of Bayesian ridge regression for genomic prediction yielded an accuracy rate exceeding 60% for all four stomatal traits in both SNP matrices, with stomatal width achieving a rate in excess of 70%. Additionally, three Kompetitive allele-specific PCR markers were developed and validated, representing a significant advancement in marker-assisted prediction. Overall, these results will contribute to a more comprehensive understanding of wheat stomatal traits and provide a valuable reference for germplasm screening and innovation in wheat germplasm with novel stomatal traits.

Genetic dissection of a major locus SC9.1 conferring seed color in broomcorn millet (Panicum miliaceum)

KEY MESSAGE

A major locus SC9.1 was identified and finely mapped into a 92.68 Kb region, and longmi004412 was identified as the casual gene regulating brown seed color in broomcorn millet. Broomcorn millet is a cereal crop with abundant genetic variations in morphology, agronomy, and yield-related traits. The diversity of seed color is among the most distinctive morphological characteristics. However, genetic determinants governing seed coloration have rarely been reported. Here, the F2 and F3 populations from a cross between Longmi12 and Zhang778 were employed to elucidate the genetic basis of seed color. Statistical analysis conducted on the seed color in F1, F2, and F3 progeny conclusively demonstrated that brown seed color was controlled by a single dominant locus in broomcorn millet. The genetic control locus, SC9.1, was preliminarily located on chromosome 9 in the 32,175,878-44,281,406 bp region through bulked segregant analysis sequencing (BSA-seq). Furthermore, SC9.1 was narrowed down to a 92.68 kb interval harboring 11 genes using fine mapping with 260 recessive individual genotypes. Combined with gene structural variation, the transcriptome profile, and functional comparison, longmi004412 was identified as the causal gene resulting in brown seed color formation in broomcorn millet. In addition, haplotype analysis of the longmi004412 gene in 516 accessions was performed to clarify the types for broomcorn millet seed color. These findings lay the foundation for precise identification of germplasm at the molecular level, molecular-assisted selection breeding, and the application of gene editing technology in broomcorn millet.

Reconciliation of wheat 660K and 90K SNP arrays and their utilization in dough rheological properties of bread wheat

INTRODUCTION

High-density Wheat 660K and 90K SNP arrays are powerful tools for understanding the genetic basis of wheat traits. However, their inconsistantly physical positions that were caused by different versions of Chinese Spring genome during developing arrays are confused and inconvenient for further application.

OBJECTIVE

With the repid development of wheat geonome sequencing, we aim to reconciliate Wheat 660K and 90K SNP arrays in modern cultivar and reveal the genetic basis of dough rheological properties in bread wheat.

METHODS

We refined physical positions of Wheat 660K and 90K SNP arrays in the currently popular wheat cultivar AK58 genome that was released more recently. We next performed genome-wide association studies (GWAS) and linkage analysis to identify important genetic loci related to quality traits using updated and un-updated arrays, respectively.

RESULTS

Refining results showed that 92.3% and 83% of SNPs in the Wheat 660K and 90K SNP arrays were precisely mapped to the AK58 genome, respective. GWAS results by the updated 660K and 90K arrays indicated that 26 intervals composed of 1032 significant SNPs were associated with 9 quality traits in multiple environments. The significant interval for stability time on 1D was narrowed into an 8.4-Mb region using the updated arrays, whereas the interval is 405 Mb using the un-updated arrays. Linkage analysis revealed an important QTL QST.henau-1D.2 for stability time with 1.64 Mb. Integration of GWAS and QTL results narrowed the significant interval into 6.46 Mb containing 35 annotation genes by collinearity analysis. After T-test, gene expression analysis, seven of them are potential candidate genes and thus favorable haplotypes are identified to benefit marker-assisted selection.

CONCLUSION

A reconciliation of Wheat 660K and 90K arrays promote their efficient applications. Important genetic loci and favorable haplotypes identified in this study provided valuable information for wheat quality breeding.

Differentiation of genome-wide DNA methylation between japonica and indica rice

Rice (Oryza sativa L.) subspecies japonica and indica show distinct morphological and genetic differentiation. However, the differences in the genome-wide DNA methylation and its effects on gene expression and metabolic levels between japonica and indica rice remain unclear. In this study, we investigated the genome-wide DNA methylation, transcriptomes and metabolomes of 12 representative japonica and indica rice accessions, to reveal the differentiation between rice subspecies. We detected 83 327 differentially methylated regions (DMRs) and 14 903 DMR-associated genes between two subspecies. Indica rice showed significantly lower levels of the CG, CHG, and CHH methylation compared with japonica rice. Subsequently, we identified 5596 differentially expressed genes between the two subspecies, predominantly enriched in pathways related to carbohydrate and amino acid metabolism. By integrating DNA methylation with transcriptomic data, a significant correlation was established between methylation patterns and the expression level of key agronomic genes in rice. Furthermore, multi-omics analyses reveal that carbohydrate metabolism pathways, especially the tricarboxylic acid (TCA) cycle metabolites, are remarkable differentiation between rice subspecies. These results provide a foundation for future studies in rice domestication and genetic improvement.

Genome-wide association study of metabolic traits in the giant duckweed Spirodela polyrhiza

The exceptionally high growth rate and high flavonoid content make the giant duckweed Spirodela polyrhiza (L.) Schleid. (Arales: Lemnaceae Martinov) an ideal organism for food production and metabolic engineering. To facilitate this, identification of the genetic basis underlying growth and metabolic traits is essential. Here, we analysed growth and content of 42 metabolites in 137 S. polyrhiza genotypes and characterized the genetics underpinning these traits using a genome-wide association (GWA) approach. We found that biomass positively correlated with the content of many free amino acids, including L-glutamine, L-tryptophan, and L-serine, but negatively correlated with specialized metabolites, such as flavonoids. GWA analysis showed that several candidate genes involved in processes such as photosynthesis, protein degradation, and organ development were jointly associated with multiple metabolic traits. The results suggest the above genes are suitable targets for simultaneous optimization of duckweed growth and metabolite levels. This study provides insights into the metabolic diversity of S. polyrhiza and its underlying genetic architecture, paving the way for industrial applications of this plant via targeted breeding or genetic engineering.

A GWAS identified loci and candidate genes associated with fiber quality traits in a new cotton MAGIC population

KEY MESSAGE

GWAS of a new MAGIC population containing alleles from five tetraploid Gossypium species identified novel fiber QTL and confirmed previously identified stable QTL. Identification of loci and underlying genes for fiber quality traits will facilitate genetic improvement in cotton fiber quality. In this research, a genome-wide association study (GWAS) was carried out for fiber quality attributes using a new multi-parent advanced generation inter-cross (MAGIC) population consisting of 372 recombinant inbred lines (RILs). Sixteen parents including 12 exotic germplasm lines derived from five tetraploid Gossypium species and 4 Upland cotton varieties were intercrossed to develop the population. Both RILs and parental lines were evaluated at three locations, College Station, Texas; Las Cruses, New Mexico; and Stoneville, Mississippi, from 2016 through 2023. Fiber length (UHM) had positive correlation with strength (STR) and length uniformity (UNI) and a negative correlation with micronaire (MIC) and elongation (ELO). By combining all the data from all locations, we identified significant SNPs for ELO, UHM, and UNI while STR and MIC were location dependent. These results suggest an important role of genotype by environment interaction in a GWAS of fiber traits. Twenty possible novel fiber QTL were identified: 10 for STR, three for UNI, and seven for MIC. The QTL for ELO (Chr.D04: 53-Mb), UHM (Chr.D11: 24-Mb), and UNI (Chr.D04: 34-Mb) were stable across multiple environments and may be useful for marker-assisted selection to improve fiber quality. For STR, we found candidate genes Gh_A07G1574 and Gh_A07G1581 to be present in the previously identified QTL region (Chr.A07: 77-Mb) on chromosome A07. Identified loci and their corresponding candidate genes will be useful to improve fiber quality via marker-assisted selection in a cotton breeding.

Using VIS-NIR hyperspectral imaging and deep learning for non-destructive high-throughput quantification and visualization of nutrients in wheat grains

High-throughput and low-cost quantification of the nutrient content in crop grains is crucial for food processing and nutritional research. However, traditional methods are time-consuming and destructive. A high-throughput and low-cost method of quantification of wheat nutrients with VIS-NIR (400-1700 nm) hyperspectral imaging is proposed in this study. Stepwise linear regression (SLR) was used to predict hundreds of nutrients accurately (R2 > 0.6); results improved when the hyperspectral data was processed with the first derivative. Knockout materials were also used to verify their practical application value. Various nutrients' characteristic wavelengths were mainly concentrated in the visible regions of 400-500 nm and 900-1000 nm. Finally, we proposed an improved pix2pix conditional generative network model to visualize the nutrients distribution and showed better results compared with the original. This research highlights the potential of hyperspectral technology in high-throughput and non-destructive determination and visualization of grain nutrients with deep learning.

A metabolic roadmap of waxy corn flavor

As well as being a popular vegetable crop worldwide, waxy corn represents an important amylopectin source, but little is known about its breeding history and flavor characteristics. In this study, through comparative-omic analyses between 318 diverse waxy corn and 507 representative field corn inbred lines we revealed that many metabolic pathways and genes exhibited selection characteristics during the breeding history of waxy corn, contributing to the divergence between waxy and field corn. We showed that waxy corn is not only altered in its glutinous property but also its sweetness, aroma, and palatability are all significantly affected. A substantial proportion (43%) of flavor-related metabolites have pleiotropic effects, affecting both flavor and yield characteristics, and 27% of these metabolites are related to antagonistic outcomes on yield and flavor. Furthermore, through multiple concrete examples, we demonstrated how yield and quality are coordinately or antagonistically regulated at the genetic level. In particular, some sweet molecules, such as DIMBOA and raffinose, which do not participate in the starch biosynthesis pathway, were identified as potential targets for breeding a new type of "sweet-waxy" corn. Taken together, our findings shed light on the historical selection of waxy corn and demonstrate the genetic and metabolic basis of waxy corn flavor, collectively providing valuable resources and knowledge for future crop breeding for improved nutritional quality.

Volatilome-based GWAS identifies OsWRKY19 and OsNAC021 as key regulators of rice aroma

Aromatic rice is globally favored for its distinctive scent, which not only increases its nutritional value but also enhances its economic importance. However, apart from 2-acetyl-1-pyrroline (2-AP), the metabolic basis of aroma remains to be clarified, and the genetic basis of the accumulation of fragrance metabolites is largely unknown. In this study, we revealed 2-AP and fatty acid-derived volatiles (FAVs) as key contributors to rice aroma by combining aroma rating with molecular docking. Using a volatilome-based genome-wide association study, we identified two regulatory genes that determine the natural variation of these fragrance metabolites. Genetic and molecular analyses showed that OsWRKY19 not only enhances fragrance by negatively regulating OsBADH2 but also improves agricultural traits in rice. Furthermore, we revealed that OsNAC021 negatively regulates FAV contents via the lipoxygenase pathway, and its knockout resulted in over-accumulation of grain FAVs without a yield penalty. Collectively, our study not only identifies two key regulators of rice aroma but also provides a compelling example about how to deciphering the genetic regulatory mechanisms that underlie rice fragrance, thereby paving the way for the creation of aromatic rice varieties.

Engineering rice genomes towards green super rice

Rice, cultivated for millennia across diverse geographical regions, has witnessed tremendous advancements in recent decades, epitomized by the emergence of Green Super Rice. These efforts aim to address challenges such as climate change, pest and disease threats, and sustainable agriculture. Driven by the advent of multiomics big data, breakthroughs in genomic tools and resources, hybrid rice breeding techniques, and the extensive utilization of green genes, rice genomes are undergoing delicate modifications to produce varieties with high yield, superior quality, enhanced nutrient efficiency, and resilience to pests and environmental stresses, leading to the development of green agriculture in China. Additionally, the utilization of wild relatives and the promotion of genomic breeding approaches have further enriched our understanding of rice improvement. In the future, international efforts to develop next-generation green rice varieties remain both challenging and imperative for the whole community.

Exploring genetics and genomics trends to understand the link between secondary metabolic genes and agronomic traits in cereals under stress

The plant metabolome is considered an important interface between the genome and its phenome, where it plays a significant role in regulating plant growth in response to various environmental cues. A wide array of specialized metabolites is produced by plants, which are essential for mediating environmental interactions and their adaptation. Notably, enhanced accumulation of these specialized metabolites, particularly plant secondary metabolites (PSMs), is a part of the chemical defense response that is directly linked to improved stress tolerance. Therefore, exploring the genetic diversity underlying the immense variation of the secondary metabolite pool could unravel the adaptation mechanisms in plants against different environmental stresses. The post-genomic profiling platforms have enabled the exploration of the link between metabolic diversity and important agronomic traits. The current review focuses on the major achievements and future challenges associated with plant secondary metabolite (PSM) research in graminaceous crops using advanced omics approaches. Given this, we briefly summarize different strategies adopted to explore the genetic diversity and evolution of PSMs in cereal crops. Further, we have discussed the recent technological advancements to integrate multi-omics approaches linking the metabolome diversity with the genome, transcriptome, and proteome of these crops under stress. Combining these data with phenomics (the omics of phenotypes) provides a holistic view of how plants respond to stress. Next, we outlined the genetic manipulation studies performed so far in cereals to engineer secondary metabolic pathways for enhanced stress tolerance. In summary, our review provides new insight into developing genetic and genomic trends in exploring the secondary metabolite diversity in graminaceous crops and discusses how this information can be utilized in designing strategies to generate future stress-resilient crops.

Compositional Analysis of Grape Berries: Mapping the Global Metabolism of Grapes

To characterize the nutrients and bioactive compounds in grape berries and to explore the real cause of the "French paradox" phenomenon, we performed metabolomic analysis of 66 grape varieties worldwide using liquid chromatography-tandem mass spectrometry (LC-MS). A nontargeted metabolomics approach detected a total of 4889 metabolite signals. From these, 964 bioactive and nutrient compounds were identified and quantified, including modified flavonoids, medicinal pentacyclic triterpenoids, vitamins, amino acids, lipids, etc. Interestingly, metabolic variations between varieties are not explained by geography or subspecies but can be significantly distinguished by grapes' color, even after excluding flavonoids and anthocyanins. In our analysis, we found that purple grape varieties had the highest levels of key bioactive components such as flavonoids, pentacyclic triterpenes, and polyphenols, which are thought to have a variety of health benefits such as antioxidant, anti-inflammatory, and antitumor properties, when compared to grapes of other colors. In addition, we found higher levels of vitamins in red and pink grapes, possibly explaining their role in preventing anemia and scurvy and protecting the skin. These findings may be a major factor in the greater health benefits of wines made from purple grapes. Our study provides comprehensive metabolic profiling data of grape berries that may contribute to future research on the French paradox.

Phylogenomic and synteny analysis of BAHD and SCP/SCPL gene families reveal their evolutionary histories in plant specialized metabolism

Plant chemical diversity is largely owing to a number of enzymes which catalyse reactions involved in the assembly, and in the subsequent chemical modifications, of the core structures of major classes of plant specialized metabolites. One such reaction is acylation. With this in mind, to study the deep evolutionary history of BAHD and the serine-carboxypeptidase-like (SCPL) acyltransferase genes, we assembled phylogenomic synteny networks based on a large-scale inference analysis of orthologues across whole-genome sequences of 126 species spanning Stramenopiles and Archaeplastida, including Arabidopsis thaliana, tomato (Solanum lycopersicum) and maize (Zea mays). As such, this study combined the study of genomic location with changes in gene sequences. Our analyses revealed that serine-carboxypeptidase (SCP)/serine-carboxypeptidase-like (SCPL) genes had a deeper evolutionary origin than BAHD genes, which expanded massively on the transition to land and with the development of the vascular system. The two gene families additionally display quite distinct patterns of copy number variation across phylogenies as well as differences in cross-phylogenetic syntenic network components. In unlocking the above observations, our analyses demonstrate the possibilities afforded by modern phylogenomic (syntenic) networks, but also highlight their current limitations, as demonstrated by the inability of phylogenetic methods to separate authentic SCPL acyltransferases from standard SCP peptide hydrolases.This article is part of the theme issue 'The evolution of plant metabolism'.

OsbHLH5 Synergically Regulates Phenolamide and Diterpenoid Phytoalexins Involved in the Defense of Rice Against Pathogens

Rice (Oryza sativa) produces phenolamides and diterpenoids as major phytoalexins. Although the biosynthetic pathways of phenolamides and diterpenoids in plants have been revealed, knowledge of their accumulation regulatory mechanisms remains limited, and, in particular, no co-regulatory factor has been identified to date. Here, using a combined co-expression and evolutionary analysis, we identified the basic helix-loop-helix (bHLH) transcription factor OsbHLH5 as a positive bifunctional regulator of phenolamide and diterpenoid biosynthesis in rice. Metabolomic analysis revealed that OsbHLH5 significantly increased the content of phenolamides (such as feruloyl tryptamine (Fer-Trm) and p-coumaroyl tyramine (Cou-Tyr)) and diterpenoid phytoalexins (such as momilactones A, momilactones B) in the overexpression lines, while their content was reduced in the OsbHLH5 knockout lines. Gene expression and dual-luciferase assays revealed that OsbHLH5 activates phenolamide biosynthetic genes (including putrescine hydroxycinnamoyltransferase 3 (OsPHT3), tyramine hydroxycinnamoyltransferases 1/2 (OsTHT1/2), and tryptamine benzoyltransferase 2 (OsTBT2)) as well as diterpenoid biosynthetic genes (including copalyl diphosphate synthase 4 (OsCPS4) and kaurene synthase-like 4/7/10/11 (OsKSL4/7/10/11)). Furthermore, we have demonstrated that OsbHLH5 is induced by jasmonic acid (JA), while pathogen inoculation assays indicated that the overexpression of OsbHLH5 in transgenic rice plants leads to enhanced resistance to Xanthomonas oryzae pv. oryzae (Xoo). Overall, we have identified a positive regulator of phenolamide and diterpenoid biosynthesis and have demonstrated that biotic stress induces phytoalexin accumulation partly in an OsbHLH5-dependent manner, providing new insights into the metabolic interactions involved in pathogen response and offering valuable gene resources for the development, through genetic improvement, of new rice varieties that are resistant to diseases.

Advances in genomic tools for plant breeding: harnessing DNA molecular markers, genomic selection, and genome editing

Conventional pre-genomics breeding methodologies have significantly improved crop yields since the mid-twentieth century. Genomics provides breeders with advanced tools for whole-genome study, enabling a direct genotype-phenotype analysis. This shift has led to precise and efficient crop development through genomics-based approaches, including molecular markers, genomic selection, and genome editing. Molecular markers, such as SNPs, are crucial for identifying genomic regions linked to important traits, enhancing breeding accuracy and efficiency. Genomic resources viz. genetic markers, reference genomes, sequence and protein databases, transcriptomes, and gene expression profiles, are vital in plant breeding and aid in the identification of key traits, understanding genetic diversity, assist in genomic mapping, support marker-assisted selection and speeding up breeding programs. Advanced techniques like CRISPR/Cas9 allow precise gene modification, accelerating breeding processes. Key techniques like Genome-Wide Association study (GWAS), Marker-Assisted Selection (MAS), and Genomic Selection (GS) enable precise trait selection and prediction of breeding outcomes, improving crop yield, disease resistance, and stress tolerance. These tools are handy for complex traits influenced by multiple genes and environmental factors. This paper explores new genomic technologies like molecular markers, genomic selection, and genome editing for plant breeding showcasing their impact on developing new plant varieties.

Uncovering fruit flavor and genetic diversity across diploid wild Fragaria species via comparative metabolomics profiling

Wild Fragaria resources exhibit extensive genetic diversity and desirable edible traits, such as high soluble solid content and flavor compounds. However, specific metabolites in different wild strawberry fruits remain unknown. In this study, we characterized 1008 metabolites covering 11 subclasses among 13 wild diploid resources representing eight species, including F. vesca, F. nilgerrensis, F. viridis, F. nubicola, F. pentaphylla, F. mandschurica, F. chinensis, and F. emeiensis. Fifteen potential metabolite biomarkers were identified to distinguish fruit flavors among the 13 diploid wild Fragaria accessions. A total of nine distinct modules were employed to explore key metabolites related to fruit quality through weighted gene co-expression module analysis, with significant enrichment in amino acid biosynthesis pathway. Notably, the identified significantly different key metabolites highlighted the close association of amino acids, sugars, and anthocyanins with flavor formation. These findings offer valuable resources for improving fruit quality through metabolome-assisted breeding.

Cellulose synthase-like OsCSLD4: a key regulator of agronomic traits, disease resistance, and metabolic indices in rice

KEY MESSAGE

Cellulose synthase-like OsCSLD4 plays a pivotal role in regulating diverse agronomic traits, enhancing resistance against bacterial leaf blight, and modulating metabolite indices based on the multi-omics analysis in rice. To delve deeper into this complex network between agronomic traits and metabolites in rice, we have compiled a dataset encompassing genome, phenome, and metabolome, including 524 diverse accessions, 11 agronomic traits, and 841 metabolites, enabling us to pinpoint eight hotspots through GWAS. We later discovered four distinct metabolite categories, encompassing 15 metabolites that are concurrently present on the QTL qC12.1, associated with leaf angle of flag and spikelet length, and finally focused the cellulose synthase-like OsCSLD4, which was pinpointed through a rigorous process encompassing sequence variation, haplotype, ATAC, and differential expression across diverse tissues. Compared to the wild type, csld4 exhibited significant reductions in the plant height, flag leaf length, leaf width, spikelet length, 1000-grain weight, grain width, grain thickness, fertility, yield per plant, and bacterial blight resistance. However, there were significant increase in tiller numbers, degree of leaf rolling, flowering period, growth period, grain length, and empty kernel rate. Furthermore, the content of four polyphenol metabolites, excluding metabolite N-feruloyltyramine (mr1268), notably rose, whereas the levels of the other three polyphenol metabolites, smiglaside C (mr1498), 4-coumaric acid (mr1622), and smiglaside A (mr1925) decreased significantly in mutant csld4. The content of amino acid L-tyramine (mr1446) exhibited a notable increase, whereas the alkaloid trigonelline (mr1188) displayed a substantial decrease among the mutants. This study offered a comprehensive multi-omics perspective to analyze the genetic mechanism of OsCSLD4, and breeders can potentially enhance rice's yield, bacterial leaf blight resistance, and metabolite content, leading to more sustainable and profitable rice production.

Integrating evolutionary genomics of forest trees to inform future tree breeding amid rapid climate change

Global climate change is leading to rapid and drastic shifts in environmental conditions, posing threats to biodiversity and nearly all life forms worldwide. Forest trees serve as foundational components of terrestrial ecosystems and play a crucial and leading role in combating and mitigating the adverse effects of extreme climate events, despite their own vulnerability to these threats. Therefore, understanding and monitoring how natural forests respond to rapid climate change is a key priority for biodiversity conservation. Recent progress in evolutionary genomics, driven primarily by cutting-edge multi-omics technologies, offers powerful new tools to address several key issues. These include precise delineation of species and evolutionary units, inference of past evolutionary histories and demographic fluctuations, identification of environmentally adaptive variants, and measurement of genetic load levels. As the urgency to deal with more extreme environmental stresses grows, understanding the genomics of evolutionary history, local adaptation, future responses to climate change, and conservation and restoration of natural forest trees will be critical for research at the nexus of global change, population genomics, and conservation biology. In this review, we explore the application of evolutionary genomics to assess the effects of global climate change using multi-omics approaches and discuss the outlook for breeding of climate-adapted trees.

Genome-wide association study and genotypic variation for the major tocopherol content in rice grain

Rice tocopherols, vitamin E compounds with antioxidant activity, play essential roles in human health. Even though the key genes involved in vitamin E biosynthetic pathways have been identified in plants, the genetic architecture of vitamin E content in rice grain remains unclear. A genome-wide association study (GWAS) on 179 genotypically diverse rice accessions with 34,323 SNP markers was conducted to detect QTLs that define total and α- tocopherol contents in rice grains. Total and α-tocopherol contents had a strong positive correlation and varied greatly across the accessions, ranging from 0.230-31.76 and 0.011-30.83 (μg/g), respectively. A total of 13 QTLs were identified, which were spread across five of the rice chromosomes. Among the 13 QTLs, 11 were considered major with phenotypic variation explained (PVE) greater than 10%. Twelve transcription factor (TF) genes, one microprotein (miP), and a transposon were found to be associated with the QTLs with putative roles in controlling tocopherol contents. Moreover, intracellular transport proteins, ABC transporters, nonaspanins, and SNARE, were identified as associated genes on chromosomes 1 and 8. In the vicinity of seven QTLs, protein kinases were identified as key signaling factors. Haplotype analysis revealed the QTLs qAlph1.1, qTot1.1, qAlph2.1, qAlph6.1, qTot6.1, and qTot8.3 to have significant haplogroups. Quantitative RT-PCR validated the expression direction and magnitude of WRKY39 (Os02g0265200), PIP5Ks (Os08g0450800), and MADS59 (Os06g0347700) in defining the major tocopherol contents. This study provides insights for ongoing biofortification efforts to breed and/or engineer vitamin E and antioxidant levels in rice and other cereals.

Amino acid permease OsAAP12 negatively regulates rice tillers and grain yield by transporting specific amino acids to affect nitrogen and cytokinin pathways

Amino acids are necessary nutrients for the growth of Oryza sativa (rice), which can be mediated by amino acid transporter; however, our understanding of these transporters is still limited. This study found that the expression levels of amino acid permease gene OsAAP12 differed between indica and japonica rice. Altered expression of OsAAP12 negatively regulated tillering and yield in transgenic rice lines. Subcellular localization revealed that OsAAP12 was primarily localized to the plasma membrane. Moreover, it was indicated that OsAAP12 transported polar neutral amino acids asparagine (Asn), threonine (Thr), and serine (Ser) through experiments involving yeast heterologous complementation, fluorescence amino acid uptake, and amino acid content determination. Additionally, exogenous application of amino acids Asn, Thr, and Ser suppressed axillary buds outgrowth in OsAAP12 overexpression lines compared with wild-type ZH11. Conversely, the opposite trend was observed in CRISPR mutant lines. RNA-seq analysis showed that the expression patterns of genes involved in the nitrogen and cytokinin pathways were generally altered in OsAAP12 modified lines. Hormone assays indicated that OsAAP12 mutant lines accumulated cytokinins in the basal part of rice, whereas overexpression lines had the opposite effect. In summary, CRISPR mutant of OsAAP12 boosted rice tillering and grain yield by coordinating the content of amino acids and cytokinins, which has potential application value in high-yield rice breeding.

The AP2/ERF transcription factor ORA59 regulates ethylene-induced phytoalexin synthesis through modulation of an acyltransferase gene expression

The gaseous ethylene (ET) and the oxylipin-derived jasmonic acid (JA) in plants jointly regulate an arsenal of pathogen responsive genes involved in defending against necrotrophic pathogens. The APETALA2 (AP2)/ETHYLENE RESPONSE FACTOR (ERF) transcription factor ORA59 is a major positive regulator of the ET/JA-mediated defense pathway in Arabidopsis thaliana. The Arabidopsis agmatine coumaroyltransferase (AtACT) catalyzes the formation of hydroxycinnamic acid amides (HCAAs) which are effective toxic antimicrobial substances known as phytoalexins and play an important role in plant defense response. However, induction and regulation of AtACT gene expression and HCAAs synthesis in plants remain less understood. Through gene coexpression network analysis, we identified a list of GCC-box cis-element containing genes that were coexpressed with ORA59 under diverse biotic stress conditions and might be potential downstream targets of this AP2/ERF-domain transcription factor. Particularly, ORA59 directly binds to AtACT gene promoter via the GCC-boxes and activates AtACT gene expression. The ET precursor 1-aminocyclopropane-1-carboxylic acid (ACC)-treatment significantly induces AtACT gene expression. Both ORA59 and members of the class II TGA transcription factors are indispensable for ACC-induced AtACT expression. Interestingly, the expression of AtACT is also subject to the signaling crosstalk of the salicylic acid- and ET/JA-mediated defense response pathways. In addition, we found that genes of the phenylpropanoid metabolism pathway were specifically induced by Botrytis cinerea. Taking together, these evidence suggest that the ET/JA signaling pathway activate the expression of AtACT to increase antimicrobial HCAAs production through the transcription factor ORA59 in response to the infection of necrotrophic plant pathogens.

Variation in WIDTH OF LEAF AND GRAIN contributes to grain and leaf size by controlling LARGE2 stability in rice

Grain and flag leaf size are two important agronomic traits that influence grain yield in rice (Oryza sativa). Many quantitative trait loci (QTLs) and genes that regulate these traits individually have been identified, however, few QTLs and genes that simultaneously control these two traits have been identified. In this study, we conducted a genome-wide association analysis in rice and detected a major locus, WIDTH OF LEAF AND GRAIN (WLG), that was associated with both grain and flag leaf width. WLG encodes a RING-domain E3 ubiquitin ligase. WLGhap.B, which possesses five single nucleotide polymophysim (SNP) variations compared to WLGhap.A, encodes a protein with enhanced ubiquitination activity that confers increased rice leaf width and grain size, whereas mutation of WLG leads to narrower leaves and smaller grains. Both WLGhap.A and WLGhap.B interact with LARGE2, a HETC-type E3 ligase, however, WLGhap.B exhibits stronger interaction with LARGE2, thus higher ubiquitination activity toward LARGE2 compared with WLGhap.A. Lysine1021 is crucial for the ubiquitination of LARGE2 by WLG. Loss-of-function of LARGE2 in wlg-1 phenocopies large2-c in grain and leaf width, suggesting that WLG acts upstream of LARGE2. These findings reveal the genetic and molecular mechanism by which the WLG-LARGE2 module mediates grain and leaf size in rice and suggest the potential of WLGhap.B in improving rice yield.

Dissecting the genetic basis of UV-B responsive metabolites in rice

BACKGROUND

UV-B, an important environmental factor, has been shown to affect the yield and quality of rice (Oryza sativa) worldwide. However, the molecular mechanisms underlying the response to UV-B stress remain elusive in rice.

RESULTS

We perform comprehensive metabolic profiling of leaves from 160 diverse rice accessions under UV-B and normal light conditions using a widely targeted metabolomics approach. Our results reveal substantial differences in metabolite accumulation between the two major rice subspecies indica and japonica, especially after UV-B treatment, implying the possible role and mechanism of metabolome changes in subspecies differentiation and the stress response. We next conduct a transcriptome analysis from four representative rice varieties under UV-B stress, revealing genes from amino acid and flavonoid pathways involved in the UV-B response. We further perform a metabolite-based genome-wide association study (mGWAS), which reveals 3307 distinct loci under UV-B stress. Identification and functional validation of candidate genes show that OsMYB44 regulates tryptamine accumulation to mediate UV-B tolerance, while OsUVR8 interacts with OsMYB110 to promote flavonoid accumulation and UV-B tolerance in a coordinated manner. Additionally, haplotype analysis suggests that natural variation of OsUVR8groupA contributes to UV-B resistance in rice.

CONCLUSIONS

Our study reveals the complex biochemical and genetic foundations that govern the metabolite dynamics underlying the response, tolerance, and adaptive strategies of rice to UV-B stress. These findings provide new insights into the biochemical and genetic basis of the metabolome underlying the crop response, tolerance, and adaptation to UV-B stress.

Variation in a Poaceae-conserved fatty acid metabolic gene cluster controls rice yield by regulating male fertility

A wide variety of metabolic gene clusters exist in eukaryotic genomes, but fatty acid metabolic gene clusters have not been discovered. Here, combining with metabolic and phenotypic genome-wide association studies, we identify a major locus containing a six-gene fatty acid metabolic gene cluster on chromosome 3 (FGC3) that controls the cutin monomer hydroxymonoacylglycerols (HMGs) contents and rice yield, possibly through variation in the transcription of FGC3 members. We show that HMGs are sequentially synthesized in the endoplasmic reticulum by OsFAR2, OsKCS11, OsGPAT6, OsCYP704B2 and subsequently transported to the apoplast by OsABCG22 and OsLTPL82. Mutation of FGC3 members reduces HMGs, leading to defective male reproductive development and a significant decrease in yield. OsMADS6 and OsMADS17 directly regulate FGC3 and thus influence male reproduction and yield. FGC3 is conserved in Poaceae and likely formed prior to the divergence of Pharus latifolius. The eukaryotic fatty acid and plant primary metabolic gene cluster we identified show a significant impact on the origin and evolution of Poaceae and has potential for application in hybrid crop breeding.

Transcriptome-referenced association study provides insights into the regulation of oil and fatty acid biosynthesis in Torreya grandis kernel

INTRODUCTION

Torreya grandis is a gymnosperm belonging to Taxodiaceae. As an economically important tree, its kernels are edible and rich in oil with high unsaturated fatty acids, such as sciadonic acid. However, the kernels from different T. grandis landraces exhibit fatty acid and oil content variations.

OBJECTIVES

As a gymnosperm, does T. grandis have special regulation mechanisms for oil biosynthesis? The aim of this study was to dissect the genetic architecture of fatty acid and oil content and the underlying mechanism in T. grandis.

METHODS

We constructed a high integrity reference sequence of expressed regions of the genome in T. grandis and performed transcriptome-referenced association study (TRAS) for 10 fatty acid and oil traits of kernels in the 170 diverse T. grandis landraces. To confirm the TRAS result, we performed functional validation and molecular biology experiments for oil significantly associated genes.

RESULTS

We identified 41 SNPs from 34 transcripts significantly associated with 7 traits by TRAS (-log10 (P) greater than 6.0). Results showed that LOB domain-containing protein 40 (LBD40) and surfeit locus protein 1 (SURF1) may be indirectly involved in the regulation of oil and sciadonic acid biosynthesis, respectively. Moreover, overexpression of TgLBD40 significantly increased seed oil content. The nonsynonymous variant in the TgLBD40 coding region discovered by TRAS could alter the oil content in plants. Pearson's correlation analysis and dual-luciferase assay indicated that TgLBD40 positively enhanced oil accumulation by affecting oil biosynthesis pathway genes, such as TgDGAT1.

CONCLUSION

Our study provides new insights into the genetic basis of oil biosynthesis in T. grandis and demonstrates that integrating RNA sequencing and TRAS is a powerful strategy to perform association study independent of a reference genome for dissecting important traits in T. grandis.

Beyond pathways: Accelerated flavonoids candidate identification and novel exploration of enzymatic properties using combined mapping populations of wheat

Although forward-genetics-metabolomics methods such as mGWAS and mQTL have proven effective in providing myriad loci affecting metabolite contents, they are somehow constrained by their respective constitutional flaws such as the hidden population structure for GWAS and insufficient recombinant rate for QTL. Here, the combination of mGWAS and mQTL was performed, conveying an improved statistical power to investigate the flavonoid pathways in common wheat. A total of 941 and 289 loci were, respectively, generated from mGWAS and mQTL, within which 13 of them were co-mapped using both approaches. Subsequently, the mGWAS or mQTL outputs alone and their combination were, respectively, utilized to delineate the metabolic routes. Using this approach, we identified two MYB transcription factor encoding genes and five structural genes, and the flavonoid pathway in wheat was accordingly updated. Moreover, we have discovered some rare-activity-exhibiting flavonoid glycosyl- and methyl-transferases, which may possess unique biological significance, and harnessing these novel catalytic capabilities provides potentially new breeding directions. Collectively, we propose our survey illustrates that the forward-genetics-metabolomics approaches including multiple populations with high density markers could be more frequently applied for delineating metabolic pathways in common wheat, which will ultimately contribute to metabolomics-assisted wheat crop improvement.

Tryptophan regulates sorghum root growth and enhances low nitrogen tolerance

Over evolutionary time, plants have developed sophisticated regulatory mechanisms to adapt to fluctuating nitrogen (N) environments, ensuring that their growth is balanced with their responses to N stress. This study explored the potential of L-tryptophan (Trp) in regulating sorghum root growth under conditions of N limitation. Here, two distinct sorghum genotypes (low-N tolerance 398B and low-N sensitive CS3541) were utilized for investigating effect of low-N stress on root morphology and conducting a comparative transcriptomics analysis. Our foundings indicated that 398B exhibited longer roots, greater root dry weights, and a higher Trp content compared to CS3541 under low-N conditions. Furthermore, transcriptome analysis revealed substantial differences in gene expression profiles related to Trp pathway and carbon (C) and N metabolism pathways between the two genotypes. Additional experiments were conducted to assess the effects of exogenous Trp treatment on the interplay between sorghum root growth and low-N tolerance. Our observations showed that Trp-treated plants developed longer root and had elevated levels of Trp and IAA under low-N conditons. Concurrently, these plants demonstrated stronger physiological activities in C and N metabolism when subjected to low-N stress. These results underscored the pivotal role of Trp on root growth and low-N stress responses by balancing IAA levels and C and N metabolism. This study not only deepens our understanding of how plants maintain growth plasticity during environmental stress but also provides valuable insights into the availability of amino acid in crops, which could be instrumental in developing strategies for promoting crop resilience to N deficiency.

Combined genome and transcriptome analysis of elite fiber quality in Gossypium barbadense

Gossypium barbadense, which is one of several species of cotton, is well known for its superior fiber quality. However, the genetic basis of its high-quality fiber remains largely unexplored. Here, we resequenced 269 G. barbadense accessions. Phylogenetic structure analysis showed that the set of accessions was clustered into 3 groups: G1 and G2 mainly included modern cultivars from Xinjiang, China, and G3 was related to widely introduced accessions in different regions worldwide. A genome-wide association study of 5 fiber quality traits across multiple field environments identified a total of 512 qtls (main-effect QTLs) and 94 qtlEs (QTL-by-environment interactions) related to fiber quality, of which 292 qtls and 57 qtlEs colocated with previous studies. We extracted the genes located in these loci and performed expression comparison, local association analysis, and introgression segment identification. The results showed that high expression of hormone-related genes during fiber development, introgressions from Gossypium hirsutum, and the recombination of domesticated elite allelic variation were 3 major contributors to improve the fiber quality of G. barbadense. In total, 839 candidate genes with encoding region variations associated with elite fiber quality were mined. We confirmed that haplotype GB_D03G0092H traced to G. hirsutum introgression, with a 1-bp deletion leading to a frameshift mutation compared with GB_D03G0092B, significantly improved fiber quality. GB_D03G0092H is localized in the plasma membrane, while GB_D03G0092B is in both the nucleus and plasma membrane. Overexpression of GB_D03G0092H in Arabidopsis (Arabidopsis thaliana) significantly improved the elongation of longitudinal cells. Our study systematically reveals the genetic basis of the superior fiber quality of G. barbadense and provides elite segments and gene resources for breeding high-quality cotton cultivars.

Identification and characterization of camptothecin tailoring enzymes in Nothapodytes tomentosa

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.

Blocking of amino acid transporter OsAAP7 promoted tillering and yield by determining basic and neutral amino acids accumulation in rice

BACKGROUND

Amino acids are not only the main form of N in rice, but also are vital for its growth and development. These processes are facilitated by amino acid transporters within the plant. Despite their significance, only a few AAP amino acid transporters have been reported.

RESULTS

In this study, we observed that there were differences in the expression of amino acid transporter OsAAP7 among 521 wild cultivated rice varieties, and it directly negatively correlated with tillering and grain yield per plant. We revealed that OsAAP7 protein was localized to the endoplasmic reticulum and had absorption and transport affinity for amino acids such as phenylalanine (Phe), lysine (Lys), leucine (Leu), and arginine (Arg) using subcellular localization, yeast substrate testing, fluorescent amino acid uptake, and amino acid content determination. Further hydroponic studies showed that exogenous application of amino acids Phe, Lys and Arg inhibited the growth of axillary buds in the overexpression lines, and promoted the elongation of axillary buds in the mutant lines. Finally, RNA-seq analysis showed that the expression patterns of genes related to nitrogen, auxin and cytokinin pathways were changed in axillary buds of OsAAP7 transgenic plants.

CONCLUSIONS

This study revealed the gene function of OsAAP7, and found that blocking of amino acid transporter OsAAP7 with CRISPR/Cas9 technology promoted tillering and yield by determining basic and neutral amino acids accumulation in rice.

OsLSC6 Regulates Leaf Sheath Color and Cold Tolerance in Rice Revealed by Metabolite Genome Wide Association Study

Plant metabolites including anthocyanins play an important role in the growth of plants, as well as in regulating biotic and abiotic stress responses to the environment. Here we report comprehensive profiling of 3315 metabolites and a further metabolic-based genome-wide association study (mGWAS) based on 292,485 SNPs obtained from 311 rice accessions, including 160 wild and 151 cultivars. We identified hundreds of common variants affecting a large number of secondary metabolites with large effects at high throughput. Finally, we identified a novel gene namely OsLSC6 (Oryza sativa leaf sheath color 6), which encoded a UDP 3-O-glucosyltransferase and involved in the anthocyanin biosynthesis of Cyanidin-3-Galc (sd1825) responsible for leaf sheath color, and resulted in significant different accumulation of sd1825 between wild (purple) and cultivars (green). The results of knockout transgenic experiments showed that OsLSC6 regulated the biosynthesis and accumulation of sd1825, controlled the purple leaf sheath. Our further research revealed that OsLSC6 also confers resistance to cold stress during the seedling stage in rice. And we identified that a SNP in OsLSC6 was responsible for the leaf sheath color and chilling tolerance, supporting the importance of OsLSC6 in plant adaption. Our study could not only demonstrate that OsLSC6 is a vital regulator during anthocyanin biosynthesis and abiotic stress responses, but also provide a powerful complementary tool based on metabolites-to-genes analysis by mGWAS for functional gene identification andpromising candidate in future rice breeding and improvement.

Integrated transcriptome and metabolome analysis reveals the regulation of phlorizin synthesis in Lithocarpus polystachyus under nitrogen fertilization

BACKGROUND

Nitrogen (N) is essential for plant growth and development. In Lithocarpus polystachyus Rehd., a species known for its medicinal and food value, phlorizin is the major bioactive compound with pharmacological activity. Research has revealed a positive correlation between plant nitrogen (N) content and phlorizin synthesis in this species. However, no study has analyzed the effect of N fertilization on phlorizin content and elucidated the molecular mechanisms underlying phlorizin synthesis in L. polystachyus.

RESULTS

A comparison of the L. polystachyus plants grown without (0 mg/plant) and with N fertilization (25, 75, 125, 175, 225, and 275 mg/plant) revealed that 75 mg N/plant fertilization resulted in the greatest seedling height, ground diameter, crown width, and total phlorizin content. Subsequent analysis of the leaves using ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) detected 150 metabolites, including 42 flavonoids, that were differentially accumulated between the plants grown without and with 75 mg/plant N fertilization. Transcriptomic analysis of the L. polystachyus plants via RNA sequencing revealed 162 genes involved in flavonoid biosynthesis, among which 53 significantly differed between the N-treated and untreated plants. Fertilization (75 mg N/plant) specifically upregulated the expression of the genes phenylalanine ammonia-lyase (PAL), 4-coumarate-CoA ligase (4CL), and phlorizin synthase (PGT1) but downregulated the expression of trans-cinnamate 4-monooxygenase (C4H), shikimate O-hydroxycinnamoyltransferase (HCT), and chalcone isomerase (CHI), which are related to phlorizin synthesis. Finally, an integrated analysis of the transcriptome and metabolome revealed that the increase in phlorizin after N fertilization was consistent with the upregulation of phlorizin biosynthetic genes. Quantitative real-time PCR (qRT‒PCR) was used to validate the RNA sequencing data. Thus, our results indicated that N fertilization increased phlorizin metabolism in L. polystachyus by regulating the expression levels of the PAL, PGT1, 5-O-(4-coumaroyl)-D-quinate 3'-monooxygenase (C3'H), C4H, and HCT genes.

CONCLUSIONS

Our results demonstrated that the addition of 75 mg/plant N to L. polystachyus significantly promoted the accumulation of flavonoids, including phlorizin, and the expression of flavonoid synthesis-related genes. Under these conditions, the genes PAL, 4CL, and PGT1 were positively correlated with phlorizin accumulation, while C4H, CHI, and HCT were negatively correlated with phlorizin accumulation. Therefore, we speculate that PAL, 4CL, and PGT1 participate in the phlorizin pathway under an optimal N environment, regulating phlorizin biosynthesis. These findings provide a basis for improving plant bioactive constituents and serve as a reference for further pharmacological studies.

Comparative metabolomics combined with genome sequencing provides insights into novel wolfberry-specific metabolites and their formation mechanisms

Wolfberry (Lycium, of the family Solanaceae) has special nutritional benefits due to its valuable metabolites. Here, 16 wolfberry-specific metabolites were identified by comparing the metabolome of wolfberry with those of six species, including maize, rice, wheat, soybean, tomato and grape. The copy numbers of the riboflavin and phenyllactate degradation genes riboflavin kinase (RFK) and phenyllactate UDP-glycosyltransferase (UGT1) were lower in wolfberry than in other species, while the copy number of the phenyllactate synthesis gene hydroxyphenyl-pyruvate reductase (HPPR) was higher in wolfberry, suggesting that the copy number variation of these genes among species may be the main reason for the specific accumulation of riboflavin and phenyllactate in wolfberry. Moreover, the metabolome-based neighbor-joining tree revealed distinct clustering of monocots and dicots, suggesting that metabolites could reflect the evolutionary relationship among those species. Taken together, we identified 16 specific metabolites in wolfberry and provided new insight into the accumulation mechanism of species-specific metabolites at the genomic level.

New Insights into the Genetic Basis of Lysine Accumulation in Rice Revealed by Multi-Model GWAS

Lysine is an essential amino acid that cannot be synthesized in humans. Rice is a global staple food for humans but has a rather low lysine content. Identification of the quantitative trait nucleotides (QTNs) and genes underlying lysine content is crucial to increase lysine accumulation. In this study, five grain and three leaf lysine content datasets and 4,630,367 single nucleotide polymorphisms (SNPs) of 387 rice accessions were used to perform a genome-wide association study (GWAS) by ten statistical models. A total of 248 and 71 common QTNs associated with grain/leaf lysine content were identified. The accuracy of genomic selection/prediction RR-BLUP models was up to 0.85, and the significant correlation between the number of favorable alleles per accession and lysine content was up to 0.71, which validated the reliability and additive effects of these QTNs. Several key genes were uncovered for fine-tuning lysine accumulation. Additionally, 20 and 30 QTN-by-environment interactions (QEIs) were detected in grains/leaves. The QEI-sf0111954416 candidate gene LOC_Os01g21380 putatively accounted for gene-by-environment interaction was identified in grains. These findings suggested the application of multi-model GWAS facilitates a better understanding of lysine accumulation in rice. The identified QTNs and genes hold the potential for lysine-rich rice with a normal phenotype.

Combining association with linkage mapping to dissect the phenolamides metabolism of the maize kernel

Phenolamides are important secondary metabolites in plant species. They play important roles in plant defense responses against pathogens and insect herbivores, protection against UV irradiation and floral induction and development. However, the accumulation and variation in phenolamides content in diverse maize lines and the genes responsible for their biosynthesis remain largely unknown. Here, we combined genetic mapping, protein regulatory network and bioinformatics analysis to further enhance the understanding of maize phenolamides biosynthesis. Sixteen phenolamides were identified in multiple populations, and they were all significantly correlated with one or several of 19 phenotypic traits. By linkage mapping, 58, 58, 39 and 67 QTLs, with an average of 3.9, 3.6, 3.6 and 4.2 QTLs for each trait were mapped in BBE1, BBE2, ZYE1 and ZYE2, explaining 9.47%, 10.78%, 9.51% and 11.40% phenotypic variation for each QTL on average, respectively. By GWAS, 39 and 36 significant loci were detected in two different environments, 3.3 and 2.8 loci for each trait, explaining 10.00% and 9.97% phenotypic variation for each locus on average, respectively. Totally, 58 unique candidate genes were identified, 31% of them encoding enzymes involved in amine and derivative metabolic processes. Gene Ontology term analysis of the 358 protein-protein interrelated genes revealed significant enrichment in terms relating to cellular nitrogen metabolism, amine metabolism. GRMZM2G066142, GRMZM2G066049, GRMZM2G165390 and GRMZM2G159587 were further validated involvement in phenolamides biosynthesis. Our results provide insights into the genetic basis of phenolamides biosynthesis in maize kernels, understanding phenolamides biosynthesis and its nutritional content and ability to withstand biotic and abiotic stress.

Genetic dissection of ten photosynthesis-related traits based on InDel- and SNP-GWAS in soybean

KEY MESSAGE

A total of 416 InDels and 112 SNPs were significantly associated with soybean photosynthesis-related traits. GmIWS1 and GmCDC48 might be related to chlorophyll fluorescence and gas-exchange parameters, respectively. Photosynthesis is one of the main factors determining crop yield. A better understanding of the genetic architecture for photosynthesis is of great significance for soybean yield improvement. Our previous studies identified 5,410,112 single nucleotide polymorphisms (SNPs) from the resequencing data of 219 natural soybean accessions. Here, we identified 634,106 insertions and deletions (InDels) from these 219 accessions and used these InDel variations to perform principal component and linkage disequilibrium analysis of this population. The genome-wide association study (GWAS) were conducted on six chlorophyll fluorescence parameters (chlorophyll content, light energy absorbed per reaction center, quantum yield for electron transport, probability that a trapped exciton moves an electron into the electron transport chain beyond primary quinone acceptor, maximum quantum yield of photosystem II primary photochemistry in the dark-adapted state, performance index on absorption basis) and four gas-exchange parameters (intercellular carbon dioxide concentration, stomatal conductance, net photosynthesis rate, transpiration rate) and revealed 416 significant InDels and 112 significant SNPs. Based on GWAS results, GmIWS1 (encoding a transcription elongation factor) and GmCDC48 (encoding a cell division cycle protein) with the highest expression in the mapping region were determined as the candidate genes responsible for chlorophyll fluorescence and gas-exchange parameters, respectively. Further identification of favorable haplotypes with higher photosynthesis, seed weight and seed yield were carried out for GmIWS1 and GmCDC48. Overall, this study revealed the natural variations and candidate genes underlying the photosynthesis-related traits based on abundant phenotypic and genetic data, providing valuable insights into the genetic mechanisms controlling photosynthesis and yield in soybean.

Natural variation in OsMYB8 confers diurnal floret opening time divergence between indica and japonica subspecies

The inter-subspecific indica-japonica hybrid rice confer potential higher yield than the widely used indica-indica intra-subspecific hybrid rice. Nevertheless, the utilization of this strong heterosis is currently hindered by asynchronous diurnal floret opening time (DFOT) of indica and japonica parental lines. Here, we identify OsMYB8 as a key regulator of rice DFOT. OsMYB8 induces the transcription of JA-Ile synthetase OsJAR1, thereby regulating the expression of genes related to cell osmolality and cell wall remodeling in lodicules to promote floret opening. Natural variations of OsMYB8 promoter contribute to its differential expression, thus differential transcription of OsJAR1 and accumulation of JA-Ile in lodicules of indica and japonica subspecies. Furthermore, introgression of the indica haplotype of OsMYB8 into japonica effectively promotes DFOT in japonica. Our findings reveal an OsMYB8-OsJAR1 module that regulates differential DFOT in indica and japonica, and provide a strategy for breeding early DFOT japonica to facilitate breeding of indica-japonica hybrids.