New Data and New Features of the FunRiceGenes (Functionally Characterized Rice Genes) Database: 2021 Update

Potential Roles of the GRF Transcription Factors in Sorghum Internodes during Post-Reproductive Stages

Growth-regulating factor (GRF) is a plant-specific family of transcription factors crucial for meristem development and plant growth. Sorghum (Sorghum bicolor L. Moench) is a cereal species widely used for food, feed and fuel. While sorghum stems are important biomass components, the regulation of stem development and the carbohydrate composition of the stem tissues remain largely unknown. Here, we identified 11 SbGRF-encoding genes and found the SbGRF expansion driven by whole-genome duplication events. By comparative analyses of GRFs between rice and sorghum, we demonstrated the divergence of whole-genome duplication (WGD)-derived OsGRFs and SbGRFs. A comparison of SbGRFs' expression profiles supports that the WGD-duplicated OsGRFs and SbGRFs experienced distinct evolutionary trajectories, possibly leading to diverged functions. RNA-seq analysis of the internode tissues identified several SbGRFs involved in internode elongation, maturation and cell wall metabolism. We constructed co-expression networks with the RNA-seq data of sorghum internodes. Network analysis discovered that SbGRF1, 5 and 7 could be involved in the down-regulation of the biosynthesis of cell wall components, while SbGRF4, 6, 8 and 9 could be associated with the regulation of cell wall loosening, reassembly and/or starch biosynthesis. In summary, our genome-wide analysis of SbGRFs reveals the distinct evolutionary trajectories of WGD-derived SbGRF pairs. Importantly, expression analyses highlight previously unknown functions of several SbGRFs in internode elongation, maturation and the potential involvement in the metabolism of the cell wall and starch during post-anthesis stages.

Divergent Roles of the Auxin Response Factors in Lemongrass (Cymbopogon flexuosus (Nees ex Steud.) W. Watson) during Plant Growth

Auxin Response Factors (ARFs) make up a plant-specific transcription factor family that mainly couples perception of the phytohormone, auxin, and gene expression programs and plays an important and multi-faceted role during plant growth and development. Lemongrass (Cymbopogon flexuosus) is a representative Cymbopogon species widely used in gardening, beverages, fragrances, traditional medicine, and heavy metal phytoremediation. Biomass yield is an important trait for several agro-economic purposes of lemongrass, such as landscaping, essential oil production, and phytoremediation. Therefore, we performed gene mining of CfARFs and identified 26 and 27 CfARF-encoding genes in each of the haplotype genomes of lemongrass, respectively. Phylogenetic and domain architecture analyses showed that CfARFs can be divided into four groups, among which groups 1, 2, and 3 correspond to activator, repressor, and ETTN-like ARFs, respectively. To identify the CfARFs that may play major roles during the growth of lemongrass plants, RNA-seq was performed on three tissues (leaf, stem, and root) and four developmental stages (3-leaf, 4-leaf, 5-leaf. and mature stages). The expression profiling of CfARFs identified several highly expressed activator and repressor CfARFs and three CfARFs (CfARF3, 18, and 35) with gradually increased levels during leaf growth. Haplotype-resolved transcriptome analysis revealed that biallelic expression dominance is frequent among CfARFs and contributes to their gene expression patterns. In addition, co-expression network analysis identified the modules enriched with CfARFs. By establishing orthologous relationships among CfARFs, sorghum ARFs, and maize ARFs, we showed that CfARFs were mainly expanded by whole-genome duplications, and that the duplicated CfARFs might have been divergent due to differential expression and variations in domains and motifs. Our work provides a detailed catalog of CfARFs in lemongrass, representing a first step toward characterizing CfARF functions, and may be useful in molecular breeding to enhance lemongrass plant growth.

Functional differentiation and genetic diversity of rice cation exchanger (CAX) genes and their potential use in rice improvement

