High Throughput Sequencing Analysis on Mirna in Root of Iris Lactea Var. Chinensis Response to Cd Stress
LIU Liangqin,SONG Aiping,ZHANG Yongxia,YUAN Haiyan,HUANG Suzhen,LIU Zhaolei,GU Chunsun
DOI: https://doi.org/10.3969/j.issn.1674-7895.2016.03.01
2016-01-01
Abstract:In order to understand the expression pattern of miRNA in root of Iris lactea var. chinensis ( Fisch.) Koidz. under Cd stress, sRNA libraries ( which is CK and Cd libraries, respectively) from root of I. lactea var. chinensis after stressed by 100 μmol·L-1 Cd for 0 ( CK) and 24 h ( Cd) were analyzed by method of high throughput sequencing analysis. Also, miRNA with significantly differential expression was screened, and functions of target genes of these miRNA were predicted. On this basis, expression pattern of some miRNA and their target genes were verified by qRT-PCR technology. The results show that in CK and Cd libraries, there are more unannotation sRNA sequences, accounting for 86. 4% and 80. 5% of total number of their own sRNA unique reads, respectively. Among annotation sRNA sequences, percentage of miRNA is the lowest with a value of 0. 3% and 0. 5%, respectively, while that of rRNA is the highest with a value of 9. 4% and 11. 8%, respectively. Length of sRNA in the two libraries is mainly 21-24 nt, and the most of them is 21 nt. Thirty-two miRNA with significantly differential expression are screened from sRNA in root of I. lactea var. chinensis under Cd stress, in which, expression of twenty miRNA is down-regulated ( belonging to miR165, miR166, miR167, miR168, miR390 and miR396 families, respectively), that of twelve miRNA is up-regulated. The result of function prediction indicates that functions of target genes of these miRNA are mainly concentrated in three aspects, i. e. biological process, cellular component and molecular function, while from KEGG pathway enrichment analysis, number of target gene of differential expression miRNA enriching in three pathways of ribosome, biosynthesis of amino acids and carbon metabolism is 122, 88 and 82, respectively. Verification result of qRT-PCR shows that between CK and Cd libaries, there are significantly differences in relative expression of eleven miRNA with differential expression and eight target genes (P<0. 05). In which, up-regulated and down-regulated trends of relative expression of eleven miRNA are consistent with above screening result, and when relative expression of miRNA is up-regulated, that of its target gene is down-regulated, and vice versa, meaning that under Cd stress condition, miRNA from root of I. lactea var. chinensis negatively regulates its target gene expression, and these target genes mainly are involved in processes of coding transcription factor, HD-ZIP protein and signaling protein, etc.