Abstract
In plants, arsenic (As) stress modulates metabolic cascades at various developmental stages by influencing the pattern of gene expressions mediated by small non-coding RNAs, especially Micro-RNAs, involved in the moderation of a myriad of cellular processes needed for plant adaptation upon oxidative stress. miR156j of miR156 gene family, involved mainly in the regulation of growth and development in plants. This study was designed to find out the role of arsenic toxicity on Osa-miR156j expression in all physiological growth stages. To better understand the functional role of Osa-miR156j in rice, we observed the expression in different developmental stages (seedlings, tillering and flowering) and various tissues of leaf, stem and root tissues (at 0, 24, 48, and 72 h) under 25 μM arsenite [As (III)] exposure. Additionally, using bioinformatic tools to target genes of Osa-miR156j and the potential co-expressed genes were explored at different development stages in the various tissues of rice under stress conditions. The expression of Osa-miR156j showed its temporal downregulation in various tissues in different developmental stages. Of note, the downregulation was more pronounced in root tissues at seedlings, tillering, and flowering stages during 0–72 h under arsenite exposure as compared to other tissues. Overall, the As stress altered the gene expression more prominently at seedlings developmental stage followed by flowering and tillering. Additionally, through the In silico approach, the target functions and presence of oxidative stress-responsive cis-acting regulatory elements/motifs also confirmed Osa-miR156j involvement in the regulation of arsenic stress in rice. The findings of this study demonstrate the prominent role of Osa-miR156j in rice under arsenite stress, which was found to modulate the metabolic activities in rice plants at different developmental stages, and thus it might be useful for the development of arsenic tolerant varieties.
Introduction
Arsenic (As) is a metalloid, which is ubiquitously distributed in soils, most plants species, and aquifers (Zhang and Selim, 2006). This metalloid can cause some severe threats to human health primarily through biomagnification in the food chain due to contaminated crops (). It’s been estimated that an average natural arsenic content level is about 5 mg/kg in the soils worldwide, which is due to leaching from sedimentary rocks, mining and smelting of the As bearing ores and minerals as well as the coal mining (; ; ). However, due to the indiscriminate use of arsenic in herbicides, insecticides, and wood preservatives, etc., its concentration is increasing at an alarming rate in the soil of South-East Asian sub-continent (). Rice (Oryza sativa) is the staple food and hence, principal crop for over half of the world. The presence of high arsenic content in the soil has led to its accumulation in paddy, which has resulted in its collection at the concentration beyond the food safety threshold (Zhu et al., 2008; Zheng et al., 2011). Arsenic accumulation in plants culminates in both direct and indirect toxicity (). Some studies on arsenic toxicity on rice seedlings showed compromised photosynthetic rate (; ), disturbed carbohydrate metabolism (), subdued nitrogen assimilation (), elevated phytochelatin synthesis (), overproduction of reactive oxygen species (ROS), marked increase in lipid peroxidation in it leading to the oxidative stress (; ; ). However, plants have an array of mechanisms to detoxify arsenic poisoning, which mainly include metal transport, chelation, and sequestration through organoarsenic compounds (; ; ; ).
MicroRNAs (miRNAs) are a discrete class of short-chain (21–24 nucleotides) small RNAs that play a central role in modulating gene expression at post-transcriptional levels (). The miRNAs have been reported to play an essential role in plant tolerance to abiotic stresses such as drought, salinity, cold, heavy metal toxicity, and nutrient deprivation (). A variety of miRNAs have been identified and characterized for various metabolic and biological processes, including tissue maturation (), signal transduction (), hormone signaling (), the phase change from vegetative to reproductive growth (), and response to environmental stresses (). The plant miR156 has been reported to target squamosa promoter binding protein-like (SPL) transcription factor family (), which are expressed more at the younger phase whose expression gradually declines with the plant age (). In Arabidopsis, miR156 plays a critical role during different developmental stages such as leaf development, shoot maturation, phase change, and flowering by targeting 11 SPLs out of 17 SPL genes. The over-expression of miR156 in transgenic Arabidopsis showed reduced apical dominance, delayed flowering time, caused dwarfism, and increased total leaf numbers and biomass (). Out of 19 SPLs, Osa-miR156 targets 11 SPLs, which influence various morphological changes during the development of rice (Xie et al., 2006). Some reports suggest its crucial roles during various abiotic stresses (; ). A recent report from Yu et al. (2012) pointed out the inevitable role of Osa-miR156j, a member of the Osa-miR156 family, due to its significant down-regulation in rice seedlings during arsenic-stress. This down-regulation suggests some important regulatory functions associated with Osa-miR156j during arsenic stress at developmental stages in rice. Therefore, this study was planned to extricate the effect of arsenic stress on the expression of Osa-miR156j at different developmental stages in different tissues for better understanding of miRNA functional genomics and the molecular role of involved genes against As-induced stress tolerance in rice.
Materials and Methods
Rice Growth and Arsenic Treatment
Rice (Oryza sativa L.) cultivar IR-64 seeds were germinated in a soil-less medium. Seeds were surface-sterilized with 3% sodium hypochlorite solution for 20 min and washed with distilled water. The sterilized seeds were soaked in distilled water for 24 h in a flask at 37°C. Soaked seeds were kept on moistened filter paper in Petri dish and incubated at 30°C under the dark condition for three days for germination. After germination, the coleoptile of germinated seeds was slotted in the styrofoam and placed in Yoshida nutrient solution (YS) (Yoshida et al., 1971) containing disposable 250 ml plastic cups. For the natural growth of the seedlings, the hydroponic system was transferred in glasshouse conditions with day/night temperatures of 30/20°C and relative humidity of at least 50% during the day. The 14 days old seedlings were treated with 25 μM arsenite (NaAsO2). For the treatments of arsenite in tillering and flowering stages, 14 days old fresh seedlings were subsequently transferred to soil-filled pots containing balanced nutrients as per requirement. After attaining the tillering stage (27 days) and flowering stage (60 days), plants were harvested with roots. The roots were washed with tap water. For the treatment of arsenite, plants were subjected to hydroponic media containing 25 μM arsenite. At all the developmental stages, plant samples were collected at 0, 24, 48, and 72 h from leaf, stem, and roots. The samples were frozen immediately in liquid nitrogen and stored at -80°C until used for further analysis.
Targets and cis-Elements/Motif Prediction
To find out targets of Osa-miR156j, the publically available tool psRNATarget1 was used by searching cDNA OSA1 release five and transcript, RAP-DB, version 1.0 as reference genome (). To identify the molecular functions of the searched targets, we used the AgriGO2 online analysis tool (). The upstream sequence 200-bp was analyzed by PlantCARE3 to search for the presence of cis-regulatory elements/motifs ().
In silico Co-expression Network and MPSS Analysis of Target Genes
The Gene co-expression network of miRNA-target genes and co-expressed genes were obtained from the publically available tool RiceFrend4, and networks were constructed using Cytoscape 3.2.1 tool (; ). Massively Parallel Signature Sequencing (MPSS) data were obtained from the rice MPSS database for different development stages and various tissues of Rice. Collected data were further analyzed by the Clustvis tool5.
RNA Extraction, Polyadenylation, and cDNA Synthesis
Total RNA in collected samples was extracted by using TRIzol® reagent (Invitrogen) according to the manufacturer’s instructions. The quality of extracted RNA was visually examined by the 28S/18S ribosomal bands on an agarose gel, and the quantity was checked by spectrophotometer. Total RNA and Short length of miRNA were treated with E. coli poly A polymerase (NEB, United Kingdom) to generate a poly-A tail at the 3′ end of each RNA molecule. The polyadenylation reaction with a final volume of 20 μl includes 250 ng of total RNA, 4.0 μl of 5 × poly-A polymerase buffer, 1.0 μl of 10 mM rATP and 1 μl of E. coli poly A polymerase (New England Biolabs, United States). The reaction was incubated at 37°C for 30 min and adenylation terminated by incubating reaction at 95°C for 5 min. The1st-strand cDNA synthesis was performed using a miRNA 1st-strand cDNA synthesis kit (Stratagene) as per manufacturer’s instructions.
Primers Design and qRT-PCR Analysis
The DNA primers used in real-time analysis of Osa-miR156j were designed as previously described by . Reference gene, β-Actin (Wu et al., 2009; ) primers were designed using Primer Express®Software v2.0 (Applied Biosystems, United States) under default parameters. The primer sequences used in this study are given in Supplementary Table S1. For quantification of miRNA, final reaction, 25 μl was comprised 2 μl of cDNA (10 ng), 12.5 μl of 2× SYBR advantage premix, 0.5 μl of 50× ROX dye (Clontech, United States), 0.5 μl of 10 μM miRNA specific primers. A negative control without a template, the reaction was also performed with the same primers. Three biological replicates for each sample were used for RT-PCR analysis. The real-time PCR was performed, employing a 7300 Real-Time PCR System (Applied Biosystems, United States). PCR cycles 95°C for 10 s, followed by 35 cycles of 95°C for 10 s and 60°C for 20 s, were performed in 96-well optical reaction plates (Applied Biosystems). The specificity of the amplicon was observed by analysis of the melting curve after 35 cycles at 60–90°C. The abundance/decline of Osa-miR156j was normalized to a geometric average of endogenous control (β-actin) for ΔCt and expressed as relative Osa-miR156j expression.
Northern Blotting Analysis
Total RNA from different samples was isolated with Trizol reagent. The total RNA, 25 μg, was resolved in the per lane on a denaturing 15% polyacrylamide gel with 7% urea at 40 mA for 2 h. The gel was stained with ethidium bromide, and 5S RNA/tRNA bands were used as loading controls. RNAs were electro-blotted on Hybond N + membranes (Amersham Biosciences, NJ) using a Trans-Blot Cell (Bio-Rad, CA, United States). Membranes were UV cross-linked and baked at 80°C for 1 h. DNA oligonucleotide probe complementary to miRNA was labeled with ATP γ-P32 using T4 polynucleotide kinase (NEB, United Kingdom). The blot was pre-hybridized for 30 min at 40°C in pre-hybridization buffer containing 7% SDS, 200 mM Na2HPO4 (pH 7.0), 5 μg/ml salmon sperm DNA. After removing the pre-hybridization buffer, the blot was hybridized with a 50 pmol/ml probe containing a hybridization buffer for 16 h at 40°C. The hybridized blot was subsequently washed three times with 1X SSC, 0.1% SDS buffer at room temperature. Finally, the blot was exposed to a phosphor screen and cross-linked with 1200 μ joules for 20 min to improve sensitivity. The bands were analyzed using the Alpha Imager documentation system. The data were expressed as mean ± standard deviation (SD) of the band density of the experiments.
Statistical Analysis
All the experiments were performed in triplicate and repeated three times, and the statistical analysis between the different groups was performed by Student t-test using GraphPad Prism 7.0 software. The data were represented as mean ± SD with p-value < 0.05 was considered significant.
Results
Identification of Potential Targets and cis-Regulatory Elements/Motifs in Osa-miR156j and Their Co-expression Gene Networks and MPSS
The potential targets of Osa-miR156j (locus IDs), mainly comprised of DNA binding transcription factors and DNA binding proteins, acetyltransferases, phospholipase A1, teosinte glume architecture 1, etc., as shown in Figure 1 and Supplementary Table S2. Furthermore, Osa-miR156j co-expression target networks are shown in Figure 2 and Supplementary Table S3. To elucidate cis-acting regulatory elements linked with the regulation of diverse metabolic and physiological functions, the immature sequence of Osa-miR156j was analyzed using the bioinformatics tool, PlantCARE. The potential cis-acting regulatory elements in the promoter region of miRNA were comprised with the presence of many common abiotic stress regulatory motifs/elements like ABRE (abscisic acid responsiveness), CCAAT-box (MYB binding site), GC-motif (enhancer like elements involved in anoxic specific inducibility), HSE (heat stress-responsive elements), MSA (involved in cell cycle regulation response) and TCA (salicylic acid responsiveness).
FIGURE 1
FIGURE 2
Using MPSS data, Heat map generated using the Clustvis tool shows in silico expression level and classification of identified genes at different development stages, various tissues under biotic and abiotic stress in rice (Figure 3 and Supplementary Table S4). Further, the expression values obtained from Rice MPSS for 17 and 20 bp sequences were used for analysis and representation. Red color designates the highest level of expression frequency at different development stages and various tissues followed by other colors for their corresponding level of expression in the heat map.
FIGURE 3
Expression Analysis by RT-PCR and Northern Blotting in Various Plant Tissues
To investigate the role of Osa-miR156j under As stress at various developmental stages, the expression patterns were studied in different tissues of seedlings, tillering, and flowering physiological stages in rice during 0–72 h. The results show that the expression of Osa-miR156j at the seedling stage in the stem was a nearly similar or slight decline in expression at 48 and 72 h exposure (Figure 4). Up to 24 h of As (III) treatment, there was no change in expression, while 1.61 and two-fold decline in expression of Osa-miR156j was observed in 48 and 72 h exposure. In the leaf of the seedling stage, a 1.8-fold decline in expression at 24 h, followed by a 2.31 and 4.05-fold decline in expression of Osa-miR156j at 48 and 72 h exposure was recorded. The expression of Osa-miR156j at the seedling stage in root followed the same pattern as the leaf. The 2.71, 4.52, and 11.17-fold decline in expression was observed at 24, 48, and 72 h As (III) exposure, respectively. At the tillering stage, the expression of the Osa-miR156j in leaf declined only by 1.26-fold, followed by 1.9-fold at 24 and 48 h exposed plants. The 4.22-fold decrease in Osa-miR156j expression was recorded after 72 h As (III) exposure. But at the tillering phase, a sharp decline in expression Osa-miR156j was observed, i.e., 8.05, followed by 40.25-fold in 48 and 72 h exposed plant stems, respectively. In the root of tillering phase 3, 9, and 45-fold decline in expression of Osa-miR156j was recorded at 24, 48, and 72 h exposed to As (III), respectively. While at flowering stages, the sharp decline in expression of the Osa-miR156j gene was recorded with an increase in As (III) exposure time. The expression of Osa-miR156j at the flowering stage in leaf tissue is lower in comparison to the stem and root cells of plants. The expression of the Osa-miR156j in leaf declined by 3.08, followed by 7.70 and 22-fold at 24, 48, and 72 h As (III) exposure, respectively. Similar results were obtained in stem and root cells as nearly 5.5-fold decline in expression of Osa-miR156j exposed to 24 h, while 16-fold decline in stem and root cells exposed to 48 h As (III) exposure.
FIGURE 4
For further validation of the expression results obtained, we used Northern blotting analysis to examine expression levels. Northern blotting analysis gave a similar analysis pattern in agreement with the RT-PCR results (Figure 5). At 0th h in the leaf, stem, and root at all developmental stages under arsenic stress, we observed Osa-miR156j significant (p < 0.05) decreased band density. This has also been seen in 24, 48, and 72 h, with a constant significant decrease in the band density levels of Osa-miR156j under As (III) with 72 h showing the least band density in the seedlings, tillering, and flowering stages. The duration of As (III) exposure was also found to affect the gene expression such that with an increase in exposure time, the downregulation of the expression of Osa-miR156j was more pronounced.
FIGURE 5
Discussion
Oryza sativa is the primary cereal crop, and its productivity is compromised under several abiotic stresses such as heavy metals, heat, salt, and drought, etc. The distribution of metalloid arsenic is ubiquitous in the environment, which is reported to cause severe stress in rice. Studies by and have shown the accumulation of As at maximum levels in roots as compared to leaves and stems, as the dosage increases As accumulation increases in a similar manner (; ). These accumulations occur maximum in the roots followed by stem than in the leaves and least in the economic production. Besides, it is said to get accumulated in the edible parts, i.e., grains (2 mg/kg) and present a noteworthy illustration of biomagnification by finally reaching the human food chain where they exhibit several deleterious effects (; ). Arsenic in the form of arsenite, i.e., As (III) enters plant by the roots through nodulin 26-like intrinsic protein (aquaporins) (; ). Studies in rice plants have shown that aquaporins (OsNIP1;1, OsNIP3;1, OsNIP3;2, OsNIP3;3, OsNIP2;1, OsNIP2;2, OsLsi2, OsPIP2;4, OsPIP2;6 etc.) help in the uptake, transport, and influx of As (III) (; ; ; ; ).
Albeit our substantial progress in understanding heavy metal(s) responses in plants, the molecular parameters responsible for physiological responses are still not well identified. Under metal toxicity, plants try to modify gene expression through miRNAs to regulate excess metals accumulation by various means such as complexation, chelation, antioxidant defense against ROS, and control multiple biological responses through signal transduction (). In response to As stress, plants modulate metabolic cascades at various developmental stages through alteration in gene expression. Therefore, it becomes prudent to understand the As-induced abiotic stress at the molecular level so that we may develop a strategy to curb this menace, or we can develop As-resistant rice variety. Plant miRNAs are involved in various stress responses that regulate the gene expression by negatively regulating the complementary target genes.
There are different microRNAs reported in plants with various functions, in which miR156 has significant agricultural importance for biotic stress tolerance and plant development (). The genomic investigation of miR156 deciphered that miR156 regulatory repertoire is highly conserved during plant evolution (; ; Zhang et al., 2015). ; Xing et al., 2010 have reported that miR156 downregulates Squamosa Promoter Binding Protein Like genes, which enrich growth from juvenile to mature, leaf formation, flowering, and vegetative phase (; Xing et al., 2010). reported that all the conserved miRNAs, such as miR156, miR171, miR172, miR395, miR397, and miR398, were down-regulated in grapevine under abiotic stress (). Yu et al., 2012 also reported that miR156, among other 12 miRNA families, plays a key role during the As stress in rice (Yu et al., 2012). also deciphered the non-coding RNA-level response to arsenic stress in rice (Oryza sativa) and enlarge the present molecular understanding of As stress response, particularly at the non-coding RNA level (). The miR156j, a member of the miR156 gene family, is highly conserved and helps in regulating the growth and development in plants. This study emphases the response of Osa-miR156j in rice cultivars exposed to As stress. Rice is an extensively studied model system in plant biotechnology; yet, the response of these miRNAs against As stress at the different growth phases of rice plant (seedling, tillering and flowering) on root leaf and stem separately has not been studied until our report. The scope of this study could be utilized for a better understanding of molecular mechanisms of arsenic stress and the roles of miRNAs. In the present study, we have investigated the expression pattern of Osa-miR156j and their associated target networks in different developmental stages in various tissues of rice during As stress.
Osa-miR156j is one of the more predominant plant microRNA families predicted to target SPL transcription factors to control developmental timing, and the phase change from vegetative to reproductive growth (). The candidates of this family are expressed mainly in the juvenile stage and gradually decrease with the advancement of the plant age (). Identification of micro RNAs and their targets in switchgrass, a model biofuel plant species, target analysis of Osa-miR156j showed that its cleavage targets are SPL9, SPL10, and SPL12 genes. During the growth phase, SPL9 acts as a transcriptional repressor (). Genes SPL10 and SPL11 play a central role in the development of lateral organs in association with shoot maturation at the reproduction stage, but SPL10 is the only transcription factor that controls shoot development regulation during vegetative growth (). Although, the expression of Osa-miR156j has been studied previously at seedlings stage in some plants including rice under arsenic stress (; Yu et al., 2012) however, its expression profile (and hence the role) is entirely unknown at other developmental stages in different plant tissues at different time intervals. Taken together, our data suggest that the regulation of Osa-miR156j expression is highly influenced by SPLs transcription factors, acetyltransferases/phospholipase A1 enzymes, and its down-regulation might be results of activation/deactivation of these transcription factors/enzymes. In previous studies, these target genes have been reported for their involvement in plant growth and development during oxidative stresses (; ). Acetyltransferases play a vital role in histone modification by acetylation during regulation of gene expression by DNA replication, DNA repair, DNA recombination, and gene transcription (; ). Histone acetylation by acetyltransferases (HATs) activates gene expression, whereas de-acetylation leads to gene repression (, ). In rice, HATs (OsHATs) are expressed constitutively, and their expression regulated by hormones and oxidative stresses (). In Arabidopsis, HATs (AtHATs) have been reported to regulate light-inducible gene expression (), root stem-cell niche maintenance () and regulation of miRNAs accumulation at transcriptional and posttranscriptional levels (). Phospholipases are phospholipid hydrolyzing enzymes, which activate plant immune responses associated with the development of hypersensitive cell death (). These enzymes are known for the production of outstanding defense signaling molecules like oxylipins, jasmonates, and phosphatidic acid in plants (). In plants, members of the phospholipases super-family, phospholipase A1 plays an essential role in maintenance and remodeling cell membrane and regulates various cellular mechanisms by the production of various lysophospholipid mediators, such as lysophosphatidylserine and lysophosphatidic acid, which involve in multiple biological functions ().
In recent years, independent research on receptor-like protein kinases (RLKs) has shown its role in biotic and abiotic stress response and its resistance besides signal transduction in plants, which can further promote agriculture development (Ye et al., 2017). Our results through In silico networking have shown that Osa-miR156j is associated with eight potential genes. These identified potential genes (Os08g0509600, Os02g0174100, Os02g0139400, Os06g0663500, Os09g0513100, Os09g0507100,and Os01g0922600) were associated with zinc finger proteins, wall-associated receptor kinase, lectin receptor-like kinase, UDP-glycosyltransferase and lipid metabolic process, etc. through co-expression networks, which were analogous to the previous studies (Zhou et al., 2011; Yu et al., 2012; ; ; ; ). The presence of cis-regulatory elements in the promoter region indicated that Osa-miR156j might have an essential role in the regulation of transcription factors, metabolisms, stress responses, signaling transductions, cellular structural components, and other cellular processes. ABRE cis-regulatory elements play a vital role in regulating oxidative stresses in rice (Yamaguchi-Shinozaki and Shinozaki, 2006). Promoter elements ABREs are known to be controlled by phytohormone abscisic acid (ABA) and salicylic acid signaling (). They contribute to a range of developmental and adaptive processes to environmental stimuli in plants (). With an increase in endogenous ABA levels in response to osmotic stresses, ABA activates the expression of many genes through ABREs in their promoter regions (). HSEs induce oxidative stress genes and play an essential role during signal transduction pathways mediating activation of high temperatures and other abiotic stresses cold and heavy metals (Yi and Liu, 2009; ). MSAs are M specific activators, and they control the M- phase-specific gene expression in plants (). Hence, the presence of these cis-acting elements in the promoter region mainly suggested that these Osa-miR156j played a critical role in the regulation of transcription factors, metabolisms, stress responses, signaling transductions, cellular structural components, and other cellular processes.
Conclusion
Under arsenic abiotic stress, the expression analysis of Osa-miR156j has successfully indicated Osa-miR156j importance at various developmental stages and tissues in rice. Additionally, the target functions and presence of oxidative stress-responsive cis-acting regulatory elements/motifs also confirmed the miR156j involvement in the regulation of arsenic abiotic stress, which is highly influenced by exposure duration and plant tissues.
Statements
Data availability statement
All datasets generated for this study are included in the article/Supplementary Material.
Author contributions
AV conceived grant-in-aid, designed, and supervised the complete study and supervised the experimental procedure and proofread the entire manuscript. MG and AP performed all the wet and dry lab experiments and statistical analysis besides draft preparation. All authors contributed to the manuscript and approved the submitted version.
Funding
A grant-in-aid supported this work under Project No. SERC/LS-504/2012 by Department of Science and Technology, Government of India to AV. AV is sincerely thankful to DST, New Delhi, for providing financial support.
Acknowledgments
AP, MG, and AV are sincerely thankful to Banaras Hindu University and U.P. Council of Sugarcane Research for providing facilities to perform the experiments.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.00752/full#supplementary-material
TABLE S1Primers and probes used in the study.
TABLE S2Osa-miR156j potential target genes, locus IDs, and their target functions during arsenic abiotic stress.
TABLE S3Osa-miR156j target gene; Os01g0922600, Os06g0663500, Os07g0505200Os09g0507100, Os02g0174100, Os09g0513100, and their potential co-expressed target genes, functions in the co-expression network.
TABLE S4Indicates the data obtained from the rice Massively Parallel Signature Sequencing database with both 17 and 20 bp sequences for different development stages and various tissues in rice.
Footnotes
1.^http://plantgrn.noble.org/psRNATarget/
2.^http://bioinfo.cau.edu.cn/agriGO/
3.^http://bioinformatics.psb.ugent.be/webtools/plantcare/
References
1
AxtellM. J.SnyderJ. A.BartelD. P. (2007). Common functions for diverse small RNAs of land plants.Plant Cell191750–1769. 10.1105/tpc.107.051706
2
BalcellsI.CireraS.BuskP. K. (2011). Specific and sensitive quantitative RT-PCR of miRNAs with DNA primers.BMC Biotechnol.11:70. 10.1186/1472-6750-11-70
3
BartelD. P. (2004). MicroRNAs: genomics, biogenesis, mechanism, and function.Cell116281–297.
4
BellandeK.BonoJ.-J.SavelliB.JametE.CanutH. (2017). Plant lectins and lectin receptor-like kinases: how do they sense the outside?Int. J. Mol. Sci.18:1164. 10.3390/ijms18061164
5
BienertG. P.ThorsenM.SchüsslerM. D.NilssonH. R.WagnerA.TamásM. J.et al (2008). A subgroup of plant aquaporins facilitate the bi-directional diffusion of As (OH) 3 and Sb (OH) 3 across membranes.BMC Biol.6:26. 10.1186/1741-7007-6-26
6
CanonneJ.Froidure-NicolasS.RivasS. (2011). Phospholipases in action during plant defense signaling.Plant Signal. Behav.613–18. 10.4161/psb.6.1.14037
7
ChenY.HanY.-H.CaoY.ZhuY.-G.RathinasabapathiB.MaL. Q. (2017). Arsenic transport in rice and biological solutions to reduce arsenic risk from rice.Front. Plant Sci.8:268. 10.3389/fpls.2017.00268
8
ChenZ. J.TianL. (2007). Roles of dynamic and reversible histone acetylation in plant development and polyploidy.Biochim. Biophys. Acta Gene Struct. Express.1769295–307. 10.1016/j.bbaexp.2007.04.007
9
ChengX.WuY.GuoJ.DuB.ChenR.ZhuL.et al (2013). A rice lectin receptor-like kinase that is involved in innate immune responses also contributes to seed germination.Plant J.76687–698. 10.1111/tpj.12328
10
ChoudhuryB.ChowdhuryS.BiswasA. K. (2011). Regulation of growth and metabolism in rice (Oryza sativa L.) by arsenic and its possible reversal by phosphate.J. Plant Interact.615–24. 10.1080/17429140903487552
11
ChowdhuryN. R.DasR.JoardarM.GhoshS.BhowmickS.RoychowdhuryT. (2018). Arsenic accumulation in paddy plants at different phases of pre-monsoon cultivation.Chemosphere210987–997. 10.1016/j.chemosphere.2018.07.041
12
DaiX.ZhaoP. X. (2011). psRNATarget: a plant small RNA target analysis server.Nucleic Acids Res.39(Suppl._2) W155–W159.
13
DasI.GhoshK.DasD.SanyalS. (2013). Assessment of arsenic toxicity in rice plants in areas of West Bengal.Chem. Spec. Bioavailabil.25201–208. 10.3184/095422913x13785717162124
14
DuZ.ZhouX.LingY.ZhangZ.SuZ. (2010). agriGO: a GO analysis toolkit for the agricultural community.Nucleic Acids Res.38(Suppl._2) W64–W70.
15
EckardtN. A. (2005). MicroRNAs regulate auxin homeostasis and plant development.Am. Soc. Plant Biol.171335–1338. 10.1105/tpc.105.033159
16
FrancesconiK.VisoottivisethP.SridokchanW.GoesslerW. (2002). Arsenic species in an arsenic hyperaccumulating fern, Pityrogramma calomelanos: a potential phytoremediator of arsenic-contaminated soils.Sci. Total Environ.28427–35. 10.1016/s0048-9697(01)00854-3
17
FuS.-F.ChenP.-Y.NguyenQ. T. T.HuangL.-Y.ZengG.-R.HuangT.-L.et al (2014). Transcriptome profiling of genes and pathways associated with arsenic toxicity and tolerance in Arabidopsis.BMC Plant Biol.14:94. 10.1186/1471-2229-14-94
18
FujitaY.FujitaM.ShinozakiK.Yamaguchi-ShinozakiK. (2011). ABA-mediated transcriptional regulation in response to osmotic stress in plants.J. Plant Res.124509–525. 10.1007/s10265-011-0412-3
19
FujitaY.YoshidaT.Yamaguchi-ShinozakiK. (2013). Pivotal role of the AREB/ABF-SnRK2 pathway in ABRE-mediated transcription in response to osmotic stress in plants.Physiol. Plant.14715–27. 10.1111/j.1399-3054.2012.01635.x
20
GautamA.PandeyA. K.DubeyR. S. (2019). Effect of arsenic toxicity on photosynthesis, oxidative stress and alleviation of toxicitywith herbal extracts in growing rice seedlings.Indian J. Agric. Biochem.32143–148.
21
GautamA.PandeyA. K.DubeyR. S. (2020). Azadirachta indica and Ocimum sanctum leaf extracts alleviate arsenic toxicity by reducing arsenic uptake and improving antioxidant system in rice seedlings.Physiol. Mol. Biol. Plants2663–81. 10.1007/s12298-019-00730-z
22
GautamA.PandeyP.PandeyA. K. (2020). “Proteomics in relation to abiotic stress tolerance in plants,” in Plant Life under Changing Environment, Responses and Management, edsTripathiD. K.SinghV. P.ChauhanD. K.SharmaS. K.PrasadS. M.DubeyN. K.RamawatN. (Elsevier: Academic Press), 513–541. 10.1016/B978-0-12-818204-8.00023-0
23
GielenH.RemansT.VangronsveldJ.CuypersA. (2012). MicroRNAs in metal stress: specific roles or secondary responses?Int. J. Mol. Sci.1315826–15847. 10.3390/ijms131215826
24
GouJ.-Y.FelippesF. F.LiuC.-J.WeigelD.WangJ.-W. (2011). Negative regulation of anthocyanin biosynthesis in Arabidopsis by a miR156-targeted SPL transcription factor.Plant Cell231512–1522. 10.1105/tpc.111.084525
25
HuijserP.SchmidM. (2011). The control of developmental phase transitions in plants.Development1384117–4129. 10.1242/dev.063511
26
ItoM.IwaseM.KodamaH.LavisseP.KomamineA.NishihamaR.et al (1998). A novel cis-acting element in promoters of plant B-type cyclin genes activates M phase–specific transcription.Plant Cell10331–341. 10.1105/tpc.10.3.331
27
JhaA.DubeyR. (2004). Carbohydrate metabolism in growing rice seedlings under arsenic toxicity.J. Plant Physiol.161867–872. 10.1016/j.jplph.2004.01.004
28
JhaA.DubeyR. (2005). Effect of arsenic on behaviour of enzymes of sugar metabolism in germinating rice seeds.Acta Physiol. Plant.27341–347. 10.1007/s11738-005-0010-x
29
JonesJ. D.DanglJ. L. (2006). The plant immune system.Nature444:323.
30
Jones-RhoadesM. W.BartelD. P. (2004). Computational identification of plant microRNAs and their targets, including a stress-induced miRNA.Mol. Cell14787–799. 10.1016/j.molcel.2004.05.027
31
KatsuharaM.SasanoS.HorieT.MatsumotoT.RheeJ.ShibasakaM. (2014). Functional and molecular characteristics of rice and barley NIP aquaporins transporting water, hydrogen peroxide and arsenite.Plant Biotechnol.31213–219. 10.5511/plantbiotechnology.14.0421a
32
KhanE.GuptaM. (2018). Arsenic–silicon priming of rice (Oryza sativa L.) seeds influence mineral nutrient uptake and biochemical responses through modulation of Lsi-1, Lsi-2, Lsi-6 and nutrient transporter genes.Sci. Rep.8:10301.
33
KhraiweshB.ZhuJ.-K.ZhuJ. (2012). Role of miRNAs and siRNAs in biotic and abiotic stress responses of plants.Biochim. Biophys. Acta Gene Regul. Mech.1819137–148. 10.1016/j.bbagrm.2011.05.001
34
KimW.BenhamedM.ServetC.LatrasseD.ZhangW.DelarueM.et al (2009). Histone acetyltransferase GCN5 interferes with the miRNA pathway in Arabidopsis.Cell Res.19899–909. 10.1038/cr.2009.59
35
KornetN.ScheresB. (2009). Members of the GCN5 histone acetyltransferase complex regulate PLETHORA-mediated root stem cell niche maintenance and transit amplifying cell proliferation in Arabidopsis.Plant Cell211070–1079. 10.1105/tpc.108.065300
36
LauterN.KampaniA.CarlsonS.GoebelM.MooseS. P. (2005). microRNA172 down-regulates glossy15 to promote vegetative phase change in maize.Proc. Natl. Acad. Sci. U.S.A.1029412–9417. 10.1073/pnas.0503927102
37
Le HéchoI.MateraV. (2001). “Arsenic behavior in contaminated soils: mobility and speciation,” in Heavy Metals Release in Soils, edsSelimH. M.SparksD. L. (CRC Press), 213–241.
38
LescotM.DéhaisP.ThijsG.MarchalK.MoreauY.Van de PeerY.et al (2002). PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences.Nucleic Acids Res.30325–327. 10.1093/nar/30.1.325
39
LiG.SantoniV.MaurelC. (2014). Plant aquaporins: roles in plant physiology.Biochim. Biophys. Acta Gen. Subj.18401574–1582. 10.1016/j.bbagen.2013.11.004
40
LiN.WangJ.SongW.-Y. (2016). Arsenic uptake and translocation in plants.Plant Cell Physiol.574–13. 10.1093/pcp/pcv143
41
LiuX.LuoM.ZhangW.ZhaoJ.ZhangJ.WuK.et al (2012). Histone acetyltransferases in rice (Oryza sativa L.): phylogenetic analysis, subcellular localization and expression.BMC Plant Biol.12:145. 10.1186/1471-2229-12-145
42
MaJ. F.YamajiN.MitaniN.XuX.-Y.SuY.-H.McGrathS. P.et al (2008). Transporters of arsenite in rice and their role in arsenic accumulation in rice grain.Proc. Natl. Acad. Sci. U.S.A.1059931–9935. 10.1073/pnas.0802361105
43
MartinR.DowlingK.PearceD.SillitoeJ.FlorentineS. (2014). Health effects associated with inhalation of airborne arsenic arising from mining operations.Geosciences4128–175. 10.3390/geosciences4030128
44
MattsJ.JagadeeswaranG.RoeB. A.SunkarR. (2010). Identification of microRNAs and their targets in switchgrass, a model biofuel plant species.J. Plant Physiol.167896–904. 10.1016/j.jplph.2010.02.001
45
MehargA. A.Hartley-WhitakerJ. (2002). Arsenic uptake and metabolism in arsenic resistant and nonresistant plant species.New Phytol.15429–43. 10.1046/j.1469-8137.2002.00363.x
46
MishraS.JhaA.DubeyR. (2011). Arsenite treatment induces oxidative stress, upregulates antioxidant system, and causes phytochelatin synthesis in rice seedlings.Protoplasma248565–577. 10.1007/s00709-010-0210-0
47
MlotshwaS.YangZ.KimY.ChenX. (2006). Floral patterning defects induced by Arabidopsis APETALA2 and microRNA172 expression in Nicotiana benthamiana.Plant Mol. Biol.61781–793. 10.1007/s11103-006-0049-0
48
MoldovanD.SpriggsA.YangJ.PogsonB. J.DennisE. S.WilsonI. W. (2009). Hypoxia-responsive microRNAs and trans-acting small interfering RNAs in Arabidopsis.J. Exp. Bot.61165–177. 10.1093/jxb/erp296
49
MoonS.KimS.-R.ZhaoG.YiJ.YooY.JinP.et al (2013). Rice glycosyltransferase1 encodes a glycosyltransferase essential for pollen wall formation.Plant Physiol.161663–675. 10.1104/pp.112.210948
50
MukhopadhyayR.BhattacharjeeH.RosenB. P. (2014). Aquaglyceroporins: generalized metalloid channels.Biochimi. Biophys. Acta Gen. Subj.18401583–1591. 10.1016/j.bbagen.2013.11.021
51
PandeyA. K.GautamA.DubeyR. S. (2019). Transport and detoxification of metalloids in plants in relation to plant-metalloid tolerance.Plant Gene17:100171. 10.1016/j.plgene.2019.100171
52
PengT.SunH.DuY.ZhangJ.LiJ.LiuY.et al (2013). Characterization and expression patterns of microRNAs involved in rice grain filling.PLoS One8:e54148. 10.1371/journal.pone.0054148
53
PetrovV.HilleJ.Mueller-RoeberB.GechevT. S. (2015). ROS-mediated abiotic stress-induced programmed cell death in plants.Front. Plant Sci.6:69. 10.3389/fpls.2015.00069
54
PraneenararatT.TakagiT.IwasakiW. (2012). Integration of interactive, multi-scale network navigation approach with Cytoscape for functional genomics in the big data era.BMC Genomics13(Suppl. 7):S24. 10.1186/1471-2164-13-S7-S24
55
RichmondG. S.SmithT. K. (2011). Phospholipases A1.Int. J. Mol. Sci.12588–612. 10.3390/ijms12010588
56
SchwabR.PalatnikJ. F.RiesterM.SchommerC.SchmidM.WeigelD. (2005). Specific effects of microRNAs on the plant transcriptome.Dev. Cell8517–527. 10.1016/j.devcel.2005.01.018
57
SchwarzS.GrandeA. V.BujdosoN.SaedlerH.HuijserP. (2008). The microRNA regulated SBP-box genes SPL9 and SPL15 control shoot maturation in Arabidopsis.Plant Mol. Biol.67183–195. 10.1007/s11103-008-9310-z
58
ServetC.e SilvaN. C.ZhouD.-X. (2010). Histone acetyltransferase AtGCN5/HAG1 is a versatile regulator of developmental and inducible gene expression in Arabidopsis.Mol. Plant3670–677. 10.1093/mp/ssq018
59
ShankarS.ShankerU. (2014). Arsenic contamination of groundwater: a review of sources, prevalence, health risks, and strategies for mitigation.Sci. World J.2014:304524.
60
ShikataM.KoyamaT.MitsudaN.Ohme-TakagiM. (2009). Arabidopsis SBP-box genes SPL10, SPL11 and SPL2 control morphological change in association with shoot maturation in the reproductive phase.Plant Cell Physiol.502133–2145. 10.1093/pcp/pcp148
61
SrivastavaS.SrivastavaA. K.SuprasannaP.D’souzaS. (2012). Identification and profiling of arsenic stress-induced microRNAs in Brassica juncea.J. Exp. Bot.64303–315. 10.1093/jxb/ers333
62
StoevaN.BerovaM.ZlatevZ. (2003). Physiological response of maize to arsenic contamination.Biol. Plant.47449–452. 10.1023/b:biop.0000023893.12939.48
63
SunX.FanG.SuL.WangW.LiangZ.LiS.et al (2015). Identification of cold-inducible microRNAs in grapevine.Front. Plant Sci.6:595. 10.3389/fpls.2015.00595
64
SunkarR.LiY.-F.JagadeeswaranG. (2012). Functions of microRNAs in plant stress responses.Trends Plant Sci.17196–203. 10.1016/j.tplants.2012.01.010
65
SunkarR.ZhuJ.-K. (2004). Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis.Plant Cell162001–2019. 10.1105/tpc.104.022830
66
TangZ.XuM.ItoH.CaiJ.MaX.QinJ.et al (2019). Deciphering the non-coding RNA-level response to arsenic stress in rice (Oryza sativa).Plant Signal. Behav.14:1629268. 10.1080/15592324.2019.1629268
67
TiwariN.KumarV.GeddaM. R.SinghA. K.SinghV. K.GannavaramS.et al (2017). Identification and characterization of miRNAs in response to leishmania donovani infection: delineation of their roles in macrophage dysfunction.Front. Microbiol.8:314. 10.3389/fmicb.2017.00314
68
VerbruggenN.HermansC.SchatH. (2009). Mechanisms to cope with arsenic or cadmium excess in plants.Curr. Opin. Plant Biol.12364–372. 10.1016/j.pbi.2009.05.001
69
VernouxT.BenfeyP. N. (2005). Signals that regulate stem cell activity during plant development.Curr. Opin. Genet. Dev.15388–394. 10.1016/j.gde.2005.06.008
70
WangY.HuZ.YangY.ChenX.ChenG. (2009). Function annotation of an SBP-box gene in Arabidopsis based on analysis of co-expression networks and promoters.Int. J. Mol. Sci.10116–132. 10.3390/ijms10010116
71
WilliamsP.PriceA.RaabA.HossainS.FeldmannJ.MehargA. A. (2005). Variation in arsenic speciation and concentration in paddy rice related to dietary exposure.Environ. Sci. Technol.395531–5540. 10.1021/es0502324
72
WuG.PoethigR. S. (2006). Temporal regulation of shoot development in Arabidopsis thaliana by miR156 and its target SPL3.Development1333539–3547. 10.1242/dev.02521
73
WuL.ZhangQ.ZhouH.NiF.WuX.QiY. (2009). Rice microRNA effector complexes and targets.Plant Cell213421–3435. 10.1105/tpc.109.070938
74
XieK.WuC.XiongL. (2006). Genomic organization, differential expression, and interaction of SQUAMOSA promoter-binding-like transcription factors and microRNA156 in rice.Plant Physiol.142280–293. 10.1104/pp.106.084475
75
XingS.SalinasM.HöhmannS.BerndtgenR.HuijserP. (2010). miR156-targeted and nontargeted SBP-box transcription factors act in concert to secure male fertility in Arabidopsis.Plant Cell223935–3950. 10.1105/tpc.110.079343
76
Yamaguchi-ShinozakiK.ShinozakiK. (2006). Transcriptional regulatory networks in cellular responses and tolerance to dehydration and cold stresses.Annu. Rev. Plant Biol.57781–803. 10.1146/annurev.arplant.57.032905.105444
77
YeY.DingY.JiangQ.WangF.SunJ.ZhuC. (2017). The role of receptor-like protein kinases (RLKs) in abiotic stress response in plants.Plant Cell Rep.36235–242. 10.1007/s00299-016-2084-x
78
YiS.-Y.LiuJ. (2009). Combinatorial interactions of two cis-acting elements, AT-rich regions and HSEs, in the expression of tomato Lehsp23. 8 upon heat and non-heat stresses.J. Plant Biol.52:560. 10.1007/s12374-009-9072-4
79
YoshidaS.FornoD. A.CockJ. H. (1971). Laboratory Manual for Physiological Studies of Rice.Los Baños, PH: International Rice Research Institute.
80
YuL. J.LuoY. F.LiaoB.XieL. J.ChenL.XiaoS.et al (2012). Comparative transcriptome analysis of transporters, phytohormone and lipid metabolism pathways in response to arsenic stress in rice (Oryza sativa).New Phytol.19597–112. 10.1111/j.1469-8137.2012.04154.x
81
ZhangH.SelimH. (2006). Modeling the transport and retention of arsenic (V) in soils.Soil Sci. Soc. Am. J.701677–1687. 10.2136/sssaj2006.0035
82
ZhangS.-D.LingL.-Z.ZhangQ.-F.XuJ.-D.ChengL. (2015). Evolutionary comparison of two combinatorial regulators of SBP-Box genes, MiR156 and MiR529, in plants.PLoS One10:e0124621. 10.1371/journal.pone.0124621
83
ZhengM. Z.CaiC.HuY.SunG. X.WilliamsP. N.CuiH. J.et al (2011). Spatial distribution of arsenic and temporal variation of its concentration in rice.New Phytol.189200–209. 10.1111/j.1469-8137.2010.03456.x
84
ZhouX.SunX.CooperK. L.WangF.LiuK. J.HudsonL. G. (2011). Arsenite interacts selectively with zinc finger proteins containing C3H1 or C4 motifs.J. Biol. Chem.28622855–22863. 10.1074/jbc.m111.232926
85
ZhuY.-G.WilliamsP. N.MehargA. A. (2008). Exposure to inorganic arsenic from rice: a global health issue?Environ. Pollut.154169–171. 10.1016/j.envpol.2008.03.015
Summary
Keywords
rice, arsenic, microRNA, Osa-miR156j, gene expression, developmental stage
Citation
Pandey AK, Gedda MR and Verma AK (2020) Effect of Arsenic Stress on Expression Pattern of a Rice Specific miR156j at Various Developmental Stages and Their Allied Co-expression Target Networks. Front. Plant Sci. 11:752. doi: 10.3389/fpls.2020.00752
Received
20 November 2019
Accepted
12 May 2020
Published
16 June 2020
Volume
11 - 2020
Edited by
Eric Ruelland, UMR 7618 Institut d’Écologie et des Sciences de l’Environnement de Paris (IEES), France
Reviewed by
Jirina Szakova, Czech University of Life Sciences Prague, Czechia; Junhong Zhang, Zhejiang Agriculture and Forestry University, China
Updates
Copyright
© 2020 Pandey, Gedda and Verma.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ashok K. Verma, vrmashok75@gmail.com
†These authors have contributed equally to this work
This article was submitted to Plant Abiotic Stress, a section of the journal Frontiers in Plant Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.