Genome‑wide identification and characterization of miR396 family members and their target genes GRF in sorghum (Sorghum bicolor (L.) moench)

Abstract

MicroRNAs (miRNAs) widely participate in plant growth and development. The miR396 family, one of the most conserved miRNA families, remains poorly understood in sorghum. To reveal the evolution and expression pattern of Sbi-miR396 gene family in sorghum, bioinformatics analysis and target gene prediction were performed on the sequences of the Sbi-miR396 gene family members. The results showed that five Sbi-miR396 members, located on chromosomes 4, 6, and 10, were identified at the whole-genome level. The secondary structure analysis showed that the precursor sequences of all five Sbi-miR396 potentially form a stable secondary stem–loop structure, and the mature miRNA sequences were generated on the 5′ arm of the precursors. Sequence analysis identified the mature sequences of the five sbi-miR396 genes were high identity, with differences only at the 1st, 9th and 21st bases at the 5’ end. Phylogenetic analysis revealed that Sbi-miR396a, Sbi-miR396b, and Sbi-miR396c were clustered into Group I, and Sbi-miR396d and Sbi-miR396e were clustered into Group II, and all five sbi-miR396 genes were closely related to those of maize and foxtail millet. Expression analysis of different tissue found that Sbi-miR396d/e and Sbi-miR396a/b/c were preferentially and barely expressed, respectively, in leaves, flowers, and panicles. Target gene prediction indicates that the growth-regulating factor family members (SbiGRF1/2/3/4/5/6/7/8/10) were target genes of Sbi-miR396d/e. Thus, Sbi-miR396d/e may affect the growth and development of sorghum by targeting SbiGRFs. In addition, expression analysis of different tissues and developmental stages found that all Sbi-miR396 target genes, SbiGRFs, were barely expressed in leaves, root and shoot, but were predominantly expressed in inflorescence and seed development stage, especially SbiGRF1/5/8. Therefore, inhibition the expression of sbi-miR396d/e may increase the expression of SbiGRF1/5/8, thereby affecting floral organ and seed development in sorghum. These findings provide the basis for studying the expression of the Sbi-mir396 family members and the function of their target genes.

Introduction

Sorghum (Sorghum bicolor [L.] Moench) is one of the world’s top five grain crops, along with maize, rice, wheat, and barley. It is a staple food for more than 500 million people worldwide [1–3]. Owing to its excellent adaptability to heat and drought stress, sorghum is a key crop for safeguarding food security in response to climate change [4]. Recently, sorghum has also been widely used for different purposes, for example, as feed and forage, in brewing, and in biomass energy production, which has a significant effect on global agricultural production [5–8]. In recent years, with the continuous development of animal husbandry and brewery industries, sorghum production has struggled to meet the market demand [9, 10]. Crop yield is influenced by various factors, of which plant architecture has a significant effect [11, 12]. An ideal plant architecture is beneficial for a crop plant to develop an optimal shape for light exposure, which enables it to fully capture and use solar energy throughout the growth period and consequently maximizes the economic yield [11]. Plant architecture studies mainly focus on tillering and stem, leaf, panicle types, which are the traits determining crop yield [13, 14].

MicroRNAs (miRNAs) are a class of endogenous non-coding small RNAs commonly found in eukaryotic organisms, generally consisting of 20–22 nucleotides [15]. In plants, miRNA bind to the complementary sites within target mRNA to repress gene expression at the posttranscriptional level, thereby regulating plant growth, development, and response to adverse stress [16–18]. The miR396 family, a conserved miRNA family in plants, has been identified in crops such as Arabidopsis, rice, wheat, tobacco, and soybean [19–23]. Most, but not all, growth-regulating factor genes (GRFs) are the target genes of miR396 [24, 25]. miR396 negatively regulates the expression of GRFs to participate in the regulation of leaf and floral organ morphogenesis, root and stem tissue elongation, and grain size and number, ultimately influencing plant architecture and crop yield. In Arabidopsis, overexpression of miR396 decreases the mRNA levels of target GRFs, resulting in phenotypes such as short roots, small leaves, reduced number of petals and stamens, pistils with a single carpel, small siliques, and abnormal growing points at the tip of roots [25–27]. In rice, either plants overexpressing miR396 or mutants carrying an insertionally inactivated GRF gene exhibit morphological changes such as retarded growth, reduced plant height, shortened roots, and narrowed leaves, whereas inhibition of miR396 expression alters the panicle structure [28–30]. Furthermore, rice miR396 (OsmiR396d) diminishes both gibberellin biosynthesis and signaling by inhibiting OsGRF6 expression, which reduces plant height; it also relieves the inhibition of brassinolide signaling through suppression of OsGRF4, thereby positively regulating leaf angle [31]. The OsmiR396–OsGRF3/4/8 module effectively regulates grain size and yield [14, 21, 32–34], and the OsmiR396–OsGRF6/10 module regulates floral organ development. For example, rice plants overexpressing OsmiR396 or osgrf6/10 double mutant exhibit open husks, long sterile lemmas, and altered floral organ morphology, whereas rice plants with suppressed expression of OsmiR396 or overexpression of OsGRF6 exhibit a considerable increase in the number of secondary branches per panicle and glumes [35].

In addition to their essential roles in plant growth and development, miR396 and its target GRF genes also enhance plant responses required to tolerate adverse stress conditions. Overexpression of miR396 increases drought tolerance in crop plants. A previous study reports that miR396a/b-overexpressing Arabidopsis plants displayed a narrow leaf phenotype, with lower stomatal density and decreased transpiration rate; their survival rate was remarkably higher than that of empty vector controls after 14 d of drought treatment [25]. In tobacco lines overexpressing miR396a-5p, expression of NtGRFs is inhibited, osmoregulation and drought tolerance are enhanced, and reactive oxygen species levels are lowered [36]. miR396 mutation may increase crop tolerance to saline–alkali and nitrogen deficiency stress [37–40]. Under nitrogen-deficient conditions, osmir396ef double mutants exhibit higher grain yield and aboveground biomass. Recent research has also found that miR396–GRF markedly increases the regeneration efficiency of explants [41, 42].

