funRiceGenes dataset for comprehensive understanding and application of rice functional genes

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.

Innovative computational tools provide new insights into the polyploid wheat genome

Bread wheat (Triticum aestivum) is an important crop and serves as a significant source of protein and calories for humans, worldwide. Nevertheless, its large and allopolyploid genome poses constraints on genetic improvement. The complex reticulate evolutionary history and the intricacy of genomic resources make the deciphering of the functional genome considerably more challenging. Recently, we have developed a comprehensive list of versatile computational tools with the integration of statistical models for dissecting the polyploid wheat genome. Here, we summarize the methodological innovations and applications of these tools and databases. A series of step-by-step examples illustrates how these tools can be utilized for dissecting wheat germplasm resources and unveiling functional genes associated with important agronomic traits. Furthermore, we outline future perspectives on new advanced tools and databases, taking into consideration the unique features of bread wheat, to accelerate genomic-assisted wheat breeding.

Superior haplotypes of key drought-responsive genes reveal opportunities for the development of climate-resilient rice varieties

Haplotype-based breeding is an emerging and innovative concept that enables the development of designer crop varieties by exploiting and exploring superior alleles/haplotypes among target genes to create new traits in breeding programs. In this regard, whole-genome re-sequencing of 399 genotypes (landraces and breeding lines) from the 3000 rice genomes panel (3K-RG) is mined to identify the superior haplotypes for 95 drought-responsive candidate genes. Candidate gene-based association analysis reveals 69 marker-trait associations (MTAs) in 16 genes for single plant yield (SPY) under drought stress. Haplo-pheno analysis of these 16 genes identifies superior haplotypes for seven genes associated with the higher SPY under drought stress. Our study reveals that the performance of lines possessing superior haplotypes is significantly higher (p ≤ 0.05) as measured by single plant yield (SPY), for the OsGSK1-H4, OsDSR2-H3, OsDIL1-H22, OsDREB1C-H3, ASR3-H88, DSM3-H4 and ZFP182-H4 genes as compared to lines without the superior haplotypes. The validation results indicate that a superior haplotype for the DREB transcription factor (OsDREB1C) is present in all the drought-tolerant rice varieties, while it was notably absent in all susceptible varieties. These lines carrying the superior haplotypes can be used as potential donors in haplotype-based breeding to develop high-yielding drought-tolerant rice varieties.

Improving yield-related traits by editing the promoter of the heading date gene Ehd1 in rice

KEY MESSAGE

We developed an efficient promoter editing method to create different weak Ehd1 alleles in elite japonica rice variety ZJ8 with slightly delayed heading and improved yield for use in breeding. Heading date is an important agronomic trait of rice (Oryza sativa) that determines the planting areas and cultivation seasons of different varieties, thus affecting final yield. Early heading date 1 (Ehd1) is a major rice integrator gene in the regulatory network of heading date whose expression level is negatively correlated with heading date and grain yield. Some elite japonica varieties such as Zhongjia 8 (ZJ8) show very early heading with poor agronomic traits when planted in South China. This problem can be addressed by downregulating the expression of Ehd1. In this study, we analyzed the cis-regulatory elements in the Ehd1 promoter region. We then used CRISPR/Cas9-mediated editing to modify the Ehd1 promoter at multiple target sites in ZJ8. We rapidly identified homozygous allelic mutations in the T2 generation via long-read sequencing. We obtained several Ehd1 promoter mutants with different degrees of lower Ehd1 expression, delayed heading date, and improved yield-related traits. We developed an efficient promoter editing method to create different weak Ehd1 alleles for breeding selection. Using this method, a series of heading date materials from elite varieties can be created to expand the planting area of rice and improve grain yields.

Deciphering genetic basis of developmental and agronomic traits by integrating high-throughput optical phenotyping and genome-wide association studies in wheat

Dissecting the genetic basis of complex traits such as dynamic growth and yield potential is a major challenge in crops. Monitoring the growth throughout growing season in a large wheat population to uncover the temporal genetic controls for plant growth and yield-related traits has so far not been explored. In this study, a diverse wheat panel composed of 288 lines was monitored by a non-invasive and high-throughput phenotyping platform to collect growth traits from seedling to grain filling stage and their relationship with yield-related traits was further explored. Whole genome re-sequencing of the panel provided 12.64 million markers for a high-resolution genome-wide association analysis using 190 image-based traits and 17 agronomic traits. A total of 8327 marker-trait associations were detected and clustered into 1605 quantitative trait loci (QTLs) including a number of known genes or QTLs. We identified 277 pleiotropic QTLs controlling multiple traits at different growth stages which revealed temporal dynamics of QTLs action on plant development and yield production in wheat. A candidate gene related to plant growth that was detected by image traits was further validated. Particularly, our study demonstrated that the yield-related traits are largely predictable using models developed based on i-traits and provide possibility for high-throughput early selection, thus to accelerate breeding process. Our study explored the genetic architecture of growth and yield-related traits by combining high-throughput phenotyping and genotyping, which further unravelled the complex and stage-specific contributions of genetic loci to optimize growth and yield in wheat.

Multiple domestications of Asian rice

The origin of domesticated Asian rice (Oryza sativa L.) has been controversial for more than half a century. The debates have focused on two leading hypotheses: a single domestication event in China or multiple domestication events in geographically separate areas. These two hypotheses differ in their predicted history of genes/alleles selected during domestication. Here we amassed a dataset of 1,578 resequenced genomes, including an expanded sample of wild rice from throughout its geographic range. We identified 993 selected genes that generated phylogenetic trees on which japonica and indica formed a monophyletic group, suggesting that the domestication alleles of these genes originated only once in either japonica or indica. Importantly, the domestication alleles of most selected genes (~80%) stemmed from wild rice in China, but the domestication alleles of a substantial minority of selected genes (~20%) originated from wild rice in South and Southeast Asia, demonstrating separate domestication events of Asian rice.

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.

The translational landscape of bread wheat during grain development

The dynamics of gene expression in crop grains has typically been investigated at the transcriptional level. However, this approach neglects translational regulation, a widespread mechanism that rapidly modulates gene expression to increase the plasticity of organisms. Here, we performed ribosome profiling and polysome profiling to obtain a comprehensive translatome data set of developing bread wheat (Triticum aestivum) grains. We further investigated the genome-wide translational dynamics during grain development, revealing that the translation of many functional genes is modulated in a stage-specific manner. The unbalanced translation between subgenomes is pervasive, which increases the expression flexibility of allohexaploid wheat. In addition, we uncovered widespread previously unannotated translation events, including upstream open reading frames (uORFs), downstream open reading frames (dORFs), and open reading frames (ORFs) in long noncoding RNAs, and characterized the temporal expression dynamics of small ORFs. We demonstrated that uORFs act as cis-regulatory elements that can repress or even enhance the translation of mRNAs. Gene translation may be combinatorially modulated by uORFs, dORFs, and microRNAs. In summary, our study presents a translatomic resource that provides a comprehensive and detailed overview of the translational regulation in developing bread wheat grains. This resource will facilitate future crop improvements for optimal yield and quality.

Mining genes regulating root system architecture in maize based on data integration analysis

KEY MESSAGE

