Plant Public RNA-seq Database: a comprehensive online database for expression analysis of ~45 000 plant public RNA-Seq libraries

Spatial accumulation of lignin monomers and cellulose underlying stalk strength in maize

Lodging largely affects yield, quality and mechanical harvesting of maize. Stalk strength is one of the major factors that affect maize lodging. Although plant cell wall components including lignin and cellulose were known to be associated with stalk strength and lodging resistance, spatial accumulation of specific lignin monomers and cellulose in different tissues and their association with stalk strength in maize was not clearly understood. In this study, we found that both G and S lignin monomers accumulate highest in root, stem rind and leaf vein. Consistently, most lignin biosynthetic genes were expressed higher in root and stem than in other tissues. However, cellulose appears to be lowest in root. There are only mild changes of G lignin and cellulose in different internodes. Instead, we noticed a dramatic decrease of S-lignin accumulation and lignin biosynthetic gene expression in 2nd to 4th internodes wherein stem breakage usually occurs, thereby revealing a few candidate lignin biosynthetic genes associated with stalk strength. Moreover, stalk strength is positively correlated with G, S lignin, and cellulose, but negatively correlated with S/G ratio based on data of maize lines with high or low stalk strength. Loss-of-function of a caffeic acid o-methyltransferase (COMT), which is involved in S lignin biosynthesis, in the maize bm3 mutant, leads to lower stalk strength. Our data collectively suggest that stalk strength is determined by tissue-specific accumulation of lignin monomers and cellulose, and manipulation of the cell wall components by genetic engineering is vital to improve maize stalk strength and lodging resistance.

A repertoire of intronic lariat RNAs reveals tissue-specific regulation and target mimicry potential in plants

Lariat RNA is concomitantly produced by excised intron during RNA splicing, which is usually debranched by DBR1, an RNA debranching enzyme. However, increasing evidence showed that some lariat RNA could escape debranching. Little is known about how and why these lariat RNAs could be retained. By comparing the atlas of lariat RNAs between the non-dividing cell (mature pollen) and three actively dividing tissues (young shoot apex, young seeds, and young roots), we identified hundreds to thousands of lariat RNA naturally retained in each tissue, and the incidence of lariat RNA retention is much less in shoot apex while much more in pollen. Many lariat RNAs derived from the same intron or different lariat RNAs from the same pre-mRNA could be retained in one tissue while degraded in the other tissues. By deciphering lariat RNA sequences, we identified an AG-rich (RAAAAVAAAR) motif and a UC-rich (UCUCUYUCUC) motif for pollen-specific and the other three tissues-retained lariat RNAs, respectively. Reconstitution of the pollen-specific AG-rich motif indeed enhanced lariat RNA retention in plants. Biologically, hundreds of lariat RNAs harbored miRNA binding sites, and dual-luciferase reporter assay showed that these natural lariat RNAs had the potential to protect expression of miRNA target genes. Collectively, our results uncover that selective retention of lariat RNA is an actively regulatory process, and provide new insights into understanding how lariat RNA metabolism may impact miRNA activity.

QTL mapping and genomic selection of stem and branch diameter in soybean (Glycine max L.)

Introduction

Soybean stem diameter (SD) and branch diameter (BD) are closely related traits, and genetic clarification of SD and BD is crucial for soybean breeding.

Methods

SD and BD were genetically analyzed by a population of 363 RIL derived from the cross between Zhongdou41 (ZD41) and ZYD02878 using restricted two-stage multi-locus genome-wide association, inclusive composite interval mapping, and three-variance component multi-locus random SNP effect mixed linear modeling. Then candidate genes of major QTLs were selected and genetic selection model of SD and BD were constructed respectively.

Results and discussion

The results showed that SD and BD were significantly correlated (r = 0.74, P < 0.001). A total of 93 and 84 unique quantitative trait loci (QTL) were detected for SD and BD, respectively by three different methods. There were two and ten major QTLs for SD and BD, respectively, with phenotypic variance explained (PVE) by more than 10%. Within these loci, seven genes involved in the regulation of phytohormones (IAA and GA) and cell proliferation and showing extensive expression of shoot apical meristematic genes were selected as candidate genes. Genomic selection (GS) analysis showed that the trait-associated markers identified in this study reached 0.47-0.73 in terms of prediction accuracy, which was enhanced by 6.56-23.69% compared with genome-wide markers. These results clarify the genetic basis of SD and BD, which laid solid foundation in regulation gene cloning, and GS models constructed could be potentially applied in future breeding programs.

Comprehensive integration of single-cell transcriptomic data illuminates the regulatory network architecture of plant cell fate specification

Plant morphogenesis relies on precise gene expression programs at the proper time and position which is orchestrated by transcription factors (TFs) in intricate regulatory networks in a cell-type specific manner. Here we introduced a comprehensive single-cell transcriptomic atlas of Arabidopsis seedlings. This atlas is the result of meticulous integration of 63 previously published scRNA-seq datasets, addressing batch effects and conserving biological variance. This integration spans a broad spectrum of tissues, including both below- and above-ground parts. Utilizing a rigorous approach for cell type annotation, we identified 47 distinct cell types or states, largely expanding our current view of plant cell compositions. We systematically constructed cell-type specific gene regulatory networks and uncovered key regulators that act in a coordinated manner to control cell-type specific gene expression. Taken together, our study not only offers extensive plant cell atlas exploration that serves as a valuable resource, but also provides molecular insights into gene-regulatory programs that varies from different cell types.

Plant genome information facilitates plant functional genomics

MAIN CONCLUSION

In this review, we give an overview of plant sequencing efforts and how this impacts plant functional genomics research. Plant genome sequence information greatly facilitates the studies of plant biology, functional genomics, evolution of genomes and genes, domestication processes, phylogenetic relationships, among many others. More than two decades of sequencing efforts have boosted the number of available sequenced plant genomes. The first plant genome, of Arabidopsis, was published in the year 2000 and currently, 4604 plant genomes from 1482 plant species have been published. Various large sequence initiatives are running, which are planning to produce tens of thousands of sequenced plant genomes in the near future. In this review, we give an overview on the status of sequenced plant genomes and on the use of genome information in different research areas.

PEO: Plant Expression Omnibus - a comparative transcriptomic database for 103 Archaeplastida

The Plant Expression Omnibus (PEO) is a web application that provides biologists with access to gene expression insights across over 100 plant species, ~60 000 manually annotated RNA-seq samples, and more than 4 million genes. The tool allows users to explore the expression patterns of genes across different organs, identify organ-specific genes, and discover top co-expressed genes for any gene of interest. PEO also provides functional annotations for each gene, allowing for the identification of genetic modules and pathways. PEO is designed to facilitate comparative kingdom-wide gene expression analysis and provide a valuable resource for plant biology research. We provide two case studies to demonstrate the utility of PEO in identifying candidate genes in pollen coat biosynthesis in Arabidopsis and investigating the biosynthetic pathway components of capsaicin in Capsicum annuum. The database is freely available at https://expression.plant.tools/.

Evolution of plant cell-type-specific cis -regulatory elements

