MiR529a controls plant height, tiller number, panicle architecture and grain size by regulating SPL target genes in rice (Oryza sativa L.)

Pleiotropic Effects of miR5504 Underlying Plant Height, Grain Yield and Quality in Rice

MicroRNAs (miRNAs) are known to play critical roles in regulating rice agronomic traits through mRNA cleavage or translational repression. Our previous study indicated that miR5504 regulates plant height by affecting cell proliferation and expansion. Here, the two independent homozygous mir5504 mutants (CR1 and CR2) and overexpression lines (OE1 and OE2) were further used to investigate the functions of miR5504. The panicle length, 1000-grain weight and grain yield per plant of miR5504-OE lines were identical to those of Nipponbare (NIP), but the 1000-grain weight of mir5504 mutants was reduced by about 10% and 9%, respectively. Meanwhile, the grain width and thickness of mir5504 mutants decreased significantly by approximately 10% and 11%, respectively. Moreover, the cytological results revealed a significant decrease in cell number along grain width direction and cell width in spikelet in mir5504, compared with those in NIP. In addition, several major storage substances of the rice seeds were measured. Compared to NIP, the amylose content of the mir5504 seeds was noticeably decreased, leading to an increase of nearly 10 mm in gel consistency (GC) in mir5504 lines. Further investigation confirmed that LOC_Os08g16914 was the genuine target of miR5504: LOC_Os08g16914 over-expression plants phenocopied the mir5504 mutants. This study provides insights into the role of miR5504 in rice seed development.

Genome-wide analysis of the SPL transcription factor family and its response to water stress in sunflower (Helianthus annuus)

SPL (SQUAMOSA promoter binding proteins-like) are plant-specific transcription factors that play essential roles in a variety of developmental processes as well as the ability to withstand biotic and abiotic stresses. To date, numerous species have been investigated for the SPL gene family, but so far, no SPL family genes have been thoroughly identified and characterized in the sunflower (Helianthus annuus). In this study, 25 SPL genes were identified in the sunflower genome and were unevenly distributed on 11 chromosomes. According to phylogeny analysis, 59 SPL genes from H. annuus, O. sativa, and A. thaliana were clustered into seven groups. Furthermore, the SPL genes in groups-I and II were demonstrated to be potential targets of miR156. Synteny analysis showed that 7 paralogous gene pairs exist in HaSPL genes and 26 orthologous gene pairs exist between sunflower and rice, whereas 21 orthologous gene pairs were found between sunflower and Arabidopsis. Segmental duplication appears to have played a vital role in the expansion processes of sunflower SPL genes, and because of selection pressure, all duplicated genes have undergone purifying selection. Tissue-specific gene expression analysis of the HaSBP genes proved their diverse spatiotemporal expression patterns, which were predominantly expressed in floral organs and differentially expressed in stem, axil, and root tissues. The expression pattern of HaSPL genes under water stress showed broad involvement of HaSPLs in the response to flood and drought stresses. This genome-wide identification investigation provides detailed information on the sunflower SPL transcription factor gene family and establishes a strong platform for future research on sunflower responses to abiotic stress tolerance.

A multi-level approach reveals key physiological and molecular traits in the response of two rice genotypes subjected to water deficit at the reproductive stage

Rice is more vulnerable to drought than maize, wheat, and sorghum because its water requirements remain high throughout the rice life cycle. The effects of drought vary depending on the timing, intensity, and duration of the events, as well as on the rice genotype and developmental stage. It can affect all levels of organization, from genes to the cells, tissues, and/or organs. In this study, a moderate water deficit was applied to two contrasting rice genotypes, IAC 25 and CIRAD 409, during their reproductive stage. Multi-level transcriptomic, metabolomic, physiological, and morphological analyses were performed to investigate the complex traits involved in their response to drought. Weighted gene network correlation analysis was used to identify the specific molecular mechanisms regulated by each genotype, and the correlations between gene networks and phenotypic traits. A holistic analysis of all the data provided a deeper understanding of the specific mechanisms regulated by each genotype, and enabled the identification of gene markers. Under non-limiting water conditions, CIRAD 409 had a denser shoot, but shoot growth was slower despite better photosynthetic performance. Under water deficit, CIRAD 409 was weakly affected regardless of the plant level analyzed. In contrast, IAC 25 had reduced growth and reproductive development. It regulated transcriptomic and metabolic activities at a high level, and activated a complex gene regulatory network involved in growth-limiting processes. By comparing two contrasting genotypes, the present study identified the regulation of some fundamental processes and gene markers, that drive rice development, and influence its response to water deficit, in particular, the importance of the biosynthetic and regulatory pathways for cell wall metabolism. These key processes determine the biological and mechanical properties of the cell wall and thus influence plant development, organ expansion, and turgor maintenance under water deficit. Our results also question the genericity of the antagonism between morphogenesis and organogenesis observed in the two genotypes.

Rice miR5504 regulates plant height by affecting cell proliferation and expansion

miRNAs play critical roles in the regulation of plant growth and development by cleaving mRNA or repressing transcription. In our previous study, miR5504 with unknown functions was captured by small RNA sequencing. Here, the function and characters of miR5504 were extensively analyzed using CRISPR/Cas9, overexpression strategy, Northern blot, cytological analysis, and transcriptomics analysis. We found that the dwarf phenotype of mir5504 mutants (mir5504-1 and mir5504-2) appeared on 35-day seedlings and became more apparent at the mature stage. The cytological results showed a substantial decrease in the vascular bundle number, cell number and cell length in the mir5504 mutant compared with NIP. In addition, we found that miR5504 regulated plant height by targeting LOC_Os08g16914. The results of RNA-seq revealed that numerous biological processes were mainly enriched, including DNA-binding transcription factor activity, transferase activity, regulation of transcription, metabolic process, and protein binding. Meanwhile, KEEG analysis showed that numerous proteins were associated with cellular processes and metabolism pathways. Taken together, miR5504 may be involved in the regulation of plant height by affecting cell expansion and division of internode in rice.