A new strategy that integrated multiple public data resources was established to construct root gene co-expression network and mine genes regulating root system architecture in maize. A root gene co-expression network, containing 13,874 genes, was constructed. A total of 53 root hub genes and 16 priority root candidate genes were identified. One priority root candidate was further functionally verified using overexpression transgenic maize lines. Root system architecture (RSA) is crucial for crops productivity and stress tolerance. In maize, few RSA genes are functionally cloned, and effective discovery of RSA genes remains a great of challenge. In this work, we established a strategy to mine maize RSA genes by integrating functionally characterized root genes, root transcriptome, weighted gene co-expression network analysis (WGCNA) and genome-wide association analysis (GWAS) of RSA traits based on public data resources. A total of 589 maize root genes were collected by searching well-characterized root genes in maize or homologous genes of other species. We performed WGCNA to construct a maize root gene co-expression network containing 13874 genes based on public available root transcriptome data, and further discovered the 53 hub genes related to root traits. In addition, by the prediction function of obtained root gene co-expression network, a total of 1082 new root candidate genes were explored. By further overlapping the obtained new root candidate gene with the root-related GWAS of RSA candidate genes, 16 priority root candidate genes were identified. Finally, a priority root candidate gene, Zm00001d023379 (encodes pyruvate kinase 2), was validated to modulate root open angle and shoot-borne roots number using its overexpression transgenic lines. Our results develop an integration analysis method for effectively exploring regulatory genes of RSA in maize and open a new avenue to mine the candidate genes underlying complex traits.

Wheat EARLY FLOWERING 3 affects heading date without disrupting circadian oscillations

Plant breeders have indirectly selected for variation at circadian-associated loci in many of the world's major crops, when breeding to increase yield and improve crop performance. Using an eight-parent Multiparent Advanced Generation Inter-Cross (MAGIC) population, we investigated how variation in circadian clock-associated genes contributes to the regulation of heading date in UK and European winter wheat (Triticum aestivum) varieties. We identified homoeologues of EARLY FLOWERING 3 (ELF3) as candidates for the Earliness per se (Eps) D1 and B1 loci under field conditions. We then confirmed a single-nucleotide polymorphism within the coding region of TaELF3-B1 as a candidate polymorphism underlying the Eps-B1 locus. We found that a reported deletion at the Eps-D1 locus encompassing TaELF3-D1 is, instead, an allele that lies within an introgression region containing an inversion relative to the Chinese Spring D genome. Using Triticum turgidum cv. Kronos carrying loss-of-function alleles of TtELF3, we showed that ELF3 regulates heading, with loss of a single ELF3 homoeologue sufficient to alter heading date. These studies demonstrated that ELF3 forms part of the circadian oscillator; however, the loss of all homoeologues was required to affect circadian rhythms. Similarly, loss of functional LUX ARRHYTHMO (LUX) in T. aestivum, an orthologue of a protein partner of Arabidopsis (Arabidopsis thaliana) ELF3, severely disrupted circadian rhythms. ELF3 and LUX transcripts are not co-expressed at dusk, suggesting that the structure of the wheat circadian oscillator might differ from that of Arabidopsis. Our demonstration that alterations to ELF3 homoeologues can affect heading date separately from effects on the circadian oscillator suggests a role for ELF3 in cereal photoperiodic responses that could be selected for without pleiotropic deleterious alterations to circadian rhythms.

Genomics and transcriptomics to protect rice (Oryza sativa. L.) from abiotic stressors: -pathways to achieving zero hunger

More over half of the world's population depends on rice as a major food crop. Rice (Oryza sativa L.) is vulnerable to abiotic challenges including drought, cold, and salinity since it grown in semi-aquatic, tropical, or subtropical settings. Abiotic stress resistance has bred into rice plants since the earliest rice cultivation techniques. Prior to the discovery of the genome, abiotic stress-related genes were identified using forward genetic methods, and abiotic stress-tolerant lines have developed using traditional breeding methods. Dynamic transcriptome expression represents the degree of gene expression in a specific cell, tissue, or organ of an individual organism at a specific point in its growth and development. Transcriptomics can reveal the expression at the entire genome level during stressful conditions from the entire transcriptional level, which can be helpful in understanding the intricate regulatory network relating to the stress tolerance and adaptability of plants. Rice (Oryza sativa L.) gene families found comparatively using the reference genome sequences of other plant species, allowing for genome-wide identification. Transcriptomics via gene expression profiling which have recently dominated by RNA-seq complements genomic techniques. The identification of numerous important qtl,s genes, promoter elements, transcription factors and miRNAs involved in rice response to abiotic stress was made possible by all of these genomic and transcriptomic techniques. The use of several genomes and transcriptome methodologies to comprehend rice (Oryza sativa, L.) ability to withstand abiotic stress have been discussed in this review.

Improving Yield and Yield Stability in Winter Rye by Hybrid Breeding

Rye is the only cross-pollinating small-grain cereal. The unique reproduction biology results in an exceptional complexity concerning genetic improvement of rye by breeding. Rye is a close relative of wheat and has a strong adaptation potential that refers to its mating system, making this overlooked cereal readily adjustable to a changing environment. Rye breeding addresses the emerging challenges of food security associated with climate change. The systematic identification, management, and use of its valuable natural diversity became a feasible option in outbreeding rye only following the establishment of hybrid breeding late in the 20th century. In this article, we review the most recent technological advances to improve yield and yield stability in winter rye. Based on recently released reference genome sequences, SMART breeding approaches are described to counterbalance undesired linkage drag effects of major restorer genes on grain yield. We present the development of gibberellin-sensitive semidwarf hybrids as a novel plant breeding innovation based on an approach that is different from current methods of increasing productivity in rye and wheat. Breeding of new rye cultivars with improved performance and resilience is indispensable for a renaissance of this healthy minor cereal as a homogeneous commodity with cultural relevance in Europe that allows for comparatively smooth but substantial complementation of wheat with rye-based diets, supporting the necessary restoration of the balance between human action and nature.

Overexpression of the WAPO-A1 gene increases the number of spikelets per spike in bread wheat

Two homoeologous QTLs for number of spikelets per spike (SPS) were mapped on chromosomes 7AL and 7BL using two wheat MAGIC populations. Sets of lines contrasting for the QTL on 7AL were developed which allowed for the validation and fine mapping of the 7AL QTL and for the identification of a previously described candidate gene, WHEAT ORTHOLOG OF APO1 (WAPO1). Using transgenic overexpression in both a low and a high SPS line, we provide a functional validation for the role of this gene in determining SPS also in hexaploid wheat. We show that the expression levels of this gene positively correlate with SPS in multiple MAGIC founder lines under field conditions as well as in transgenic lines grown in the greenhouse. This work highlights the potential use of WAPO1 in hexaploid wheat for further yield increases. The impact of WAPO1 and SPS on yield depends on other genetic and environmental factors, hence, will require a finely balanced expression level to avoid the development of detrimental pleiotropic phenotypes.

The genetic basis for panicle trait variation in switchgrass (Panicum virgatum)

We investigate the genetic basis of panicle architecture in switchgrass in two mapping populations across a latitudinal gradient, and find many stable, repeatable genetic effects and limited genetic interactions with the environment. Grass species exhibit large diversity in panicle architecture influenced by genes, the environment, and their interaction. The genetic study of panicle architecture in perennial grasses is limited. In this study, we evaluate the genetic basis of panicle architecture including panicle length, primary branching number, and secondary branching number in an outcrossed switchgrass QTL population grown across ten field sites in the central USA through multi-environment mixed QTL analysis. We also evaluate genetic effects in a diversity panel of switchgrass grown at three of the ten field sites using genome-wide association (GWAS) and multivariate adaptive shrinkage. Furthermore, we search for candidate genes underlying panicle traits in both of these independent mapping populations. Overall, 18 QTL were detected in the QTL mapping population for the three panicle traits, and 146 unlinked genomic regions in the diversity panel affected one or more panicle trait. Twelve of the QTL exhibited consistent effects (i.e., no QTL by environment interactions or no QTL × E), and most (four of six) of the effects with QTL × E exhibited site-specific effects. Most (59.3%) significant partially linked diversity panel SNPs had significant effects in all panicle traits and all field sites and showed pervasive pleiotropy and limited environment interactions. Panicle QTL co-localized with significant SNPs found using GWAS, providing additional power to distinguish between true and false associations in the diversity panel.

