Next Article in Journal
Agrobacterium-Mediated Genetic Transformation of the Medicinal Plant Veratrum dahuricum
Previous Article in Journal
Modulation of Cadmium Tolerance in Rice: Insight into Vanillic Acid-Induced Upregulation of Antioxidant Defense and Glyoxalase Systems
Previous Article in Special Issue
The Use of Nitrogen and Its Regulation in Cereals: Structural Genes, Transcription Factors, and the Role of miRNAs
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 9
Issue 2
10.3390/plants9020189
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transgenerational Effects of Water-Deficit and Heat Stress on Germination and Seedling Vigour—New Insights from Durum Wheat microRNAs

by
Haipei Liu
*,
Amanda J. Able
and
Jason A. Able
School of Agriculture, Food & Wine, Waite Research Institute, The University of Adelaide, Urrbrae SA5064, Australia
*
Author to whom correspondence should be addressed.
Plants 2020, 9(2), 189; https://doi.org/10.3390/plants9020189
Submission received: 14 January 2020 / Revised: 3 February 2020 / Accepted: 3 February 2020 / Published: 4 February 2020
(This article belongs to the Special Issue Small RNAs in Crop Improvement and Breeding)
Browse Figures
Review Reports Versions Notes

Abstract

:
Water deficiency and heat stress can severely limit crop production and quality. Stress imposed on the parents during reproduction could have transgenerational effects on their progeny. Seeds with different origins can vary significantly in their germination and early growth. Here, we investigated how water-deficit and heat stress on parental durum wheat plants affected seedling establishment of the subsequent generation. One stress-tolerant and one stress-sensitive Australian durum genotype were used. Seeds were collected from parents with or without exposure to stress during reproduction. Generally, stress on the previous generation negatively affected seed germination and seedling vigour, but to a lesser extent in the tolerant variety. Small RNA sequencing utilising the new durum genome assembly revealed significant differences in microRNA (miRNA) expression in the two genotypes. A bioinformatics approach was used to identify multiple miRNA targets which have critical molecular functions in stress adaptation and plant development and could therefore contribute to the phenotypic differences observed. Our data provide the first confirmation of the transgenerational effects of reproductive-stage stress on germination and seedling establishment in durum wheat. New insights gained on the epigenetic level indicate that durum miRNAs could be key factors in optimising seed vigour for breeding superior germplasm and/or varieties.

1. Introduction

Global food demand is constantly increasing with the rise in our world population. However, even with a finite amount of arable land, agricultural production is still expected to provide more food and fibre, as well as to adapt to the environmental challenges (i.e., climate resilience). High quality seeds or grains are therefore essential for sustainable agriculture and food security. Cereal crop production starts with seed germination to establish vigorous healthy seedlings in the field. Successful germination and robust seedling growth are critical early stages of crop life to ensure reliable canopy establishment and production efficiency [1,2]. Seeds with high vigour are expected to germinate rapidly and uniformly, thus warranting field emergence with minimum delay after sowing [3]. Seedlings with high vigour are expected to be more resource-use-efficient and to have competitive advantages over unwanted weed species [1,2].
For cereal crops, successful germination and vigorous seedling emergence relies on high quality seeds. Durum wheat (Triticum turgidum L. ssp. durum) is a vital cereal crop that provides substantial economic output in the Mediterranean region [4] due to its unique use in multiple food products such as pasta and couscous. Durum is better adapted than common hexaploid wheat to semi-arid environments, but durum wheat production and quality suffer greatly from several abiotic stress constraints including water deficiency and heat [5,6]. Many studies have investigated the impact of water-deficit stress and heat stress on yield and quality, with emphasis on grain traits, such as protein content and carbohydrate composition, that are directly associated with end use (e.g., pasta-making) [7,8,9]. However, there is no evidence of research to date that has focused on investigating the transgenerational effects of stress on seed vigour and seedling establishment in durum wheat. Previous studies in other crops with agronomic significance have shown that the effects of environmental stresses can be transmitted to the offspring with significant impact [10,11,12,13]. A recent study in soybean reported that exposing parental plants to drought stress significantly reduces germination and seedling growth of the F1 generation grown under optimum conditions [11]. In rice, transient heat stress shortly after fertilisation can either promote or supress seed germination rate, depending on the timing of the stress [12]. In canola, drought stress on the mother plants actually increased the progeny seedling vigour, independent of the maternal yield performance [10]. Thus, establishing the transgenerational effects of abiotic stress in durum wheat is of relevance, particular for the broader wheat breeding programs that exist globally, where such information could be translated to enhance breeding strategies of the future.
Seeds originating from different parents not only inherit their full genetic components, but also their epigenetic marks [14]. Rapidly accumulating evidence shows that epigenetic modification can have transgenerational effects which have an impact on adaptive traits in plants [15,16]. Small RNAs (sRNAs) are a class of non-coding RNAs with essential regulatory roles in the epigenetic phenomena [17]. In plants, the major sRNA types include microRNAs (miRNAs) and small-interfering RNAs (siRNAs) [18]. Plant miRNAs can precisely reprogram the expression of their protein-coding target genes by inducing mRNA degradation or translational inhibition through post-transcriptional gene silencing [17,18]. In crops, miRNAs can rapidly respond to environmental stress and act as mobile signal molecules to help restore cellular homeostasis [19]. Numerous studies have shown that miRNAs are important for crop plants in growth, development and stress adaptation [20,21,22,23,24,25,26,27,28]. The involvement of miRNAs in stress memory has also been demonstrated as part of the within- or cross-generation phenotypic plasticity [16,29]. Next generation sequencing (NGS) or microarray technologies have been used to study the miRNA profile in durum wheat under various types of abiotic stress [20,21,22,30,31,32]. However, most efforts have focused on the same generation of stress occurrence. We lack the knowledge of the transgenerational impact of stress on the miRNA transcriptome (miRNAome) in progeny, and its possible implication associated with early seedling establishment.
In this study, we therefore aimed to elucidate the transgenerational impact of water-deficit and heat stress on germination and seedling vigour of stress-tolerant and stress-sensitive durum genotypes. To understand the molecular mechanism underlying altered genotypic performance, sRNA sequencing was used to profile the miRNAome along with in silico characterisation of their functional target genes. To our knowledge, this is the first description of the transgenerational effects of abiotic stress in durum wheat integrated with genome-wide sRNA sequencing. Moreover, the current study is the first to employ the new durum genome assembly (Svevo.v1, available 2019 [33]) for epigenetic studies, providing a more accurate, in-depth analysis on the whole-genome scale.

2. Results

2.1. Evaluation of Germination and Seedling Establishment

Two Australian durum wheat genotypes were used in this study in controlled environment facilities. DBA Aurora is a leading Australian commercial variety that has shown high tolerance to pre-anthesis water-deficit stress and post-anthesis heat stress [8]. Line 6 (L6) is an advanced University of Adelaide breeding line (UAD1301020-8) that is sensitive to the combination of water-deficit and heat stress [8]. The seeds used in this study were collected from the parents that were in the control group and the water-deficit and heat stress group of both genotypes [8]. In total, there were four different seed groups: AuCG (seeds from DBA Aurora control group), AuWH (seeds from DBA Aurora water-deficit and heat stress group), L6CG (seeds from L6 control group), and L6WH (seeds from L6 water-deficit and heat stress group).
For seed germination measurements (Table 1), genotype had a significant impact (p < 0.05) on the germination potential (Gp), germination rate (Gr), and germination index (Gi) without interaction with the parent treatment factor. Overall, DBA Aurora had better germination performance than L6. In particular, AuCG and AuWH almost had double the germination potential of L6CG and L6WH. Parent treatment only had a significant impact (p < 0.05) on Gp and Gi without interaction with the genotype. Overall, water-deficit and heat stress treatment on the parents negatively affected the germination performance of the next generation. Coleoptile length was not significantly affected by any of the factors. Interestingly, a significant interaction (p = 0.005) was detected between genotype and parent treatment for young root length. For DBA Aurora, the parent stress treatment did not significantly affect the young root of the progenies. However, for L6, stress on the parents significantly reduced the length of the young roots.
For the seedling vigour index I (SVI-I) and seedling vigour index II (SVI-II) (Table 2), both genotype and parent treatment had a significant impact (p < 0.05) without interaction with each other. Overall DBA Aurora outperformed L6. Water-deficit and heat stress (WH) treatment of the parents negatively impacted both genotypes, but to a greater extent in L6. For other seedling measurements (Table 2), genotype had a significant impact (p < 0.05) on shoot fresh weight, shoot dry weight, and seedling height, with no interaction with the parent treatment factor. Overall, DBA Aurora again displayed a better performance for these traits when compared to L6. Parent treatment also had a significant impact (p < 0.05) on shoot fresh weight and seedling height without interaction with the genotype factor. The WH treatment of the parents significantly reduced the shoot fresh weight and seedling height of the next generation (p < 0.05). Interestingly, a significant interaction (p < 0.05) was detected between genotype and parent treatment for all root traits (root fresh weight, root dry weight and root length). For DBA Aurora, the WH treatment of the parents had no impact on the root fresh weight and dry weight of the next generation seedlings. For L6, the WH treatment of the parents significantly reduced the root fresh weight and dry weight. The WH treatment of the parents promoted the root length in both genotypes, but to a greater extent in DBA Aurora.

2.2. The microRNAome in the Progenies