Revitalizing miRNAs mediated agronomical advantageous traits improvement in rice

One of the key enigmas in conventional and modern crop improvement programmes is how to introduce beneficial traits without any penalty impairment. Rice (Oryza sativa L.), among the essential staple food crops grown and utilized worldwide, needs to improve genotypes in multifaceted ways. With the global view to feed ten billion under the climatic perturbation, only a potent functional master regulator can withstand with hope for the next green revolution and food security. miRNAs are such, miniature, fine tuners for crop improvement and provide a value addition in emerging technologies, namely large-scale genotyping, phenotyping, genome editing, marker-assisted selection, and genomic selection, to make rice production feasible. There has been surplus research output generated since the last decade on miRNAs in rice, however, recent functional knowledge is limited to reaping the benefits for conventional and modern improvements in rice to avoid ambiguity and redundancy in the generated data. Here, we present the latest functional understanding of miRNAs in rice. In addition, their biogenesis, intra- and inter-kingdom signaling and communication, implication of amiRNAs, and consequences upon integration with CRISPR-Cas9. Further, highlights refer to the application of miRNAs for rice agronomical trait improvements, broadly classified into three functional domains. The majority of functionally established miRNAs are responsible for growth and development, followed by biotic and abiotic stresses. Tabular cataloguing reveals and highlights two multifaceted modules that were extensively studied. These belong to miRNA families 156 and 396, orchestrate multifarious aspects of advantageous agronomical traits. Moreover, updated and exhaustive functional aspects of different supplemental miRNA modules that would strengthen rice improvement are also being discussed.

Rice microRNA156/529-SQUAMOSA PROMOTER BINDING PROTEIN-LIKE7/14/17 modules regulate defenses against bacteria

Rice (Oryza sativa L.) microRNA156/529-SQUAMOSA PROMOTER BINDING PROTEIN-LIKE7/14/17 (miR156/529-SPL7/14/17) modules have pleiotropic effects on many biological pathways. OsSPL7/14 can interact with DELLA protein SLENDER RICE1 (SLR1) to modulate gibberellin acid (GA) signal transduction against the bacterial pathogen Xanthomonas oryzae pv. oryzae. However, whether the miR156/529-OsSPL7/14/17 modules also regulate resistance against other pathogens is unclear. Notably, OsSPL7/14/17 functioning as transcriptional activators, their target genes, and the corresponding downstream signaling pathways remain largely unexplored. Here, we demonstrate that miR156/529 play negative roles in plant immunity and that miR156/529-regulated OsSPL7/14/17 confer broad-spectrum resistance against 2 devastating bacterial pathogens. Three OsSPL7/14/17 proteins directly bind to the promoters of rice Allene Oxide Synthase 2 (OsAOS2) and NONEXPRESSOR OF PATHOGENESIS-RELATED GENES1 (OsNPR1) and activate their transcription, regulating jasmonic acid (JA) accumulation and the salicylic acid (SA) signaling pathway, respectively. Overexpression of OsAOS2 or OsNPR1 impairs the susceptibility of the osspl7/14/17 triple mutant. Exogenous application of JA enhances resistance of the osspl7/14/17 triple mutant and the miR156 overexpressing plants. In addition, genetic evidence confirms that bacterial pathogen-activated miR156/529 negatively regulate pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) responses, such as pattern recognition receptor Xa3/Xa26-initiated PTI. Our findings demonstrate that bacterial pathogens modulate miR156/529-OsSPL7/14/17 modules to suppress OsAOS2-catalyzed JA accumulation and the OsNPR1-promoted SA signaling pathway, facilitating pathogen infection. The uncovered miR156/529-OsSPL7/14/17-OsAOS2/OsNPR1 regulatory network provides a potential strategy to genetically improve rice disease resistance.

Cytological, transcriptome and miRNome temporal landscapes decode enhancement of rice grain size

BACKGROUND

Rice grain size (GS) is an essential agronomic trait. Though several genes and miRNA modules influencing GS are known and seed development transcriptomes analyzed, a comprehensive compendium connecting all possible players is lacking. This study utilizes two contrasting GS indica rice genotypes (small-grained SN and large-grained LGR). Rice seed development involves five stages (S1-S5). Comparative transcriptome and miRNome atlases, substantiated with morphological and cytological studies, from S1-S5 stages and flag leaf have been analyzed to identify GS proponents.

RESULTS

Histology shows prolonged endosperm development and cell enlargement in LGR. Stand-alone and comparative RNAseq analyses manifest S3 (5-10 days after pollination) stage as crucial for GS enhancement, coherently with cell cycle, endoreduplication, and programmed cell death participating genes. Seed storage protein and carbohydrate accumulation, cytologically and by RNAseq, is shown to be delayed in LGR. Fourteen transcription factor families influence GS. Pathway genes for four phytohormones display opposite patterns of higher expression. A total of 186 genes generated from the transcriptome analyses are located within GS trait-related QTLs deciphered by a cross between SN and LGR. Fourteen miRNA families express specifically in SN or LGR seeds. Eight miRNA-target modules display contrasting expressions amongst SN and LGR, while 26 (SN) and 43 (LGR) modules are differentially expressed in all stages.

