Diverse set of microRNAs are responsive to powdery mildew infection and heat stress in wheat (Triticum aestivum L.)
- Mingming Xin†1, 2,
- Yu Wang†1, 2,
- Yingyin Yao†1, 2,
- Chaojie Xie1, 2, 3Email author,
- Huiru Peng1, 2,
- Zhongfu Ni1, 2 and
- Qixin Sun1, 2, 3Email author
© Xin et al; licensee BioMed Central Ltd. 2010
Received: 10 July 2009
Accepted: 24 June 2010
Published: 24 June 2010
MicroRNAs (miRNAs) are a class of small non-coding regulatory RNAs that regulate gene expression by guiding target mRNA cleavage or translational inhibition. MiRNAs can have large-scale regulatory effects on development and stress response in plants.
To test whether miRNAs play roles in regulating response to powdery mildew infection and heat stress in wheat, by using Solexa high-throughput sequencing we cloned the small RNA from wheat leaves infected by preponderant physiological strain Erysiphe graminis f. sp. tritici (Egt) or by heat stress treatment. A total of 153 miRNAs were identified, which belong to 51 known and 81 novel miRNA families. We found that 24 and 12 miRNAs were responsive to powdery mildew infection and heat stress, respectively. We further predicted that 149 target genes were potentially regulated by the novel wheat miRNA.
Our results indicated that diverse set of wheat miRNAs were responsive to powdery mildew infection and heat stress and could function in wheat responses to both biotic and abiotic stresses.
Wheat (Triticum aestivum, AABBDD, 2n = 42) is the most widely grown crop plant, occupying 17% of all the cultivated land, provide approximately 55% of carbohydrates for world human consumption , Biotic and abiotic stresses are important limiting factors for yield and grain quality in wheat production. For instance, powdery mildew, caused by the obligate biotrophic fungus Blumeria graminis f. sp. tritici (Bgt), is one of the most devastating diseases of wheat in China and worldwide, and causing significant yield losses [2, 3]. High temperature, often combined with drought stress, causes yield loss and reduces the grain quality . To reduce the damage caused by biotic and abiotic stresses, plants have evolved sophisticated adaptive response mechanisms to reprogram gene expression at the transcriptional, post-transcriptional and post-translational levels . Therefore, transcript profiling have been successfully employed to determine the transcriptional responses to powdery mildew infection and heat stress in wheat, and the results revealed that a number of genes were significantly induced or repressed in response to these stresses [6, 7]. In addition, much progress has been made in unraveling the complex post-transcriptional regulation mechanisms in response to stress. Recently discovered microRNAs (miRNAs) and endogenous small interfering RNAs (siRNAs) have emerged as important posttranscriptional regulators in plant stress responses . However, no report has been published on the roles of small RNAs in wheat response to abiotic and biotic stresses.
MicroRNAs (miRNAs) are a class of small RNAs that serve as posttranscriptional negative regulators of gene expression in plants and animals [9–14]. Up to date, plant miRNAs have been shown to function in many plant processes, including developmental transitions [15, 16], leaf growth , organ polarity, auxin signaling  and RNA metabolism [20–22]. Importantly, increasing evidence indicated that miRNAs also play important roles in plant response to abiotic and biotic stresses . The observation that some of the miRNAs are up- or down- regulated in response to stress implies that these miRNAs could play important roles in stress tolerance . For example, in Arabidopsis, miR393 and other miRNAs are induced by cold stress [22, 23]. In rice, miR169g and miR393 are up-regulated under drought stress . Arabidopsis miR398 directs the cleavage of CSD1 and CSD2 mRNA under normal conditions, and down-regulation of miR398 by oxidative stress results in accumulation of CDS1 and CSD2 mRNAs . Recent studies in Arabidopsis have also established that miR399, miR395 and miR398 are induced in response to phosphate-, sulfate- and Cu2+-deprived conditions, respectively [8, 26–29]. In addition, among the 42 Populus miRNA families, expression of some miRNAs are altered in response to cold, heat, salinity, dehydration, and mechanical stresses . Pathogenic bacteria, fungi, viruses, insect pests and nematodes cause severe damage to plants . Recent discovery also revealed that miR393-guided post-transcriptional regulation plays a crucial role in the plant defense against pathogens through targeting transport inhibitor response 1 (TIR1), an auxin receptor . Transgenic Arabidopsis overexpressing miR393a shows enhanced resistance to virulent P. syringae pv. tomato . Another report in loblolly pine indicated that expressions of 10 miRNAs are decreased in response to the rust fungus . Functional analyses have demonstrated that several plant miRNAs play vital roles in plant resistance to abiotic as well as biotic stresses [4, 33–37].
Although some of the stress-responsive miRNA families are deeply conserved among various plant species, including Arabidopsis, rice and Populus, their species-specific function may be the results of adaptation to long-term growth and survival in stressful environment . Moreover, some species-specific miRNAs could also play possible roles in the regulatory networks associated with the stress resistance, and individual miRNAs of a family response differentially to stress, suggesting that different members in same miRNA family may have different functions . Recently, a total of 58 miRNAs comprising 43 miRNA families have been cloned in wheat , however, their expression level in response to abiotic and biotic stress are still unknown. With the development of high-throughput sequencing technology, it became possible to discover several species-specific or lowly expressed miRNAs and different members in the same miRNA family. In this study, by using Solexa high-throughput sequencing, we indentified a diverse set of wheat small RNAs which are responsive to powdery mildew infection and heat stress. A total of 51 known conserved miRNAs and 81 new identified miRNAs were obtained, increasing the number of wheat miRNA families from 43 to 132. Moreover we also found that many of these wheat miRNAs showed differential expression in response to powdery mildew infection and heat stress. In addition, 149 genes were predicted as potential targets for novel wheat miRNAs, which included transcription factors implicated in development as well as genes involved in other physiological processes, such as stress responses.
High-throughput sequencing of wheat small RNAs
To investigate the role of wheat small RNA in powdery mildew resistance and heat stress, six small RNA libraries were constructed and sequenced. Disease susceptible wheat cultivar Jingdong8 (JD8) and its near-isogenic resistant line Jingdong8-Pm30 (JD8-Pm30) with single resistant gene Pm30 were used for powdery mildew infection and the preponderant physiological strain Erysiphe graminis f. sp. tritici (Egt) (Isolate E09) was used for powdery mildew infection. Non-inoculated and inoculated seedlings of JD8 and JD8-Pm30 were used to construct four small RNA libraries, designated JD8, JD8-Egt, JD8-Pm30 and JD8-Pm30-Egt, respectively. As for heat stress treatment, seedlings of heat tolerant cultivar TAM107 that undergone 40°C for 1 hour and those grown in normal condition were selected to build small RNA libraries, designated TAM107 and TAM107-heat, respectively.
These six wheat small RNA libraries were sequenced by Solexa high-throughput sequencing, and the numbers of raw reads for JD8, JD8-Egt, JD8-Pm30, JD8-Pm30-Egt, TAM107 and TAM107-heat were 3101553, 3848063, 2664384, 2560124, 4036437 and 4802099, respectively. To simplify the sequencing data, the identical sequence reads among the 6 small RNA libraries were grouped and converted into unique sequences. After removing the low quality reads, a total of 3839724 unique sequences ranging from 18-30 nt in length were obtained. Theoretically, these sequences could cover most of the small RNAs in our six samples used in this study.
Identification of known miRNAs
We searched for known miRNAs in the six small RNA libraries by use of homolog analysis to find miRNA sequences (orthologs/paralogs) matching at least 18 nt and leaving 3 nt for possible sequence variations. A total of 72 known miRNAs from 51 families were identified (Additional files 1 and 2). Our previous miRNA cloning by 454 sequencing in wheat have identified 43 miRNA families , among which 37 miRNA families were confirmed in current Solexa data set, however, other 6 families, including TamiR479, TamiR503, TamiR508, TamiR513, TamiR518, and TamiR522, were absent in our current Solexa data set. This is probably due to that the small RNA library for 454 sequencing was constructed from the pooled wheat RNA isolated from leaves, roots and spikes. In addition, 14 miRNA families, which were not found in the previous 454 sequencing data, were identified in present study with lower frequency. This indicated that these miRNAs were not abundant in wheat and could be detected only in the larger sequencing data set. Seven known rice miRNAs, including miR528, miR530, miR535, miR812, miR818, miR820, and miR827, which were not found in other species like Arabidopsis, were also represented in our data set, suggesting that they might belong to monocot-specific miRNAs.
Identification of novel wheat miRNAs
For the identification of novel wheat miRNAs, we firstly rely on wheat EST sequences as miRNA surrounding sequences in prediction, since the information of wheat genome sequence is limited. We found that a total of 330,925 small RNA sequences can be perfectly matched to at least one EST. The matched ESTs were then searched against Rfam and protein database to eliminate non-coding RNAs such as rRNA and tRNA as well as the degradation products from protein-coding sequences. With this analysis, we obtained 38167 sequences which are used to predict for fold-back RNA secondary structure. And we evaluated reads that fell within potential miRNA-like hairpins, considering the following five criteria: (1) the pairing characteristics of the hairpin; (2) the expression of the candidate, as measured by the abundance of sequence reads sharing the same 5' terminus; 3) presence in no less than two independent libraries, 4) the presence of miRNA* for several new microRNAs; 5) the absence of annotation suggesting non-miRNA biogenesis. These hairpins were further checked manually to ensure that they were accorded with the new criteria for annotation of plant miRNAs, provided recently by Meyers et al . Based on such analysis, we identified a total of 81 novel microRNA candidates (Additional files 3, 4 and 5). The length of the newly identified miRNAs range from 21 bp and 24 bp in length, and the negative folding free energies vary from -98.4 to -34.5 kcal mol-1 (with an average of -59.5 kcal mol-1) according to MFOLD, which is similar to the free energy values of other plant miRNA precursors (-72.4 kcal mol-1 in wheat, -71.0 kcal mol-1 in rice and -59.5 kcal mol-1 in Arabidopsis respectively). We also found the presence of 17 miRNA* in our libraries. Moreover, we took advantage of approximate 4× depth genomic sequences of Brachypodium distachyon, a grass species related to wheat, for further discovery of novel miRNAs. A total of 20 novel miRNAs were identified (Additional file 6), among which 16 (Additional files 6 and Additional file 4, underlined) can not be matched to wheat EST sequences due to the limitation of available wheat EST sequences.
In order to understand the conservation of new identified wheat miRNA, we compared wheat novel miRNAs to genomes of species representing important lineages, including maize, Arabidopsis, rice, Brachypodium distachyon and soybean. Among the 81 new miRNAs, 67 were wheat specific, and the others were conserved in different species. For example, miR2019 was conserved in monocots such as wheat, Brachypodium distachyon, rice and maize, suggesting that it might be monocot-specific miRNA. miR2072 and miR2050 were conserved in Brachypodium distachyon, wheat and rice; miR2070 was conserved in Brachypodium distachyon and wheat, miR2014 was conserved in rice and wheat.
A total of 149 genes were predicted as potential targets for these 81 new identified miRNAs (Additional file 7). Up to 62 new identified wheat miRNAs have predicted targets, whereas 10 miRNAs do not, which could be due to limitation of wheat EST sequences. We found that 52 new miRNAs have 139 targets with definite function, among which 16 encode proteins responding to disease or stress, 6 for transcription factor, 117 encode other proteins and most of them are involved in cellular metabolism. This is quite different as compared to the conserved miRNAs, as most of their targets are transcription factors. Interesting, we found that 28 target genes were putatively regulated by more than 2 miRNAs. It has also been widely suggested that microRNAs, similarly to transcription factors, can act in combination (or cooperatively) by binding to the same mRNA in a concentration-dependent manner .
Response of wheat miRNAs to powdery mildew infection and heat stress
Thirdly, we identified that 9 miRNAs are co-regulated by both powdery mildew infection and heat stress (Figure 3), among which two (miR827, miR2005) were up-regulated both in powdery mildew infection and heat stress, indicating that they might play important roles in both abiotic and biotic stress response in wheat. The others exhibited opposite expression pattern in response to powdery mildew infection and heat stress.
Expression patterns of potential target genes in response to powdery mildew infection and heat stress
Discussion and Conclusions
Identification of entire set of miRNAs and their targets will lay the foundation to unravel the complex miRNA-mediated regulatory networks controlling development and other physiological processes. Several computational studies estimated that organisms probably contain about 1% miRNA genes of the total protein-coding genes [34, 35, 45]. Recently, a large number of miRNA have been found in various species. For example, the identified number of miRNA in Arabidopsis, rice and maize were 187, 353 and 96, respectively. A total of 58 wheat miRNAs were identified by 454 sequencing, . In present study, our Solexa high-throughput sequencing of wheat small RNAs revealed a diverse and complex small RNA population and expression of the 51 known plant miRNA families (Additional files 1 and 2) and 81 novel miRNAs (Additional files 3 and 4) were determined. This increased the number of known wheat miRNA families from 43 to 132. Although the number of deeply conserved miRNA families in wheat largely remains the same as in Arabidopsis, all of the newly identified wheat miRNAs do not appear to be conserved in Arabidopsis and have predicted target genes with more diverse functions than those of conserved miRNAs. Recent deep sequencing of plant small RNA libraries also demonstrated that plants express more non-conserved than conserved miRNAs [8, 42, 46]. Only 5 miRNAs were conserved in different species, suggesting that wheat appears to have evolved species-specific miRNAs because of the functional diversification among the various species. Therefore, our study revealed that wheat genome encoded more non-conserved miRNA families than conserved miRNA families. It has been proposed that these non-conserved miRNAs presumably emerged and dissipated in short evolutionary time scales [42, 47], and such rapid emergence of new genes is likely facilitated by the small size and simple architecture of miRNA genes [47, 48].
By using Northern blotting, at least 3 and 5 wheat miRNAs families were found to be responsive to powdery mildew infection and heat stress, respectively. It is interesting to note that the members of the same miRNA family were differentially regulated in the response to powdery mildew infection and heat stress. This is consistent with results obtained from cold-stressed Populus , drought-stressed rice  and UV-B-treated Arabidopsis . Therefore, the functions of plant miRNAs can be dissimilar even if they share a high degree of sequence similarity and belong to the same family . In addition, several reported stress-responsive miRNAs in Arabidopsis and Populus are also involved in wheat response to powdery mildew infection and heat stress. For example, miR156 and miR160 were significantly repressed in the galled loblolly pine stem infected with the fungus C. quercuum f.sp.fusiforme , and miR156, miR164 and miR160 were induced in tobacco after plant virus's infection . In Arabidopsis, miR156 and miR164 were induced by infection with plant virus TYMV p69, and also induced in transgenic Arabidopsis plants expressing the viral silencing suppressor P1/HC-Pro . In this study, we also found that miR156 was significantly repressed by powdery mildew infection, whereas, miR156 was induced in response to heat stress. It was reported that miR393 was induced by a bacterial PAMP peptide flg22 at 20 to 60 min post-treatment . And in our study, expression of miR393 was also strongly increased in JD8 under infection of Egt. Moreover, we also found that 8 and 2 novel wheat miRNAs families were responsive to powdery mildew infection and heat stress. Taken together, these stress-related miRNAs identified in this study may function in one of the most critical defense systems for wheat biotic and abiotic stress tolerance, but the underlying molecular mechanism remains an area to be elucidated.
Auxin is an important hormone in plants, and involved in many plant processes including cell elongation and thus growth. Report shows that susceptible barley cultivar infected by powdery mildew have a remarkably increase of IAA level after inoculation  and simultaneously, IAA level decreased after the inoculation by parasite in the resistant barley cultivar "Professor Schiemann" . More recently, it was reported that miR393 was involved in auxin signaling pathway and played an important role in basal defense. Previous studies revealed that miR160, miR166, miR167 and miR393 may involve in auxin signaling pathway by regulating different transcription factors. One interesting observation in this study is that the expression patterns of miR393 related to auxin signaling pathway was significantly different between near isogenic JD8 and JD8-Pm30. Therefore, we speculated that auxin pathways affected by miRNAs could play important roles in powdery mildew disease resistance.
Hexaploid wheat (Triticum aestivum L.) line JD8 (susceptible to wheat powdery mildew) and its near-isogenic resistant line JD8-Pm30 were grown in a growth chamber at a relative humidity of 75% and 26/20°C day and night temperature with light intensity of 3000 lx. Seven-day-old plants were used for all experiments. One isolate of the powdery mildew fungus Erysiphe graminis f. sp. tritici (Egt) (Isolate E09) was maintained on the wheat cultivar Fidel by weekly transfer to new plants. Inoculations were performed at a density of 100-150 conidia/mm2 by brushing plants. Leaves were collected at 12 hours after infection.
Two wheat genotypes, heat susceptible 'Chinese Spring' (CS) and heat tolerant 'TAM107' were used in this study. Seeds were surface-sterilized in 1% sodium hypochlorite for 15 min, rinsed in distilled water, and soaked in dark overnight at room temperature. The germinated seeds were transferred into the pots (25 seedlings per pot) containing vermiculite. The seedlings were subjected to 40°C (heat stress) for 0.5, 1 and 2 hours and returned to normal conditions for 2 hours. At the end of heat treatments, the leaves of the seedlings were frozen immediately in the liquid nitrogen, and stored at -80°C for further use.
RNA extraction and small RNA cloning
Total RNA was isolated from the frozen leaves by using Trizol agent (TaKaRa, Inc., Dalian, China) according to the manufacturer's instructions. Cloning of the small RNAs was performed as described by Sunkar and Zhu . Briefly, low molecular weight RNA was enriched by 0.5 M NaCl and 10% PEG8000 precipitation. About 100 μg low molecular weight RNA was separated on a denaturing 15% polyacrylamide gel. Labeled RNA oligonucleotides with 18 and 26 nt were used as size standards. The nucleotides from 18 to 26 nt were excised, and RNA was eluted overnight with 0.4 M NaCl at 4°C. The RNA was dephosphorylated by alkaline phosphatase (New England Biolabs Inc, Beijing China) and recovered by ethanol precipitated. The small RNAs were then ligated sequentially to RNA/DNA chimeric oligonucleotide adapters, and then reverse transcription was preformed, followed by PCR amplification. The resulting PCR products were sequenced using Solexa technology.
Automated base calling of the raw sequence and vector removal were performed with PHRED and CROSS MATCH programs [23, 54]. Trimmed 3' and 5' adapters sequences, removed RNAs less than 17 nt and polyA, only sequences longer than 17 nt with a unique ID were used for further analysis. These sequences were used to search for the Rfam database  with BLASTN  to remove most non-siRNA and non-miRNA sequences. Putative origins for the remaining sequences were identified by BLASTN search against wheat EST database from NCBI. The protein-coding EST sequences were removed and the remaining non-coding candidate wheat ESTs with perfect matches with small RNA sequences were used for fold back secondary structure prediction with MFOLD program . In NCBI Unigene database, closely related wheat ESTs have been assembled to Unigene cluster, therefore the Unigene accessions were selected and recorded. Based on these analyses, putative miRNAs were then searched against NCBI NT database to check whether these miRNAs exist in other species. We also map the small RNA sequences to Brachypodium distachyon genomic sequences, which were downloaded from http://www.brachypodium.org/.
Target gene prediction was carried out as described by Jones-Rhoades et al [19, 58]. It was performed by searching the wheat EST database and NCBI NT database for miRNA complementary sequences. These criteria include allowing one mismatch in the region complementary to nucleotide positions 2 to 12 of the miRNA but not at the position 10/11 which is predicted cleavage site, and three additional mismatches were permitted between 12 and 22 nucleotide positions, but no more than two continuous mismatches within this region.
Differential expression analysis of miRNAs based on high-throughput sequencing
The miRNAs with TPM less than 20 showed the less expression level, which did not be included for the differential expression analysis.
RNA gel blot analysis
Low molecular weight RNA was loaded and resolved on a denaturing 15% polyacrylamide gel, and transferred electrophoretically to Hybond-N+ membranes (Amersham Biosciences, Buckinghamshire, UK). Membranes were UV cross-linked and baked for 2 hours at 80°C. DNA oligonucleotides complementary to miRNA sequences were end labeled with γ-32P-ATP using T4 polynucleotide kinase (TaKaRa, Dalian, China). Membranes were prehybridized for more than 8 hours and hybridized overnight using Church buffer at 38°C. Blots were washed three times (two times with 2 × SSC + 1% SDS and one time with 1 × SSC + 0.5% SDS) at 50°C. The membranes were briefly air dried and then exposed to X-ray films for autograph at -80°C.
Infiltration of Agrobacterium tumefaciens intoN. Benthamiana
Precursor sequences of miR156 including the hairpin structures, the Ta3711 and mTa3711 were cloned to downstream of 35 S promoter, respectively. A mutated version of the Ta3711 transgene (mTa3711) was generated by PCR. The mutated mTa3711 primers used were as follows: 5'-ATCTTCAGCACGCACTGTCACTACTCTCTAGCAACCCAGTGG-3 'and 5'-CCACTGGGTTGCTAGAGAGTAGTGACAGTGCGTGCTGAAGAT-3'. The nucleotide sequence in miR156 binding domain of the mTa3711 was changed from CATGCTCTCTCTCTTCTGTCA to CACGCACTGTCACTACTCTCT.
35S::miR156, 35S::Ta3711, 35S::mTa3711 were introduced into A. tumefaciens and the bacteria injected into N. benthamiana leaves with a syringe according to the method of Llave . For co-injections of two different constructs, bacteria were resuspended in infiltration medium (0.5' Murashige and Skoog salts, 5% sucrose, 0.5 g/l Mes) at OD600 = 1, and incubated for 3 hours at room temperature with 150 mM acetosyringone. Zones of infiltration were harvested at 4 and 5 days after injection for RNA isolations.
Semi-quantitative RT-PCR analysis
Total RNA was isolated by using Trizol (Invitrogen, USA) according to the manufacturer's instructions and treated with RNase-free DNase I (Promega, Madison, USA). Two microgram of total RNA from each sample was used for first-strand cDNA synthesis in 20 μl reactions containing 50 mM Tris-HCl (pH 8.3), 75 mM KCl, 3 mM MgCl2, 10 mM DTT, 50 μM dNTPs, 200 U M-MLV reverse transcriptase (Promega, Madison, USA) and 50 pmol oligonucleotides T15. Reverse transcription was performed at 37°C for 60 min with a final denaturation at 95°C for 5 min. Gene-specific RT-PCR primers for 5 miRNA targets were designed according to the EST sequences.
Three RT-PCR replications were conducted using independently isolated total RNAs with the following thermal cycling parameters: 94°C for 30 sec, 57°C for 30 sec, and 72°C for 30 sec. A 350 bp β-actin gene fragment was amplified as a positive control using the primer pair 5'-CAGCAACTGGGATGATATGG-3' and 5'-ATTTCGCTTTCAGCAGTGGT-3'. The RT-PCR products were sequenced to verify the specificity of PCR amplifications. The primers used for RT-PCR were listed in the Additional file 9.
Realtime PCR analysis
The TaGAMYB1 gene specific primer was designed using DNAMAN software as followings: 5'-CCTGAATTGAGCGACACC' and 5'-ACACACTCCGACTT-CACTG-3'. The TaGAMYB2 gene specific primer was designed using DNAMAN software as followings: 5'-GTGGGGCGATTTCATTGAT-3' and 5'-TTGAGTA-TGCGGACCAGTTG-3'. PCR reactions were performed in a volume of 20 μl containing 10 pmol primers 10 mmol/l Tris-HCl)pH8.5), 50 mmol/l KCl, 2 mmol/l MgCl2, 0.4 μl DMSO, 200 mmol/l dNTPs, 1 U Taq DNA Polymerase (TaKaRa, Dalian), 0.5 μl SYBR GREEN I. PCR amplification protocol was followed by 95°C (3 min) and 40 cycles of amplification cycle (95°C (30 s), 55°C (30 s), and 72°C (1 min)) using Opticon PTC200 system (MJ Research, USA).
All reactions were run in triplicate and included no template and no reverse transcription controls. Quantification results were expressed in terms of the cycle threshold (CT) value according to the baseline adjusted to 0.04. The comparative CT method (PE Applied Biosystems) was used to quantify relative gene expression compared with actin. Briefly, the CT values were averaged for each triplicate. Differences between the mean CT values of TaGAMYB and those of β-actin were calculated as ΔCTsample = CT TaGAMYB - CTβ-actin. Final results, the sample relative expression level were determined as 2 -ΔCtsample. Statistical significance was tested using the Student's t-test (P < 0.05).
This work was financially supported by National Basic Research Program of China (2007CB109000), 863 Project of China (2007AA10Z138,2006AA10A104) and National Natural Science Foundation of China (30871529,30871528).
- Gill BS, Appels R, Botha-Oberholster AM, Buell CR, Bennetzen JL, Chalhoub B, Chumley F, Dvorak J, Iwanaga M, Keller B, Li W, McCombie WR, Ogihara Y, Quetier F: A workshop report on wheat genome sequencing: International Genome Research on Wheat Consortium. Genetics. 2004, 168: 1087-1096. 10.1534/genetics.104.034769.PubMedPubMed CentralView ArticleGoogle Scholar
- Griffey CA, Das MK, Stromberg EL: ffectiveness ofadult-plant resistance in reducing grain yield loss to powderymildew in winter wheat. Plant Dis. 1993, 77: 618-622.View ArticleGoogle Scholar
- Leath S, Bowen KL: Effects of powdery mildew, triadimenolseed treatment, and triadimefon foliar sprays on yieldof winter wheat in North Carolina. Phytopathology. 1989, 79: 152-155. 10.1094/Phyto-79-152.View ArticleGoogle Scholar
- Wang X, Zhang J, Li F, Gu J, He T, Zhang X, Li Y: MicroRNA identification based on sequence and structure alignment. Bioinformatics. 2005, 21: 3610-3614. 10.1093/bioinformatics/bti562.PubMedView ArticleGoogle Scholar
- Shukla LI, Chinnusamy V, Sunkar R: The role of microRNAs and other endogenous small RNAs in plant stress responses. Biochim Biophys Acta. 2008, 1779: 743-748.PubMedView ArticleGoogle Scholar
- Bruggmann R, Abderhalden O, Reymond P, Dudler R: Analysis of epidermis- and mesophyll-specific transcript accumulation in powdery mildew-inoculated wheat leaves. Plant Mol Biol. 2005, 58: 247-267. 10.1007/s11103-005-3099-9.PubMedView ArticleGoogle Scholar
- Qin D, Wu H, Peng H, Yao Y, Ni Z, Li Z, Zhou C, Sun Q: Heat stress-responsive transcriptome analysis in heat susceptible and tolerant wheat (Triticum aestivum L.) by using Wheat Genome Array. BMC Genomics. 2008, 9: 432-10.1186/1471-2164-9-432.PubMedPubMed CentralView ArticleGoogle Scholar
- Sunkar R, Chinnusamy V, Zhu J, Zhu JK: Small RNAs as big players in plant abiotic stress responses and nutrient deprivation. Trends Plant Sci. 2007, 12: 301-309. 10.1016/j.tplants.2007.05.001.PubMedView ArticleGoogle Scholar
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.PubMedView ArticleGoogle Scholar
- Baulcombe D: RNA silencing in plants. Nature. 2004, 431: 356-363. 10.1038/nature02874.PubMedView ArticleGoogle Scholar
- Carrington JC, Ambros V: Role of microRNAs in plant and animal development. Science. 2003, 301: 336-338. 10.1126/science.1085242.PubMedView ArticleGoogle Scholar
- Chen PY, Meister G: microRNA-guided posttranscriptional gene regulation. Biol Chem. 2005, 386: 1205-1218. 10.1515/BC.2005.139.PubMedView ArticleGoogle Scholar
- Lee RC, Ambros V: An extensive class of small RNAs in Caenorhabditis elegans. Science. 2001, 294: 862-864. 10.1126/science.1065329.PubMedView ArticleGoogle Scholar
- Matzke M, Aufsatz W, Kanno T, Daxinger L, Papp I, Mette MF, Matzke AJ: Genetic analysis of RNA-mediated transcriptional gene silencing. Biochim Biophys Acta. 2004, 1677: 129-141.PubMedView ArticleGoogle Scholar
- Aukerman MJ, Sakai H: Regulation of flowering time and floral organ identity by a MicroRNA and its APETALA2-like target genes. Plant Cell. 2003, 15: 2730-2741. 10.1105/tpc.016238.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen CZ, Li L, Lodish HF, Bartel DP: MicroRNAs modulate hematopoietic lineage differentiation. Science. 2004, 303: 83-86. 10.1126/science.1091903.PubMedView ArticleGoogle Scholar
- Palatnik JF, Allen E, Wu X, Schommer C, Schwab R, Carrington JC, Weigel D: Control of leaf morphogenesis by microRNAs. Nature. 2003, 425: 257-263. 10.1038/nature01958.PubMedView ArticleGoogle Scholar
- Kidner CA, Martienssen RA: Spatially restricted microRNA directs leaf polarity through ARGONAUTE1. Nature. 2004, 428: 81-84. 10.1038/nature02366.PubMedView ArticleGoogle Scholar
- Jones-Rhoades MW, Bartel DP: Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol Cell. 2004, 14: 787-799. 10.1016/j.molcel.2004.05.027.PubMedView ArticleGoogle Scholar
- Vaucheret H, Vazquez F, Crete P, Bartel DP: The action of ARGONAUTE1 in the miRNA pathway and its regulation by the miRNA pathway are crucial for plant development. Genes Dev. 2004, 18: 1187-1197. 10.1101/gad.1201404.PubMedPubMed CentralView ArticleGoogle Scholar
- Xie Z, Kasschau KD, Carrington JC: Negative feedback regulation of Dicer-Like1 in Arabidopsis by microRNA-guided mRNA degradation. Curr Biol. 2003, 13: 784-789. 10.1016/S0960-9822(03)00281-1.PubMedView ArticleGoogle Scholar
- Zhou X, Wang G, Sutoh K, Zhu JK, Zhang W: Identification of cold-inducible microRNAs in plants by transcriptome analysis. Biochim Biophys Acta. 2008, 1779: 780-788.PubMedView ArticleGoogle Scholar
- Sunkar R, Zhu JK: Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis. Plant Cell. 2004, 16: 2001-2019. 10.1105/tpc.104.022830.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhao Y, Srivastava D: A developmental view of microRNA function. Trends Biochem Sci. 2007, 32: 189-197. 10.1016/j.tibs.2007.02.006.PubMedView ArticleGoogle Scholar
- Sunkar R, Kapoor A, Zhu JK: Posttranscriptional induction of two Cu/Zn superoxide dismutase genes in Arabidopsis is mediated by downregulation of miR398 and important for oxidative stress tolerance. Plant Cell. 2006, 18: 2051-2065. 10.1105/tpc.106.041673.PubMedPubMed CentralView ArticleGoogle Scholar
- Aung K, Lin SI, Wu CC, Huang YT, Su CL, Chiou TJ: pho2, a phosphate overaccumulator, is caused by a nonsense mutation in a microRNA399 target gene. Plant Physiol. 2006, 141: 1000-1011. 10.1104/pp.106.078063.PubMedPubMed CentralView ArticleGoogle Scholar
- Bari R, Datt PB, Stitt M, Scheible WR: PHO2, microRNA399, and PHR1 define a phosphate-signaling pathway in plants. Plant Physiol. 2006, 141: 988-999. 10.1104/pp.106.079707.PubMedPubMed CentralView ArticleGoogle Scholar
- Chiou TJ, Aung K, Lin SI, Wu CC, Chiang SF, Su CL: Regulation of phosphate homeostasis by MicroRNA in Arabidopsis. Plant Cell. 2006, 18: 412-421. 10.1105/tpc.105.038943.PubMedPubMed CentralView ArticleGoogle Scholar
- Fujii H, Chiou TJ, Lin SI, Aung K, Zhu JK: A miRNA involved in phosphate-starvation response in Arabidopsis. Curr Biol. 2005, 15: 2038-2043. 10.1016/j.cub.2005.10.016.PubMedView ArticleGoogle Scholar
- Lu YD, Gan QH, Chi XY, Qin S: Roles of microRNA in plant defense and virus offense interaction. Plant Cell Rep. 2008, 27: 1571-1579. 10.1007/s00299-008-0584-z.PubMedView ArticleGoogle Scholar
- Navarro L, Dunoyer P, Jay F, Arnold B, Dharmasiri N, Estelle M, Voinnet O, Jones JD: A plant miRNA contributes to antibacterial resistance by repressing auxin signaling. Science. 2006, 312: 436-439. 10.1126/science.1126088.PubMedView ArticleGoogle Scholar
- Lu S, Sun YH, Amerson H, Chiang VL: MicroRNAs in loblolly pine (Pinus taeda L.) and their association with fusiform rust gall development. Plant J. 2007, 51: 1077-1098. 10.1111/j.1365-313X.2007.03208.x.PubMedView ArticleGoogle Scholar
- Bonnet E, Wuyts J, Rouze P, Van de PY: Detection of 91 potential conserved plant microRNAs in Arabidopsis thaliana and Oryza sativa identifies important target genes. Proc Natl Acad Sci USA. 2004, 101: 11511-11516. 10.1073/pnas.0404025101.PubMedPubMed CentralView ArticleGoogle Scholar
- Lim LP, Glasner ME, Yekta S, Burge CB, Bartel DP: Vertebrate microRNA genes. Science. 2003, 299: 1540-10.1126/science.1080372.PubMedView ArticleGoogle Scholar
- Lim LP, Lau NC, Weinstein EG, Abdelhakim A, Yekta S, Rhoades MW, Burge CB: The microRNAs of Caenorhabditis elegans. Genes Dev. 2003, 17: 991-1008. 10.1101/gad.1074403.PubMedPubMed CentralView ArticleGoogle Scholar
- Mallory AC, Vaucheret H: Functions of microRNAs and related small RNAs in plants. Nat Genet. 2006, 38 (Suppl): S31-S36. 10.1038/ng1791.PubMedView ArticleGoogle Scholar
- Weber MJ: New human and mouse microRNA genes found by homology search. FEBS J. 2005, 272: 59-73. 10.1111/j.1432-1033.2004.04389.x.PubMedView ArticleGoogle Scholar
- Yao Y, Guo G, Ni Z, Sunkar R, Du J, Zhu JK, Sun Q: Cloning and characterization of microRNAs from wheat (Triticum aestivum L.). Genome Biol. 2007, 8: R96-10.1186/gb-2007-8-6-r96.PubMedPubMed CentralView ArticleGoogle Scholar
- Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, Green PJ, Griffiths-Jones S, Jacobsen SE, Mallory AC, Martienssen RA, Poethig RS, Qi Y, Vaucheret H, Voinnet O, Watanabe Y, Weigel D, Zhu JK: Criteria for Annotation of Plant MicroRNAs. Plant Cell. 2008, 19:Google Scholar
- Ambros V, Bartel B, Bartel DP, Burge CB, Carrington JC, Chen X, Dreyfuss G, Eddy SR, Griffiths-Jones S, Marshall M, Matzke M, Ruvkun G, Tuschl T: A uniform system for microRNA annotation. RNA. 2003, 9: 277-279. 10.1261/rna.2183803.PubMedPubMed CentralView ArticleGoogle Scholar
- Grun D, Wang YL, Langenberger D, Gunsalus KC, Rajewsky N: microRNA target predictions across seven Drosophila species and comparison to mammalian targets. PLoS Comput Bio. 2005, 1: e13-10.1371/journal.pcbi.0010013.View ArticleGoogle Scholar
- Fahlgren N, Howell MD, Kasschau KD, Chapman EJ, Sullivan CM, Cumbie JS, Givan SA, Law TF, Grant SR, Dangl JL, Carrington JC: High-throughput sequencing of Arabidopsis microRNAs: evidence for frequent birth and death of MIRNA genes. PLoS ONE. 2007, 2: e219-10.1371/journal.pone.0000219.PubMedPubMed CentralView ArticleGoogle Scholar
- Johnson C, Bowman L, Adai AT, Vance V, Sundaresan V: CSRDB: a small RNA integrated database and browser resource for cereals. Nucleic Acids Res. 2007, 35: D829-D833. 10.1093/nar/gkl991.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhu QH, Spriggs A, Matthew L, Fan L, Kennedy G, Gubler F, Helliwell C: A diverse set of microRNAs and microRNA-like small RNAs in developing rice grains. Genome Res. 2008, 18: 1456-1465. 10.1101/gr.075572.107.PubMedPubMed CentralView ArticleGoogle Scholar
- Lai EC, Tomancak P, Williams RW, Rubin GM: Computational identification of Drosophila microRNA genes. Genome Biol. 2003, 4: R42-10.1186/gb-2003-4-7-r42.PubMedPubMed CentralView ArticleGoogle Scholar
- Rajagopalan R, Vaucheret H, Trejo J, Bartel DP: A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes Dev. 2006, 20: 3407-3425. 10.1101/gad.1476406.PubMedPubMed CentralView ArticleGoogle Scholar
- Sunkar R, Zhou X, Zheng Y, Zhang W, Zhu JK: Identification of novel and candidate miRNAs in rice by high throughput sequencing. BMC Plant Biol. 2008, 8: 25-10.1186/1471-2229-8-25.PubMedPubMed CentralView ArticleGoogle Scholar
- Allen E, Xie Z, Gustafson AM, Sung GH, Spatafora JW, Carrington JC: Evolution of microRNA genes by inverted duplication of target gene sequences in Arabidopsis thaliana. Nat Genet. 2004, 36: 1282-1290. 10.1038/ng1478.PubMedView ArticleGoogle Scholar
- Zhao B, Liang R, Ge L, Li W, Xiao H, Lin H, Ruan K, Jin Y: Identification of drought-induced microRNAs in rice. Biochem Biophys Res Commun. 2007, 354: 585-590. 10.1016/j.bbrc.2007.01.022.PubMedView ArticleGoogle Scholar
- Zhou X, Wang G, Zhang W: UV-B responsive microRNA genes in Arabidopsis thaliana. Mol Syst Biol. 2007, 3: 103-10.1038/msb4100143.PubMedPubMed CentralView ArticleGoogle Scholar
- Bazzini AA, Hopp HE, Beachy RN, Asurmendi S: Infection and coaccumulation of tobacco mosaic virus proteins alter microRNA levels, correlating with symptom and plant development. Proc Natl Acad Sci USA. 2007, 104: 12157-12162. 10.1073/pnas.0705114104.PubMedPubMed CentralView ArticleGoogle Scholar
- Kasschau KD, Xie Z, Allen E, Llave C, Chapman EJ, Krizan KA, Carrington JC: P1/HC-Pro, a viral suppressor of RNA silencing, interferes with Arabidopsis development and miRNA unction. Dev Cell. 2003, 4: 205-217. 10.1016/S1534-5807(03)00025-X.PubMedView ArticleGoogle Scholar
- Vizarova G: Study of the level of the indole type of growth regulators in barley during the pathogenesis of the powdery mildew (Erysiphe graminis f. sp. hordei Marchal). Biologia (Bratisl). 1968, 23: 89-94.Google Scholar
- Sunkar R, Girke T, Jain PK, Zhu JK: Cloning and characterization of microRNAs from rice. Plant Cell. 2005, 17: 1397-1411. 10.1105/tpc.105.031682.PubMedPubMed CentralView ArticleGoogle Scholar
- Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen Z, Dewell SB, Du L, Fierro JM, Gomes XV, Godwin BC, He W, Helgesen S, Ho CH, Irzyk GP, Jando SC, Alenquer ML, Jarvie TP, Jirage KB, Kim JB, Knight JR, Lanza JR, Leamon JH, Lefkowitz SM, Lei M, Li J, Lohman KL, Lu H, Makhijani VB, McDade KE, McKenna MP, Myers EW, Nickerson E, Nobile JR, Plant R, Puc BP, Ronan MT, Roth GT, Sarkis GJ, Simons JF, Simpson JW, Srinivasan M, Tartaro KR, Tomasz A, Vogt KA, Volkmer GA, Wang SH, Wang Y, Weiner MP, Begley RF, Rothberg JM: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437: 376-380.PubMedPubMed CentralGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31: 3406-3415. 10.1093/nar/gkg595.PubMedPubMed CentralView ArticleGoogle Scholar
- Rhoades MW, Reinhart BG, Lim LP, Burge CB, Bartel B, Bartel DP: Prediction of plant miRNA targets. Cell. 2002, 110: 513-520. 10.1016/S0092-8674(02)00863-2.PubMedView ArticleGoogle Scholar
- Llave C, Xie Z, Kasschau KD, Carrington JC: Cleavage of Scarecrow-like mRNA targets directed by a class of Arabidopsis miRNA. Science. 2002, 297: 2053-2056. 10.1126/science.1076311.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.