Four sRNA libraries (TCG, TWH, SCG and SWH) were constructed and sequenced from the four seedling groups (AuCG, AuWH, L6CG and L6WH), respectively. A total of 83,991,830 raw reads were obtained (Table 3). Approximately 52.04 million clean sRNA reads within the 18–25 nt range remained after data processing. Among these, 26.16 million non-redundant unique sRNA reads were obtained. The size distribution of sRNA reads exhibited a two-peak pattern (a weak peak at 21 nt and a strong peak at 24 nt) (Figure 1), a classic profile expected for cleaved products of plant Dicer-like enzymes [20,22,34].
A total of 2202 MIR-miRNA entries (different combinations of MIR loci origin and mature miRNA products) were identified (Table S1). All MIR-miRNA entries were categorised into five groups (G1–G5), with G1–G4 representing conserved miRNAs and G5 representing novel miRNAs (Table 4). The sequence reads that mapped to the durum genome and ESTs (expressed sequence tags) but their extended genome sequences could not form secondary hairpins (G3) had the highest number of unique miRNAs. The group that had neither aligned reference pre-miRNAs nor sequencing reads that could be further mapped to the genome (G4) had the lowest number of unique miRNAs (Table 4). In Table S1, MIR-miRNA entries with the same pre-miRNA index number within each group share the same pre-miRNA sequence. As one pre-miRNA can potentially produce multiple unique mature miRNAs, each unique miRNA is labelled with its own miRNA index number within each group. The pre-miRNA and miRNA index numbers in Table S1 can be linked to Table 4. For Group 2 to Group 4, the total number of pre-miRNAs and the total number of unique mature miRNAs correspond to the highest pre-miRNA index number and the highest miRNA index number within each group in Table S1 (e.g., the last MIR-miRNA entry for Group 2 in Table S1 is indexed as 446–483). However, the highest miRNA index number in Group 1 (Table S1) is different to the total number of unique miRNAs in Group 1 (Table 4) as the same mature miRNA (unique miRNA) can be produced from different pre-miRNAs with different index numbers.
Among all entries (Table S1), 200 miRNAs were highly expressed (normalised reads number of the reported miRNA is higher than the average of the dataset), including 190 conserved miRNAs and 10 novel miRNAs. A total of 1435 miRNAs showed a medium expression level (normalised reads number of the reported miRNA is higher than 10 but less than the average of the dataset), including 1002 conserved miRNAs and 433 novel miRNAs. A total of 567 miRNA entries were lowly expressed (normalised reads number of the reported miRNA is less than 10), including 525 conserved miRNAs 42 novel miRNAs. In total, these 2202 MIR-miRNA entries represent 2177 unique mature miRNA products among four libraries, as shown in the Venn diagram (Figure 2). Collectively, the conserved miRNAs identified in this study represent 70 MIR families (Table S2). The conservation profiles of all miRNAs across the reference plant species (Table S3) in miRBase are also shown (Figure 3).

2.3. Differentially Expressed miRNAs (DEMs)

Differentially expressed miRNAs (DEMs) were identified based on their normalised reads number. A total of 1368 miRNAs (Table S4) showed differential expression across four libraries (TCG vs TWH vs SCG vs SWH). A heatmap was compiled to show DEMs across four libraries using their log10 (normalised count) value (Figure S1). Some miRNAs exhibited a genotype-specific pattern (Table S5). A total of 144 miRNA were found to be exclusively expressed in DBA Aurora (TCG and/or TWH), while 139 miRNAs were only expressed in L6 (SCG and/or SWH). A number of miRNAs also showed specificity to parent treatment (Table S5). A total of 129 miRNAs were only detected in the CG libraries (TCG and/or SCG) while 72 miRNAs were only found in the WH libraries (TWH and/or SWH).
When comparisons were made between two libraries within the genotype or within the treatment source (TWH vs TCG, SWH vs SCG, SCG vs TCG, SWH vs TWH), DEMs were identified using Chi-squared 2 × 2 test and Fisher’s test (see Materials and Methods). The number of up- and down-regulated DEMs identified in each comparison using different p values as thresholds are shown in Figure 4. Interestingly, a genotypic pattern can be observed here for the number of DEMs. For the stress-tolerant genotype DBA Aurora (TWH vs TCG), the number of miRNAs that were up-regulated was higher than the ones that were down-regulated, irrespective of the p value used. However, for the stress-sensitive genotype L6 (SWH vs SCG), an opposite trend was observed, where there were more miRNAs that were down-regulated than up-regulated. When the same parent treatment source was compared between two genotypes (SCG vs TCG, SWH vs TWH), there are more down-regulated miRNAs than up-regulated miRNAs (except for SCG vs TCG at p < 0.1). The expression levels of all DEMs are shown in Table S6. Highly expressed DEMs with a |log2 (fold change)| > 1 in each comparison (p < 0.05) are listed in Table 5.

2.4. Target Prediction and Enrichment Analysis

The target genes were predicted for highly expressed DEMs using in silico analysis with the recently constructed durum wheat genome (see Materials and Methods). A total of 311,964 transcripts (target genes) were identified as potential targets for 189 DEMs (Table S7). Gene ontology (GO) analysis was conducted to classify the biochemical functions of the predicted target genes. GO annotation grouped all targets in the following three GO categories (i) Biological Process (BP), (ii) Molecular Function (MF), and (iii) Cellular Component (CC). GO enrichment analysis of the predicted targets was performed (Table S8). The most represented GO terms are shown in Figure 5. Notably, a few biological processes with importance are highly represented, including the regulation of transcription, DNA-templating (GO:0006355), protein phosphorylation (GO:0006468), the oxidation-reduction process (GO:0055114), transmembrane transport (GO:0055085) and signal transduction (GO:0007165). Highly represented cellular components include the plasma membrane (GO:0005886), cytoplasm (GO:0005737), integral components of the membrane (GO:0016021), cytosol (GO:0005829) and chloroplast (GO:0009507). In the molecular function category, binding activities are highly represented, including protein binding (GO:0005515), ATP binding (GO:0005524), DNA binding (GO:0003677) and metal ion binding (GO:0046872).
KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis was performed as previously described [35] to identify the biochemical pathways with which durum target genes were associated (Table S9). The top-ranked enriched pathways are shown in Figure 6. In particular, a total of 2406 targets were matched to protein processing in endoplasmic reticulum (ko04141), 2248 targets were matched to starch and sucrose metabolism (ko00500), 2171 targets were matched to the MAPK signalling pathway (ko04010), 1191 target were matched to the spliceosome (ko03040), and 1720 targets were matched to ABC transporters (Table S9).

3. Discussion