CONCLUSIONS

Integration of all analyses concludes in a "Domino effect" model for GS regulation highlighting chronology and fruition of each event. This study delineates the essence of GS regulation, providing scope for future exploits. The rice grain development database (RGDD) ( www.nipgr.ac.in/RGDD/index.php ; https://doi.org/10.5281/zenodo.7762870 ) has been developed for easy access of data generated in this paper.

A Novel miRNA in Rice Associated with the Low Seed Setting Rate Symptom of Rice Stripe Virus

MicroRNAs play key regulatory roles in plant development. The changed pattern of miRNA expression is involved in the production of viral symptoms. Here, we showed that a small RNA, Seq119, a putative novel microRNA, is associated with the low seed setting rate, a viral symptom of rice stripe virus (RSV)-infected rice. The expression of Seq 119 was downregulated in RSV-infected rice. The overexpression of Seq119 in transgenic rice plants did not cause any obvious phenotypic changes in plant development. When the expression of Seq119 was suppressed in rice plants either by expressing a mimic target or by CRISPR/Cas editing, seed setting rates were extremely low, similar to the effects of RSV infection. The putative targets of Seq119 were then predicted. The overexpression of the target of Seq119 in rice caused a low seed setting rate, similar to that in Seq119-suppressed or edited rice plants. Consistently, the expression of the target was upregulated in Seq119-suppressed and edited rice plants. These results suggest that downregulated Seq119 is associated with the low seed setting rate symptom of the RSV in rice.

Transcriptome and Gene Co-Expression Network Analysis Identifying Differentially Expressed Genes and Signal Pathways Involved in the Height Development of Banana (Musa spp.)

Plant height is an important and valuable agronomic trait associated with yield and resistance to abiotic and biotic stresses. Dwarfism has positive effects on plant development and field management, especially for tall monocotyledon banana (Musa spp.). However, several key genes and their regulation mechanism of controlling plant height during banana development are unclear. In the present study, the popular cultivar 'Brazilian banana' ('BX') and its dwarf mutant ('RK') were selected to identify plant height-related genes by comparing the phenotypic and transcriptomic data. Banana seedlings with 3-4 leaves were planted in the greenhouse and field. We found that the third and fourth weeks are the key period of plant height development of the selected cultivars. A total of 4563 and 10507 differentially expressed genes (DEGs) were identified in the third and fourth weeks, respectively. Twenty modules were produced by the weighted gene co-expression network analysis (WGCNA). Eight modules were positively correlated with the plant height, and twelve other modules were negatively correlated. Combining with the analysis of DEGs and WGCNA, 13 genes in the signaling pathway of gibberellic acid (GA) and 7 genes in the signaling pathway of indole acetic acid (IAA) were identified. Hub genes related to plant height development were obtained in light of the significantly different expression levels (|log2FC| ≥ 1) at the critical stages. Moreover, GA3 treatment significantly induced the transcription expressions of the selected candidate genes, suggesting that GA signaling could play a key role in plant height development of banana. It provides an important gene resource for the regulation mechanism of banana plant development and assisted breeding of ideal plant architecture.

Molecular characterization of SPL gene family during flower morphogenesis and regulation in blueberry

The SPL gene is a plant-specific transcription factor involved in the regulation of plant growth and development, which have been identified in woody plants. The process of floral bud differentiation affects the timing of flowering and fruit set and regulates plant growth, however, the mechanism of regulation of flower development by SPL genes is less studied. In this study, 56 VcSPL genes were identified in the tetraploid blueberry. The VcSPL gene family was classified into six subfamilies, and analysis of cis-elements showed that VcSPL genes were regulated by light, phytohormones (abscisic acid, MeJA), and low temperature. In the evolutionary analysis, segmental replication may play an important role in VcSPL gene amplification. Interestingly, we also studied diploid blueberry (Bilberry), in which 24 SPL genes were identified, and 36 homologous pairs were found, suggesting a high degree of convergence in the syntenic relationship between blueberry (Vaccinium corymbosum L) and bilberry (Vaccinium darrowii). Based on the expression profile, VcSPL genes were expressed at high levels in flowers, shoots, and roots, indicating a diversity of gene functions. Then we selected 20 differentially-expressed SPL genes to further investigate the role of VcSPL in floral induction and initiation. It showed that the genes VcSPL40, VcSPL35, VcSPL45, and VcSPL53 may play a crucial role in the blueberry floral transition phase (from vegetative growth to flower initiation). These results provided important information for understanding and exploring the role of VcSPLs in flower morphogenesis and plant growth.

Identification of Alfalfa SPL gene family and expression analysis under biotic and abiotic stresses

The SQUAMOSA promoter binding-like protein (SPL) is a specific transcription factor that affects plant growth and development. The SPL gene family has been explored in various plants, but information about these genes in alfalfa is limited. This study, based on the whole genome data of alfalfa SPL, the fundamental physicochemical properties, phylogenetic evolution, gene structure, cis-acting elements, and gene expression of members of the MsSPL gene family were analyzed by bioinformatics methods. We identified 82 SPL sequences in the alfalfa, which were annotated into 23 genes, including 7 (30.43%) genes with four alleles, 10 (43.47%) with three, 3 (13.04%) with two, 3 (13.04%) with one allele. These SPL genes were divided into six groups, that are constructed from A. thaliana, M. truncatula and alfalfa. Chromosomal localization of the identified SPL genes showed arbitary distribution. The subcellular localization predictions showed that all MsSPL proteins were located in the nucleus. A total of 71 pairs of duplicated genes were identified, and segmental duplication mainly contributed to the expansion of the MsSPL gene family. Analysis of the Ka/Ks ratios indicated that paralogs of the MsSPL gene family principally underwent purifying selection. Protein-protein interaction analysis of MsSPL proteins were performed to predict their roles in potential regulatory networks. Twelve cis-acting elements including phytohormone and stress elements were detected in the regions of MsSPL genes. We further analyzed that the MsSPLs had apparent responses to abiotic stresses such as drought and salt and the biotic stress of methyl jasmonate. These results provide comprehensive information on the MsSPL gene family in alfalfa and lay a solid foundation for elucidating the biological functions of MsSPLs. This study also provides valuable on the regulation mechanism and function of MsSPLs in response to biotic and abiotic stresses.

Genome-Wide Analysis of SQUAMOSA-Promoter-Binding Protein-like Family in Flowering Pleioblastus pygmaeus

SQUAMOSA Promoter-Binding Protein-Like (SPL) family is well-known for playing an important role in plant growth and development, specifically in the reproductive process. Bamboo plants have special reproductive characteristics with a prolonged vegetative phase and uncertain flowering time. However, the underlying functions of SPL genes in reproductive growth are undisclosed in bamboo plants. In the study, a total of 28 SPLs were screened from an ornamental dwarf bamboo species, Pleioblastus pygmaeus. Phylogenetic analysis indicates that 183 SPLs from eight plant species can be classified into nine subfamilies, and the 28 PpSPLs are distributed among eight subfamilies. Homologous analysis shows that as many as 32 pairs of homologous genes were found between P. pygmaeus and rice, and 83 pairs were found between P. pygmaeus and Moso bamboo, whose Ka/Ks values are all <1. MiRNA target prediction reveals that 13 out of the 28 PpSPLs have recognition sites complementary to miRNA156. To screen the SPLs involved in the reproductive growth of bamboo plants, the mRNA abundance of the 28 PpSPLs was profiled in the different tissues of flowering P. pygmaeus and non-flowering plants by RNA-Seq. Moreover, the relative expression level of eight PpSPLs is significantly higher in flowering P. pygmaeus than that in non-flowering plants, which was also validated by RT-qPCR. Combined with phylogenetic analysis and homologous analysis, the eight significant, differentially expressed PpSPLs were identified to be associated with the reproductive process and flower organ development. Among them, there are four potential miRNA156-targeting PpSPLs involved in the flowering process. Of significant interest in the study is the identification of 28 SPLs and the exploration of four key flowering-related SPLs from P. pygmaeus, which provides a theoretic basis for revealing the underlying functions of SPLs in the reproductive growth of bamboo plants.

Non-coding RNAs fine-tune the balance between plant growth and abiotic stress tolerance

To survive in adverse environmental conditions, plants have evolved sophisticated genetic and epigenetic regulatory mechanisms to balance their growth and abiotic stress tolerance. An increasing number of non-coding RNAs (ncRNAs), including small RNAs (sRNAs) and long non-coding RNAs (lncRNAs) have been identified as essential regulators which enable plants to coordinate multiple aspects of growth and responses to environmental stresses through modulating the expression of target genes at both the transcriptional and posttranscriptional levels. In this review, we summarize recent advances in understanding ncRNAs-mediated prioritization towards plant growth or tolerance to abiotic stresses, especially to cold, heat, drought and salt stresses. We highlight the diverse roles of evolutionally conserved microRNAs (miRNAs) and small interfering RNAs (siRNAs), and the underlying phytohormone-based signaling crosstalk in regulating the balance between plant growth and abiotic stress tolerance. We also review current discoveries regarding the potential roles of ncRNAs in stress memory in plants, which offer their descendants the potential for better fitness. Future ncRNAs-based breeding strategies are proposed to optimize the balance between growth and stress tolerance to maximize crop yield under the changing climate.

Molecular evolution and functional modification of plant miRNAs with CRISPR

Gene editing using clustered regularly interspaced short palindromic repeat/CRISPR-associated proteins (CRISPR/Cas) has revolutionized biotechnology and provides genetic tools for medicine and life sciences. However, the application of this technology to miRNAs, with the function as negative gene regulators, has not been extensively reviewed in plants. Here, we summarize the evolution, biogenesis, and structure of miRNAs, as well as their interactions with mRNAs and computational models for predicting target genes. In addition, we review current advances in CRISPR/Cas for functional analysis and for modulating miRNA genes in plants. Extending our knowledge of miRNAs and their manipulation with CRISPR will provide fundamental understanding of the functions of plant miRNAs and facilitate more sustainable and publicly acceptable genetic engineering of crops.

Genetic and molecular factors in determining grain number per panicle of rice

It was suggested that the most effective way to improve rice grain yield is to increase the grain number per panicle (GN) through the breeding practice in recent decades. GN is a representative quantitative trait affected by multiple genetic and environmental factors. Understanding the mechanisms controlling GN has become an important research field in rice biotechnology and breeding. The regulation of rice GN is coordinately controlled by panicle architecture and branch differentiation, and many GN-associated genes showed pleiotropic effect in regulating tillering, grain size, flowering time, and other domestication-related traits. It is also revealed that GN determination is closely related to vascular development and the metabolism of some phytohormones. In this review, we summarize the recent findings in rice GN determination and discuss the genetic and molecular mechanisms of GN regulators.

OsSPL14 acts upstream of OsPIN1b and PILS6b to modulate axillary bud outgrowth by fine-tuning auxin transport in rice

As a multigenic trait, rice tillering can optimize plant architecture for the maximum agronomic yield. SQUAMOSA PROMOTER BINDING PROTEIN-LIKE14 (OsSPL14) has been demonstrated to be necessary and sufficient to inhibit rice branching, but the underlying mechanism remains largely unclear. Here, we demonstrated that OsSPL14, which is cleaved by miR529 and miR156, inhibits tillering by fine-tuning auxin transport in rice. RNA interference of OsSPL14 or miR529 and miR156 overexpression significantly increased the tiller number, whereas OsSPL14 overexpression decreased the tiller number. Histological analysis revealed that the OsSPL14-overexpressing line had normal initiation of axillary buds but inhibited outgrowth of tillers. Moreover, OsSPL14 was found to be responsive to indole-acetic acid and 1-naphthylphthalamic acid, and RNA interference of OsSPL14 reduced polar auxin transport and increased 1-naphthylphthalamic acid sensitivity of rice plants. Further analysis revealed that OsSPL14 directly binds to the promoter of PIN-FORMED 1b (OsPIN1b) and PIN-LIKE6b (PILS6b) to regulate their expression positively. OsPIN1b and PILS6b were highly expressed in axillary buds and proved involved in bud outgrowth. Loss of function of OsPIN1b or PILS6b increased the tiller number of rice. Taken together, our findings suggested that OsSPL14 could control axillary bud outgrowth and tiller number by activating the expression of OsPIN1b and PILS6b to fine-tune auxin transport in rice.

Osa-miR1320 targets the ERF transcription factor OsERF096 to regulate cold tolerance via JA-mediated signaling

MicroRNAs play key roles in abiotic stress response. Rice (Oryza sativa L.) miR1320 is a species-specific miRNA that contributes to miR168-regulated immunity. However, it is still unknown whether miR1320 is involved in rice response to abiotic stress. In this study, we illustrated that the miR1320 precursor generated two mature miR1320s, miR1320-3p, and miR1320-5p, and they both displayed decreased expression under cold stress. Genetic evidence showed that miR1320 overexpression resulted in increased cold tolerance, while miR1320 knock down (KD) reduced cold tolerance. Furthermore, an APETALA2/ethylene-responsive factor (ERF) transcription factor OsERF096 was identified as a target of miR1320 via 5'-RACE and dual luciferase assays. OsERF096 expression was altered by miR1320 overexpression and KD and exhibited an opposite pattern to that of miR1320 in different tissues and under cold stress. Consistently, OsERF096 negatively regulated cold stress tolerance. Furthermore, we suggested that OsERF096 could bind to the GCC and DRE cis-elements and act as a transcriptional activator in the nucleus. Based on RNA-sequencing and targeted metabolomics assays, we found that OsERF096 modified hormone content and signaling pathways. Finally, phenotypic and reverse transcription-quantitative PCR assays showed that jasmonic acid (JA) methyl ester application recovered the cold-sensitive phenotype and JA-activated expression of three Dehydration Responsive Element Binding genes in the OsERF096-OE line. Taken together, our results strongly suggest that the miR1320-OsERF096 module regulates cold tolerance by repressing the JA-mediated cold signaling pathway.

Effect of Panicle Morphology on Grain Filling and Rice Yield: Genetic Control and Molecular Regulation

The demand for rice is likely to increase approximately 1.5 times by the year 2050. In contrast, the rice production is stagnant since the past decade as the ongoing rice breeding program is unable to increase the production further, primarily because of the problem in grain filling. Investigations have revealed several reasons for poor filling of the grains in the inferior spikelets of the compact panicle, which are otherwise genetically competent to develop into well-filled grains. Among these, the important reasons are 1) poor activities of the starch biosynthesizing enzymes, 2) high ethylene production leading to inhibition in expressions of the starch biosynthesizing enzymes, 3) insufficient division of the endosperm cells and endoreduplication of their nuclei, 4) low accumulation of cytokinins and indole-3-acetic acid (IAA) that promote grain filling, and 5) altered expressions of the miRNAs unfavorable for grain filling. At the genetic level, several genes/QTLs linked to the yield traits have been identified, but the information so far has not been put into perspective toward increasing the rice production. Keeping in view the genetic competency of the inferior spikelets to develop into well-filled grains and based on the findings from the recent research studies, improving grain filling in these spikelets seems plausible through the following biotechnological interventions: 1) spikelet-specific knockdown of the genes involved in ethylene synthesis and overexpression of β-CAS (β-cyanoalanine) for enhanced scavenging of CN⁻ formed as a byproduct of ethylene biosynthesis; 2) designing molecular means for increased accumulation of cytokinins, abscisic acid (ABA), and IAA in the caryopses; 3) manipulation of expression of the transcription factors like MYC and OsbZIP58 to drive the expression of the starch biosynthesizing enzymes; 4) spikelet-specific overexpression of the cyclins like CycB;1 and CycH;1 for promoting endosperm cell division; and 5) the targeted increase in accumulation of ABA in the straw during the grain filling stage for increased carbon resource remobilization to the grains. Identification of genes determining panicle compactness could also lead to an increase in rice yield through conversion of a compact-panicle into a lax/open one. These efforts have the ability to increase rice production by as much as 30%, which could be more than the set production target by the year 2050.

