A dynamic gene expression atlas covering the entire life cycle of rice

Mass spectrometry-based proteomic landscape of rice reveals a post-transcriptional regulatory role of N6-methyladenosine

Rice is one of the most important staple food and model species in plant biology, yet its quantitative proteomes are largely uncharacterized. Here we quantify the relative protein levels of over 15,000 genes across major rice tissues using a tandem mass tag strategy followed by intensive fractionation and mass spectrometry. We identify tissue-specific and tissue-enriched proteins that are linked to the functional specificity of individual tissues. Proteogenomic comparison of rice and Arabidopsis reveals conserved proteome expression, which differs from mammals in that there is a strong separation of species rather than tissues. Notably, profiling of N6-methyladenosine (m6A) across the rice major tissues shows that m6A at untranslated regions is negatively correlated with protein abundance and contributes to the discordance between RNA and protein levels. We also demonstrate that our data are valuable for identifying novel genes required for regulating m6A methylation. Taken together, this study provides a paradigm for further research into rice proteogenome.

Genome-wide characterization and analysis of rice DUF247 gene family

BACKGROUND

The domain of unknown function 247 (DUF247) proteins is involved in plant development and stress response. Rice is an important worldwide cereal crop, although an increasing number of DUF proteins have been identified, the understanding of DUF proteins is still very limited in rice.

RESULTS

In this study, we identified 69 genes that encode DUF247 proteins in the rice (Oryza sativa) genome by homology searches and domain prediction. All the OsDUF247 proteins were classified into four major groups (I, II, III and IV) by phylogenetic analysis. Remarkably, OsDUF247 genes clustered on the chromosomes solely show close phylogenetic relationships, suggesting that gene duplications have driven the expansion of the DUF247 gene family in the rice genome. Tissue profile analysis showed that most DUF247 genes expressed at constitutive levels in seedlings, roots, stems, and leaves, except for seven genes (LOC_Os01g21670, LOC_Os03g19700, LOC_Os05g04060, LOC_Os08g26820, LOC_Os08g26840, LOC_Os08g26850 and LOC_Os09g13410) in panicles. These seven genes were induced by various abiotic stress, including cold, drought, heat, hormone treatment, and especially salt, as demonstrated by further experimental analysis. DUF247 proteins contain transmembrane domains located on the membrane, suggesting their significant roles in rice development and adaptation to the environment.

CONCLUSIONS

These findings lay the foundation for functional characterizations of DUF247 genes to unravel their exact role in rice cultivars.

Transcriptome-wide association analyses reveal the impact of regulatory variants on rice panicle architecture and causal gene regulatory networks

Panicle architecture is a key determinant of rice grain yield and is mainly determined at the 1-2 mm young panicle stage. Here, we investigated the transcriptome of the 1-2 mm young panicles from 275 rice varieties and identified thousands of genes whose expression levels were associated with panicle traits. Multimodel association studies suggested that many small-effect genetic loci determine spikelet per panicle (SPP) by regulating the expression of genes associated with panicle traits. We found that alleles at cis-expression quantitative trait loci of SPP-associated genes underwent positive selection, with a strong preference for alleles increasing SPP. We further developed a method that integrates the associations of cis- and trans-expression components of genes with traits to identify causal genes at even small-effect loci and construct regulatory networks. We identified 36 putative causal genes of SPP, including SDT (MIR156j) and OsMADS17, and inferred that OsMADS17 regulates SDT expression, which was experimentally validated. Our study reveals the impact of regulatory variants on rice panicle architecture and provides new insights into the gene regulatory networks of panicle traits.

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.

Transcriptome and DNA methylome divergence of inflorescence development between two ecotypes in Panicum hallii

The morphological diversity of the inflorescence determines flower and seed production, which is critical for plant adaptation. Hall's panicgrass (Panicum hallii, P. hallii) is a wild perennial grass that has been developed as a model to study perennial grass biology and adaptive evolution. Highly divergent inflorescences have evolved between the two major ecotypes in P. hallii, the upland ecotype (P. hallii var hallii, HAL2 genotype) with compact inflorescence and large seed and the lowland ecotype (P. hallii var filipes, FIL2 genotype) with an open inflorescence and small seed. Here we conducted a comparative analysis of the transcriptome and DNA methylome, an epigenetic mark that influences gene expression regulation, across different stages of inflorescence development using genomic references for each ecotype. Global transcriptome analysis of differentially expressed genes (DEGs) and co-expression modules underlying the inflorescence divergence revealed the potential role of cytokinin signaling in heterochronic changes. Comparing DNA methylome profiles revealed a remarkable level of differential DNA methylation associated with the evolution of P. hallii inflorescence. We found that a large proportion of differentially methylated regions (DMRs) were located in the flanking regulatory regions of genes. Intriguingly, we observed a substantial bias of CHH hypermethylation in the promoters of FIL2 genes. The integration of DEGs, DMRs, and Ka/Ks ratio results characterized the evolutionary features of DMRs-associated DEGs that contribute to the divergence of the P. hallii inflorescence. This study provides insights into the transcriptome and epigenetic landscape of inflorescence divergence in P. hallii and a genomic resource for perennial grass biology.

Identification and co-expression network analysis of plumule-preferentially expressed genes in Oryza sativa

BACKGROUND

The seedling establishment is controlled by the programmed expression of sets of genes at the specific tissues of seed, abundance and environment. Plumule is an important part of the seed embryo and expresses the suits of genes to exert distinct functions during seed germination. Although rice genomic resources are available and developed rapidly, thousands of transcripts have not previously been located in the plumule of rice.

OBJECTIVE

This study was performed to identify plumule-preferentially expressed (OsPluP) genes in rice and determine the expression profiles and functions of OsPluP genes.

METHODS

We identified the OsPluP genes through Affymetrix microarray data. Meanwhile, qRT-PCR was performed to validate the expression pattern, also found that OsPluP genes were regulated by dark/light treatment. The cis-acting regulatory elements were analyzed in the promoters' regions of OsPluP genes. The T-DNA mutant of the OsPluP seed was used to reveal the function in seed germination.

RESULTS

In this study, a genomic survey of OsPluP genes was performed, and we identified 88 OsPluP genes based on Affymetrix microarray data. The expression profiles of 88 OsPluP members in 24 representative tissues covering rice whole life cycle can be roughly classified into three major groups, suggesting functional divergence of OsPluP genes in seed germination. The microarray data, qRT-PCR, and promoter analysis results demonstrated that transcripts of more than half OsPluPs (54 genes) could be enhanced in the darkness and respond to phytohormone. Gene Ontology (GO)and Kyoto encyclopedia of genes and genomes (KEGG) analysis demonstrated that OsPluP and their co-expressed genes were highly enriched in fatty acid metabolism. Moreover, OsPluP82 T-DNA mutant seeds displayed short plumule length and storage lipid accumulation.

CONCLUSION

This study would enable the functions of OsPluP genes during seed germination and contribute to the goal of molecular regulatory networks that lay the foundation for further studies of seedling growth.

Overexpression of maize GOLDEN2 in rice and maize calli improves regeneration by activating chloroplast development

Golden2 (G2), a member of the GARP transcription factor superfamily, regulates several biological processes and phytohormone signaling pathways in plants. In this study, we used a rice codon-optimized maize G2 gene (rZmG2) to improve the regeneration efficiency of rice and maize calli for genetic transformation. We isolated a promoter driving strong and callus-specific expression from rice to drive rZmG2 transcription from a transgene after transformation of two indica and two japonica rice cultivars. The resulting rZmG2 transgenic calli turned green in advance at the differentiation stage, thus significantly raising the regeneration rates of the transgenic indica and japonica rice plants relative to control transformations. Similar effect of this gene on improving maize transformation was also observed. Transcriptome sequencing and RT-qPCR analyses showed that many rice genes related to chloroplast development and phytohormones are upregulated in rZmG2-transgenic calli. These results demonstrate that rZmG2 can promote embryogenic callus differentiation and improve regeneration efficiency by activating chloroplast development and phytohormone pathways. We also established a heat-inducible Cre/loxP-based gene-excision system to remove rZmG2 and the antibiotic selectable gene after obtaining the transgenic plants. This study provides a useful tool for functional genomics work and biotechnology in plants.

