Genome-wide transcriptome profiling provides insights into panicle development of rice (Oryza sativa L.)

Genome-wide association study of rice (Oryza sativa L.) inflorescence architecture

Rice yield strongly depends on panicle size and architecture but the genetics underlying these traits and their coordination with environmental cues through various signaling pathways have remained elusive. A genome-wide association study (GWAS) was performed to pinpoint the underlying genetic determinants for rice panicle architecture by analyzing 20 panicle-related traits using a data set consisting of 44,100 SNPs. We defined QTL windows around significant SNPs by the rate of LD decay for each chromosome and used these windows to identify putative candidate genes associated with the trait. Using a publicly available RNA-seq data set we performed analyses to identify the differentially expressed genes between stem and panicle with putative functions in panicle architecture. In total, 52 significant SNPs were identified, corresponding to 41 unique QTLs across the 12 rice chromosomes, with the most signals appearing on chromosome 1 (nine associated SNPs), and seven significant SNPs for each of chromosomes 8 and 12. Some novel genes such as Ankyrin, Duf, Kinesin and Brassinosteroid insensitive were found to be associated with panicle size. A haplotype analysis showed that genetic variation in haplotypes qMIL2 and qNSBBH21 were related to two traits, MIL, the greatest distance between two nodes on the rachis, and NSBBH, the number of primary branches in the bottom half of a panicle, respectively. Analysis of epistatic interactions revealed a marker affecting clustered traits. Several QTLs were identified on different chromosomes for the first time which may explain the phenotypic diversity of rice panicle architecture we observe in our collection of accessions. The identified candidate genes and haplotypes could be used in marker-assisted selection to improve rice yield through gene pyramiding.

Virulence effectors encoded in the rice yellow dwarf phytoplasma genome participate in pathogenesis

Bacteria-like phytoplasmas alternate between plant and insect hosts, secreting proteins that disrupt host development. In this study, we sequenced the complete genome of "Candidatus Phytoplasma oryzae" strain HN2022, associated with rice yellow dwarf (RYD) disease, using PacBio HiFi technology. The strain was classified within the 16Sr XI-B subgroup. Through SignalP v5.0 for prediction and subsequent expression analysis of secreted proteins in Nicotiana benthamiana and rice (Oryza sativa L.), we identified the key virulence effector proteins RY348 and RY378. RY348, a homolog of Secreted Aster Yellows Phytoplasma Effector 54 (SAP54), targets and degrades the MADS-box transcription factors MADS1 and MADS15, causing pollen sterility. Meanwhile, RY378 impacts the strigolactone and auxin signaling pathways, substantially increasing tillering. These findings offer insights into the interactions between plants and phytoplasmas.

GASA Proteins: Review of Their Functions in Plant Environmental Stress Tolerance

Gibberellic acid-stimulated Arabidopsis (GASA) gene family is a class of functional cysteine-rich proteins characterized by an N-terminal signal peptide and a C-terminal-conserved GASA domain with 12 invariant cysteine (Cys) residues. GASA proteins are widely distributed among plant species, and the majority of them are involved in the signal transmission of plant hormones, the regulation of plant development and growth, and the responses to different environmental constraints. To date, their action mechanisms are not completely elucidated. This review reports an overview of the diversity, structure, and subcellular localization of GASA proteins, their involvement in hormone crosstalk and redox regulation during development, and plant responses to abiotic and biotic stresses. Knowledge of this complex regulation can be a contribution to promoting multiple abiotic stress tolerance with potential agricultural applications through the engineering of genes encoding GASA proteins and the production of transgenic plants.

Ribosome profiling reveals the translational landscape and allele-specific translational efficiency in rice

Translational regulation is a critical step in the process of gene expression and governs the synthesis of proteins from mRNAs. Many studies have revealed translational regulation in plants in response to various environmental stimuli. However, there have been no studies documenting the comprehensive landscape of translational regulation and allele-specific translational efficiency in multiple plant tissues, especially those of rice, a main staple crop that feeds nearly half of the world's population. Here we used RNA sequencing and ribosome profiling data to analyze the transcriptome and translatome of an elite hybrid rice, Shanyou 63 (SY63), and its parental varieties Zhenshan 97 and Minghui 63. The results revealed that gene expression patterns varied more among tissues than among varieties at the transcriptional and translational levels. We identified 3392 upstream open reading frames (uORFs), and the uORF-containing genes were enriched in transcription factors. Only 668 of 13 492 long non-coding RNAs could be translated into peptides. Finally, we discovered numerous genes with allele-specific translational efficiency in SY63 and demonstrated that some cis-regulatory elements may contribute to allelic divergence in translational efficiency. Overall, these findings may improve our understanding of translational regulation in rice and provide information for molecular breeding research.