Genome-wide identification and characterization of the SPL gene family and its expression in the various developmental stages and stress conditions in foxtail millet (Setaria italica)

Background

Among the major transcription factors, SPL plays a crucial role in plant growth, development, and stress response. Foxtail millet (Setaria italica), as a C4 crop, is rich in nutrients and is beneficial to human health. However, research on the foxtail millet SPL (SQUAMOSA PROMOTER BINDING-LIKE) gene family is limited. 

Results

In this study, a total of 18 SPL genes were identified for the comprehensive analysis of the whole genome of foxtail millet. These SiSPL genes were divided into seven subfamilies (I, II, III, V, VI, VII, and VIII) according to the classification of the Arabidopsis thaliana SPL gene family. Structural analysis of the SiSPL genes showed that the number of introns in subfamilies I and II were much larger than others, and the promoter regions of SiSPL genes were rich in different cis-acting elements. Among the 18 SiSPL genes, nine genes had putative binding sites with foxtail millet miR156. No tandem duplication events were found between the SiSPL genes, but four pairs of segmental duplications were detected. The SiSPL genes expression were detected in different tissues, which was generally highly expressed in seeds development process, especially SiSPL6 and SiSPL16, which deserve further study. The results of the expression levels of SiSPL genes under eight types of abiotic stresses showed that many stress responsive genes, especially SiSPL9, SiSPL10, and SiSPL16, were highly expressed under multiple stresses, which deserves further attention.