Growing Maize Root: Lectins Involved in Consecutive Stages of Cell Development

Proteins that carry specific carbohydrate-binding lectin domains have a great variety and are ubiquitous across the plant kingdom. In turn, the plant cell wall has a complex carbohydrate composition, which is subjected to constant changes in the course of plant development. In this regard, proteins with lectin domains are of great interest in the context of studying their contribution to the tuning and monitoring of the cell wall during its modifications in the course of plant organ development. We performed a genome-wide screening of lectin motifs in the Zea mays genome and analyzed the transcriptomic data from five zones of primary maize root with cells at different development stages. This allowed us to obtain 306 gene sequences encoding putative lectins and to relate their expressions to the stages of root cell development and peculiarities of cell wall metabolism. Among the lectins whose expression was high and differentially regulated in growing maize root were the members of the EUL, dirigent–jacalin, malectin, malectin-like, GNA and Nictaba families, many of which are predicted as cell wall proteins or lectin receptor-like kinases that have direct access to the cell wall. Thus, a set of molecular players was identified with high potential to play important roles in the early stages of root morphogenesis.

The Triticum ispahanicum elongated glume locus P2 maps to chromosome 6A and is associated with the ectopic expression of SVP-A1

We propose the MADS-box transcription factor SVP-A1 as a promising candidate gene for the elongated glume locus P2, which maps to chromosome 6A instead of the previously proposed chromosome 7B. In rice and wheat, glume and floral organ length are positively correlated with grain size, making them an important target to increase grain size and potentially yield. The wheat subspecies Triticum ispahanicum is known to develop elongated glumes and floral organs as well as long grains. These multiple phenotypic effects are controlled by the P2 locus, which was previously mapped to wheat chromosome 7B. Using three mapping populations, we show that the long glume locus P2 does not map to chromosome 7B, but instead maps to a 1.68 Mbp interval on chromosome 6A. Within this interval, we identified SVP-A1, a MADS box transcription factor which is the direct ortholog of the maize gene underlying the 'pod corn' Tunicate locus and is a paralog to the T. polonicum elongated glume P1 gene. In T. ispahanicum, we identified a unique allele which has a 482-bp deletion in the SVP-A1 promoter and is associated with ectopic and higher expression of SVP-A1 in the elongated glumes and floral organs. We used near-isogenic lines (NILs) to show that P2 has a consistent positive effect on the length of glume, lemma, palea, spike and grain. Based on the mapping data, natural variation, biological function of SVP genes in cereals and expression analyses, we propose the MADS-box transcription factor SVP-A1 as a promising candidate for P2.

Transcriptome profiling of two super hybrid rice provides insights into the genetic basis of heterosis

BACKGROUND

Heterosis is a phenomenon that hybrids show superior performance over their parents. The successful utilization of heterosis has greatly improved rice productivity, but the molecular basis of heterosis remains largely unclear.

RESULTS

Here, the transcriptomes of young panicles and leaves of the two widely grown two-line super hybrid rice varieties (Jing-Liang-You-Hua-Zhan (JLYHZ) and Long-Liang-You-Hua-Zhan (LLYHZ)) and their parents were analyzed by RNA-seq. Transcriptome profiling of the hybrids revealed 1,778 ~ 9,404 differentially expressed genes (DEGs) in two tissues, which were identified by comparing with their parents. GO, and KEGG enrichment analysis showed that the pathways significantly enriched in both tissues of two hybrids were all related to yield and resistance, like circadian rhythm (GO:0,007,623), response to water deprivation (GO:0,009,414), and photosynthetic genes (osa00196). Allele-specific expression genes (ASEGs) were also identified in hybrids. The ASEGs were most significantly enriched in ionotropic glutamate receptor signaling pathway, which was hypothesized to be potential amino acid sensors in plants. Moreover, the ASEGs were also differentially expressed between parents. The number of variations in ASEGs is higher than expected, especially for large effect variations. The DEGs and ASEGs are the potential reasons for the formation of heterosis in the two elite super hybrid rice.

CONCLUSIONS

Our results provide a comprehensive understanding of the heterosis of two-line super hybrid rice and facilitate the exploitation of heterosis in hybrid rice breeding with high yield heterosis.

Genetic Analysis of a Collection of Rice Germplasm (Oryza sativa L.) through High-Density SNP Array Provides Useful Information for Further Breeding Practices

Traditional breeding strategies mainly focus on the evaluation of trait performance, but pay less attention to the changing genetic background. A comprehensive understanding of the genetic diversity of germplasms is crucial for the deliberate improvement of specific traits. A collection of 154 highland rice varieties were collected as the initial genetic resource in our breeding program to improve the pathogen resistance and eating and cooking qualities. These varieties were analyzed using a whole-genome SNP array and were clustered into three groups. Further analysis revealed that the favorable alleles of pathogen resistance genes are mostly absent in our collected varieties. However, it showed that most varieties possess favorable alleles of Waxy (Wx) and ALKALI DEGENERATION (ALK), which are able to enhance the eating and cooking qualities. Moreover, only about one fifth of all varieties harbors favorable the allele of fragrance gene Betainealdehyde dehydrogenase (BADH2). Together, these results give an overall view of the genetic constitution of the target traits, which provide useful information for future genetic improvement in breeding practices.

Comparative transcriptome analysis of pollen and anther wall reveals novel insights into the regulatory mechanisms underlying anther wall development and its dehiscence in rice

To further understand the regulatory mechanism for anther dehiscence in rice, we carried out transcriptome analysis for the following two tissues: the anther wall and pollen at the anthesis stage. With the anatomical meta-expression data, in addition to these tissues, the differentially expressed genes (DEGs) between the two tissues were further refined to identify 1,717 pollen-preferred genes and 534 anther wall-preferred genes. A GUS transgenic line and RT-qPCR analysis for anther wall-preferred genes supported the fidelity of our gene candidates for further analysis. The refined DEGs were functionally classified through Gene Ontology (GO) enrichment and MapMan analyses. Through the analysis of cis-acting elements and alternative splicing variants, we also suggest the feature of regulatory sequences in promoter regions for anther wall-preferred expression and provide information of the unique splicing variants in anther wall. Subsequently, it was found that hormone signaling and the resulting transcriptional regulation pathways may play an important role in anther dehiscence and anther wall development. Our results could provide useful insights into future research to broaden the molecular mechanism of anther dehiscence or anther wall development in rice.

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

As a major food crop and model organism, rice has been mostly studied with the largest number of functionally characterized genes among all crops. We previously built the funRiceGenes database including ~ 2800 functionally characterized rice genes and ~ 5000 members of different gene families. Since being published, the funRiceGenes database has been accessed by more than 54,400 users with over 540,000 pageviews. The funRiceGenes database has been continuously updated with newly cloned rice genes and newly published literature, based on the progress of rice functional genomics studies. Up to Nov 2021, ~ 4100 functionally characterized rice genes and ~ 6000 members of different gene families were collected in funRiceGenes, accounting for 22.3% of the 39,045 annotated protein-coding genes in the rice genome. Here, we summarized the update of the funRiceGenes database with new data and new features in the last 5 years.

Combining Genome-Wide Association Study and Gene-Based Haplotype Analysis to Identify Candidate Genes for Alkali Tolerance at the Germination Stage in Rice

