High-throughput sequencing reveals microRNAs and their targets in response to drought stress in wheat (Triticum aestivum L.)

Figures & data

Table 1. Overview of small RNA sequences for two wheat (Triticum aestivum L. cv. XF 20) libraries.

Data type	Control		Drought
	Total reads	Unique reads	Total reads	Unique reads
Clean reads	11723790	1508601	12034102	1357647
Genome	9365098	904379	10151515	918905
Rfam	3935209	166060	3590058	114691
Known miRNA	814254	18198	1287494	27292
Novel miRNA	12972	1025	12512	1124

Note: This table summarizes the number of total and unique reads for both the control and drought-stressed plant cDNA libraries. Clean reads, reads mapped to the genome of Brachypodium distachyon, reads mapped to Rfam, known miRNA and novel miRNA read numbers are listed.

Table 2. Differentially expressed miRNA genes in wheat (Triticum aestivum L. cv. XF 20) in response to drought stress.

miRNA name	Control	Drought	Fold	p value
miR165a-5p	3551.50	1688.12	−1.07	**
miR419	167.01	1.00	−7.38	**
miR482a-5p	1192.36	1.00	−10.22	**
miR530b	50.15	1.00	−5.65	**
miR533e	158.40	1.00	−7.31	**
miR858a	1.00	21.52	4.43	**
miR1023a-3p	47.08	20.28	−1.22	**
miR1120c-5p	1.00	11.72	3.55	**
miR1127b-3p	1.00	29.25	4.87	**
miR1140	649.02	1.00	−9.34	**
miR1918	23.54	1.00	−4.56	**
miR2111c	192.69	1.00	−7.59	**
miR4391	10.41	1.00	−3.38	**
miR5037c	766.99	102.29	−2.91	**
miR5054	90.33	27.17	−1.73	**
miR5060	52.80	1.00	−5.72	**
miR5374-5p	237.12	1.00	−7.89	**
miR5377	127.26	1.00	−6.99	**
miR5384-3p	24.99	7.23	−1.79	**
miR5662	935.11	105.87	−3.14	**
miR5800	35.40	13.13	−1.43	**
miR6117	1.00	125.81	6.98	**
miR6135f	1003.69	1.00	−9.97	**
miR7513	64.31	27.34	−1.23	**
miR7696a-3p	1.00	7367.56	12.85	**
miR7731-3p	37.53	1.00	−5.23	**
miR7810	541.72	181.57	−1.58	**
miR8651	1.00	305.63	8.26	**
miR8714	27.21	6.15	−2.15	**
miR9674a-3p	14.59	1.00	−3.87	**
miR9674b-3p	1.00	9.64	3.27	**
miR9763	1.00	1733.16	10.76	**
NmiRNA07	33.52	1.00	−5.07	**
NmiRNA08	20.73	1.00	−4.37	**
NmiRNA15	19.45	1.00	−4.28	**
NmiRNA23	1.00	71.38	6.16	**
NmiRNA25	1.00	11.30	3.50	**

Note The list includes all the miRNAs that are significantly differentially expressed between control and drought-stressed plants. Mean expression levels are listed in transcripts per kilobase million for each group, and the fold-change in expression in drought-stressed plants is shown relative to control plants. The table is ordered by miRNA name. The chi-square test was used to calculate p values.

** p < 0.01 indicates significant differences.

Figure 1. Gene ontology analysis of miRNA target genes identified in drought stress of wheat (Triticum aestivum L. cv. XF 20). Note: The digits on the left y-axis and right y-axis show the percentage and number of miRNA targets, respectively. Targets of all expressed miRNAs are shown in cyan; targets of all differentially expressed miRNAs are shown in pink.

Data availability

The Triticum aestivum (cv. XF 20) small RNA sequences, degradome sequences, and transcriptome sequences have all been deposited in the NCBI SRA database with accession number SRP 166912.