水稻CYP81A6基因干扰对其他基因表达的影响

Knockdown of CYP81A6 Gene Affects Expression of Other Genes in Rice

为了明确抑制基因CYP81A6的表达对水稻的自身功能基因表达的影响,对亲本水稻秀水110和转基因水稻450-3的mRNA进行深度测序,并将获得的序列信息通过水稻基因数据库和生物信息学手段进行注释和功能基因比较。结果表明,在3529个差异表达的基因中,上调表达基因1018个,下调表达基因2511个;秀水110与450-3中共同表达的基因有2533个,仅在450-3中表达的基因有363个,仅在秀水110表达的基因有633个。与水稻基因数据库功能已知的296个抗性基因对比,发现450-3中差异表达的基因有64个功能已知且与水稻抗性相关,包括2个抗虫相关基因,17个抗病相关基因和45个抗逆相关基因,预示着基因CYP81A6的抑制表达对水稻抗性会产生一定的影响。GO功能和通路显著性富集分析显示,差异表达基因显著富集在38个GO条目和13个通路上,在结合能力(GO:0020037)和光合作用(KO 00195)通路上富集最为显著。通过对转基因水稻的转录组分析,阐明转基因水稻基因与亲本相比发生的变化,进而分析这些变化对转基因水稻产生的影响,从基因生态学角度为转基因水稻的安全性评价提供科学依据,可作为一种新的方法应用于转基因生物的安全性评价。

To study the effects of knockdown the gene CYP81A6 on gene expression of transgenic rice, we report the transcriptome profiling analysis of transgenic rice 450-3 and normal control rice Xiushui 110. In this study, mRNA of transgenic rice and its parent （control） were deeply sequenced, sequences information obtained and functional genes comparison against the rice gene database were analyzed. Results showed that 3529 genes in total were significantly differentialexpressed, among which 1018 genes were upregulated and 2511 genes were downregulated. The differential expressed genes distribution of the two samples showed that there were 2533 genes expressed in both Xiushui 110 and 450-3, 633 genes only expressed in Xiushui 110, while 363 genes expressed just in 450-3. Comparative analysis with 296 resistance genes in the rice database showed that 64 resistance genes were involved in rice resistance, including two insect-resistant genes, 17 disease-resistant genes and 45 stress-resistant genes, indicating that resistance of transgenic plant 450-3 would be affected. Gene enrichment studies and further analysis showed that differential expressed genes were significantly enriched in 38 GO terms and 13 KO terms in 450-3, ＇heme binding＇ （GO： 0020037） and ＇photosynthesis＇ （KO 00195） were enriched most dramatically. Comparative transcriptome analysis of the transgenic rice 450-3 and Xiushui 110 was conducted to clarify the changes of gene expression level in 450-3. Knockdown of gene CYP81A6 may result in extensive and variable change in the transcriptome of 450-3. Viewed from the gene evology point, these observations may provide valuable information to address the safety of genetically modified crops, and it could be used as a new method for the safety evaluation of genetically modified organisms.