Investigations on Regulation of MicroRNAs in Rice Reveal [Ca2+]cyt Signal Transduction Regulated MicroRNAs

MicroRNAs (miRNAs) are critical components of the multidimensional regulatory networks in eukaryotic systems. Given their diverse spectrum of function, it is apparent that the transcription, processing, and activity of the miRNAs themselves, is very dynamically regulated. One of the most important and universally implicated signaling molecules is [Ca2+]cyt. It is known to regulate a plethora of developmental and metabolic processes in both plants and animals; however, its impact on the regulation of miRNA expression is relatively less explored. The current study employed a combination of internal and external calcium channel inhibitors to establishing that [Ca2+]cyt signatures actively regulate miRNA expression in rice. Involvement of [Ca2+]cyt in the regulation of miRNA expression was further confirmed by treatment with calcimycin, the calcium ionophore. Modulation of the cytosolic calcium levels was also found to regulate the drought-responsive expression as well as ABA-mediated response of miRNA genes in rice seedlings. The study further establishes the role of calmodulins and Calmodulin-binding Transcription Activators (CAMTAs) as important components of the signal transduction schema that regulates miRNA expression. Yeast one-hybrid assay established that OsCAMTA4 & 6 are involved in the transcriptional regulation of miR156a and miR167h. Thus, the study was able to establish that [Ca2+]cyt is actively involved in regulating the expression of miRNA genes both under control and stress conditions.