Simultaneous alleviation of Cd availability in contaminated soil and accumulation in rice (Oryza sativa L.) by Fe-Mn oxide-modified biochar

Fe-Mn oxide-modified biochar (BC-FM) was used to remediate Cd-contaminated soil and mitigate Cd accumulation in rice. The roles of Fe and Mn in soil Cd immobilization and in controlling Cd uptake by rice were investigated via X-ray photoelectron spectroscopy (XPS) characterization and chemical analysis. Fe and Mn loaded on BC-FM increased the removal efficiencies of CaCl₂ extractable Cd in soil and Cd in pore water compared to those in only biochar (BC)-treated soil, with maximum removal rates at 67.9 % and 77.8 %, respectively. The XPS results indicated that the redox reactions of the Fe-Mn oxides on BC-FM surface affected Cd immobilization in the soil. The Fe (II/III) components on BC-FM were primarily converted to Fe₃O₄ in the soil system, which may form stable complexes with Cd²⁺ (Fe-O-Cd) during the entire rice growth period, and Cd may be bound to MnO or Mn₂O₃ in the form of CdMn₂O₄. The excellent adsorption performance of BC-FM enhanced by Fe-Mn oxides reduced the available Cd in the soil and stimulated Fe and Mn transport in rice, thereby inhibiting Cd accumulation in the aerial parts of rice. Cd concentrations in brown rice under BC-FM treatments reached the national safety standard (0.2 mg/kg, GB2762-2017). And BC-FM significantly increased the biomass of brown rice with a maximum rate of 26.8 %. These findings suggest that BC-FM could be used as an efficient material for Cd-contaminated soil remediation, and Fe-Mn plays important role in immobilizing Cd in soil and reducing Cd transport in rice.

Systematic annotation of conservation states provides insights into regulatory regions in rice

Plant genomes contain a large fraction of noncoding sequences. The discovery and annotation of conserved noncoding sequences (CNSs) in plants is an ongoing challenge. Here we report the application of comparative genomics to systematically identify CNSs in 50 well-annotated Gramineae genomes using rice (Oryza sativa) as the reference. We conduct multiple-way whole-genome alignments to the rice genome. The rice genome is annotated as 20 conservation states (CSs) at single-nucleotide resolution using a multivariate hidden Markov model (ConsHMM) based on the multiple-genome alignments. Different states show distinct enrichments for various genomic features, and the conservation scores of CSs are highly correlated with the level of associated chromatin accessibility. We find that at least 33.5% of the rice genome is highly under selection, with more than 70% of the sequence lying outside of coding regions. A catalog of 855,366 regulatory CNSs is generated, and they significantly overlapped with putative active regulatory elements such as promoters, enhancers, and transcription factor binding sites. Collectively, our study provides a resource for elucidating functional noncoding regions of the rice genome and an evolutionary aspect of regulatory sequences in higher plants.

Improvement of Rice Agronomic Traits by Editing Type-B Response Regulators

Type-B response regulator proteins in rice contain a conserved receiver domain, followed by a GARP DNA binding domain and a longer C-terminus. Some type-B response regulators such as RR21, RR22 and RR23 are involved in the development of rice leaf, root, flower and trichome. In this study, to evaluate the application potential of type-B response regulators in rice genetic improvement, thirteen type-B response regulator genes in rice were respectively knocked out by using CRISPR/Cas9 genome editing technology. Two guide RNAs (gRNAs) were simultaneously expressed on a knockout vector to mutate one gene. T₀ transformed plants were used to screen the plants with deletion of large DNA fragments through PCR with specific primers. The mutants of CRISPR/Cas9 gene editing were detected by Cas9 specific primer in the T₁ generation, and homozygous mutants without Cas9 were screened, whose target regions were confirmed by sequencing. Mutant materials of 12 OsRRs were obtained, except for RR24. Preliminary phenotypic observation revealed variations of various important traits in different mutant materials, including plant height, tiller number, tillering angle, heading date, panicle length and yield. The osrr30 mutant in the T₂ generation was then further examined. As a result, the heading date of the osrr30 mutant was delayed by about 18 d, while the yield was increased by about 30%, and the chalkiness was significantly reduced compared with those of the wild-type under field high temperature stress. These results indicated that osrr30 has great application value in rice breeding. Our findings suggest that it is feasible to perform genetic improvement of rice by editing the type-B response regulators.

4D genetic networks reveal the genetic basis of metabolites and seed oil-related traits in 398 soybean RILs

Background

The yield and quality of soybean oil are determined by seed oil-related traits, and metabolites/lipids act as bridges between genes and traits. Although there are many studies on the mode of inheritance of metabolites or traits, studies on multi-dimensional genetic network (MDGN) are limited.

Results

In this study, six seed oil-related traits, 59 metabolites, and 107 lipids in 398 recombinant inbred lines, along with their candidate genes and miRNAs, were used to construct an MDGN in soybean. Around 175 quantitative trait loci (QTLs), 36 QTL-by-environment interactions, and 302 metabolic QTL clusters, 70 and 181 candidate genes, including 46 and 70 known homologs, were previously reported to be associated with the traits and metabolites, respectively. Gene regulatory networks were constructed using co-expression, protein–protein interaction, and transcription factor binding site and miRNA target predictions between candidate genes and 26 key miRNAs. Using modern statistical methods, 463 metabolite–lipid, 62 trait–metabolite, and 89 trait–lipid associations were found to be significant. Integrating these associations into the above networks, an MDGN was constructed, and 128 sub-networks were extracted. Among these sub-networks, the gene–trait or gene–metabolite relationships in 38 sub-networks were in agreement with previous studies, e.g., oleic acid (trait)–GmSEI–GmDGAT1a–triacylglycerol (16:0/18:2/18:3), gene and metabolite in each of 64 sub-networks were predicted to be in the same pathway, e.g., oleic acid (trait)–GmPHS–d-glucose, and others were new, e.g., triacylglycerol (16:0/18:1/18:2)–GmbZIP123–GmHD-ZIPIII-10–miR166s–oil content.

Conclusions

This study showed the advantages of MGDN in dissecting the genetic relationships between complex traits and metabolites. Using sub-networks in MGDN, 3D genetic sub-networks including pyruvate/threonine/citric acid revealed genetic relationships between carbohydrates, oil, and protein content, and 4D genetic sub-networks including PLDs revealed the relationships between oil-related traits and phospholipid metabolism likely influenced by the environment. This study will be helpful in soybean quality improvement and molecular biological research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13068-022-02191-1.

In-Silico Study of Brassinosteroid Signaling Genes in Rice Provides Insight Into Mechanisms Which Regulate Their Expression

Brassinosteroids (BRs) regulate a diverse spectrum of processes during plant growth and development and modulate plant physiology in response to environmental fluctuations and stress factors. Thus, the BR signaling regulators have the potential to be targeted for gene editing to optimize the architecture of plants and make them more resilient to environmental stress. Our understanding of the BR signaling mechanism in monocot crop species is limited compared to our knowledge of this process accumulated in the model dicot species - Arabidopsis thaliana. A deeper understanding of the BR signaling and response during plant growth and adaptation to continually changing environmental conditions will provide insight into mechanisms that govern the coordinated expression of the BR signaling genes in rice (Oryza sativa) which is a model for cereal crops. Therefore, in this study a comprehensive and detailed in silico analysis of promoter sequences of rice BR signaling genes was performed. Moreover, expression profiles of these genes during various developmental stages and reactions to several stress conditions were analyzed. Additionally, a model of interactions between the encoded proteins was also established. The obtained results revealed that promoters of the 39 BR signaling genes are involved in various regulatory mechanisms and interdependent processes that influence growth, development, and stress response in rice. Different transcription factor-binding sites and cis-regulatory elements in the gene promoters were identified which are involved in regulation of the genes’ expression during plant development and reactions to stress conditions. The in-silico analysis of BR signaling genes in O. sativa provides information about mechanisms which regulate the coordinated expression of these genes during rice development and in response to other phytohormones and environmental factors. Since rice is both an important crop and the model species for other cereals, this information may be important for understanding the regulatory mechanisms that modulate the BR signaling in monocot species. It can also provide new ways for the plant genetic engineering technology by providing novel potential targets, either cis-elements or transcriptional factors, to create elite genotypes with desirable traits.

