Beyond gene ontology (GO): using biocuration approach to improve the gene nomenclature and functional annotation of rice S-domain kinase subfamily

Plant Reactome Knowledgebase: empowering plant pathway exploration and OMICS data analysis

Plant Reactome (https://plantreactome.gramene.org) is a freely accessible, comprehensive plant pathway knowledgebase. It provides curated reference pathways from rice (Oryza sativa) and gene-orthology-based pathway projections to 129 additional species, spanning single-cell photoautotrophs, non-vascular plants, and higher plants, thus encompassing a wide-ranging taxonomic diversity. Currently, Plant Reactome houses a collection of 339 reference pathways, covering metabolic and transport pathways, hormone signaling, genetic regulations of developmental processes, and intricate transcriptional networks that orchestrate a plant's response to abiotic and biotic stimuli. Beyond being a mere repository, Plant Reactome serves as a dynamic data discovery platform. Users can analyze and visualize omics data, such as gene expression, gene-gene interaction, proteome, and metabolome data, all within the rich context of plant pathways. Plant Reactome is dedicated to fostering data interoperability, upholding global data standards, and embracing the tenets of the Findable, Accessible, Interoperable and Re-usable (FAIR) data policy.

Reference genome of the nutrition-rich orphan crop chia (Salvia hispanica) and its implications for future breeding

Chia (Salvia hispanica L.) is one of the most popular nutrition-rich foods and pseudocereal crops of the family Lamiaceae. Chia seeds are a rich source of proteins, polyunsaturated fatty acids (PUFAs), dietary fibers, and antioxidants. In this study, we present the assembly of the chia reference genome, which spans 303.6 Mb and encodes 48,090 annotated protein-coding genes. Our analysis revealed that ~42% of the chia genome harbors repetitive content, and identified ~3 million single nucleotide polymorphisms (SNPs) and 15,380 simple sequence repeat (SSR) marker sites. By investigating the chia transcriptome, we discovered that ~44% of the genes undergo alternative splicing with a higher frequency of intron retention events. Additionally, we identified chia genes associated with important nutrient content and quality traits, such as the biosynthesis of PUFAs and seed mucilage fiber (dietary fiber) polysaccharides. Notably, this is the first report of in-silico annotation of a plant genome for protein-derived small bioactive peptides (biopeptides) associated with improving human health. To facilitate further research and translational applications of this valuable orphan crop, we have developed the Salvia genomics database (SalviaGDB), accessible at https://salviagdb.org.

Genetic dissection of monosaccharides contents in rice whole grain using genome-wide association study

The simplest form of carbohydrates are monosaccharides which are the building blocks for the synthesis of polymers or complex carbohydrates. Monosaccharide contents of 197 rice accessions were quantified by HPAEC-PAD in rice (Oryza sativa L.) whole grain (RWG). A genome-wide association study (GWAS) was carried out using 33,812 single nucleotide polymorphisms (SNPs) to identify corresponding genomic regions influencing neutral monosaccharides contents. In total, 49 GWAS signals contained in 17 genomic regions (quantitative trait loci [QTLs]) on seven chromosomes of rice were determined to be associated with monosaccharides contents of whole grain. The QTLs were found for fucose (1), mannose (1), xylose (2), arabinose (2), galactose (4), and rhamnose (7) contents, all of which are novel. Based on co-location of annotated rice genes in the vicinity of GWAS signals, the constituents of the whole grain were associated with the following candidate genes: arabinose content with α-N-arabinofuranosidase, pectinesterase inhibitor, and glucosamine-fructose-6-phosphate aminotransferase 1; xylose content with ZOS1-10 (a C2H2 zinc finger transcription factor [TF]); mannose content with aldose 1-epimerase-like protein and a MYB family TF; galactose content with a GT8 family member (galacturonosyltransferase-like 3), a GRAS family TF, and a GH16 family member (xyloglucan endotransglucosylase/hydrolase xyloglucan 23); fucose content with gibberellin 20 oxidase and a lysine-rich arabinogalactan protein 19, and finally rhamnose content with myo-inositol-1-phosphate synthase, UDP-arabinopyranose mutase, and COBRA-like protein precursor. The results of this study should improve our understanding of the genetic basis of the factors that might be involved in the biosynthesis, regulation, and turnover of monosaccharides in RWG, aiming to enhance the nutritional value of rice grain and impact the related industries.

Comparison of Transcriptome between Tolerant and Susceptible Rice Cultivar Reveals Positive and Negative Regulators of Response to Rhizoctonia solani in Rice

Rice (Oryza sativa L.) is one of the world's most crucial food crops, as it currently supports more than half of the world's population. However, the presence of sheath blight (SB) caused by Rhizoctonia solani has become a significant issue for rice agriculture. This disease is responsible for causing severe yield losses each year and is a threat to global food security. The breeding of SB-resistant rice varieties requires a thorough understanding of the molecular mechanisms involved and the exploration of immune genes in rice. To this end, we conducted a screening of rice cultivars for resistance to SB and compared the transcriptome based on RNA-seq between the most tolerant and susceptible cultivars. Our study revealed significant transcriptomic differences between the tolerant cultivar ZhengDao 22 (ZD) and the most susceptible cultivar XinZhi No.1 (XZ) in response to R. solani invasion. Specifically, the tolerant cultivar showed 7066 differentially expressed genes (DEGs), while the susceptible cultivar showed only 60 DEGs. In further analysis, we observed clear differences in gene category between up- and down-regulated expression of genes (uDEGs and dDEGs) based on Gene Ontology (GO) classes in response to infection in the tolerant cultivar ZD, and then identified uDEGs related to cell surface pattern recognition receptors, the Ca2+ ion signaling pathway, and the Mitogen-Activated Protein Kinase (MAPK) cascade that play a positive role against R. solani. In addition, DEGs of the jasmonic acid and ethylene signaling pathways were mainly positively regulated, whereas DEGs of the auxin signaling pathway were mainly negatively regulated. Transcription factors were involved in the immune response as either positive or negative regulators of the response to this pathogen. Furthermore, our results showed that chloroplasts play a crucial role and that reduced photosynthetic capacity is a critical feature of this response. The results of this research have important implications for better characterization of the molecular mechanism of SB resistance and for the development of resistant cultivars through molecular breeding methods.

Exploring Pan-Genomes: An Overview of Resources and Tools for Unraveling Structure, Function, and Evolution of Crop Genes and Genomes

The availability of multiple sequenced genomes from a single species made it possible to explore intra- and inter-specific genomic comparisons at higher resolution and build clade-specific pan-genomes of several crops. The pan-genomes of crops constructed from various cultivars, accessions, landraces, and wild ancestral species represent a compendium of genes and structural variations and allow researchers to search for the novel genes and alleles that were inadvertently lost in domesticated crops during the historical process of crop domestication or in the process of extensive plant breeding. Fortunately, many valuable genes and alleles associated with desirable traits like disease resistance, abiotic stress tolerance, plant architecture, and nutrition qualities exist in landraces, ancestral species, and crop wild relatives. The novel genes from the wild ancestors and landraces can be introduced back to high-yielding varieties of modern crops by implementing classical plant breeding, genomic selection, and transgenic/gene editing approaches. Thus, pan-genomic represents a great leap in plant research and offers new avenues for targeted breeding to mitigate the impact of global climate change. Here, we summarize the tools used for pan-genome assembly and annotations, web-portals hosting plant pan-genomes, etc. Furthermore, we highlight a few discoveries made in crops using the pan-genomic approach and future potential of this emerging field of study.

Biocuration of a Transcription Factors Network Involved in Submergence Tolerance during Seed Germination and Coleoptile Elongation in Rice (Oryza sativa)

Modeling biological processes and genetic-regulatory networks using in silico approaches provides a valuable framework for understanding how genes and associated allelic and genotypic differences result in specific traits. Submergence tolerance is a significant agronomic trait in rice; however, the gene-gene interactions linked with this polygenic trait remain largely unknown. In this study, we constructed a network of 57 transcription factors involved in seed germination and coleoptile elongation under submergence. The gene-gene interactions were based on the co-expression profiles of genes and the presence of transcription factor binding sites in the promoter region of target genes. We also incorporated published experimental evidence, wherever available, to support gene-gene, gene-protein, and protein-protein interactions. The co-expression data were obtained by re-analyzing publicly available transcriptome data from rice. Notably, this network includes OSH1, OSH15, OSH71, Sub1B, ERFs, WRKYs, NACs, ZFP36, TCPs, etc., which play key regulatory roles in seed germination, coleoptile elongation and submergence response, and mediate gravitropic signaling by regulating OsLAZY1 and/or IL2. The network of transcription factors was manually biocurated and submitted to the Plant Reactome Knowledgebase to make it publicly accessible. We expect this work will facilitate the re-analysis/re-use of OMICs data and aid genomics research to accelerate crop improvement.

Comparative transcriptomics and co-expression networks reveal cultivar-specific molecular signatures associated with reproductive-stage cold stress in rice

The resistance of Huaidao5 results from the high constitutive expression of tolerance genes, while that of Huaidao9 is due to the cold-induced resistance in flag leaves and panicles. The regulation mechanism of rice seedlings' cold tolerance is relatively clear, and knowledge of its underlying mechanisms at the reproductive stage is limited. We performed differential expression and co-expression network analyses to transcriptomes from panicle and flag leaf tissues of a cold-tolerant cultivar (Huaidao5), and a sensitive cultivar (Huaidao9), under reproductive-stage cold stress. The results revealed that the expression levels of genes in stress-related pathways such as MAPK signaling pathway, diterpenoid biosynthesis, glutathione metabolism, plant-pathogen interaction and plant hormone signal transduction were constitutively highly expressed in Huaidao5, especially in panicles. Moreover, the Hudaidao5's panicle sample-specific (under cold) module contained some genes related to rice yield, such as GW5L, GGC2, SG1 and CTPS1. However, the resistance of Huaidao9 was derived from the induced resistance to cold in flag leaves and panicles. In the flag leaves, the responses included a series of stress response and signal transduction, while in the panicles nitrogen metabolism was severely affected, especially 66 endosperm-specific genes. Through integrating differential expression with co-expression networks, we predicted 161 candidate genes (79 cold-responsive genes common to both cultivars and 82 cold-tolerance genes associated with differences in cold tolerance between cultivars) potentially affecting cold response/tolerance, among which 85 (52.80%) were known to be cold-related genes. Moreover, 52 (65.82%) cold-responsive genes (e.g., TIFY11C, LSK1 and LPA) could be confirmed by previous transcriptome studies and 72 (87.80%) cold-tolerance genes (e.g., APX5, OsFbox17 and OsSTA109) were located within QTLs associated with cold tolerance. This study provides an efficient strategy for further discovery of mechanisms of cold tolerance in rice.

Meta-analysis of QTLs and candidate genes associated with seed germination in rice (Oryza sativa L.)

Seed germination is one of the critical stages of plant life, and many quantitative trait loci (QTLs) control this complex trait. Meta-analysis of QTLs is a powerful computational technique for estimating the most stable QTLs regardless of the population's genetic background. Besides, this analysis effectively narrows down the confidence interval (CI) to identify candidate genes (CGs) and marker development. In the current study, a comprehensive genome-wide meta-analysis was performed on QTLs associated with germination in rice. This analysis was conducted based on the data reported over the last two decades. In this case, various analyses were performed, including seed germination rate, plumule length, radicle length, germination percentage, coleoptile length, coleorhiza length, radicle fresh weight, germination potential, and germination index. A total of 67 QTLs were projected onto a reference map for these traits and then integrated into 32 meta-QTLs (MQTLs) to provide a genetic framework for seed germination. The average CI of MQTLs was considerably reduced from 15.125 to 8.73 cM compared to the initial QTLs. This situation identified 728 well-known functionally characterized genes and novel putative CGs for investigated traits. The fold change calculation demonstrated that 155 CGs had significant changes in expression analysis. In this case, 112 and 43 CGs were up-regulated and down-regulated during germination, respectively. This study provides an overview and compares genetic loci controlling traits related to seed germination in rice. The findings can bridge the gap between QTLs and CGs for seed germination.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01232-1.

First Plant Cell Atlas symposium report

The Plant Cell Atlas (PCA) community hosted a virtual symposium on December 9 and 10, 2021 on single cell and spatial omics technologies. The conference gathered almost 500 academic, industry, and government leaders to identify the needs and directions of the PCA community and to explore how establishing a data synthesis center would address these needs and accelerate progress. This report details the presentations and discussions focused on the possibility of a data synthesis center for a PCA and the expected impacts of such a center on advancing science and technology globally. Community discussions focused on topics such as data analysis tools and annotation standards; computational expertise and cyber-infrastructure; modes of community organization and engagement; methods for ensuring a broad reach in the PCA community; recruitment, training, and nurturing of new talent; and the overall impact of the PCA initiative. These targeted discussions facilitated dialogue among the participants to gauge whether PCA might be a vehicle for formulating a data synthesis center. The conversations also explored how online tools can be leveraged to help broaden the reach of the PCA (i.e., online contests, virtual networking, and social media stakeholder engagement) and decrease costs of conducting research (e.g., virtual REU opportunities). Major recommendations for the future of the PCA included establishing standards, creating dashboards for easy and intuitive access to data, and engaging with a broad community of stakeholders. The discussions also identified the following as being essential to the PCA's success: identifying homologous cell-type markers and their biocuration, publishing datasets and computational pipelines, utilizing online tools for communication (such as Slack), and user-friendly data visualization and data sharing. In conclusion, the development of a data synthesis center will help the PCA community achieve these goals by providing a centralized repository for existing and new data, a platform for sharing tools, and new analytical approaches through collaborative, multidisciplinary efforts. A data synthesis center will help the PCA reach milestones, such as community-supported data evaluation metrics, accelerating plant research necessary for human and environmental health.

The Rice Serine/Arginine Splicing Factor RS33 Regulates Pre-mRNA Splicing during Abiotic Stress Responses

Abiotic stresses profoundly affect plant growth and development and limit crop productivity. Pre-mRNA splicing is a major form of gene regulation that helps plants cope with various stresses. Serine/arginine (SR)-rich splicing factors play a key role in pre-mRNA splicing to regulate different biological processes under stress conditions. Alternative splicing (AS) of SR transcripts and other transcripts of stress-responsive genes generates multiple splice isoforms that contribute to protein diversity, modulate gene expression, and affect plant stress tolerance. Here, we investigated the function of the plant-specific SR protein RS33 in regulating pre-mRNA splicing and abiotic stress responses in rice. The loss-of-function mutant rs33 showed increased sensitivity to salt and low-temperature stresses. Genome-wide analyses of gene expression and splicing in wild-type and rs33 seedlings subjected to these stresses identified multiple splice isoforms of stress-responsive genes whose AS are regulated by RS33. The number of RS33-regulated genes was much higher under low-temperature stress than under salt stress. Our results suggest that the plant-specific splicing factor RS33 plays a crucial role during plant responses to abiotic stresses.

Plant lectins and their many roles: Carbohydrate-binding and beyond

Lectins are ubiquitous proteins that reversibly bind to specific carbohydrates and, thus, serve as readers of the sugar code. In photosynthetic organisms, lectin family proteins play important roles in capturing and releasing photosynthates via an endogenous lectin cycle. Often, lectin proteins consist of one or more lectin domains in combination with other types of domains. This structural diversity of lectins is the basis for their current classification, which is consistent with their diverse functions in cell signaling associated with growth and development, as well as in the plant's response to biotic, symbiotic, and abiotic stimuli. Furthermore, the lectin family shows evolutionary expansion that has distinct clade-specific signatures. Although the function(s) of many plant lectin family genes are unknown, studies in the model plant Arabidopsis thaliana have provided insights into their diverse roles. Here, we have used a biocuration approach rooted in the critical review of scientific literature and information available in the public genomic databases to summarize the expression, localization, and known functions of lectins in Arabidopsis. A better understanding of the structure and function of lectins is expected to aid in improving agricultural productivity through the manipulation of candidate genes for breeding climate-resilient crops, or by regulating metabolic pathways by applications of plant growth regulators.