Next Article in Journal
Effect of Prohydrojasmon on the Growth of Eggplant and Komatsuna
Next Article in Special Issue
Metabolic Responses to Waterlogging Differ between Roots and Shoots and Reflect Phloem Transport Alteration in Medicago truncatula
Previous Article in Journal
Correction of Potassium Fertigation Rate of Apple Tree (Malus domestica Borkh.) in Central Russia during the Growing Season
Previous Article in Special Issue
The Molecular Regulatory Pathways and Metabolic Adaptation in the Seed Germination and Early Seedling Growth of Rice in Response to Low O2 Stress
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Differential Expression of Maize and Teosinte microRNAs under Submergence, Drought, and Alternated Stress

1
Laboratorio de Biotecnología Vegetal, Instituto de Biotecnología, Universidad del Papaloapan, Tuxtepec 68301, Mexico
2
Division de Estudios de Posgrado, Universidad del Papaloapan, Tuxtepec 68301, Mexico
3
Cátedras CONACyT-UNPA, Universidad del Papaloapan, Tuxtepec 68301, Mexico
4
Center for Agroforestry Mega Data Science, Haixia Institute of Science and Technology, Fujian Agriculture and Forestry University, Fuzhou 350002, China
*
Author to whom correspondence should be addressed.
Plants 2020, 9(10), 1367; https://doi.org/10.3390/plants9101367
Submission received: 15 August 2020 / Revised: 1 October 2020 / Accepted: 11 October 2020 / Published: 15 October 2020
(This article belongs to the Special Issue Plant Responses to Hypoxia)

Abstract

:
Submergence and drought stresses are the main constraints to crop production worldwide. MicroRNAs (miRNAs) are known to play a major role in plant response to various stresses. In this study, we analyzed the expression of maize and teosinte miRNAs by high-throughput sequencing of small RNA libraries in maize and its ancestor teosinte (Zea mays ssp. parviglumis), under submergence, drought, and alternated stress. We found that the expression patterns of 67 miRNA sequences representing 23 miRNA families in maize and other plants were regulated by submergence or drought. miR159a, miR166b, miR167c, and miR169c were downregulated by submergence in both plants but more severely in maize. miR156k and miR164e were upregulated by drought in teosinte but downregulated in maize. Small RNA profiling of teosinte subject to alternate treatments with drought and submergence revealed that submergence as the first stress attenuated the response to drought, while drought being the first stress did not alter the response to submergence. The miRNAs identified herein, and their potential targets, indicate that control of development, growth, and response to oxidative stress could be crucial for adaptation and that there exists evolutionary divergence between these two subspecies in miRNA response to abiotic stresses.

1. Introduction

As a consequence of global warming, hydrological fluctuation events such as excessive rainfall and droughts are common and projected to continue in the future, affecting economic activities [1]. Alterations in water availability in the field, either caused by water surplus or deficit, produce water stress in plants that negatively impacts the growth and productivity of crops worldwide [2].
Flooding affects the properties of soil and the composition of associated microbial communities [3] and reduces the availability of nutrients and flux of oxygen to the plant. This leads to hypoxic stress at the cellular level, especially when the column of water exceeds the length of the stem (i.e., submergence) [4,5]. In response to flooding stress, plants express genes known as hypoxia core genes (HCG) that promote anaerobic metabolism [6]. Additionally, some plants use the escape strategy by redirecting growth to stem elongation to overpass the column of water and maintain the oxygen flux. Plants may also delay growth, flowering, and other processes that are energetically expensive, as part of a strategy known as quiescence to save resources useful to reestablish normal growth when the water level recedes [7]. On the contrary, when water availability is limited, as in the case of drought, plants undergo a similar series of changes to adapt to the adverse conditions, which include stomatal closure with limitations in CO2 uptake and reduction in photosynthetic activity, increase in root elongation rate, early or delayed flowering time, and other physiological and morphological changes [8].
Cellular dehydration caused by drought, and hypoxia driven by submergence, trigger a cascade of adaptive responses that are regulated at the molecular level in plants and are directed to maintain vital functions and protect structures from damage. Since cellular protection is a priority during stress, it has been shown that some responses to dehydration and excess water have common molecular effectors acting as nodes in the crosstalk between responses to drought and submergence [9,10,11]. This may have allowed wild plants to evolve tolerance mechanisms to cope with alternate drought and flooding events, which commonly occur consecutively in some regions of the world [12,13].
Maize (Zea mays ssp. mays) is one of the main dietary cereals in the world [14]. In the United States, the world’s main maize producer, extreme drought and excessive rainfall are the first and second major causes of losses in maize production, respectively, with both affecting yields to a comparable extent [15]. According to archeological, botanical, and genetic evidence, maize was domesticated from its single ancestor teosinte (Zea mays ssp. parviglumis) around 10,000 years ago [16,17]. It is estimated that maize retained 83% of the nucleotide diversity from its ancestor [18]. During domestication and breeding, maize has been subject to selection for traits that are beneficial to improve feeding quality and yield, while other traits like tolerance to stress could have been lost. For this reason, teosinte is considered a living reservoir of genes and mechanisms that could be of importance to improve tolerance to stress in maize [19,20].
MicroRNAs (miRNAs) have been proposed as promising targets for developing plants with improved tolerance to multiple abiotic stresses [21]. miRNAs are small RNA molecules with pivotal roles in the response to environmental stresses in plants that usually act as negative regulators of gene expression via silencing target genes, which are recognized by sequence complementarity and subsequently targeted for degradation and/or translation inhibition [22]. With the development of new platforms for the high-throughput sequencing of small RNAs, several miRNAs involved in the response to drought or submergence have been described in plants such as Arabidopsis thaliana [23,24], Brachypodium distachyon [25,26], Oryza sativa [27], Nelumbo nucifera, [28] and Zea mays [29,30].
According to miRBase, 325 mature miRNAs from 174 miRNA precursors have been identified in maize [31], and some of them have been reported to be responsive to drought [32,33], submergence, [34] or waterlogging [35]. For example, miR159, miR164, miR167, miR393, miR408, and miR528 are upregulated by short-term waterlogging in roots of maize lines with high tolerance and suggest the involvement of hormonal control in the response mediated by miRNAs [30]. Another study of two maize lines, with contrasting tolerance to drought, observed the downregulation of miR164 and upregulation of miR159, miR390, and miR398 in the tolerant line compared with the sensitive line [29]. However, it is not clear how miRNAs are regulated in teosinte under submergence or drought stress and whether miRNA response to these stresses is conserved during maize domestication and breeding.
Therefore, this study aimed to identify microRNAs expressed in maize and teosinte under submergence and drought stress in order to analyze the conserved and differential mechanisms of response in these plants. For this, we used high-throughput small RNA sequencing to profile maize and teosinte under submergence, drought, and consecutive drought and submergence stresses. We found a group of miRNAs that were regulated by drought or submergence, including some miRNAs that were differentially regulated in the two species. Analyses of small RNA data from teosinte subject to alternate treatments of drought and submergence indicated that miRNA response to drought was attenuated when submergence was the first stress treatment, but miRNA response to submergence did not change when drought stress was applied first. The potential targets of these miRNAs indicate that control of development, growth, and response to oxidative stress could be crucial for understanding the conservation and divergence of stress tolerance between maize and teosinte.

2. Results

2.1. Physiological Response of Maize and Teosinte to Submergence and Drought

With the objective to determine a time period of treatment for RNA collection, maize seedlings with two vegetative leaves and 14 days after sowing (stage V2, 14 DAS) were exposed to submergence for 2, 4, and 6 days (Figure 1A). All treatments induced growth reduction in maize, but plants submerged for 4 or 6 days developed necrotic spots on the leaves, and plants treated for 6 days were not able to recover from the stress. Therefore, a 2 day submergence treatment was selected for subsequent assays. Similarly, drought assays were performed on maize seedlings that were deprived of water until they reached half their initial weight and were maintained under this level of water limitation for 2, 4, 6, and 8 days (Figure 1B).
Drought caused a reduction in the growth of maize plants when compared with control plants, but this reduction was stabilized, which may be explained by the fact that pots were maintained at a constant, intermediate water deficit condition, and not exposed to progressive drought stress. Based on these observations, maize and teosinte plants were exposed to submergence for 2 days and drought for 8 days, and the submergence tolerance coefficient (STC) and relative water content (RWC) were registered for the respective treatments. Submergence affected the growth of maize and teosinte seedlings to a similar extent, since STC was not significantly different (Figure 1C). There were also no significant differences in RWC between maize and teosinte in control, drought, and 24 h of recovery after the stress (Figure 1D). These results suggest that both plants are able to sense and respond to stress by adjusting water status and growth to a similar extent at this stage. This is expected because both maize CML496 and teosinte grow in conditions prevalent in the tropics, where drought, high temperature, and extreme seasonal variations in monthly rainfall are commonly observed. Therefore, this constitutes a solid system to investigate the intra-specific conservation and divergence in the response to stress between maize and its direct ancestor that lead to adaptation in both plants.

2.2. Overall Analysis of Small RNA Sequencing Data