Characterization and expression analysis of the glycosyltransferase 64 familyin rice (Oryza sativa)

The glycosyltransferase 64 (GT64) family is widely conserved in many species, including animals and plants. The functions of GT64 family genes in animals have been well characterized in the biosynthesis of extracellular heparan sulfate, whereas two GT64 members in Arabidopsis thaliana are involved in the glycosylation of plasma membrane glycosylinositol phosphorylceramides (GIPCs). GIPCs are the main components of plant sphingolipids and serve as important signal molecules in various developmental processes and stress responses. Rice (Oryza sativa), a model monocot plant, contains four GT64 members in its genome. Using phylogenetic analysis, 73 GT64s from 19 plant species were divided into three main groups. Each group can be represented by the three members in Arabidopsis and show a trend of monocot-eudicot divergences. A promoter and genomic variation analysis of GT64s in rice showed that various stress-related regulatory elements exist in their promoters, and many sequence variations were found between the two main rice subspecies, japonica and indica. Additionally, the transmembrane domain and subcellular localization analyses revealed that these genes all encode membrane-bound glycosyltransferases and localize to the Golgi apparatus. Finally, expression analysis of the four GT64 genes in rice, as assessed by real-time quantitative PCR, showed that they have distinct tissue-specific expression patterns and respond to different hormone treatments or abiotic stresses. Our results indicated that this family of genes may play a role in different stress responses and hormone signaling pathways in rice, and therefore provides fundamental information for the further investigation of their function in rice.

Systems seed biology to understand and manipulate rice grain quality and nutrition

Rice is one of the most essential crops since it meets the calorific needs of 3 billion people around the world. Rice seed development initiates upon fertilization, leading to the establishment of two distinct filial tissues, the endosperm and embryo, which accumulate distinct seed storage products, such as starch, storage proteins, and lipids. A range of systems biology tools deployed in dissecting the spatiotemporal dynamics of transcriptome data, methylation, and small RNA based regulation operative during seed development, influencing the accumulation of storage products was reviewed. Studies of other model systems are also considered due to the limited information on the rice transcriptome. This review highlights key genes identified through a holistic view of systems biology targeted to modify biochemical composition and influence rice grain quality and nutritional value with the target of improving rice as a functional food.

Gene fusion as an important mechanism to generate new genes in the genus Oryza

BACKGROUND

Events of gene fusion have been reported in several organisms. However, the general role of gene fusion as part of new gene origination remains unknown.

RESULTS

We conduct genome-wide interrogations of four Oryza genomes by designing and implementing novel pipelines to detect fusion genes. Based on the phylogeny of ten plant species, we detect 310 fusion genes across four Oryza species. The estimated rate of origination of fusion genes in the Oryza genus is as high as 63 fusion genes per species per million years, which is fixed at 16 fusion genes per species per million years and much higher than that in flies. By RNA sequencing analysis, we find more than 44% of the fusion genes are expressed and 90% of gene pairs show strong signals of purifying selection. Further analysis of CRISPR/Cas9 knockout lines indicates that newly formed fusion genes regulate phenotype traits including seed germination, shoot length and root length, suggesting the functional significance of these genes.

CONCLUSIONS

We detect new fusion genes that may drive phenotype evolution in Oryza. This study provides novel insights into the genome evolution of Oryza.

Twenty years of rice genomics research: From sequencing and functional genomics to quantitative genomics

Since the completion of the rice genome sequencing project in 2005, we have entered the era of rice genomics, which is still in its ascendancy. Rice genomics studies can be classified into three stages: structural genomics, functional genomics, and quantitative genomics. Structural genomics refers primarily to genome sequencing for the construction of a complete map of rice genome sequence. This is fundamental for rice genetics and molecular biology research. Functional genomics aims to decode the functions of rice genes. Quantitative genomics is large-scale sequence- and statistics-based research to define the quantitative traits and genetic features of rice populations. Rice genomics has been a transformative influence on rice biological research and contributes significantly to rice breeding, making rice a good model plant for studying crop sciences.

The rice heavy-metal transporter OsNRAMP1 regulates disease resistance by modulating ROS homoeostasis

Crop diseases threaten food security and sustainable agriculture. Consumption of crops containing nonessential toxic metals leads to health risks for humans. Therefore, cultivation of disease-resistant and toxic metal-safe crops is a double-gain strategy that can contribute to food security. Here, we show that rice heavy-metal transporter OsNRAMP1 plays an important role in plant immunity by modulating metal ion and reactive oxygen species (ROS) homoeostasis. OsNRAMP1 expression was induced after pathogenic bacteria and fungi infections. The osnramp1 mutants had an increased content of H₂ O₂ and activity of superoxide dismutase, but decreased activity of catalase, showing enhanced broad-spectrum resistance against bacterial and fungal pathogens. RNA-seq analysis identified a number of differentially expressed genes that were involved in metal ion and ROS homoeostasis. Altered expression of metal ion-dependent ROS-scavenging enzymes genes and lower accumulation of cations such as Mn together induced compromised metal ion-dependent enzyme-catalysing activity and modulated ROS homoeostasis, which together contributed towards disease resistance in osnramp1 mutants. Furthermore, the osnramp1 mutants contained lower levels of toxic heavy metals Cd and Pb and micronutrients Ni and Mn in leaves and grains. Taken together, a proof of concept was achieved that broad-spectrum disease-resistant and toxic heavy-metal-safe rice was engineered by removal of the OsNRAMP1 gene.

Understanding the genetic and molecular constitutions of heterosis for developing hybrid rice

The wide adoption of hybrid rice has greatly increased rice yield in the last several decades. The utilization of heterosis facilitated by male sterility has been a common strategy for hybrid rice development. Here, we summarize our efforts in the genetic and molecular understanding of heterosis and male sterility together with the related progress from other research groups. Analyses of F₁ diallel crosses show that strong heterosis widely exists in hybrids of diverse germplasms, and inter-subspecific hybrids often display higher heterosis. Using the elite hybrid Shanyou 63 as a model, an immortalized F₂ population design is conducted for systematic characterization of the biological mechanism of heterosis, with identification of loci controlling heterosis of yield and yield component traits. Dominance, overdominance, and epistasis all play important roles in the genetic basis of heterosis; quantitative assessment of these components well addressed the three classical genetic hypotheses for heterosis. Environment-sensitive genic male sterility (EGMS) enables the development of two-line hybrids, and long noncoding RNAs often function as regulators of EGMS. Inter-subspecific hybrids show greatly reduced fertility; the identification and molecular characterization of hybrid sterility genes offer strategies for overcoming inter-subspecific hybrid sterility. These developments have significant implications for future hybrid rice improvement using genomic breeding.

Isolation and functional analysis of OsAOS1 promoter for resistance to Nilaparvata lugens Stål infestation in rice