Conclusions

In this research, 18 SPL genes were identified in foxtail millet, and their phylogenetic relationships, gene structural features, duplication events, gene expression and potential roles in foxtail millet development were studied. The findings provide a new perspective for the mining of the excellent SiSPL gene and the molecular breeding of foxtail millet.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08633-2.

Osa-miR535 targets SQUAMOSA promoter binding protein-like 4 to regulate blast disease resistance in rice

Many rice microRNAs have been identified as fine-tuning factors in the regulation of agronomic traits and immunity. Among them, Osa-miR535 targets SQUAMOSA promoter binding protein-like 14 (OsSPL14) to positively regulate tillers but negatively regulate yield and immunity. Here, we uncovered that Osa-miR535 targets another SPL gene, OsSPL4, to suppress rice immunity against Magnaporthe oryzae. Overexpression of Osa-miR535 significantly decreased the accumulation of the fusion protein SPL4_TBS -YFP that contains the target site of Osa-miR535 in OsSPL4. Consistently, Osa-miR535 mediated the cleavage of OsSPL4 mRNA between the 10th and 11th base pair of the predicted binding site at the 3' untranslated region. Transgenic rice lines overexpressing OsSPL4 (OXSPL4) displayed enhanced blast disease resistance accompanied by enhanced immune responses, including increased expression of defense-relative genes and up-accumulated H₂ O₂ . By contrast, the knockout mutant osspl4 exhibited susceptibility. Moreover, OsSPL4 binds to the promoter of GH3.2, an indole-3-acetic acid-amido synthetase, and promotes its expression. Together, these data indicate that Os-miR535 targets OsSPL4 and OsSPL4-GH3.2, which may parallel the OsSPL14-WRKY45 module in rice blast disease resistance.