Cereal crops are often challenged by multiple biotic and abiotic factors during their life cycle, affecting both production volume and quality [8,36,37,38,39]. The detrimental effects on the crop plants are often exacerbated when several stress factors are combined, simultaneously or sequentially [8,40]. Increasing evidence has shown that the effects of environmental stress on the parental plants can be carried over to the subsequent generation(s). In the attempt to achieve sustainable crop production, it is essential that we understand the impact of abiotic stress within and across generations, as well as the fundamental networks that underpin crop performance. The influence of the maternal environment on the offspring early growth has rarely been explored in durum wheat before the current study.
For the two genotypes used in this study (DBA Aurora and L6), the parental exposure of water-deficit and heat stress during the reproductive stage significantly reduced the germination potential and germination index. These traits reflect the time-course and the cumulative status of seed germination. Similarly, in soybean, water shortage during the reproductive stage had a pronounced transgenerational effect on seed germination performance [11]. The cumulative percentage of germination negatively correlated with the severity of the drought stress that was imposed on the parents. Germination rate was also significantly reduced in the F1 generation, with reduced seed number and individual seed weight in both soybean cultivars. However, in our study, the germination rate was not affected by the parental stress treatment factor. Moreover, the genotype had a significant impact on all germination traits measured (except for the coleoptile length), with superior performance observed for DBA Aurora. In our previous study, the effect of stress on the parental plants, including physiological responses, yield component and particularly grain traits, was well documented [8]. DBA Aurora produced bigger grains that were full-looking under water-deficit stress and heat stress. The 1000-grain weight of DBA Aurora water-deficit plus heat stress group was significantly higher than that of the control group. On the other hand, L6 had significantly smaller seeds with reduced weight under the effect of such stress (that is, significantly lower 1000-grain weight compared with its control group). For cereal species, high temperature and water shortage during reproduction often leads to an accelerated grain growth rate and a shortened grain development period [40]. As a result, smaller grains with poor grain filling are often produced. However, elite stress-tolerant varieties like DBA Aurora could mitigate the adverse impact of stress and still manage to produce bigger grains, which could (at least in part) contribute to its superior germination performance over L6. It is also worth noting that there is no significant difference between the 1000-grain weight of L6WH seeds and AuWH seeds, suggesting that other factors could contribute to the transgenerational effect of stress rather than grain size.
Similar to germination traits, early durum seedling growth of the offspring was also significantly affected by the parental stress exposure. Genotype seems to play a significant role here. Although seedling vigour was significantly reduced under the transgenerational effect of stress in both genotypes, a greater impact was always observed in L6. In particular, a significant interaction was found between genotype and parent treatment for all root traits measured in the young seedlings. Only L6 has shown reduced root weight due to parental stress exposure. Such genotypic differences in the transgenerational effects of stress have also been noted before [41]. The findings suggest that breeders and growers should take care in selecting varieties under different parental environments for crossing or seed production (this is often overlooked, with a ‘mother stock’ source used for such purposes). Furthermore, the stress exposure of the parents increased the young seedling root length for both genotypes (with a higher increase in DBA Aurora), suggesting that root elongation could be promoted through transgenerational stress exposure irrespective of genotypes. Establishing a longer, deeper root system at earlier developmental stages might be an adaptive mechanism to prepare the offspring to improve water and nutrient uptake in any possible future stress occurrence [11]. Previous research has suggested that the impact of parental stress on early seedling growth is often associated with seed storage reserves such as carbohydrates and proteins. However, in our case, both the DBA Aurora and the L6 parents produced grains with increased grain protein content and reduced starch content under the same stress without any significant genotypic difference [8]. In the miRNA-target analysis, starch and sucrose metabolism came second for the most-enriched KEGG pathways (Table S9). The 1639 targets enriched for this pathway included many enzymes, such as sucrose synthase, starch synthase, beta-glucosidase, sucrose phosphatase, and isoamylase (Table S8). Altered enzyme activities regulated by differentially expressed miRNAs can affect how storage reserves were utilised during germination and early growth. The results suggest that, apart from seed quality traits being affected by stressful paternal environments, epigenetic mechanisms such as miRNA-target modules could be in play during the transmission across generations [11].
As a central class of epigenetic regulators, miRNAs take part in a wide range of cellular biological processes [17,18,42]. Hence, the identification and annotation of the miRNAome enables the epigenetic networks underlying the transgenerational effect of stress to be illustrated. The combination of NGS and bioinformatics analysis using genome reference has been the most effective approach to characterise the miRNA population on a genome-wide scale [17]. Such technologies have led to the complete discovery of small RNAs in many plant species on an unprecedented scale. However, despite a few reports in durum wheat (including our previous studies) [20,22,30,34], the discovery of novel miRNAs and a fully comprehensive miRNAome description in this important cereal has not been previously achieved due to incomplete genome information. Only recently, the ground-breaking genome assembly of the modern durum cultivar Svevo was made available [33]. The complete Triticum turgidum genome presented us with a timely opportunity to harness additional information from sRNA sequencing datasets. For the first time, to our knowledge, a systematic and complete archive of durum miRNAs has now been provided. In total, we successfully identified 1717 conserved and 485 novel miRNA entries matched with their MIR loci origins. This number exceeded all previous reports in durum wheat. Interestingly, among all the conserved miRNA entries, the top plant species with the highest miRBase conservation match was Glycine max, followed by Oryza sativa. Hexaploid wheat (Triticum aestivum), with tetraploid wheat as one of its genome progenitors, only came third in conservation ranking. This is likely due to soybean miRNAs being very well-documented (684 precursors, 756 mature) in the miRBase when compared to bread wheat (122 precursors, 125 mature), which is not surprising considering that the Glycine max genome was well-annotated ahead of all wheat species sequenced. Notably, there is currently only one Triticum turgidum miRNA entry (ttu-miR160, MI0016454) deposited in the miRBase [43]. It could be foreseen that this number will increase rapidly with the first durum genome assembly now being available to the scientific community.
Plant miRNAs can act rapidly in response to both environmental and developmental signals [17]. In particular, various studies have described the functional roles of miRNA-target modules in combatting water deficiency and extreme temperatures in cereal crops [8,23,24]. In our previous research, we have demonstrated that the expression of durum miRNAs can change differently with different timings of stress adaption within the life cycle [23,24]. Such patterns were also genotype-dependent among stress-tolerant and -sensitive varieties, suggesting that miRNA-mediated regulatory modules could contribute to the genotypic differences in the physiological performance, ultimately affecting yield production [23]. The present study further illustrated that durum miRNAs behave in a genotypic manner under the effect of transgenerational stress. Overall, there were more up-regulated miRNAs (275, p < 0.05) than down-regulated miRNAs (181, p < 0.05) in the stress-tolerant genotype, but more down-regulated miRNAs (242, p < 0.05) than up-regulated (153, p < 0.05) miRNAs in the stress-sensitive genotype. Durum miRNAs that were induced by the transgenerational effects could repress the expression of protein-coding genes with negative regulatory roles, while reduced miRNA expression could promote the expression of target genes to increase protein products with positive functions. A higher number of induced miRNAs in the stress-tolerant genotype and higher number of repressed miRNAs in the stress-sensitive genotype suggest that contrasting mechanisms were preferably deployed in response to transgenerational stress, resulting in the difference in young seedling growth.
Investigations into the expression pattern of specific durum miRNAs, together with the molecular functions of their predicted targets, provide further information on the possible pathways activated in different genotypes. As an example, tae-miR1847-5p was significantly down-regulated in DBA Aurora (−0.90 fold) under the transgenerational effects of stress (Table S6). In silico analysis revealed that tae-miR1847 targets several CREB-binding protein-like family members located on the 6A and 6B chromosomes (Table S7). CREB-binding proteins (CBPs) are a type of transcriptional coactivator that can switch on the cAMP-dependent signalling pathways, which are key components employed by many hormone-regulated activities to modulate gene expression [44,45]. CBPs initiate the function of many classes of transcription factors through their intrinsic transferase activity to acetylate histones, which is a key process associated with chromatin remodelling [44]. Stress-related chromatin modification plays a critical role in transgenerational memory in plants [15]. The repressed expression of tae-miR1847 in the stress-tolerant genotype would increase the expression of its target CBPs, which could ultimately promote histone acetyltransferase-regulated transduction pathways. However, this pathway might not be in play in the stress-sensitive genotype, as tae-miR1847-tp was up-regulated (0.23-fold) under the transgenerational effects of stress. Another example of a miRNA-regulated module that could contribute to the differential seedling performance involved mitogen-activated protein kinase (MAPK). MAPK signalling is one of the top KEGG pathways enriched among the predicted targets. MAPK cascades are a universal signal transduction mechanism that play diverse roles in connecting intra- and extra-cellular crosstalk in eukaryotes [46]. MAPKs are involved in virtually all the important aspects of plant development and response to environmental signals. In particular, MAPKs facilitate plant adaptation and survival under many abiotic stresses, including drought, heat and salt stress, via coordinating the transduction of secondary messengers and hormone signals [46,47]. Here, a few durum miRNAs that target MAPKs exhibited genotype-dependent expression patterns. For example, conserved miRNAs such as ttu-miR160, hvu-miR444b, and tae-miR1847-5p were all down-regulated (−0.42-fold, −0.45-fold and −0.90-fold, respectively) in DBA Aurora under the effect of transgenerational stress (Table S6). This promotes the expression of MAPKs in the stress-tolerant variety to facilitate signal transduction. However, the expression of these miRNAs was either unchanged (ttu-miR160 and hvu-miR444b) or increased (tae-miR1847-5p) in L6, suggesting that miRNA-MAPK modules were not activated in the stress-sensitive genotype. Research in bread wheat has also demonstrated the involvement of MAPK in transgenerational tolerance to heat stress [48]. Increased expression of MAPK in the progenies of heat-stressed parents contributed to the enhanced activation of Ca2+ signalling pathways to facilitate stress adaptation [48]. Moreover, stress-induced MAPK could promote root tip growth through controlling the organisation of the actin cytoskeleton [49]. Durum miRNA-promoted MAPK expression could possibly contribute to the adaptive root traits observed in the stress-tolerant genotype. Future research should investigate the co-expression of durum miRNAs and MAPK family members in different tissues types to further our understanding of miRNA-MAPK modules in transgenerational stress adaption.
The spliceosome (ko03040) is another major KEGG pathway enriched for the targets of differentially expressed durum miRNAs. The pre-mRNA splicing machinery is governed by the spliceosome, a large ribonucleoprotein complex that recognises the pre-mRNA splice sites to remove introns and join exons [50]. Research has shown that alternative splicing (AS), which enables the production of multiple mRNA isoforms from one single gene, is closely associated with abiotic stress adaption [51,52]. The AS process provides a fine-tuning mechanism to regulate the levels of active and inactive isoforms of stress-associated proteins. The specificity of splice-site selection in AS largely relies on serine/arginine-rich (SR) proteins (splicing enhancer) and heterogeneous nuclear ribonucleoproteins (hnRNPs, splicing suppressor) [53]. In our study, a large number of SR and hnRNP family members were targeted by differentially expressed durum miRNAs. Some miRNAs exhibited a genotype-dependent expression, while others seem to have universal regulatory roles. For example, ttu-miR160 targets several hnRNP-encoding genes located on the 4A and 4B chromosomes (Table S7). It also targets thirteen SR splicing factor-coding genes located on chromosome 6A and 6B. ttu-miR160 was down-regulated in DBA Aurora (−0.42-fold) in response to transgenerational stress, indicating that the expression of SR proteins and hnRNPs were possibly increased to coordinate the AS process. ttu-miR160-regulated AS was not in play in the sensitive genotype, based on the unaltered miRNA expression. However, a few other miRNA-regulated modules were active in both genotypes. For example, tae-miR9772, ata-miR396c-5p, and a miR396c variant (sbi-miR396c_R+1) all target hnRNP 1-like genes. The expression of these miRNAs was down-regulated in both genotypes, suggesting that hnRNP-modulated splicing events were activated by miRNAs without genotypic specificity. Furthermore, stress-induced AS can regulate the processing of non-coding RNA transcripts, so as to control the production of different pre-miRNA hairpins [54]. In barley, heat stress induced the splicing of introns hosting the miRNA primary transcripts (pri-miRNAs) of miR160a and miR5175a, which was correlated with the accumulation of mature miRNA products [55]. The results suggest that the durum miRNA-modulated AS mechanism could possibly provide a feedback loop to fine-tune the post-transcriptional regulation of miRNA precursor processing. Looking forward, the complexities in durum miRNA expression, spliceosomal protein activities, and the production of pri-miRNA isoforms require further exploration. In addition, it would be desirable for future work to study the upstream regulatory mechanisms of durum miRNAs so as to understand the underlying reasons for miRNA expression changes between stress-tolerant and stress-sensitive genotypes. Genomic analysis of the surrounding regions (such as transcription start sites and cis-regulatory motifs) of major stress-responsive miRNAs can compare the promoter features of the stress-tolerant and stress-sensitive durum genotypes to provide valuable information on transcriptional regulation of MIR genes. Our current research has paved the way for future investigation, having provided in-depth details of the MIR gene (genome ID, length, pre-miRNA sequence, secondary structure etc., Table S1).
In conclusion, the present study provides a valuable resource for an improved understanding of the transgenerational effects of water-deficit and heat stress on germination and seedling growth in durum wheat, as well as new insights of conserved and novel miRNAs that are involved in the epigenetic regulatory networks. Such knowledge provides a new foundation for durum researchers and breeders to use effective strategies for developing and selecting elite varieties that can mitigate or adapt to the negative impacts of stress that are transmitted across generations.