Insect pests have a great impact on the yield and quality of crops. Insecticide applications are an effective method of pest control, however, they also have adverse effects on the environment. Using insect-inducible promoters to drive insect-resistant genes in transgenic crops is a potential sustainable pest management strategy, but insect-inducible promoters have been rarely reported. In this study, we found rice allene oxide synthase gene (AOS, LOC_Os03g12500) can be highly upregulated following brown planthopper (Nilaparvata lugens Stål, BPH) infestation. Then, we amplified the promoter of OsAOS1 and the β- glucuronidase reporter gene was used to analyze the expression pattern of the promoter. Through a series of 5' truncated assays, three positive regulatory regions in response to BPH infestation in the promoter were identified. The transgenic plants, P1R123-min 35S and P1TR1-min 35S promoter-driven snowdrop lectin (Galanthus nivalis agglutinin, GNA) gene, demonstrated the highest expression levels of GNA and lowest BPH survival. Our work identified a BPH-inducible promoter and three positive regions within it. Transgenic rice with GNA driven by OsAOS1 promoter and positive regions exhibited an expected lethal effect on BPH. This study proved the application potential of BPH-inducible promoter and provided a novel path for the selection of insect-resistant tools in the future.

Epigenetic regulation of seed-specific gene expression by DNA methylation valleys in castor bean

BACKGROUND

Understanding the processes governing angiosperm seed growth and development is essential both for fundamental plant biology and for agronomic purposes. Master regulators of angiosperm seed development are expressed in a seed-specific manner. However, it is unclear how this seed specificity of transcription is established. In some vertebrates, DNA methylation valleys (DMVs) are highly conserved and strongly associated with key developmental genes, but comparable studies in plants are limited to Arabidopsis and soybean. Castor bean (Ricinus communis) is a valuable model system for the study of seed biology in dicots and source of economically important castor oil. Unlike other dicots such as Arabidopsis and soybean, castor bean seeds have a relatively large and persistent endosperm throughout seed development, representing substantial structural differences in mature seeds. Here, we performed an integrated analysis of RNA-seq, whole-genome bisulfite sequencing, and ChIP-seq for various histone marks in the castor bean.

RESULTS

We present a gene expression atlas covering 16 representative tissues and identified 1162 seed-specific genes in castor bean (Ricinus communis), a valuable model for the study of seed biology in dicots. Upon whole-genome DNA methylation analyses, we detected 32,567 DMVs across five tissues, covering ~33% of the castor bean genome. These DMVs are highly hypomethylated during development and conserved across plant species. We found that DMVs have the potential to activate transcription, especially that of tissue-specific genes. Focusing on seed development, we found that many key developmental regulators of seed/endosperm development, including AGL61, AGL62, LEC1, LEC2, ABI3, and WRI1, were located within DMVs. ChIP-seq for five histone modifications in leaves and seeds clearly showed that the vast majority of histone modification peaks were enriched within DMVs, and their remodeling within DMVs has a critical role in the regulation of seed-specific gene expression. Importantly, further experiment analysis revealed that distal DMVs may act as cis-regulatory elements, like enhancers, to activate downstream gene expression.

CONCLUSIONS

Our results point to the importance of DMVs and special distal DMVs behaving like enhancers, in the regulation of seed-specific genes, via the reprogramming of histone modifications within DMVs. Furthermore, these results provide a comprehensive understanding of the epigenetic regulator roles in seed development in castor bean and other important crops.

Global Identification and Characterization of C2 Domain-Containing Proteins Associated with Abiotic Stress Response in Rice (Oryza sativa L.)

C2 domain-containing proteins (C2DPs) have been identified in different genomes that contain single or multiple C2 domains in their C- or N-terminal. It possesses higher functional activity in the transmembrane regions. The identification of C2 domains were reported in a previous study, such as multiple C2 domains and transmembrane-region proteins (MCTPs) and N-terminal-TM-C2 domain proteins (NTMC2s) of rice, Arabidopsis thaliana, and cotton, whereas the C2DP gene family in rice has not been comprehensively studied, and the role of the C2DP gene in rice in response to abiotic stress is not yet fully understood. In this study, we identified 82 C2DPs in the rice genome and divided them into seven groups through phylogenetic analysis. The synteny analysis revealed that duplication events were either exhibited within the genome of rice or between the genomes of rice and other species. Through the analysis of cis-acting elements in promoters, expression profiles, and qRT-PCR results, the functions of OsC2DPs were found to be widely distributed in diverse tissues and were extensively involved in phytohormones-related and abiotic stresses response in rice. The prediction of the microRNA (miRNA) targets of OsC2DPs revealed the possibility of regulation by consistent miRNAs. Notably, OsC2DP50/51/52 as a co-tandem duplication exhibited similar expression variations and involved the coincident miRNA-regulation pathway. Moreover, the results of the genotypic variation and haplotype analysis revealed that OsC2DP17, OsC2DP29, and OsC2DP49 were associated with cold stress responses. These findings provided comprehensive insights for characterizations of OsC2DPs in rice as well as for their roles for abiotic stress.

Rice metabolic regulatory network spanning the entire life cycle

As one of the most important crops in the world, rice (Oryza sativa) is a model plant for metabolome research. Although many studies have focused on the analysis of specific tissues, the changes in metabolite abundance across the entire life cycle have not yet been determined. In this study, combining both targeted and nontargeted metabolite profiling methods, a total of 825 annotated metabolites were quantified in rice samples from different tissues covering the entire life cycle. The contents of metabolites in different tissues of rice were significantly different, with various metabolites accumulating in the plumule and radicle during seed germination. Combining these data with transcriptome data obtained from the same time period, we constructed the Rice Metabolic Regulation Network. The metabolites and co-expressed genes were further divided into 12 clusters according to their accumulation patterns, with members within each cluster displaying a uniform and clear pattern of abundance across development. Using this dataset, we established a comprehensive metabolic profile of the rice life cycle and used two independent strategies to identify novel transcription factors-namely the use of known regulatory genes as bait to screen for new networks underlying lignin metabolism and the unbiased identification of new glycerophospholipid metabolism regulators on the basis of tissue specificity. This study thus demonstrates how guilt-by-association analysis of metabolome and transcriptome data spanning the entire life cycle in cereal crops provides novel resources and tools to aid in understanding the mechanisms underlying important agronomic traits.

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 T0 plants and 83.2% of T1 plants. A high proportion (55.6%) of homozygous marker-excised plants were obtained in T1 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.

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 atlas of Phalaenopsis equestris

The vast diversity of Orchidaceae together with sophisticated adaptations to pollinators and other unique features make this family an attractive model for evolutionary and functional studies. The sequenced genome of Phalaenopsis equestris facilitates Orchidaceae research. Here, we present an RNA-seq-based transcriptome map of P. equestris that covers 19 organs of the plant, including leaves, roots, floral organs and the shoot apical meristem. We demonstrated the high quality of the data and showed the similarity of the P. equestris transcriptome map with the gene expression atlases of other plants. The transcriptome map can be easily accessed through our database Transcriptome Variation Analysis (TraVA) for visualizing gene expression profiles. As an example of the application, we analyzed the expression of Phalaenopsis “orphan” genes–those that do not have recognizable similarity with the genes of other plants. We found that approximately half of these genes were not expressed; the ones that were expressed were predominantly expressed in reproductive structures.

Toward Integrated Multi-Omics Intervention: Rice Trait Improvement and Stress Management

Rice (Oryza sativa) is an imperative staple crop for nearly half of the world's population. Challenging environmental conditions encompassing abiotic and biotic stresses negatively impact the quality and yield of rice. To assure food supply for the unprecedented ever-growing world population, the improvement of rice as a crop is of utmost importance. In this era, "omics" techniques have been comprehensively utilized to decipher the regulatory mechanisms and cellular intricacies in rice. Advancements in omics technologies have provided a strong platform for the reliable exploration of genetic resources involved in rice trait development. Omics disciplines like genomics, transcriptomics, proteomics, and metabolomics have significantly contributed toward the achievement of desired improvements in rice under optimal and stressful environments. The present review recapitulates the basic and applied multi-omics technologies in providing new orchestration toward the improvement of rice desirable traits. The article also provides a catalog of current scenario of omics applications in comprehending this imperative crop in relation to yield enhancement and various environmental stresses. Further, the appropriate databases in the field of data science to analyze big data, and retrieve relevant information vis-à-vis rice trait improvement and stress management are described.