OsSPLs Regulate Male Fertility in Response to Different Temperatures by Flavonoid Biosynthesis and Tapetum PCD in PTGMS Rice

Photoperiod and thermo-sensitive genic male sterile (PTGMS) rice is an important resource for two line hybrid rice production. The SQUAMOSA–promoter binding, such as the (SPL) gene family, encode the plant specific transcription factors that regulate development and defense responses in plants. However, the reports about SPLs participating in male fertility regulation are limited. Here, we identified 19 OsSPL family members and investigated their involvement in the fertility regulation of the PTGMS rice lines, PA2364S and PA2864S, with different fertility transition temperatures. The results demonstrated that OsSPL2, OsSPL4, OsSPL16 and OsSPL17 affect male fertility in response to temperature changes through the MiR156-SPL module. WGCNA (weighted gene co-expression network analysis) revealed that CHI and APX1 were co-expressed with OsSPL17. Targeted metabolite and flavonoid biosynthetic gene expression analysis revealed that OsSPL17 regulates the expression of flavonoid biosynthesis genes CHI, and the up regulation of flavanones (eriodictvol and naringenin) and flavones (apigenin and luteolin) content contributed to plant fertility. Meanwhile, OsSPL17 negatively regulates APX1 to affect APX (ascorbate peroxidase) activity, thereby regulating ROS (reactive oxygen species) content in the tapetum, controlling the PCD (programmed cell death) process and regulating male fertility in rice. Overall, this report highlights the potential role of OsSPL for the regulation of male fertility in rice and provides a new insight for the further understanding of fertility molecular mechanisms in PTGMS rice.

Agrobacterium-Mediated High-Efficiency Genetic Transformation and Genome Editing of Chaling Common Wild Rice (Oryza rufipogon Griff.) Using Scutellum Tissue of Embryos in Mature Seeds