INTRODUCTION
Plants being sessile in nature need to face all the environmental and biotic challenges in one place. Plants modulate their internal system to cope with these challenges. Changes in the growth or environment of the plants are conveyed to the cell via primary messengers such as hormones while it is relayed to the nucleus via the secondary messengers such as cAMP, cGMP, Inositol triphosphate, diacylglycerol, and the most widely studied [Ca 2+ ] cyt . Apart from these, players like hydrogen sulfide, hydrogen peroxide, carbon monoxide, and nitric oxide, collectively termed as "gasotransmitters, " have also emerged as potential signaling messengers in cells (Gaupels et al., 2011;Mittler et al., 2011;Lisjak et al., 2013). Their signaling cascade is often routed via [Ca 2+ ] cyt or MAPK signal transduction components. [Ca 2+ ] cyt as a second messenger of internal and external signals has been well-studied and documented in a plethora of studies and reviews, describing its role in signaling in various plant physiological functions such as abiotic stress response (Dodd et al., 2010), stomatal aperture (McAinsh et al., 1990), self-incompatibility (Franklin-Tong et al., 1993), pathogenic and symbiotic interactions (Ma and Berkowitz, 2007), and growth of pollen tube and root tips (Hepler et al., 2001). Upon perceiving an environmental stimulus or developmental cue, rapid and precise [Ca 2+ ] oscillations are generated within the cytoplasm through the influx of [Ca 2+ ] from internal and extracellular stores. Besides cytosolic [Ca 2+ ] signatures and signaling, the plant cell nucleus has also been proved to contain within an independent Ca 2+ -signaling system (Pauly et al., 2000(Pauly et al., , 2001Xiong et al., 2004) to regulate gene transcription. Calmodulin (CaM) is said to be the primary decoder of the signal along with the support of CaM-binding proteins in the nucleus. For instance, nuclear apyrase in pea binds to CaM and is activated by Ca 2+ /CaM complex (Hsieh et al., 2000); nuclear localized potato CaM-binding protein (PCBP) identified from developing potato tubers binds to CaM in a Ca 2+dependent manner (Reddy et al., 2002). The strongest evidence for the theory is the identification of nuclear transcription factors that are regulated directly or indirectly by CaM (Yang and Poovaiah, 2003). For instance, CAM7 in Arabidopsis has been demonstrated to regulate photomorphogenic growth and lightresponsive gene expression by binding to Z-/G-box elements in their promoters, including CAB1 and RBCS1A. CAM7 was also shown to work independently yet in an interdependent pathway with HY5 to regulate the photomorphogenic growth of Arabidopsis (Kushwaha et al., 2008). WRKY class of TFs has also been shown to interact with Ca 2+ /CaM complex, especially the WRKY group IId . Another separate class of transcription factors known as CAMTAs (Calmodulinbinding transcription activators), described to be evolutionarily conserved from humans to plants (Finkler et al., 2006), responds rapidly to multiple abiotic stresses such as drought, salinity, heat, cold, and UV, thus regulating the signaling required (Yang and Poovaiah, 2003). Despite the name, CAMTAs can act as transcriptional activators as well as repressors; as CAMTA3 was demonstrated for activating the expression of CBF2-imparting cold tolerance (Doherty et al., 2009), as well as repressing the expression of EDS1, a regulator of the salicylic acid level (Du et al., 2009). Recently, AtCAMTA3 has been shown to repress the genes involved in SA-biosynthesis and SA-mediated immunity in healthy plants (Kim et al., 2017).
In biological systems, microRNAs (miRNAs) represent a significant and critical layer of regulation of protein-coding genes at the post-transcriptional level. In our previous studies and others as well, we have seen that several miRNAs can express differentially to a similar cue in different rice cultivars, thus emphasizing the evolution of distinct regulatory mechanisms controlling miRNA expression in different cultivars (Mutum et al., 2013;Kansal et al., 2015;Balyan et al., 2017). Such unique regulation would have a cascading effect on plant biology since a single miRNA regulates several protein-coding genes. Thus, the study of the regulation of miRNA expression is of critical importance. MiRNAs have been shown to be differentially regulated by different stages of growth and development, biotic and abiotic stresses, etc. At the molecular level, miRNAs have been shown to respond to secondary messengers of signaling in various publications. A study showed seven miRNAs, including miR169, miR397, miR1425, miR408-5p, and miR827 that were upregulated, while miR319a-3p and miR408-5p were downregulated by H 2 O 2 stress (Li et al., 2011). Similarly, miRNAs were also shown responsive to exogenously supplied H 2 S, and it was established that their expression under drought simulation was mediated via H 2 S (Shen et al., 2013). Analysis in the lcd (knockdown mutant of H 2 S-producing enzyme L-CYSTEINE DESULFHYDRASE) mutants under control and droughtsimulated conditions showed lower expression levels of some miRNAs like miR167a/c/d, miR393a, miR396a, miR398a/b/c, which could be rescued by resupplying H 2 S. This confirmed that H 2 S regulates miRNA activity to modulate the drought response of Arabidopsis. Recently, miRNAs have been shown to be differentially regulated during early embryo development under conditions of calcium deficiency in the growth medium in peanuts (Chen et al., 2019). The targets of such miRNAs included seed/embryo development-related genes, cell-division and proliferation-related genes, and phytohormone signalingrelated genes. However, to the best of our knowledge, the demonstration that miRNAs are responsive to [Ca 2+ ] cyt has not been made till now.
In the current study, we attempt to establish a stronger link between calcium signaling and its control over the regulation of miRNAs in rice. The miRNAs were confirmed to be calcium responsive using miRNA expression tools [small RNA libraries and quantitative real-time PCR (qRT-PCR)] under conditions of calcium scarcity (created using Ca 2+ channel inhibitors) and excess of [Ca 2+ ] cyt using a Ca 2+ ionophore A23187. Furthermore, the involvement of Ca 2+ signaling in differential expressions of miRNAs under abiotic stress such as dehydration and hormonal stress such as ABA was also shown. Additionally, we show here the participation of calmodulins and CAMTAs in the regulation of miRNA expression. Furthermore, our results show physical interaction between OsCAMTA4 and OsCAMTA6 and the promoters of miR156a and miR167h, thereby strengthening the evidence of hold of calcium signaling over miRNA regulation.

Rice Plant Growth and Stress Treatment
Rice seeds [O. sativa subsp. indica var Nagina22 (N22)] were obtained from the Indian Agricultural Research Institute (IARI) Pusa, New Delhi. Seeds were sterilized and planted on a muslin cloth tied to a tray filled with a rice growth medium (as described in Kansal et al., 2015). These were grown for 1 week in a culture room with 28 • C ± 2 • C temperature and 16-h/8h day/night settings and then subjected to various stresses. For seedling stress experiments, the seedlings were transferred to a solution of calcium channel inhibitors [mix of LaCl 3 (Lanthanum chloride, 5 mM), Verapamil (50 µM), LiCl (Lithium chloride, 5 mM), Ruthenium Red (100 µM) for 3 h], calcium ionophore A23187 (10 µM for 5 h), TriFluoPerazine (200 µM for 6 h) or ABA (Abscisic acid, 100 µM for 3 h). All these solutions and control (MQ water) contained the surfactant Silwett at 0.01% final concentration for better and uniform absorption. All the chemicals were obtained from Sigma Aldrich.
For dehydration stress, the seedlings were taken out of the respective solution, dabbed with tissue paper to remove excess liquid sticking to the roots, and kept open in the air for air-drying until leaf rolling was observed. Whole seedlings were frozen in liquid nitrogen and stored at −80 • C.
Drought stress at the mature stage was given as per Mutum et al. (2016). Briefly, drought stress was given by withdrawing water supply to the field bed 10 days before the "mean heading date" of the cultivar. The samples were collected from plants that had attained heading, displayed leaf rolling, and had soil moisture content dropped below 15% (as measured using Hydra Probe Soil Moisture Sensor). Tissue (flag leaf and inflorescence) was frozen in liquid nitrogen in situ and stored at −80 • C. Induction of drought was checked by estimating transcript levels of drought stress marker genes such as Rubisco small subunit (RBCS) and OsbZIP23 in the collected tissue as described in Mutum et al. (2016).

RNA Isolation and cDNA Synthesis
RNA was isolated using TriReagent (Sigma) as per the protocol of the manufacturer. It was subjected to DNaseI treatment followed by small RNA enrichment by the 4M LiCl method as given in Kansal et al. (2015). For cDNA synthesis, 2 µg of RNA was incubated with 100 pmoles of either oligodT primer (Eurofins Genomics) for total RNA or miR_oligodT primer (Eurofins Genomics) for small RNA using the SuperscriptII cDNA synthesis kit (Invitrogen). All the primers have been enlisted in Supplementary Table 1.

Small RNA Library Preparation, Quality Check, and Analysis
Small RNA libraries were prepared with 2 µg of enriched small RNA using NEBNext R Small RNA Library Prep Set for Illumina R (Multiplex Compatible) (New England Biolabs) according to the protocol of the manufacturer in biological duplicates. Each biological replicate consisted of a set of 5-6 seedlings individually subjected to each step of stress, RNA isolation, library preparation, and sequencing. Briefly describing the library preparation, 3 ′ and 5 ′ small RNA adapters were ligated to the small RNA population, sequentially followed by First Strand cDNA synthesis. After this, the Small RNA primer for Illumina along with index primers was added, and the sample was subjected to PCR. The cDNA prepared was thus run on 6% PAGE gel for size excision of miRNA band (∼140 nts) and purified further. The libraries prepared were thus checked for quality using 2100 Bioanalyzer High Sensitivity DNA chip (Agilent) according to the protocol of the manufacturer. The libraries were sequenced using the Illumina platform by Sandor LifeSciences Pvt Ltd. Small RNA Library Analysis was done using CLC Genomics Workbench version 7.0.4 (Qiagen) with the help of the "Small RNA Analysis" tool. For a further detailed analysis of differentially expressed genes, the "Statistical Analysis" tool, in particular, "Empirical Analysis of DGE" which uses "Edge" as the tool for calculating differentially expressed genes, was used. The statistically significant DEGs were filtered at having p ≤ 0.05 and expression fold change of 2X in stress condition vs. the control.

Quantitative Real-Time PCR
The first-strand cDNA synthesized (section RNA isolation and cDNA synthesis) was used for qRT-PCR with appropriate dilution. The reaction was set up by adding TaqMan Master Mix chemistry (2x, Applied Biosystems). The reaction was run in the ABI StepOnePlus Real-Time system as per its default cycling for TaqMan chemistry. Rice actin and 5S rRNA were used as endogenous controls for protein-coding gene and miRNA profiling, respectively. For each sample, at least three technical and three biological replicates were set up. The expression fold change was calculated as stress vs. control with 5S rRNA used as the internal control for normalization purposes. The use of * above a bar in the bar graphs depicts its significance as * p ≤ 0.05 as recorded by the t-test. A two-tailed t-test (with the assumption of unequal variance) has been used to determine if there is a significant difference between the means of control and treatment groups.

Mirna-Target Analysis and Classification
Putative targets of miRNAs were selected from in-house analyzed publicly available rice degradome libraries (GSM455939, GSM455938, GSM434596, GSM960648, and GSM476257; (Wu L. et al., 2009;Li et al., 2010;Wang et al., 2012) and in-house generated three additional degradome libraries from spikelets at the heading and anthesis stage, flag leaf at the heading stage of N22 (Mutum et al., 2016). The targets falling within the criteria of category 0-2 along with read no ≥10 and p ≤ 0.05 were selected as putative targets. Category "0" is defined as >1 raw read at the position, with abundance at a position equal to the maximum on the transcript, and with only one maximum on the transcript. Category "1" is described as >1 raw read at the position, with abundance at the position equal to the maximum on the transcript, and more than one maximum position on the transcript. Category "2" includes >1 raw read at the position, and abundance at the position less than the maximum but higher than the median for the transcript (Yang et al., 2013). The target loci were then used to extract their GO annotations from the RGAP database and classified according to their GO term of "molecular function." Similarly, the loci were checked for their presence in the RiceCyc database available at the Gramene database.
The promoter sequences (2 kb upstream of the precursor start site) were retrieved from IRDb (Indica Rice Database, www.genomeindia.du.ac.in/irdb) and subjected to a sequencebased motif search for the calcium-responsive motifs enlisted in literature (Galon et al., 2010).

Yeast One-Hybrid Assay
Yeast one-hybrid assay was performed using the EZ yeast transformation kit (MP Biomedicals) as per the instructions. For the Y1H promoter construct, 110 and 200 bp regions were selected from 2-kb promoter regions of MIR156 and MIR167, respectively, that contained at least one CAMTA-binding site (Supplementary Figure 1). These mentioned promoter fragments were PCR amplified from the genomic DNA of N22 and cloned into the pAbAI vector as bait, harboring the aureobasidin (AUR1-C) reporter gene. Full-length CDS of OsCAMTA4 and OsCAMTA6 were PCR amplified from genomic DNA of N22 and cloned into the pGADT7-AD vector as prey. The miR156a-pAbAi, miR167h-pAbAi, OsCAMTA4-pGADT7, and OsCAMTA6-pGADT7 plasmids were linearized and cotransformed into the Y1H gold yeast strain. The transformants were selected on SD medium-lacking Uracil (U) and Leucine (L) (SD/-L-U), respectively, as per the protocol of the manufacturer. Colonies obtained were screened for a positive insert via colony PCR. Furthermore, the positive clones were inoculated as primary culture and subsequently used in the secondary culture (3 ml) and grown at 30 • C, 200 rpm for 3 h till the OD = 1 at 600 nm. Each culture was then serially diluted (10 −1 , 10 −2 , and 10 −3 ), and droplets of 10 µl of each dilution, including the undiluted culture, were placed on the selection media (SD/-L-U) and incubated at 30 • C for 3-5 days till the formation of colonies. The interaction was determined based on the ability of cotransformed yeast cells to grow on dropout media-lacking Uracil and Leucine with 100, 200, and 300 ng/ml of Aureobasidin A.

Accessions
Accession IDs for camta mutants used in the study are as follows:

RESULTS miRNAs Are Responsive to [Ca 2+ ] cyt in N22 Rice
Cellular calcium levels are primarily regulated by the activity of the calcium channels present either on the plasma membrane (blocked by LaCl 3 and Verapamil) and/or the mitochondrial (inhibited by Ruthenium Red) or endoplasmic reticulum membrane (blocked by LiCl). These are well-characterized inhibitors and have been previously experimentally shown to alter the cellular calcium levels in plants (Knight et al., 1996;Sedbrook et al., 1996;Legué et al., 1997;Hu et al., 2004;Song et al., 2011;Gao and Zhang, 2019;Zheng et al., 2020). In order to identify miRNA genes that are regulated by [Ca 2+ ] cyt, a cocktail of Ca 2+ channel blockers (viz. LaCl 3 , LiCl, Ruthenium Red, and Verapamil) was used to treat young rice seedlings, followed by sequencing of the sRNA population. The efficiency of the inhibitor treatment was confirmed by monitoring the expression of a calcium-responsive "Ca 2+ -dependent protein kinase" OsCPK6 (Wan et al., 2007) (Supplementary Figure 2). Briefly describing, 293 and 303 miRNAs were detected in control libraries, while 258 and 279 miRNAs were detected from the tissue treated with calcium channel inhibitors (Supplementary Table 2). Furthermore, comparative analysis of the datasets revealed 17 differentially expressed miRNAs (fold change ≥ 2 and p ≤ 0.05) (Supplementary Table 3); out of which, 10 were found to be downregulated, while 7 were upregulated in the presence of Ca 2+ channel inhibitors ( Figure 1A). Subsequently, qRT-PCR analysis reaffirmed the calcium inhibitor-dependent differential expression of 11 miRNAs in rice (miR1425-5p, miR156a, miR159b, miR166g-3p, miR167h-5p, miR1862a,d, miR1876, miR1878, miR396c-3p, and miR444b.1) ( Figure 1B).
The [Ca 2+ ] cyt inducibility of these miRNAs was further verified with the help of a Ca 2+ ionophore, i.e., A23187 which is expected to reverse the expression pattern obtained with inhibitors. A23187, also known as calcimycin, is an ionophore that binds to Ca 2+ ions and acts as their carrier by inserting itself in the plasma membrane, thereby acting as a channel for the divalent. Quantitative real-time PCR profiling under control and ionophore treatment confirmed nine miRNAs (miR1425-5p, miR156a, miR159b, miR167h-5p, miR1862a,d, miR1876, miR1878, and miR396c-3p) that responded to channel inhibitors as well as to the ionophore but in inverse fashion, thereby confirming their [Ca 2+ ] cyt inducibility ( Figure 1C).
To identify the source of [Ca 2+ ] cyt inducibility of these miRNAs, subsequent analysis was done by treating rice seedlings with only internal or extracellular channel inhibitors to decipher the independent involvement of internal and external calcium stores in this regulation ( Figure 1D). Interestingly, miR156a, miR159b, miR167h-5p, and miR396c-3p were significantly downregulated when internal channels were inhibited but showed no remarkable change with the blocking of plasma membrane channels. On the other hand, miR166g-3p was downregulated in the presence of external inhibitors only. Thus, it goes to show that specific miRNAs are uniquely regulated by the internal or external calcium reserves indicative of the complexity of the regulatory phenomenon.
However, from these data, it appears that the effect of [Ca 2+ ] cyt is not global since only a small proportion of miRNAs was found to be responsive, indicative of a transcriptional route (rather than miRNA processing where a higher proportion of miRNAs would have been differentially changed). Thus, dissecting closely, the [Ca 2+ ] cyt responsiveness of the three main miRNA-processing enzymes, namely DCL1, HYL1, and SE, was checked, and the response was found neither significant nor confirmed by ionophore treatment (Supplementary Figure 3).
Besides, we also found various calcium-responsive ciselements in the promoter regions of miRNAs. The promoter sequences of some miRNAs were subjected to a search for the calcium-responsive motifs enlisted in literature (Galon et al., 2010). These motifs are ABA-related ABRE elements [(T/C) GT-box [GGTAATT], and Rapid Stress Responsive Elements [CGCGTT]. For instance, MIR1425 has two ABRE, three CAMTA, and five E-box elements; MIR156A has one CARE and two CAMTA sites, and MIR159B has nine E-box sites.

Targets of [Ca 2+ ] cyt -Regulated miRNAs in N22 Are Functionally Widespread
To analyze the impact of these [Ca 2+ ] cyt regulated miRNAs, data from eight degradome libraries, including both in-house generated and publicly available, were analyzed to identify Frontiers in Plant Science | www.frontiersin.org miRNA targets (Mutum et al., 2016). Targets identified under the category "0-2" and ≥10 degradome reads were considered for the analysis ( Table 1).
These targets included some well-known genes, such as certain transcription factors like OsSPL (miR156a), MYB [miR159b and START domain containing (miR166), auxin response factors (miR160, miR167). A gene ontology classification of these targets (based on molecular function) shows transcription factors being the major category along with protein-binding proteins and catalytic activity. The targets are widespread along almost all major categories like lipid and RNA-binding, kinase activity, transporters, structural molecule activity, etc. (Figure 2). Similarly, the majority of targets are related to the "biological process" terms like nucleic acid metabolism, biosynthetic process, multicellular organismal development, flower development, and metabolic processes.

Dehydration and ABA Response of miRNAs Are Mediated by [Ca 2+ ] cyt
Since it is a well-known fact that miRNAs respond to various abiotic stresses and that [Ca 2+ ] cyt plays a significant role in relaying these abiotic stress cues, we investigated whether the dehydration response of miRNAs in rice is mediated by [Ca 2+ ] cyt . Thus, 7-day-old rice seedlings were subjected to air dehydration in the presence or absence of the calcium channel inhibitor cocktail as described previously. Subsequently, sRNA libraries were generated and sequenced (Supplementary Tables 2, 4). Comparison of both the datasets gives us an insight into [Ca 2+ ] cyt -mediated dehydration response of miRNAs. Consequently, through qRT-PCR confirmation, it was possible to identify that dehydration response of miR156a, miR167h-5p, miR168a-5p, miR5083, and miR5788 is being mediated by [Ca 2+ ] cyt levels (Figures 3A-C). Interestingly, miR156a and miR167h-5p also showed [Ca 2+ ] cyt responsiveness under control conditions.
Another important signaling molecule involved in the regulation of gene expression is abscisic acid (ABA), aka the stress hormone that mediates several developmental processes in plants, including dormancy, abscission, seed germination, as well as the abiotic and biotic stress response. ABA levels increase during certain stress conditions and mediate adaptive response. Thus, experiments were conducted to assess the response of miRNAs to elevated ABA levels by treating 7-dayold seedlings with ABA (100 µM). Based on previous literature, a selected set of miRNAs was profiled, and the results indicated that miR1425-5p, miR159b, miR168a-5p, and miR529b are significantly downregulated, while miR319b and miR530-5p are upregulated in the presence of ABA (Figure 4).
Expression of miR156a, miR162b, miR1320-5p, miR396a, miR528-5p did not deviate much from control. The correlation of the ABA and dehydration responsiveness indicates that miR1425-5p and miR168a-5p are downregulated during dehydration and ABA. miR156a is downregulated by dehydration but is not ABA responsive. On the other hand, miR162b, miR396a-5p, miR528-5p, and miR1320-5p turn out neither dehydration nor ABA responsive. [Ca 2+ ] cyt was found to mediate the ABA response of miRNAs as well, since pre-blocking [Ca 2+ ] cyt channels alleviated the hormone-induced response of some miRNAs, such as miR159b, miR319b, and miR530-5p (Figure 4). Notably, miR159b is responsive to [Ca 2+ ] cyt in resting the state as well; however, its response to dehydration was not significant. Hence, here is a sample of demonstration of the versatile roles that [Ca 2+ ] cyt plays in regulating miRNAs under variable environmental conditions.
Subsequent to finding [Ca 2+ ] cyt and calmodulin-responsive miRNAs in rice, further components involved in mediating the [Ca 2+ ] signal transduction of these miRNAs were explored. "Calmodulin-binding transcriptional activators" or CAMTAs are calcium-binding transcription factors that play an important role in mediating Ca 2+ /CaM-mediated gene expression. Five miRNAs viz. miR156a, miR160a-5p, miR166a-3p, miR167h-5p, and miR168a-5p are orthologous in both rice and Arabidopsis plus their promoter regions were found to have CAMTAbinding sites. To investigate whether CAMTAs regulate their expression, these miRNAs were profiled in the Arabidopsis camta mutants (available at the ABRC seed stock center). There are six CAMTA genes in Arabidopsis, and mutants specific to CAMTA loci camta1, camta3, camta4, camta5, and camta6 were procured (Supplementary Figure 4; for accessions, refer to section Accessions).
Quantitative real-time profiling for these miRNAs in 10day-old mutant seedlings revealed that miR156a is significantly reduced in camta4 and camta6 (Figure 6). Expression of miR160a-5p is reduced in camta4 while enhanced in camta5 and camta6. Similarly, miR168a-5p is reduced in camta3,4,5 but not in camta6. On the other hand, miR167h-5p is reduced in camta5&6 but over accumulated in camta1 mutant. Thus, these results indicate that CAMTAs are, indeed, involved in the regulation of miRNA expression. The promoter of miR167h has one CAMTA-binding site, while that of miR156a harbors two sites. Furthermore, CAMTA4 appears to have a broader and major role, as it appears to regulate the expression of all the five miRNAs. In the case of miR156a, miR160a-5p, miR167h-5p, and miR168a-5p, it has a positive influence, while, for miR166a-3p, it acts to negatively regulate it. Taking this cue further, we ventured to find if the orthologs of these AtCAMTA4 and 6, i.e., OsCAMTA4 (LOC_Os04g31900) and OsCAMTA6 (LOC_Os07g43030) are involved in transcriptional regulation of miR156a and miR167h.
To find out whether OsCAMTA4 binds to the promoter of miR156a and miR167h, Y1H was performed using their promoter regions as bait. The promoter fragment containing the CAMTA-binding sites was cloned into the pABAi and transformed into Y1HGold to generate reporter strains. The basal expression of the bait reporter strain in the absence of prey was checked and found to be null. On the other hand, the fulllength-coding sequence of OsCAMTA4 and 6 was cloned into pGAD7 to generate a prey vector. Both the prey vector and empty pGAD7 (as negative control) were transformed into Y1H gold strain-harboring promoter segments of miR156a and miR167h in separate experiments. Diploid yeast colonies containing miR156a promoter as bait and OsCAMTA4 and 6 as prey were able to grow strongly on the selective auxotrophic media-containing aureobasidin A. But miR167h shows a strong interaction with OsCAMTA6 and not with OsCAMTA4 ( Figure 6B). Thus, hereby, the physical interaction of the CAMTA-binding sites residing in the promoters of miRNA156a, miR167h, with OsCAMTAs, was confirmed in rice.
Furthermore, expression patterns of OsCAMTA4 and 6 were checked in different tissues under control and stress conditions. Both the CAMTAs respond to inhibitor treatment in seedlings by the reduction in expression; however, their dehydration response could not be confirmed (Figure 7A). In the mature droughtstressed rice plant, both the miRNAs, miR156a, and miR167h-5p are downregulated greatly in Flag Leaf while remaining close to control in a spikelet (Balyan et al., 2017). Under the same drought conditions, both the CAMTAs are downregulated significantly in the Flag Leaf with only OsCAMTA4 reducing in inflorescence ( Figure 7B). Notably, these two miRNAs (miR156a and miR167h-5p) appear to co-regulate with OsCAMTA4 during seedling dehydration stress as well as FL drought stress. Thus, OsCAMTA4 appears to be a major regulator of these miRNAs under different growth and environmental stimuli.

DISCUSSION
These regulatory small molecules named miRNAs are involved in almost all the critical biological processes in plants, including growth and development, as well as combating various stresses. While a significant number of studies have been done to understand the molecular and biochemical processes that are regulated by miRNAs, few studies attempt to study the regulation of miRNA genes themselves. To address this area, we investigated the molecular nature of the schema that regulates miRNA genes by looking into the involvement of one of the most important signaling molecules, i.e., cytosolic levels of calcium in regulating miRNA gene expression. In the past, Kaplan et al., 2006 studied the effect of [Ca 2+ ] cyt bursts induced by calmodulin inhibitors on gene expression of protein-coding genes in Arabidopsis. Over 200 genes were found to be differentially regulated (162 upregulated and 68 downregulated significantly), and it was summarized that [Ca 2+ ] cyt bursts in the cytosol affect the transcriptome to a certain extent and are responsible for the induction of signaling-related genes. In another study, where different types of [Ca 2+ ] cyt bursts were given (single peak, oscillations, and prolonged [Ca 2+ ] cyt elevation) by electrical stimulus, changes in gene expression were observed in Arabidopsis (Whalley et al., 2011). Herein again, the authors found that more genes were upregulated than downregulated in response to [Ca 2+ ] cyt single FIGURE 4 | [Ca 2+ ] cyt mediates the ABA response of miRNAs in N22. The qRT expression profile of miRNAs in seedlings under control conditions, treated with ABA (100 µM) or treated with calcium channel inhibitors, followed by ABA (100 µM), has been shown. The asterisks denote significant change as observed by the t-test. peak burst (104 up vs. 30 downregulated genes) and oscillations (256 upregulated vs. 97 downregulated genes) with a large number of overlapping genes among the upregulated ones. In our study as well, we found that miRNAs, indeed, respond to [Ca 2+ ] cyt levels. The proportion of such miRNAs in N22 (∼2.8%) is pretty close to the estimate of protein-coding genes found by Kaplan et al. in Arabidopsis (∼3%) and by Feske et al. (2001) in human T lymphocyte cells (∼2.1%). Additionally, a higher proportion of miRNAs was upregulated by [Ca 2+ ] cyt than downregulated, which is similar to the findings mentioned above. Among the [Ca 2+ ] cyt -responsive miRNAs as well, a big proportion of miRNAs (such as miR156, miR159, miR166, miR167, miR319, and miR396) is also responsive to, at least, one or multiple abiotic stresses, such as cold, drought, salt, UV, and mechanical .
The data described above identified several miRNAs whose expression is regulated by cytosolic calcium levels. On the other hand, it is also possible that miRNAs target genes involved in the calcium-signaling cascade. Although prediction-based models have specified [Ca 2+ ] cyt -signaling components as miRNA targets for a long time now (Sunkar et al., 2008), our study identifies those targets on the basis of degradome data. Data analysis identified genes such as IQ CaM-binding motif family protein (miR531b), OsCML16 (miR5517), OsCML22 (miR164a, c, d, and e), CAMK_CAMK_like.9 (miR2090), CAMK_CAMK_like.11 (miR5151), and MAPKKK (miR1428e-5p, miR2876-3p). In an earlier study in tomatoes, an autoinhibited [Ca 2+ ] cyt ATPase ACA10 has been shown to be targeted by miR4376. This miRtarget module is involved in proper fruit development since the overexpression of either the miRNA or the miR-resistant ACA10 resulted in elongated stamen filaments and a reduced number of mature fruits . Thus, it appears that miRNAs play an important role in the calcium-mediated regulatory schema of the plant cell.
[Ca 2+ ] cyt is known to be the mediator of several biotic and abiotic stress responses like heat, drought, touch, cold, and salinity (Knight et al., 1996;Albrecht et al., 2003;Choi et al., 2005;Wu et al., 2012;Schulz et al., 2013). In this study, we show that [Ca 2+ ] cyt also mediates the dehydration response of several miRNA genes as well. These included miRNAs like miR156a, miR167h-5p, miR168a-5p, and miR5788. Among these, the former two are also responsive to resting [Ca 2+ ] cyt under nonstress conditions, while the latter two are not. The trend exhibited by the latter group could be speculated by the involvement of some other factors such as dehydration-specific transcription factors that might be controlled by [Ca 2+ ] cyt signaling. Earlier, H 2 S, which is a potent-signaling messenger in various plant physiological processes, including drought (Jin et al., 2011, has been shown to mediate drought response of some miRNAs like miR167a/c/d, miR393a, miR396a, and miR398a/b/c in Arabidopsis (Shen et al., 2013). Another significant hormone known to be involved in many abiotic stresses, especially drought/dehydration stress-the abscisic acid and its influence on miRNA expression-has also been discussed in this study. In general, [Ca 2+ ] cyt is known to mediate many ABA responses, such a stomatal closure and seedling growth (McAinsh et al., 1990;Guo et al., 2002), and, thus, there is a significant overlap of ABA and calcium-mediated signaling. Indeed, the expression of several miRNAs was found to be regulated by ABA. Expression trends of miR1425-5p, miR156a, miR162b, miR168a, miR319b, miR528-5p, and miR530-5p match with previously known knowledge, while the rest had no previous information on their ABA response (Shen et al., 2010;Tian et al., 2015). Further analysis conducted to decipher whether modulation in the cytosolic calcium levels mediates ABA inducibility of the miRNAs resulted in the finding that cytosolic calcium levels do, indeed, mediate the ABA response of miR159b, miR319b, and miR530-5p.
Further investigation into the signaling cascade regulating the miRNA expression revealed the involvement of calmodulins (miR156a, miR1878, miR396c-3p, miR166g-3p, miR167h-5p, miR1425-5p, and miR1876). Calmodulins are signal sensor-relay proteins that directly bind to calcium and relay the signal to other proteins, such as CAMTAs (calmodulin-binding transcription activators), that bind to particular motifs in promoter sequences of genes and either activate or suppress transcription. The use of calmodulin inhibitors and camta knockdown mutants affected the expression of certain miRNAs, namely, miR156a, miR160a-5p, miR166a-3p, miR167h-5p, and miR168a-5p. The effect of knockdown of camtas on the miRNA expression in the mutants can be a result of direct interaction between the TF and the miRNA promoter or it can be routed via some other TFs/proteins that might be regulated by the CAMTAs. CAMTAs have been previously shown to mediate cold response as well as salicylic acid-mediated immunity in Arabidopsis. However, their involvement in miRNA regulation is not wellknown in plants. Besides, several calcium-responsive promoter motifs have been identified in the promoter sequences of calcium-responsive genes (Finkler et al., 2006;Kaplan et al., 2006;Galon et al., 2010). Accordingly, we found the presence of these calcium-responsive motifs, namely, CARE (calciumresponsive elements), ABRE (ABA-responsive elements and related motifs), E-box, G-box, GT-box, Z-box, and CAMTAbinding sites in the promoter sequences of calcium-responsive miRNAs. The Y1H shows strong physical interaction between OsCAMTA4 and miR156a, but this CAMTA does not interact with miR167h. However, OsCAMTA6 shows strong interaction with the promoters of miR156a and miR167h. Thus, hereby, we reveal novel and key players in the transcriptional regulation of miR156a and miR167h. While miR156a is seen to be responsive to [Ca 2+ ] cyt in resting conditions, positively by calmodulin as well as [Ca 2+ ]-mediated dehydration, miR167h-5p has a slightly different story. Although [Ca 2+ ] cyt affects its expression under control conditions as well as dehydration, miR167h-5p shows negative regulation by calmodulins. Both of these miRNAs show binding and trans-activation by OsCAMTA4, which appears to be co-regulated in seedling dehydration as well as FL drought stress. This appears as a significant finding since miRNA expression has to date not been associated with CAMTAs. miR156 has been studied in quite a detail, and there is plenty of knowledge available about its function and regulation. The miRNA is a floral repressor, and promoter of the juvenile phase in Arabidopsis (Wu G. et al., 2009) is involved in lateral root development and leaf morphology in Arabidopsis (Gao et al., 2018). Regarding its regulation, it has been shown that AGL15 and AGL18 act in cooperation to promote its transcription by binding to the CArG motifs in its promoter (Serivichyaswat et al., 2015). The phytochrome-interacting factors or PIF1, PIF3, PIF4, and PIF5 have also been shown to directly bind and repress the expression of miR156b/d/e/f/h to enhance the shade avoidance syndrome in Arabidopsis. Another protein called DOG1 (DELAY OF GERMINATION 1) plays a role in the efficient processing of primary miR156 to its active mature form by regulating the processing proteins-DCL1, HYL1, SE, TGH, and CDC5 (Huo et al., 2016). Recently, another player has been found in miR159 that targets MYB33 that, in turn, binds and promotes the transcription of miR156a&c during the young seedling stage (Guo et al., 2017). Furthermore, clues behind its temporal expression pattern have been found in epigenetic regulation. The transcription-activating mark H3K4me3 is seen abundantly at the miR156a and miR156c loci during the early seedling stages wherein it contributes to its high expression (Xu et al., 2018), while, during the vegetative phase change, there is an increase in the level of histone H3K27me3 with a simultaneous decrease in H3K4me4 and H3K27ac at regions upstream and immediately downstream of its TSS, resulting in its decline in abundance (Xu et al., 2016). Besides, a cycling DOF transcription factor CDF2 has also been shown to be its transcriptional activator (and repressor of miR172) acting in the same signaling pathway to control flowering in Arabidopsis (Sun et al., 2015). Regarding the regulation of miR156 in rice, drought has already been shown to downregulate the miRNA in inflorescence tissue . Besides, it is also responsive to several hormones, such as auxin in Arabidopsis [downregulation; (Marin et al., 2010)] and ethylene in tomato [downregulation; (Zuo et al., 2012)]. In rice, the miRNA has been shown to be independent of any regulation by gibberellin during juvenile to the adult phase transition (Tanaka, 2012). Our data bring more regulators of miR156 into light, which are [Ca 2+ ] and OsCAMTA4 and 6. Since these are located pretty high in the signaling hierarchy that generates any responses inside a cell, this new mode of regulation might help to explain the various responses miR156 displays under the various abiotic stresses and different stages of the lifecycle of the plant. miR167 is also a conserved miRNA across monocots and dicots and is known to target the ARF6/8 genes, thereby acting as a node in auxin signaling (Barik et al., 2015). It is known to be involved in adventitious root development in rice (Meng et al., 2011), modulation of auxin signaling during bacterial infection in tomato (Jodder et al., 2017), in regenerating calli in rice (Sinha et al., 2019), in blue-light signaling in Arabidopsis (Pashkovskiy et al., 2016), and abiotic stress, such as salinity (Liu et al., 2008;Ding et al., 2009;Frazier et al., 2011) and drought in rice (Balyan et al., 2017). For such an evolutionarily conserved and functionally significant miRNA, the information about its regulation was lacking. In our study, we showed its regulation via [Ca 2+ ] cyt , calmodulin, and OsCAMTA6. Its dehydration response is also mediated via [Ca] 2+ cyt . Thus, it is a remarkable discovery for a miRNA that acts as a node in a critical signaling pathway such as the auxin signaling, ultimately regulating a plethora of plant functions.

CONCLUSIONS
The understanding of the regulatory pathways governing the expression of another class of regulators, i.e., the miRNAs is critical in order to be able to manipulate them for advantageous traits in plants. With this view, the study explored the involvement of cytosolic calciumin-regulating expression of miRNA genes under control and drought stress conditions that we could establish with the help of calcium channel inhibitors and ionophores. The fact that [Ca] 2+ cyt mediates several abiotic responses was also explored and it was demonstrated, it is the same for miRNAs in rice as well. The further dissection revealed the involvement of calmodulins and CAMTAs. Through yeast-onehybrid experiments, OsCAMTA4 and 6 were proved to bind the CAMTA-binding sites of the very critical miR156a and miR167h. Both these CAMTAs were found to be coregulated with these two miRNAs at various developmental and stress stages indicative of a possible regulatory schema for these miRNAs in rice.

DATA AVAILABILITY STATEMENT
The data presented in the study are deposited in NCBI under the accession number PRJEB47136.

AUTHOR CONTRIBUTIONS
SR conceived the concept. SR and SK designed the experiments and prepared the manuscript. SK constructed the NGS libraries and performed the data analysis and qRT-PCR. VP performed qRT-PCR and yeast one-hybrid experiments. RM performed the degradome analysis. All the authors read and approved the final manuscript.

FUNDING
The study was funded by Science and Engineering Research Board, Department of Science and Technology under the Grant No: EMR/2016/006081 awarded to the corresponding author.

ACKNOWLEDGMENTS
SK is grateful to CSIR for the fellowship, and RM and VP thank Delhi University for the fellowships during the research tenure.