The Divergent Roles of the Rice bcl-2 Associated Athanogene (BAG) Genes in Plant Development and Environmental Responses

Bcl-2-associated athanogene (BAG), a group of proteins evolutionarily conserved and functioned as co-chaperones in plants and animals, is involved in various cell activities and diverse physiological processes. However, the biological functions of this gene family in rice are largely unknown. In this study, we identified a total of six BAG members in rice. These genes were classified into two groups, OsBAG1, -2, -3, and -4 are in group I with a conserved ubiquitin-like structure and OsBAG5 and -6 are in group Ⅱ with a calmodulin-binding domain, in addition to a common BAG domain. The BAG genes exhibited diverse expression patterns, with OsBAG4 showing the highest expression level, followed by OsBAG1 and OsBAG3, and OsBAG6 preferentially expressed in the panicle, endosperm, and calli. The co-expression analysis and the hierarchical cluster analysis indicated that the OsBAG1 and OsBAG3 were co-expressed with primary cell wall-biosynthesizing genes, OsBAG4 was co-expressed with phytohormone and transcriptional factors, and OsBAG6 was co-expressed with disease and shock-associated genes. β-glucuronidase (GUS) staining further indicated that OsBAG3 is mainly involved in primary young tissues under both primary and secondary growth. In addition, the expression of the BAG genes under brown planthopper (BPH) feeding, N, P, and K deficiency, heat, drought and plant hormones treatments was investigated. Our results clearly showed that OsBAGs are multifunctional molecules as inferred by their protein structures, subcellular localizations, and expression profiles. BAGs in group I are mainly involved in plant development, whereas BAGs in group II are reactive in gene regulations and stress responses. Our results provide a solid basis for the further elucidation of the biological functions of plant BAG genes.

A genome-wide association study reveals that the 2-oxoglutarate/malate translocator mediates seed vigor in rice

Seed vigor is an important trait for the direct seeding of rice (Oryza sativa L.). In this study, we examined the genetic architecture of variation in the germination rate using a diverse panel of rice accessions. Four quantitative trait loci for germination rate were identified using a genome-wide association study during early germination. One candidate gene, encoding the 2-oxoglutarate/malate translocator (OsOMT), was validated for qGR11. Disruption of this gene (Osomt mutants) reduced seed vigor, including seed germination and seedling growth, in rice. Functional analysis revealed that OsOMT influences seed vigor mainly by modulating amino acid levels and glycolysis and tricarboxylic acid cycle processes. The levels of most amino acids, including the Glu family (Glu, Pro, Arg, and GABA), Asp family (Asp, Thr, Lys, Ile, and Met), Ser family (Ser, Gly, and Cys), and others (His, Ala, Leu, and Val), were significantly reduced in the mature grains and the early germinating seeds of Osomt mutants compared to wild type (WT). The glucose and soluble sugar contents, as well as adenosine triphosphate levels, were significantly decreased in germinating seeds of Osomt mutants compared to WT. These results provide important insights into the role of OsOMT in seed vigor in rice.

An inferred functional impact map of genetic variants in rice