Genetic transformation is an important strategy for revealing gene function, and it is used extensively in both functional genomics study and molecular breeding of rice. Demand for its application in wild Oryza species is rising for their extensive genetic diversity. However, genetic transformation of wild Oryza accessions with AA genome using calli induced from scutellum tissue of embryos in mature seeds has not been successfully established. In the present study, we used Chaling common wild rice (CLCWR) (Oryza rufipogon Griff.) with AA genome to successfully establish an Agrobacterium-mediated genetic transformation system based on scutellum tissue of embryos in mature seeds. The calli from embryos in mature seeds of CLCWR were easy to be induced and regenerated. The callus induction rate and texture were optimum under 2.5 mg/L 2,4-D. The optimal hormone combination used for regeneration was 2 mg/L ZT + 0.1 mg/L NAA. Studies on genetic transformation and genome editing showed that the transformation efficiency was 87-94%, the efficiency of single genome editing and multiplex genome editing were about 60-70% and 20-40%, respectively. Compared with Nipponbare (Nip), CLCWR had higher Hygromycin-resistant callus frequency and transformation efficiency. Taken together, our study establishes a highly efficient transformation system for common wild rice with AA genome and provides a good rice material for de novo domestication by genome editing in the future.

Genes and Their Molecular Functions Determining Seed Structure, Components, and Quality of Rice

With the improvement of people's living standards and rice trade worldwide, the demand for high-quality rice is increasing. Therefore, breeding high quality rice is critical to meet the market demand. However, progress in improving rice grain quality lags far behind that of rice yield. This might be because of the complexity of rice grain quality research, and the lack of consensus definition and evaluation standards for high quality rice. In general, the main components of rice grain quality are milling quality (MQ), appearance quality (AQ), eating and cooking quality (ECQ), and nutritional quality (NQ). Importantly, all these quality traits are determined directly or indirectly by the structure and composition of the rice seeds. Structurally, rice seeds mainly comprise the spikelet hull, seed coat, aleurone layer, embryo, and endosperm. Among them, the size of spikelet hull is the key determinant of rice grain size, which usually affects rice AQ, MQ, and ECQ. The endosperm, mainly composed of starch and protein, is the major edible part of the rice seed. Therefore, the content, constitution, and physicochemical properties of starch and protein are crucial for multiple rice grain quality traits. Moreover, the other substances, such as lipids, minerals, vitamins, and phytochemicals, included in different parts of the rice seed, also contribute significantly to rice grain quality, especially the NQ. Rice seed growth and development are precisely controlled by many genes; therefore, cloning and dissecting these quality-related genes will enhance our knowledge of rice grain quality and will assist with the breeding of high quality rice. This review focuses on summarizing the recent progress on cloning key genes and their functions in regulating rice seed structure and composition, and their corresponding contributions to rice grain quality. This information will facilitate and advance future high quality rice breeding programs.

Micro-RNA-Regulated SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) Gene Expression and Cytokinin Accumulation Distinguish Early-Developing Male and Female Inflorescences in Oil Palm (Elaeis guineensis)

Sexual differentiation of inflorescences and flowers is important for reproduction and affects crop plant productivity. We report here on a molecular study of the process of sexual differentiation in the immature inflorescence of oil palm (Elaeis guineensis). This species is monoecious and exhibits gender diphasy, producing male and female inflorescences separately on the same plant in alternation. Three main approaches were used: small RNA-seq to characterise and study the expression of miRNA genes; RNA-seq to monitor mRNA accumulation patterns; hormone quantification to assess the role of cytokinins and auxins in inflorescence differentiation. Our study allowed the characterisation of 30 previously unreported palm MIRNA genes. In differential gene and miRNA expression studies, we identified a number of key developmental genes and miRNA-mRNA target modules previously described in relation to their developmental regulatory role in the cereal panicle, notably the miR156/529/535-SQUAMOSA PROMOTER-BINDING PROTEIN-LIKE (SPL) gene regulatory module. Gene enrichment analysis highlighted the importance of hormone-related genes, and this observation was corroborated by the detection of much higher levels of cytokinins in the female inflorescence. Our data illustrate the importance of branching regulation within the developmental window studied, during which the female inflorescence, unlike its male counterpart, produces flower clusters on new successive axes by sympodial growth.

Conservation and Divergence of SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) Gene Family between Wheat and Rice

The SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) gene family affects plant architecture, panicle structure, and grain development, representing key genes for crop improvements. The objective of the present study is to utilize the well characterized SPLs’ functions in rice to facilitate the functional genomics of TaSPL genes. To achieve these goals, we combined several approaches, including genome-wide analysis of TaSPLs, comparative genomic analysis, expression profiling, and functional study of TaSPL3 in rice. We established the orthologous relationships of 56 TaSPL genes with the corresponding OsSPLs, laying a foundation for the comparison of known SPL functions between wheat and rice. Some TaSPLs exhibited different spatial–temporal expression patterns when compared to their rice orthologs, thus implicating functional divergence. TaSPL2/6/8/10 were identified to respond to different abiotic stresses through the combination of RNA-seq and qPCR expression analysis. Additionally, ectopic expression of TaSPL3 in rice promotes heading dates, affects leaf and stem development, and leads to smaller panicles and decreased yields per panicle. In conclusion, our work provides useful information toward cataloging of the functions of TaSPLs, emphasized the conservation and divergence between TaSPLs and OsSPLs, and identified the important SPL genes for wheat improvement.

Molecular identification and functional verification of SPL9 and SPL15 of Lilium