Cation exchanger (CAX) genes play an important role in plant growth/development and response to biotic and abiotic stresses. Here, we tried to obtain important information on the functionalities and phenotypic effects of CAX gene family by systematic analyses of their expression patterns, genetic diversity (gene CDS haplotypes, structural variations, gene presence/absence variations) in 3010 rice genomes and nine parents of 496 Huanghuazhan introgression lines, the frequency shifts of the predominant gcHaps at these loci to artificial selection during modern breeding, and their association with tolerances to several abiotic stresses. Significant amounts of variation also exist in the cis-regulatory elements (CREs) of the OsCAX gene promoters in 50 high-quality rice genomes. The functional differentiation of OsCAX gene family were reflected primarily by their tissue and development specific expression patterns and in varied responses to different treatments, by unique sets of CREs in their promoters and their associations with specific agronomic traits/abiotic stress tolerances. Our results indicated that OsCAX1a and OsCAX2 as general signal transporters were in many processes of rice growth/development and responses to diverse environments, but they might be of less value in rice improvement. OsCAX1b, OsCAX1c, OsCAX3 and OsCAX4 was expected to be of potential value in rice improvement because of their associations with specific traits, responsiveness to specific abiotic stresses or phytohormones, and relatively high gcHap and CRE diversity. Our strategy was demonstrated to be highly efficient to obtain important genetic information on genes/alleles of specific gene family and can be used to systematically characterize the other rice gene families.

Development of activation-tagged gain-of-functional mutants in indica rice line (BPT 5204) for sheath blight resistance

BACKGROUND

The development of sheath blight (ShB) resistance varieties has been a challenge for scientists for long time in rice. Activation tagging is an efficient gain-of-function mutation approach to create novel phenotypes and to identify their underlying genes. In this study, a mutant population was developed employing activation tagging in the recalcitrant indica rice (Oryza sativa L.) cv. BPT 5204 (Samba Mahsuri) through activation tagging.

METHODS AND RESULTS

In this study, we have generated more than 1000 activation tagged lines in indica rice, from these mutant population 38 (GFP- RFP+) stable Ds plants were generated through germinal transposition at T2 generation based on molecular analysis and seeds selected on hygromycin (50 mg/L) containing medium segregation analyses confirmed that the transgene inherited as mendelian segregation ratio of 3:1 (3 resistant: 1 susceptible). Of them, five stable activation tagged Ds lines (M-Ds-1, M-Ds-2, M-Ds-3, M-Ds-4 and M-Ds-5) were selected based on phenotypic observation through screening for sheath blight (ShB) resistance caused by fungal pathogen Rhizoctonia solani (R. solani),. Among them, M-Ds-3 and M-Ds-5 lines showed significant resistance for ShB over other tagged lines and wild type (WT) plants. Furthermore, analysed for launch pad insertion through TAIL-PCR results and mapped on corresponding rice chromosomes. Flanking sequence and gene expression analysis revealed that the upregulation of glycoside hydrolase-OsGH or similar to Class III chitinase homologue (LOC_Os08g40680) in M-Ds-3 and a hypothetical protein gene (LOC_Os01g55000) in M-Ds-5 are potential candidate genes for sheath blight resistance in rice.

CONCLUSION

In the present study, we developed Ac-Ds based ShB resistance gain-of-functional mutants through activation tagging in rice. These activation tagged mutant lines can be excellent sources for the development of ShB resistant cultivars in rice.

Molecular analysis of broad-spectrum induced resistance in rice by the green leaf volatile Z-3-hexenyl acetate