Salinity-alkalinity 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 the tolerance to salinity-alkalinity of rice varieties at the germination stage. Although we have a more comprehensive understanding of salt tolerance in rice, the genetic basis of alkali tolerance in rice is still poorly understood. In this study, we measured seven germination-related traits under alkali stress and control conditions using 428 diverse rice accessions. The alkali tolerance levels of rice germplasms varied considerably during germination. Xian/indica accessions had generally higher tolerance to alkali stress than Geng/japonica accessions at the germination stage. Using genome-wide association analysis, 90 loci were identified as significantly associated with alkali tolerance. Eight genes (LOC_Os01g12000, LOC_Os03g60240, LOC_Os03g08960, LOC_Os04g41410, LOC_Os09g25060, LOC_Os11g35350, LOC_Os12g09350, and LOC_Os12g13300) were selected as important candidate genes for alkali tolerance based on the gene functional annotation and gene-CDS-haplotype analysis. According to the expression levels of LOC_Os09g25060 (OsWRKY76), it is likely to play a negative regulatory role in alkali tolerance during rice germination. An effective strategy for improving rice alkali tolerance may be to pyramid alkali-tolerant haplotypes of multiple candidate genes to obtain the optimal haplotype combination. Our findings may provide valuable genetic information and expand the use of alkali tolerance germplasm resources in rice molecular breeding to improve the alkali tolerance at the germination stage.

Genome-wide association study and genomic prediction for yield and grain quality traits of hybrid rice

Genomic selection is an efficient tool for breeding selection, especially for quantitative traits controlled by multiples genes with low heritability. To validate the application of genomic selection in hybrid rice breeding, the yield and grain quality traits of 404 hybrid rice breeding lines were investigated, and the same accessions were genotyped by using a 56 K SNP chip. There were wide variances among the tested accessions for all the measured traits, and most of the traits were correlated. A total of 67 significant loci were identified for the yield-related traits, and 123 significant loci were identified for the grain quality traits by GWAS. Two of these loci associated with increasing grain yield but decreasing grain quality. The GEBVs of all the yield and grain quality traits were calculated by using 15 different prediction algorithms. The plant height, panicle length, thousand grain weight, grain length and width ratio, amylose content, and alkali value have higher predictability than other traits. However, the predictive accuracy of different GS models is different for different traits. This study provided useful information for genomic selection of specific trait using proper markers and prediction models.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01289-6.

Grain dispersal mechanism in cereals arose from a genome duplication followed by changes in spatial expression of genes involved in pollen development

Grain disarticulation in wild progenitor of wheat and barley evolved through a local duplication event followed by neo-functionalization resulting from changes in location of gene expression. One of the most critical events in the process of cereal domestication was the loss of the natural mode of grain dispersal. Grain dispersal in barley is controlled by two major genes, Btr1 and Btr2, which affect the thickness of cell walls around the disarticulation zone. The barley genome also encodes Btr1-like and Btr2-like genes, which have been shown to be the ancestral copies. While Btr and Btr-like genes are non-redundant, the biological function of Btr-like genes is unknown. We explored the potential biological role of the Btr-like genes by surveying their expression profile across 212 publicly available transcriptome datasets representing diverse organs, developmental stages and stress conditions. We found that Btr1-like and Btr2-like are expressed exclusively in immature anther samples throughout Prophase I of meiosis within the meiocyte. The similar and restricted expression profile of these two genes suggests they are involved in a common biological function. Further analysis revealed 141 genes co-expressed with Btr1-like and 122 genes co-expressed with Btr2-like, with 105 genes in common, supporting Btr-like genes involvement in a shared molecular pathway. We hypothesize that the Btr-like genes play a crucial role in pollen development by facilitating the formation of the callose wall around the meiocyte or in the secretion of callase by the tapetum. Our data suggest that Btr genes retained an ancestral function in cell wall modification and gained a new role in grain dispersal due to changes in their spatial expression becoming spike specific after gene duplication.

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.

A FLASH pipeline for arrayed CRISPR library construction and the gene function discovery of rice receptor-like kinases

Clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9)-mediated gene editing is revolutionizing plant research and crop breeding. Here, we present an effective and streamlined pipeline for arrayed CRISPR library construction and demonstrate it is suitable for small- to large-scale genome editing in plants. This pipeline introduces artificial PCR fragment-length markers for distinguishing guide RNAs (gRNAs) (FLASH), and a group of 12 constructs harboring different FLASH tags are co-transformed into plants each time. The identities of gRNAs in Agrobacterium mixtures and transgenic plants can therefore be read out by detecting the FLASH tags, a process that requires only conventional PCR and gel electrophoresis rather than sequencing. We generated an arrayed CRISPR library targeting all 1,072 members of the receptor-like kinase (RLK) family in rice. One-shot transformation generated a mutant population that covers gRNAs targeting 955 RLKs, and 74.3% (710/955) of the target genes had three or more independent T₀ lines. Our results indicate that the FLASH tags act as bona fide surrogates for the gRNAs and are tightly (92.1%) associated with frameshift mutations in the target genes. In addition, the FLASH pipeline allows for rapid identification of unintended editing events without corresponding T-DNA integrations and generates high-order mutants of closely related RLK genes. Furthermore, we showed that the RLK mutant library enables rapid discovery of defense-related RLK genes. This study introduces an effective pipeline for arrayed CRISPR library construction and provides genome-wide rice RLK mutant resources for functional genomics.

The genome of the rice variety LTH provides insight into its universal susceptibility mechanism to worldwide rice blast fungal strains

The widely used rice variety Lijiangxintuanheigu (LTH) shows a universal susceptibility to thousands of Magnaporthe oryzae isolates, the causal agent of devastating rice blast, making LTH an ideal line in resistance (R) gene cloning. However, the underlying genetic mechanism of the universal susceptibility has not been fully revealed because of the lack of a high-quality genome. Here, we took a genomic approach together with experimental assays to investigate LTH’s universal susceptibility to rice blast. Using Nanopore long reads, we assembled a chromosome-level genome. Millions of genomic variants were detected by comparing LTH with 10 other rice varieties, of which large-effect variants could affect plant immunity. Gene family analyses show that the number of R genes and leucine-rich repeat receptor-like protein kinase (LRR-RLK)-encoding genes decrease significantly in LTH. Rice blast resistance genes called Pi genes are either absent or disrupted by genomic variations. Additionally, residual R genes of LTH are likely under weak pathogen selection pressure, and other plant defense-related genes are weakly induced by rice blast. In contrast, the pattern-triggered immunity (PTI) of LTH is normal, as demonstrated by experimental assays. Therefore, we conclude that weak effector-trigger immunity (ETI)-mediated primarily by Pi genes but not PTI results in the universal susceptibility of LTH to rice blast. The attenuated ETI of LTH may be also associated with reduced numbers of R genes and LRR-RLKs, and minimally functional residual defense-related genes. Finally, we demonstrate the use of the LTH genome by rapid cloning of the Pi gene Piak from a resistant variety.

An Efficient Marker Gene Excision Strategy Based on CRISPR/Cas9-Mediated Homology-Directed Repair in Rice

In order to separate transformed cells from non-transformed cells, antibiotic selectable marker genes are usually utilized in genetic transformation. After obtaining transgenic plants, it is often necessary to remove the marker gene from the plant genome in order to avoid regulatory issues. However, many marker-free systems are time-consuming and labor-intensive. Homology-directed repair (HDR) is a process of homologous recombination using homologous arms for efficient and precise repair of DNA double-strand breaks (DSBs). The clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein-9 (Cas9) system is a powerful genome editing tool that can efficiently cause DSBs. Here, we isolated a rice promoter (Pssi) of a gene that highly expressed in stem, shoot tip and inflorescence, and established a high-efficiency sequence-excision strategy by using this Pssi to drive CRISPR/Cas9-mediated HDR for marker free (PssiCHMF). In our study, PssiCHMF-induced marker gene deletion was detected in 73.3% of T₀ plants and 83.2% of T₁ plants. A high proportion (55.6%) of homozygous marker-excised plants were obtained in T₁ progeny. The recombinant GUS reporter-aided analysis and its sequencing of the recombinant products showed precise deletion and repair mediated by the PssiCHMF method. In conclusion, our CRISPR/Cas9-mediated HDR auto-excision method provides a time-saving and efficient strategy for removing the marker genes from transgenic plants.