Cis -regulatory elements (CREs) are critical in regulating gene expression, and yet our understanding of CRE evolution remains a challenge. Here, we constructed a comprehensive single-cell atlas of chromatin accessibility in Oryza sativa , integrating data from 104,029 nuclei representing 128 discrete cell states across nine distinct organs. We used comparative genomics to compare cell-type resolved chromatin accessibility between O. sativa and 57,552 nuclei from four additional grass species ( Zea mays, Sorghum bicolor, Panicum miliaceum , and Urochloa fusca ). Accessible chromatin regions (ACRs) had different levels of conservation depending on the degree of cell-type specificity. We found a complex relationship between ACRs with conserved noncoding sequences, cell-type specificity, conservation, and tissue-specific switching. Additionally, we found that epidermal ACRs were less conserved compared to other cell types, potentially indicating that more rapid regulatory evolution has occurred in the L1 epidermal layer of these species. Finally, we identified and characterized a conserved subset of ACRs that overlapped the repressive histone modification H3K27me3, implicating them as potentially critical silencer CREs maintained by evolution. Collectively, this comparative genomics approach highlights the dynamics of cell-type-specific CRE evolution in plants.

Molecular evolution and interaction of 14-3-3 proteins with H+-ATPases in plant abiotic stresses

Environmental stresses severely affect plant growth and crop productivity. Regulated by 14-3-3 proteins (14-3-3s), H+-ATPases (AHAs) are important proton pumps that can induce diverse secondary transport via channels and co-transporters for the abiotic stress response of plants. Many studies demonstrated the roles of 14-3-3s and AHAs in coordinating the processes of plant growth, phytohormone signaling, and stress responses. However, the molecular evolution of 14-3-3s and AHAs has not been summarized in parallel with evolutionary insights across multiple plant species. Here, we comprehensively review the roles of 14-3-3s and AHAs in cell signaling to enhance plant responses to diverse environmental stresses. We analyzed the molecular evolution of key proteins and functional domains that are associated with 14-3-3s and AHAs in plant growth and hormone signaling. The results revealed evolution, duplication, contraction, and expansion of 14-3-3s and AHAs in green plants. We also discussed the stress-specific expression of those 14-3-3and AHA genes in a eudicotyledon (Arabidopsis thaliana), a monocotyledon (Hordeum vulgare), and a moss (Physcomitrium patens) under abiotic stresses. We propose that 14-3-3s and AHAs respond to abiotic stresses through many important targets and signaling components of phytohormones, which could be promising to improve plant tolerance to single or multiple environmental stresses.

Uncovering Novel QTLs and Candidate Genes for Salt Tolerance at the Bud Burst Stage in Rice through Genome-Wide Association Study

Salt stress is one of the most important factors limiting rice growth and yield increase. Salt tolerance of rice at the bud burst (STB) stage determines whether germinated seeds can grow normally under salt stress, which is very important for direct seeding. However, reports on quantitative trait loci (QTLs) and candidate genes for STB in rice are very limited. In this study, a natural population of 130 indica and 81 japonica rice accessions was used to identify STB-related QTLs and candidate genes using a genome-wide association study (GWAS). Nine QTLs, including five for relative shoot length (RSL), two for relative root length (RRL), and two for relative root number (RRN), were identified. Five of these STB-related QTLs are located at the same site as the characterized salt tolerance genes, such as OsMDH1, OsSRFP1, and OsCDPK7. However, an important QTL related to RSL, qRSL1-2, has not been previously identified and was detected on chromosome 1. The candidate region for qRSL1-2 was identified by linkage disequilibrium analysis, 18 genes were found to have altered expression levels under salt stress through the RNA-seq database, and 10 of them were found to be highly expressed in the shoot. It was also found that, eight candidate genes (LOC_Os01g62980, LOC_Os01g63190, LOC_Os01g63230, LOC_Os01g63280, LOC_Os01g63400, LOC_Os01g63460, and LOC_Os01g63580) for qRSL1-2 carry different haplotypes between indica and japonica rice, which exactly corresponds to the significant difference in RSL values between indica and japonica rice in this study. Most of the accessions with elite haplotypes were indica rice, which had higher RSL values. These genes with indica-japonica specific haplotypes were identified as candidate genes. Rice accessions with elite haplotypes could be used as important resources for direct seeding. This study also provides new insights into the genetic mechanism of STB.

Navigating rice seedling cold resilience: QTL mapping in two inbred line populations and the search for genes

Due to global climate change resulting in extreme temperature fluctuations, it becomes increasingly necessary to explore the natural genetic variation in model crops such as rice to facilitate the breeding of climate-resilient cultivars. To uncover genomic regions in rice involved in managing cold stress tolerance responses and to identify associated cold tolerance genes, two inbred line populations developed from crosses between cold-tolerant and cold-sensitive parents were used for quantitative trait locus (QTL) mapping of two traits: degree of membrane damage after 1 week of cold exposure quantified as percent electrolyte leakage (EL) and percent low-temperature seedling survivability (LTSS) after 1 week of recovery growth. This revealed four EL QTL and 12 LTSS QTL, all overlapping with larger QTL regions previously uncovered by genome-wide association study (GWAS) mapping approaches. Within the QTL regions, 25 cold-tolerant candidate genes were identified based on genomic differences between the cold-tolerant and cold-sensitive parents. Of those genes, 20% coded for receptor-like kinases potentially involved in signal transduction of cold tolerance responses; 16% coded for transcription factors or factors potentially involved in regulating cold tolerance response effector genes; and 64% coded for protein chaperons or enzymes potentially serving as cold tolerance effector proteins. Most of the 25 genes were cold temperature regulated and had deleterious nucleotide variants in the cold-sensitive parent, which might contribute to its cold-sensitive phenotype.

Integrating omics databases for enhanced crop breeding

Crop plant breeding involves selecting and developing new plant varieties with desirable traits such as increased yield, improved disease resistance, and enhanced nutritional value. With the development of high-throughput technologies, such as genomics, transcriptomics, and metabolomics, crop breeding has entered a new era. However, to effectively use these technologies, integration of multi-omics data from different databases is required. Integration of omics data provides a comprehensive understanding of the biological processes underlying plant traits and their interactions. This review highlights the importance of integrating omics databases in crop plant breeding, discusses available omics data and databases, describes integration challenges, and highlights recent developments and potential benefits. Taken together, the integration of omics databases is a critical step towards enhancing crop plant breeding and improving global food security.

NRTPredictor: identifying rice root cell state in single-cell RNA-seq via ensemble learning

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) measurements of gene expression show great promise for studying the cellular heterogeneity of rice roots. How precisely annotating cell identity is a major unresolved problem in plant scRNA-seq analysis due to the inherent high dimensionality and sparsity.

RESULTS

To address this challenge, we present NRTPredictor, an ensemble-learning system, to predict rice root cell stage and mine biomarkers through complete model interpretability. The performance of NRTPredictor was evaluated using a test dataset, with 98.01% accuracy and 95.45% recall. With the power of interpretability provided by NRTPredictor, our model recognizes 110 marker genes partially involved in phenylpropanoid biosynthesis. Expression patterns of rice root could be mapped by the above-mentioned candidate genes, showing the superiority of NRTPredictor. Integrated analysis of scRNA and bulk RNA-seq data revealed aberrant expression of Epidermis cell subpopulations in flooding, Pi, and salt stresses.

CONCLUSION

Taken together, our results demonstrate that NRTPredictor is a useful tool for automated prediction of rice root cell stage and provides a valuable resource for deciphering the rice root cellular heterogeneity and the molecular mechanisms of flooding, Pi, and salt stresses. Based on the proposed model, a free webserver has been established, which is available at https://www.cgris.net/nrtp .

Deciphering the evolution and complexity of wheat germplasm from a genomic perspective

