Abstract
Phenotypic plasticity is crucial for plants to survive in changing environments. Discovering microRNAs, identifying their targets and further inferring microRNA functions in mediating plastic developmental responses to environmental changes have been a critical strategy for understanding the underlying molecular mechanisms of phenotypic plasticity. In this study, the dynamic expression patterns of microRNAs under contrasting hydrological habitats in the amphibious species Alternanthera philoxeroides were identified by time course expression profiling using high-throughput sequencing technology. A total of 128 known and 18 novel microRNAs were found to be differentially expressed under contrasting hydrological habitats. The microRNA:mRNA pairs potentially associated with plastic internode elongation were identified by integrative analysis of microRNA and mRNA expression profiles, and were validated by qRT-PCR and 5′ RLM-RACE. The results showed that both the universal microRNAs conserved across different plants and the unique microRNAs novelly identified in A. philoxeroides were involved in the responses to varied water regimes. The results also showed that most of the differentially expressed microRNAs were transiently up-/down-regulated at certain time points during the treatments. The fine-scale temporal changes in microRNA expression highlighted the importance of time-series sampling in identifying stress-responsive microRNAs and analyzing their role in stress response/tolerance.
Introduction
The ability of organisms to produce distinct phenotypes under different environmental conditions is referred to as phenotypic plasticity (; ; ), which is crucial for organisms to survive in a world that is constantly changing, especially when the environment changes quickly (). It is usually considered that plants have greater plasticity than animals for their sessile lifestyle (; ; ). Plants are forced to cope with local environmental fluctuations by plastic modifications in morphology, physiology, behavior or life history traits, to overcome the challenges from changing environments and increase their niche breadth (; ). Over the past two decades, the molecular and cellular mechanisms of environmental perception and gene regulation underlying various plastic responses in plants have been an explicit focus of the ecological developmental (Eco-Devo) research, with the primary goal of characterizing the developmental genetic pathways supporting the ability to respond plastically to environment changes (, ; ). Understanding the regulatory events informed by both external and internal environmental signals and identifying genetic elements implicated in environmentally contingent development will provide key insights to the immediate tolerance and potential evolutionary resilience of organisms to environmental fluctuations (, ; ).
Investigating the development of plastic traits in ecological context has revealed an association between environmentally contingent phenotypic expression and the environmental sensitivity of gene expression (). This environmental sensitivity of gene expression is, in turn, driven by the flexibility of epigenetic gene regulation in fluctuating environments (; ). Epigenetic variations triggered by environmental perturbations can lead to inducible phenotypic changes by altering transcriptional profiles and eventually altering ontogenetic trajectory in response to environmental variation, even in the complete absence of genetic variability (; ), and thus serve as a mechanistic link between genes and environment (; ). A suite of epigenetic mechanisms, including DNA methylation, histone modifications and small RNAs, have been shown to be able to control the temporal, spatial and abundance patterns of gene expression and RNA-translation under varied environments (; ). These epigenetic mechanisms may act separately or concomitantly (; ).
microRNAs, a type of endogenous small-non-coding RNAs, have emerged as key post-transcriptional regulators of gene expression through base pairing with their complementary mRNA targets (). Functional genomic studies have shown the involvement of plant microRNAs in a broad range of developmental processes (; ; ). miR165/166 have been proved to participate in SAM development, acting as mobile signal molecules for stem cell maintenance (). miR319 mediates the change of plant leaf shape via targeting TCP transcription factors (), while miR160 and miR393 play key roles in root architecture regulation via post-transcriptional modification of the auxin signal pathway genes (; ). Both miR156 and miR172 are involved in regulation of flowering time by age-dependent pathway (; ). Multiple microRNAs, including miR169, miR164, miR319, miR159, and miR167 function to control floral organ identity fate during flower development by regulating the spatial boundaries of expression of target genes, and to specify particular cell types during later stages of flower development (). It has also been revealed that a single miRNA may have diverse regulatory roles. For example, miR159 is involved in flowering time control, leaf development, seed size and shape determination, respectively (; ; ).
In addition to control various developmental processes, plant microRNAs also act as node molecules to coordinate the balance between plant development and environmental clues, playing critical roles in response to biotic and abiotic stresses (; ; ; ), microRNAs involved in responses to water-logging/flooding have been demonstrated in different plants. For instance, and identified microRNAs responding to waterlogging in maize roots. predicted the associated microRNAs under waterlogged conditions in sugarcane. found significant changes in the expression of seven conserved microRNA families and five novel microRNAs in response to flooding stress in Populus tomentosa. revealed the microRNA-mediated regulatory networks in lotus by identifying submergence-responsive microRNAs and their targets. Flooding exposes plants to severe hypoxia. Thus, microRNAs mediating hypoxic stress responses in plants were also thought to facilitate the survival and adaptation to the flooding condition by modifying gene expression, metabolic networks, and developmental processes (; ).
Alternanthera philoxeroides is an exotic amphibious weed that is native to South America but has now invaded into the temperate and tropical regions across the world (; ). A. philoxeroides rarely produces viable seeds in its introduced range, and propagates mainly via vegetative regeneration (; ; ). Plants produced by vegetative propagation are genetically identical. However, these plants can grow in diverse habitats, from purely aquatic environments (ponds, rivers) to dry lands (). Individuals growing in aquatic and terrestrial habitats exhibited notable morphological differences, with the plants growing in aquatic environments showing significantly longer internode length and greater stem diameter (; ). Plasticity in internode elongation, aerenchyma formation, and other morphological and physiological traits allows the amphibious plants of A. philoxeroides to withstand alternative wetting and drying conditions and to proliferate in both aquatic and terrestrial environments (, ; ).
As the most conspicuous responses of plants subject to submergence, enhanced internode/petiole elongation and secondary aerenchyma formation have been found in many species (; ; ; ). The molecular and cellular mechanisms underlying these responses, such as the external and internal environmental cues, signaling pathways, and genetic elements implicated in developmental outcomes, have been explored in deepwater rice (), Rumex (), and other species (; ; ). However, the epigenetic mechanisms underpinning submergence-induced plastic development, including the potential roles of microRNA-mediated regulation of gene expression, remain elusive. Since the transition from aquatic to terrestrial environments, or vice versa, requires timely and flexible adjustments in morphological traits, physiological processes, and developmental machinery to extremely different conditions (; ), the epigenetic control system is most likely to provide an effective short-term strategy for amphibious plants to plastically respond to environmental fluctuations by reversibly controlling effector gene expression states depending on environmental conditions.
The genetically identical but phenotypically distinct plants of A. philoxeroides provide a suitable model for us to rule out genetic contributions to observed differences, to explore the correlations between environmental stimuli, epigenetic modification and phenotypic variation. Previously, we have investigated the gene expression reaction norms of A. philoxeroides in contrasting hydrological conditions, and identified 11 transcriptionally coordinated gene groups associated with environment-induced phenotypic variation in A. philoxeroides (). We have also examined genome-wide DNA methylation patterns of A. philoxeroides in natural and manipulated habitats, and found that plants of different source populations not only underwent significant morphological changes in common garden environments, but also underwent genome-wide DNA methylation alterations in response to different water treatments, demonstrating a correlation between epigenetic reprogramming and the reversible phenotypic response of A. philoxeroides to particular environmental factors (; ). In this study, we aim to detect microRNAs potentially involved in the regulation of the plastic responses of A. philoxeroides to contrasting hydrological conditions, especially microRNAs associated with enhanced internode elongation growth under submergence. By preliminarily profiling the time-dependent changes of A. philoxeroides miRNome in alternate wetting and drying conditions using deep sequencing technology, and validating the most differentially expressed microRNAs by independent quantitative reverse transcription-PCR (qRT-PCR), we identified the submergence-sensitive microRNAs of A. philoxeroides, and determined microRNA:mRNA functional pairs and their temporal expression dynamics possibly associated with submergence-induced internode elongation in A. philoxeroides. This expanded investigation enhanced our understanding of the role of microRNAs in mediating adaptive developmental response to environmental changes, promoting the study of epigenetic regulation and phenotypic variation in an ecological context.
Materials and Methods
Plant Materials, Common Garden Treatments, and Phenotypic Evaluation
Plants of A. philoxeroides used in this study were collected from the wild populations in Zhuji, Zhejiang Province (E120°29′ N29°40′) in 2006, and were maintained in the green house of Fudan University, Shanghai (E121°29′ N31°14′). Healthy storage roots from genetically uniform plants were used to produce asexual seedlings. The storage roots were sowed in plastic plates. After the first two new leaves appeared, seedlings of similar sizes were individually transplanted into pots (20 cm in diameter, 15 cm in depth) containing a 1:1:1 mix of peat, vermiculite and sand. After 2 months, the seedlings were subjected to two different types of treatments: pond (submergence) and upland treatments. Two treatments were conducted simultaneously, during the summer months of July and August with about 12 h of daytime, and under identical climatic conditions. In the pond treatment, the potted plants were completely submerged in a series of plastic tanks containing tap water to a depth of 50 cm. Plants were supplied with 1 L water per day in the upland treatment and the soil was kept well-drained (). The length and diameter of the fifth internode from top of 10 individuals were measured under each treatment, at 0, 1, 3, 6, 12, and 120 h from the start of treatments. Tukey’s HSD (honest significant difference) test was performed based on raw data to test differences among sample means for significance. Internode tissues were also harvested and stored in RNALater (Ambion) for RNA extraction. Five replicates from separate plants were collected for each time point and treatment.
RNA Isolation and Small RNA Sequencing
Total RNAs of each sample were isolated using TRIzol reagent (Invitrogen, Carlsbad, CA, United States), following manufacturer’s instructions. The quality and quantity of obtained RNAs were evaluated on a Bioanalyzer 2100 (Agilent Technologies). RNAs of five replicates for each sample were equivalently mixed according to different time points and treatments. The small RNA libraries were constructed following the standard protocol. In brief, small RNAs with 18–30 nt were separated from the total RNA by polyacrylamide gel electrophoresis. The selected fragments were then ligated to a 5′ RNA adapter and a 3′ RNA adapter, followed by reverse transcription and PCR amplification. PCR products were collected by gel purification and sequenced using the Illumina HiSeq 2000 platform at the Beijing Genomics Institute (BGI), Shenzhen, China. A total of eleven libraries were constructed and sequenced, including an untreated sample library (the control library G0), five upland libraries (D1, D3, D6, D12, and D120) and five pond treatment libraries (W1, W3, W6, W12, and W120). The sequence data were deposited in the United States National Center for Biotechnology Information (NCBI) Sequence Read Archive under accession number SRR5891543 to SRR5891553.
Identification of Known and Novel microRNAs
Raw reads produced by sequencing were first gone through a data cleaning pipeline to remove low quality reads, contaminating sequences, and reads less than 18 nt. The clean reads were further filtered to remove reads with less than five counts in each library (). Sequences matching non-coding rRNAs, tRNAs, snRNAs, and snoRNAs in the NCBI GenBank1 and Rfam2 databases were eliminated. The remaining small RNA sequences were aligned against precursor/mature microRNA sequences in miRBase 213 to identify known (conserved) microRNAs in A. philoxeroides allowing two mismatches and three gaps. The highest expression microRNA of each mature microRNA family was selected to create a temporary microRNA database. Then clean reads was aligned to this temporary microRNA database and the frequency of each microRNA was generated by summing the count of tags which can align to the temporary microRNA database within two mismatches.
Clean reads that could not be assigned to the aforementioned RNA classes were used to predict potential novel microRNAs in A. philoxeroides. Due to the non-availability of the genome sequence of A. philoxeroides, the transcriptome of A. philoxeroides () was used to predict novel microRNAs. MIREAP4 was used to examine the hairpin structure, the Dicer cleavage site and the minimum free energy of candidate sequences according to the following parameters: (i) minimal microRNA reference sequence length of 20 nt and maximal microRNA reference sequence length of 24 nt; (ii) maximal copy number of microRNAs on the reference was 20; (iii) maximal free energy allowed for a microRNA precursor was -18 kcal/mol; (iv) maximal space between the microRNA and microRNA∗ was 300 nt; (v) minimal base pairs of microRNA and microRNA∗ was 16; (vi) maximal bulge of microRNA and microRNA∗ was 4; (vii) maximal asymmetry of the microRNA/microRNA∗ duplex was 4; and (viii) flank sequence length of microRNA precursor was 20.
Differential Expression Analysis of microRNAs
The frequency of known and novel microRNAs in each library was then normalized to reads per million (RPM) by dividing the read counts of a microRNA by the total number of clean reads of each library, and multiplying this number by 105. The temporal expression pattern of each microRNAs under different treatments was detected by calculating the fold-changes of five upland-treated (D1, D3, D6, D12, and D120) vs. control (G0) and five pond-treated (W1, W3, W6, W12, and W120) samples vs. control (G0), respectively. We considered a fold change of at least 2 (| Log2 (treatment/control)| >= 1, P-value < 0.05) as an indication of significant change. If the read count of a microRNA was 0 in one sample, its expression value was revised to 0.01; If the RPM of a microRNA was less than 1 in both samples compared, this microRNA was excluded from differentially expression analysis, owing to its low expression (). Pearson’s r was used to measure the correlations of microRNA expression in upland- and pond-treated samples, using the R software package for statistical programming, version 3.1.16. To identify differentially expressed microRNAs induced by distinct treatments, the fold-changes of microRNA at different time points of treatment were also calculated between contrast upland- and pond- treated samples. MicroRNAs with the log2 fold-change above 1 and a p-value < 0.05 (except microRNAs expressed only in upland- treated samples) were defined as submergence-responsive microRNAs.
Prediction and Functional Annotation of microRNA Targets
To identify the potential targets, all known and novel microRNAs identified in this study were used as queries to search against the A. philoxeroides transcriptome () using the psRNATarget program with default parameter setting (). The BGI microRNA target prediction pipeline which based on the criteria suggested by the two previous reports (; ) was also used for microRNA target prediction. A pool of potential target genes of identified microRNAs was constructed by combining the results of two target prediction tools. By filtering with gene expression profiles (), the target genes that were differentially expressed in contrasting hydrological conditions were identified. GO enrichment analysis was performed using the Bioconductor package topGO6. The results of GO enrichment analysis were summarized using WEGO (). The target genes of known and novel submergence-responsive microRNAs were also subjected to GO annotation separately. The web tool REVIGO was used to summarize and visualize the results, which could cluster long, unintelligible lists of GO terms into several representative subsets (). Additionally, not all microRNAs were expressed in all eleven libraries. A part of microRNAs were only expressed in one library. To explore the functional difference of target genes between all-library-expressed microRNAs and one-library-expressed microRNAs, multi-GOEAST analysis was used to compare target genes GO enrichment status between these two categories of microRNAs ().
qRT-PCR Validation of Transient Co-expression of microRNAs with Their Target Genes
Real-time qPCR was used to validate the co-expression patterns of submergence-responsive microRNAs and their target genes. Four conserved and four novel submergence-responsive microRNAs, were selected and subjected to qRT-PCR analysis. Plant materials used for qRT-PCR validation were collected from another independent common garden experiment. Total RNAs of different samples were extracted using the method mentioned above. For quantification of microRNA expression, the total RNA (1 μg) was first polyadenylated and then reverse transcribed into cDNA with an oligo-dT adapter primer, using the miRcute microRNA First-Strand cDNA Synthesis Kit (TIANGEN, Beijing, China). qRT-PCR was carried out using the miRcute microRNA qPCR Detection kit (SYBR Green) (TIANGEN, Beijing, China) following manufacturer’s instructions. The specific forward primer designed for each microRNA and the universal reverse primer provided by kit were used for amplification. For quantification of target mRNAs by qRT-PCR, the total RNA of each sample was reverse transcribed into cDNA using PrimeScriptTM RT Master Mix (Perfect Real Time) (Takara, Dalian, China). qRT-PCR was performed with LightCycler® 480 SYBR Green I Master (Roche, United States) following manufacturer’s instructions. Ubiquitin C (UBC) was used as an internal control. The primers used for qRT-PCR were listed in Supplementary Table S1. All reactions were performed in triplicate for each sample on LightCycler® 96 System. The relative time-course changes in microRNA and mRNA expression levels between pond and upland treatments were quantified using the 2-ΔΔCt method ().
Target Validation by 5′ RLM-RACE
The 5′ RLM-RACE (RNA ligase-mediated 5′ rapid amplification of cDNA ends) method was used for target validation (). The total RNAs of pond-treated samples were isolated using TRIzol reagent (Invitrogen, Carlsbad, CA, United States) following manufacturer’s instructions. mRNAs were separated from the total RNA using PolyA Tract mRNA Isolation System III (Promega, United States), and were then used to perform 5′ RLM-RACE using GeneRacerTM Kit (Invitrogen, Carlsbad, CA, United States), according to the instructions of the manufacturer. Gene-specific nested reverse primers (listed in Supplementary Table S1) were used in combination with the 5′ primer provided by the kit to amplify the cleaved transcripts. Amplicons were purified using the AxyPrep DNA gel extraction kit (Axygen, United States), cloned into the pEASY-T3 vector (Transgene, Beijing, China), and sequenced.
Results
Phenotypical Variation of A. philoxeroides under Pond and Upland Treatments
Internode elongation and increases in stem diameter were not significant between plants under distinct treatments at the first four time points (1, 3, 6, and 12 h). One hundred and twenty hours later, significant differences were found between plants subjected to pond and upland treatments, respectively, with the plants growing in pond environment exhibiting significantly enhanced internode elongation, accompanied by increases in stem diameter (Figure 1).
FIGURE 1
Conserved and Novel microRNAs of A. philoxeroides
Illumina HiSeq was used to sequence the small RNA libraries. On average, 12.98 million (11719363 to 14236593) reads were obtained for each library. After stringent quality check and data cleaning, approximately 12.73 million (11596962 to 14093418) clean reads were left for each library (Supplementary Table S2). Clean reads were classified as different types of non-coding transcripts, such as rRNA, tRNA, snRNA, snoRNA, microRNA, and unannotated reads. The numbers of each class in each library were shown in Supplementary Table S3. The 24 nt class was the most abundant group of sRNAs in all libraries followed by 21 nt sequences (Supplementary Figure S1). Conserved microRNAs were identified by BLAST searching against the miRBase database. A total of 220 known microRNAs belonging to 76 families were identified. The number of members varied in different microRNA families. The MIR166 family had the most abundant members in A. philoxeroides (nine members), followed by MIR159 (six members). Four families (MIR167_1, MIR169_1, MIR172, and MIR396) contained five members. Fifty-two families (68.4%) had only one representative in A. philoxeroides. Investigation of the distribution of the 220 known microRNAs among different libraries revealed that 113 (51.4%) microRNAs were expressed not only in the pond- and upland-treated samples, but also in the control library G0. Among them, 66 microRNAs were constantly expressed in all libraries across different time points. These microRNAs came from 32 families, and about a half of these families were well-conserved across the plant kingdom, such as MIR159 (three members) and MIR166 (six members). Some known microRNAs were found to be specific to the control (3), the pond-treated (37) and the upland-treated (18) samples, respectively (Supplementary Figure S2A). The microRNAs specifically expressed in the pond-treated samples mostly belonged to the less-conserved microRNA families such as MIR824 (one member), MIR837 (one member) and MIR7982 (1 member), yet some of them came from well-conserved families such as MIR169_1 (miR169d-5p, miR169n-3p, and miR169i).
The MIREAP prediction software was used to identify potentially novel microRNAs in A. philoxeroides. A total of 81 novel microRNAs were identified. The length of the novel microRNA precursors varied from 72 to 345 nt, with an average of 131 nt. The average minimum free energy (MFE) ranged from -107.9 to -19 kcal/mol, with an average of -37.93 kcal/mol. The sequences and hairpin secondary structures of novel microRNAs were shown in Supplementary Table S4. Among 81 newly identified microRNAs, 12 (14.8%) were expressed in both the control and the pond- and upland-treated samples. 16 (19.75%) were shared by the pond- and upland-treated samples. 5 (6.17%), 21 (25.92%), and 24 (29.62%) were specific to the control, pond- and upland-treated samples, respectively (Supplementary Figure S2B).
Expression Dynamics of Identified microRNAs
Time course expression profiles revealed that, among 301 known and novel microRNAs identified in A. philoxeroides, 70 (66 known and 4 novel) microRNAs were constantly expressed in all libraries across different time points of the treatments. The rest of microRNAs showed varied expression between different libraries. Ninety-four (47 known and 47 novel) microRNAs were expressed only in a single library. Comparison of the expression profiles led to the identification of 179 (153 known and 26 novel) microRNAs that were differentially expressed between the pond- and upland-treated samples. Of them, 24 and 40 microRNAs were specifically expressed in the upland- (Figure 2A and Supplementary Table S5) and pond-treated (Figure 2B and Supplementary Table S5) samples, respectively. The rest 115 microRNAs were expressed in both pond- and upland-treated samples but with different expression patterns (Figure 2C and Supplementary Table S5). Twelve microRNAs showed consistent expression patterns between the upland- and pond-treated samples (r > 0.9, p < 0.05, Figure 2D and Supplementary Table S5). To identify submergence-responsive microRNAs, the upland-treated samples were used as controls to calculate the relative expression value of microRNAs between treatments. A total of 146 (128 known and 18 novel) microRNAs were identified as submergence-responsive microRNAs (Supplementary Table S6), and their expression levels were significantly altered at different time points of the treatments (Figure 2E). Eighty-four microRNAs (e.g., miR1023a-3p, miR3630-3p, and novel_mir_25) responded to the pond treatment at the early stage of the treatment (1–12 h), and 19 microRNAs (e.g., miR1439, miR7542 and novel_mir_48) responded at the later stage (120 h). Forty-three microRNAs (e.g., miR167b-3p, miR5488, and novel_mir_73) exhibited varied patterns of differential expression between the upland- and pond-treated samples at different stages of the treatments.
FIGURE 2
Predicted Targets of Known and Novel microRNAs
Based on in silico target prediction, a total of 2503 genes were identified as the potential target genes of the known and novel A. philoxeroides microRNAs, in which 1439 were found to be differentially expressed in pond and upland environments. Among these differentially expressed genes, 732 were predicted as the targets of known submergence-responsive microRNAs, and 132 were targeted by novel submergence-responsive microRNAs. The cleavage sites of two target genes (Contig16589_AHP4 and Contig18934_UDP-Glycosyltransferase superfamily protein) recognized by novel_mir_19 and novel_mir_25, respectively, were confirmed by 5′ RLM-RACE experiments. The cleavage site was concentrated at the 10th nucleotide from the 5′ end of each microRNA (Figure 3).
FIGURE 3
The results of the GO term enrichment analysis showed that the microRNAs ubiquitously expressed in all samples were more likely to regulate genes involved in basic life processes, whereas the targets of transiently expressed microRNAs were much more versatile than those of constantly expressed microRNAs (Supplementary Figure S3). The target genes regulated by submergence-responsive microRNAs were involved in a variety of cellular and physiological processes, such as anatomical structure formation, response to stimulus and developmental or metabolic process, and the distributions of the different GO categories are shown in Supplementary Figure S4. Interestingly, genes targeted by known microRNAs seem more likely to be engaged in cell wall metabolism and remodeling in response to flooding stress, including genes participating in cell wall organization, lignin catabolism, and UDP-D-xylose biosynthesis (Supplementary Figure S5A). Instead, targets of novel microRNAs mainly participated in plant hormone response and epigenetic processes, such as cellular response to gibberellin and regulation of DNA methylation (Supplementary Figure S5B).
Co-expression Patterns between microRNAs and Their Target mRNAs
Contrasting hydrologic treatments led to global changes in microRNA and mRNA expression patterns in A. philoxeroides. By integrating expression profiles of microRNA and mRNA in conjunction with the predicted microRNA targets, we identified eight microRNA:mRNA pairs potentially involved in the regulation of internode elongation under submergence. Quantitative analysis using qRT-PCR confirmed the inverse expression patterns between microRNAs and their target mRNAs, with the exception of novel_mir_29:Contig23974 (Figure 4). Mir167d-5p was dramatically induced after 12 h of submergence following a subtle up-regulation at 3 h of the pond treatment, and then down-regulated. Its target mRNA cotig45839 showed the expected negative regulation pattern (Figure 4A). mir3630-3p, novel_mir_25 and novel_mir_26 were transiently up-regulated at 6 h of submergence and returned to normal levels at 12 h of the pond treatment, with the expression levels of their target mRNAs (contig7264, contig18934/contig49354, and contig539) being negatively correlated with those of the corresponding microRNAs (Figures 4D,F,G). miR169u, miR396a-3p and novel_mir_19 exhibited fluctuate expression patterns. They were up- and down-regulated at the early stage (1–12 h) of the treatment, and then up again during the late phase (120 h) (Figures 4B,C,E). The co-expression pattern of novel_mir_29 with its target gene Contig23974 was distinct from others. Although the expression of novel_mir_29 was differentially regulated at different time points of the treatment, the expression level of the putative target transcript was not changed accordingly in an inversely correlated pattern (Figure 4H).
FIGURE 4
Discussion
Due to their sessile lifestyle, plants have developed a remarkable plasticity for shaping their body plan to adapt to changes in their immediate environment. The high sensitivity of microRNAs to environmental perturbations and their significant role in readjusting the expression levels of genes and phenotypic profiles under changing environments have made them the best candidates as dynamic regulators of phenotypic plasticity ().
Genome-wide perspective of miRNAome in response to pond and upland treatments revealed 146 submergence-responsive microRNAs belonging to 49 microRNA families in A. philoxeroides. The members from 19 families (such as MIR159, MIR167_1, MIR169, MIR396, MIR482, MIR530, and MIR3630) were not only involved in responses to water-logging/flooding in A. philoxeroides, but also in other plant species. For instance, gene expression profiling and qRT-PCR analyses revealed the varied expression patterns of miR169u and miR396a-3p in A. philoxeroides during submergence. The similar patterns were also found in Nelumbo nucifera () and Zea mays (; ). Three members of the miR167_1 family (miR167b-3p, miR167d-5p, and miR167d-3p) showed significantly altered expression between the pond- and upland-treated samples in A. philoxeroides. Similarly, the members from this family were differentially expressed under flooding condition in P. tomentosa (). These results suggest the common roles of these microRNAs in different plant species. The submergence-responsive microRNAs identified in A. philoxeroides can also be induced by components of submergence, such as hypoxic stress (e.g., miR159, miR169, miR319, miR394, and miR395) (), elevated CO2 concentration (e.g., miR156, miR157, miR172, and miR160) (), darkness (e.g., miR160, miR167, miR171, miR396, and miR398) () and ethylene entrapment (e.g., miR399, miR408, miR482, miR4414, and miR5139) ().
Target prediction and functional annotation of the submergence-responsive microRNAs recognized microRNAs potentially associated with plastic internode elongation in A. philoxeroides upon submergence. The target of novel_mir_26 was homeobox protein 22 (HB22), which has been shown to be associated with the synthesis of gibberellic acid (GA) in Arabidopsis (). The primary action of GA in stem growth was on cell elongation (). The expression of A. philoxeroides HB22 (Contig539) was significantly up-regulated at 12 h of the pond treatment accompanied by down-regulation of novel_mir_26, that might promote internode elongation in submerged plants by enhancing GA biosynthesis (). Both cytokinin and brassinosteroid (BR) were also implicated in stem elongation growth (; ). The target of novel_mir_19 in A. philoxeroides was contig16589, which was homologous to Arabidopsis AHP4 and potentially involved in cytokinin signaling pathway (). Novel_mir_19 and its target exhibited a significant inverse expression pattern under the pond treatment (Figure 4E), and the target cleavage site was verified by 5′ RLM-RACE. A. philoxeroides contig2589, a homolog of Arabidopsis CDL1 that positively regulates BR signaling () and stem elongation (), was predicted as the target of miR398a-5p. The expression levels of contig2589 and miR398a-5p were also changed accordingly in an inversely correlated pattern according to the expression profile.
Abscisic acid (ABA) antagonizes the effects of GA and other phytohormones in numerous processes during plant growth and development (; ). Growth-regulating factors (GRFs), a conserved class of plant-specific transcription factors, participate in ABA signaling pathway as a transcriptional repressor of ABA- and osmotic-stress responsive genes (). The expression of several GRFs is controlled by miR396 (), functioning as a regulator of hormone homeostasis and signaling under flooding stress (). In our study, miR396-3p presented fluctuate expression patterns opposite to contig34390 (GRF6) during the pond treatment (Figure 4C), indicating the interplay between ABA and other phytohormones in A. philoxeroides upon submergence. Previous studies have shown that interactions between different plant hormones play a central role in generating plastic phenotypic responses to environmental variation (; ), including in promoting shoot elongation upon submergence (; ; , ; ; ). The detection of dynamic expressions of microRNAs and their corresponding targets associated with different hormone signaling pathways in this study possibly suggests the potential role of A. philoxeroides microRNAs in coordinating multiple phytohormone responses during the submergence treatment. The crosstalk between microRNAs and hormone signaling cascades deserves better attention while studying submergence-induced elongation growth in flooding-tolerant plants.
Previous gene expression profiling showed that a large number of genes associated with different processes of cell wall strengthening, such as the xylan biosynthetic process, the lignin biosynthetic process, secondary cell wall biogenesis, etc., exhibited a significant “up-down” pattern in the early stage of the pond treatment, suggesting a “Just-in-time” cell wall thickening in response to submergence though its function remaining unclear (). In coincidence with the transient cell wall strengthening, we found in this study that NF-YA3 (A. philoxeroides contig18972), a transcript factor which has been shown to impair cell elongation (), was transiently up-regulated in the early stage of the pond treatment (Figure 4B), inversely correlated with the expression of the corresponding miR169u. In contrast, the expression of SKS5 (A. philoxeroides contig45839), a gene that has been identified as a cell wall-related gene () and was expressed most strongly in expanding tissues (), was highly repressed in the early stage of the pond treatment accompanied by induced and transient up-regulation of the corresponding miR167d-5p (Figure 4A). The transient down-regulation of the cell wall expansion-promoting gene SKS5 was consistent with the transient up-regulation of cell wall strengthening genes. microRNAs, while functioning as a trigger of hormone response to external environment, may also coordinate the expression of genes associated with different processes of cell wall modification under submergence.
However, the relationships between microRNAs and their target genes was not always straightforward. One microRNA seemed to be able to regulate more than one target gene and, in turn, one gene might be regulated by more than one microRNA. For example, both contig49354 (DUF579) and contig18934 (UDP-Glycosyltransferase superfamily protein) were identified as the targets of novel_mir_25 and showed similar inverse relationships with novel_mir_25 in expression (Figure 4F), while the expression of contig2773 (SPL3) was regulated by miR157a and miR4387b, respectively, at different time points of the treatment. We also found a positive relation between the expression of novel_mir_29 and its target Contig23974 (vacuolar invertase, VI2), a gene that has long been considered as a major player in cell expansion () (Figure 4H). The most striking feature revealed by time course expression profiling was the fluctuating pattern in microRNA expression. Most of the differentially expressed microRNAs were transiently up-/down-regulated at certain time points during the treatments. More than half of the submergence-responsive microRNAs (e.g., miR1023a-3p, miR3630-3p, and novel_mir_25) responded to the pond treatment at the early stage of the treatment (1–12 h). About one-third of the microRNAs were differentially expressed at the 1st hour from the start of treatments, while 19 microRNAs (e.g., miR1439, miR7542, and novel_mir_48) responded to the treatment 120 h later. The fine-scale temporal changes in microRNA expression highlight the importance of time-series sampling in identifying stress-responsive microRNAs and analyzing their role in stress response/tolerance. Incomplete sampling might give a biased conclusion about the microRNA’s role in readjusting the expression levels of genes and phenotypic profiles under changing environments.
Statements
Author contributions
GL performed the wet lab work and data analysis, and drafted the manuscript. YD, YG, and CZ helped to set up experiments. YW, WZ, and ZS participated in data analysis. LG and JY conceived the idea, participated in the design of the study and finalized the manuscript. All of the authors read and approved the final manuscript.
Funding
This study was supported by National Key R&D Program of China (2016YFC1200603), the National Natural Science Foundation of China (31670222 and 31300189), the Science and Technology Commission of Shanghai Municipality (14DZ2260400) and Chenshan Key Scientific Research Projects (G142431).
Acknowledgments
The authors are grateful to Dr. Luoyan Zhang (Shandong Normal University) and Dr. Jian Lu (Peking University) for their helpful advice on data analysis.
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. The reviewer OV-L and handling Editor declared their shared affiliation.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2017.02078/full#supplementary-material
FIGURE S1Length distribution of small RNAs in each library.
FIGURE S2microRNAs expressed in different sample libraries. (A) Known microRNAs; (B) novel microRNAs.
FIGURE S3Results of multi-GOEAST analysis, showing that targets of ubiquitously expressed microRNAs were more likely to be involved in basic life processes (red), whereas the targets of transiently expressed microRNAs were much more versatile (green). Yellow boxes represent common GO terms enriched in both categories of targets. The degree of color saturation indicates the significance of enrichment for corresponding GO terms. Arrows represent connections between different GO terms. Red arrows represent direct relationships between two GO terms, and black arrows (solid and dashed) represent indirect relationships between GO terms.
FIGURE S4Gene Ontology classifications of the target genes based on cellular component, molecular function, and biological process.
FIGURE S5Summarized biological process of known (A) and novel (B) microRNA target genes. Loosely related terms were jointed into “super-clusters” by different colors with a representative showed in middle. Size of the rectangles was adjusted to p-value of GO annotation.
Footnotes
1.^https://www.ncbi.nlm.nih.gov/genbank/
2.^http://www.sanger.ac.uk/Software/Rfam
4.^https://sourceforge.net/projects/mireap/
5.^http://bioconductor.org/packages/release/bioc/html/topGO.html
References
1
AchardP.HerrA.BaulcombeD. C.HarberdN. P. (2004). Modulation of floral development by a gibberellin-regulated microRNA.Development1313357–3365. 10.1242/dev.01206
2
AllenE.XieZ.GustafsonA. M.CarringtonJ. C. (2005). microRNA-Directed phasing during trans-acting siRNA biogenesis in plants.Cell121207–221. 10.1016/j.cell.2005.04.004
3
ArneyS.MancinelliP. (1966). The basic action of gibberellic acid in elongation of ‘Meteor’ pea stems.New Phytol.65161–175. 10.1111/j.1469-8137.1966.tb06349.x
4
Aubin-HorthN.RennS. C. (2009). Genomic reaction norms: using integrative biology to understand molecular mechanisms of phenotypic plasticity.Mol. Eol.183763–3780. 10.1111/j.1365-294X.2009.04313.x
5
AzpirozR.WuY.LoCascioJ. C.FeldmannK. A. (1998). An Arabidopsis brassinosteroid-dependent mutant is blocked in cell elongation.Plant Cell10219–230. 10.1105/tpc.10.2.219
6
Bailey-SerresJ.VoesenekL. A. (2010). Life in the balance: a signaling network controlling survival of flooding.Curr. Opin. Plant Biol.13489–494. 10.1016/j.pbi.2010.08.002
7
BianH.XieY.GuoF.HanN.MaS.ZengZ.et al (2012). Distinctive expression patterns and roles of the miRNA393/TIR1 homolog module in regulating flag leaf inclination and primary and crown root growth in rice (Oryza sativa).New Phytol.196149–161. 10.1111/j.1469-8137.2012.04248.x
8
BilichakA.IlnystkyyY.HollunderJ.KovalchukI. (2012). The progeny of Arabidopsis thaliana plants exposed to salt exhibit changes in DNA methylation, histone modifications and gene expression.PLOS ONE7:e30515. 10.1371/journal.pone.0030515
9
BlazquezM. A.GreenR.NilssonO.SussmanM. R.WeigelD. (1998). Gibberellins promote flowering of Arabidopsis by activating the LEAFY promoter.Plant Cell10791–800. 10.1105/tpc.10.5.791
10
BoegerM. R. T.PoulsonM. E. (2003). Morphological adaptations and photosynthetic rates of amphibious Veronica anagallis-aquatica L. (Scrophulariaceae) under different flow regimes.Aquat. Bot.75123–135. 10.1016/S0304-3770(02)00174-2
11
BorgesR. M. (2008). Plasticity comparisons between plants and animals: concepts and mechanisms.Plant Signal. Behav.3367–375. 10.4161/psb.3.6.5823
12
BossdorfO.RichardsC. L.PigliucciM. (2008). Epigenetics for ecologists.Ecol. Lett.11106–115. 10.1111/j.1461-0248.2007.01130.x
13
BradshawA. D. (1965). Evolutionary significance of phenotypic plasticity in plants.Adv. Genet.13115–155. 10.1016/S0065-2660(08)60048-6
14
BuesoE.Munoz-BertomeuJ.CamposF.BrunaudV.MartinezL.SayasE.et al (2014). ARABIDOPSIS THALIANA HOMEOBOX25 uncovers a role for Gibberellins in seed longevity.Plant Physiol.164999–1010. 10.1104/pp.113.232223
15
ChinnusamyV.ZhuJ. K. (2009). Epigenetic regulation of stress responses in plants.Curr. Opin. Plant Biol.12133–139. 10.1016/j.pbi.2008.12.006
16
CoxM. C. H.BenschopJ. J.VreeburgR. A.WagemakerC. A.MoritzT.PeetersA. J.et al (2004). The roles of ethylene, auxin, abscisic acid, and gibberellin in the hyponastic growth of submerged Rumex palustris petioles.Plant Physiol.1362948–2960. 10.1104/pp.104.049197
17
CoxM. C. H.PeetersA. J. M.VoesenekL. A. C. J. (2006). The stimulating effects of ethylene and auxin on petiole elongation and on hyponastic curvature are independent processes in submerged Rumex palustris.Plant Cell Environ.29282–290. 10.1111/j.1365-3040.2005.01420.x
18
DaiX.ZhaoP. X. (2011). psRNATarget: a plant small RNA target analysis server.Nucleic Acids Res.39W155–W159. 10.1093/nar/gkr319
19
DengY.GaoL.ZhuZ.YangJ. (2014). Differential expression of DNA methylation regulating factors and dynamic methylation patterns of Alternanthera philoxeroides under different water treatments.Plant Sci. J.32475–486. 10.11913/PSJ.2095-0837.2014.50475
20
EldemV.OkayS.ÜnverT. (2013). Plant microRNAs: new players in functional genomics.Turk. J. Agric. For.371–21. 10.3906/sag-1210-51
21
FinkelsteinR. (2013). Abscisic acid synthesis and response.Arabidopsis Book11:e0166. 10.1199/tab.0166
22
GaoL.GengY.LiB.ChenJ.YangJ. (2010). Genome-wide DNA methylation alterations of Alternanthera philoxeroides in natural and manipulated habitats: implications for epigenetic regulation of rapid responses to environmental fluctuation and phenotypic variation.Plant Cell Environ.331820–1827. 10.1111/j.1365-3040.2010.02186.x
23
GaoL.GengY.YangH.HuY.YangJ. (2015). Gene expression reaction norms unravel the molecular and cellular processes underpinning the plastic phenotypes of Alternanthera philoxeroides in contrasting hydrological conditions.Front. Plant Sci.6:991. 10.3389/fpls.2015.00991
24
GaoZ.ShiT.LuoX.ZhangZ.ZhuangW.WangL. (2012). High-throughput sequencing of small RNAs and analysis of differentially expressed microRNAs associated with pistil development in Japanese apricot.BMC Genomics13:371. 10.1186/1471-2164-13-371
25
GengY. P.PanX. Y.XuC. Y.ZhangW. J.LiB.ChenJ. K. (2006). Phenotypic plasticity of invasive Alternanthera philoxeroides in relation to different water availability, compared to its native congener.Acta Oecol.30380–385. 10.1016/j.actao.2006.07.002
26
GengY. P.PanX. Y.XuC. Y.ZhangW. J.LiB.ChenJ. K.et al (2007). Phenotypic plasticity rather than locally adapted ecotypes allows the invasive alligator weed to colonize a wide range of habitats.Biol. Invasions9245–256. 10.1007/s10530-006-9029-1
27
GhalamborC. K.McKayJ. K.CarrollS. P.ReznickD. N. (2007). Adaptive versus non-adaptive phenotypic plasticity and the potential for contemporary adaptation in new environments.Funct. Ecol.21394–407. 10.1111/j.1365-2435.2007.01283.x
28
GilbertS. F. (2016). Chapter twenty-two-developmental plasticity and developmental symbiosis: the return of eco-devo.Cur. Top. Dev. Biol.116415–433. 10.1016/bs.ctdb.2015.12.006
29
HansenM.ChaeH. S.KieberJ. J. (2009). Regulation of ACS protein stability by cytokinin and brassinosteroid.Plant J.57606–614. 10.1111/j.1365-313X.2008.03711.x
30
HuoX.WangC.TengY.LiuX. (2015). Identification of miRNAs associated with dark-induced senescence in Arabidopsis.BMC Plant Biol.15:266. 10.1186/s12870-015-0656-5
31
HutchisonC. E.LiJ.ArguesoC.GonzalezM.LeeE.LewisM. W.et al (2006). The Arabidopsis histidine phosphotransfer proteins are redundant positive regulators of cytokinin signaling.Plant Cell183073–3087. 10.1105/tpc.106.045674
32
JacksonM. B. (2008). Ethylene-promoted elongation: an adaptation to submergence stress.Ann. Bot.101229–248. 10.1093/aob/mcm237
33
JacksonM.ArmstrongW. (1999). Formation of aerenchyma and the processes of plant ventilation in relation to soil flooding and submergence.Plant Biol.1274–287. 10.1055/s-2007-978516
34
JinD.WangY.ZhaoY.ChenM. (2013). MicroRNAs and their cross-talks in plant development.J. Genet. Genomics40161–170. 10.1016/j.jgg.2013.02.003
35
JinQ.XuY.MattsonN.LiX.WangB.ZhangX.et al (2017). Identification of submergence-responsive microRNAs and their targets reveals complex miRNA-mediated regulatory networks in lotus (Nelumbo nucifera Gaertn).Front. Plant Sci.8:6. 10.3389/fpls.2017.00006
36
Jones-RhoadesM. W.BartelD. P.BartelB. (2006). MicroRNAs and their regulatory roles in plants.Annu. Rev. Plant Biol.5719–53. 10.1146/annurev.arplant.57.032905.105218
37
JulienM. H.SkarrattB.MaywaldG. (1995). Potential geographical distribution of alligator weed and its biological control by Agasicles hygrophila.J. Aquat. Plant Manage.3355–60.
38
JulienM. H.StanleyJ. N. (1999). “The management of alligator weed, a challenge for the new millennium,” inProceedings of the 10th Biennial Noxious Weed Conference (Ballina, NSW: New South Wales Department of Agriculture).
39
KaliszS.KramerE. (2008). Variation and constraint in plant evolution and development.Heredity100171–177. 10.1038/sj.hdy.6800939
40
KhanM. S.KhraiweshB.PugalenthiG.GuptaR. S.SinghJ.DuttamajumderS. K.et al (2014). Subtractive hybridization-mediated analysis of genes and in silico prediction of associated microRNAs under waterlogged conditions in sugarcane (Saccharum spp.).FEBS Open Bio4533–541. 10.1016/j.fob.2014.05.007
41
KimJ. S.MizoiJ.KidokoroS.MaruyamaK.NakajimaJ.NakashimaK.et al (2012). Arabidopsis growth-regulating factor7 functions as a transcriptional repressor of abscisic acid- and osmotic stress-responsive genes, including DREB2A.Plant Cell243393–3405. 10.1105/tpc.112.100933
42
KimT. W.GuanS.BurlingameA. L.WangZ. Y. (2011). The CDG1 kinase mediates brassinosteroid signal transduction from BRI1 receptor kinase to BSU1 phosphatase and GSK3-like kinase BIN2.Mol. Cell.43561–571. 10.1016/j.molcel.2011.05.037
43
LemaS. C. (2014). Hormones and phenotypic plasticity in an ecological context: linking physiological mechanisms to evolutionary processes.Integr. Comp. Biol.54850–863. 10.1093/icb/icu019
44
Leyva-GonzálezM. A.Ibarra-LacletteE.Cruz-RamírezA.Herrera-EstrellaL. (2012). Functional and transcriptome analysis reveals an acclimatization strategy for abiotic stress tolerance mediated by Arabidopsis NF-YA family members.PLOS ONE7:e48138. 10.1371/journal.pone.0048138
45
LiC.ZhangB. (2016). MicroRNAs in control of plant development.J. Cell. Physiol.231303–313. 10.1002/jcp.25125
46
LicausiF.WeitsD. A.PantB. D.ScheibleW. R.GeigenbergerP.van DongenJ. T. (2011). Hypoxia responsive gene expression is mediated by various subsets of transcription factors and miRNAs that are determined by the actual oxygen availability.New Phytol.190442–456. 10.1111/j.1469-8137.2010.03451.x
47
ListerR.GregoryB. D.EckerJ. R. (2009). Next is now: new technologies for sequencing of genomes, transcriptomes, and beyond.Curr. Opin. Plant Biol.12107–118. 10.1016/j.pbi.2008.11.004
48
LiuZ.KumariS.ZhangL.ZhengY.WareD. (2012). Characterization of miRNAs in response to short-term waterlogging in three inbred lines of Zea mays.PLOS ONE7:e39786. 10.1371/journal.pone.0039786
49
LivakK. J.SchmittgenT. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method.Methods25402–408. 10.1006/meth.2001.1262
50
LoretiE.Van VeenH.PerataP. (2016). Plant responses to flooding stress.Curr. Opin. Plant Biol.3364–71. 10.1016/j.pbi.2016.06.005
51
MayP.LiaoW.WuY.ShuaiB.McCombieW. R.ZhangM. Q.et al (2013). The effects of carbon dioxide and temperature on microRNA expression in Arabidopsis development.Nat. Commun.4:2145. 10.1038/ncomms3145
52
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
53
MorrisM.RogersS. M. (2014). Integrating phenotypic plasticity within an ecological genomics framework: recent insights from the genomics, evolution, ecology, and fitness of plasticity.Adv. Exp. Med. Biol.78173–105. 10.1007/978-94-007-7347-9_5
54
MutwilM.RuprechtC.GiorgiF. M.BringmannM.UsadelB.PerssonS. (2009). Transcriptional wiring of cell wall-related genes in Arabidopsis.Mol. Plant21015–1024. 10.1093/mp/ssp055
55
NagA.JackT. (2010). Chapter twelve-sculpting the flower; the role of microRNAs in flower development.Curr. Top. Dev. Biol.91349–378. 10.1016/S0070-2153(10)91012-0
56
NielsenS. L.NielsenH. D. (2006). Pigments, photosynthesis and photoinhibition in two amphibious plants: consequences of varying carbon availability.New Phytol.170311–319. 10.1111/j.1469-8137.2006.01670.x
57
NishiuchiS.YamauchiT.TakahashiH.KotulaL.NakazonoM. (2012). Mechanisms for coping with submergence and waterlogging in rice.Rice5:2. 10.1186/1939-8433-5-2
58
NomanA.FahadS.AqeelM.AliU.AnwarS.BalochS. K.et al (2017). miRNAs: major modulators for crop growth and development under abiotic stresses.Biotechnol. Lett.391–16. 10.1007/s10529-017-2302-9
59
OmidbakhshfardM. A.ProostS.FujikuraU.Mueller-RoeberB. (2015). Growth-regulating factors (GRFs): a small transcription factor family with important functions in plant biology.Mol. Plant8998–1010. 10.1016/j.molp.2015.01.013
60
PalatnikJ. F.AllenE.WuX.SchommerC.SchwabR.CarringtonJ. C.et al (2003). Control of leaf morphogenesis by microRNAs.Nature425257–263. 10.1038/nature01958
61
PanX.GengY.SosaA.ZhangW.LiB.ChenJ. (2007). Invasive Alternanthera philoxeroides: biology, ecology and management.Acta Phytotaxon. Sin.45884–900. 10.1360/aps06134
62
ParkJ. H.ShinC. (2014). MicroRNA-directed cleavage of targets: mechanism and experimental approaches.BMB Rep.47417–423. 10.5483/BMBRep.2014.47.8.109
63
PeetersA. J.CoxM. C.BenschopJ. J.VreeburgR. A.BouJ.VoesenekL. A. (2002). Submergence research using Rumex palustris as a model; looking back and going forward.J. Exp. Bot.53391–398. 10.1093/jexbot/53.368.391
64
PeiH.MaN.ChenJ.ZhengY.TianJ.LiJ.et al (2013). Integrative analysis of miRNA and mRNA profiles in response to ethylene in rose petals during flower opening.PLOS ONE8:e64290. 10.1371/journal.pone.0064290
65
PigliucciM. (2005). Evolution of phenotypic plasticity: Where are we going now?Trends Ecol. Evol.20481–486. 10.1016/j.tree.2005.06.001
66
RenY.ChenL.ZhangY.KangX.ZhangZ.WangY. (2012). Identification of novel and conserved Populus tomentosa microRNA as components of a response to water stress.Funct. Integr. Genomics12327–339. 10.1007/s10142-012-0271-6
67
RichardsC. L.AlonsoC.BeckerC.BossdorfO.BucherE.Colome-TatcheM.et al (2017). Ecological plant epigenetics: evidence from model and non-model species, and the way forward. 201576–1590. 10.1111/ele.12858
68
RichardsC. L.BossdorfO.MuthN. Z.GurevitchJ.PigliucciM. (2006). Jack of all trades, master of some? On the role of phenotypic plasticity in plant invasions.Ecol. Lett.9981–993. 10.1111/j.1461-0248.2006.00950.x
69
Rubio-SomozaI.WeigelD. (2011). MicroRNA networks and developmental plasticity in plants.Trends Plant Sci.16258–264. 10.1016/j.tplants.2011.03.001
70
SchlichtingC. D. (1986). The evolution of phenotypic plasticity in plants.Annu. Rev. Ecol. Evol. Syst.17667–693. 10.1146/annurev.es.17.110186.003315
71
SchlichtingC. D.SmithH. (2002). Phenotypic plasticity: linking molecular mechanisms with evolutionary outcomes.Evol. Ecol.16189–211. 10.1023/A:1019624425971
72
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
73
SedbrookJ. C.CarrollK. L.HungK. F.MassonP. H.SomervilleC. R. (2002). The Arabidopsis SKU5 gene encodes an extracellular glycosyl phosphatidylinositol–anchored glycoprotein involved in directional root growth.Plant Cell141635–1648. 10.1105/tpc.002360
74
SultanS. E. (2000). Phenotypic plasticity for plant development, function and life history.Trends Plant Sci.5537–542. 10.1016/S1360-1385(00)01797-0
75
SultanS. E. (2007). Development in context: the timely emergence of eco-devo.Trends Ecol. Evol.22575–582. 10.1016/j.tree.2007.06.014
76
SultanS. E. (2010). Plant developmental responses to the environment: eco-devo insights.Curr. Opin. Plant Biol.1396–101. 10.1016/j.pbi.2009.09.021
77
SunkarR.LiY. F.JagadeeswaranG. (2012). Functions of microRNAs in plant stress responses.Trends Plant Sci.17196–203. 10.1016/j.tplants.2012.01.010
78
SupekF.BošnjakM.ŠkuncaN.ŠmucT. (2011). REVIGO summarizes and visualizes long lists of gene ontology terms.PLOS ONE6:e21800. 10.1371/journal.pone.0021800
79
Van Der StraetenD.ZhouZ.PrinsenE.Van OnckelenH. A.Van MontaguM. C. (2001). A comparative molecular-physiological study of submergence response in lowland and deepwater rice.Plant Physiol.125955–968. 10.1104/pp.125.2.955
80
Van VeenH.MustrophA.BardingG. A.EijkM. H. V.WelschenevertmanR. A. M.PedersenO.et al (2013). Two Rumex species from contrasting hydrological niches regulate flooding tolerance through distinct mechanisms.Plant Cell254691–4707. 10.1105/tpc.113.119016
81
VoesenekL. A. C. J.Bailey-SerresJ. (2015). Flood adaptive traits and processes: an overview.New Phytol.20657–73. 10.1111/nph.13209
82
VoesenekL. A. C. J.BenschopJ. J.BouJ.CoxM. C. H.GroeneveldH. W.MillenaarF. F.et al (2003). Interactions between plant hormones regulate submergence-induced shoot elongation in the flooding-tolerant dicot Rumex palustris.Ann. Bot.91205–211. 10.1093/aob/mcf116
83
VoesenekL. A. C. J.van VeenH.SasidharanR. (2014). Learning from nature: the use of non-model species to identify novel acclimations to flooding stress.AoB Plants6:plu016. 10.1093/aobpla/plu016
84
WangJ. W. (2014). Regulation of flowering time by the miR156-mediated age pathway.J. Exp. Bot.654723–4730. 10.1093/jxb/eru246
85
WangL.LiX. R.LianH.NiD. A.HeY. K.ChenX. Y.et al (2010). Evidence that high activity of vacuolar invertase is required for cotton fiber and Arabidopsis root elongation through osmotic dependent and independent pathways, respectively.Plant Physiol.154744–756. 10.1104/pp.110.162487
86
WangJ. W.WangL. J.MaoY. B.CaiW. J.XueH. W.ChenX. Y. (2005). Control of root cap formation by microRNA-targeted auxin response factors in Arabidopsis.Plant Cell172204–2216. 10.1105/tpc.105.033076
87
YamauchiT.ShimamuraS.NakazonoM.MochizukiT. (2013). Aerenchyma formation in crop species: a review.Field Crops Res.1528–16. 10.1016/j.fcr.2012.12.008
88
YeJ.FangL.ZhengH.ZhangY.ChenJ.ZhangZ.et al (2006). WEGO: a web tool for plotting GO annotations.Nucleic Acids Res.34W293–W297. 10.1093/nar/gkl031
89
ZhaiL.LiuZ.ZouX.JiangY.QiuF.ZhengY.et al (2013). Genome-wide identification and analysis of microRNA responding to long-term waterlogging in crown roots of maize seedlings.Physiol. Plant.147181–193. 10.1111/j.1399-3054.2012.01653.x
90
ZhangB. (2015). MicroRNA: a new target for improving plant tolerance to abiotic stress.J. Exp. Bot.661749–1761. 10.1093/jxb/erv013
91
ZhangZ.WeiL.ZouX.TaoY.LiuZ.ZhengY. (2008). Submergence-responsive microRNAs are potentially involved in the regulation of morphological and metabolic adaptations in maize root cells.Ann. Bot.102509–519. 10.1093/aob/mcn129
92
ZhangZ.ZhangX. (2012). Argonautes compete for miR165/166 to regulate shoot apical meristem development.Curr. Opin. Plant Biol.15652–658. 10.1016/j.pbi.2012.05.007
93
ZhengQ.WangX. J. (2008). GOEAST: a web-based software toolkit for Gene Ontology enrichment analysis.Nucleic Acids Res.36W358–W363. 10.1093/nar/gkn276
Summary
Keywords
microRNA, dynamic expression, time course analysis, phenotypic plasticity, Alternanthera philoxeroides
Citation
Li G, Deng Y, Geng Y, Zhou C, Wang Y, Zhang W, Song Z, Gao L and Yang J (2017) Differentially Expressed microRNAs and Target Genes Associated with Plastic Internode Elongation in Alternanthera philoxeroides in Contrasting Hydrological Habitats. Front. Plant Sci. 8:2078. doi: 10.3389/fpls.2017.02078
Received
09 August 2017
Accepted
21 November 2017
Published
05 December 2017
Volume
8 - 2017
Edited by
Alejandra A. Covarrubias, Universidad Nacional Autónoma de México, Mexico
Reviewed by
Oswaldo Valdes-Lopez, Universidad Nacional Autónoma de México, Mexico; Francesco Licausi, Sant’Anna School of Advanced Studies, Italy
Updates
Copyright
© 2017 Li, Deng, Geng, Zhou, Wang, Zhang, Song, Gao and Yang.
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) or licensor 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: Lexuan Gao, lexuangao@hotmail.com Ji Yang, jiyang@fudan.edu.cn
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.