Identification and Validation of Quantitative Trait Loci Mapping for Spike-Layer Uniformity in Wheat

Spike-layer uniformity (SLU), the consistency of the spike distribution in the vertical space, is an important trait. It directly affects the yield potential and appearance. Revealing the genetic basis of SLU will provide new insights into wheat improvement. To map the SLU-related quantitative trait loci (QTL), 300 recombinant inbred lines (RILs) that were derived from a cross between H461 and Chinese Spring were used in this study. The RILs and parents were tested in fields from two continuous years from two different pilots. Phenotypic analysis showed that H461 was more consistent in the vertical spatial distribution of the spike layer than in Chinese Spring. Based on inclusive composite interval mapping, four QTL were identified for SLU. There were two major QTL on chromosomes 2BL and 2DL and two minor QTL on chromosomes 1BS and 2BL that were identified. The additive effects of QSlu.sicau-1B, Qslu.sicau-2B-2, and QSlu.sicau-2D were all from the parent, H461. The major QTL, QSlu.sicau-2B-2 and QSlu.sicau-2D, were detected in each of the conducted trials. Based on the best linear unbiased prediction values, the two loci explained 23.97% and 15.98% of the phenotypic variation, respectively. Compared with previous studies, the two major loci were potentially novel and the two minor loci were overlapped. Based on the kompetitive allele-specific PCR (KASP) marker, the genetic effects for QSlu.sicau-2B-2 were validated in an additional RIL population. The genetic effects ranged from 26.65% to 32.56%, with an average value of 30.40%. In addition, QSlu.sicau-2B-2 showed a significant (p < 0.01) and positive influence on the spike length, spikelet number, and thousand kernel weight. The identified QTL and the developed KASP marker will be helpful for fine-mapping these loci, finally contributing to wheat breeding programs in a marker-assisted selection way.

From Green Super Rice to green agriculture: Reaping the promise of functional genomics research

Producing sufficient food with finite resources to feed the growing global population while having a smaller impact on the environment has always been a great challenge. Here, we review the concept and practices of Green Super Rice (GSR) that have led to a paradigm shift in goals for crop genetic improvement and models of food production for promoting sustainable agriculture. The momentous achievements and global deliveries of GSR have been fueled by the integration of abundant genetic resources, functional gene discoveries, and innovative breeding techniques with precise gene and whole-genome selection and efficient agronomic management to promote resource-saving, environmentally friendly crop production systems. We also provide perspectives on new horizons in genomic breeding technologies geared toward delivering green and nutritious crop varieties to further enhance the development of green agriculture and better nourish the world population.

CG and CHG Methylation Contribute to the Transcriptional Control of OsPRR37-Output Genes in Rice

Plant circadian clock coordinates endogenous transcriptional rhythms with diurnal changes of environmental cues. OsPRR37, a negative component in the rice circadian clock, reportedly regulates transcriptome rhythms, and agronomically important traits. However, the underlying regulatory mechanisms of OsPRR37-output genes remain largely unknown. In this study, whole genome bisulfite sequencing and high-throughput RNA sequencing were applied to verify the role of DNA methylation in the transcriptional control of OsPRR37-output genes. We found that the overexpression of OsPRR37 suppressed rice growth and altered cytosine methylations in CG and CHG sequence contexts in but not the CHH context (H represents A, T, or C). In total, 35 overlapping genes were identified, and 25 of them showed negative correlation between the methylation level and gene expression. The promoter of the hexokinase gene OsHXK1 was hypomethylated at both CG and CHG sites, and the expression of OsHXK1 was significantly increased. Meanwhile, the leaf starch content was consistently lower in OsPRR37 overexpression lines than in the recipient parent Guangluai 4. Further analysis with published data of time-course transcriptomes revealed that most overlapping genes showed peak expression phases from dusk to dawn. The genes involved in DNA methylation, methylation maintenance, and DNA demethylation were found to be actively expressed around dusk. A DNA glycosylase, namely ROS1A/DNG702, was probably the upstream candidate that demethylated the promoter of OsHXK1. Taken together, our results revealed that CG and CHG methylation contribute to the transcriptional regulation of OsPRR37-output genes, and hypomethylation of OsHXK1 leads to decreased starch content and reduced plant growth in rice.

Transcriptome Analysis of Triple Mutant for OsMADS62, OsMADS63, and OsMADS68 Reveals the Downstream Regulatory Mechanism for Pollen Germination in Rice (Oryza sativa)

The MADS (MCM1-AGAMOUS-DEFFICIENS-SRF) gene family has a preserved domain called MADS-box that regulates downstream gene expression as a transcriptional factor. Reports have revealed three MADS genes in rice, OsMADS62, OsMADS63, and OsMADS68, which exhibits preferential expression in mature rice pollen grains. To better understand the transcriptional regulation of pollen germination and tube growth in rice, we generated the loss-of-function homozygous mutant of these three OsMADS genes using the CRISPR-Cas9 (clustered regularly interspaced short palindromic repeats-CRISPR associated protein 9) system in wild-type backgrounds. Results showed that the triple knockout (KO) mutant showed a complete sterile phenotype without pollen germination. Next, to determine downstream candidate genes that are transcriptionally regulated by the three OsMADS genes during pollen development, we proceeded with RNA-seq analysis by sampling the mature anther of the mutant and wild-type. Two hundred and seventy-four upregulated and 658 downregulated genes with preferential expressions in the anthers were selected. Furthermore, downregulated genes possessed cell wall modification, clathrin coat assembly, and cellular cell wall organization features. We also selected downregulated genes predicted to be directly regulated by three OsMADS genes through the analyses for promoter sequences. Thus, this study provides a molecular background for understanding pollen germination and tube growth mediated by OsMADS62, OsMADS63, and OsMADS68 with mature pollen preferred expression.

CRISPR-Cas9 mediated mutation in GRAIN WIDTH and WEIGHT2 (GW2) locus improves aleurone layer and grain nutritional quality in rice

Enhancing crop productivity and their nutritional quality are the key components and primary focus of crop improvement strategy for fulfilling future food demand and improving human health. Grain filling and endosperm development are the key determinants of grain yield and nutritional quality. GRAIN WIDTH and WEIGHT2 (GW2) gene encodes a RING-type E3 ubiquitin ligase and determines the grain weight in cereal crops. Here we report GW2 knockout (KO) mutants in Indica (var. MTU1010) through CRISPR/Cas9 genome editing. The endosperm of GW2-KO mutant seed displays a thick aleurone layer with enhanced grain protein content. Further the loss of function of OsGW2 results in improved accumulation of essential dietary minerals (Fe, Zn, K, P, Ca) in the endosperm of rice grain. Additionally, the mutants displayed an early growth vigour phenotype with an improved root and shoot architecture. The hull morphology of GW2-KO lines also showed improved, grain filling thereby promoting larger grain architecture. Together, our findings indicate that GW2 may serve as a key regulator of improved grain architecture, grain nutritional quality and an important modulator of plant morphology. The study offers a strategy for the development of improved rice cultivars with enriched nutritional quality and its possible implementation in other cereals as well.

A Genome-Wide Association Study Pinpoints Quantitative Trait Genes for Plant Height, Heading Date, Grain Quality, and Yield in Rye (Secale cereale L.)