Green leaf volatiles (GLVs), volatile organic compounds released by plants upon tissue damage, are key signaling molecules in plant immunity. The ability of exogenous GLV application to trigger an induced resistance (IR) phenotype against arthropod pests has been widely reported, but its effectiveness against plant pathogens is less well understood. In this study, we combined mRNA sequencing-based transcriptomics and phytohormone measurements with multispectral imaging-based precision phenotyping to gain insights into the molecular basis of Z-3-hexenyl acetate-induced resistance (Z-3-HAC-IR) in rice. Furthermore, we evaluated the efficacy of Z-3-HAC-IR against a panel of economically significant rice pathogens: Pyricularia oryzae, Rhizoctonia solani, Xanthomonas oryzae pv. oryzae, Cochliobolus miyabeanus, and Meloidogyne graminicola. Our data revealed rapid induction of jasmonate metabolism and systemic induction of plant immune responses upon Z-3-HAC exposure, as well as a transient allocation cost due to accelerated chlorophyll degradation and nutrient remobilization. Z-3-HAC-IR proved effective against all tested pathogens except for C. miyabeanus, including against the (hemi)biotrophs M. graminicola, X. oryzae pv. oryzae, and P. oryzae. The Z-3-HAC-IR phenotype was lost in the jasmonate (JA)-deficient hebiba mutant, which confirms the causal role of JA in Z-3-HAC-IR. Together, our results show that GLV exposure in rice induces broad-spectrum, JA-mediated disease resistance with limited allocation costs, and may thus be a promising alternative crop protection approach.

Genomic introgressions from African rice (Oryza glaberrima) in Asian rice (O. sativa) lead to the identification of key QTLs for panicle architecture

BACKGROUND

Developing high yielding varieties is a major challenge for breeders tackling the challenges of climate change in agriculture. The panicle (inflorescence) architecture of rice is one of the key components of yield potential and displays high inter- and intra-specific variability. The genus Oryza features two different crop species: Asian rice (Oryza sativa L.) and the African rice (O. glaberrima Steud.). One of the main morphological differences between the two independently domesticated species is the structure (or complexity) of the panicle, with O. sativa displaying a highly branched panicle, which in turn produces a larger number of grains than that of O. glaberrima. The gene regulatory network that governs intra- and interspecific panicle diversity is still under-studied.

RESULTS

To identify genetic factors linked to panicle architecture diversity in the two species, we used a set of 60 Chromosome Segment Substitution Lines (CSSLs) issued from third generation backcross (BC3DH) and carrying genomic segments from O. glaberrima cv. MG12 in the genetic background of O. sativa Tropical Japonica cv. Caiapó. Phenotypic data were collected for rachis and primary branch length, primary, secondary and tertiary branch number and spikelet number. A total of 15 QTLs were localized on chromosomes 1, 2, 3, 7, 11 and 12, QTLs associated with enhanced secondary and tertiary branch numbers were detected in two CSSLs. Furthermore, BC4F3:5 lines carrying different combinations of substituted segments were produced to decipher the effects of the identified QTL regions on variations in panicle architecture. A detailed analysis of phenotypes versus genotypes was carried out between the two parental genomes within these regions in order to understand how O. glaberrima introgression events may lead to alterations in panicle traits.

CONCLUSION

Our analysis led to the detection of genomic variations between O. sativa cv. Caiapó and O. glaberrima cv. MG12 in regions associated with enhanced panicle traits in specific CSSLs. These regions contain a number of key genes that regulate panicle development in O. sativa and their interspecific genomic variations may explain the phenotypic effects observed.

Inositolphosphorylceramide synthases, OsIPCSs, regulate plant height in rice

Inositolphosphorylceramide synthase (IPCS) catalyses ceramides and phosphatidylinositol (PI) into inositolphosphorylceramide (IPC), which is involved in the regulation of plant growth and development. A total of three OsIPCS family genes have been identified in rice. However, most of their functions remain unknown. Here, the functions of OsIPCSs were analyzed by CRISPR/Cas9 technology, lipidomics analysis, and transcriptomics analysis. Single-gene mutation of OsIPCSs resulted in dwarf phenotype. Among them, the phenotype of osipcs3 mutant was more severe. Multi-gene mutation of OsIPCS genes led to more severe phenotypes, indicating the additive effects of OsIPCSs. We further determined that a significant decrease in epidermal cell elongation of internode in the mutants. There was a significant decrease in the content of IPC detected in the osipcs2/3 and osipcs1/2/3 mutants. The contents of glycosyl inositol phosphoryl ceramide (GIPC) were also decreased by 20% and 10% in osipcs2/3 and osipcs1/2/3, respectively. The results of RNA-seq showed that numerous DEGs found to be associated with cellular component organization, anatomical structure morphogenesis, and cell growth in the osipcs2, osipcs2/3, and osipcs1/2/3. Taken together, OsIPCSs may be involved in the regulation of plant height through affecting cell growth and sphingolipid metabolism in rice.