Interpreting the functional impacts of genetic variants (GVs) is an important challenge for functional genomic studies in crops and next-generation breeding. Previous studies in rice (Oryza sativa) have focused mainly on the identification of GVs, whereas systematic functional annotation of GVs has not yet been performed. Here, we present a functional impact map of GVs in rice. We curated haplotype information for 17 397 026 GVs from sequencing data of 4726 rice accessions. We quantitatively evaluated the effects of missense mutations in coding regions in each haplotype based on the conservation of amino acid residues and obtained the effects of 918 848 non-redundant missense GVs. Furthermore, we generated high-quality chromatin accessibility (CA) data from six representative rice tissues and used these data to train deep convolutional neural network models to predict the impacts of 5 067 405 GVs for CA in regulatory regions. We characterized the functional properties and tissue specificity of the GV effects and found that large-effect GVs in coding and regulatory regions may be subject to selection in different directions. Finally, we demonstrated how the functional impact map could be used to prioritize causal variants in mapping populations. This impact map will be a useful resource for accelerating gene cloning and functional studies in rice, and can be freely queried in RiceVarMap V2.0 (http://ricevarmap.ncpgr.cn).

Genome-Wide Association Study in Rice Revealed a Novel Gene in Determining Plant Height and Stem Development, by Encoding a WRKY Transcription Factor

Semi-dwarfism is a main agronomic trait in crop breeding. In this study, we performed genome-wide association study (GWAS) and identified a new quantitative trait nucleotide (QTN) for rice shoot length. The peak QTN (C/T) was located in the first coding region of a group III WRKY transcription factor OsWRKY21 (LOC_Os01g60640). Interestingly, further haplotype analysis showed that C/T difference only existed in the indica group but not in the japonica group, resulting in significant differences in plant height among the different indica rice varieties. OsWRKY21 was expressed in embryo, radicle, shoots, leaves, and stems. Most notably, overexpressing OsWRKY21 resulted in the semi-dwarf phenotype, early heading date and short internodes compared to the wild type, while the knockout mutant plants by CRISPR/Cas9 technology yielded the opposite. The overexpressing lines exhibited the decreased length of the cells near sclerenchyma epidermis, accompanied with the lower levels of indole-3-acetic acid (IAA) and gibberellin 3 (GA₃), but increased levels of the abscisic acid (ABA) and salicylic acid (SA) in the internodes at heading stage. Moreover, the semi-dwarf phenotype could be fully rescued by exogenous GA₃ application at seedling stage. The RNA-seq and qRT-PCR analysis confirmed the differential expression levels of genes in development and the stress responses in rice, including GA metabolism (GA20ox2, GA2ox6, and YABY1) and cell wall biosynthesis (CesA4, 7, and 9) and regulation (MYB103L). These data suggest the essential role of OsWRKY21 in regulation of internode elongation and plant height in rice.

Homoploid F1 hybrids and segmental allotetraploids of japonica and indica rice subspecies show similar and enhanced tolerance to nitrogen deficiency than parental lines

It remains unclear whether the merger of two divergent genomes by hybridization at the homoploid level or coupled with whole-genome duplication (WGD; allopolyploidy) can result in plants having better tolerance to stress conditions. In this study, we compared phenotypic performance and gene expression in the two diploid subspecies of rice (Oryza sativa subsp. japonica and indica), their reciprocal F1 hybrids, and in segmental allotetraploids under normal and nitrogen (N)-deficient conditions. We found that F1 hybrids and tetraploids showed higher and similar levels of tolerance to N deficiency than either parent. In parallel, total expression levels of 18 relevant functional genes were less perturbed by N deficiency in the F1 hybrids and tetraploids than in the parents. This was consistent with stable intrinsic partitioning of allelic/homoeologous expression defined by parental legacy in the homoploid F1 hybrids/tetraploids between the two conditions. The results suggest that genetic additivity at both the homoploid and allopolyploidy level might lead to similar beneficial phenotypic responses to nitrogen stress compared with the parents. The lack of synergistic responses to N limitation concomitant with WGD, relative to that exhibited by F1 hybrids, adds new empirical evidence in support of the emerging hypothesis that hybridization by itself can play a significant role in plant adaptive evolution in times of stress.

Bract suppression regulated by the miR156/529-SPLs-NL1-PLA1 module is required for the transition from vegetative to reproductive branching in rice

Reproductive transition of grasses is characterized by switching the pattern of lateral branches, featuring the suppression of outgrowth of the subtending leaves (bracts) and rapid formation of higher-order branches in the inflorescence (panicle). However, the molecular mechanisms underlying such changes remain largely unknown. Here, we show that bract suppression is required for the reproductive branching in rice. We identified a pathway involving the intrinsic time ruler microRNA156/529, their targets SQUAMOSA PROMOTER BINDING PROTEIN LIKE (SPL) genes, NECK LEAF1 (NL1), and PLASTOCHRON1 (PLA1), which regulates the bract outgrowth and thus affects the pattern switch between vegetative and reproductive branching. Suppression of the bract results in global reprogramming of transcriptome and chromatin accessibility following the reproductive transition, while these processes are largely dysregulated in the mutants of these genes. These discoveries contribute to our understanding of the dynamic plant architecture and provide novel insights for improving crop yields.

Systems biology for crop improvement

In recent years, generation of large-scale data from genome, transcriptome, proteome, metabolome, epigenome, and others, has become routine in several plant species. Most of these datasets in different crop species, however, were studied independently and as a result, full insight could not be gained on the molecular basis of complex traits and biological networks. A systems biology approach involving integration of multiple omics data, modeling, and prediction of the cellular functions is required to understand the flow of biological information that underlies complex traits. In this context, systems biology with multiomics data integration is crucial and allows a holistic understanding of the dynamic system with the different levels of biological organization interacting with external environment for a phenotypic expression. Here, we present recent progress made in the area of various omics studies-integrative and systems biology approaches with a special focus on application to crop improvement. We have also discussed the challenges and opportunities in multiomics data integration, modeling, and understanding of the biology of complex traits underpinning yield and stress tolerance in major cereals and legumes.

Identification and candidate gene screening of qCIR9.1, a novel QTL associated with anther culturability in rice (Oryza sativa L.)

A novel QTL, qCIR9.1, that controls callus induction rate in anther culture was identified on chromosome 9 in rice, and based on RNA-seq data, Os09g0551600 was the most promising candidate gene. Anther culture, a doubled haploid (DH) technique, has become an important technology in many plant-breeding programmes. Although anther culturability is the key factor in this technique, its genetic mechanisms in rice remain poorly understood. In this study, we mapped quantitative trait loci (QTLs) responsible for anther culturability by using 192 recombinant inbred lines (RILs) derived from YZX (Oryza sativa ssp. indica) × 02428 (Oryza sativa ssp. japonica) and a high-density bin map. A total of eight QTLs for anther culturability were detected in three environments. Among these QTLs, a novel major QTL for callus induction rate (CIR) named qCIR9.1 was repeatedly mapped to a ~ 100 kb genomic interval on chromosome 9 and explained 8.39-14.14% of the phenotypic variation. Additionally, RNA sequencing (RNA-seq) was performed for the parents (YZX and 02428), low- (L-Pool) and high-CIR RILs (H-Pool) after 16 and 26 days of culture. By using the RNA of the bulked RILs for background normalization, the number of differentially expressed genes (DEGs) both between the parents and between the bulked RILs after 26 days of culture was drastically reduced to only 78. Among these DEGs, only one gene, Os09g0551600, encoding a high-mobility group (HMG) protein, was located in the candidate region of qCIR9.1. qRT-PCR analysis of Os09g0551600 showed the same results as RNA-seq, and the expression of this gene was decreased in the low-callus-induction parent (YZX) and L-Pool. Our results provide a foundational step for further cloning of qCIR9.1 and will be very useful for improving anther culturability in rice.

OsTPR boosts the superior grains through increase in upper secondary rachis branches without incurring a grain quality penalty

To address the future food security in Asia, we need to improve the genetic gain of grain yield while ensuring the consumer acceptance. This study aimed to identify novel genes influencing the number of upper secondary rachis branches (USRB) to elevate superior grains without compromising grain quality by studying the genetic variance of 310 diverse O. sativa var. indica panel using single- and multi-locus genome-wide association studies (GWAS), gene set analyses and gene regulatory network analysis. GWAS of USRB identified 230 significant (q-value < 0.05) SNPs from chromosomes 1 and 2. GWAS targets narrowed down using gene set analyses identified large effect association on an important locus LOC_Os02g50790/LOC_Os02g50799 encoding a nuclear-pore anchor protein (OsTPR). The superior haplotype derived from non-synonymous SNPs identified in OsTPR was specifically associated with increase in USRB with superior grains being low chalk. Through haplotype mining, we further demonstrated the synergy of offering added yield advantage due to superior allele of OsTPR in elite materials with low glycaemic index (GI) property. We further validated the importance of OsTPR using recombinant inbred lines (RILs) population by introgressing a superior allele of OsTPR into elite materials resulted in raise in productivity in high amylose background. This confirmed a critical role for OsTPR in influencing yield while maintaining grain and nutritional quality.

Identification of stable pollen development related reference genes for accurate qRT-PCR analysis and morphological variations in autotetraploid and diploid rice

Autotetraploid rice exhibited hybrid vigor and greater genetic variation compared to diploid rice, but low pollen fertility is a major hindrance for its utilization. Our previous analysis revealed that large number of pollen fertility genes were exhibited down-regulation in autotetraploid rice. Hence, it is of utmost importance to reveal the expression patterns of pollen fertility genes with high accuracy. To find stable reference genes for autotetraploid rice, we compared the pollen development stages between diploid and autotetraploid rice, and 14 candidate genes were selected based on transcriptome analysis to evaluate their expression levels. Autotetraploid rice (i.e. Taichung65-4x) displayed lower seed set (40.40%) and higher percentage of abnormalities during the pollen development process than its diploid counterpart. To detect the candidate reference genes for pollen development of autotetraploid and diploid rice, we used five different algorithms, including NormFinder, BestKeeper, ΔCt method, geNorm and Re-Finder to evaluate their expression patterns stability. Consequently, we identified two genes, Cytochrome b5 and CPI, as the best candidate reference genes for qRT-PCR normalization in autotetraploid and diploid rice during pre-meiosis, meiosis, single microspore and bicellular pollen development stages. However, Cytochrome b5 was found to be the most stably expressed gene during different pollen development stages in autotetraploid rice. The results of our study provide a platform for subsequent gene expression analyses in autotetraploid rice, which could also be used in other polyploid plants.

MircroRNA Profiles of Early Rice Inflorescence Revealed a Specific miRNA5506 Regulating Development of Floral Organs and Female Megagametophyte in Rice

The developmental process of inflorescence and gametophytes is vital for sexual reproduction in rice. Multiple genes and conserved miRNAs have been characterized to regulate the process. The changes of miRNAs expression during the early development of rice inflorescence remain unknown. In this study, the analysis of miRNAs profiles in the early stage of rice inflorescence development identified 671 miRNAs, including 67 known and 44 novel differentially expressed miRNAs (DEMs). Six distinct clusters of miRNAs expression patterns were detected, and Cluster 5 comprised 110 DEMs, including unconserved, rice-specific osa-miR5506. Overexpression of osa-miR5506 caused pleiotropic abnormalities, including over- or under-developed palea, various numbers of floral organs and spikelet indeterminacy. In addition, the defects of ovaries development were frequently characterized by multiple megasporocytes, ovule-free ovary, megasporocyte degenerated and embryo sac degenerated in the transgenic lines. osa-miR5506 targeted REM transcription factor LOC_Os03g11370. Summarily, these results demonstrated that rice-specific osa-miR5506 plays an essential role in the regulation of floral organ number, spikelet determinacy and female gametophyte development in rice.

Designing Future Crops: Genomics-Assisted Breeding Comes of Age

Over the past decade, genomics-assisted breeding (GAB) has been instrumental in harnessing the potential of modern genome resources and characterizing and exploiting allelic variation for germplasm enhancement and cultivar development. Sustaining GAB in the future (GAB 2.0) will rely upon a suite of new approaches that fast-track targeted manipulation of allelic variation for creating novel diversity and facilitate their rapid and efficient incorporation in crop improvement programs. Genomic breeding strategies that optimize crop genomes with accumulation of beneficial alleles and purging of deleterious alleles will be indispensable for designing future crops. In coming decades, GAB 2.0 is expected to play a crucial role in breeding more climate-smart crop cultivars with higher nutritional value in a cost-effective and timely manner.

OsGRETCHENHAGEN3-2 modulates rice seed storability via accumulation of abscisic acid and protective substances

Seed storability largely determines the vigor of seeds during storage and is significant in agriculture and ecology. However, the underlying genetic basis remains unclear. In the present study, we report the cloning and characterization of the rice (Oryza sativa) indole-3-acetic acid (IAA)-amido synthetase gene GRETCHEN HAGEN3-2 (OsGH3-2) associated with seed storability. OsGH3-2 was identified by performing a genome-wide association study in rice germplasms with linkage mapping in chromosome substitution segment lines, contributing to the wide variation of seed viability in the populations after long periods of storage and artificial ageing. OsGH3-2 was dominantly expressed in the developing seeds and catalyzed IAA conjugation to amino acids, forming inactive auxin. Transgenic overexpression, knockout, and knockdown experiments demonstrated that OsGH3-2 affected seed storability by regulating the accumulation level of abscisic acid (ABA). Overexpression of OsGH3-2 significantly decreased seed storability, while knockout or knockdown of the gene enhanced seed storability compared with the wild-type. OsGH3-2 acted as a negative regulator of seed storability by modulating many genes related to the ABA pathway and probably subsequently late embryogenesis-abundant proteins at the transcription level. These findings shed light on the molecular mechanisms underlying seed storability and will facilitate the improvement of seed vigor by genomic breeding and gene-editing approaches in rice.

A heat stress responsive NAC transcription factor heterodimer plays key roles in rice grain filling

High temperature often leads to failure of grain filling in rice (Oryza sativa) causing yield loss, but the underlying mechanisms are still not elucidated. Here, we report that two genes encoding seed-specific NAM/ATAF/CUC (NAC) domain transcription factors, ONAC127 and ONAC129, are responsive to heat stress and involved in the grain filling process of rice. ONAC127 and ONAC129 are dominantly expressed in the pericarp and can form a heterodimer during rice grain filling. CRISPR/Cas9 induced mutants and overexpression lines were then generated to investigate the function of these two transcription factors. Interestingly, both knock-out and overexpression plants showed incomplete grain filling and shrunken grains, which became more severe under heat stress. Transcriptome analysis revealed that ONAC127 and ONAC129 mainly regulate stimulus response and nutrient transport. ChIP-seq analysis identified that the direct target genes of ONAC127 and ONAC129 in developing rice seeds include monosaccharide transporter gene OsMST6, sugar transporter gene OsSWEET4, calmodulin-like protein gene OsMSR2 and AP2/ERF factor gene OsEATB. These results suggest that ONAC127 and ONAC129 regulate grain filling by affecting sugar transportation and abiotic stress responses. Overall, this study demonstrates a transcriptional regulatory network with ONAC127 and ONAC129 coordinating multiple pathways to modulate seed development and heat stress responses at rice reproductive stages.

SWEET Transporters and the Potential Functions of These Sequences in Tea (Camellia sinensis)

Tea (Camellia sinensis) is an important economic beverage crop. Its flowers and leaves could be used as healthcare tea for its medicinal value. SWEET proteins were recently identified in plants as sugar transporters, which participate in diverse physiological processes, including pathogen nutrition, seed filling, nectar secretion, and phloem loading. Although SWEET genes have been characterized and identified in model plants, such as Arabidopsis thaliana and Oryza sativa, there is very little knowledge of these genes in C. sinensis. In this study, 28 CsSWEETs were identified in C. sinensis and further phylogenetically divided into four subfamilies with A. thaliana. These identified CsSWEETs contained seven transmembrane helixes (TMHs) which were generated by an ancestral three-TMH unit with an internal duplication experience. Microsynteny analysis revealed that the large-scale duplication events were the main driving forces for members from CsSWEET family expansion in C. sinensis. The expression profiles of the 28 CsSWEETs revealed that some genes were highly expressed in reproductive tissues. Among them, CsSWEET1a might play crucial roles in the efflux of sucrose, and CsSWEET17b could control fructose content as a hexose transporter in C. sinensis. Remarkably, CsSWEET12 and CsSWEET17c were specifically expressed in flowers, indicating that these two genes might be involved in sugar transport during flower development. The expression patterns of all CsSWEETs were differentially regulated under cold and drought treatments. This work provided a systematic understanding of the members from the CsSWEET gene family, which would be helpful for further functional studies of CsSWEETs in C. sinensis.

Arachis hypogaea gene expression atlas for fastigiata subspecies of cultivated groundnut to accelerate functional and translational genomics applications

Spatio-temporal and developmental stage-specific transcriptome analysis plays a crucial role in systems biology-based improvement of any species. In this context, we report here the Arachis hypogaea gene expression atlas (AhGEA) for the world's widest cultivated subsp. fastigiata based on RNA-seq data using 20 diverse tissues across five key developmental stages. Approximately 480 million paired-end filtered reads were generated followed by identification of 81 901 transcripts from an early-maturing, high-yielding, drought-tolerant groundnut variety, ICGV 91114. Further, 57 344 genome-wide transcripts were identified with ≥1 FPKM across different tissues and stages. Our in-depth analysis of the global transcriptome sheds light into complex regulatory networks namely gravitropism and photomorphogenesis, seed development, allergens and oil biosynthesis in groundnut. Importantly, interesting insights into molecular basis of seed development and nodulation have immense potential for translational genomics research. We have also identified a set of stable expressing transcripts across the selected tissues, which could be utilized as internal controls in groundnut functional genomics studies. The AhGEA revealed potential transcripts associated with allergens, which upon appropriate validation could be deployed in the coming years to develop consumer-friendly groundnut varieties. Taken together, the AhGEA touches upon various important and key features of cultivated groundnut and provides a reference for further functional, comparative and translational genomics research for various economically important traits.

Genome-wide association studies of ionomic and agronomic traits in USDA mini core collection of rice and comparative analyses of different mapping methods

BACKGROUND

Rice is an important human staple food vulnerable to heavy metal contamination leading to serious concerns. High yield with low heavy metal contamination is a common but highly challenging goal for rice breeders worldwide due to lack of genetic knowledge and markers.

RESULTS

To identify candidate QTLs and develop molecular markers for rice yield and heavy metal content, a total of 191 accessions from the USDA Rice mini-core collection with over 3.2 million SNPs were employed to investigate the QTLs. Sixteen ionomic and thirteen agronomic traits were analyzed utilizing two univariate (GLM and MLM) and two multivariate (MLMM and FarmCPU) GWAS methods. 106, 47, and 97 QTLs were identified for ionomics flooded, ionomics unflooded, and agronomic traits, respectively, with the criterium of p-value < 1.53 × 10- 8, which was determined by the Bonferroni correction for p-value of 0.05. While 49 (~ 20%) of the 250 QTLs were coinciding with previously reported QTLs/genes, about 201 (~ 80%) were new. In addition, several new candidate genes involved in ionomic and agronomic traits control were identified by analyzing the DNA sequence, gene expression, and the homologs of the QTL regions. Our results further showed that each of the four GWAS methods can identify unique as well as common QTLs, suggesting that using multiple GWAS methods can complement each other in QTL identification, especially by combining univariate and multivariate methods.

CONCLUSIONS

While 49 previously reported QTLs/genes were rediscovered, over 200 new QTLs for ionomic and agronomic traits were found in the rice genome. Moreover, multiple new candidate genes for agronomic and ionomic traits were identified. This research provides novel insights into the genetic basis of both ionomic and agronomic variations in rice, establishing the foundation for marker development in breeding and further investigation on reducing heavy-metal contamination and improving crop yields. Finally, the comparative analysis of the GWAS methods showed that each method has unique features and different methods can complement each other.

Analyzing Differentially Expressed Genes and Pathways Associated with Pistil Abortion in Japanese Apricot via RNA-Seq

Reproduction is a critical stage in the flower development process, and its failure causes serious problems affecting fruit quality and yield. Pistil abortion is one of the main factors in unsuccessful reproduction and occurs in many fruit plants. In Japanese apricot, the problem of pistil abortion is very common and affects fruit quality and plant yield; however, its molecular mechanism is not clearly understood. Therefore, in the current study, we used RNA-Seq to identify the differentially expressed genes (DEGs) and pathways actively involved in pistil abortion. A total of 3882 differentially expressed genes were found after cutoff and pairwise comparison analysis. According to KEGG pathway analysis, plant hormone signaling transduction and metabolic pathways were found most significantly enriched in this study. A total of 60 transcription factor families such as MADS-box, NAC and TCP showed their role in this process. RT-qPCR assays confirmed that the expression levels were consistent with RNA-Seq results. This study provides an alternative to be considered for further studies and understanding of pistil abortion processes in Japanese apricot, and it provides a reference related to this issue for other deciduous fruit crops.

Enhanced photosynthetic activity in pak choi hybrids is associated with increased grana thylakoids in chloroplasts

Increased photosynthetic activity is closely linked to heterosis in plants, but the underlying molecular mechanisms remain elusive. Pak choi (Brassica rapa ssp. chinensis) is a widely grown vegetable in Asia, and the most commercial cultivars are F₁ hybrids. Here, the inbred pak choi lines WTC and 2Q, and their reciprocal F₁ hybrids WQ and QW, were used to characterize the increased photosynthetic activity in these hybrids at the physiological, cellular and molecular levels. We found that the hybrids had larger leaves, with more grana thylakoids. Additionally, these hybrids had significantly increased net photosynthetic rates (P_n ) under both saturating and low irradiance conditions. These data indicate that the increased photosynthetic activity in pak choi hybrids was associated with an improved photosynthetic mechanism and larger leaves. Next, we obtained genome-wide data using transcriptome and bisulfite sequencing. Gene ontology (GO) analysis showed that the differentially expressed genes among the parents and hybrids were mostly enriched in the 'photosynthesis', 'thylakoid', and 'chloroplast' categories, indicating that the increased number of grana thylakoids contributes to the enhanced photosynthetic capacity in hybrids. Furthermore, we found that the increased number of grana thylakoids was associated with the upregulation of light-harvesting complex of photosystem II 1 (BrLhcb1). Yeast one-hybrid and transient assay showed that the BrLhcb1 promoter was directly bound by CIRCADIAN CLOCK ASSOCIATED 1 (BrCCA1), resulting in increased BrLhcb1 expression and enhanced carbon fixation in hybrids. Finally, our findings provide new insight into molecular mechanisms underlying enhanced photosynthesis in pak choi hybrids.

New insight into comprehensive analysis of INDETERMINATE DOMAIN (IDD) gene family in rice

The INDETERMINATE DOMAIN (IDD) transcription factor (TF), as a family of plant-specific zinc-finger proteins, regulates a variety of development processes and abiotic stresses in plants. IDD genes have been identified and characterized in other plants, however, the rice IDD family genes have not been investigated at genome-wide. In this study, 15 OsIDD genes were identified in rice genome and phylogenetically classified into two groups. Conserved motifs and potential interaction protein analysis about OsIDD proteins were carried out. Exon-intron structures, cis-acting elements and expression profiles of OsIDD genes were also examined. Exon-intron structures analysis revealed that overall structures of OsIDD genes were relatively conserved although they contained different numbers of introns. Cis-acting elements analysis suggested that most OsIDD gene transcripts could be induced by various abiotic stresses and phytohormones. The expression patterns of OsIDD genes were detected by qRT-PCR under cold and drought conditions, and by exogenous auxin (2,4-D), gibberellin (GA₃), and abscisic acid (ABA) treatments, respectively. The results showed that the OsIDDs might play essential roles under abiotic stresses and hormone responses. Distinct expression profiles in tissues/organs suggested that OsIDDs might be involved in different development processes in rice. More interestingly, the prediction of protein-protein interactions (PPIs) revealed OsIDDs could cooperate with some histone modifiers. Yeast two-hybrid assays were performed and confirmed it. Collectively, these results provide a foundation for further elucidation on the molecular mechanisms of OsIDD genes and advance our understanding of their biological function in rice.

Cloning and functional identification of a Chilo suppressalis-inducible promoter of rice gene, OsHPL2

BACKGROUND

Promoters play a key role in driving insect-resistant genes during breeding of transgenic plants. In current transgenic procedures for breeding rice resistance to striped stem borer (Chilo suppressalis Walker, SSB), the constitutive promoter is used to drive the insect-resistant gene. To reduce the burden of constitutive promoters on plant growth, isolation and identification of insect-inducible promoters are particularly important. However, few promoters are induced specifically by insect feeding.

RESULTS

We found rice hydroperoxide lyase gene (OsHPL2) (LOC_Os02g12680) was upregulated after feeding by SSB. We subsequently cloned the promoter of OsHPL2 and analysed its expression pattern using the β-glucuronidase (GUS) reporter gene. Histochemical assays and quantitative analyses of GUS activity confirmed that P _HPL2 :GUS was activated by SSB, but did not respond to brown planthopper (Nilaparvata lugens Stål, BPH) infestation, mechanical wounding or phytohormone treatments. A series of 5' truncated assays were conducted and three positive regulatory regions (-1452 to -1213, -903 to -624, and -376 to -176) induced by SSB infestation were identified. P2R123-min 35S and P2TR2-min 35S promoters linked with cry1C of transgenic plants showed the highest levels of Cry1C protein expression and SSB larval mortality.

CONCLUSION

We identified an SSB-inducible promoter and three positive internal regions. Transgenic rice plants with the OsHPL2 promoter and its positive regions driving cry1C exhibited the expected larvicidal effect on SSB. Our study is the first report of an SSB-inducible promoter that could be used as a potential resource for breeding insect-resistant transgenic crops. © 2020 Society of Chemical Industry.

Dynamic transcriptome and metabolome analyses of two types of rice during the seed germination and young seedling growth stages

Background

Seed germination and young seedling growth are important agricultural traits for developing populations of both irrigated and directly seeded rice. Previous studies have focused on the identification of QTLs. However, there are few studies on the metabolome or transcriptome of germination and young seedling growth in rice.

Results

Here, an indica rice and a japonica rice were used as materials, and the transcripts and metabolites were detected during the germination and young seedling growth periods on a large scale by using RNA sequencing and a widely targeted metabolomics method, respectively. Fourteen shared transcripts and 15 shared metabolites that were continuously differentially expressed in the two materials were identified and may be essential for seed germination and young seedling growth. Enrichment analysis of differentially expressed genes in transcriptome expression profiles at different stages indicated that cell wall metabolism, lipid metabolism, nucleotide degradation, amino acid, etc., were enriched at 0–2 days, and most of the results are consistent with those of previous reports. Specifically, phenylpropanoid biosynthesis and glutathione metabolism were continuously enriched during the seed germination and young seedling growth stages. Next, KO enrichment analysis was conducted by using the differentially expressed genes of the two materials at 2, 3 and 4 days. Fourteen pathways were enriched. Additionally, 44 differentially expressed metabolites at 2, 3 and 4 days were identified. These metabolites may be responsible for the differences in germination and young seedling growth between the two materials. Further attention was focused on the ascorbate–glutathione pathway, and it was found that differences in ROS-scavenging abilities mediated by some APX, GPX and GST genes may be directly involved in mediating differences in the germination and young seedling growth speed of the two materials.

Conclusions

In summary, these results may enhance the understanding of the overall mechanism of seed germination and young seedling growth, and the outcome of this study is expected to facilitate rice breeding for direct seeding.