Rye is the only cross-pollinating Triticeae crop species. Knowledge of rye genes controlling complex-inherited traits is scarce, which, currently, largely disables the genomics assisted introgression of untapped genetic variation from self-incompatible germplasm collections in elite inbred lines for hybrid breeding. We report on the first genome-wide association study (GWAS) in rye based on the phenotypic evaluation of 526 experimental hybrids for plant height, heading date, grain quality, and yield in 2 years and up to 19 environments. We established a cross-validated NIRS calibration model as a fast, effective, and robust analytical method to determine grain quality parameters. We observed phenotypic plasticity in plant height and tiller number as a resource use strategy of rye under drought and identified increased grain arabinoxylan content as a striking phenotype in osmotically stressed rye. We used DArTseq™ as a genotyping-by-sequencing technology to reduce the complexity of the rye genome. We established a novel high-density genetic linkage map that describes the position of almost 19k markers and that allowed us to estimate a low genome-wide LD based on the assessed genetic diversity in elite germplasm. We analyzed the relationship between plant height, heading date, agronomic, as well as grain quality traits, and genotype based on 20k novel single-nucleotide polymorphism markers. In addition, we integrated the DArTseq™ markers in the recently established 'Lo7' reference genome assembly. We identified cross-validated SNPs in 'Lo7' protein-coding genes associated with all traits studied. These include associations of the WUSCHEL-related homeobox transcription factor DWT1 and grain yield, the DELLA protein gene SLR1 and heading date, the Ethylene overproducer 1-like protein gene ETOL1 and thousand-grain weight, protein and starch content, as well as the Lectin receptor kinase SIT2 and plant height. A Leucine-rich repeat receptor protein kinase and a Xyloglucan alpha-1,6-xylosyltransferase count among the cross-validated genes associated with water-extractable arabinoxylan content. This study demonstrates the power of GWAS, hybrid breeding, and the reference genome sequence in rye genetics research to dissect and identify the function of genes shaping genetic diversity in agronomic and grain quality traits of rye. The described links between genetic causes and phenotypic variation will accelerate genomics-enabled rye improvement.

Low Temperature Antioxidant Activity QTL Associate with Genomic Regions Involved in Physiological Cold Stress Tolerance Responses in Rice (Oryza sativa L.)

Boosting cold stress tolerance in crop plants can minimize stress-mediated yield losses. Asian rice (Oryza sativa L.), one of the most consumed cereal crops, originated from subtropical regions and is generally sensitive to low temperature environments. Within the two subspecies of rice, JAPONICA, and INDICA, the cold tolerance potential of its accessions is highly variable and depends on their genetic background. Yet, cold stress tolerance response mechanisms are complex and not well understood. This study utilized 370 accessions from the Rice Diversity Panel 1 (RDP1) to investigate and correlate four cold stress tolerance response phenotypes: membrane damage, seedling survivability, and catalase and anthocyanin antioxidative activity. Most JAPONICA accessions, and admixed accessions within JAPONICA, had lower membrane damage, higher antioxidative activity, and overall, higher seedling survivability compared to INDICA accessions. Genome-wide association study (GWAS) mapping was done using the four traits to find novel quantitative trait loci (QTL), and to validate and fine-map previously identified QTL. A total of 20 QTL associated to two or more traits were uncovered by our study. Gene Ontology (GO) term enrichment analyses satisfying four layers of filtering retrieved three potential pathways: signal transduction, maintenance of plasma membrane and cell wall integrity, and nucleic acids metabolism as general mechanisms of cold stress tolerance responses involving antioxidant activity.

Mapping causal genes and genetic interactions for agronomic traits using a large F2 population in rice

Dissecting the genetic mechanisms underlying agronomic traits is of great importance for crop breeding. Agronomic traits are usually controlled by multiple quantitative trait loci (QTLs) and genetic interactions, and mapping the underlying causal genes is still labor-intensive and time-consuming. Here, we present a genetic tool for directly targeting the specific causal genes by using a single-gene resolution linkage map that was constructed from 3756 F2 rice plants via targeted sequencing technology and Tukey-Kramer multiple comparisons test. Three large- and moderate-effect QTLs, qHD6-2, qGL3-1, and qGW5-2, were successfully mapped to their specific causal genes, Hd1, GS3, and GW5, respectively. A complex genetic interaction network containing 30 QTL-QTL interactions was constructed, revealing that the alternative allele of hub QTL, qHD6-2, can hide or release the genetic contributions of the alleles at interacting loci. Moreover, arranging genetic interactions in the models lead to more accurate phenotypic predictions. These results provide a community resource and new feasible strategy for deciphering the genetic mechanisms of complex agronomic traits and accelerating crop breeding.

A new comprehensive annotation of leucine-rich repeat-containing receptors in rice

Oryza sativa (rice) plays an essential food security role for more than half of the world's population. Obtaining crops with high levels of disease resistance is a major challenge for breeders, especially today, given the urgent need for agriculture to be more sustainable. Plant resistance genes are mainly encoded by three large leucine-rich repeat (LRR)-containing receptor (LRR-CR) families: the LRR-receptor-like kinase (LRR-RLK), LRR-receptor-like protein (LRR-RLP) and nucleotide-binding LRR receptor (NLR). Using lrrprofiler, a pipeline that we developed to annotate and classify these proteins, we compared three publicly available annotations of the rice Nipponbare reference genome. The extended discrepancies that we observed for LRR-CR gene models led us to perform an in-depth manual curation of their annotations while paying special attention to nonsense mutations. We then transferred this manually curated annotation to Kitaake, a cultivar that is closely related to Nipponbare, using an optimized strategy. Here, we discuss the breakthrough achieved by manual curation when comparing genomes and, in addition to 'functional' and 'structural' annotations, we propose that the community adopts this approach, which we call 'comprehensive' annotation. The resulting data are crucial for further studies on the natural variability and evolution of LRR-CR genes in order to promote their use in breeding future resilient varieties.

Reciprocal adaptation of rice and Xanthomonas oryzae pv. oryzae: cross-species 2D GWAS reveals the underlying genetics

A 1D/2D genome-wide association study strategy was adopted to investigate the genetic systems underlying the reciprocal adaptation of rice (Oryza sativa) and its bacterial pathogen, Xanthomonas oryzae pv. oryzae (Xoo) using the whole-genome sequencing and large-scale phenotyping data of 701 rice accessions and 23 diverse Xoo strains. Forty-seven Xoo virulence-related genes and 318 rice quantitative resistance genes (QR-genes) mainly located in 41 genomic regions, and genome-wide interactions between the detected virulence-related genes and QR genes were identified, including well-known resistance genes/virulence genes plus many previously uncharacterized ones. The relationship between rice and Xoo was characterized by strong differentiation among Xoo races corresponding to the subspecific differentiation of rice, by strong shifts toward increased resistance/virulence of rice/Xoo populations and by rich genetic diversity at the detected rice QR-genes and Xoo virulence genes, and by genome-wide interactions between many rice QR-genes and Xoo virulence genes in a multiple-to-multiple manner, presumably resulting either from direct protein-protein interactions or from genetic epistasis. The observed complex genetic interaction system between rice and Xoo likely exists in other crop-pathogen systems that would maintain high levels of diversity at their QR-loci/virulence-loci, resulting in dynamic coevolutionary consequences during their reciprocal adaptation.

Endoplasmic reticulum stress pathway mediates the early heat stress response of developing rice seeds