Transcriptome and metabolome profiling of the medicinal plant Veratrum mengtzeanum reveal key components of the alkaloid biosynthesis

Veratrum mengtzeanum is the main ingredient for Chinese folk medicine known as "Pimacao" due to its unique alkaloids. A diverse class of plant-specific metabolites having key pharmacological activities. There are limited studies on alkaloid synthesis and its metabolic pathways in plants. To elucidate the alkaloid pathway and identify novel biosynthetic enzymes and compounds in V. mengtzeanum, transcriptome and metabolome profiling has been conducted in leaves and roots. The transcriptome of V. mengtzeanum leaves and roots yielded 190,161 unigenes, of which 33,942 genes expressed differentially (DEGs) in both tissues. Three enriched regulatory pathways (isoquinoline alkaloid biosynthesis, indole alkaloid biosynthesis and tropane, piperidine and pyridine alkaloid biosynthesis) and a considerable number of genes such as AED3-like, A4U43, 21 kDa protein-like, 3-O-glycotransferase 2-like, AtDIR19, MST4, CASP-like protein 1D1 were discovered in association with the biosynthesis of alkaloids in leaves and roots. Some transcription factor families, i.e., AP2/ERF, GRAS, NAC, bHLH, MYB-related, C3H, FARI, WRKY, HB-HD-ZIP, C2H2, and bZIP were also found to have a prominent role in regulating the synthesis of alkaloids and steroidal alkaloids in the leaves and roots of V. mengtzeanum. The metabolome analysis revealed 74 significantly accumulated metabolites, with 55 differentially accumulated in leaves compared to root tissues. Out of 74 metabolites, 18 alkaloids were highly accumulated in the roots. A novel alkaloid compound viz; 3-Vanilloylygadenine was discovered in root samples. Conjoint analysis of transcriptome and metabolome studies has also highlighted potential genes involved in regulation and transport of alkaloid compounds. Here, we have presented a comprehensive metabolic and transcriptome profiling of V. mengtzeanum tissues. In earlier reports, only the roots were reported as a rich source of alkaloid biosynthesis, but the current findings revealed both leaves and roots as significant manufacturing factories for alkaloid biosynthesis.

Transcriptomics and metabolomics reveal tolerance new mechanism of rice roots to Al stress