The transformation of plants from juveniles to adults is a key process in plant growth and development, and the main regulatory factors are miR156 and SQUAMOSA promoter binding protein-like (SPL) transcription factors. Lilium is an ornamental bulb, but it has a long maturation time. In this experiment, Lilium bulbs were subjected to a temperature treatment of 15 °C for 4 weeks to initiate vegetative phase change. Transmission electron microscopy indicated the cell wall of bud core tissue undergoing vegetative phase change became thinner, the starch grains were reduced, and the growth of the juvenile stage was accelerated. The key transcription factors LbrSPL9 and LbrSPL15 were cloned, and the phylogenetic analysis showed they possessed high homology with other plant SPLs. Subcellular localization and transcription activation experiments confirmed LbrSPL9 and LbrSPL15 were mainly located in the nucleus and exhibited transcriptional activity. The results of in situ hybridization showed the expression levels of LbrSPL9 and LbrSPL15 were increased after temperature change treatment. The functional verification experiment of the transgenic plants confirmed that the overexpression of LbrSPL9 and LbrSPL15 could shorten maturation time. These findings help elucidate the regulatory mechanisms of phase transition in Lilium and provide a reference for breeding research in other bulbous flowers.

The regulation of DKGA2ox1 and miR171f_3 in scion dwarfing with Diospyros kaki Thunb. cv. 'Nan-tong-xiao-fang-shi' as interstocks

High-throughput sequencing and yeast one and two-hybrid library screening reveal that DKGA2ox1 and miR171f_3 are involved in the regulation of scion dwarfing with 'Nan-tong-xiao-fang-shi' as interstocks. Diospyros kaki Thunb. cv. Nan-tong-xiao-fang-shi ('Nan-tong-xiao-fang-shi') interstocks play a critical role in the scion dwarfing. However, the understanding of the molecular signaling pathways that regulate the scion dwarfing with 'Nan-tong-xiao-fang-shi' as interstocks remain unclear. In this work, the regulatory network in the scion dwarfing with 'Nan-tong-xiao-fang-shi' as interstocks was identified. Using a yeast one-hybrid library screening, luciferase activity analysis, luciferase complementation imaging assays and GFP signal detection, a SPL transcription factor named Diospyros kaki SPL (DKSPL), potentially functioning as a transcriptional activator of the Diospyros kaki GA2ox1 (DKGA2ox1) gene, was identified as a key stimulating factor in the persimmon growth and development. The DKSPL was found in the nucleus, and might play a role in the transcriptional regulation system. A microRNA named miR171f_3 was identified, which might act as a negative regulator of Diospyros kaki SCR (DKSCR) in persimmon. The interactions between DKSCR and seven proteins were experimentally validated with a yeast two-hybrid library screening. Compared to the non-grafted wildtype persimmon, the tissue section of graft union healed well due to the increased expression of cinnamyl-alcohol dehydrogenase. These results indicate that DKGA2ox1 and miR171f_3 may co-promote the scion dwarfing by plant hormone signal transduction pathways.

Genome-Wide Characterization and Expression Analysis of the SBP-Box Gene Family in Sweet Orange (Citrus sinensis)

SBP-box is an important plant-specific transcription factor family and is involved in diverse biological processes. Here, we identified a total of 15 SBP-BOX genes in the important fruit crop sweet orange (Citrus sinensis) and characterized their gene structures, conserved domain and motif, chromosomal location, and cis-acting regulatory elements. SBP genes were classified into four subfamilies based on the amino acid sequence homology, and the classification is equally strongly supported by the gene and protein structures. Our analysis revealed that segmental duplication events were the main driving force in the evolution of CsSBP genes, and gene pairs might undergo extensive purifying selection. Further synteny analysis of the SBP members among sweet orange and other plant species provides valuable information for clarifying the CsSBP family evolutionary relationship. According to publicly available RNA-seq data and qRT-PCR analysis from various sweet orange tissues, CsSBP genes may be expressed in different tissues and developmental stages. Gene expression analysis showed variable expression profiles of CsSBP genes under various abiotic stresses, such as high and low-temperature, salt, and wound treatments, demonstrating the potential role of SBP members in sweet orange response to abiotic stress. Noticeably, all CsSBP genes were also downregulated in sweet orange upon the infection of an important fungal pathogen Diaporthe citri. Our results provide valuable information for exploring the role of SBP-Box in sweet orange.

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.

Photocontrol of Axillary Bud Outgrowth by MicroRNAs: Current State-of-the-Art and Novel Perspectives Gained From the Rosebush Model

Shoot branching is highly dependent on environmental factors. While many species show some light dependence for branching, the rosebush shows a strict requirement for light to allow branching, making this species an excellent model to further understand how light impinges on branching. Here, in the first part, we provide a review of the current understanding of how light may modulate the complex regulatory network of endogenous factors like hormones (SL, IAA, CK, GA, and ABA), nutrients (sugar and nitrogen), and ROS to control branching. We review the regulatory contribution of microRNAs (miRNAs) to branching in different species, highlighting the action of such evolutionarily conserved factors. We underline some possible pathways by which light may modulate miRNA-dependent regulation of branching. In the second part, we exploit the strict light dependence of rosebush for branching to identify putative miRNAs that could contribute to the photocontrol of branching. For this, we first performed a profiling of the miRNAs expressed in early light-induced rosebush buds and next tested whether they were predicted to target recognized regulators of branching. Thus, we identified seven miRNAs (miR156, miR159, miR164, miR166, miR399, miR477, and miR8175) that could target nine genes (CKX1/6, EXPA3, MAX4, CYCD3;1, SUSY, 6PFK, APX1, and RBOHB1). Because these genes are affecting branching through different hormonal or metabolic pathways and because expression of some of these genes is photoregulated, our bioinformatic analysis suggests that miRNAs may trigger a rearrangement of the regulatory network to modulate branching in response to light environment.