The alterations in growth, STC, and RWC observed in maize and teosinte seedlings in the assays suggest that the times and stress intensities can help to study molecular responses in both plants. Accordingly, maize and teosinte seedlings exposed to 2 days of submergence or 6 days of constant drought were selected for profiling the microRNA populations involved in the corresponding responses. In addition, teosinte seedlings exposed to alternated stress of drought-submergence and submergence-drought were also selected.
High-throughput sequencing of small RNA libraries prepared from teosinte and maize under submergence, drought, and alternated stress rendered from 16 to 46 million raw reads (Table 1), with PHRED score above 30 (Figure S1). After removing low-quality reads and adapters, the libraries were analyzed for size distribution of reads, revealing an enrichment of 24 nt sRNAs across the libraries (Figure S2). Reads less than 18 or longer than 35 nt were removed, resulting in a number of clean reads ranging from 14.8 to 41.9 million (Table 1).
On average, 99.7% of the clean reads were mapped to the maize genome. From the reads mapped to the maize genome, less than 3% represented plant miRNAs reported in the miRBase (Table 1). Interestingly, libraries from teosinte showed the lowest percentage of miRNAs (approximately 0.74% for teosinte and approximately 2.19% for maize), suggesting a potential underrepresentation effect caused by the absence of teosinte miRNAs in the miRBase. However, the percentage of teosinte reads that were mapped to the maize genome (approximately 99.75%) was the same as that for maize (approximately 99.77%), indicating that differences are likely due to lower levels of miRNA expression in teosinte. In addition, a reduced abundance of plant miRNAs under submergence stress was observed in both maize and teosinte, and in teosinte treated with drought-submergence stress compared with the corresponding controls (Table 1). After removing low frequency reads and normalization to tags per million (TPM), 178 miRNA sequences in total (118 mapping to maize and 60 to miRNAs of other plants) were expressed across the libraries (Table S1).

2.3. Global Changes in Expression of miRNA Sequences in Response to Stress

Further analysis revealed that 164 miRNA sequences (110 from maize and 54 annotated in other plants) were differentially expressed in maize and/or teosinte in response to at least one of the treatments (Table S1).
Regarding submergence stress, maize and teosinte showed common expression profiles of 10 upregulated and 72 downregulated miRNA sequences (Figure 2A). A total of 42 miRNA sequences (14 upregulated and 28 downregulated) were regulated exclusively in maize, while 27 miRNA sequences (4 upregulated and 23 downregulated) were exclusive of teosinte (Figure 2A). A similar trend in the number of miRNAs responsive to stress, with lesser miRNAs being regulated in teosinte than in maize, was observed for drought treatments (Figure 2B). The analysis revealed 70 miRNA sequences (21 upregulated, 49 downregulated) regulated exclusively in maize, while 27 were differentially expressed only in teosinte (17 upregulated, 10 downregulated).
The overlap between submergence and drought responses in maize and teosinte showed that in maize, most of the miRNAs (68 sequences) were downregulated in both submergence and drought (Figure 2C), whereas in teosinte, the number of sequences downregulated exclusively in submergence was the highest (69 sequences; Figure 2D).
In addition, small RNA libraries from teosinte exposed to alternated water stress were analyzed to assess the effects of consecutive treatments of drought and submergence on miRNA expression. Figure 2E shows that 18 miRNA sequences (5 upregulated and 13 downregulated) were differentially expressed in a similar pattern in response to submergence, drought, and alternated drought-submergence treatments. The number of miRNA sequences (69) overlapping between drought-submergence and submergence (representing 63% of total sequences in submergence) was higher compared with only 22 miRNA sequences between drought-submergence and drought (39% of total sequences in drought; Figure 2E), indicating that when submergence is the second stress, it has the strongest effect on miRNA expression.
When submergence is followed by drought, 21 out of 56 miRNA sequences regulated under drought overlapped between submergence-drought and drought treatments (37% of sequences in drought), indicating that when submergence is the first stress, it affects the subsequent response to drought in teosinte (Figure 2F).

2.4. Redundancy of miRNA Reads Annotated to miRBASE and Grouping of Representative Sequences

The data showed a high level of redundancy of reads in the libraries annotated to miRBASE, mainly caused by positional variants of mature miRNAs and other sequences matching miRNA precursors. Regrouping of mature miRNAs, positional variants derived from the same mature sequence, and sequences that exactly match those from other plants but with a close similarity to known maize miRNAs (i.e., differing by one or two nucleotides in the 5′ or 3′ ends), rendered 67 representative sequences belonging to 23 miRNA families (Table S2).

2.5. Search for Novel miRNAs

In order to identify new variants of known miRNAs in teosinte, a search was performed against the miRBASE, allowing up to two mismatches in the mature or precursor sequences. This rendered 102 polymorphic variants that showed the same expression profiles as those of their corresponding exact sequences (Table S3). However, these sequences failed to provide a perfect hit to the reference genome, which did not allow the analysis of a putative precursor sequence. We also searched for potential novel miRNAs in our datasets and generated several candidates, but after a detailed inspection, they did not comply with the criteria for annotation of novel miRNAs in terms of secondary structure and/or distribution of reads along the potential precursor [36].

2.6. Analysis of miRNAs Responsive to Submergence and Drought in Maize and Teosinte

In the group of 67 representative non-redundant miRNA sequences, 13 were identified as responsive to at least one of the stress conditions, with fold change values ranging from −4.39 for the most downregulated miRNA (miR166bd) to 2.17 for the most upregulated miRNA (miR319b) (Table 2). Besides these miRNAs, others showing differential expression such as miR1511, miR2916, miR482, and miR4995 (Table S2), were discarded from the analysis due to failure to hit the maize genome or originating from ribosomal RNA or transposons.
Overall, the grouping of miRNAs presented a clear picture of the stress response. Changes in miRNA expression were more marked in maize than in teosinte for most miRNAs downregulated under submergence, suggesting that maize plants could be more sensitive or reactive to submergence than teosinte in terms of miRNA responses. miRNAs that were specifically downregulated in maize were miR159ab, miR164e, miR166bd, miR167cdeg, miR169cr, miR396cd, and miR398b (Table 2).
In maize, drought stress repressed the expression of five miRNAs (miR156k, miR164e, miR166bd, miR167cdeg, and miR528a). In teosinte, only miR166bd and miR408 were downregulated, while five miRNAs were upregulated under drought (miR156k, miR159ab, miR164e, miR319b, and miR396cd) (Table 2).

2.7. Analysis of miRNAs Responsive to Alternated Stress in Teosinte

When plants were exposed to drought-submergence conditions, the miRNA expression profile resembled that of submergence, suggesting that at the end of the drought treatment, plants were not compromised to respond to a submergence event (Table 3). The changes between drought-submergence and submergence stress were limited only to the expression of miR408. This miRNA target gene encodes plastocyanin (PLC) and laccase (LAC) involved in copper homeostasis and the formation of lignin.
When plants were exposed to submergence-drought conditions, most miRNAs were regulated in a way similar to that in plants under drought alone, but the response was attenuated. For example, miR159ab, miR319b, and miR396cd were upregulated in drought as a single stress, but when drought was preceded by submergence, regulation of these miRNAs was maintained within the same trend but at low levels (Table 3). These results suggest that plants are plastic in response to environmental cues and have the capacity to oscillate to adapt to different stresses.

2.8. Assessment and Validation of miRNA Expression by Quantitative RT-PCR

Quantitative RT-PCR assays were performed using a modified SL-RTPCR method for seven miRNAs regulated by stress. Based on the fold change values of miRNAs expressed in the libraries, miR166c-5P was selected as a constitutive control for qRT-PCR assays, since it showed intermediate abundance levels and stable regulation in single and alternated stress conditions (Table 2 and Table 3). SL-RTPCR assays showed changes in expression similar to those observed by analysis of sequencing data for miR156k, miR159ab, miR167cdeg, miR396cd, miR398ab, miR408b, and miR528ab (Figure 3A–G).
For example, based on sequencing data, miR167 had a fold-change value of 0.12 and 0.17 in maize under submergence and drought, respectively. The values obtained by SL-RTPCR were 0.39 (±0.03) and 0.40 (±0.06), respectively, which confirms that miR167 is downregulated by these treatments and thus the expression trends for upregulated or downregulated miRNAs were maintained in both methods.

2.9. Other Small RNAs

Mitochondria and chloroplasts are the main compartments for energy production that are affected by low water and oxygen availability and may be a source of signals to integrate the cellular response [37,38]. A search was performed in the libraries with the objective of identifying small RNAs potentially originating from maize mitochondria or plastids. For both mitochondria and plastids, the highest number of reads mapped to the subunit 2 of NADH dehydrogenase (Table S4). The reads were derived from a region spanning 53 bp of NADH dehydrogenase subunit 2 with an abundance of up to 820 TPM, comparable to miR408 abundance (Table S2). The second highest number of reads mapped to the trnL-trnF intergenic spacer and tRNA-Phe (trnF) gene from the maize chloroplast. Interestingly, these reads derived from a 32 bp region of the gene did not show a clear expression trend in response to stress; however, a high accumulation of reads in teosinte compared with that in maize was evident, with up to 2273 TPM, comparable to miR398 abundance (Tables S2 and S5). The causes and implications of the existence of these small RNA signatures need to be investigated in the future.

3. Discussion