The prevalence of soluble aluminum (Al) ions is one of the major limitations to crop production worldwide on acid soils. Therefore, understanding the Al tolerance mechanism of rice and applying Al tolerance functional genes in sensitive plants can significantly improve Al stress resistance. In this study, transcriptomics and metabolomics analyses were performed to reveal the mechanism of Al tolerance differences between two rice landraces (Al-tolerant genotype Shibanzhan (KR) and Al-sensitive genotype Hekedanuo (MR) with different Al tolerance. The results showed that DEG related to phenylpropanoid biosynthesis was highly enriched in KR and MR after Al stress, indicating that phenylpropanoid biosynthesis may be closely related to Al tolerance. E1.11.1.7 (peroxidase) was the most significant enzyme of phenylpropanoid biosynthesis in KR and MR under Al stress and is regulated by multiple genes. We further identified that two candidate genes Os02g0770800 and Os06g0521900 may be involved in the regulation of Al tolerance in rice. Our results not only reveal the resistance mechanism of rice to Al stress to some extent, but also provide a useful reference for the molecular mechanism of different effects of Al poisoning on plants.

Comprehensive Transcriptome Analysis of GS3 Near-Isogenic Lines During Panicle Development in Rice (Oryza sativa L.)

Panicle architecture is an important agronomic trait in rice that affects rice yields and quality. The GRAIN SIZE 3 (GS3) locus has been identified as a major quantitative trait locus (QTL) affecting grain length and weight. The current understanding of the function of the GS3 gene, especially concerning the regulatory mechanism of panicle development, is still in its infancy. In this study, we generated GS3 near-isogenic lines (NILs) by successive crossing and backcrossing of TD70 (large grain) with Kasalath (small grain), using Kasalath as the recurrent parent. To identify potential transcription dynamic changes in rice panicle formation and grain shape, we deeply analyzed transcriptional profiles for the NILs (NIL-GS3 and NIL-gs3) at three different panicle developmental stages (S, M, and L). A total of 887, 1,768, and 1,478 differentially expressed genes (DEGs) were identified at stages S, M, and L, respectively. We also found 542 differential expressed long non-coding RNAs (lncRNAs). Co-expression analysis further revealed significant clusters associated with different development periods in NIL-gs3 lines. Gene Ontology and KEGG enrichment analysis revealed G-protein signaling and hormones pathway were successively activated at the M and L stages of NIL-gs3, which indicated activation of the G-protein signaling pathway might trigger the down-streaming hormone signaling transduction. we found that other hormones such ABA, Auxin, CK were significantly enriched in the L stage in the NIL-gs3. We highlighted the synergistic interplay of G-protein and multiple hormones signaling pathways and their essential roles in regulating rice panicle formation and the grain shape. Our study provides an invaluable resource for further molecular mechanistic studies that affect rice grain size and provide new insight for directed selection by marker-assisted backcross breeding.

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.

Diverse Roles of MAX1 Homologues in Rice

Cytochrome P450 enzymes encoded by MORE AXILLARY GROWTH1 (MAX1)-like genes produce most of the structural diversity of strigolactones during the final steps of strigolactone biosynthesis. The diverse copies of MAX1 in Oryza sativa provide a resource to investigate why plants produce such a wide range of strigolactones. Here we performed in silico analyses of transcription factors and microRNAs that may regulate each rice MAX1, and compared the results with available data about MAX1 expression profiles and genes co-expressed with MAX1 genes. Data suggest that distinct mechanisms regulate the expression of each MAX1. Moreover, there may be novel functions for MAX1 homologues, such as the regulation of flower development or responses to heavy metals. In addition, individual MAX1s could be involved in specific functions, such as the regulation of seed development or wax synthesis in rice. Our analysis reveals potential new avenues of strigolactone research that may otherwise not be obvious.

The Candidate Genes Underlying a Stably Expressed QTL for Low Temperature Germinability in Rice (Oryza sativa L.)

BACKGROUND

Direct seeding is an efficient cultivation technique in rice. However, poor low temperature germinability (LTG) of modern rice cultivars limits its application. Identifying the genes associated with LTG and performing molecular breeding is the fundamental way to address this issue. However, few LTG QTLs have been fine mapped and cloned so far.

RESULTS

In the present study, the LTG evaluation of 375 rice accessions selected from the Rice Diversity Panel 2 showed that there were large LTG variations within the population, and the LTG of Indica group was significantly higher than that of Japonica and Aus groups (p < 0.01). In total, eleven QTLs for LTG were identified through genome-wide association study (GWAS). Among them, qLTG_sRDP2-3/qLTG_JAP-3, qLTG_AUS-3 and qLTG_sRDP2-12 are first reported in the present study. The QTL on chromosome 10, qLTG_sRDP2-10a had the largest contribution to LTG variations in 375 rice accessions, and was further validated using single segment substitution line (SSSL). The presence of qLTG_sRDP2-10a could result in 59.8% increase in LTG under 15 °C low temperature. The expression analysis of the genes within qLTG_sRDP2-10a region indicated that LOC_Os10g22520 and LOC_Os10g22484 exhibited differential expression between the high and low LTG lines. Further sequence comparisons revealed that there were insertion and deletion sequence differences in the promoter and intron region of LOC_Os10g22520, and an about 6 kb variation at the 3' end of LOC_Os10g22484 between the high and low LTG lines, suggesting that the sequence variations of the two genes could be the cause for their differential expression in high and low LTG lines.

CONCLUSION

Among the 11 QTLs identified in this study, qLTG_sRDP2-10a could also be detected in other three studies using different germplasm under different cold environments. Its large effect and stable expression make qLTG_sRDP2-10a particularly valuable in rice breeding. The two genes, LOC_Os10g22484 and LOC_Os10g22520, were considered as the candidate genes underlying qLTG_sRDP2-10a. Our results suggest that integrating GWAS and SSSL can facilitate identification of QTL for complex traits in rice. The identification of qLTG_sRDP2-10a and its candidate genes provide a promising source for gene cloning of LTG and molecular breeding for LTG in rice.

Transcriptome Analysis Provides Insights into Grain Filling in Foxtail Millet (Setaria italica L.)

Grain filling is an importantly developmental process which is associated with the yield and quality of foxtail millet (Setaria italic L.). However, the molecular mechanisms of grain filling are rarely reported in foxtail millet. In our study, RNA-seq was performed to investigate the transcriptional dynamics and identify the key genes involved in grain filling in foxtail millet at five different developmental stages. A total of 11,399 differentially expressed genes (DEGs), including 902 transcription factors (TFs), were identified. Certain important genes involved in grain filling were discovered through a function annotation and temporal expression patterns analysis. These genes included genes associated with starch biosynthesis, cell-wall invertases, hormone signal transduction, and polyamine metabolism pathways. The expression levels of seven randomly selected DEGs were validated by a quantitative real-time polymerase chain reaction (qRT-PCR). This study provides the first insight into the changes in the gene expression of grain filling at different developmental stages in foxtail millet. These results could help understand the complex molecular mechanisms of the panicle formation in foxtail millet and other cereal crops.

Metabolomics and transcriptomics reveal defense mechanism of rice (Oryza sativa) grains under stress of 2,2',4,4'-tetrabromodiphenyl ether

2,2',4,4'-Tetrabromodiphenyl ether (BDE-47), a predominant polybrominated diphenyl ether (PBDE), has received extensive attention for its potential environmental impact. An integrated study of metabolomics and transcriptomics was conducted on two rice (Oryza sativa) cultivars, Lianjing-7 (LJ-7) and Yongyou-9 (YY-9), which have been identified as tolerant and sensitive cultivars to BDE-47, respectively. The objective was to investigate the molecular mechanisms of their different ability to tolerate BDE-47. Both rice plants were cultivated to maturity in soils containing three concentrations of BDE-47 (10, 20, and 50 mg/kg). Metabolomic analyses of rice grains identified 65 metabolites in LJ-7 and 45 metabolites in YY-9, including amino acids, saccharides, organic acids, fatty acids, and secondary metabolites. In the tolerant cultivar LJ-7 exposed to 50 mg/kg BDE-47, concentrations of most of the metabolites increased significantly, with α-ketoglutaric acid increased by 20-fold and stigmastanol increased by 12-fold. In the sensitive cultivar YY-9, the concentrations of most metabolites increased after the plant was exposed to 1 and 10 mg/kg BDE-47 but decreased after the plant was exposed to 50 mg/kg BDE-47. Transcriptomic data demonstrated that regulation of gene expressions was affected most in LJ-7 exposed to 50 mg/kg BDE-47 (966 genes up-regulated and 620 genes down-regulated) and in YY-9 exposed to 10 mg/kg BDE-47 (85 genes up-regulated and 291 genes down-regulated), in good accordance with the observed metabolic alternation in the two cultivars. Analyses of metabolic pathways and KEGG enrichment revealed that many biological processes, including energy consumption and biosynthesis, were perturbed in the two rice cultivars by BDE-47. A majority of metabolites and genes involved in dominating pathways of energy consumption (e.g., tricarboxylic acid cycle) and the biosynthesis (e.g., metabolism of saccharides and amino acids) were enhanced in LJ-7 by BDE-47. In contrast, energy consumption was increased while biosynthetic processes were inhibited in YY-9 by BDE-47, which could lead to the sensitivity of YY-9 to BDE-47. The combined results suggest that the different defensive abilities of these two rice cultivars in response to BDE-47 could be attributed to their differences in energy-consumption strategy and biosynthesis of nutritional components in grains. This study provides a useful reference for rice cultivation in PBDE-polluted areas.

A systematic in silico prediction of gibberellic acid stimulated GASA family members: A novel small peptide contributes to floral architecture and transcriptomic changes induced by external stimuli in rice

The GASA (GA-stimulated Arabidopsis) gene family is highly specific to plants, signifying a crucial role in plant growth and development. Herein, we retrieved 119 GASA genes in 10 different plant species in two major lineages (monocots and eudicots). Further, in the phylogenetic tree we classified these genes into four well-conserved subgroups. All the proteins contain a conserved GASA domain with similar characteristics and a highly specific 12-cysteine residue of the C-terminus position. According to the global microarray data and qRT-PCR based analysis, the OsGASA gene family was dominantly expressed in the seedling and transition phase of floral stages. Despite this, OsGASA genes profoundly contribute to rice grain size and length, whereas the highest abundance of transcript level was noticed in stage-2 (Inf 6, 3.0-cm-long spikelet) and stage-3 (Inf 7, 5.0-cm-long spikelet) under GA treatment during panicle formation. Additionally, the maximum expression level of these genes was recorded in response to GA and ABA in young seedlings. Further, in response to abiotic stresses, OsGASA1/8/10 was up- regulated by salt, OsGASA2/5/7 by drought, OsGASA3/6 by cold, and OsGASA4/9 by heat stress. With the exception of OsGASA4, the higher transcription levels of all the other GASA genes were induced by Cd and Cr metal stresses (8-10 fold changes) at various time points. Finally, the GO ontology analysis of GASAs revealed the biological involvement in the GA-mediated signaling pathway and abiotic stresses. Prominently, most of these proteins are localized in cellular components such as the cell wall and extracellular region, where the molecular functions such as ATP binding and protein binding were observed. These results imply that GASAs are significantly involved in rice panicle developmental stages, responses to external stimuli, and hormones.