Genome-Wide Identification and Analysis of Collar Region-Preferential Genes in Rice

The collar region plays a crucial role in leaf angle formation and plant architecture, which is important for improving crop yield given the challenges of diminishing arable land and changing environmental conditions. To determine collar region-preferential genes (CRPGs) affecting plant architecture and crop yield, we conducted genome-wide transcriptomic analysis. By integrating our RNA sequencing data with public rice anatomical expression data, we identified 657 CRPGs. Verification involved testing six randomly selected CRPGs, all of which exhibited collar-preferential expression. The functional significance of CRPGs was assessed via Gene Ontology enrichment analysis, utilizing MapMan and KEGG, and literature analysis provided additional information for characterized CRPGs. Our findings revealed links between manipulating leaf angle and phytohormone-related pathways and stress responses. Moreover, based on the CRPGs, five transcription factors downstream of the liguleless 1 (LG1) gene were identified. Overall, the identified CRPGs provide potential targets for further research and breeding applications aimed at improving crop productivity by manipulating leaf architecture.

Integrated omic analysis provides insights into the molecular regulation of stress tolerance by partial root-zone drying in rice

Partial root-zone drying (PRD) is an effective water-saving irrigation strategy that improves stress tolerance and facilitates efficient water use in several crops. It has long been considered that abscisic acid (ABA)-dependent drought resistance may be involved during partial root-zone drying. However, the molecular mechanisms underlying PRD-mediated stress tolerance remain unclear. It's hypothesized that other mechanisms might contribute to PRD-mediated drought tolerance. Here, rice seedlings were used as a research model and the complex transcriptomic and metabolic reprogramming processes were revealed during PRD, with several key genes involved in osmotic stress tolerance identified by using a combination of physiological, transcriptome, and metabolome analyses. Our results demonstrated that PRD induces transcriptomic alteration mainly in the roots but not in the leaves and adjusts several amino-acid and phytohormone metabolic pathways to maintain the balance between growth and stress response compared to the polyethylene glycol (PEG)-treated roots. Integrated analysis of the transcriptome and metabolome associated the co-expression modules with PRD-induced metabolic reprogramming. Several genes encoding the key transcription factors (TFs) were identified in these co-expression modules, highlighting several key TFs, including TCP19, WRI1a, ABF1, ABF2, DERF1, and TZF7, involved in nitrogen metabolism, lipid metabolism, ABA signaling, ethylene signaling, and stress regulation. Thus, our work presents the first evidence that molecular mechanisms other than ABA-mediated drought resistance are involved in PRD-mediated stress tolerance. Overall, our results provide new insights into PRD-mediated osmotic stress tolerance, clarify the molecular regulation induced by PRD, and identify genes useful for further improving water-use efficiency and/or stress tolerance in rice.

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.

Traces of Introgression from cAus into Tropical Japonica Observed in African Upland Rice Varieties

BACKGROUND

Asian rice Oryza sativa, first domesticated in East Asia, has considerable success in African fields. When and where this introduction occurred is unclear. Rice varieties of Asian origin may have evolved locally during and after migration to Africa, resulting in unique adaptations, particularly in relation to upland cultivation as frequently practiced in Africa.

METHODS