Bread wheat provides an essential fraction of the daily calorific intake for humanity. Due to its huge and complex genome, progress in studying on the wheat genome is substantially trailed behind those of the other two major crops, rice and maize, for at least a decade. With rapid advances in genome assembling and reduced cost of high-throughput sequencing, emerging de novo genome assemblies of wheat and whole-genome sequencing data are leading to a paradigm shift in wheat research. Here, we review recent progress in dissecting the complex genome and germplasm evolution of wheat since the release of the first high-quality wheat genome. New insights have been gained in the evolution of wheat germplasm during domestication and modern breeding progress, genomic variations at multiple scales contributing to the diversity of wheat germplasm, and complex transcriptional and epigenetic regulations of functional genes in polyploid wheat. Genomics databases and bioinformatics tools meeting the urgent needs of wheat genomics research are also summarized. The ever-increasing omics data, along with advanced tools and well-structured databases, are expected to accelerate deciphering the germplasm and gene resources in wheat for future breeding advances.

Long noncoding RNA study: Genome-wide approaches

Long noncoding RNAs (lncRNAs) have been confirmed to play a crucial role in various biological processes across several species. Though many efforts have been devoted to the expansion of the lncRNAs landscape, much about lncRNAs is still unknown due to their great complexity. The development of high-throughput technologies and the constantly improved bioinformatic methods have resulted in a rapid expansion of lncRNA research and relevant databases. In this review, we introduced genome-wide research of lncRNAs in three parts: (i) novel lncRNA identification by high-throughput sequencing and computational pipelines; (ii) functional characterization of lncRNAs by expression atlas profiling, genome-scale screening, and the research of cancer-related lncRNAs; (iii) mechanism research by large-scale experimental technologies and computational analysis. Besides, primary experimental methods and bioinformatic pipelines related to these three parts are summarized. This review aimed to provide a comprehensive and systemic overview of lncRNA genome-wide research strategies and indicate a genome-wide lncRNA research system.

ZmEREB92 plays a negative role in seed germination by regulating ethylene signaling and starch mobilization in maize

Rapid and uniform seed germination is required for modern cropping system. Thus, it is important to optimize germination performance through breeding strategies in maize, in which identification for key regulators is needed. Here, we characterized an AP2/ERF transcription factor, ZmEREB92, as a negative regulator of seed germination in maize. Enhanced germination in ereb92 mutants is contributed by elevated ethylene signaling and starch degradation. Consistently, an ethylene signaling gene ZmEIL7 and an α-amylase gene ZmAMYa2 are identified as direct targets repressed by ZmEREB92. OsERF74, the rice ortholog of ZmEREB92, shows conserved function in negatively regulating seed germination in rice. Importantly, this orthologous gene pair is likely experienced convergently selection during maize and rice domestication. Besides, mutation of ZmEREB92 and OsERF74 both lead to enhanced germination under cold condition, suggesting their regulation on seed germination might be coupled with temperature sensitivity. Collectively, our findings uncovered the ZmEREB92-mediated regulatory mechanism of seed germination in maize and provide breeding targets for maize and rice to optimize seed germination performance towards changing climates.

A new family of proteins is required for tethering of Casparian strip membrane domain and nutrient homoeostasis in rice

Cell-cell junctions are essential for multicellular organisms to maintain nutrient homoeostasis. A plant-type tight junction, the Casparian strip (CS)-Casparian strip membrane domain (CSD) that seals the paracellular space between adjacent endodermal cells, has been known for more than one hundred years. However, the molecular basis of this structure remains unknown. Here we report that a new family of proteins containing a glycine/alanine/proline-rich domain, a lectin domain and a secretory signal peptide (GAPLESS) mediates tethering of the plasma membrane to the CS in rice. The GAPLESS proteins are specifically localized in the CS of root endodermal cells, and loss of their functions results in a disabled cell-cell junction and disrupted nutrient homoeostasis. The GAPLESS protein forms a tight complex with OsCASP1 in the plasma membrane, thereby mediating the CS-CSD junction. This study provides valuable insights into the junctional complex of plant endodermal cells, shedding light on our understanding of nutrient homoeostasis in crops and the cell junctions of eukaryotes.

JGI Plant Gene Atlas: an updateable transcriptome resource to improve functional gene descriptions across the plant kingdom

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

Spontaneous movement of a retrotransposon generated genic dominant male sterility providing a useful tool for rice breeding

Male sterility in plants provides valuable breeding tools in germplasm innovation and hybrid crop production. However, genetic resources for dominant genic male sterility, which hold great promise to facilitate breeding processes, are extremely rare in natural germplasm. Here we characterized the Sanming Dominant Genic Male Sterility in rice and identified the gene SDGMS using a map-based cloning approach. We found that spontaneous movement of a 1978-bp long terminal repeat (LTR) retrotransposon into the promoter region of the SDGMS gene activates its expression in anther tapetum, which causes abnormal programmed cell death of tapetal cells resulting in dominant male sterility. SDGMS encodes a ribosome inactivating protein showing N-glycosidase activity. The activation of SDGMS triggers transcription reprogramming of genes responsive to biotic stress leading to a hypersensitive response which causes sterility. The results demonstrate that an ectopic gene activation by transposon movement can give birth to a novel trait which enriches phenotypic diversity with practical utility.

ZmMS1/ZmLBD30-orchestrated transcriptional regulatory networks precisely control pollen exine development

Because of its significance for plant male fertility and, hence, direct impact on crop yield, pollen exine development has inspired decades of scientific inquiry. However, the molecular mechanism underlying exine formation and thickness remains elusive. In this study, we identified that a previously unrecognized repressor, ZmMS1/ZmLBD30, controls proper pollen exine development in maize. Using an ms1 mutant with aberrantly thickened exine, we cloned a male-sterility gene, ZmMs1, which encodes a tapetum-specific lateral organ boundary domain transcription factor, ZmLBD30. We showed that ZmMs1/ZmLBD30 is initially turned on by a transcriptional activation cascade of ZmbHLH51-ZmMYB84-ZmMS7, and then it serves as a repressor to shut down this cascade via feedback repression to ensure timely tapetal degeneration and proper level of exine. This activation-feedback repression loop regulating male fertility is conserved in maize and sorghum, and similar regulatory mechanism may also exist in other flowering plants such as rice and Arabidopsis. Collectively, these findings reveal a novel regulatory mechanism of pollen exine development by which a long-sought master repressor of upstream activators prevents excessive exine formation.

Unraveling genomic regions and candidate genes for multiple disease resistance in upland cotton using meta-QTL analysis

Verticillium wilt (VW), Fusarium wilt (FW) and Root-knot nematode (RKN) are the main diseases affecting cotton production. However, many reported quantitative trait loci (QTLs) for cotton resistance have not been used for agricultural practices because of inconsistencies in the cotton genetic background. The integration of existing cotton genetic resources can facilitate the discovery of important genomic regions and candidate genes involved in disease resistance. Here, an improved and comprehensive meta-QTL analysis was conducted on 487 disease resistant QTLs from 31 studies in the last two decades. A consensus linkage map with genetic overall length of 3006.59 cM containing 8650 markers was constructed. A total of 28 Meta-QTLs (MQTLs) were discovered, among which nine MQTLs were identified as related to resistance to multiple diseases. Candidate genes were predicted based on public transcriptome data and enriched in pathways related to disease resistance. This study used a method based on the integration of Meta-QTL, known genes and transcriptomics to reveal major genomic regions and putative candidate genes for resistance to multiple diseases, providing a new basis for marker-assisted selection of high disease resistance in cotton breeding.

CropGF: a comprehensive visual platform for crop gene family mining and analysis