Drought and submergence represent complex conditions that can be dissected into several components such as hypoxia, nutrient, light, osmotic, temperature, and oxidative stress [39]. Understanding how the different components act to integrate the response to stress at all levels of regulation is fundamental to underpin the improvement of tolerance to drought and flooding stress in crops.
miRNAs are a class of small RNAs that act as regulators of gene expression at the transcriptional level. In this study, high-throughput sequencing and analysis of small RNA populations from maize and teosinte under submergence and drought revealed that submergence caused a global reduction in the representation of miRNA expression compared with control or drought conditions. This was observed in roots of wild tomato treated with hypoxia, where only 1.45% of total reads were expressed in hypoxia-treated roots against 2.45% expressed in control roots [40].
Further exploration of the expression of specific miRNA family members allowed the identification of 13 miRNAs with differential expression profiles in response to submergence and drought. The analysis of their expression patterns and functions of their predicted, or previously confirmed targets, allowed the recognition of two main components of the responses to submergence and drought in maize: transcriptional regulation and antioxidant activity.

3.1. Transcriptional Regulation

The growth and development processes are driven by hormonal signaling pathways that are fundamental in the response to submergence and drought in plants [41], where transcription factors and miRNAs participate. Several differentially expressed miRNAs identified in this work are known to act over target genes encoding transcription factors such as miR156 (squamosa promoter-binding protein like, SPL) [42,43,44], miR159 (giberellic acid-MYB transcription factor, GAMYB [28,29]), miR164 (NAM-ATAF-CUC domain transcription factor NAC, MYB [45]), miR166 (basic leucine zipper transcription factor, HD-ZIP III [46,47]), miR167 (auxin-responsive factor, ARF [30,35,48]), miR169 (nuclear transcription factor Y subunit alpha, NFYA [49,50,51,52]), miR319 (teosinte branched1/cycloidea/proliferating cell nuclear antigen factor, TCP [53]), and miR396 (growth-regulating factor, GRF [54,55]). Other target genes predicted in this study are listed in Table S6.
Most miRNAs in this category have been related to responses to abiotic stress. miR156 controls developmental transitions and flowering through the activity of their targets (SPL) [42,43,44] which is upregulated in Arabidopsis roots under hypoxia and in lotus under submergence [23,28]. Other miRNAs such as miR166 were upregulated in maize roots in response to waterlogging [34] and downregulated in Arabidopsis roots in response to hypoxia [23]. miR164 participates in the development of lateral roots [45] and the response to waterlogging in maize roots [30]. miR159 controls petiole elongation and flowering [56] and is upregulated in maize roots exposed to waterlogging [29]. miR167 is upregulated by waterlogging in maize roots [30,34], and downregulated in submerged lotus seedlings [28].
In this study, miRNAs in the category of transcriptional regulation overlapped between maize and teosinte in response to submergence and between submergence and drought in maize (Figure 4), indicating that the control of growth is of central importance for adaptation in maize and teosinte. Interestingly, all these miRNAs were downregulated in the overlaps. Downregulation of this set of miRNAs means that their targets may be released from post transcriptional control, which may activate hormone-responsive genes to initiate elongation (miR159), cell elongation and differentiation (miR167), transition to flowering (miR156), and other responses. In addition, it is possible that changes observed in miRNA expression at the transcriptional level are not only directed to increase growth, but also to mediate feedback in hormonal pathways, or to prepare for reoxygenation or rehydration after stress. Interestingly, the response to drought in teosinte was regulated by five miRNAs that were upregulated (miR156, miR164, miR159, miR319, and miR396). The upregulation of these miRNAs suggests that growth processes may be restricted in teosinte during drought stress (Figure 4).

3.2. Antioxidant Activity

Submergence and drought stress result in excessive production of reactive oxygen species (ROS) and efficient mechanisms to cope with oxidative injury are determinants of tolerance [57]. miR398, miR408, and miR528 are included in this category. miR398 target genes encode Cu/Zn superoxide dismutases (CSD1 and CSD2) that act in the defense against toxic ROS in Arabidopsis [58]. In Phaseolus vulgaris, miR398 is downregulated under drought and submergence [59] and is also downregulated in response to submergence in Arabidopsis [60]. We predicted other miR398 targets including a selenium-binding protein (SBP) (Table S6). An SBP was recently reported to be induced in Arabidopsis in response to submergence [60] and confirmed as a miR398 target by degradome analysis in maize [61]. In this study, miR398 was downregulated in maize and teosinte under submergence, thereby implicating a potential increase in antioxidant activity.
Other members, such as miR408 and miR528, target genes encoding the cupredoxins PLC and LAC [62,63,64], and other genes controlling circadian clock and flowering time [65,66]. PLC functions as an electron transporter and LAC participates in the formation of lignin by oxidation. miR408 and miR528 have been previously reported to respond to submergence in Lotus seedlings and maize roots [28,30]. In this study, miR408 and miR528 were downregulated in maize and teosinte under submergence, suggesting the maintenance of electron flux, oxidation homeostasis, and lignin synthesis (Figure 4).

3.3. Comparison of miRNA Expression in Maize and Teosinte

Previous work demonstrated a high level of conservation of miRNA mature sequences between maize and teosinte [67]. Thus, divergence in tolerance to stress is most likely the result of differences in gene expression between these plants. Overall, in this work maize and teosinte showed similar responses to submergence, with a clear trend toward downregulation of miRNAs involved in the control of growth, development, antioxidant response, and copper homeostasis. However, the response to submergence was more intense in maize than in teosinte for most of the miRNAs responsive to stress, especially for those miRNAs controlling growth. This suggests that maize is more reactive to submergence and responds by releasing the control on their target genes, which in turn may accelerate growth to escape from stress. In teosinte, these responses were less dramatic, suggesting that a finer tuning of expression could be a key for tolerance in this plant. A similar pattern was observed in a study comparing miRNA expression in inbred maize lines with different levels of tolerance to waterlogging, where the stress sensitive line reacted by downregulating most of the miRNAs and induced target genes to accelerate growth. Instead, the tolerant line responded with a moderate level of downregulation, or even upregulation of some miRNAs to repress growth under stress [30].
In the case of drought, maize responded in a similar way as that in submergence by downregulating miRNAs to induce growth, while teosinte upregulated miRNAs to restrict growth and probably save energy resources. Regarding defense against oxidative stress, oxidation homeostasis, and reinforcement of cell structure, both plants maintained the downregulation of miRNAs, suggesting the possibility that their target genes were actively protecting the plants from the oxidative damage caused by stress.

3.4. Effects of Alternated Stress on Teosinte

Alternated events of drought and submergence are common in some regions of the world where plants experience the succession of these events throughout their life cycle. In this study, teosinte plants treated with alternated stress revealed that drought as the first stress had no effect on the capacity to respond to a subsequent event of submergence (Figure 5). In fact, plants treated with alternated stress, and plants treated only with submergence had resembling replicates for most of the miRNAs. This suggests that teosinte was able to rapidly reverse the effects of drought and redirect the response towards submergence, even though drought had induced clear symptoms of stress in plants. The only miRNA that did not respond to submergence as the second stress was miR408, but a functional overlap with miR528 may have supplied the activity of LAC and PLC.
In the opposite order, submergence followed by drought produced an attenuated response to drought. The effects of submergence as the first stress had a long-term influence, considering that between the end of submergence and the beginning of drought, there was a gap of 8 days. This effect also suggests the plastic nature of plants to adapt to adverse factors and oscillate between different stresses.
Research has shown that recovery from submergence represents an additional stress by itself. Upon de-submergence, plants suffer reoxygenation stress conducive to the accumulation of ROS, reillumination stress, dehydration, and senescence [68]. During these processes, hormonal signals, ROS, protective proteins, and other molecular effectors could be induced in plants over a period of time that could be useful to survive any other event of stress such as drought, establishing an overlap between submergence and drought tolerance, as described in rice harboring the SUB1A gene [9].
Several aspects in the involvement of miRNAs in the responses to stress still need to be addressed for constructing a solid system that improves stress tolerance in plants. These include a detailed study of expression dynamics between miRNA and target genes during exposure to combined or alternated stress and during recovery in lines with contrasting tolerance, implications of other factors in the response such as intensity and quality of light during stress, circadian cycle on miRNA expression, plant age, the function of other sRNAs involved, and additional layers of regulation of gene expression. The present study constitutes a base for information that could be useful in further efforts to broaden the knowledge on the roles of miRNAs in stress response and explore miRNAs as tools for improving tolerance to multiple stresses.

4. Materials and Methods

4.1. Plant Materials

Maize seeds (Zea mays ssp. mays) of tropical homozygous line CML496, and teosinte (Zea mays ssp. parviglumis) were obtained from Centro de Investigación y de Estudios Avanzados (CINVESTAV) and Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias (INIFAP) seed collections, respectively. The seeds were germinated directly in pots with Sunshine Mix 3 (N: 36 ppm from NO3/9 ppm from NH4; P: 7 ppm; K: 49 ppm; Ca: 45 ppm; Mg: 30 ppm; S: 39 ppm; B: 0.046 ppm; Cu: 0.01 ppm; Fe: 0.478 ppm; Mn: 0.251 ppm; Mo: 0.035 ppm; Na: 10 ppm; Cl: 8 ppm; Al: 0.48 ppm) (Sun Gro Horticulture Distribution Inc., Agawam, MA, USA), three seeds per pot. In the case of teosinte, the seeds were mechanically scarified to remove the coat before sowing.

4.2. Growth Conditions and Treatments