We investigated the genetic differentiation between Asian and African varieties using the 3000 Rice Genomes SNP dataset. African upland cultivars were first characterized using principal component analysis among 292 tropical Japonica accessions from Africa and Asia. The particularities of African accessions were then explored using two inference techniques, PCA-KDE for supervised classification and chromosome painting, and ELAI for individual allelic dosage monitoring.

KEY RESULTS

Ambiguities of local differentiation between Japonica and other groups pointed at genomic segments that potentially resulted from genetic exchange. Those specific to West African upland accessions were concentrated on chromosome 6 and featured several cAus introgression signals, including a large one between 17.9 and 21.7 Mb. We found iHS statistics in support of positive selection in this region and we provide a list of candidate genes enriched in GO terms that have regulatory functions involved in stress responses that could have facilitated adaptation to harsh upland growing conditions.

Genome-wide association study to identify genomic loci associated with early vigor in bread wheat under simulated water deficit complemented with quantitative trait loci meta-analysis

A genome-wide association study (GWAS) was used to identify associated loci with early vigor under simulated water deficit and grain yield under field drought in a diverse collection of Iranian bread wheat landraces. In addition, a meta-quantitative trait loci (MQTL) analysis was used to further expand our approach by retrieving already published quantitative trait loci (QTL) from recombinant inbred lines, double haploids, back-crosses, and F2 mapping populations. In the current study, around 16%, 14%, and 16% of SNPs were in significant linkage disequilibrium (LD) in the A, B, and D genomes, respectively, and varied between 5.44% (4A) and 21.85% (6A). Three main subgroups were identified among the landraces with different degrees of admixture, and population structure was further explored through principal component analysis. Our GWAS identified 54 marker-trait associations (MTAs) that were located across the wheat genome but with the highest number found in the B sub-genome. The gene ontology (GO) analysis of MTAs revealed that around 75% were located within or closed to protein-coding genes. In the MQTL analysis, 23 MQTLs, from a total of 215 QTLs, were identified and successfully projected onto the reference map. MQT-YLD4, MQT-YLD9, MQT-YLD13, MQT-YLD17, MQT-YLD18, MQT-YLD19, and MQTL-RL1 contributed to the highest number of projected QTLs and were therefore regarded as the most reliable and stable QTLs under water deficit conditions. These MQTLs greatly facilitate the identification of putative candidate genes underlying at each MQTL interval due to the reduced confidence of intervals associated with MQTLs. These findings provide important information on the genetic basis of early vigor traits and grain yield under water deficit conditions and set the foundation for future investigations into adaptation to water deficit in bread wheat.

Integrating GWAS and transcriptomics to identify candidate genes conferring heat tolerance in rice

INTRODUCTION

Rice (Oryza sativa L.) production is being challenged by global warming. Identifying new loci and favorable alleles associated with heat tolerance is crucial to developing rice heat-tolerant varieties.

METHODS

We evaluated the heat tolerance at the seedling stage using 620 diverse rice accessions. A total of six loci associated with heat tolerance were identified by a genome-wide association study (GWAS) with ~2.8 million single nucleotide polymorphisms (SNPs).

RESULTS

Among the six detected loci, qHT7 harbored the strongest association signal and the most associated SNPs. By comparing the transcriptomes of two representative accessions with contrasting heat tolerance, LOC_Os07g48710 (OsVQ30) was selected as a promising candidate gene in qHT7 due to the significant difference in its expression level between the two accessions. Haplotype 4 (Hap4) of LOC_Os07g48710 was determined as the favorable haplotype for heat tolerance via the gene-based haplotype analysis. The heat-tolerant haplotype LOC_Os07g48710Hap4 is highly enriched in the tropical Geng/Japonica accessions, and its frequency has decreased significantly during the improvement process of rice varieties.

DISCUSSION

Based on the GWAS and transcriptomics integrated results, a hypothetical model modulated by qHT7 in response to heat stress was proposed. Our results provide valuable candidate genes for improving rice heat tolerance through molecular breeding.