The miR396-mediated GRF regulatory module (miR396–GRF) has shown potential application value in improving plant type, crop yield, and stress resistance and increasing the efficiency of plant genetic transformation. However, to date, members of the sorghum miR396 (Sbi-miR396) family have not been identified and functionally analyzed. In the present study, bioinformatics methods were used to analyze the chromosomal distribution, phylogenetic relationships, base conservation, secondary structure, and target gene expression of the Sbi-miR396 family. The findings provide molecular evidence for the evolution and biological functions of the miR396 family in sorghum.

Materials and methods

Genome-wide identification and sequence analysis of Sbi-miR396 and its target genes

The precursor sequences, mature sequences, and chromosomal localization information of the Sbi-miR396 family members were downloaded from the miRBase database (https://www.mirbase.org/ftp.shtml). Base conservation of these sequences were analyzed using WebLogo 3 (http://weblogo.threeplusone.com/create.cgi) [40]. The sorghum whole-genome gff file (Sorghum_bicolor_NCBIv3.54.gff3) was downloaded from EnsemblPlants (http://plants.ensembl.org/Sorghum_bicolor/Info/Index). Physical position of each Sbi-miR396 family gene was extracted from this file, and chromosomal localization map was constructed using Gene Location Visualize from GTF/GFF in TBtools software [43]. The secondary stem–loop structure of the Sbi-miR396 precursors was predicted using RNAfold (http://rna.tbi.univie.ac.at./cgi-bin/RNAfold.cgi). A folding algorithm with basic options was employed to select the minimum free energy and partition function.

The target genes of Sbi-miR396 were predicted using psRNATarget (http://plantgrn.noble.org/psRNATarget/) with default parameter values. Genes with the maximum expectation value ≤1.5 were selected as target genes, and the sites of interaction between Sbi-miR396 and candidate target genes were analyzed. In addition, we estimated the molecular weight (MW) and isoelectric point (pI) of the candidate target proteins using ExPASy (https://web.expasy.org/protparam/), the functional annotation was performed using the National Center for Biotechnology Information (NCBI) database (http://www.ncbi.nlm.nih.gov/), and multiple alignment of the amino acid sequences were performed using DNAMAN 6.0 software.

Phylogenetic analysis of Sbi-miR396 and its target gene SbiGRFs

The miR396 precursor sequences of maize (Zea mays L.), foxtail millet (Setaria italica [L.] P. Beauv), and rice (Oryza sativa L.) were downloaded from the miRBase database. The GRF protein sequences of maize, foxtail millet, rice, and Arabidopsis were downloaded from PlantTFDB database (http://planttfdb.gao-lab.org/index.php?sp=Sbi). Multiple alignment of miR396 precursor sequences and GRF protein sequences were performed using ClustalW with default parameters. Phylogenetic trees were then constructed using the neighbor-joining method of MEGA v7.0 software with 1000 bootstraps and other default parameters.

Sequence and structural analysis of SbiGRF gene family targeted by Sbi-miR396

Conserved motifs of SbiGRFs were analyzed using the multiple expectation maximization for motif elicitation (MEME, http://meme-suite.org/). Gene structure information was extracted from the sorghum whole-genome gff file. These identified motifs and gene structures were further visualized by TBtools.

Expression analysis of Sbi-miR396 and its target gene SbiGRFs

Published RNA-seq data from different tissues or developmental stages of sorghum were downloaded from the PmiREN database (https://www.pmiren.com/browsehref4?wzid=209&mbt=5&type=sRNA-Seq#zqnav) and the Expression Atlas database (https://www.ebi.ac.uk/gxa/home) for the expression analysis of Sbi-miR396 and its target gene SbiGRFs, respectively. The PmiREN database contained sorghum leaves, flowers, and panicles. The Expression Atlas database included different tissues and developmental stages, roots, shoots, leaves at seedling developmental stage; pistils and pollen at booting stage; inflorescences during the developmental process; and seeds during the developmental process. The expression level of Sbi-miR396 and SbiGRFs were estimated using the fragments per kilobase per million reads value (FPKM), and heatmaps were displayed by Tbtools.

RNA isolation and qRT-PCR analysis

Three materials were used for qRT-PCR analysis, namely, 10480A, L17R and Jinnuo 3. Total RNA was extracted from young spikes of the inflorescence developmental stage and early development seeds using the RNAprep Pure Plant Kit (TIANGEN, Beijing, China) following the manufacturer’s instructions. The purity and content of total RNA were detected by 1.0% agarose gel electrophoresis and NanoDrop 2000 (Thermo Fisher Scientific, USA). First-strand cDNA of miRNA and mRNA were synthesized with the Mir-X miRNA First-Strand Synthesis Kit and the PrimeScript^TM RT reagent Kit with gDNA Eraser (Perfect Real Time) (Takara, Beijing, China), respectively. Quantitative real-time PCR (qRT-PCR) of Sbi-miR396 and SbiGRFs were performed with three technical and three biological replicates in the LightCyclerTM 96 (Roche, Basel) using the Mir-X miRNA qRT-PCR TB Green Kit and the TB Green Premix Ex Taq^TM II (Tli RNaseH Plus) (Takara, Beijing, China), respectively. The U6 and SbiGAPDH gene were used as internal reference. The qRT-PCR reaction system of Sbi-mi396 contained 10 μL 2×TB Green Advantage Premix, 0.4 μL (10 μM) RT-primer, 0.4 μL (10 μM) 5’ primer (S1 Table), 2 μL cDNA template, 0.4 μL 50× Rox Reference Dye II, and 6.8 μL ddH₂O. The reaction procedure was as follows: denaturation at 95°C for 10 s, followed by 40 cycles at 95°C for 5 s and 60°C for 20 s. The qRT-PCR reaction system of SbiGRFs contained 10 μL 2× TB Green Premix Ex Taq, 0.8 μL (10 μM) forward primer, 0.8 μL (10 μM) reverse primer (S1 Table), 2 μL cDNA template, 0.4 μL 50× Rox Reference Dye II, and 6.0 μL ddH₂O. The reaction procedure was as follows: denaturation at 95°C for 3 min, followed by 40 cycles at 95°C for 10 s and 60°C for 30 s. After the reaction was completed, relative expression levels were calculated using the△CT method. Statistical significance was calculated using SPSS 19.0 software (SPSS Corp., Chicago, IL).

Results

Sequences and chromosomal localization of the Sbi-miR396 family members

According to the miRBase database, the miR396 family comprises five members in the sorghum genome, designated as Sbi-miR396a, b, c, d, and e (Table 1). The precursor sequences of the five Sbi-miR396 members consisted of 88–186 bases. Bases at positions 1–22 at the 5′ end were highly conserved in the precursor sequences, whereas bases at other positions were less conserved (Fig 1A).

[Figure omitted. See PDF.]

(A) Precursor sequence; (B) Mature sequence; (C) Star sequence.

[Figure omitted. See PDF.]

Mature Sbi-miR396s were all produced from the highly conserved bases at positions 1–22 of the precursor sequences. The mature sequences of Sbi-miR396a, b, and c each contained 21 bases, and Sbi-miR396d and e each contained 22 bases. All 21 bases were highly conserved between Sbi-miR396a and Sbi-miR396b. Excluding a dissimilar base at position 1 at the 5′ end, 21 bases were highly conserved between Sbi-miR396d and Sbi-miR396e. Excluding a dissimilar base at position 1 at the 3′ end, Sbi-miR396c were highly conserved with Sbi-miR396a and Sbi-miR396b for 20 bases. Sbi-miR396d and Sbi-miR396e contained an additional base, G, at position 9 from the 5′ end, compared with Sbi-miR396a and Sbi-miR396b (Fig 1B, Table 1). The star sequences of all five Sbi-miR396 members comprised 21 bases. Of these, bases at positions 1, 3–6, 8–11, and 17–20 from the 5′ end were highly conserved; bases at positions 2, 12, 15, and 21 were moderately conserved; and bases at the remaining positions were poorly conserved (Fig 1C, Table 1).

The chromosomal localization results indicated that Sbi-miR396a, Sbi-miR396c, and Sbi-miR396d were all located on chromosome 4, whereas Sbi-miR396b and Sbi-miR396e were located on chromosomes 10 and 6, respectively (Fig 2, Table 1).

[Figure omitted. See PDF.]

Secondary structure of the Sbi-miR396 precursor sequences

The precursor sequences of the five Sbi-miR396 family members differed in length; Sbi-miR396e was the longest at 186 bp, and Sbi-miR396a was the shortest at 88 bp. The secondary structure prediction results showed that a stable secondary stem–loop structure could be formed in all precursor sequences. All mature miRNAs were located on the 5′ arm of the corresponding precursors and were highly conserved (Fig 3).

[Figure omitted. See PDF.]

Phylogenetic evolution of the Sbi-miR396 family members

Phylogenetic analysis of the precursor sequences revealed that the miR396 family members of sorghum, maize, foxtail millet, and rice were divided into three major groups (Groups I–III). Sbi-miR396a, Sbi-miR396b, and Sbi-miR396c clustered in Group I, whereas Sbi-miR396d and Sbi-miR396e clustered in Group II. Specifically, the miR396 family members in sorghum were closely related to those in maize and foxtail millet, exhibiting clusters such as Sbi-miR396a/Zma-miR396g/Sit-miR396c/Zma-miR396a and Sbi-miR396e/Sti-miR396b/Zma-miR396c, but distantly related to those in rice (Fig 4).

[Figure omitted. See PDF.]

Gene prefixtion Sbi, Zma, Sit, and Osa represent Sorghum bicolor, Zea mays, Setaria italica, and Oryza sativa, respectively.

Tissue-specific expression of Sbi-miR396 family members

Tissue-specific expression analysis revealed that Sbi-miR396a and Sbi-miR396b were barely expressed in young leaves, flowers, and young panicles of sorghum. Sbi-miR396c was slightly expressed in young leaves but not expressed in flowers and young panicles. In contrast, Sbi-miR396d and Sbi-miR396e were preferentially expressed in young leaves, flowers, and young panicles; their highest expression levels were observed in young leaves, followed by those in the flowers (Fig 5).

[Figure omitted. See PDF.]

Prediction and physicochemical features of target genes of Sbi-miR396

Taking into account the minimal expression of Sbi-miR396a/b/c and preferential expression of Sbi-miR396d/e in various tissues of sorghum, subsequent analyses were performed only for Sbi-miR396d and Sbi-miR396e. Table 2 lists the predicted target genes of Sbi-miR396. Sbi-miR396d and Sbi-miR396e were found to target the same nine SbiGRF family members, SbiGRF1, 2, 3, 4, 5, 6, 8, 9, and 10. The target sites of the SbiGRFs were all located in the third or second exon. The protein sequences of the nine SbiGFRs ranged from 271 to 601 amino acids in length, with a theoretical isoelectric point of 4.73–10.03 and molecular weight of 28.31–62.48 kDa.

[Figure omitted. See PDF.]

Phylogeny and classification of SbiGRFs

The GRF families of maize, foxtail millet, rice, and Arabidopsis had 15, 10, 12, and 9 members, respectively. The protein sequences of these GRFs and the nine GRFs of sorghum were used to construct a phylogenetic tree to analyze the evolutionary relationships among the GRF families. The 55 GRFs clustered into five groups (Group I–V; Fig 6). Sorghum, maize, foxtail millet, and rice GRFs were distributed in Group I, II, IV, and V. Arabidopsis GRFs were distributed in Group I, II, III, and IV. All GRF members in Group V were derived from dicotyledons. The results suggest that Group I, II, and IV already existed prior to differentiation of monocotyledonous and dicotyledonous plants, whereas Group V emerged after differentiation of monocotyledonous and dicotyledonous plants. The phylogenetic tree suggested that the GRF family of sorghum was most closely related to that of maize and foxtail millet, followed by rice, and most distantly related to that of the dicotyledon, Arabidopsis.

[Figure omitted. See PDF.]

Gene prefixtion Sbi, Zma, Sit, Osa, and At represent Sorghum bicolor, Zea mays, Setaria italica, Oryza sativa, and Arabidopsis thaliana, respectively.

Gene structure of SbiGRFs

The conserved protein sequences of the GRF family genes, putative targets of Sbi-miR396, were analyzed using MEME, which identified a total of five motifs. The N terminus of nine SbiGFR proteins contained intact conserved domains, WRC (Motif1) and QLQ (Motif2); in all cases, the QLQ domain was located in front of the WRC domain (Fig 7). In addition to the conserved amino acids Q/Gln-L/Leu-Q/GLn, the QLQ domain also comprised the hydrophobic, acidic amino acid residues F/Phe, Y/Tyr, L/Leu, E/Glu, and P/Pro. The WRC domain contained the conserved amino acid residues W/Trp-R/Arg-C/Cys and the C3H motif consisting of three C/Cys and one H/His, immediately followed by the conserved RSRK-VE region (Fig 8). Furthermore, SbiGRF3–5 contained Motif3, 4, and 5; SbiGRF1 and 2 contained Motif3 and 4; SbiGRF6 and 8 contained Motif4; and SbiGRF10 contained Motif3.

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

(A) WRC domain; (B) QLQ domain.

Gene structure analysis revealed that the nine SbiGRF genes differed in the number of exons. SbiGRF2 and 10 contained two exons, SbiGRF1 and 5 contained three exons each, SbiGRF4 contained five exons, whereas the remaining four SbiGRF genes contained four exons each (Fig 7).

Expression analysis of SbiGRFs in different tissues and developmental stages

miR396 affects plant growth and development by regulating the expression of its target genes SbiGRFs. We used online available transcriptome data of three tissues (leaf, shoot and root) and eight developmental stages of inflorescence and seed (pistil and pollen at booting stage, Early inflorescence, Emerging inflorescence, Inflorescence size 1 to 5 millimeter, Inflorescence size 1 to 10 millimeter, Inflorescence size 1 to 2 centimeter, seed 5 days after pollination, seed 10 days after pollination) to investigate the expression patterns of SbiGRFs genes (Fig 9). It was found that SbiGRFs were predominantly expressed in inflorescence and seed. The expression patterns of SbiGRFs exhibited remarkable tissue specificity, indicating that they are mainly involved in regulating the growth and development of floral organs and seeds in sorghum. Specifically, at the seedling developmental stage, the expression levels of SbiGRFs were considerably low across in all tissues, especially in seedling leaves. At booting stage, SbiGRFs were barely expressed in pollen but preferentially expressed in pistils. During inflorescence development stage, the expression levels of SbiGRFs were higher after than before inflorescence emergence. Before inflorescence emergence (I1), SbiGRF6/8 were expressed at slightly higher levels, whereas after inflorescence emergence(I3, I4 and I5), SbiGRF5/8 were expressed at relatively higher levels, followed by SbiGRF1/3/4/6, and SbiGRF2/10 showed low expression levels. During seed development stage, the expression levels of SbiGRFs in early developing seeds (P5S) were higher than those in intermediate developing seeds (P10S). In addition, SbiGRF1/3/5/6/8 showed higher expression levels than other SbiGRFs.

[Figure omitted. See PDF.]

SR, SL, and SS represent root, leaf, and shoot at seedling developmental stage, respectively; Bpi and Bpo represent pistil and pollen at booting stage; I1, Early inflorescence; I2, Emerging inflorescence; I3, Inflorescence size 1 to 5 millimeter; I4, Inflorescence size 1 to 10 millimeter; I5, Inflorescence size 1 to 2 centimeter; P5S and P10S represent seed at 5 and 10 days after pollination.

To further clarify the role of Sbi-miR396 and their target gene SbiGRFs in inflorescence and seed development, qRT-PCR was performed from young spikes of the inflorescence developmental stage and early development seeds in three sorghum varieties, and the results are shown in Fig 10. Consistent with the RNA-seq results, Sbi-miR396a, Sbi-miR396b, and Sbi-miR396c were expressed at low levels or not expressed in the young spikes of the inflorescence developmental stage and early development seeds, while Sbi-miR396d and Sbi-miR396e were preferentially expressed (Fig 10A). Moreover, Sbi-miR396d was expressed slightly higher than Sbi-miR396e at both developmental stages. Nine Sbi-miR396 target genes, SbiGRFs, were expressed at both developmental stages, and most SbiGRFs were expressed slightly higher in the inflorescence development stage than in the seed development stage (Fig 10B). In sorghum varieties, SbiGRF8 exhibited the highest expression level, followed by SbiGRF1 and SbiGRF5. SbiGRF9 and SbiGRF4 showed the lowest expression level. SbiGRF2, SbiGRF3, SbiGRF6, and SbiGRF10 showed intermediate expression level. Furthermore, SbiGRF8 and SbiGRF1 were expressed at about two times that of SbiGRF5 and about ten times that of SbiGRF9 and SbiGRF4. These findings suggest that Sbi-miR396d/e and SbiGRF8/1/5 play a major role in inflorescence and seed development stage in sorghum.

[Figure omitted. See PDF.]

Expression analysis of Sbi-miR396 (A) and SbiGRFs (B) in inflorescence and seed development. 1, young spikes of the inflorescence developmental stage; 2, early development seeds.

Discussion

Sorghum exhibits resistance to environmental stresses such as drought stress, flood stress, saline-alkali stress, nutrient deficiency, and heat stress. Consequently, it has comparative advantages when grown in fields with poor ecological and site conditions. However, compared with major crops such as rice and wheat or model crops such as Arabidopsis, sorghum has received less research attention regarding its molecular biology, especially miRNAs. To date, the miRBase database includes miR396s present in 52 species, with 221 mature sequences and 170 precursor sequences. Five plant species exhibit a high number of mature miR396 sequences, namely soybean (n = 15), rice (n = 13), maize (n = 13), Brachypodium distachyon (n = 10), and Aegilops tauschii (n = 10). Other plants such as wheat, sugarcane, and papaya have only one mature miR396 sequence each. These findings indicate that the number of mature miR396 sequences does not correspond exactly to the number of precursor sequences and some precursor sequences can produce miR396 at both ends of the hairpin. Different mature sequences may originate from the same precursor sequence; in contrast, the same mature sequence may be derived from various precursor sequences. In the present study, it was found that the Sbi-miR396 family comprises five mature and five precursor sequences in sorghum. Each mature Sbi-miR396 is derived from the highly conserved positions 1–22 at the 5′ end of the corresponding precursor sequence. The five mature Sbi-miR396 sequences are highly similar to with each other; except for three base differences between Sbi-miR396c and Sbi-miR396d, only one or no base difference was observed between any two Sbi-miR396 sequences.

The miR396 family shows high evolutionary conservation across species [44]. The present findings indicate that the miR396 precursor sequences of sorghum have a close evolutionary relationship with those of foxtail millet and maize, suggesting that the evolutionary relationship of the miR396 gene families is consistent with the phylogenetic relationship among the species themselves. The Arabidopsis miR396 family includes four members, Ath-miR396a-5p, Ath-miR396a-3p, Ath-miR396b-5p, and Ath-miR396b-3p. The sequence alignment analysis demonstrated that the five Sbi-miR396 sequences of sorghum have the same forms as the two mature sequences in Arabidopsis, Ath-miR396a-5p and Ath-miR396b-5p, indicating that the mature miR396 sequences from the conserved region at the 5′ end are consistent between species, in agreement with the findings reported by Shan et al. [45]. Tissue-specific expression analysis revealed extremely low expression levels of Sbi-miR396a, Sbi-miR396b, and Sbi-miR396c in all tissues of sorghum, in contrast to the preferential expression of Sbi-miR396d and Sbi-miR396e in young leaves, flowers, and young panicles. The phylogenetic evolution analysis showed that Sbi-miR396d and Sbi-miR396e and Osa-miR396e and Osa-miR396f were clustered in the same group (Group II), so we speculated that they may have similar functions. Zhang et al. [39] reported that Osa-miR396f double mutant significantly increased panicle number, grain length, grain width, and yield compared with WT. In this study, qRT-PCR results showed Sbi-miR396a, Sbi-miR396b, and Sbi-miR396c were expressed at low levels or not expressed in the young spikes of the inflorescence developmental stage and early development seeds, while Sbi-miR396d and Sbi-miR396e were highly expressed at both developmental stages. These results suggest that only Sbi-miR396d and Sbi-miR396e are extensively involved in the growth and development of organs and tissues, especially the panicle development in sorghum. Therefore, it is possible to change the panicle and grain development and increase the yield by regulating the expression level of Sbi-miR396d and Sbi-miR396e in sorghum.

GRFs are plant-specific transcription factors that play a key role in regulating various stages of plant growth and development. The N-terminus of GRF proteins commonly contains two conserved functional domains, QLQ and WRC [46]. In the present study, all nine GRFs targeted by Sbi-miR396 also exhibited the intact conserved domains WRC (Motif1) and QLQ (Motif2), with QLQ located in front of WRC in all cases. The conserved RSRK-VE region present at the end of the conserved domain WRC corresponds to a highly conserved nucleotide coding sequence and is nearly complementary to the mature Sbi-miR396 sequence. It has been reported that miR396 cuts the mRNA sequence of the GRF gene between the sequences “CGUUCAAGAA” and “AGCNUGUGGAA” at this conserved site, thereby negatively regulating GRF gene expression [47]. The motif and gene structure analyses of SbiGRFs suggested that the nine SbiGRFs have different numbers of motifs and introns, implying that they perform differential functions in the regulation of the growth and development of sorghum.

In sorghum, Sbi-miR396 influences the growth and development processes by targeting SbiGRFs and regulating their expression levels. Therefore, analysis of the expression patterns of the sorghum GRF genes facilitates the understanding of the function of the sorghum miR396–GRF module. In rice, most GRF genes are expressed in several tissues and participate in the regulation of multiple physiological processes during plant growth and development. However, expression analysis of SbiGRFs in different tissues and developmental stages revealed that they are mainly expressed in young inflorescences and early developing seeds, indicating remarkable tissue specificity. These results indicate that SbiGRFs mainly participate in the regulation of floral organ and seed growth and development in sorghum. The phylogenetic analysis demonstrated that SbiGRF8/6/9 were clustered with OsGRF6/7/8/9 in Group I, highlighting their close phylogenetic relationship. Research in rice has demonstrated that OsGRF6 plays a vital role in regulating plant height, stem elongation, and floral organ development [31, 35]; OsGRF7 mainly controls plant height, leaf length, and tiller number [13, 29]; OsGRF8 regulates grain size and yield [33]; and OsGRF9 is not involved in rice growth and development regulation but regulates rice blast resistance [48]. Regarding expression analysis of SbiGRFs in different tissues and developmental stages, the expression level of SbiGRF8 was higher than that of other SbiGRFs in the young spikes of the inflorescence developmental stage and early development seeds. Thus, SbiGRF8 may play a major role in regulating the development of floral organs and seeds in sorghum. SbiGRF9 exhibited the lowest expression level at both developmental stages, suggesting that this gene is not associated with the regulation of sorghum growth and development. It may be functionally consistent with the homologous gene OsGRF9 in rice and plays a role in disease resistance. SbiGRF5/3/4 were clustered with OsGRF3/4/5 in Group V. In rice, OsGRF3/4 regulate grain size [21, 32, 34], and OsGRF5 regulates plant height and inflorescence development [49]. In the present study, SbiGRF5 were found to be highly expressed in inflorescence developmental and early development seeds of sorghum. In particular, SbiGRF5 was preferentially expressed after inflorescences emergence, whereas SbiGRF3/4 showed relatively low expression levels in various tissues. These results suggest that in sorghum, SbiGRF5 was associated with the regulation of seed development and floral organ development. Furthermore, SbiGRF2 and SbiGRF10 also showed relatively low expression levels, whereas SbiGRF1 exhibited very high expression levels at both developmental stages. These results imply that SbiGRFs may be both redundant and functionally specialized in regulating the growth and development of sorghum flower organs and seeds. Therefore, increasing the expression level of SbiGRFs, especially SbiGRF8, SbiGRF1 and SbiGRF5, by means of genetic engineering may affect the panicle and grain development, and increase the yield in sorghum.

Conclusions

In this study, a total of five miR396 members were identified from sorghum genome. Notably, only Sbi-miR396d and Sbi-miR396e could regulate sorghum growth and development. Sbi-miR396d/e mainly target the growth-regulating factor genes SbiGRF1/2/3/4/5/6/8/9/10. However, only SbiGRF1/5/8 were highly expressed in floral organ and early developing seeds. Thus, the miR396d/e–GRF1/5/8 modules may be involved in regulation of sorghum floral organ and seed development, thus changing the sorghum spike structure and subsequently affecting its yield. Taken together, our study is the first comprehensive characterization of miR396–GRF in sorghum. All these findings provide the foundation to elucidate the molecular mechanism of miR396–GRF in sorghum floral organ and seed development.

Supporting information

S1 Table. Primers of qRT-PCR.

https://doi.org/10.1371/journal.pone.0285494.s001

(DOCX)

Citation: Wang H, Zhang Y, Liang D, Zhang X, Fan X, Guo Q, et al. (2023) Genome‑wide identification and characterization of miR396 family members and their target genes GRF in sorghum (Sorghum bicolor (L.) moench). PLoS ONE 18(5): e0285494. https://doi.org/10.1371/journal.pone.0285494

About the Authors:

Huiyan Wang

Roles: Conceptualization, Data curation, Visualization, Writing – original draft, Writing – review & editing

Affiliations: Shanxi Key Laboratory of Sorghum Genetic and Germplasm Innovation, Sorghum Research Institute, Shanxi Agricultural University, Yuci, Shanxi Province, China, Shanxi Key Laboratory of Minor Crops Germplasm Innovation and Molecular Breeding, State Key Laboratory of Sustainable Dryland Agriculture, Shanxi Agricultural University, Taiyuan, Shanxi Province, China

Yizhong Zhang

Roles: Conceptualization, Visualization, Writing – review & editing

E-mail: [email protected] (YZ); [email protected] (QL)

Affiliations: Shanxi Key Laboratory of Sorghum Genetic and Germplasm Innovation, Sorghum Research Institute, Shanxi Agricultural University, Yuci, Shanxi Province, China, Shanxi Key Laboratory of Minor Crops Germplasm Innovation and Molecular Breeding, State Key Laboratory of Sustainable Dryland Agriculture, Shanxi Agricultural University, Taiyuan, Shanxi Province, China

ORICD: https://orcid.org/0000-0001-5176-0040

Du Liang

Roles: Data curation

Affiliations: Shanxi Key Laboratory of Sorghum Genetic and Germplasm Innovation, Sorghum Research Institute, Shanxi Agricultural University, Yuci, Shanxi Province, China, Shanxi Key Laboratory of Minor Crops Germplasm Innovation and Molecular Breeding, State Key Laboratory of Sustainable Dryland Agriculture, Shanxi Agricultural University, Taiyuan, Shanxi Province, China

Xiaojuan Zhang

Roles: Data curation

Affiliations: Shanxi Key Laboratory of Sorghum Genetic and Germplasm Innovation, Sorghum Research Institute, Shanxi Agricultural University, Yuci, Shanxi Province, China, Shanxi Key Laboratory of Minor Crops Germplasm Innovation and Molecular Breeding, State Key Laboratory of Sustainable Dryland Agriculture, Shanxi Agricultural University, Taiyuan, Shanxi Province, China

Xinqi Fan

Roles: Data curation

Affiliations: Shanxi Key Laboratory of Sorghum Genetic and Germplasm Innovation, Sorghum Research Institute, Shanxi Agricultural University, Yuci, Shanxi Province, China, Shanxi Key Laboratory of Minor Crops Germplasm Innovation and Molecular Breeding, State Key Laboratory of Sustainable Dryland Agriculture, Shanxi Agricultural University, Taiyuan, Shanxi Province, China

Qi Guo

Roles: Visualization

Affiliations: Shanxi Key Laboratory of Sorghum Genetic and Germplasm Innovation, Sorghum Research Institute, Shanxi Agricultural University, Yuci, Shanxi Province, China, Shanxi Key Laboratory of Minor Crops Germplasm Innovation and Molecular Breeding, State Key Laboratory of Sustainable Dryland Agriculture, Shanxi Agricultural University, Taiyuan, Shanxi Province, China

Linfang Wang

Roles: Visualization

Affiliations: Shanxi Key Laboratory of Sorghum Genetic and Germplasm Innovation, Sorghum Research Institute, Shanxi Agricultural University, Yuci, Shanxi Province, China, Shanxi Key Laboratory of Minor Crops Germplasm Innovation and Molecular Breeding, State Key Laboratory of Sustainable Dryland Agriculture, Shanxi Agricultural University, Taiyuan, Shanxi Province, China

Jingxue Wang

Roles: Visualization

Affiliation: School of Life Science, Shanxi University, Taiyuan, Shanxi Province, China

Qingshan Liu

Roles: Conceptualization, Funding acquisition

E-mail: [email protected] (YZ); [email protected] (QL)

Affiliation: Shanxi Key Laboratory of Minor Crops Germplasm Innovation and Molecular Breeding, State Key Laboratory of Sustainable Dryland Agriculture, Shanxi Agricultural University, Taiyuan, Shanxi Province, China

References

1. Rooney W, Blumenthal J, Bean B, Mullet J. Designing sorghum as a dedicated bioenergy feedstock. Biofuel Bioprod Bior. 2007;1:147–57.

2. Chen J, Zhu M, Ruixiang L, Zhang M, Lv Y, Liu Y, et al. BIOMASS YIELD 1 regulates sorghum biomass and grain yield via the shikimate pathway. J Exp Bot. 2020;71(18):5506–20. pmid:32497182

3. Paterson A, Bowers J, Bruggmann Ré, Dubchak I, Grimwood J, Gundlach H, et al. The Sorghum bicolor genome and the diversification of grasses. Nature. 2009;457:551–6. pmid:19189423

4. Tao Y, Zhao X, Wang X, Hathorn A, Hunt C, Cruickshank AW, et al. Large-scale GWAS in sorghum reveals common genetic control of grain size among cereals. Plant Biotechnol J. 2020;18(4):1093–105. pmid:31659829

5. Xie Q, Xu Z. Sustainable agriculture: From sweet sorghum planting, ensiling to ruminant feeding. Mol Plant. 2019;12. pmid:31002980

6. Williams G, Capps O. Generic promotion of sorghum for food and industrial uses. J Int Food Agribus Market. 2018;32:1–17.

7. Takaki M, Tan L, Murakami T, Tang Y-Q, Sun Z-Y, Morimura S, et al. Production of biofuels from sweet sorghum juice via ethanol–methane two-stage fermentation. Ind Crop Prod. 2015;63:329–36.

8. Upadhyaya HD, Wang L, Prakash CS, Liu Y, Gao L, Meng R, et al. Genome-wide association mapping identifies an SNF4 ortholog that impacts biomass and sugar yield in sorghum and sugarcane. J Exp Bot. 2022:erac110. pmid:35290448

9. Liu C, Zhang H, Xin X. Analysis of the development characteristics and trends of sorghum industry in China. J Agr Sci Tech. 2020;22(10):1–9.

10. Li S, Liu M, Liu F, Zou J, Lu X, Diao X. Current status and future prospective of sorghum production and seed industry in China. Sci Agr Sinica. 2021;54(03):471–82.

11. Wang Y, Li J. Molecular basis of plant architecture. Annu Rev Plant Biol. 2008;59:253–79. pmid:18444901

12. Khush G. Productivity improvements in rice. Nutr Rev. 2003;61:S114–6. pmid:12908742

13. Yunping C, Dan Z, Feng G, Chen P, fan F, Li S. Rice GROWTH-REGULATING FACTOR7 modulates plant architecture through regulating GA and IAA metabolism. Plant Physiol. 2020;184:pp.00302.2020. pmid:32581114

14. Miao C, Wang D. Mutations in MIR396e and MIR396f increase grain size and modulate shoot architecture in rice. Plant Biotechnol J. 2020;18(2):491–501. pmid:31336020

15. Zhang S, Wu Y, Huang X, Wu W, Lyu L, Li W. Research progress about microRNAs involved in plant secondary metabolism. Int J Biol Macromol. 2022;216:820–9. pmid:35921910

16. Tang J, Chu C. MicroRNAs in crop improvement: fine-tuners for complex traits. Nat Plants. 2017;3:17077. pmid:28665396

17. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116(2):281–97. pmid:14744438

18. Zhang H, Zhang J, Yan J, Gou F, Mao Y, Tang G, et al. Short tandem target mimic rice lines uncover functions of miRNAs in regulating important agronomic traits. PNAS. 2017;114(20):5277–82. pmid:28461499

19. Pegler J, Nguyen DQ, Oultram J, Grof C, Eamens A. Molecular manipulation of the MiR396/GRF expression module alters the salt stress response of Arabidopsis thaliana. Agronomy. 2021;11:1751.

20. Yu Y, Sun F, Chen N, Sun G, Wang C-Y, Wu D-X. MiR396 regulatory network and its expression during grain development in wheat. Protoplasma. 2021;258:103–13. pmid:32929630

21. Li S, Gao F, Xie K, Zeng X, Cao Y, Zeng J, et al. The OsmiR396c-OsGRF4-OsGIF1 regulatory module determines grain size and yield in Rice. Plant Biotechnol J. 2016;14:2134–46. pmid:27107174

22. Subramanian S, Fu Y, Sunkar R, Barbazuk W, Zhu J-K, Yu O. Novel and nodulation-regulated microRNAs in soybean roots. BMC Genomics. 2008;9:160. pmid:18402695

23. Yang F, Liang G, Liu D, Yu D. Arabidopsis MiR396 mediates the development of leaves and flowers in transgenic tobacco. J Plant Biol. 2009;52:475–81.

24. Hewezi T, Baum TJ. Complex feedback regulations govern the expression of miRNA396 and its GRF target genes. Plant Signal Behav. 2012;7(7):749–51. pmid:22751317

25. Liu D, Song Y, Chen Z, Yu D. Ectopic expression of miR396 suppresses GRF target gene expression and alters leaf growth in Arabidopsis. Physiol Plant. 2009;136(2):223–36. pmid:19453503

26. Liang G, He H, Li Y, Wang F, Yu D. Molecular mechanism of microRNA396 mediating pistil development in Arabidopsis. Plant Physiol. 2013;164:249–58. pmid:24285851

27. Rodriguez R, Mecchia M, Debernardi J, Schommer C, Weigel D, Palatnik J. Control of cell proliferation in Arabidopsis thaliana by microRNA miR396. Development. 2010;137:103–12. pmid:20023165

28. Lu Y, Feng Z, Meng Y, Bian L, Xie H, Mysore K, et al. SLENDER RICE1 and Oryza sativa INDETERMINATE DOMAIN2 regulating OsmiR396 are involved in stem elongation. Plant Physiol. 2020;182(4):2213–27. pmid:31953375

29. Jathar V, Saini K, Chauhan A, Rani R, Ichihashi Y, Ranjan A. Spatial control of cell division by GA-OsGRF7/8 module in a leaf explaining the leaf length variation between cultivated and wild rice. New Phytol. 2022;234(3):867–83. pmid:35152411

30. Feng G, Wang K, Liu Y, Yunping C, Chen P, Shi Z, et al. Blocking miR396 increases rice yield by shaping inflorescence architecture. Nat Plants. 2015;2:15196. pmid:27250748

31. Yongyan T, Huanhuan L, Guo S, Wang B, Li Z, Chong K, et al. OsmiR396d affects gibberellin and brassinosteroid signaling to regulate plant architecture in rice. Plant Physiol. 2017;176:946–56. pmid:29180380

32. Hu J, Wang Y, Fang Y, Zeng L, Xu J, Yu H, et al. A rare allele of GS2 enhances grain size and grain yield in rice. Mol Plant. 2015;8(10):1455–65. pmid:26187814

33. Yang X, Zhao X, Dai Z, Ma F, Miao X, Shi Z. OsmiR396/Growth Regulating Factor modulate rice grain size through direct regulation of embryo-specific miR408. Plant Physiol. 2021;186. pmid:33620493

34. Duan P, Ni S, Wang J, Zhang B, Xu R, Wang Y, et al. Regulation of OsGRF4 by OsmiR396 controls grain size and yield in rice. Nat Plants. 2015;2:15203. pmid:27250749

35. Liu H, Guo S, Xu Y, Li C, Zhang Z, Zhang D, et al. OsmiR396d-regulated OsGRFs function in floral organogenesis in rice through binding to their targets OsJMJ706 and OsCR4. Plant Physiol. 2014;165(1):160–74. pmid:24596329

36. Chen L, Luan Y, Zhai J. Sp-miR396a-5p acts as a stress-responsive genes regulator by conferring tolerance to abiotic stresses and susceptibility to Phytophthora nicotianae infection in transgenic tobacco. Plant Cell Rep. 2015;34(12):2013–25. pmid:26242449

37. Gao P, Bai X, Yang L, Lv D, Li Y, Cai H, et al. Over-expression of osa-MIR396c decreases salt and alkali stress tolerance. Planta. 2010;231:991–1001. pmid:20135324

38. Yang F-X, Yu D-Q. Overexpression of Arabidopsis MiR396 enhances drought tolerance in transgenic tobacco plants. Acta Botanica Yunnanica. 2010;31:421–6.

39. Zhang J, Zhou Z, Bai J, Tao X, Wang L, Zhang H, et al. Disruption of MIR396e and MIR396f improves rice yield under nitrogen-deficient conditions. Natl Sci Rev. 2020;7:102–12. pmid:34692021

40. Crooks GE, Hon G, Chandonia J-M, Brenner SE. WebLogo: a sequence logo generator. Genome Res. 2004;14(6):1188–90. pmid:15173120

41. Kong J, Martin-Ortigosa S, Finer J, Orchard N, Gunadi A, Batts LA, et al. Overexpression of the transcription factor GROWTH-REGULATING FACTOR5 improves transformation of dicot and monocot species. Front Plant Sci. 2020;11:572319. pmid:33154762

42. Qiu F, Xing S, Xue C, Liu J, Chen K, Chai T, et al. Transient expression of a TaGRF4-TaGIF1 complex stimulates wheat regeneration and improves genome editing. Sci China Life Sci. 2022;65(4):731–8. pmid:34406572

43. Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194–202. pmid:32585190

44. Noon JB, Hewezi T, Baum TJ. Homeostasis in the soybean miRNA396-GRF network is essential for productive soybean cyst nematode infections. J Exp Bot. 2019;70(5):1653–68. pmid:30715445

45. Shan Q, Jia H, Yao W, Liu w, Li H. Research progress in biological functions of miR396-GRF module in plants and its potential application values. Biotechnol Bull. 2022;38(11):28–38.

46. Liebsch D, Palatnik JF. MicroRNA miR396, GRF transcription factors and GIF co-regulators: a conserved plant growth regulatory module with potential for breeding and biotechnology. Curr Opin Plant Biol. 2020;53:31–42. pmid:31726426

47. Liu W, Zhou Y, Li X, Wang X, Dong Y, Wang N, et al. Tissue-specific regulation of Gma-miR396 family on coordinating development and low water availability responses. Front Plant Sci. 2017;8:1112. pmid:28694817

48. Chandran V, Wang H, Gao F, Cao XL, Chen YP, Li GB, et al. miR396-OsGRFs module balances growth and rice blast disease-resistance. Front Plant Sci. 2018;9:1999. pmid:30693011

49. Kuijt SJ, Greco R, Agalou A, Shao J, t Hoen CC, Overnäs E, et al. Interaction between the GROWTH-REGULATING FACTOR and KNOTTED1-LIKE HOMEOBOX families of transcription factors. Plant Physiol. 2014;164(4):1952–66. pmid:24532604