A gene family refers to a group of genes that share a common ancestry and encode proteins or RNA molecules with similar functions or structural features. Gene families play a crucial role in determining the traits of plants and can be utilized to develop new crop varieties. Therefore, a comprehensive database of gene family is significant for gaining deep insight into crops. To address this need, we have developed CropGF (https://bis.zju.edu.cn/cropgf), a comprehensive visual platform that encompasses six important crops (rice, wheat, maize, barley, sorghum and foxtail millet) and one model plant (Arabidopsis), as well as genomics, transcriptomics and proteomics data for gene family mining and analysis, covering a total of 314 611 genes and 4399 types of domains. CropGF provides a versatile search system that allows for the identification of gene families and their members in a single crop or multiple crops. Users can customize their search based on gene family domains and/or homology using keywords or BLAST. To enhance usability, we have collected the corresponding ID information from various public databases for both genes and domains. Furthermore, CropGF comprises numerous downstream analysis modules, such as ka/ks analysis, phylogenetic tree construction, subcellular localization analysis and more. These visually-displayed modules provide intuitive insights into gene expression patterns, gene family expansion and functional relationships across different molecular levels and different species. We believe that CropGF will be a valuable resource for deep mining and analysis in future studies of crop gene families. Database URL https://bis.zju.edu.cn/cropgf.

The maize ZmVPS23-like protein relocates the nucleotide-binding leucine-rich repeat protein Rp1-D21 to endosomes and suppresses the defense response

Plants often utilize nucleotide-binding leucine-rich repeat (NLR) proteins to perceive pathogen infections and trigger a hypersensitive response (HR). The endosomal sorting complex required for transport (ESCRT) machinery is a conserved multisubunit complex that is essential for the biogenesis of multivesicular bodies and cargo protein sorting. VPS23 is a key component of ESCRT-I and plays important roles in plant development and abiotic stresses. ZmVPS23L, a homolog of VPS23-like in maize (Zea mays), was previously identified as a candidate gene in modulating HR mediated by the autoactive NLR protein Rp1-D21 in different maize populations. Here, we demonstrate that ZmVPS23L suppresses Rp1-D21-mediated HR in maize and Nicotiana benthamiana. Variation in the suppressive effect of HR by different ZmVPS23L alleles was correlated with variation in their expression levels. ZmVPS23 also suppressed Rp1-D21-mediated HR. ZmVPS23L and ZmVPS23 predominantly localized to endosomes, and they physically interacted with the coiled-coil domain of Rp1-D21 and mediated the relocation of Rp1-D21 from the nucleo-cytoplasm to endosomes. In summary, we demonstrate that ZmVPS23L and ZmVPS23 are negative regulators of Rp1-D21-mediated HR, likely by sequestrating Rp1-D21 in endosomes via physical interaction. Our findings reveal the role of ESCRT components in controlling plant NLR-mediated defense responses.

Genome-wide association and epistasis studies reveal the genetic basis of saline-alkali tolerance at the germination stage in rice

Introduction

Saline-alkali stress is one of the main abiotic factors limiting rice production worldwide. With the widespread use of rice direct seeding technology, it has become increasingly important to improve rice saline-alkali tolerance at the germination stage.

Methods

To understand the genetic basis of saline-alkali tolerance and facilitate breeding efforts for developing saline-alkali tolerant rice varieties, the genetic basis of rice saline-alkali tolerance was dissected by phenotyping seven germination-related traits of 736 diverse rice accessions under the saline-alkali stress and control conditions using genome-wide association and epistasis analysis (GWAES).

Results

Totally, 165 main-effect quantitative trait nucleotides (QTNs) and 124 additional epistatic QTNs were identified as significantly associated with saline-alkali tolerance, which explained a significant portion of the total phenotypic variation of the saline-alkali tolerance traits in the 736 rice accessions. Most of these QTNs were located in genomic regions either harboring saline-alkali tolerance QTNs or known genes for saline-alkali tolerance reported previously. Epistasis as an important genetic basis of rice saline-alkali tolerance was validated by genomic best linear unbiased prediction in which inclusion of both main-effect and epistatic QTNs showed a consistently better prediction accuracy than either main-effect or epistatic QTNs alone. Candidate genes for two pairs of important epistatic QTNs were suggested based on combined evidence from the high-resolution mapping plus their reported molecular functions. The first pair included a glycosyltransferase gene LOC_Os02g51900 (UGT85E1) and an E3 ligase gene LOC_Os04g01490 (OsSIRP4), while the second pair comprised an ethylene-responsive transcriptional factor, AP59 (LOC_Os02g43790), and a Bcl-2-associated athanogene gene, OsBAG1 (LOC_Os09g35630) for salt tolerance. Detailed haplotype analyses at both gene promoter and CDS regions of these candidate genes for important QTNs identified favorable haplotype combinations with large effects on saline-alkali tolerance, which can be used to improve rice saline-alkali tolerance by selective introgression.

Discussion

Our findings provided saline-alkali tolerant germplasm resources and valuable genetic information to be used in future functional genomic and breeding efforts of rice saline-alkali tolerance at the germination stage.

Genome-Wide Identification and Expression Pattern Analysis of Dirigent Members in the Genus Oryza

Dirigent (DIR) members have been shown to play essential roles in plant growth, development and adaptation to environmental changes. However, to date, there has been no systematic analysis of the DIR members in the genus Oryza. Here, 420 genes were identified from nine rice species to have the conserved DIR domain. Importantly, the cultivated rice species Oryza sativa has more DIR family members than the wild rice species. DIR proteins in rice could be classified into six subfamilies based on phylogeny analysis. Gene duplication event analysis suggests that whole genome/segmental duplication and tandem duplication are the primary drivers for DIR genes' evolution in Oryza, while tandem duplication is the main mechanism of gene family expansion in the DIR-b/d and DIR-c subfamilies. Analysis of the RNA sequencing data indicates that OsjDIR genes respond to a wide range of environmental factors, and most OsjDIR genes have a high expression level in roots. Qualitative reverse transcription PCR assays confirmed the responsiveness of OsjDIR genes to the undersupply of mineral elements, the excess of heavy metals and the infection of Rhizoctonia solani. Furthermore, there exist extensive interactions between DIR family members. Taken together, our results shed light on and provide a research foundation for the further exploration of DIR genes in rice.

BarleyExpDB: an integrative gene expression database for barley

BACKGROUND

RNA-sequencing (RNA-seq) has been widely used to study the dynamic expression patterns of transcribed genes, which can lead to new biological insights. However, processing and analyzing these huge amounts of histological data remains a great challenge for wet labs and field researchers who lack bioinformatics experience and computational resources.

RESULTS

We present BarleyExpDB, an easy-to-operate, free, and web-accessible database that integrates transcriptional profiles of barley at different growth and developmental stages, tissues, and stress conditions, as well as differential expression of mutants and populations to build a platform for barley expression and visualization. The expression of a gene of interest can be easily queried by searching by known gene ID or sequence similarity. Expression data can be displayed as a heat map, along with functional descriptions as well as Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Proteins Families Database, and Simple Modular Architecture Research Tool annotations.

CONCLUSIONS

BarleyExpDB will serve as a valuable resource for the barley research community to leverage the vast publicly available RNA-seq datasets for functional genomics research and crop molecular breeding.

LMTdb: A comprehensive transcriptome database for climate-resilient, nutritionally rich little millet (Panicum sumatrense)

Little millet (Panicum sumatrense) a native of Chhattisgarh, belongs to the minor millet group and is primarily known as a climate-resilient and nutritionally rich crop. However, due to the lack of enough Omic studies on the crop, the scientific community has largely remained unaware of the potential of this crop, resulting in less scope for its utilization in crop improvement programs. Looking at global warming, erratic climate change, nutritional security, and limited genetic information available, the Little Millet Transcriptome Database (LMT^db) (https://igkv.ac.in/xenom/index.aspx) was conceptualized upon completion of the transcriptome sequencing of little millet with the aim of deciphering the genetic signatures of this largely unknown crop. The database was developed with the view of providing information about the most comprehensive part of the genome, the 'Transcriptome'. The database includes transcriptome sequence information, functional annotation, microsatellite markers, DEGs, and pathway information. The database is a freely available resource that provides breeders and scientists a portal to search, browse, and query data to facilitate functional and applied Omic studies in millet crops.

Plant Intron-Splicing Efficiency Database (PISE): exploring splicing of ∼1,650,000 introns in Arabidopsis, maize, rice, and soybean from ∼57,000 public RNA-seq libraries

Intron retention is the most common alternative splicing event in plants and plays a crucial role in the responses of plants to environmental signals. Despite a large number of RNA-seq libraries from different treatments and genetic mutants stored in public domains, a resource for querying the intron-splicing ratio of individual intron is still required. Here, we established the first-ever large-scale splicing efficiency database in any organism. Our database includes over 57,000 plant public RNA-seq libraries, comprising 25,283 from Arabidopsis, 17,789 from maize, 10,710 from rice, and 3,974 from soybean, and covers a total of 1.6 million introns in these four species. In addition, we manually curated and annotated all the mutant- and treatment-related libraries as well as their matched controls included in our library collection, and added graphics to display intron-splicing efficiency across various tissues, developmental stages, and stress-related conditions. The result is a large collection of 3,313 treatment conditions and 3,594 genetic mutants for discovering differentially regulated splicing efficiency. Our online database can be accessed at https://plantintron.com/ .

From Classical to Modern Computational Approaches to Identify Key Genetic Regulatory Components in Plant Biology

The selection of plant genotypes with improved productivity and tolerance to environmental constraints has always been a major concern in plant breeding. Classical approaches based on the generation of variability and selection of better phenotypes from large variant collections have improved their efficacy and processivity due to the implementation of molecular biology techniques, particularly genomics, Next Generation Sequencing and other omics such as proteomics and metabolomics. In this regard, the identification of interesting variants before they develop the phenotype trait of interest with molecular markers has advanced the breeding process of new varieties. Moreover, the correlation of phenotype or biochemical traits with gene expression or protein abundance has boosted the identification of potential new regulators of the traits of interest, using a relatively low number of variants. These important breakthrough technologies, built on top of classical approaches, will be improved in the future by including the spatial variable, allowing the identification of gene(s) involved in key processes at the tissue and cell levels.

Characterization and evolutionary analysis of phosphate starvation response genes in wheat and other major gramineous plants

Developing cultivars with improved Pi use efficiency is essential for the sustainability of agriculture as well as the environment. Phosphate starvation response (PHR) regulators have not yet been systematically studied in wheat. This study provides the detailed characteristics of PHRs in hexaploid wheat as well as other major gramineous plants at the genome-wide level. The identified PHR proteins were divided into six subfamilies through phylogeny analysis, and a total of 63 paralogous TaPHR pairs were designated as arising from duplication events, with strong purifying selection. The promoters of TaPHRs were identified as stations for many transcription factors. Protein-protein interaction network and gene ontology enrichment analysis indicated a core biological process of cellular response to phosphate starvation. The three-dimensional structures of core PHR proteins showed a high phylogenetic relationship, but amino acid deletions in core protein domains may cause functional differentiation between rice and wheat. TaPHR3 could interact with TaSPX1 and TaSPX5 proteins, which is regarded as a novel interaction mode. Under different Pi gradient treatments, TaPHRs showed low inducible expression patterns among all subfamilies. Our study is the first to comprehensively clarify the basic properties of TaPHR proteins and might accumulate basic data for improving grain yield and environmental homeostasis.

Comparative analysis of physiological variations and genetic architecture for cold stress response in soybean germplasm

Soybean (Glycine max L.) is susceptible to low temperatures. Increasing lines of evidence indicate that abiotic stress-responsive genes are involved in plant low-temperature stress response. However, the involvement of photosynthesis, antioxidants and metabolites genes in low temperature response is largely unexplored in Soybean. In the current study, a genetic panel of diverse soybean varieties was analyzed for photosynthesis, chlorophyll fluorescence and leaf injury parameters under cold stress and control conditions. This helps us to identify cold tolerant (V100) and cold sensitive (V45) varieties. The V100 variety outperformed for antioxidant enzymes activities and relative expression of photosynthesis (Glyma.08G204800.1, Glyma.12G232000.1), GmSOD (GmSOD01, GmSOD08), GmPOD (GmPOD29, GmPOD47), trehalose (GmTPS01, GmTPS13) and cold marker genes (DREB1E, DREB1D, SCOF1) than V45 under cold stress. Upon cold stress, the V100 variety showed reduced accumulation of H₂O₂ and MDA levels and subsequently showed lower leaf injury compared to V45. Together, our results uncovered new avenues for identifying cold tolerant soybean varieties from a large panel. Additionally, we identified the role of antioxidants, osmo-protectants and their posttranscriptional regulators miRNAs such as miR319, miR394, miR397, and miR398 in Soybean cold stress tolerance.

Large-scale analyses of angiosperm Flowering Locus T genes reveal duplication and functional divergence in monocots

FLOWERING LOCUS T (FT) are well-known key genes for initiating flowering in plants. Delineating the evolutionary history and functional diversity of FT genes is important for understanding the diversification of flowering time and how plants adapt to the changing surroundings. We performed a comprehensive phylogenetic analysis of FT genes in 47 sequenced flowering plants and the 1,000 Plant Transcriptomes (1KP) database with a focus on monocots, especially cereals. We revealed the evolutionary history of FT genes. The FT genes in monocots can be divided into three clades (I, II, and III), whereas only one monophyletic group was detected in early angiosperms, magnoliids, and eudicots. Multiple rounds of whole-genome duplications (WGD) events followed by gene retention contributed to the expansion and variation of FT genes in monocots. Amino acid sites in the clade II and III genes were preferentially under high positive selection, and some sites located in vital domain regions are known to change functions when mutated. Clade II and clade III genes exhibited high variability in important regions and functional divergence compared with clade I genes; thus, clade I is more conserved than clade II and III. Genes in clade I displayed higher expression levels in studied organs and tissues than the clade II and III genes. The co-expression modules showed that some of the FT genes might have experienced neofunctionalization and subfunctionalization, such as the acquisition of environmental resistance. Overall, FT genes in monocots might form three clades by the ancient gene duplication, and each clade was subsequently subjected to different selection pressures and amino acid substitutions, which eventually led to different expression patterns and functional diversification. Our study provides a global picture of FT genes' evolution in monocots, paving a road for investigating FT genes' function in future.

Genome-Wide Association Studies of Seven Root Traits in Soybean (Glycine max L.) Landraces

Soybean [Glycine max (L.) Merr.], an important oilseed crop, is a low-cost source of protein and oil. In Southeast Asia and Africa, soybeans are widely cultivated for use as traditional food and feed and industrial purposes. Given the ongoing changes in global climate, developing crops that are resistant to climatic extremes and produce viable yields under predicted climatic conditions will be essential in the coming decades. To develop such crops, it will be necessary to gain a thorough understanding of the genetic basis of agronomic and plant root traits. As plant roots generally lie beneath the soil surface, detailed observations and phenotyping throughout plant development present several challenges, and thus the associated traits have tended to be ignored in genomics studies. In this study, we phenotyped 357 soybean landraces at the early vegetative (V2) growth stages and used a 180 K single-nucleotide polymorphism (SNP) soybean array in a genome-wide association study (GWAS) conducted to determine the phenotypic relationships among root traits, elucidate the genetic bases, and identify significant SNPs associated with root trait-controlling genomic regions/loci. A total of 112 significant SNP loci/regions were detected for seven root traits, and we identified 55 putative candidate genes considered to be the most promising. Our findings in this study indicate that a combined approach based on SNP array and GWAS analyses can be applied to unravel the genetic basis of complex root traits in soybean, and may provide an alternative high-resolution marker strategy to traditional bi-parental mapping. In addition, the identified SNPs, candidate genes, and diverse variations in the root traits of soybean landraces will serve as a valuable basis for further application in genetic studies and the breeding of climate-resilient soybeans characterized by improved root traits.

Genetic dissection reveals the complex architecture of amino acid composition in soybean seeds

Five loci related to soybean protein and amino acid contents were colocated by performing linkage mapping and GWAS. The haplotype analysis showed that Glyma.08G109100 may be useful to improve the soybean seed composition. Soybean (Glycine max (L.) Merr.) seeds are good protein sources. Although genetic variation is abundant, natural variation in seed amino acids and their derived traits is lacking across soybean accessions. Here, we determined the contents of protein and 17 amino acids, obtained 36 derived traits based on the protein and total amino acid contents, and derived 34 traits based on seven amino acid family groups. Furthermore, we performed a linkage analysis of the contents of 17 amino acids and 73 amino acid-derived traits based on the recombinant inbred line (RIL)-derived Kefeng No. 1 × Nannong 1138-2. Six hundred thirty-nine quantitative trait loci (QTLs) were identified, explaining 6.07-39.00% of the phenotypic variation. Among these loci, five were detected in diverse soybean accessions using a genome-wide association study. A network analysis revealed that some loci that were significantly associated with multiple amino acids were tightly linked on chromosome 8 based on linkage disequilibrium values, which also further confirmed the results of the correlation analysis among amino acid traits. Through a combination of a genome-wide association study, linkage analysis, qRT-PCR, and genomic polymorphism comparison, Glyma.08G109100 on chromosome 8, which may affect amino acid contents, was selected. The haplotype analysis showed that Hap-T of Glyma.08G109100 may be useful to improve the contents of protein and 16 amino acids in soybean. This study provides new insights into the genetic basis of the amino acid composition in soybean seeds and may facilitate marker-based breeding of soybean with improved nutritional value.

Wildsoydb DataHub: a platform for accessing soybean multiomic datasets across multiple reference genomes

The Wildsoydb DataHub is an integrated interface for biologists and breeders to access soybean genomic resources easily, allowing them to fully utilize the results of genomic research.

From single- to multi-omics: future research trends in medicinal plants

Medicinal plants are the main source of natural metabolites with specialised pharmacological activities and have been widely examined by plant researchers. Numerous omics studies of medicinal plants have been performed to identify molecular markers of species and functional genes controlling key biological traits, as well as to understand biosynthetic pathways of bioactive metabolites and the regulatory mechanisms of environmental responses. Omics technologies have been widely applied to medicinal plants, including as taxonomics, transcriptomics, metabolomics, proteomics, genomics, pangenomics, epigenomics and mutagenomics. However, because of the complex biological regulation network, single omics usually fail to explain the specific biological phenomena. In recent years, reports of integrated multi-omics studies of medicinal plants have increased. Until now, there have few assessments of recent developments and upcoming trends in omics studies of medicinal plants. We highlight recent developments in omics research of medicinal plants, summarise the typical bioinformatics resources available for analysing omics datasets, and discuss related future directions and challenges. This information facilitates further studies of medicinal plants, refinement of current approaches and leads to new ideas.

Genome-wide association studies for soybean epicotyl length in two environments using 3VmrMLM

Germination of soybean seed is the imminent vital process after sowing. The status of plumular axis and radicle determine whether soybean seed can emerge normally. Epicotyl, an organ between cotyledons and first functional leaves, is essential for soybean seed germination, seedling growth and early morphogenesis. Epicotyl length (EL) is a quantitative trait controlled by multiple genes/QTLs. Here, the present study analyzes the phenotypic diversity and genetic basis of EL using 951 soybean improved cultivars and landraces from Asia, America, Europe and Africa. 3VmrMLM was used to analyze the associations between EL in 2016 and 2020 and 1,639,846 SNPs for the identification of QTNs and QTN-by-environment interactions (QEIs)".A total of 180 QTNs and QEIs associated with EL were detected. Among them, 74 QTNs (ELS_Q) and 16 QEIs (ELS_QE) were identified to be associated with ELS (epicotyl length of single plant emergence), and 60 QTNs (ELT_Q) and 30 QEIs (ELT_QE) were identified to be associated with ELT (epicotyl length of three seedlings). Based on transcript abundance analysis, GO (Gene Ontology) enrichment and haplotype analysis, ten candidate genes were predicted within nine genic SNPs located in introns, upstream or downstream, which were supposed to be directly or indirectly involved in the process of seed germination and seedling development., Of 10 candidate genes, two of them (Glyma.04G122400 and Glyma.18G183600) could possibly affect epicotyl length elongation. These results indicate the genetic basis of EL and provides a valuable basis for specific functional studies of epicotyl traits.

Genome-wide association studies reveal novel QTLs, QTL-by-environment interactions and their candidate genes for tocopherol content in soybean seed

Genome-wide association studies (GWAS) is an efficient method to detect quantitative trait locus (QTL), and has dissected many complex traits in soybean [Glycine max (L.) Merr.]. Although these results have undoubtedly played a far-reaching role in the study of soybean biology, environmental interactions for complex traits in traditional GWAS models are frequently overlooked. Recently, a new GWAS model, 3VmrMLM, was established to identify QTLs and QTL-by-environment interactions (QEIs) for complex traits. In this study, the GLM, MLM, CMLM, FarmCPU, BLINK, and 3VmrMLM models were used to identify QTLs and QEIs for tocopherol (Toc) content in soybean seed, including δ-Tocotrienol (δ-Toc) content, γ-Tocotrienol (γ-Toc) content, α-Tocopherol (α-Toc) content, and total Tocopherol (T-Toc) content. As a result, 101 QTLs were detected by the above methods in single-environment analysis, and 57 QTLs and 13 QEIs were detected by 3VmrMLM in multi-environment analysis. Among these QTLs, some QTLs (Group I) were repeatedly detected three times or by at least two models, and some QTLs (Group II) were repeatedly detected only by 3VmrMLM. In the two Groups, 3VmrMLM was able to correctly detect all known QTLs in group I, while good results were achieved in Group II, for example, 8 novel QTLs were detected in Group II. In addition, comparative genomic analysis revealed that the proportion of Glyma_max specific genes near QEIs was higher, in other words, these QEIs nearby genes are more susceptible to environmental influences. Finally, around the 8 novel QTLs, 11 important candidate genes were identified using haplotype, and validated by RNA-Seq data and qRT-PCR analysis. In summary, we used phenotypic data of Toc content in soybean, and tested the accuracy and reliability of 3VmrMLM, and then revealed novel QTLs, QEIs and candidate genes for these traits. Hence, the 3VmrMLM model has broad prospects and potential for analyzing the genetic structure of complex quantitative traits in soybean.

PlantExp: a platform for exploration of gene expression and alternative splicing based on public plant RNA-seq samples

Over the last decade, RNA-seq has produced a massive amount of plant transcriptomic sequencing data deposited in public databases. Reanalysis of these public datasets can generate additional novel hypotheses not included in original studies. However, the large data volume and the requirement for specialized computational resources and expertise present a barrier for experimental biologists to explore public repositories. Here, we introduce PlantExp (https://biotec.njau.edu.cn/plantExp), a database platform for exploration of plant gene expression and alternative splicing profiles based on 131 423 uniformly processed publicly available RNA-seq samples from 85 species in 24 plant orders. In addition to two common retrieval accesses to gene expression and alternative splicing profiles by functional terms and sequence similarity, PlantExp is equipped with four online analysis tools, including differential expression analysis, specific expression analysis, co-expression network analysis and cross-species expression conservation analysis. With these online analysis tools, users can flexibly customize sample groups to reanalyze public RNA-seq datasets and obtain new insights. Furthermore, it offers a wide range of visualization tools to help users intuitively understand analysis results. In conclusion, PlantExp provides a valuable data resource and analysis platform for plant biologists to utilize public RNA-seq. datasets.

Identifying Soybean Pod Borer (Leguminivora glycinivorella) Resistance QTLs and the Mechanism of Induced Defense Using Linkage Mapping and RNA-Seq Analysis

The soybean pod borer (Leguminivora glycinivorella) (SPB) is a major cause of soybean (Glycine max L.) yield losses in northeast Asia, thus it is desirable to elucidate the resistance mechanisms involved in soybean response to the SPB. However, few studies have mapped SPB-resistant quantitative trait loci (QTLs) and deciphered the response mechanism in soybean. Here, we selected two soybean varieties, JY93 (SPB-resistant) and K6 (SPB-sensitive), to construct F₂ and F_2:3 populations for QTL mapping and collected pod shells before and after SPB larvae chewed on the two parents to perform RNA-Seq, which can identify stable QTLs and explore the response mechanism of soybean to the SPB. The results show that four QTLs underlying SPB damage to seeds were detected on chromosomes 4, 9, 13, and 15. Among them, qESP-9-1 was scanned in all environments, hence it can be considered a stable QTL. All QTLs explained 0.79 to 6.09% of the phenotypic variation. Meanwhile, 2298 and 3509 DEGs were identified for JY93 and K6, respectively, after the SPB attack, and most of these genes were upregulated. Gene Ontology enrichment results indicated that the SPB-induced and differently expressed genes in both parents are involved in biological processes such as wound response, signal transduction, immune response, and phytohormone pathways. Interestingly, secondary metabolic processes such as flavonoid synthesis were only significantly enriched in the upregulated genes of JY93 after SPB chewing compared with K6. Finally, we identified 18 candidate genes related to soybean pod borer resistance through the integration of QTL mapping and RNA-Seq analysis. Seven of these genes had similar expression patterns to the mapping parents in four additional soybean germplasm after feeding by the SPB. These results provide additional knowledge of the early response and induced defense mechanisms against the SPB in soybean, which could help in breeding SPB-resistant soybean accessions.

Lateral transfers lead to the birth of momilactone biosynthetic gene clusters in grass

Momilactone A, an important plant labdane-related diterpenoid, functions as a phytoalexin against pathogens and an allelochemical against neighboring plants. The genes involved in the biosynthesis of momilactone A are found in clusters, i.e., momilactone A biosynthetic gene clusters (MABGCs), in the rice and barnyardgrass genomes. In addition, we know little about the origin and evolution of MABGCs. Here, we integrated results from comprehensive phylogeny and comparative genomic analyses of the core genes of MABGC-like clusters and MABGCs in 40 monocot plant genomes, providing convincing evidence for the birth and evolution of MABGCs in grass species. The MABGCs found in the PACMAD clade of the core grass lineage (including Panicoideae and Chloridoideae) originated from a MABGC-like cluster in Triticeae (BOP clade) via lateral gene transfer (LGT) and followed by recruitment of MAS1/2 and CYP76L1 genes. The MABGCs in Oryzoideae originated from PACMAD through another LGT event and lost CYP76L1 afterwards. The Oryza MABGC and another Oryza diterpenoid cluster c2BGC are two distinct clusters, with the latter originating from gene duplication and relocation within Oryzoideae. Further comparison of the expression patterns of the MABGC genes between rice and barnyardgrass in response to pathogen infection and allelopathy provides novel insights into the functional innovation of MABGCs in plants. Our results demonstrate LGT-mediated origination of MABGCs in grass and shed lights into the evolutionary innovation and optimization of plant biosynthetic pathways.

Characteristics of members of IGT family genes in controlling rice root system architecture and tiller development

Root system architecture (RSA) and tiller are important agronomic traits. However, the mechanisms of the IGT family genes regulate RSA and tiller development in different rice varieties remain unclear. In this study, we demonstrated that 38 rice varieties obtained from Yuanyang Hani's terraced fields with different RSA and could be classified into six groups based on the ratio of root length and width. We found a positive correlation between RSA (including root width, length, and area) and tiller number in most of rice varieties. Furthermore, the IGT family genes Deeper Rooting 1 (DRO1), LAZY1, TAC1, and qSOR1 showed different expression patterns when rice grown under irrigation and drought conditions. Moreover, the qSOR1 gene had higher levels in the roots and tillers, and accompanied with higher levels of PIN1b gene in roots when rice grown under drought environmental condition. DRO1 gene had two single nucleotide polymorphisms (SNPs) in the exon 3 sequences and showed different expression patterns in the roots and tillers of the 38 rice varieties. Overexpression of DRO1 with a deletion of exon 5 caused shorter root length, less lateral roots and lower levels of LAZY1, TAC1, and qSOR1. Further protein interaction network, microRNA targeting and co-expression analysis showed that DRO1 plays a critical role in the root and tiller development associated with auxin transport. These data suggest that the RSA and tiller development are regulated by the IGT family genes in an intricate network way, which is tightly related to rice genetic background in rice adapting to different environmental conditions.

GhWRKY33 Interacts with GhTIFY10A to Synergistically Modulate Both Ageing and JA-Mediated Leaf Senescence in Arabidopsis

WRKY transcription factors play critical roles in the modulation of transcriptional changes during leaf senescence, but the underlying mechanisms controlled by them in this progress still remain enigmatic. In this study, Gossypium hirsutum WRKY DNA-binding protein 33 (GhWRKY33) was characterized as a negative regulator of both ageing and JA-mediated leaf senescence. The overexpression of GhWRKY33 in Arabidopsis greatly delayed leaf senescence, as determined by elevated chlorophyll content, lower H₂O₂ content, and reduced expression of several senescence-associated genes (SAGs). An electrophoretic mobility shift assay (EMSA) and transient dual–luciferase reporter assay revealed that GhWRKY33 could bind to the promoters of both AtSAG12 and Ghcysp and suppress their expression. Yeast two-hybrid (Y2H) and firefly luciferase complementation imaging (LUC) assays showed that GhWRKY33 could interact with GhTIFY10A. Similarly, the overexpression of GhTIFY10A in Arabidopsis also dramatically delayed leaf senescence. Furthermore, both GhWRKY33 and GhTIFY10A negatively regulate JA-mediated leaf senescence. In addition, a transientdual-luciferase reporter assay indicated that GhWRKY33 and GhTIFY10A could function synergistically to inhibit the expression of both AtSAG12 and Ghcysp. Thus, our work suggested that GhWRKY33 may function as a negative regulator to modulate both ageing and JA-mediated leaf senescence and also contributes to a basis for further functional studies on cotton leaf senescence.

Genome-Wide Association Study Reveals a Genetic Mechanism of Salt Tolerance Germinability in Rice (Oryza sativa L.)

Salt stress is one of the factors that limits rice production, and an important task for researchers is to cultivate rice with strong salt tolerance. In this study, 211 rice accessions were used to determine salt tolerance germinability (STG) indices and conduct a genome-wide association study (GWAS) using 36,727 SNPs. The relative germination energy (RGE), relative germination index (RGI), relative vigor index (RVI), relative mean germination time (RMGT), relative shoot length (RSL), and relative root length (RRL) were used to determine the STG indices in rice. A total of 43 QTLs, including 15 for the RGE, 6 for the RGI, 7 for the RVI, 3 for the RMGT, 1 for the RSL, and 11 for the RRL, were identified on nine chromosome regions under 60 and 100 mM NaCl conditions. For these STG-related QTLs, 18 QTLs were co-localized with previous studies, and some characterized salt-tolerance genes, such as OsCOIN, OsHsp17.0, and OsDREB2A, are located in these QTL candidates. Among the 25 novel QTLs, qRGE60-1-2 co-localized with qRGI60-1-1 on chromosome 1, and qRGE60-3-1 and qRVI60-3-1 co-localized on chromosome 3. According to the RNA-seq database, 16 genes, including nine for qRGE60-1-2 (qRGI60-1-1) and seven for qRGE60-3-1 (qRVI60-3-1), were found to show significant differences in their expression levels between the control and salt treatments. Furthermore, the expression patterns of these differentially expressed genes were analyzed, and nine genes (five for qRGE60-1-2 and four for qRGE60-3-1) were highly expressed in embryos at the germination stage. Haplotype analysis of these nine genes showed that the rice varieties with elite haplotypes in the LOC_Os03g13560, LOC_Os03g13840, and LOC_Os03g14180 genes had high STG. GWAS validated the known genes underlying salt tolerance and identified novel loci that could enrich the current gene pool related to salt tolerance. The resources with high STG and significant loci identified in this study are potentially useful in breeding for salt tolerance.

An ORFeome of rice E3 ubiquitin ligases for global analysis of the ubiquitination interactome

BACKGROUND

Ubiquitination is essential for many cellular processes in eukaryotes, including 26S proteasome-dependent protein degradation, cell cycle progression, transcriptional regulation, and signal transduction. Although numerous ubiquitinated proteins have been empirically identified, their cognate ubiquitin E3 ligases remain largely unknown.

RESULTS

Here, we generate a complete ubiquitin E3 ligase-encoding open reading frames (UbE3-ORFeome) library containing 98.94% of the 1515 E3 ligase genes in the rice (Oryza sativa L.) genome. In the test screens with four known ubiquitinated proteins, we identify both known and new E3s. The interaction and degradation between several E3s and their substrates are confirmed in vitro and in vivo. In addition, we identify the F-box E3 ligase OsFBK16 as a hub-interacting protein of the phenylalanine ammonia lyase family OsPAL1-OsPAL7. We demonstrate that OsFBK16 promotes the degradation of OsPAL1, OsPAL5, and OsPAL6. Remarkably, we find that overexpression of OsPAL1 or OsPAL6 as well as loss-of-function of OsFBK16 in rice displayed enhanced blast resistance, indicating that OsFBK16 degrades OsPALs to negatively regulate rice immunity.

CONCLUSIONS

The rice UbE3-ORFeome is the first complete E3 ligase library in plants and represents a powerful proteomic resource for rapid identification of the cognate E3 ligases of ubiquitinated proteins and establishment of functional E3-substrate interactome in plants.

Genome-Wide Analysis of Type-III Polyketide Synthases in Wheat and Possible Roles in Wheat Sheath-Blight Resistance

The enzymes in the chalcone synthase family, also known as type-III polyketide synthases (PKSs), play important roles in the biosynthesis of various plant secondary metabolites and plant adaptation to environmental stresses. There have been few detailed reports regarding the gene and tissue expression profiles of the PKS (TaPKS) family members in wheat (Triticum aestivum L.). In this study, 81 candidate TaPKS genes were identified in the wheat genome, which were designated as TaPKS1–81. Phylogenetic analysis divided the TaPKS genes into two groups. TaPKS gene family expansion mainly occurred via tandem duplication and fragment duplication. In addition, we analyzed the physical and chemical properties, gene structures, and cis-acting elements of TaPKS gene family members. RNA-seq analysis showed that the expression of TaPKS genes was tissue-specific, and their expression levels differed before and after infection with Rhizoctonia cerealis. The expression levels of four TaPKS genes were also analyzed via qRT-PCR after treatment with methyl jasmonate, salicylic acid, abscisic acid, and ethylene. In the present study, we systematically identified and analyzed TaPKS gene family members in wheat, and our findings may facilitate the cloning of candidate genes associated with resistance to sheath blight in wheat.

Characterization of two constitutive promoters RPS28 and EIF1 for studying soybean growth, development, and symbiotic nodule development

Native promoters that can drive high and stable transgene expression are important tools for modifying plant traits. Although several such promoters have been reported in soybean (Glycine max), few of them function at multiple growth and development stages and during nodule development. Here, we report that the promoters of 40S RIBOSOMAL PROTEIN SMALL SUBUNIT S28 (RPS28) and EUKARYOTIC TRANSLATION INITIATION FACTOR 1 (EIF1) are ideal for high expression of transgene. Through bioinformatic analysis, we determined that RPS28 and EIF1 were highly expressed during soybean growth and development, nodule development, and various biotic and abiotic stresses. Fusion of both RPS28 and EIF1 promoters, with or without their first intron, with the reporter gene β-GLUCURONIDASE (uidA) in transgenic soybean, resulted in high GUS activity in seedlings, seeds, and nodules. Fluorimetric GUS assays showed that the RPS28 promoter and the EIF1 promoter yielded high expression, comparable to the soybean Ubiquitin (GmUbi) promoter. RPS28 and EIF1 promoters were also highly expressed in Arabidopsis thaliana and Nicotiana benthamiana. Our results indicate the potential of RPS28 and EIF1 promoters to facilitate future genetic engineering and breeding to improve the quality and yield of soybean, as well as in a wide variety of other plant species.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-022-00073-6.

Deciphering Haplotypic Variation and Gene Expression Dynamics Associated with Nutritional and Cooking Quality in Rice

Nutritional quality improvement of rice is the key to ensure global food security. Consequently, enormous efforts have been made to develop genomics and transcriptomics resources for rice. The available omics resources along with the molecular understanding of trait development can be utilized for efficient exploration of genetic resources for breeding programs. In the present study, 80 genes known to regulate the nutritional and cooking quality of rice were extensively studied to understand the haplotypic variability and gene expression dynamics. The haplotypic variability of selected genes were defined using whole-genome re-sequencing data of ~4700 diverse genotypes. The analytical workflow identified 133 deleterious single-nucleotide polymorphisms, which are predicted to affect the gene function. Furthermore, 788 haplotype groups were defined for 80 genes, and the distribution and evolution of these haplotype groups in rice were described. The nucleotide diversity for the selected genes was significantly reduced in cultivated rice as compared with that in wild rice. The utility of the approach was successfully demonstrated by revealing the haplotypic association of chalk5 gene with the varying degree of grain chalkiness. The gene expression atlas was developed for these genes by analyzing RNA-Seq transcriptome profiling data from 102 independent sequence libraries. Subsequently, weighted gene co-expression meta-analyses of 11,726 publicly available RNAseq libraries identified 19 genes as the hub of interactions. The comprehensive analyses of genetic polymorphisms, allelic distribution, and gene expression profiling of key quality traits will help in exploring the most desired haplotype for grain quality improvement. Similarly, the information provided here will be helpful to understand the molecular mechanism involved in the development of nutritional and cooking quality traits in rice.