A transient heat stress occurring during early seed development in rice (Oryza sativa) reduces seed size by altering endosperm development. However, the relationship between the timing of the stress and specific developmental stage on heat sensitivity is not well-understood. To address this, we imposed a series of non-overlapping heat stress treatments and found that young seeds are most sensitive during the first two days after flowering. Temporal transcriptome analysis of developing, heat stressed (35°C) seeds during this window shows that Inositol-requiring enzyme 1 (IRE1)-mediated endoplasmic reticulum (ER) stress response and jasmonic acid (JA) pathways are the early (1-3 h) drivers of heat stress response. We propose that increased JA levels under heat stress may precede ER stress response as JA application promotes the spliced form of OsbZIP50, an ER response marker gene linked to IRE1-specific pathway. This study presents temporal and mechanistic insights into the role of JA and ER stress signalling during early heat stress response of rice seeds that impact both grain size and quality. Modulating the heat sensitivity of the early sensing pathways and downstream endosperm development genes can enhance rice resilience to transient heat stress events.

The bundle sheath of rice is conditioned to play an active role in water transport as well as sulfur assimilation and jasmonic acid synthesis

Leaves comprise multiple cell types but our knowledge of the patterns of gene expression that underpin their functional specialization is fragmentary. Our understanding and ability to undertake the rational redesign of these cells is therefore limited. We aimed to identify genes associated with the incompletely understood bundle sheath of C₃ plants, which represents a key target associated with engineering traits such as C₄ photosynthesis into Oryza sativa (rice). To better understand the veins, bundle sheath and mesophyll cells of rice, we used laser capture microdissection followed by deep sequencing. Gene expression of the mesophyll is conditioned to allow coenzyme metabolism and redox homeostasis, as well as photosynthesis. In contrast, the bundle sheath is specialized in water transport, sulphur assimilation and jasmonic acid biosynthesis. Despite the small chloroplast compartment of bundle sheath cells, substantial photosynthesis gene expression was detected. These patterns of gene expression were not associated with the presence or absence of specific transcription factors in each cell type, but were instead associated with gradients in expression across the leaf. Comparative analysis with C₃ Arabidopsis identified a small gene set preferentially expressed in the bundle sheath cells of both species. This gene set included genes encoding transcription factors from 14 orthogroups and proteins allowing water transport, sulphate assimilation and jasmonic acid synthesis. The most parsimonious explanation for our findings is that bundle sheath cells from the last common ancestor of rice and Arabidopsis were specialized in this manner, and as the species diverged these patterns of gene expression have been maintained.

High-quality genome assembly of Huazhan and Tianfeng, the parents of an elite rice hybrid Tian-you-hua-zhan

High-quality rice reference genomes have accelerated the comprehensive identification of genome-wide variations and research on functional genomics and breeding. Tian-you-hua-zhan has been a leading hybrid in China over the past decade. Here, de novo genome assembly strategy optimization for the rice indica lines Huazhan (HZ) and Tianfeng (TF), including sequencing platforms, assembly pipelines and sequence depth, was carried out. The PacBio and Nanopore platforms for long-read sequencing were utilized, with the Canu, wtdbg2, SMARTdenovo, Flye, Canu-wtdbg2, Canu-SMARTdenovo and Canu-Flye assemblers. The combination of PacBio and Canu was optimal, considering the contig N50 length, contig number, assembled genome size and polishing process. The assembled contigs were scaffolded with Hi-C data, resulting in two "golden quality" rice reference genomes, and evaluated using the scaffold N50, BUSCO, and LTR assembly index. Furthermore, 42,625 and 41,815 non-transposable element genes were annotated for HZ and TF, respectively. Based on our assembly of HZ and TF, as well as Zhenshan97, Minghui63, Shuhui498 and 9311, comprehensive variations were identified using Nipponbare as a reference. The de novo assembly strategy for rice we optimized and the "golden quality" rice genomes we produced for HZ and TF will benefit rice genomics and breeding research, especially with respect to uncovering the genomic basis of the elite traits of HZ and TF.

Using Network-Based Machine Learning to Predict Transcription Factors Involved in Drought Resistance

Gene regulatory networks underpin stress response pathways in plants. However, parsing these networks to prioritize key genes underlying a particular trait is challenging. Here, we have built the Gene Regulation and Association Network (GRAiN) of rice (Oryza sativa). GRAiN is an interactive query-based web-platform that allows users to study functional relationships between transcription factors (TFs) and genetic modules underlying abiotic-stress responses. We built GRAiN by applying a combination of different network inference algorithms to publicly available gene expression data. We propose a supervised machine learning framework that complements GRAiN in prioritizing genes that regulate stress signal transduction and modulate gene expression under drought conditions. Our framework converts intricate network connectivity patterns of 2160 TFs into a single drought score. We observed that TFs with the highest drought scores define the functional, structural, and evolutionary characteristics of drought resistance in rice. Our approach accurately predicted the function of OsbHLH148 TF, which we validated using in vitro protein-DNA binding assays and mRNA sequencing loss-of-function mutants grown under control and drought stress conditions. Our network and the complementary machine learning strategy lends itself to predicting key regulatory genes underlying other agricultural traits and will assist in the genetic engineering of desirable rice varieties.

PyRice: a Python package for querying Oryza sativa databases

SUMMARY

Currently, gene information available for Oryza sativa species is located in various online heterogeneous data sources. Moreover, methods of access are also diverse, mostly web-based and sometimes query APIs, which might not always be straightforward for domain experts. The challenge is to collect information quickly from these applications and combine it logically, to facilitate scientific research. We developed a Python package named PyRice, a unified programing API to access all supported databases at the same time with consistent output. PyRice design is modular and implements a smart query system, which fits the computing resources to optimize the query speed. As a result, PyRice is easy to use and produces intuitive results.

AVAILABILITY AND IMPLEMENTATION

https://github.com/SouthGreenPlatform/PyRice.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Pan-genome analysis of 33 genetically diverse rice accessions reveals hidden genomic variations

Structural variations (SVs) and gene copy number variations (gCNVs) have contributed to crop evolution, domestication, and improvement. Here, we assembled 31 high-quality genomes of genetically diverse rice accessions. Coupling with two existing assemblies, we developed pan-genome-scale genomic resources including a graph-based genome, providing access to rice genomic variations. Specifically, we discovered 171,072 SVs and 25,549 gCNVs and used an Oryza glaberrima assembly to infer the derived states of SVs in the Oryza sativa population. Our analyses of SV formation mechanisms, impacts on gene expression, and distributions among subpopulations illustrate the utility of these resources for understanding how SVs and gCNVs shaped rice environmental adaptation and domestication. Our graph-based genome enabled genome-wide association study (GWAS)-based identification of phenotype-associated genetic variations undetectable when using only SNPs and a single reference assembly. Our work provides rich population-scale resources paired with easy-to-access tools to facilitate rice breeding as well as plant functional genomics and evolutionary biology research.

Genome-wide association mapping for resistance to bacterial blight and bacterial leaf streak in rice

Using genome-wide SNP association mapping, a total of 77 and 7 loci were identified for rice bacterial blight and bacterial leaf streak resistance, respectively, which may facilitate rice resistance improvement. Bacterial blight (BB) and bacterial leaf streak (BLS) caused by Gram-negative bacteria Xanthomonas oryzae pv. oryzae (Xoo) and X. oryzae pv. oryzicola (Xoc), respectively, are two economically important diseases negatively affecting rice production. To mine new sources of resistance, a set of rice germplasm collection consisting of 895 re-sequenced accessions from the 3000 Rice Genomes Project (3 K RGP) were screened for BB and BLS resistance under field conditions. Higher levels of BB resistance were observed in aus/boro subgroup, whereas the japonica, temperate japonica and tropical japonica subgroups possessed comparatively high levels of resistance to BLS. A genome-wide association study (GWAS) mined 77 genomic loci significantly associated with BB and 7 with BLS resistance. The phenotypic variance (R²) explained by these loci ranged from 0.4 to 30.2%. Among the loci, 7 for BB resistance were co-localized with known BB resistance genes and one for BLS resistance overlapped with a previously reported BLS resistance QTL. A search for the candidates in other novel loci revealed several defense-related genes that may be involved in resistance to BB and BLS. High levels of phenotypic resistance to BB or BLS could be attributed to the accumulation of the resistance (R) alleles at the associated loci, indicating their potential value in rice resistance breeding via gene pyramiding. The GWAS analysis validated the known genes underlying BB and BLS resistance and identified novel loci that could enrich the current resistance gene pool. The resources with strong resistance and significant SNPs identified in this study are potentially useful in breeding for BB and BLS resistance.