Maize and teosinte plants were grown in a plant growth room at 23 °C and 16 h/8 h light/dark cycle in order to lengthen the vegetative growth phase. The drought stress applied in this study was set to maintain plants at a constant intermediate intensity of drought rather than a progressive drought stress. For drought stress, pots with plants in vegetative stage V2, 12 days after sowing (12 DAS) were water-deprived for 6 days by interrupting water supply until half of their initial weight (water deficit) was reached. Then, water was added in small amounts on a daily basis to constantly maintain the plants in a water deficit condition at half their weight for 6 days. Submergence stress was applied on plants at 26 DAS, with pots in empty 32 cm deep submergence tanks that were filled with filtered tap water to increase the water column by 6 cm every 2 h for 10 h, until the maximum capacity was reached, and plants were fully submerged. The plants were maintained in this condition for 48 h after the first addition of water. For all the treatments, growth conditions were maintained at 23 °C and 16 h/8 h light/dark cycle.
For drought followed by submergence (drought-submergence), teosinte plants (14 DAS) were water-deprived for 6 days followed by 6 days under constant drought stress and then submerged for 48 h as previously described. For submergence followed by drought (submergence-drought), teosinte plants (12 DAS) were submerged for 48 h and removed from the water tank followed by depriving them of water until the plants reached half their initial weight (8 days). Then, small amounts of water were added to maintain a constant water deficit condition for 6 days. The single and alternated stress experiments with maize and teosinte were designed to end simultaneously at 28 DAS. At the end of the experiments, the shoots were immediately collected, frozen in liquid nitrogen, and stored at −80 °C for further processing. All the assays, including control and treatments, were performed with three biological replicates, with each replicate containing six plants.

4.3. Assessment of Physiological and Morphological Effects of Stress in Maize and Teosinte

For drought stress, the relative water content (RWC) was measured in the second leaf of maize plants under drought stress and control conditions, as previously described [69]. Briefly, square leaf sections (2 × 2 cm) were excised from plants and weighed to obtain the fresh weight, then submerged in deionized water for 24 h to obtain the turgid weight, and finally dried at 70 °C until constant dry weight was achieved. Then, RWC was calculated as [(fresh weight-dry weight)/(turgid weight-dry weight) × 100].
For submergence stress, the submergence tolerance coefficient (STC) was adapted from the coefficient applied for waterlogging [70]. Briefly, the length of shoots and roots of plants exposed to submergence or control conditions were noted at the end of experiments. The plants were then dried at 70 °C until a constant dry weight was reached. The sums of length and dry weight of plants exposed to the stress and control conditions were used to calculate the STC as [(sum of treated plants)/(sum of control plants)]. RWC and STC measurements were performed in three biological replicates, each consisting of six plants per treatment.

4.4. Total RNA Extraction and Construction of Small RNA Libraries

Total RNA was extracted from treated and control samples (three biological replicates, each consisting of a pool of six plants) with TRIzol reagent (Life Technologies, Carlsbad, CA, USA) according to the manufacturer’s instructions. Total RNA was quantified using a Nanodrop ND-1000 Spectrophotometer (Thermo Fisher Scientific Inc., Waltham, MA, USA) and visualized by 1% agarose gel electrophoresis under denaturing conditions. Total RNA from triplicates of the same treatment was pooled (4 μg from each biological replicate in a total of 12 μg per treatment). The small RNA fraction (20–30 nt) was purified from total RNA in a denaturing 15% acrylamide gel. Next, a pre-adenylated adapter (linker1, Integrated DNA Technologies, San Diego, CA, USA) was ligated to the 3’ end of small RNAs with T4 RNA ligase II (New England Biolabs, Ipswich, MA, USA) in the absence of ATP. Then, the RNA adapter Illumina RA5 was ligated to the 5´ end of small RNAs with T4 RNA ligase 1 (Ambion, Austin, TX, USA). The ligated products were used as templates in a cDNA synthesis reaction with the primer RT Bridge and the retrotranscriptase MMuLV (Thermo Fisher Scientific Inc., Waltham, MA, USA), followed by a 15-cycle PCR reaction with Phusion High-Fidelity DNA polymerase (New England Biolabs). The PCR products were purified, quantified, and sequenced in the Illumina HiSEQ2500 at the Genomics Core of the University of California-Riverside (Riverside, CA, USA). Adapters and primers are listed in Table S7.

4.5. Bioinformatics Analysis of Small RNA Libraries

Raw sequencing datasets from small RNA libraries were processed with the CLC Genomics Workbench 12 (QIAGEN, Hilden, Germany) to remove adapters, low-quality reads, and reads shorter than 18 nt and longer than 35 nt. Clean reads were mapped to the B73 RefGen_v5 maize genome [71] and classified into the different non-coding RNA categories (rRNA, tRNA, snRNA, snoRNA, lncRNA, and repeats). In addition, clean reads were mapped to the maize chloroplast and mitochondria DNA (B73 RefGen_v4). The clean reads were searched against the plant miRNAs deposited in the miRBASE (Release 22.1) [31], allowing (1) zero mismatches, or (2) up to two mismatches to detect polymorphisms. miRNA sequence counts were normalized as TPM relative to the total clean read counts in the corresponding library. Individual miRNA sequences were then analyzed for differential expression. Fold change values were calculated as log2 (normalized counts in treatment/normalized counts in control) (Table S1). Only miRNA sequences with raw counts ≥ 50 in at least one of the libraries were considered for differential expression analysis, and fold change values ≥1 or ≤−1 were considered as upregulated or downregulated, respectively.
miRNA sequences with a high level of similarity were then grouped. All sequences that differed from the mature sequence by one or two nucleotides in the 5′ or 3′ ends were considered positional variants and were then grouped with the mature sequence to constitute a representative miRNA group. The sum of reads and fold change values were recalculated for grouped miRNA sequences (Table S2).
Small RNA sequencing data were deposited into the NCBI/GEO database with accession number GSE155050.

4.6. Assessment of miRNA Differential Expression by Quantitative RT-PCR (qRT-PCR)

A modified stem-loop method based on that previously described [72,73] was used to assess miRNA expression. Briefly, total RNA from treated and control plants were treated with DNase I (Thermo Fisher Scientific, Inc, Waltham, MA, USA.). Subsequently, 0.1 μg of treated RNA was used in 10 μL-cDNA synthesis reactions using the TaqMan miRNA reverse transcription kit (Applied Biosystems, Foster City, CA, USA) and miRNA specific stem-loop primers. For quantitative PCR, 20 μL-PCR reactions were prepared using Maxima Probe qPCR Master Mix (Thermo Fisher Scientific Inc, Waltham, MA, USA) with a miRNA-specific direct primer, PCR stem-loop reverse primer, and a universal fluorescein amidite (FAM)-labeled probe, in a PikoReal real-time thermal cycler (Thermo Fisher Scientific Inc.). Differential expression values were calculated using the ΔΔCt method [74] with the use of miR166c as an internal control to normalize the relative abundance of each miRNA. All primers and probes used are listed in Table S7.

4.7. Prediction of Target Genes