4. Materials and Methods

4.1. Plant Materials

Seeds of two Australian genotypes were used in this study. The seeds of both genotypes were collected from the parental plants in the control group and the water-deficit and heat stress group from the previous year’s experiment [8]. Briefly, water-deficit stress was applied to the parents from the booting stage by maintaining the soil water content at half of the field capacity. Heat stress was applied post-anthesis by exposing the plants to 37/27 °C for 24 h at 5, 15, 25, 35, and 45 days after flowering (a total of five stress episodes on the same plants). In total, four different seed groups were used in this study—AuCG (seeds from DBA Aurora control group), AuWH (seeds from DBA Aurora stress group), L6CG (seeds from L6 control group), and L6WH (seeds from L6 stress group).

4.2. Germination Test

Germination tests of the four seed groups were performed according to the International Rules for Seed Testing 2019 (https://www.seedtest.org). The seeds used for germination were randomly selected. The seeds were sterilised by rinsing with 70% ethanol for 5 min, followed by a wash with sterile distilled water. The seeds were germinated on filter paper in petri dishes containing sterile distilled water for 96 h at 22 °C. Three replicates per group were used with 30 seeds per replicate. germination potential (Gp, %) was determined as number of seeds germinated within 24 h/total number of seeds × 100%. The germination rate (Gr, %) was determined as number of seeds germinated within 96 h/total number of seeds × 100%. The germination index (Gi) was determined as ∑(nt/dt), where nt represents the number of germinated seeds on the t day, and dt represented germination days [56]. The length of coleoptile and young roots was recorded at 96 h.

4.3. Measurement of Seedling Establishment

For each group, 16 randomly selected germinated seeds were transferred into pots (one seedling per pot). Each pot contained 1.2 kg of N40 sand with 0.5% CaCO3, and nutrient solution was provided as previously described [8]. The seedlings were grown under controlled glasshouse conditions with day/night temperature of 24/18 °C and a photoperiod of 16/8 h. Seedling height, shoot fresh weight and root fresh weight were measured on 12-day old seedlings. Shoot dry weight and root dry weight were measured after overnight oven drying at 65 °C. Fresh shoots were collected and snap-frozen in liquid nitrogen for RNA extraction. Eight biological replicates were used for seedling measurements and eight biological replicates were used for RNA extraction. Seedling vigour index I (SVI-I) was determined as germination rate × seedling fresh weight in mg [57]. Seedling vigour index II (SVI-II) was determined as germination rate × seedling dry weight in mg [57]. Two-way-ANOVA was used to analyse the statistical significance (p < 0.05) for germination and seedling data using Genstat (19th Edition, VSN International).

4.4. Small RNA Sequencing

Total RNA was extracted from young seedlings using Tri reagent (Sigma-Aldrich, North Ryde, Australia) and treated with TURBO DNase (ThermoFisher Scientific, Scoresby, Australia) following the manufacturer’s instructions. The concentration and quality of the RNA samples were measured on a NanoDrop Lite spectrophotometer (ThermoFisher Scientific, Scoresby, Australia). RNA integrity was assessed by agarose gel electrophoresis and a Bioanalyzer (RIN > 9). Total RNA from eight biological replicates in each group was pooled for small RNA sequencing. Four small RNA libraries (TCG, TWH, SCG, SWH) were constructed using the TruSeq Small RNA Sample Prep Kit from the four groups AuCG, AuWH, L6CG, and L6WH, respectively. Single-end sequencing (50 bp) was performed on an Illumina HiSeq2500 by LC-Bio (Hangzhou, China). The sequencing datasets were deposited in NCBI GEO database (accession number GSE143454).

4.5. Analysis of Conserved and Novel miRNAs

Sequencing analysis was performed using the in-house ACGT101-miR program (LC Sciences, Houston, TX, USA). Raw sequencing reads were first processed by removing low-quality reads and adapter sequences to obtain clean sequences. The reads with a nucleotide (nt) length <18 nt or > 25 nt were removed. Common non-coding RNA families (rRNA, tRNA, snRNA and snoRNA), repeats and mRNA sequences were then discarded using RFam, Repbase and durum wheat NCBI mRNA entries as references. Unique sRNA sequences were then obtained for each library. Conserved mature miRNAs and 3p- or 5p-derived miRNA variants were identified by BLAST search against miRNA precursors in the miRBase (Release 22.1). MIR and miRNA sequences from common plant species in the miRBase were used as references (Table S3). Single mismatch within the sequence and length variation at both 5p and 3p were allowed in the alignment. Sequences mapped to the mature miRNA in the hairpin were identified as conserved mature miRNAs. Sequences mapped to the opposite arm of mature miRNA in the hairpin were identified as 5p- or 3p-derived variants. All the mapped miRNAs were aligned to the durum wheat genome (NCBI UID 3439611, assembly Svevo.v1) to determine their genomic location.
The remaining unmatched sRNA sequences were used to identify novel durum miRNAs. Sequences were BLASTed to the durum wheat genome. To identify miRNA precursors, secondary hairpin structures containing matched sequences were predicted using RNAfold (http://rna.tbi.univie.ac.at/cgi-bin/RNAfold.cgi) from 120 nt of flanking genome sequences. The following criteria were used for the prediction of the secondary structure as previously described [35,58]: (1) the length of the hairpin (two stems and the terminal loop) is ≥50 nt; (2) the length of the hairpin loop region is ≤200 nt; (3) the cut-off of free energy is ≤−15 kCal/mol; (4) the stem region of the predicted hairpin is ≥16 bp in length; (5) for each bulge in the stem, the nucleotide length is ≤12 nt; (6) for each bulge in the mature region of the predicted hairpin, the nucleotide length is ≤4 nt; (7) the number of biased matches in each bulge in the mature region is ≤2; (8) the number of biased bulges in the mature region is ≤2; (9) the number of mismatches in the mature region is ≤4; (10) the mature region is ≥12 bp in length; and (11) the mature region in the stem is ≥80%.
All miRNAs were categorised into five groups (G1 to G5). To analyse the miRNA expression profiles, data normalisation of the reads number was carried out as previously described [59]. Differentially expressed miRNAs (DEMs) were identified based on the normalised count. Chi-squared n × n test was used to detect DEMs when comparing TWH vs TCG vs SWH vs SCG. Chi-squared 2 × 2 test and Fisher’s test was used to detect DEMs when comparisons were made between two libraries (TWH vs TCG, SWH vs SCG, SCG vs TCG, SWH vs TWH). The cut-off for the p value was set as 0.01, 0.05 and 0.1 to generate Figure 4 and provide a total evaluation of the number of miRNAs that were up-/down-regulated. The cut-off for the p value was set as 0.05 for Tables S4 and S6 to identify miRNAs with significant differential expression. DEMs with a |log2 (Fold change)| > 1 in each comparison were identified and highlighted if they were up/down-regulated in Table S6.

4.6. Target Prediction and Functional GO Analysis

To predict protein-coding genes targeted by durum miRNAs, in silico analysis using the GSTAr was employed. Gene Ontology terms of the miRNAs targets were annotated according to their biological process, cellular component and molecular function. KEGG Pathways (http://www.genome.jp/kegg) of miRNA targets were also analysed. The statistical enrichment of GO terms and KEGG pathways was performed as previously described [35,60].

Supplementary Materials

The following are available online at https://www.mdpi.com/2223-7747/9/2/189/s1. Figure S1: Heatmap of differentially expressed miRNAs (DEMs) identified in comparison across four libraries (TCG vs TWH vs SCG vs SWH). The colour key (ranging from blue to red) represents the miRNA abundance using the log10 (normalised count) value, Table S1: Summary of 2202 MIR-miRNA entries identified in this study, Table S2: Conserved MIR families identified in this study, Table S3: List of plant species that were used as references in the miRBase, Table S4: List of miRNAs that were differentially expressed across TCG vs TWH vs SCG vs SWH, Table S5: List of miRNAs that were specifically expressed in DBA Aurora, L6, CG and WH seed sources, Table S6: List of differentially expressed miRNAs in comparisons made between two libraries (TWH vs TCG, SWH vs SCG, SCG vs TCG, SWH vs TWH), Table S7: Target genes predicted for differentially expressed miRNAs, Table S8: GO enrichment of predicted target genes in three categories: Biological Process (BP), Molecular Function (MF), and Cellular Component (CC), Table S9: KEGG pathway enrichment analysis of predicted target genes.

Author Contributions

Experimental design and conceptualization, H.L., A.J.A. and J.A.A.; Performing experiments, H.L.; Data analysis, H.L.; Writing—Original draft preparation, H.L.; Writing—Review & Editing, H.L., A.J.A. and J.A.A.; Funding Acquisition, H.L. and J.A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Australian Research Council, Project Number DE180100784.

Acknowledgments

We thank Nick Smart (South Australian Research and Development Corporation) for maintaining the glasshouse environment, and Dr Anh Pham (University of Adelaide) for setting up the pots.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Richards, R.A.; Rebetzke, G.J.; Condon, A.G.; van Herwaarden, A.F. Breeding opportunities for increasing the efficiency of water use and crop yield in temperate cereals. Crop Sci. 2002, 42, 111–121. [Google Scholar] [CrossRef] [PubMed]
  2. Liao, M.; Fillery, I.R.; Palta, J.A. Early vigorous growth is a major factor influencing nitrogen uptake in wheat. Funct. Plant Biol. 2004, 31, 121–129. [Google Scholar] [CrossRef]
  3. Perveen, S.; Nigar, S.; Khalil, S.; Zubair, M. Vigor tests used to rank seed lot quality and predict field emergence in wheat. Pak. J. Bot. 2010, 42, 3147–3155. [Google Scholar]
  4. Sissons, M. Role of durum wheat composition on the quality of pasta and bread. Food 2008, 2, 75–90. [Google Scholar]
  5. Li, Y.-F.; Wu, Y.; Hernandez-Espinosa, N.; Peña, R.J. Heat and drought stress on durum wheat: Responses of genotypes, yield, and quality parameters. J. Cereal Sci. 2013, 57, 398–404. [Google Scholar] [CrossRef]
  6. Liu, H.; Bruce, D.R.; Sissons, M.; Able, A.J.; Able, J.A. Genotype-dependent changes in the phenolic content of durum under water-deficit stress. Cereal Chem. 2018, 95, 59–78. [Google Scholar] [CrossRef]
  7. Guzmán, C.; Autrique, J.E.; Mondal, S.; Singh, R.P.; Govindan, V.; Morales-Dorantes, A.; Posadas-Romano, G.; Crossa, J.; Ammar, K.; Peña, R.J. Response to drought and heat stress on wheat quality, with special emphasis on bread-making quality, in durum wheat. Field Crop. Res. 2016, 186, 157–165. [Google Scholar] [CrossRef]
  8. Liu, H.; Able, A.J.; Able, J.A. Genotypic performance of Australian durum under single and combined water-deficit and heat stress during reproduction. Sci. Rep. 2019, 9, 1–17. [Google Scholar] [CrossRef] [Green Version]
  9. Liu, H.; Searle, I.R.; Mather, D.E.; Able, A.J.; Able, J.A. Morphological, physiological and yield responses of durum wheat to pre-anthesis water-deficit stress are genotype-dependent. Crop Pasture Sci. 2015, 66, 1024–1038. [Google Scholar] [CrossRef]
  10. Hatzig, S.V.; Nuppenau, J.-N.; Snowdon, R.J.; Schießl, S.V. Drought stress has transgenerational effects on seeds and seedlings in winter oilseed rape (Brassica napus L.). BMC Plant Biol. 2018, 18, 297. [Google Scholar] [CrossRef]
  11. Wijewardana, C.; Reddy, K.R.; Krutz, L.J.; Gao, W.; Bellaloui, N. Drought stress has transgenerational effects on soybean seed germination and seedling vigor. PLoS ONE 2019, 14, e0214977. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Begcy, K.; Sandhu, J.; Walia, H. Transient heat stress during early seed development primes germination and seedling establishment in rice. Front. Plant Sci. 2018, 9, 1768. [Google Scholar] [CrossRef]
  13. Carvalho, M.E.A.; Piotto, F.A.; Nogueira, M.L.; Gomes-Junior, F.G.; Chamma, H.M.C.P.; Pizzaia, D.; Azevedo, R.A. Cadmium exposure triggers genotype-dependent changes in seed vigor and germination of tomato offspring. Protoplasma 2018, 255, 989–999. [Google Scholar] [CrossRef] [PubMed]
  14. Wollmann, H.; Berger, F. Epigenetic reprogramming during plant reproduction and seed development. Curr. Opin. Plant Biol. 2012, 15, 63–69. [Google Scholar] [CrossRef] [PubMed]
  15. Lämke, J.; Bäurle, I. Epigenetic and chromatin-based mechanisms in environmental stress adaptation and stress memory in plants. Genome Biol. 2017, 18, 124. [Google Scholar] [CrossRef]
  16. Holeski, L.M.; Jander, G.; Agrawal, A.A. Transgenerational defense induction and epigenetic inheritance in plants. Trends Ecol. Evol. 2012, 27, 618–626. [Google Scholar] [CrossRef]
  17. Liu, H.; Able, A.J.; Able, J.A. SMARTER de-stressed cereal breeding. Trends Plant Sci. 2016, 21, 909–925. [Google Scholar] [CrossRef]
  18. Budak, H.; Akpinar, B.A. Plant miRNAs: Biogenesis, organization and origins. Funct. Integr. Genom. 2015, 15, 523–531. [Google Scholar] [CrossRef]
  19. Budak, H.; Kantar, M.; Bulut, R.; Akpinar, B.A. Stress responsive miRNAs and isomiRs in cereals. Plant Sci. 2015, 235, 1–13. [Google Scholar] [CrossRef]
  20. Giusti, L.; Mica, E.; Bertolini, E.; De Leonardis, A.M.; Faccioli, P.; Cattivelli, L.; Crosatti, C. microRNAs differentially modulated in response to heat and drought stress in durum wheat cultivars with contrasting water use efficiency. Funct. Integr. Genom. 2017, 17, 293–309. [Google Scholar] [CrossRef]
  21. Zuluaga, D.L.; De Paola, D.; Janni, M.; Curci, P.L.; Sonnante, G. Durum wheat miRNAs in response to nitrogen starvation at the grain filling stage. PLoS ONE 2017, 12, e0183253. [Google Scholar] [CrossRef]
  22. Liu, H.; Searle, I.R.; Watson-Haigh, N.S.; Baumann, U.; Mather, D.E.; Able, A.J.; Able, J.A. Genome-wide identification of microRNAs in leaves and the developing head of four durum genotypes during water deficit stress. PLoS ONE 2015, 10, e0142799. [Google Scholar] [CrossRef] [PubMed]
  23. Liu, H.; Able, A.J.; Able, J.A. Genotypic water-deficit stress responses in durum wheat: Association between physiological traits, microRNA regulatory modules and yield components. Funct. Plant Biol. 2017, 44, 538–551. [Google Scholar] [CrossRef] [Green Version]
  24. Liu, H.; Able, A.J.; Able, J.A. Water-deficit stress responsive microRNAs and their targets in four durum wheat genotypes. Funct. Integr. Genom. 2017, 17, 237–251. [Google Scholar] [CrossRef] [PubMed]
  25. Wang, Y.; Li, K.; Chen, L.; Zou, Y.; Liu, H.; Tian, Y.; Li, D.; Wang, R.; Zhao, F.; Ferguson, B.J.; et al. microRNA167-directed regulation of the auxin response factors, GmARF8a and GmARF8b, is required for soybean nodulation and lateral root development. Plant Physiol. 2015, 168, 984–999. [Google Scholar] [CrossRef] [Green Version]
  26. Zuluaga, D.L.; Sonnante, G. The use of nitrogen and its regulation in cereals: Structural genes, transcription factors, and the role of miRNAs. Plants 2019, 8, 294. [Google Scholar] [CrossRef] [Green Version]
  27. Ragupathy, R.; Ravichandran, S.; Mahdi, M.S.R.; Huang, D.; Reimer, E.; Domaratzki, M.; Cloutier, S. Deep sequencing of wheat sRNA transcriptome reveals distinct temporal expression pattern of miRNAs in response to heat, light and UV. Sci. Rep. 2016, 6, 39373. [Google Scholar] [CrossRef] [Green Version]
  28. Ravichandran, S.; Ragupathy, R.; Edwards, T.; Domaratzki, M.; Cloutier, S. MicroRNA-guided regulation of heat stress response in wheat. BMC Genom. 2019, 20, 488. [Google Scholar] [CrossRef] [Green Version]
  29. Yakovlev, I.A.; Fossdal, C.G. In silico analysis of small RNAs suggest roles for novel and conserved miRNAs in the formation of epigenetic memory in somatic embryos of Norway spruce. Front Physiol. 2017, 8, 674. [Google Scholar] [CrossRef]
  30. Zuluaga, D.L.; Liuzzi, V.; Curci, P.L.; Sonnante, G. MicroRNAs in durum wheat seedlings under chronic and short-term nitrogen stress. Funct. Integr. Genom. 2018, 18, 645–657. [Google Scholar] [CrossRef]
  31. Kantar, M.; Lucas, S.J.; Budak, H. miRNA expression patterns of Triticum dicoccoides in response to shock drought stress. Planta 2011, 233, 471–484. [Google Scholar] [CrossRef] [PubMed]
  32. Akpinar, B.A.; Kantar, M.; Budak, H. Root precursors of microRNAs in wild emmer and modern wheats show major differences in response to drought stress. Funct. Integr. Genom. 2015, 15, 587–598. [Google Scholar] [CrossRef] [PubMed]
  33. Maccaferri, M.; Harris, N.S.; Twardziok, S.O.; Pasam, R.K.; Gundlach, H.; Spannagl, M.; Ormanbekova, D.; Lux, T.; Prade, V.M.; Milner, S.G. Durum wheat genome highlights past domestication signatures and future improvement targets. Nat. Genet. 2019, 51, 885–895. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. De Paola, D.; Zuluaga, D.L.; Sonnante, G. The miRNAome of durum wheat: Isolation and characterisation of conserved and novel microRNAs and their target genes. BMC Genom. 2016, 17, 505. [Google Scholar] [CrossRef] [Green Version]
  35. Jeyaraj, A.; Zhang, X.; Hou, Y.; Shangguan, M.; Gajjeraman, P.; Li, Y.; Wei, C. Genome-wide identification of conserved and novel microRNAs in one bud and two tender leaves of tea plant (Camellia sinensis) by small RNA sequencing, microarray-based hybridization and genome survey scaffold sequences. BMC Plant Biol. 2017, 17, 212. [Google Scholar] [CrossRef] [Green Version]
  36. Changhai, S.; Baodi, D.; Yunzhou, Q.; Yuxin, L.; Lei, S.; Mengyu, L.; Haipei, L. Physiological regulation of high transpiration efficiency in winter wheat under drought conditions. Plant Soil Environ. 2010, 56, 340–347. [Google Scholar] [CrossRef] [Green Version]
  37. Zheng, X.; Liu, H.; Ji, H.; Wang, Y.; Dong, B.; Qiao, Y.; Liu, M.; Li, X. The wheat GT factor TaGT2L1D negatively regulates drought tolerance and plant development. Sci. Rep. 2016, 6, 27042. [Google Scholar] [CrossRef] [Green Version]
  38. Dong, B.; Zheng, X.; Liu, H.; Able, J.A.; Yang, H.; Zhao, H.; Zhang, M.; Qiao, Y.; Wang, Y.; Liu, M. Effects of drought stress on pollen sterility, grain yield, abscisic acid and protective enzymes in two winter wheat cultivars. Front. Plant Sci. 2017, 8, 1008. [Google Scholar] [CrossRef]
  39. Dong, B.; Yang, H.; Liu, H.; Qiao, Y.; Zhang, M.; Wang, Y.; Xie, Z.; Liu, M. Effects of shading stress on grain number, yield, and photosynthesis during early reproductive growth in wheat. Crop Sci. 2019, 59, 363–378. [Google Scholar] [CrossRef]
  40. Barnabás, B.; Jäger, K.; Fehér, A. The effect of drought and heat stress on reproductive processes in cereals. Plant Cell Environ. 2008, 31, 11–38. [Google Scholar] [CrossRef]
  41. Grass, L.; Burris, J. Effect of heat stress during seed development and maturation on wheat (Triticum durum) seed quality. I. Seed germination and seedling vigor. Can. J. Plant Sci. 1995, 75, 821–829. [Google Scholar] [CrossRef] [Green Version]
  42. Alptekin, B.; Langridge, P.; Budak, H. Abiotic stress miRNomes in the Triticeae. Funct. Integr. Genom. 2016, 117, 145–170. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Kozomara, A.; Birgaoanu, M.; Griffiths-Jones, S. miRBase: From microRNA sequences to function. Nucleic Acids Res. 2018, 47, D155–D162. [Google Scholar] [CrossRef] [PubMed]
  44. Kwok, R.P.; Lundblad, J.R.; Chrivia, J.C.; Richards, J.P.; Bächinger, H.P.; Brennan, R.G.; Roberts, S.G.; Green, M.R.; Goodman, R.H. Nuclear protein CBP is a coactivator for the transcription factor CREB. Nature 1994, 370, 223. [Google Scholar] [CrossRef]
  45. Gehring, C. Adenyl cyclases and cAMP in plant signaling-past and present. Cell Commum Signal. 2010, 8, 15. [Google Scholar] [CrossRef] [Green Version]
  46. Colcombet, J.; Krysan, P.J. Cellular complexity in MAPK signaling in plants: Questions and emerging tools to answer them. Front. Plant Sci. 2018, 9, 1674. [Google Scholar]
  47. Smékalová, V.; Doskočilová, A.; Komis, G.; Šamaj, J. Crosstalk between secondary messengers, hormones and MAPK modules during abiotic stress signalling in plants. Biotechnol. Adv. 2014, 32, 2–11. [Google Scholar] [CrossRef]
  48. Wang, X.; Xin, C.; Cai, J.; Zhou, Q.; Dai, T.; Cao, W.; Jiang, D. Heat priming induces trans-generational tolerance to high temperature stress in wheat. Front. Plant Sci. 2016, 7, 501. [Google Scholar] [CrossRef] [Green Version]
  49. Šamaj, J.; Ovecka, M.; Hlavacka, A.; Lecourieux, F.; Meskiene, I.; Lichtscheidl, I.; Lenart, P.; Salaj, J.; Volkmann, D.; Bögre, L. Involvement of the mitogen-activated protein kinase SIMK in regulation of root hair tip growth. EMBO J. 2002, 21, 3296–3306. [Google Scholar] [CrossRef] [Green Version]
  50. Chen, W.; Moore, M.J. Spliceosomes. Curr. Biol. 2015, 25, R181–R183. [Google Scholar] [CrossRef] [Green Version]
  51. Ding, F.; Cui, P.; Wang, Z.; Zhang, S.; Ali, S.; Xiong, L. Genome-wide analysis of alternative splicing of pre-mRNA under salt stress in Arabidopsis. BMC Genom. 2014, 15, 431. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Ling, Y.; Alshareef, S.; Butt, H.; Lozano-Juste, J.; Li, L.; Galal, A.A.; Moustafa, A.; Momin, A.A.; Tashkandi, M.; Richardson, D.N. Pre-mRNA splicing repression triggers abiotic stress signaling in plants. Plant J. 2017, 89, 291–309. [Google Scholar] [CrossRef] [PubMed]
  53. Laloum, T.; Martín, G.; Duque, P. Alternative splicing control of abiotic stress responses. Trends Plant Sci. 2018, 23, 140–150. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Yan, K.; Liu, P.; Wu, C.-A.; Yang, G.-D.; Xu, R.; Guo, Q.-H.; Huang, J.-G.; Zheng, C.-C. Stress-induced alternative splicing provides a mechanism for the regulation of microRNA processing in Arabidopsis thaliana. Mol. Cell 2012, 48, 521–531. [Google Scholar] [CrossRef] [Green Version]
  55. Kruszka, K.; Pacak, A.; Swida-Barteczka, A.; Nuc, P.; Alaba, S.; Wroblewska, Z.; Karlowski, W.; Jarmolowski, A.; Szweykowska-Kulinska, Z. Transcriptionally and post-transcriptionally regulated microRNAs in heat stress response in barley. J. Exp. Bot. 2014, 65, 6123–6135. [Google Scholar] [CrossRef] [Green Version]
  56. Meng, Y.; Qu, G.; Wang, T.; Sun, Q.; Liang, D.; Hu, S. Enhancement of germination and seedling growth of wheat seed using dielectric barrier discharge plasma with various gas sources. Plasma Chem. Plasma Process. 2017, 37, 1105–1119. [Google Scholar] [CrossRef]
  57. Zahoranová, A.; Henselová, M.; Hudecová, D.; Kaliňáková, B.; Kováčik, D.; Medvecká, V.; Černák, M. Effect of cold atmospheric pressure plasma on the wheat seedlings vigor and on the inactivation of microorganisms on the seeds surface. Plasma Chem. Plasma Process. 2016, 36, 397–414. [Google Scholar] [CrossRef]
  58. Feng, H.; Zhang, Q.; Wang, B.; Fu, Y.; Huang, L.; Wang, X.; Kang, Z. Exploration of microRNAs and their targets engaging in the resistance interaction between wheat and stripe rust. Frontiers in Plant Science 2015, 6, 469. [Google Scholar] [CrossRef] [Green Version]
  59. Li, X.; Shahid, M.Q.; Wu, J.; Wang, L.; Liu, X.; Lu, Y. Comparative small RNA analysis of pollen development in autotetraploid and diploid rice. Int. J. Mol. 2016, 17, 499. [Google Scholar] [CrossRef]
  60. Li, W.; Jia, Y.; Liu, F.; Wang, F.; Fan, F.; Wang, J.; Zhu, J.; Xu, Y.; Zhong, W.; Yang, J. Integration analysis of small RNA and degradome sequencing reveals microRNAs responsive to Dickeya zeae in resistant rice. Int. J. Mol. 2019, 20, 222. [Google Scholar] [CrossRef] [Green Version]
Plants 09 00189 g001 550
Figure 1. Small RNA size distribution of four durum wheat libraries. TCG, constructed from DBA Aurora CG seedlings. TWH, constructed from DBA Aurora WH seedlings. SCG, constructed from L6 CG seedlings. SWH, constructed from L6 WH seedlings.
Figure 1. Small RNA size distribution of four durum wheat libraries. TCG, constructed from DBA Aurora CG seedlings. TWH, constructed from DBA Aurora WH seedlings. SCG, constructed from L6 CG seedlings. SWH, constructed from L6 WH seedlings.
Plants 09 00189 g001
Plants 09 00189 g002 550
Figure 2. Venn diagram of all mature miRNAs in four libraries. TCG, constructed from DBA Aurora CG seedlings. TWH, constructed from DBA Aurora WH seedlings. SCG, constructed from L6 CG seedlings. SWH, constructed from L6 WH seedlings.
Figure 2. Venn diagram of all mature miRNAs in four libraries. TCG, constructed from DBA Aurora CG seedlings. TWH, constructed from DBA Aurora WH seedlings. SCG, constructed from L6 CG seedlings. SWH, constructed from L6 WH seedlings.
Plants 09 00189 g002
Plants 09 00189 g003 550
Figure 3. Ranking of reference plant species based on the number of miRNAs conserved. gma, Glycine max. osa, Oryza sativa. tae, Triticum aestivum. bdi, Brachypodium distachyon. ata, Aegilops tauschii. zma, Zea mays. mdm, Malus domestica. sbi, Sorghum bicolor. mes, Manihot esculenta. ptc, Populus trichocarpa. lus, Linum usitatissimum. cme, Cucumis melo. vvi, Vitis vinifera. ath, Arabidopsis thaliana. aly, Arabidopsis lyrata. csi, Citrus sinensis. mtr, Medicago truncatula. cas, Camelina sativa. ppe, Prunus persica. bna, Brassica napus. nta, Nicotiana tabacum. rco, Ricinus communis. hvu, Hordeum vulgare. sly, Solanum lycopersicum. cpa, Carica papaya. stu, Solanum tuberosum. bra, Brassica rapa. ghr, Gossypium hirsutum. sof, Saccharum officinarum. lja, Lotus japonicas. ssp, Saccharum ssp. aqc, Aquilegia caerulea. gra, Gossypium raimondii. far, Festuca arundinacea. hbr, Hevea brasiliensis. cca, Cynara cardunculus. rgl, Rehmannia glutinosa. Ttu, Triticum turgidum.
Figure 3. Ranking of reference plant species based on the number of miRNAs conserved. gma, Glycine max. osa, Oryza sativa. tae, Triticum aestivum. bdi, Brachypodium distachyon. ata, Aegilops tauschii. zma, Zea mays. mdm, Malus domestica. sbi, Sorghum bicolor. mes, Manihot esculenta. ptc, Populus trichocarpa. lus, Linum usitatissimum. cme, Cucumis melo. vvi, Vitis vinifera. ath, Arabidopsis thaliana. aly, Arabidopsis lyrata. csi, Citrus sinensis. mtr, Medicago truncatula. cas, Camelina sativa. ppe, Prunus persica. bna, Brassica napus. nta, Nicotiana tabacum. rco, Ricinus communis. hvu, Hordeum vulgare. sly, Solanum lycopersicum. cpa, Carica papaya. stu, Solanum tuberosum. bra, Brassica rapa. ghr, Gossypium hirsutum. sof, Saccharum officinarum. lja, Lotus japonicas. ssp, Saccharum ssp. aqc, Aquilegia caerulea. gra, Gossypium raimondii. far, Festuca arundinacea. hbr, Hevea brasiliensis. cca, Cynara cardunculus. rgl, Rehmannia glutinosa. Ttu, Triticum turgidum.
Plants 09 00189 g003
Plants 09 00189 g004 550
Figure 4. Number of up-regulated and down-regulated miRNAs in comparisons made between libraries at different p values (p < 0.01, p < 0.05 and p < 0.1). TCG, constructed from DBA Aurora CG seedlings. TWH, constructed from DBA Aurora WH seedlings. SCG, constructed from L6 CG seedlings. SWH, constructed from L6 WH seedlings.
Figure 4. Number of up-regulated and down-regulated miRNAs in comparisons made between libraries at different p values (p < 0.01, p < 0.05 and p < 0.1). TCG, constructed from DBA Aurora CG seedlings. TWH, constructed from DBA Aurora WH seedlings. SCG, constructed from L6 CG seedlings. SWH, constructed from L6 WH seedlings.
Plants 09 00189 g004
Plants 09 00189 g005 550
Figure 5. The most represented Gene Ontology terms among the predicted targets in three categories. GO terms were shown for each category based on the number of genes matched for each GO term (from higher to lower). Percent of genes (%) indicates the percentage of the gene number matched for each GO among the total number of genes within each category.
Figure 5. The most represented Gene Ontology terms among the predicted targets in three categories. GO terms were shown for each category based on the number of genes matched for each GO term (from higher to lower). Percent of genes (%) indicates the percentage of the gene number matched for each GO among the total number of genes within each category.
Plants 09 00189 g005
Plants 09 00189 g006 550
Figure 6. KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis of predicted target genes. p-value (red to blue colour) represents the significance of the matched gene ratio (rich factor). Dot size is proportionate to the number of target genes matched to the pathway term.
Figure 6. KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis of predicted target genes. p-value (red to blue colour) represents the significance of the matched gene ratio (rich factor). Dot size is proportionate to the number of target genes matched to the pathway term.
Plants 09 00189 g006
Table
Table 1. Germination potential (Gp, %), germination rate (Gr, %), germination index (Gi), coleoptile length, and young root length of the four seed groups. AuCG, seeds from DBA Aurora control group. AuWH, seeds from DBA Aurora stress group. L6CG, seeds from L6 control group. L6WH, seeds from L6 stress group.
Table 1. Germination potential (Gp, %), germination rate (Gr, %), germination index (Gi), coleoptile length, and young root length of the four seed groups. AuCG, seeds from DBA Aurora control group. AuWH, seeds from DBA Aurora stress group. L6CG, seeds from L6 control group. L6WH, seeds from L6 stress group.
Seed SourceGp, %Gr, %GiColeoptile Length (mm)Young Root Length (mm)
AuCG32.2% ± 1.1%98.9% ± 1.1%35.1 ± 0.831.88 ± 0.5248.56 ± 1.06
AuWH26.7% ± 1.9%97.8% ± 1.1%32.2 ± 0.931.06 ± 0.6750.44 ± 1.19
L6CG17.8% ± 1.1%93.3% ± 1.9%27.6 ± 0.230.38 ± 0.7538.88 ± 1.28
L6WH14.4% ± 1.1%92.2% ± 1.1%24.4 ± 0.630.50 ± 0.7533.75 ± 1.23
F pr.
Genotype		<0.0010.004<0.0010.134<0.001
F pr.
Parent treatment		0.0110.4380.0020.6140.178
F pr.
Genotype × Parent treatment		0.4381.0000.8420.4930.005
l.s.d
Genotype		3.1%3.1%1.6n.a2.39
l.s.d
Parent treatment		3.1%n.a1.6n.an.a
l.s.d
Genotype × Parent treatment		n.an.an.an.a3.38
Table
Table 2. Root traits, shoot traits and seedling vigour indices of the four seed groups. AuCG, seeds from DBA Aurora control group. AuWH, seeds from DBA Aurora stress group. L6CG, seeds from L6 control group. L6WH, seeds from L6 stress group.
Table 2. Root traits, shoot traits and seedling vigour indices of the four seed groups. AuCG, seeds from DBA Aurora control group. AuWH, seeds from DBA Aurora stress group. L6CG, seeds from L6 control group. L6WH, seeds from L6 stress group.
Seed SourceRoot Fresh Weight (g)Root Dry Weight (g)Root Length (cm)Shoot Fresh Weight (g)Shoot Dry Weight (g)Seedling Height (cm)Seedling Vigour Index ISeedling Vigour Index II
AuCG0.718 ± 0.0160.0807 ± 0.001625.86 ± 0.180.878 ± 0.0250.1244 ± 0.003027.33 ± 0.261577.28 ± 30.48202.80 ± 2.42
AuWH0.736 ± 0.0140.0812 ± 0.002228.46 ± 0.220.836 ± 0.0220.1189 ± 0.002326.36 ± 0.231537.56 ± 28.01195.63 ± 3.70
L6CG0.710 ± 0.0130.0721 ± 0.001725.10 ± 0.190.844 ± 0.0200.1141 ± 0.002526.98 ± 0.241450.17 ± 21.74173.78 ± 2.92
L6WH0.658 ± 0.0140.0645 ± 0.001426.56 ± 0.210.778 ± 0.0160.1116 ± 0.003325.39 ± 0.311323.39 ± 16.81162.32 ± 3.79
F pr.
Genotype		0.005<0.001<0.0010.0340.0040.018<0.001<0.001
F pr.
Parent treatment		0.2420.055<0.0010.0150.160<0.0010.0020.008
F pr.
Genotype × Parent treatment		0.0180.0300.0080.5510.6080.2450.0910.516
l.s.d
Genotype		0.0290.00360.410.0430.00570.5450.906.67
l.s.d
Parent treatment		n.an.a0.410.043n.a0.5450.906.67
l.s.d
Genotype × Parent treatment		0.0410.00510.58n.an.a0.76n.an.a
Table
Table 3. Overview of sequencing read counts in four sRNA libraries. TCG, constructed from DBA Aurora CG seedlings. TWH, constructed from DBA Aurora WH seedlings. SCG, constructed from L6 CG seedlings. SWH, constructed from L6 WH seedlings. Clean reads may not be equal to raw reads minus all the filter type reads, as there could be overlapped sequences between mRNA, Rfam and Repeats filters.
Table 3. Overview of sequencing read counts in four sRNA libraries. TCG, constructed from DBA Aurora CG seedlings. TWH, constructed from DBA Aurora WH seedlings. SCG, constructed from L6 CG seedlings. SWH, constructed from L6 WH seedlings. Clean reads may not be equal to raw reads minus all the filter type reads, as there could be overlapped sequences between mRNA, Rfam and Repeats filters.
Read TypeTCGTWHSCGSWH
Total ReadsUnique ReadsTotal ReadsUnique ReadsTotal ReadsUnique ReadsTotal ReadsUnique Reads
Raw reads20,042,8338,293,42520,900,1868,593,35724,262,3438,865,72418,786,4687,221,909
3ADT & Length filter6,000,3641,418,8935,651,8351,475,2357,669,0101,832,9616,589,9331,676,502
Junk reads92,53868,473101,06673,192111,33979,08876,37156,322
Rfam - rRNA215,7165440228,7155378308,3126809266,0897444
Rfam - tRNA230,8952119215,3402207248,0912853310,9373964
Rfam - snoRNA9600393981339610,8863638009304
Rfam - snRNA2364181258119432502643580215
Other Rfam RNA13,68073314,27670619,44588617,165893
mRNA878,03923,905915,46825,1701,090,07626,584798,62920,880
Repeats6104186607319393122657704243
Clean reads12,631,6566,773,72213,795,2487,011,28814,830,8726,916,49310,781,6725,456,127
Table
Table 4. Categorisation of all MIR-miRNA entries identified in this study. G1-G4 represent conserved miRNAs while G5 represents novel miRNAs. Unique miRNAs represent distinct mature miRNA products. G1, pre-miRNAs can be mapped to the durum genome and ESTs (expressed sequence tags). G2, sequencing reads mapped to the durum genome and ESTs with their extended genome sequences forming secondary hairpins. G3, sequencing reads mapped to the durum genome and ESTs, but their extended genome sequences do not form secondary hairpins. G4, neither aligned reference pre-miRNAs nor sequencing reads can be further mapped to the genome. G5, novel miRNAs where reads cannot be mapped to the miRBase but can be mapped to the genome, and secondary hairpins can be formed from extended genome locations. TCG, constructed from DBA Aurora CG seedlings. TWH, constructed from DBA Aurora WH seedlings. SCG, constructed from L6 CG seedlings. SWH, constructed from L6 WH seedlings.
Table 4. Categorisation of all MIR-miRNA entries identified in this study. G1-G4 represent conserved miRNAs while G5 represents novel miRNAs. Unique miRNAs represent distinct mature miRNA products. G1, pre-miRNAs can be mapped to the durum genome and ESTs (expressed sequence tags). G2, sequencing reads mapped to the durum genome and ESTs with their extended genome sequences forming secondary hairpins. G3, sequencing reads mapped to the durum genome and ESTs, but their extended genome sequences do not form secondary hairpins. G4, neither aligned reference pre-miRNAs nor sequencing reads can be further mapped to the genome. G5, novel miRNAs where reads cannot be mapped to the miRBase but can be mapped to the genome, and secondary hairpins can be formed from extended genome locations. TCG, constructed from DBA Aurora CG seedlings. TWH, constructed from DBA Aurora WH seedlings. SCG, constructed from L6 CG seedlings. SWH, constructed from L6 WH seedlings.
GroupTotalTCGTWHSCGSWH
Pre-miRNAUnique miRNAPre-miRNAUnique miRNAPre-miRNAUnique miRNAPre-miRNAUnique miRNAPre-miRNAUnique miRNA
G1161270158252158256158259158257
G2446483368392372405382410300322
G3838890663694701736722763611644
G448493434313342433434
G5437485426446429445429455421440
Table
Table 5. Highly expressed differentially expressed miRNAs (DEMs) with a |log2 (fold change)| > 1 in each comparison made between two libraries (TWH vs TCG, SWH vs SCG, SCG vs TCG, SWH vs TWH). TCG, constructed from DBA Aurora CG seedlings. TWH, constructed from DBA Aurora WH seedlings. SCG, constructed from L6 CG seedlings. SWH, constructed from L6 WH seedlings.
Table 5. Highly expressed differentially expressed miRNAs (DEMs) with a |log2 (fold change)| > 1 in each comparison made between two libraries (TWH vs TCG, SWH vs SCG, SCG vs TCG, SWH vs TWH). TCG, constructed from DBA Aurora CG seedlings. TWH, constructed from DBA Aurora WH seedlings. SCG, constructed from L6 CG seedlings. SWH, constructed from L6 WH seedlings.
miR_namemiR_sequp/downlog2(FC)TCGTWH
osa-miR168a-3p_L-3CCCGCCTTGCACCAAGTGAATup2.04357114,660
ata-miR169a-3pTGGGCAAGTCACCCTGGCTACCup1.115331154
miR_namemiR_sequp/downlog2(FC)SCGSWH
tae-miR9661-5p_1ss14GATGAAGTAGAGCAGAGACCTCAdown−1.0635131689
tae-miR9774_L+2AACAAGATATTGGGTATTTCTGTCup1.2715893841
miR_namemiR_sequp/downlog2(FC)TCGSCG
ata-miR167b-5pTGAAGCTGCCAGCATGATCTGAdown−1.2120,0898687
tae-MIR7757-p3_1ss19AGTGGATAGTTTGAGGTTTTGTTTdown−1.0465243162
ata-miR167b-3pAGGTCATGCTGGAGTTTCATCdown−1.1157462653
ttu-miR160TGCCTGGCTCCCTGTATGCCAdown−1.0447392304
mtr-miR398a-5p_1ss7ACGGAGTGCCACTGAGAACACAAGdown−1.362354917
csi-miR166c-3p_R+1TCGGACCAGGCTTCATTCCCAdown−1.521062370
tae-miR9654b-3p_1ss19GATTCCGAAAGGCTTGAAGCAAATdown−1.101542721
ata-miR396e-3pGTTCAATAAAGCTGTGGGAAAdown−1.10737345
sbi-miR396c_R+1TTCCACAGCTTTCTTGAACTTTdown−1.07758361
PC-5p-1206_1154TTTGCGCAGAAGGGAGAAATCupinf02228
tae-MIR5048-p3_1ss18TGAATATATTTGCAGGTTTGAGGup1.0846779898
ata-miR169f-3p_R-1GGCAAGTCCGTCCTTGGCTACup1.22460910,772
osa-miR168a-3p_L-3CCCGCCTTGCACCAAGTGAATup1.2635718570
ata-MIR169d-p3TTGTCCTTGGCTACACCTAGTup1.5017424936
tae-miR9661-5p_1ss14GATGAAGTAGAGCAGAGACCTCAup3.493123513
ata-MIR396d-p3AAGCCCATGGAAACCATGCCCup4.6629735
ata-miR169i-5p_1ss15TCTAGCCAAGGATGACCTGCCTGup3.3953553
ata-miR169a-3pTGGGCAAGTCACCCTGGCTACCup1.495331494
PC-5p-726_2034TACGGCAAAGCCGTCGGCATAup2.03194793
ata-miR169h-3p_L-1CAAGTTGTTCTTGGCTAGCup1.254841154
tae-MIR164-p3CATGTGCCTTTCTTCTCCACCup1.02419848
miR_namemiR_sequp/downlog2(FC)TWHSWH
mtr-miR398a-5p_1ss7ACGGAGTGCCACTGAGAACACAAGdown−1.691516471
tae-miR9659-3pTCCAATGGTTGTTCACGGCATCdown−1.17910404
osa-miR160f-5p_L-2R+2CCTGGCTCCCTGAATGCCATCdown−1.14817370
sbi-miR396c_R+1TTCCACAGCTTTCTTGAACTTTdown−1.12609280
PC-5p-1206_1154TTTGCGCAGAAGGGAGAAATCupinf02473
ata-miR169f-3p_R-1GGCAAGTCCGTCCTTGGCTACup1.1441039018
tae-miR9774_L+2AACAAGATATTGGGTATTTCTGTCup1.5712963841
tae-miR9661-5p_1ss14GATGAAGTAGAGCAGAGACCTCAup2.552881689
ata-MIR396d-p3AAGCCCATGGAAACCATGCCCup4.9723710
ata-miR169i-5p_1ss15TCTAGCCAAGGATGACCTGCCTGup2.5282469

Share and Cite

MDPI and ACS Style

Liu, H.; Able, A.J.; Able, J.A. Transgenerational Effects of Water-Deficit and Heat Stress on Germination and Seedling Vigour—New Insights from Durum Wheat microRNAs. Plants 2020, 9, 189. https://doi.org/10.3390/plants9020189

AMA Style

Liu H, Able AJ, Able JA. Transgenerational Effects of Water-Deficit and Heat Stress on Germination and Seedling Vigour—New Insights from Durum Wheat microRNAs. Plants. 2020; 9(2):189. https://doi.org/10.3390/plants9020189

Chicago/Turabian Style

Liu, Haipei, Amanda J. Able, and Jason A. Able. 2020. "Transgenerational Effects of Water-Deficit and Heat Stress on Germination and Seedling Vigour—New Insights from Durum Wheat microRNAs" Plants 9, no. 2: 189. https://doi.org/10.3390/plants9020189

APA Style

Liu, H., Able, A. J., & Able, J. A. (2020). Transgenerational Effects of Water-Deficit and Heat Stress on Germination and Seedling Vigour—New Insights from Durum Wheat microRNAs. Plants, 9(2), 189. https://doi.org/10.3390/plants9020189

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