Meta-QTL and ortho-MQTL analyses identified genomic regions controlling rice yield, yield-related traits and root architecture under water deficit conditions

Meta-QTL (MQTL) analysis is a robust approach for genetic dissection of complex quantitative traits. Rice varieties adapted to non-flooded cultivation are highly desirable in breeding programs due to the water deficit global problem. In order to identify stable QTLs for major agronomic traits under water deficit conditions, we performed a comprehensive MQTL analysis on 563 QTLs from 67 rice populations published from 2001 to 2019. Yield and yield-related traits including grain weight, heading date, plant height, tiller number as well as root architecture-related traits including root dry weight, root length, root number, root thickness, the ratio of deep rooting and plant water content under water deficit condition were investigated. A total of 61 stable MQTLs over different genetic backgrounds and environments were identified. The average confidence interval of MQTLs was considerably refined compared to the initial QTLs, resulted in the identification of some well-known functionally characterized genes and several putative novel CGs for investigated traits. Ortho-MQTL mining based on genomic collinearity between rice and maize allowed identification of five ortho-MQTLs between these two cereals. The results can help breeders to improve yield under water deficit conditions.

Orphan Crops and their Wild Relatives in the Genomic Era

More than half of the calories consumed by humans are provided by three major cereal crops (rice, maize, and wheat). Orphan crops are usually well adapted to low-input agricultural conditions, and they not only play vital roles in local areas but can also contribute to food and nutritional needs worldwide. Interestingly, many wild relatives of orphan crops are important weeds of major crops. Although orphan crops and their wild relatives have received little attentions from researchers for many years, genomic studies have recently been performed on these plants. Here, we provide an overview of genomic studies on orphan crops, with a focus on orphan cereals and their wild relatives. The genomes of at least 12 orphan cereals and/or their wild relatives have been sequenced. In addition to genomic benefits for orphan crop breeding, we discuss the potential ways for mutual utilization of genomic data from major crops, orphan crops, and their wild relatives (including weeds) and provide perspectives on genetic improvement of both orphan and major crops (including de novo domestication of orphan crops) in the coming genomic era.

Gene Expression Dynamics in Rice Peduncles at the Heading Stage

Improving grain yield in the staple food crop rice has been long sought goal of plant biotechnology. One of the traits with significant impact on rice breeding programs is peduncle elongation at the time of heading failing which leads to significant reduction in grain yield due to incomplete panicle exsertion. To decipher transcriptional dynamics and molecular players underlying peduncle elongation, we performed RNA sequencing analysis of elongating and non-elongating peduncles in two Indian cultivars, Swarna and Pokkali, at the time of heading. Along with genes associated with cell division and cell wall biosynthesis, we observed significant enrichment of genes associated with auxins, gibberellins, and brassinosteroid biosynthesis/signaling in the elongating peduncles before heading in both the genotypes. Similarly, genes associated with carbohydrate metabolism and mobilization, abiotic stress response along with cytokinin, abscisic acid, jasmonic acid, and ethylene biosynthesis/signaling were enriched in non-elongating peduncles post heading. Significant enrichment of genes belonging to key transcription factor families highlights their specialized roles in peduncle elongation and grain filling before and after heading, respectively. A comparison with anther/pollen development-related genes provided 76 candidates with overlapping roles in anther/pollen development and peduncle elongation. Some of these are important for carbohydrate remobilization to the developing grains. These can be engineered to combat with incomplete panicle exsertion in male sterile lines and manipulate carbohydrate dynamics in grasses. Overall, this study provides baseline information about potential target genes for engineering peduncle elongation with implications on plant height, biomass composition and grain yields in rice.

Two nuclear effectors of the rice blast fungus modulate host immunity via transcriptional reprogramming

Pathogens utilize multiple types of effectors to modulate plant immunity. Although many apoplastic and cytoplasmic effectors have been reported, nuclear effectors have not been well characterized in fungal pathogens. Here, we characterize two nuclear effectors of the rice blast pathogen Magnaporthe oryzae. Both nuclear effectors are secreted via the biotrophic interfacial complex, translocated into the nuclei of initially penetrated and surrounding cells, and reprogram the expression of immunity-associated genes by binding on effector binding elements in rice. Their expression in transgenic rice causes ambivalent immunity: increased susceptibility to M. oryzae and Xanthomonas oryzae pv. oryzae, hemibiotrophic pathogens, but enhanced resistance to Cochliobolus miyabeanus, a necrotrophic pathogen. Our findings help remedy a significant knowledge deficiency in the mechanism of M. oryzae–rice interactions and underscore how effector-mediated manipulation of plant immunity by one pathogen may also affect the disease severity by other pathogens.

Plant pathogens secrete various effectors to manipulate host immunity. Here, Kim et al. describe two Magnaporthe oryzae effectors that translocate into the nuclei of infected rice cells and reprogram expression of immunity-associated genes, increasing susceptibility to hemibiotrophic pathogens.

Comparative population genomic analysis provides insights into breeding of modern indica rice in China

Comparative genomic analysis within Asian cultivated rice (Oryza sativa L.) populations has greatly enriched our knowledge regarding rice domestication and the divergence of the indica and japonica subspecies, while study on genomic regions associated with improvement within the indica subspecies is still limited. Here, through combined investigation of 2,429 indica cultivar genomes from public sequencing projects, we depict the improvement of modern indica rice in China. We identify three subgroups within indica populations: two geographically distinct, historical subgroups indica I (Ind_I) and indica III (Ind_III) and a modern subgroup indica II (Ind_II). The modern indica subgroup Ind_II shows admixture of the other two subgroups and enrichment of alleles that had been low-frequency in the other two subgroups. The Chinese indica cultivars exhibit a strong subgroup component change from Ind_I to Ind_II in the 1980s. Through haplotype-based comparative analysis, we detect 187 regions associated with separation of Ind_II compared to Ind_I or Ind_III. Within those regions we find strong representation of beneficial agricultural production-related alleles in Ind_II and a positive correlation between grain yield and number of differentiated haplotypes. Phenotypic features of long and slender grain, small tiller angle and decreased flowering time were detected for Ind_II. Through haplotype-based comparative analysis between rice subpopulations and subspecies, we find differentiated haplotypes not only from indica itself but also from japonica and aus, suggesting that introgression from other rice sub-populations has substantially contributed to modern indica rice breeds. These results help clarify the evolutionary landscape of modern indica rice in China and provide useful targets for future improvement.

Decoding the plant genome: From epigenome to 3D organization

The linear genome of eukaryotes is partitioned into diverse chromatin states and packaged into a three-dimensional (3D) structure, which has functional implications in DNA replication, DNA repair, and transcriptional regulation. Over the past decades, research on plant functional genomics and epigenomics has made great progress, with thousands of genes cloned and molecular mechanisms of diverse biological processes elucidated. Recently, 3D genome research has gradually attracted great attention of many plant researchers. Herein, we briefly review the progress in genomic and epigenomic research in plants, with a focus on Arabidopsis and rice, and summarize the currently used technologies and advances in plant 3D genome organization studies. We also discuss the relationships between one-dimensional linear genome sequences, epigenomic states, and the 3D chromatin architecture. This review provides basis for future research on plant 3D genomics.