The psRNATarget analysis server (http://plantgrn.noble.org/psRNATarget/) was used for identification of target genes of miRNAs responsive to submergence and drought [75]. DPMIND [76] was used to search for miRNA targets in nine maize degradomes.

Supplementary Materials

The following are available online at https://www.mdpi.com/2223-7747/9/10/1367/s1, Figure S1: PHRED score distribution of small RNA reads in maize and teosinte libraries; Figure S2: Size distribution of small RNA reads in maize and teosinte libraries; Table S1: Known miRNA sequences expressed in maize and teosinte libraries; Table S2: Grouped sequences of known microRNAs expressed in maize and teosinte libraries; Table S3: Sequence variants of microRNAs expressed in maize and teosinte libraries; Table S4: The most abundant reads that mapped to NADH dehydrogenase subunit 2 from the Zea mays organelle DNA; Table S5: The most abundant reads that mapped to the trnL-trnF intergenic spacer and tRNA-Phe (trnF) gene from the maize chloroplast; Table S6: List of predicted targets of submergence and drought-responsive miRNAs; Table S7: List of adapters and primers used in this study.

Author Contributions

Conceptualization, R.L. and B.E.B.-F.; formal analysis, E.B.S.-G., J.F.P.-B. and B.E.B.-F.; funding acquisition, J.M.P.-C., R.L. and B.E.B.-F.; investigation, E.B.S.-G., J.M.P.-C., R.L. and B.E.B.-F.; methodology, E.B.S.-G., J.F.P.-B. and A.A.H.-H.; resources, J.M.P.-C. and B.E.B.-F.; supervision, B.E.B.-F.; writing—original draft, B.E.B.-F.; writing—review & editing, E.B.S.-G., A.A.H.-H., J.M.P.-C. and R.L. All authors have read and agreed to the published version of the manuscript.

Funding

This material is based upon work supported by a grant from the University of California Institute for Mexico and the United States (UC MEXUS) and the Consejo Nacional de Ciencia y Tecnología de México (CONACYT) to B.E.B.-F. and R.L., grants from Consejo Nacional de Ciencia y Tecnología CONACYT—México (Ciencia Básica 287137) to J.M.P.-C., and (Ciencia Básica 169619) to B.E.B.-F., CONACYT postdoctoral fellowship to E.B.S.-G. and A.A.H.-H., and CONACYT postgraduate fellowship to J.F.P.-B.

Acknowledgments

We wish to thank Axel Thiessen-Favier (CINVESTAV-Irapuato; in memoriam) for providing maize seeds; Flavio Aragon Cuevas (INIFAP Oaxaca Valles Centrales) for providing teosinte seeds; Mikael Roose (UCR) for academic advice and administrative assistance; Glenn Hicks and John Robert Weger (Genomics Core, UCR) for technical advice and performing sRNA sequencing; Fabiola Hernández, L.E. María de Jesús Hernández, L.C. Yenny Reyes Roque and Francisco Campos (Universidad del Papaloapan) for administrative assistance.

Conflicts of Interest

The authors declare no conflict of interest. The sponsors had no role in the design, execution, interpretation, or writing of the study.

References

  1. He, X.; Pan, M.; Wei, Z.; Wood, E.F.; Sheffield, J. A Global Drought and Flood Catalogue from 1950 to 2016. Bull. Am. Meteorol. Soc. 2020, 101, E508–E535. [Google Scholar] [CrossRef] [Green Version]
  2. Lesk, C.; Rowhani, P.; Ramankutty, N. Influence of Extreme Weather Disasters on Global Crop Production. Nature 2016, 529, 84–87. [Google Scholar] [CrossRef]
  3. Nguyen, L.T.T.; Osanai, Y.; Anderson, I.C.; Bange, M.P.; Tissue, D.T.; Singh, B.K. Flooding and Prolonged Drought Have Differential Legacy Impacts on Soil Nitrogen Cycling, Microbial Communities and Plant Productivity. Plant Soil 2018, 431, 371–387. [Google Scholar] [CrossRef]
  4. Voesenek, L.A.C.J.; Bailey-Serres, J. Flood Adaptive Traits and Processes: An Overview. New Phytol. 2015, 206, 57–73. [Google Scholar] [CrossRef] [PubMed]
  5. Fukao, T.; Barrera-Figueroa, B.E.; Juntawong, P.; Peña-Castro, J.M. Submergence and Waterlogging Stress in Plants: A Review Highlighting Research Opportunities and Understudied Aspects. Front. Plant Sci. 2019, 10, 340. [Google Scholar] [CrossRef] [PubMed]
  6. Loreti, E.; Valeri, M.C.; Novi, G.; Perata, P. Gene Regulation and Survival under Hypoxia Requires Starch Availability and Metabolism. Plant Physiol. 2018, 176, 1286–1298. [Google Scholar] [CrossRef] [PubMed]
  7. Bailey-Serres, J.; Lee, S.C.; Brinton, E. Waterproofing Crops: Effective Flooding Survival Strategies. Plant Physiol. 2012, 160, 1698–1709. [Google Scholar] [CrossRef] [Green Version]
  8. Osakabe, Y.; Osakabe, K.; Shinozaki, K.; Tran, L.-S.P. Response of Plants to Water Stress. Front. Plant Sci. 2014, 5. [Google Scholar] [CrossRef] [Green Version]
  9. Fukao, T.; Yeung, E.; Bailey-Serres, J. The Submergence Tolerance Regulator SUB1A Mediates Crosstalk between Submergence and Drought Tolerance in Rice. Plant Cell 2011, 23, 412–427. [Google Scholar] [CrossRef] [Green Version]
  10. Liang, S.; Xiong, W.; Yin, C.; Xie, X.; Jin, Y.; Zhang, S.; Yang, B.; Ye, G.; Chen, S.; Luan, W. Overexpression of OsARD1 Improves Submergence, Drought, and Salt Tolerances of Seedling Through the Enhancement of Ethylene Synthesis in Rice. Front. Plant Sci. 2019, 10, 1088. [Google Scholar] [CrossRef] [PubMed]
  11. Bin Rahman, A.N.M.R.; Zhang, J. Flood and Drought Tolerance in Rice: Opposite but May Coexist. Food Energy Secur. 2016, 5, 76–88. [Google Scholar] [CrossRef] [Green Version]
  12. Amissah, L.; Mohren, G.M.J.; Kyereh, B.; Agyeman, V.K.; Poorter, L. Rainfall Seasonality and Drought Performance Shape the Distribution of Tropical Tree Species in Ghana. Ecol. Evol. 2018, 8, 8582–8597. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. He, X.; Sheffield, J. Lagged Compound Occurrence of Droughts and Pluvials Globally over the Past Seven Decades. Geophys. Res. Lett. 2020. [Google Scholar] [CrossRef]
  14. Kim, W.; Iizumi, T.; Nishimori, M. Global Patterns of Crop Production Losses Associated with Droughts from 1983 to 2009. J. Appl. Meteorol. Climatol. 2019, 58, 1233–1244. [Google Scholar] [CrossRef]
  15. Li, Y.; Guan, K.; Schnitkey, G.D.; DeLucia, E.; Peng, B. Excessive Rainfall Leads to Maize Yield Loss of a Comparable Magnitude to Extreme Drought in the United States. Glob. Chang. Biol. 2019, 25, 2325–2337. [Google Scholar] [CrossRef] [Green Version]
  16. González, J.D.J.S.; Corral, J.A.R.; García, G.M.; Ojeda, G.R.; Larios, L.D.L.C.; Holland, J.B.; Medrano, R.M.; Romero, G.E.G. Ecogeography of Teosinte. PLoS ONE 2018, 13, e0192676. [Google Scholar] [CrossRef] [Green Version]
  17. Matsuoka, Y.; Vigouroux, Y.; Goodman, M.M.; Sanchez, G.J.; Buckler, E.; Doebley, J. A Single Domestication for Maize Shown by Multilocus Microsatellite Genotyping. Proc. Natl. Acad. Sci. USA 2002, 99, 6080–6084. [Google Scholar] [CrossRef] [Green Version]
  18. Hufford, M.B.; Xu, X.; van Heerwaarden, J.; Pyhäjärvi, T.; Chia, J.-M.; Cartwright, R.A.; Elshire, R.J.; Glaubitz, J.C.; Guill, K.E.; Kaeppler, S.M.; et al. Comparative Population Genomics of Maize Domestication and Improvement. Nat. Genet. 2012, 44, 808–811. [Google Scholar] [CrossRef] [Green Version]
  19. Liu, J.; Fernie, A.R.; Yan, J. The Past, Present, and Future of Maize Improvement: Domestication, Genomics, and Functional Genomic Routes toward Crop Enhancement. Plant Commun. 2020, 1, 100010. [Google Scholar] [CrossRef]
  20. Mano, Y.; Omori, F. Breeding for Flooding Tolerant Maize Using “Teosinte” as a Germplasm Resource. Plant Roots 2007, 1, 17–21. [Google Scholar] [CrossRef] [Green Version]
  21. Shriram, V.; Kumar, V.; Devarumath, R.M.; Khare, T.S.; Wani, S.H. MicroRNAs As Potential Targets for Abiotic Stress Tolerance in Plants. Front. Plant Sci. 2016, 7, 817. [Google Scholar] [CrossRef] [PubMed]
  22. Song, X.; Li, Y.; Cao, X.; Qi, Y. MicroRNAs and Their Regulatory Roles in Plant–Environment Interactions. Annu. Rev. Plant Biol. 2019, 70, 489–525. [Google Scholar] [CrossRef] [PubMed]
  23. Moldovan, D.; Spriggs, A.; Yang, J.; Pogson, B.J.; Dennis, E.S.; Wilson, I.W. Hypoxia-Responsive MicroRNAs and Trans-Acting Small Interfering RNAs in Arabidopsis. J. Exp. Bot. 2010, 61, 165–177. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Pegler, J.; Oultram, J.; Grof, C.; Eamens, A. Profiling the Abiotic Stress Responsive MicroRNA Landscape of Arabidopsis Thaliana. Plants 2019, 8, 58. [Google Scholar] [CrossRef] [Green Version]
  25. Jeong, D.-H.; Schmidt, S.A.; Rymarquis, L.A.; Park, S.; Ganssmann, M.; German, M.A.; Accerbi, M.; Zhai, J.; Fahlgren, N.; Fox, S.E.; et al. Parallel Analysis of RNA Ends Enhances Global Investigation of MicroRNAs and Target RNAs of Brachypodium Distachyon. Genome Biol. 2013, 14, R145. [Google Scholar] [CrossRef] [Green Version]
  26. Bertolini, E.; Verelst, W.; Horner, D.S.; Gianfranceschi, L.; Piccolo, V.; Inzé, D.; Pè, M.E.; Mica, E. Addressing the Role of MicroRNAs in Reprogramming Leaf Growth during Drought Stress in Brachypodium Distachyon. Mol. Plant 2013, 6, 423–443. [Google Scholar] [CrossRef] [Green Version]
  27. Barrera-Figueroa, B.E.; Gao, L.; Wu, Z.; Zhou, X.; Zhu, J.; Jin, H.; Liu, R.; Zhu, J.-K. High Throughput Sequencing Reveals Novel and Abiotic Stress-Regulated MicroRNAs in the Inflorescences of Rice. BMC Plant Biol. 2012, 12, 132. [Google Scholar] [CrossRef] [Green Version]
  28. Jin, Q.; Xu, Y.; Mattson, N.; Li, X.; Wang, B.; Zhang, X.; Jiang, H.; Liu, X.; Wang, Y.; Yao, D. Identification of Submergence-Responsive MicroRNAs and Their Targets Reveals Complex MiRNA-Mediated Regulatory Networks in Lotus (Nelumbo Nucifera Gaertn). Front. Plant Sci. 2017, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Liu, X.; Zhang, X.; Sun, B.; Hao, L.; Liu, C.; Zhang, D.; Tang, H.; Li, C.; Li, Y.; Shi, Y.; et al. Genome-Wide Identification and Comparative Analysis of Drought-Related MicroRNAs in Two Maize Inbred Lines with Contrasting Drought Tolerance by Deep Sequencing. PLoS ONE 2019, 14, e0219176. [Google Scholar] [CrossRef] [PubMed]
  30. Liu, Z.; Kumari, S.; Zhang, L.; Zheng, Y.; Ware, D. Characterization of MiRNAs in Response to Short-Term Waterlogging in Three Inbred Lines of Zea Mays. PLoS ONE 2012, 7, e39786. [Google Scholar] [CrossRef]
  31. Kozomara, A.; Birgaoanu, M.; Griffiths-Jones, S. MiRBase: From MicroRNA Sequences to Function. Nucleic Acids Res. 2019, 47, D155–D162. [Google Scholar] [CrossRef] [PubMed]
  32. Li, J.; Fu, F.; An, M.; Zhou, S.; She, Y.; Li, W. Differential Expression of MicroRNAs in Response to Drought Stress in Maize. J. Integr. Agric. 2013, 12, 1414–1422. [Google Scholar] [CrossRef]
  33. Seeve, C.M.; Sunkar, R.; Zheng, Y.; Liu, L.; Liu, Z.; McMullen, M.; Nelson, S.; Sharp, R.E.; Oliver, M.J. Water-Deficit Responsive MicroRNAs in the Primary Root Growth Zone of Maize. BMC Plant Biol. 2019, 19, 447. [Google Scholar] [CrossRef] [Green Version]
  34. Zhang, Z.; Wei, L.; Zou, X.; Tao, Y.; Liu, Z.; Zheng, Y. Submergence-Responsive MicroRNAs Are Potentially Involved in the Regulation of Morphological and Metabolic Adaptations in Maize Root Cells. Ann. Bot. 2008, 102, 509–519. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Zhai, L.; Liu, Z.; Zou, X.; Jiang, Y.; Qiu, F.; Zheng, Y.; Zhang, Z. Genome-Wide Identification and Analysis of MicroRNA Responding to Long-Term Waterlogging in Crown Roots of Maize Seedlings. Physiol. Plant. 2013, 147, 181–193. [Google Scholar] [CrossRef]
  36. Axtell, M.J.; Meyers, B.C. Revisiting Criteria for Plant MicroRNA Annotation in the Era of Big Data. Plant Cell 2018, 30, 272–284. [Google Scholar] [CrossRef] [Green Version]
  37. Meng, X.; Li, L.; Narsai, R.; De Clercq, I.; Whelan, J.; Berkowitz, O. Mitochondrial Signalling Is Critical for Acclimation and Adaptation to Flooding in Arabidopsis Thaliana. Plant J. 2020, 103, 227–247. [Google Scholar] [CrossRef]
  38. Mommer, L.; Pons, T.L.; Wolters-Arts, M.; Venema, J.H.; Visser, E.J.W. Submergence-Induced Morphological, Anatomical, and Biochemical Responses in a Terrestrial Species Affect Gas Diffusion Resistance and Photosynthetic Performance. Plant Physiol. 2005, 139, 497–508. [Google Scholar] [CrossRef] [Green Version]
  39. Schmidt, R.R.; Weits, D.A.; Feulner, C.F.J.; van Dongen, J.T. Oxygen Sensing and Integrative Stress Signaling in Plants. Plant Physiol. 2018, 176, 1131–1142. [Google Scholar] [CrossRef] [Green Version]
  40. Hou, Y.; Jiang, F.; Zheng, X.; Wu, Z. Identification and Analysis of Oxygen Responsive MicroRNAs in the Root of Wild Tomato (S. Habrochaites). BMC Plant Biol. 2019, 19, 100. [Google Scholar] [CrossRef] [Green Version]
  41. Voesenek, L.A.C.J. Interactions Between Plant Hormones Regulate Submergence-Induced Shoot Elongation in the Flooding-Tolerant Dicot Rumex Palustris. Ann. Bot. 2003, 91, 205–211. [Google Scholar] [CrossRef] [Green Version]
  42. Chuck, G.S.; Tobias, C.; Sun, L.; Kraemer, F.; Li, C.; Dibble, D.; Arora, R.; Bragg, J.N.; Vogel, J.P.; Singh, S.; et al. Overexpression of the Maize Corngrass1 MicroRNA Prevents Flowering, Improves Digestibility, and Increases Starch Content of Switchgrass. Proc. Natl. Acad. Sci. USA 2011, 108, 17550–17555. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Mao, H.-D.; Yu, L.-J.; Li, Z.-J.; Yan, Y.; Han, R.; Liu, H.; Ma, M. Genome-Wide Analysis of the SPL Family Transcription Factors and Their Responses to Abiotic Stresses in Maize. Plant Gene 2016, 6, 1–12. [Google Scholar] [CrossRef] [Green Version]
  44. Wu, G. Temporal Regulation of Shoot Development in Arabidopsis Thaliana by MiR156 and Its Target SPL3. Development 2006, 133, 3539–3547. [Google Scholar] [CrossRef] [Green Version]
  45. Guo, H.-S.; Xie, Q.; Fei, J.-F.; Chua, N.-H. MicroRNA Directs MRNA Cleavage of the Transcription Factor NAC1 to Downregulate Auxin Signals for Arabidopsis Lateral Root Development. Plant Cell 2005, 17, 1376–1386. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Zhang, J.; Long, Y.; Xue, M.; Xiao, X.; Pei, X. Identification of MicroRNAs in Response to Drought in Common Wild Rice (Oryza Rufipogon Griff.) Shoots and Roots. PLoS ONE 2017, 12, e0170330. [Google Scholar] [CrossRef] [PubMed]
  47. Singh, A.; Roy, S.; Singh, S.; Das, S.S.; Gautam, V.; Yadav, S.; Kumar, A.; Singh, A.; Samantha, S.; Sarkar, A.K. Phytohormonal Crosstalk Modulates the Expression of MiR166/165s, Target Class III HD-ZIPs, and KANADI Genes during Root Growth in Arabidopsis Thaliana. Sci. Rep. 2017, 7, 3408. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Li, S.-B.; Xie, Z.-Z.; Hu, C.-G.; Zhang, J.-Z. A Review of Auxin Response Factors (ARFs) in Plants. Front. Plant Sci. 2016, 7. [Google Scholar] [CrossRef] [Green Version]
  49. Li, W.-X.; Oono, Y.; Zhu, J.; He, X.-J.; Wu, J.-M.; Iida, K.; Lu, X.-Y.; Cui, X.; Jin, H.; Zhu, J.-K. The Arabidopsis NFYA5 Transcription Factor Is Regulated Transcriptionally and Posttranscriptionally to Promote Drought Resistance. Plant Cell 2008, 20, 2238–2251. [Google Scholar] [CrossRef] [Green Version]
  50. Luan, M.; Xu, M.; Lu, Y.; Zhang, Q.; Zhang, L.; Zhang, C.; Fan, Y.; Lang, Z.; Wang, L. Family-Wide Survey of MiR169s and NF-YAs and Their Expression Profiles Response to Abiotic Stress in Maize Roots. PLoS ONE 2014, 9, e91369. [Google Scholar] [CrossRef] [Green Version]
  51. Xu, M.Y.; Zhang, L.; Li, W.W.; Hu, X.L.; Wang, M.-B.; Fan, Y.L.; Zhang, C.Y.; Wang, L. Stress-Induced Early Flowering Is Mediated by MiR169 in Arabidopsis Thaliana. J. Exp. Bot. 2014, 65, 89–101. [Google Scholar] [CrossRef] [Green Version]
  52. Sorin, C.; Declerck, M.; Christ, A.; Blein, T.; Ma, L.; Lelandais-Brière, C.; Njo, M.F.; Beeckman, T.; Crespi, M.; Hartmann, C. A MiR169 Isoform Regulates Specific NF-YA Targets and Root Architecture in Arabidopsis. New Phytol. 2014, 202, 1197–1211. [Google Scholar] [CrossRef]
  53. Zhou, M.; Luo, H. Role of MicroRNA319 in Creeping Bentgrass Salinity and Drought Stress Response. Plant Signal. Behav. 2014, 9, e28700. [Google Scholar] [CrossRef] [Green Version]
  54. Wang, L.; Gu, X.; Xu, D.; Wang, W.; Wang, H.; Zeng, M.; Chang, Z.; Huang, H.; Cui, X. MiR396-Targeted AtGRF Transcription Factors Are Required for Coordination of Cell Division and Differentiation during Leaf Development in Arabidopsis. J. Exp. Bot. 2011, 62, 761–773. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Beltramino, M.; Ercoli, M.F.; Debernardi, J.M.; Goldy, C.; Rojas, A.M.L.; Nota, F.; Alvarez, M.E.; Vercruyssen, L.; Inzé, D.; Palatnik, J.F.; et al. Robust Increase of Leaf Size by Arabidopsis Thaliana GRF3-like Transcription Factors under Different Growth Conditions. Sci. Rep. 2018, 8, 13447. [Google Scholar] [CrossRef] [PubMed]
  56. Gocal, G.F.W.; Sheldon, C.C.; Gubler, F.; Moritz, T.; Bagnall, D.J.; MacMillan, C.P.; Li, S.F.; Parish, R.W.; Dennis, E.S.; Weigel, D.; et al. GAMYB-like Genes, Flowering, and Gibberellin Signaling in Arabidopsis. Plant Physiol. 2001, 127, 1682–1693. [Google Scholar] [CrossRef]
  57. Rivera-Contreras, I.K.; Zamora-Hernández, T.; Huerta-Heredia, A.A.; Capataz-Tafur, J.; Barrera-Figueroa, B.E.; Juntawong, P.; Peña-Castro, J.M. Transcriptomic Analysis of Submergence-Tolerant and Sensitive Brachypodium Distachyon Ecotypes Reveals Oxidative Stress as a Major Tolerance Factor. Sci. Rep. 2016, 6, 27686. [Google Scholar] [CrossRef] [Green Version]
  58. Sunkar, R.; Kapoor, A.; Zhu, J.-K. Posttranscriptional Induction of Two Cu/Zn Superoxide Dismutase Genes in Arabidopsis Is Mediated by Downregulation of MiR398 and Important for Oxidative Stress Tolerance. Plant Cell 2006, 18, 2051–2065. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. De la Rosa, C.; Covarrubias, A.A.; Reyes, J.L. A Dicistronic Precursor Encoding MiR398 and the Legume-Specific MiR2119 Coregulates CSD1 and ADH1 MRNAs in Response to Water Deficit: MiR398 and MiR2119 Co-Regulate CSD1 and ADH1 MRNAs. Plant Cell Environ. 2019, 42, 133–144. [Google Scholar] [CrossRef]
  60. Loreti, E.; Betti, F.; Ladera-Carmona, M.J.; Fontana, F.; Novi, G.; Valeri, M.C.; Perata, P. ARGONAUTE1 and ARGONAUTE4 Regulate Gene Expression and Hypoxia Tolerance. Plant Physiol. 2020, 182, 287–300. [Google Scholar] [CrossRef] [Green Version]
  61. Gong, S.; Ding, Y.; Huang, S.; Zhu, C. Identification of MiRNAs and Their Target Genes Associated with Sweet Corn Seed Vigor by Combined Small RNA and Degradome Sequencing. J. Agric. Food Chem. 2015, 63, 5485–5491. [Google Scholar] [CrossRef]
  62. Abdel-Ghany, S.E.; Pilon, M. MicroRNA-Mediated Systemic Down-Regulation of Copper Protein Expression in Response to Low Copper Availability in Arabidopsis. J. Biol. Chem. 2008, 283, 15932–15945. [Google Scholar] [CrossRef] [Green Version]
  63. Sun, Q.; Liu, X.; Yang, J.; Liu, W.; Du, Q.; Wang, H.; Fu, C.; Li, W.-X. MicroRNA528 Affects Lodging Resistance of Maize by Regulating Lignin Biosynthesis under Nitrogen-Luxury Conditions. Mol. Plant 2018, 11, 806–814. [Google Scholar] [CrossRef] [Green Version]
  64. Zhang, H.; Zhao, X.; Li, J.; Cai, H.; Deng, X.W.; Li, L. MicroRNA408 Is Critical for the HY5-SPL7 Gene Network That Mediates the Coordinated Response to Light and Copper. Plant Cell 2014, 26, 4933–4953. [Google Scholar] [CrossRef] [Green Version]
  65. Zhao, X.Y.; Hong, P.; Wu, J.Y.; Chen, X.B.; Ye, X.G.; Pan, Y.Y.; Wang, J.; Zhang, X.S. The Tae-MiR408-Mediated Control of TaTOC1 Genes Transcription Is Required for the Regulation of Heading Time in Wheat. Plant Physiol. 2016, 170, 1578–1594. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Yang, R.; Li, P.; Mei, H.; Wang, D.; Sun, J.; Yang, C.; Hao, L.; Cao, S.; Chu, C.; Hu, S.; et al. Fine-Tuning of MiR528 Accumulation Modulates Flowering Time in Rice. Mol. Plant 2019, 12, 1103–1113. [Google Scholar] [CrossRef] [Green Version]
  67. Zhang, L.; Chia, J.-M.; Kumari, S.; Stein, J.C.; Liu, Z.; Narechania, A.; Maher, C.A.; Guill, K.; McMullen, M.D.; Ware, D. A Genome-Wide Characterization of MicroRNA Genes in Maize. PLoS Genet. 2009, 5, e1000716. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Yeung, E.; Bailey-Serres, J.; Sasidharan, R. After The Deluge: Plant Revival Post-Flooding. Trends Plant. Sci. 2019, 24, 443–454. [Google Scholar] [CrossRef] [PubMed]
  69. Barrs, H.; Weatherley, P. A Re-Examination of the Relative Turgidity Technique for Estimating Water Deficits in Leaves. Aust. J. Biol. Sci. 1962, 15, 413. [Google Scholar] [CrossRef] [Green Version]
  70. Liu, Y.; Tang, B.; Zheng, Y.; Ma, K.; Xu, S.; Qiu, F. Screening Methods for Waterlogging Tolerance at Maize (Zea mays L.) Seedling Stage. Agric. Sci. China 2010, 9, 362–369. [Google Scholar] [CrossRef]
  71. Portwood, J.L.; Woodhouse, M.R.; Cannon, E.K.; Gardiner, J.M.; Harper, L.C.; Schaeffer, M.L.; Walsh, J.R.; Sen, T.Z.; Cho, K.T.; Schott, D.A.; et al. MaizeGDB 2018: The Maize Multi-Genome Genetics and Genomics Database. Nucleic Acids Res. 2019, 47, D1146–D1154. [Google Scholar] [CrossRef]
  72. Chen, C. Real-Time Quantification of MicroRNAs by Stem-Loop RT-PCR. Nucleic Acids Res. 2005, 33, e179. [Google Scholar] [CrossRef]
  73. Jung, U.; Jiang, X.; Kaufmann, S.H.E.; Patzel, V. A Universal TaqMan-Based RT-PCR Protocol for Cost-Efficient Detection of Small Noncoding RNA. RNA 2013, 19, 1864–1873. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  75. Dai, X.; Zhuang, Z.; Zhao, P.X. PsRNATarget: A Plant Small RNA Target Analysis Server (2017 Release). Nucleic Acids Res. 2018, 46, W49–W54. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Fei, Y.; Wang, R.; Li, H.; Liu, S.; Zhang, H.; Huang, J. DPMIND: Degradome-Based Plant MiRNA–Target Interaction and Network Database. Bioinformatics 2018, 34, 1618–1620. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Effects of submergence and drought on maize and teosinte plants. (A) Plants exposed to submergence for 2, 4, and 6 days, and (B) plants exposed to drought for 2, 4, 6, and 8 days. Response of teosinte and maize plants to (C) submergence for 48 h, expressed as submergence tolerance coefficient (STC), and (D) drought for 6 days and recovery 24 h after rehydration, expressed as relative water content (RWC) in the leaves. Data are means, with error bars representing ±SD of three biological replicates, each replicate consisting of six plants. Different letters above the error bars indicate statistically significant differences between samples (p < 0.05) in a Student’s t-test.
Figure 1. Effects of submergence and drought on maize and teosinte plants. (A) Plants exposed to submergence for 2, 4, and 6 days, and (B) plants exposed to drought for 2, 4, 6, and 8 days. Response of teosinte and maize plants to (C) submergence for 48 h, expressed as submergence tolerance coefficient (STC), and (D) drought for 6 days and recovery 24 h after rehydration, expressed as relative water content (RWC) in the leaves. Data are means, with error bars representing ±SD of three biological replicates, each replicate consisting of six plants. Different letters above the error bars indicate statistically significant differences between samples (p < 0.05) in a Student’s t-test.
Plants 09 01367 g001
Figure 2. Venn diagrams showing differentially expressed miRNA reads that overlap in maize and teosinte under (A) submergence and (B) drought treatments. Plant miRNAs overlapping in submergence and drought treatments in (C) maize, and (D) teosinte. Plant miRNAs overlapping in teosinte under single and alternated treatments of (E) drought followed by submergence, and (F) submergence followed by drought. Numbers are based on fold change (FC) (log2) values (treatment/control). Upward arrows represent upregulated miRNAs (FC ≥ 1) and downward arrows represent downregulated miRNAs (FC ≤ −1).
Figure 2. Venn diagrams showing differentially expressed miRNA reads that overlap in maize and teosinte under (A) submergence and (B) drought treatments. Plant miRNAs overlapping in submergence and drought treatments in (C) maize, and (D) teosinte. Plant miRNAs overlapping in teosinte under single and alternated treatments of (E) drought followed by submergence, and (F) submergence followed by drought. Numbers are based on fold change (FC) (log2) values (treatment/control). Upward arrows represent upregulated miRNAs (FC ≥ 1) and downward arrows represent downregulated miRNAs (FC ≤ −1).
Plants 09 01367 g002
Figure 3. qRT-PCR validation of miRNAs response to submergence, drought, and alternated stress by SL-RTPCR for (A) miR156, (B) miR159, (C) miR167, (D) miR396, (E) miR398, (F) miR408, and (G) miR528. Relative fold change values were calculated as ∆∆Ct for SL-RTPCR (blue bars) using miR166c as the constitutive control or as the treatment/control ratio from normalized reads for sRNA-SEQ data (black bars). Results are mean with error bars representing ±SD of three biological replicates, each replicate consisting of a pool of six plants and three technical replicates. Different letters above the error bars indicate statistically significant differences between samples (p < 0.05).
Figure 3. qRT-PCR validation of miRNAs response to submergence, drought, and alternated stress by SL-RTPCR for (A) miR156, (B) miR159, (C) miR167, (D) miR396, (E) miR398, (F) miR408, and (G) miR528. Relative fold change values were calculated as ∆∆Ct for SL-RTPCR (blue bars) using miR166c as the constitutive control or as the treatment/control ratio from normalized reads for sRNA-SEQ data (black bars). Results are mean with error bars representing ±SD of three biological replicates, each replicate consisting of a pool of six plants and three technical replicates. Different letters above the error bars indicate statistically significant differences between samples (p < 0.05).
Plants 09 01367 g003
Figure 4. Hypothetical model of the regulation of miRNAs in the response to submergence and drought in maize and teosinte. The importance of antioxidant enzymatic activity, copper homeostasis, and control of growth and development as the main lines of response regulated by miRNAs in both plants is highlighted. Interrupted lines represent overlap between plants and/or treatments. Green, downregulated miRNAs; orange, upregulated miRNAs. GRF: growth-regulating factor; GAMYB: gibberellic acid-MYB transcription factor; ARF: auxin-responsive factor; NFYA: nuclear transcription factor Y subunit alpha; PLC: plastocyanin; LAC: laccase; SOD: superoxide dismutase; SBP1: selenium-binding protein1; bZIP: basic leucine zipper transcription factor; SPL: squamosa promoter-binding protein like; NAC: NAM-ATAF-CUC domain transcription factor; TCP: teosinte branched1/cycloidea/proliferating cell nuclear antigen factor; SBP1: selenium-binding protein 1.
Figure 4. Hypothetical model of the regulation of miRNAs in the response to submergence and drought in maize and teosinte. The importance of antioxidant enzymatic activity, copper homeostasis, and control of growth and development as the main lines of response regulated by miRNAs in both plants is highlighted. Interrupted lines represent overlap between plants and/or treatments. Green, downregulated miRNAs; orange, upregulated miRNAs. GRF: growth-regulating factor; GAMYB: gibberellic acid-MYB transcription factor; ARF: auxin-responsive factor; NFYA: nuclear transcription factor Y subunit alpha; PLC: plastocyanin; LAC: laccase; SOD: superoxide dismutase; SBP1: selenium-binding protein1; bZIP: basic leucine zipper transcription factor; SPL: squamosa promoter-binding protein like; NAC: NAM-ATAF-CUC domain transcription factor; TCP: teosinte branched1/cycloidea/proliferating cell nuclear antigen factor; SBP1: selenium-binding protein 1.
Plants 09 01367 g004
Figure 5. Hypothetical model of the response to alternated stress in teosinte. The influence of the first stress and the consequences in miRNA responsiveness to the second stress are indicated. Green, downregulated miRNAs; orange, upregulated miRNAs.
Figure 5. Hypothetical model of the response to alternated stress in teosinte. The influence of the first stress and the consequences in miRNA responsiveness to the second stress are indicated. Green, downregulated miRNAs; orange, upregulated miRNAs.
Plants 09 01367 g005
Table 1. Summary of small RNA sequencing data from eight maize and teosinte libraries. Distribution of small RNAs in the corresponding non-coding RNA category.
Table 1. Summary of small RNA sequencing data from eight maize and teosinte libraries. Distribution of small RNAs in the corresponding non-coding RNA category.
MaizeTeosinte
ControlSubmergDroughtControlSubmergDroughtDrought-SubmergSubmerg-Drought
Raw16,139,35426,621,44831,642,21020,672,41328,656,75629,063,58817,997,03746,522,229
Clean14,897,05824,568,43929,744,16118,955,89226,296,92126,960,31315,760,25841,915,673
Maize14,870,94724,488,35929,696,34118,927,43426,217,67226,912,77715,727,75741,761,974
99.82%99.67%99.84%99.85%99.70%99.82%99.79%99.63%
mt87,401106,625134,420114,74173,980105,999106,451145,709
0.59%0.43%0.45%0.61%0.28%0.39%0.68%0.35%
cp182,414452,301459,718475,121444,434835,675973,9491,289,213
1.22%1.84%1.55%2.51%1.69%3.10%6.18%3.08%
mRNA195,418258,447342,833207,230202,191256,739238,774309,551
1.31%1.05%1.15%1.09%0.77%0.95%1.52%0.74%
rRNA288,257601,629778,148687,038591,868588,8191,056,2431,240,514
1.94%2.45%2.62%3.62%2.25%2.18%6.71%2.97%
tRNA103,451161,288200,294569,294480,975553,939953,7941,426,872
0.70%0.66%0.67%3.00%1.83%2.05%6.06%3.41%
snRNA516,999838,1511,071,157576,299868,349890,657421,4041,220,748
3.47%3.42%3.60%3.04%3.31%3.30%2.67%2.92%
snoRNA105,518168,088228,389141,006174,443190,297111,528256,675
0.70%68.00%0.77%0.74%0.66%0.70%0.70%0.61%
lncRNA258,027416,896519,449362,527530,505528,963231,375719,402
1.73%1.70%1.74%1.91%2.02%1.96%1.47%1.72%
repeats2,219,6563,623,1994,820,4555,064,7985,305,4815,343,3144,734,8799,267,896
14.92%14.79%16.23%26.75%20.23%19.85%30%22.19%
Maize397,473248,193774,680180,09198,923310,47847,836350,678
miRNAs2.67%1.01%2.60%0.95%0.38%1.15%0.30%0.84%
Other24,83721,52639,108649858627393647013,887
miRNAs0.17%0.09%0.13%0.03%0.02%0.03%0.04%0.03%
Submerg: Submergence; mt: mitochondria; cp: chloroplast; mRNA: messenger RNA; rRNA: ribosomal RNA; tRNA: transfer RNA; snRNA: small nuclear RNA; snoRNA: small nucleolar RNA; lncRNA: long non-coding RNA.
Table 2. Representative plant miRNAs responsive to submergence and drought in maize and teosinte.
Table 2. Representative plant miRNAs responsive to submergence and drought in maize and teosinte.
miRNA NameFold Change (FC) 1
MaizeTeosinte
SubmergenceDroughtSubmergenceDrought
miR156k−2.94−2.99−0.861.09
miR159ab−2.170.80−1.461.06
miR164e−4.24−1.49−0.291.20
miR166bd−4.39−3.79−2.77−1.49
miR167cdeg−2.99−2.48−1.62−0.33
miR169cr−3.13−0.95−1.050.11
miR319b−0.750.80−0.662.17
miR396cd−1.04−0.17−0.391.08
miR398ab−0.38−0.29−1.090.05
miR398b−2.83−0.94−0.690.29
miR408−1.34−0.51−1.78−1.76
miR408b−2.14−0.89−2.570.17
miR528ab−1.30−1.65−2.22−0.57
Constitutive
miR166c0.00−0.01−0.21−0.19
1 Fold change (treatment/control) values ≥ 1 indicate upregulated miRNAs (orange), and values ≤ −1 indicate downregulated miRNAs (green).
Table 3. Representative plant miRNAs responsive to single (submergence or drought) and alternated treatments (drought followed by submergence and submergence followed by drought) in teosinte plants.
Table 3. Representative plant miRNAs responsive to single (submergence or drought) and alternated treatments (drought followed by submergence and submergence followed by drought) in teosinte plants.
miRNA NameFold Change (FC) 1
SubmergenceDroughtDrought-SubmergenceSubmergence-Drought
miR156k−0.861.09−0.541.17
miR159ab−1.461.06−1.430.55
miR164e−0.291.200.061.09
miR166bd−2.77−1.49−1.95−0.77
miR167cdeg−1.62−0.33−1.870.71
miR169cr−1.050.11−1.600.11
miR319b−0.662.17−0.410.16
miR396cd−0.391.08−0.940.53
miR398ab−1.090.05−1.56−0.58
miR398b−0.690.29−1.92−0.57
miR408−1.78−1.76−0.880.55
miR408b−2.570.17−0.43−0.87
miR528ab−2.22−0.57−3.16−0.99
Constitutive
miR166c−0.21−0.19−0.27−0.19
1 Fold change (treatment/control) values ≥ 1 indicate upregulated miRNAs (orange), and values ≤ −1 indicate downregulated miRNAs (green).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sepúlveda-García, E.B.; Pulido-Barajas, J.F.; Huerta-Heredia, A.A.; Peña-Castro, J.M.; Liu, R.; Barrera-Figueroa, B.E. Differential Expression of Maize and Teosinte microRNAs under Submergence, Drought, and Alternated Stress. Plants 2020, 9, 1367. https://doi.org/10.3390/plants9101367

AMA Style

Sepúlveda-García EB, Pulido-Barajas JF, Huerta-Heredia AA, Peña-Castro JM, Liu R, Barrera-Figueroa BE. Differential Expression of Maize and Teosinte microRNAs under Submergence, Drought, and Alternated Stress. Plants. 2020; 9(10):1367. https://doi.org/10.3390/plants9101367

Chicago/Turabian Style

Sepúlveda-García, Edgar Baldemar, José Francisco Pulido-Barajas, Ariana Arlene Huerta-Heredia, Julián Mario Peña-Castro, Renyi Liu, and Blanca Estela Barrera-Figueroa. 2020. "Differential Expression of Maize and Teosinte microRNAs under Submergence, Drought, and Alternated Stress" Plants 9, no. 10: 1367. https://doi.org/10.3390/plants9101367

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop