長(zhǎng)白山人參microRNA鑒定及靶基因分析研究
發(fā)布時(shí)間：2019-05-23 09:20
【摘要】：目的:以15年生石柱參、6年生園參為研究對(duì)象,運(yùn)用高通量測(cè)序技術(shù)構(gòu)建人參small RNA、轉(zhuǎn)錄組及降解組數(shù)據(jù)庫(kù)。鑒定兩種人參中所含miRNA的種類,篩選差異表達(dá)miRNA,并進(jìn)行靶基因預(yù)測(cè)、驗(yàn)證及分析。方法:1、選用Trizol法抽提園參、石柱參總RNA,經(jīng)瓊脂糖凝膠電泳、TGem微量分光光度計(jì)、Agilent 2100 Bioanalyzer檢測(cè)總RNA的質(zhì)量、純度及完整性。2、利用Illumina HiSeqTM 2000測(cè)序儀進(jìn)行園參、石柱參轉(zhuǎn)錄組測(cè)序,經(jīng)短序列組裝拼接獲得兩樣本的Non-All-Unigene,將其與Nr、KOG、KEGG、GO等數(shù)據(jù)庫(kù)比對(duì)分析,得到Unigene的蛋白功能注釋及代謝通路信息。篩選樣本間差異表達(dá)mRNA,進(jìn)行GO、KEGG富集分析,利用Q-PCR對(duì)差異mRNA進(jìn)行驗(yàn)證。3、構(gòu)建園參、石柱參small RNA文庫(kù),以Unigene為參考基因組序列。將獲得的Clean reads與Rfam、Gene、Repbase及miRBase 21.0數(shù)據(jù)庫(kù)比對(duì),得到已知miRNA并統(tǒng)計(jì)其表達(dá)量;未比對(duì)上的序列用于新miRNA的預(yù)測(cè)。篩選樣本間差異miRNA,進(jìn)行靶基因預(yù)測(cè)及KEGG、GO富集分析,采用Q-PCR對(duì)差異miRNA進(jìn)行驗(yàn)證。4、構(gòu)建兩樣本降解組數(shù)據(jù)庫(kù),獲得原始序列與Rfam、Genbank、PolyN等數(shù)據(jù)庫(kù)比對(duì),得到cDNA-sense。cDNA-sense與Unigene、已知miRNA進(jìn)行miRNA降解位點(diǎn)的鑒定與分析后,再與Nr、KEGG和GO數(shù)據(jù)庫(kù)比對(duì),獲得降解基因的功能注釋信息。結(jié)果:1、瓊脂糖凝膠電泳顯示,兩樣本總RNA的28S、18S條帶清晰完整,亮度接近2倍;Agilent 2100 Bioanalyzer質(zhì)檢可見(jiàn)OD260nm/OD280nm比值在1.8~2.1范圍,OD260nm/OD280nm比值在2.0~2.3,RIN7.0,RNA濃度900 ng/μL,28S/18S1.0,即總RNA質(zhì)量滿足高通量建庫(kù)的標(biāo)準(zhǔn)。2、兩樣本轉(zhuǎn)錄組測(cè)序獲得5千多萬(wàn)條reads,經(jīng)De Novo拼接獲得63875條去冗余Unigene。對(duì)Unigene進(jìn)行功能注釋,有19340條Unigene被歸類到25個(gè)組蛋白直系KOG類別,20778條Unigene被歸類為64個(gè)GO功能類別,7268條Unigene被歸類到207個(gè)KEGG Pathway代謝通路。根據(jù)基因表達(dá)差異倍數(shù)FoldChange及P-value篩選出3216條差異mRNA,并選取差異顯著的18條基因進(jìn)行Q-PCR驗(yàn)證,驗(yàn)證結(jié)果與轉(zhuǎn)錄本表達(dá)量基本一致。3、通過(guò)small RNA高通量測(cè)序,在石柱參中鑒定235種已知miRNA,園參有246種已知miRNA,分屬于71個(gè)miRNA家族;未被任何數(shù)據(jù)庫(kù)比對(duì)上的序列進(jìn)行新miRNA的預(yù)測(cè),得到39個(gè)新miRNA。對(duì)已知和新miRNA進(jìn)行靶基因預(yù)測(cè),得到17個(gè)miRNA家族的169個(gè)靶基因,對(duì)這些靶基因進(jìn)行富集分析,按其功能注釋主要包括:蛋白轉(zhuǎn)錄因子、編碼與植物生長(zhǎng)發(fā)育有關(guān)的蛋白、生長(zhǎng)素響應(yīng)因子、泛素連接酶及非特征性蛋白等。根據(jù)miRNA表達(dá)差異倍數(shù)及P值,選取13條差異miRNA進(jìn)行Q-PCR驗(yàn)證,驗(yàn)證結(jié)果與miRNA表達(dá)量基本一致。4、通過(guò)構(gòu)建降解組文庫(kù)及與數(shù)據(jù)庫(kù)比對(duì)分析,在石柱參和園參中分別得到1千余萬(wàn)條cDNA-sense,用以降解位點(diǎn)的鑒定。根據(jù)cleaveland峰值,石柱參、園參中分別檢測(cè)到812個(gè)、724個(gè)剪切位點(diǎn)。結(jié)合生物信息學(xué)miRNA靶基因預(yù)測(cè)方法,對(duì)cDNAsense的降解位點(diǎn)進(jìn)行注釋,共得到8個(gè)miRNA家族的13個(gè)靶基因,其中兩個(gè)屬于新miRNA的靶基因。結(jié)論:1、成功構(gòu)建人參轉(zhuǎn)錄基因組、small RNA及降解組文庫(kù),鑒定石柱參、園參所含microRNA的種類,證實(shí)新鮮人參中富含microRNA。2、不同生長(zhǎng)年限人參miRNA存在差異表達(dá),顯著差異miRNA主要調(diào)控DCL、E2-UBC、ACA8、AP2等靶基因的調(diào)節(jié)。3、差異miRNA靶基因分析表明所涉及的8個(gè)miRNA家族的13個(gè)靶基因功能注釋為轉(zhuǎn)錄因子、F-Box生長(zhǎng)素應(yīng)答因子、DNA結(jié)合蛋白、Dicer酶、Ca 2+-ATP酶抑制劑、TIR-NBS-LRR抗病蛋白等,主要與人參自身能量代謝、生物脅迫及抗病免疫相關(guān)。
[Abstract]:Objective: To study the expression of ginseng small RNA, transcription and degradation group by high-throughput sequencing technique. The types of miRNAs contained in the two kinds of ginseng are identified, the differentially expressed miRNAs are screened, and the target gene prediction, the verification and the analysis are carried out. Methods:1. Trizol was used to extract the total RNA of the whole RNA. The quality, purity and integrity of total RNA were detected by agarose gel electrophoresis, TGem microspectrophotometer and Agilent 2100 Bioanalyzer. The Non-All-Uniene of two samples was obtained by short-sequence assembly and then compared with the database of Nr, KOG, KEGG, and GO to obtain the protein functional notes and the metabolic pathway information of Uniene. The differentially expressed mRNA between samples was selected for GO and KEGG enrichment analysis, and the differential mRNA was verified by Q-PCR. The obtained Clean threads are compared to the Rfam, Gene, Rebase, and the miRase 21.0 database to obtain known miRNAs and to count their amounts of expression; the sequence that is not above the pair is used for the prediction of the new miRNAs. In that invention, the difference miRNA between the sample is screened, the target gene prediction and the KEGG and GO enrichment analysis are carried out, and the differential miRNA is verified by using the Q-PCR.4, the database of the two degradation groups is constructed to obtain a database specific pair of the original sequence and the Rfam, Genbank, PolyN and the like to obtain the cDNA-sense. After the miRNA is known to perform the identification and analysis of the miRNA degradation site, the functional annotation information of the degradation gene is obtained by comparing with the Nr, the KEGG and the GO database. Results:1. The agarose gel electrophoresis showed that the 28S and 18S bands of the two total RNA were clear and intact, and the brightness was close to 2 times; the ratio of the OD260nm/ OD280nm to the Agilent 2100 Bioanalyzer was in the range of 1.8-2.1, the OD260nm/ OD280nm ratio was 2.0-2.3, the ratio of the OD260nm/ OD280nm was in the range of 2.0-2.3, RIN7.0, the RNA concentration of 900 ng/. mu.L, 28S/ 18S1.0, that is, the total RNA quality met the high-throughput building standard. Five thousand reads were obtained by sequencing two transcripts, and 63875 strips were obtained by De Novo. The functional notes to Uniene, with 19340 Uniene classified to 25 histone-line KOG categories, were classified as 64 GO functional categories, and 7268 Uniene was classified into 207 KEGG Pathways metabolic pathways. According to the gene expression difference multiple FoldChange and P-value,3216 difference mRNAs were selected, and 18 genes with significant difference were selected for Q-PCR verification, and the result of the verification was consistent with that of the transcript.3,235 known miRNAs were identified by small RNA high-throughput sequencing, and 246 known miRNAs were identified in the park. The fragment belongs to the 71 miRNA family, and the prediction of the new miRNA is not carried out by any database in comparison with the sequence, so that 39 new miRNAs are obtained. carrying out target gene prediction on the known and new miRNAs to obtain 169 target genes of 17 miRNA families, carrying out enrichment analysis on the target genes, Ubiquitin ligase and non-characteristic protein, etc. According to the difference of the expression of the miRNA and the P-value,13 different miRNAs were selected for Q-PCR, and the result of the validation was consistent with that of the miRNA. According to the clear and peak,812 and 724 shear sites were detected. According to the method for predicting the target gene of the bioinformatics miRNA, the degradation site of the cDNAsense is annotated to obtain 13 target genes of the 8 miRNA families, wherein the two target genes belonging to the new miRNA are obtained. Conclusion:1. The selection of microRNA contained in ginseng and garden ginseng is successfully constructed, and the microRNA in ginseng and garden is identified, and the microRNA is rich in fresh ginseng. The differential expression of ginseng miRNAs in different growth years is found, and the significant difference of the miRNA is mainly regulated by DCL, E2-UBC and ACA8. 3. The target gene analysis of the differentially expressed miRNA target showed that the 13 target gene functional notes of the 8 miRNA families involved were transcription factor, F-Box auxin response factor, DNA binding protein, Dicer enzyme, Ca 2 +-ATPase inhibitor, TIR-NBS-LRR disease-resistant protein, etc., mainly related to the energy metabolism of ginseng, And the biological stress and the disease-resistant immunity are related.
【學(xué)位授予單位】：廣東藥科大學(xué)
【學(xué)位級(jí)別】：碩士
【學(xué)位授予年份】：2017
【分類號(hào)】：S567.51
本文編號(hào)：2483781
[Abstract]:Objective: To study the expression of ginseng small RNA, transcription and degradation group by high-throughput sequencing technique. The types of miRNAs contained in the two kinds of ginseng are identified, the differentially expressed miRNAs are screened, and the target gene prediction, the verification and the analysis are carried out. Methods:1. Trizol was used to extract the total RNA of the whole RNA. The quality, purity and integrity of total RNA were detected by agarose gel electrophoresis, TGem microspectrophotometer and Agilent 2100 Bioanalyzer. The Non-All-Uniene of two samples was obtained by short-sequence assembly and then compared with the database of Nr, KOG, KEGG, and GO to obtain the protein functional notes and the metabolic pathway information of Uniene. The differentially expressed mRNA between samples was selected for GO and KEGG enrichment analysis, and the differential mRNA was verified by Q-PCR. The obtained Clean threads are compared to the Rfam, Gene, Rebase, and the miRase 21.0 database to obtain known miRNAs and to count their amounts of expression; the sequence that is not above the pair is used for the prediction of the new miRNAs. In that invention, the difference miRNA between the sample is screened, the target gene prediction and the KEGG and GO enrichment analysis are carried out, and the differential miRNA is verified by using the Q-PCR.4, the database of the two degradation groups is constructed to obtain a database specific pair of the original sequence and the Rfam, Genbank, PolyN and the like to obtain the cDNA-sense. After the miRNA is known to perform the identification and analysis of the miRNA degradation site, the functional annotation information of the degradation gene is obtained by comparing with the Nr, the KEGG and the GO database. Results:1. The agarose gel electrophoresis showed that the 28S and 18S bands of the two total RNA were clear and intact, and the brightness was close to 2 times; the ratio of the OD260nm/ OD280nm to the Agilent 2100 Bioanalyzer was in the range of 1.8-2.1, the OD260nm/ OD280nm ratio was 2.0-2.3, the ratio of the OD260nm/ OD280nm was in the range of 2.0-2.3, RIN7.0, the RNA concentration of 900 ng/. mu.L, 28S/ 18S1.0, that is, the total RNA quality met the high-throughput building standard. Five thousand reads were obtained by sequencing two transcripts, and 63875 strips were obtained by De Novo. The functional notes to Uniene, with 19340 Uniene classified to 25 histone-line KOG categories, were classified as 64 GO functional categories, and 7268 Uniene was classified into 207 KEGG Pathways metabolic pathways. According to the gene expression difference multiple FoldChange and P-value,3216 difference mRNAs were selected, and 18 genes with significant difference were selected for Q-PCR verification, and the result of the verification was consistent with that of the transcript.3,235 known miRNAs were identified by small RNA high-throughput sequencing, and 246 known miRNAs were identified in the park. The fragment belongs to the 71 miRNA family, and the prediction of the new miRNA is not carried out by any database in comparison with the sequence, so that 39 new miRNAs are obtained. carrying out target gene prediction on the known and new miRNAs to obtain 169 target genes of 17 miRNA families, carrying out enrichment analysis on the target genes, Ubiquitin ligase and non-characteristic protein, etc. According to the difference of the expression of the miRNA and the P-value,13 different miRNAs were selected for Q-PCR, and the result of the validation was consistent with that of the miRNA. According to the clear and peak,812 and 724 shear sites were detected. According to the method for predicting the target gene of the bioinformatics miRNA, the degradation site of the cDNAsense is annotated to obtain 13 target genes of the 8 miRNA families, wherein the two target genes belonging to the new miRNA are obtained. Conclusion:1. The selection of microRNA contained in ginseng and garden ginseng is successfully constructed, and the microRNA in ginseng and garden is identified, and the microRNA is rich in fresh ginseng. The differential expression of ginseng miRNAs in different growth years is found, and the significant difference of the miRNA is mainly regulated by DCL, E2-UBC and ACA8. 3. The target gene analysis of the differentially expressed miRNA target showed that the 13 target gene functional notes of the 8 miRNA families involved were transcription factor, F-Box auxin response factor, DNA binding protein, Dicer enzyme, Ca 2 +-ATPase inhibitor, TIR-NBS-LRR disease-resistant protein, etc., mainly related to the energy metabolism of ginseng, And the biological stress and the disease-resistant immunity are related.
【學(xué)位授予單位】：廣東藥科大學(xué)
【學(xué)位級(jí)別】：碩士
【學(xué)位授予年份】：2017
【分類號(hào)】：S567.51
【參考文獻(xiàn)】
相關(guān)期刊論文 前10條
1 王穎芳;陳艷琳;王文娟;;適用于轉(zhuǎn)錄組測(cè)序的人參根總RNA提取方法的篩選[J];廣東藥科大學(xué)學(xué)報(bào);2017年01期
2 吳亞運(yùn);張紹鵬;趙小龍;宋佳;陳平;;竹節(jié)參根狀莖總RNA提取方法的研究與優(yōu)化[J];武漢輕工大學(xué)學(xué)報(bào);2015年03期
3 張濤;韓梅;劉翠晶;胥苗苗;;人參總RNA提取方法及反轉(zhuǎn)錄酶的比較[J];江蘇農(nóng)業(yè)科學(xué);2015年08期
4 楊玖霞;張浩;王志龍;王旭麗;王國(guó)梁;;E3泛素連接酶調(diào)控植物抗病分子機(jī)理研究進(jìn)展[J];植物保護(hù);2015年04期
5 王杰;王全;田娜;王瑜;王愛(ài)香;張克中;崔金騰;;不同植物組織RNA提取方法的比較分析[J];北京農(nóng)學(xué)院學(xué)報(bào);2015年01期
6 方杰;胡毅翔;應(yīng)華忠;俞文英;金曉音;余陳歡;;米碎花藥效微小核糖核酸的鑒定及其靶基因功能分析[J];中草藥;2015年01期
7 張少會(huì);沈元月;;草莓果實(shí)cDNA文庫(kù)的構(gòu)建及鑒定[J];北京農(nóng)學(xué)院學(xué)報(bào);2014年02期
8 侯雙利;劉翠晶;楊利民;韓梅;;人參鯊烯環(huán)氧酶基因的克隆與序列分析[J];吉林農(nóng)業(yè)大學(xué)學(xué)報(bào);2014年02期
9 陳靜;高飛;周宜君;張孜宸;李章磊;曹玉震;張至瑋;;改良Trizol法提取高質(zhì)量蒙古沙冬青總RNA[J];生物技術(shù)通報(bào);2013年10期
10 劉笑夢(mèng);曹姣玲;臧玉柱;;人參皂甙對(duì)人白血病-60細(xì)胞凋亡的影響[J];中國(guó)當(dāng)代兒科雜志;2012年10期
相關(guān)碩士學(xué)位論文 前1條
1 肖新?lián)Q;甘蔗應(yīng)答黑穗病菌侵染的降解組分析及抗性相關(guān)miRNA靶標(biāo)挖掘[D];福建農(nóng)林大學(xué);2016年
，本文編號(hào)：2483781
本文鏈接：http://sikaile.net/kejilunwen/jiyingongcheng/2483781.html
最近更新
教材專著
