墨西哥灣扇貝轉(zhuǎn)錄組分析及其與扇貝“渤海紅”雜交F1代分子標(biāo)記研究
發(fā)布時間:2020-09-15 14:59
墨西哥灣扇貝(Argopecten irradians concentricus)是中國南方海域養(yǎng)殖的重要經(jīng)濟(jì)雙殼貝,1991年引進(jìn)至今,經(jīng)歷了20多年小群體繁育養(yǎng)殖,種質(zhì)退化嚴(yán)重。引進(jìn)新的遺傳資源,取得墨西哥灣扇貝種質(zhì)遺傳改良研究的突破性進(jìn)展是實現(xiàn)其養(yǎng)殖業(yè)健康可持續(xù)發(fā)展的根本。扇貝“渤海紅”為中國北方海域重點推廣養(yǎng)殖的新品種,課題組已于2015年9月從山東引進(jìn)扇貝“渤海紅”,并與本課題組選育的“墨西哥灣扇貝大閉殼肌新品系”進(jìn)行了雜交,小試獲得成功,其雜交可行性和形態(tài)雜種優(yōu)勢均已被證實,但相關(guān)分子遺傳學(xué)基礎(chǔ)還比較薄弱。因此,本研究圍繞墨西哥灣扇貝及其與扇貝“渤海紅”雜交F1代的分子遺傳學(xué)鑒定和評估,開展了以下研究:1、墨西哥灣扇貝轉(zhuǎn)錄組Illumina測序及生物信息學(xué)分析基于Illumina Hi Seq~(TM)2000平臺,對2組生長速率不同的5月齡墨西哥灣扇(普通養(yǎng)殖群體)進(jìn)行RNA-seq測序,獲得109 911條高質(zhì)量Unigenes,平均長度為1452 bp,其中50 433條(45.89%)Unigenes獲得注釋。根據(jù)FDR≤0.001且|log2Ratio|≥1的原則,篩選出9 901個差異表達(dá)基因,隨機選擇9個差異基因進(jìn)行q RT-PCR驗證,結(jié)果顯示,q PCR基因表達(dá)的整體趨勢與RNA-Seq測序結(jié)果呈線性相關(guān)(R2=0.9193),初步說明轉(zhuǎn)錄組測序?qū)τ诨虮磉_(dá)趨勢的預(yù)測是準(zhǔn)確的。對差異表達(dá)基因進(jìn)行KEGG代謝通路及GO功能富集分析,篩選到了TGF-β、MAPK信號轉(zhuǎn)導(dǎo)等幾個參與細(xì)胞增殖和分化的重要通路;差異基因在GO分類中的覆蓋面較廣,說明墨西哥灣扇貝的生長是生物學(xué)過程的全方位響應(yīng)。研究結(jié)果為墨西哥灣扇貝生長相關(guān)基因挖掘、生長發(fā)育機制揭示及分子標(biāo)記開發(fā)等提供了理論支撐。2、基于轉(zhuǎn)錄組測序的墨西哥灣扇貝生長相關(guān)基因及SNPs篩選生長相關(guān)基因是指在細(xì)胞分化及成熟或組織生長及重構(gòu)方面起促進(jìn)或抑制作用的一類基因。本研究基于Blast比對、GO和KEGG多種數(shù)據(jù)庫注釋,從墨西哥灣扇貝比較轉(zhuǎn)錄組差異表達(dá)基因中篩選出38個生長相關(guān)基因,涉及肌肉生長蛋白、消化酶和生長因子及其受體等18個類別。通過PCR擴(kuò)增直接測序法對墨西哥灣扇貝轉(zhuǎn)錄組Unigene23250和CL4903.Contig2序列上預(yù)測到的潛在SNP位點進(jìn)行驗證,共檢測到5個中度多態(tài)SNPs(0.25PIC0.50),包含4個轉(zhuǎn)換和1個顛換。CL4903.Contig2序列上的3個完全連鎖位點(A-6167-T、G-6178-A、G-6185-A)與Unigene23250序列上的G-1142-A位點發(fā)生非同義突變。對已鑒定的SNP與體質(zhì)量、閉殼肌重、殼長、殼寬和殼高的相關(guān)性進(jìn)行單因素方差分析,結(jié)果顯示,Unigene23250 G-1142-A位點AA基因型在體質(zhì)量、閉殼肌重及殼長、殼寬和殼高上的均值都顯著高于AG和GG基因型(P0.05),但AG和CG基因型在這5個性狀上的差異不顯著(P0.05)。CL4903.Contig2的G-6085-A位點AA基因型的體質(zhì)量及閉殼肌重均值顯著高于AG和GG基因型(P0.05),但這兩個位點不同基因型間形態(tài)性狀差異不顯著(P0.05)。推斷AA為優(yōu)勢基因型。研究結(jié)果表明,基于墨西哥灣扇貝轉(zhuǎn)錄組篩選到的這幾個SNP位點具有用于生長性狀分子輔助選育的潛力。3、墨西哥灣扇貝與扇貝“渤海紅”SSR開發(fā)及其雜交F1代驗證為對墨西哥灣扇貝與扇貝“渤海紅”雜交F1代進(jìn)行分子遺傳學(xué)鑒定和評估,本研究以墨西哥灣扇貝普通養(yǎng)殖群體40個個體對其轉(zhuǎn)錄組查找到的117個潛在SSR進(jìn)行篩選。結(jié)果有26個為多態(tài)位點(22.2%),其中21個(80.77%)在扇貝“渤海紅”中通用;多態(tài)位點主要為二堿基、6重復(fù),擴(kuò)增效率最高的是AT。隨機選擇8個多態(tài)SSR對扇貝“渤海紅”和墨西哥灣扇貝兩群體的遺傳多樣性進(jìn)行評估。結(jié)果顯示,墨西哥灣扇貝群體的遺傳多樣性低于扇貝“渤海紅”群體,兩個群體的遺傳距離為0.348,兩者雜交能產(chǎn)生雜種優(yōu)勢。開發(fā)3個扇貝“渤海紅”與墨西哥灣扇貝特異SSR,利用這3個特異SSR對40個扇貝“渤海紅”和墨西哥灣扇貝雜交F1代個體進(jìn)行驗證,3個特異SSR在F1代中的擴(kuò)增片段均不重疊,且同時存在來自墨西哥灣扇貝和扇貝“渤海紅”的特異等位基因,也未出現(xiàn)非親本之外的條帶,證明這些子代均為墨西哥灣扇貝和扇貝“渤海紅”的雜交子代。
【學(xué)位單位】:廣東海洋大學(xué)
【學(xué)位級別】:碩士
【學(xué)位年份】:2018
【中圖分類】:S917.4
【部分圖文】:
圖 2-1 RNA 提取步驟Fig.2-1 The procedure RNA extraction構(gòu)建、Illumina 測序和 De novo 組裝格后,2 組樣品各取 1 μg 用于轉(zhuǎn)錄組文庫構(gòu)建。具體方re Beads富集真核生物mRNA,采用高溫和金屬離子作用mRNA為模板,用六堿基隨機引物(random hexamers)合A 為模板合成二鏈 cDNA[78],再用 VAHTSTM DNAC(3)純化的雙鏈 cDNA 先進(jìn)行末端修復(fù)及加 A 尾操作,然TSTMDNA Clean Beads 進(jìn)行片段大小分選[79];(4)最后STM DNA Clean Beads 純化 PCR 產(chǎn)物,得到最終的文gilent 2100 Bioanalyzer 進(jìn)行質(zhì)檢。格后在 Illumina HiSeq 2000 平臺上進(jìn)行 6G 數(shù)據(jù)量測的原始數(shù)據(jù)(Raw reads),刪除接頭和低質(zhì)量 Reads, Trinity 軟件對高質(zhì)量 Clean reads 進(jìn)行 De novo 組裝[8gicl將其去冗余和進(jìn)一步拼接、同源轉(zhuǎn)錄本聚類,得到最
圖 2-2 De novo 組裝流程Fig.2-2 The process of De novo assemblynigene 功能注釋組裝得到的 Unigene 在六大數(shù)據(jù)庫進(jìn)行比對注釋,并分別對注釋到每個注釋上的 Unigene 數(shù)目進(jìn)行統(tǒng)計。六大數(shù)據(jù)庫分別為:非冗余蛋白序列、非冗余核酸序列數(shù)據(jù)庫 (Nt)、高質(zhì)量非冗余的蛋白數(shù)據(jù)庫 (Swiss-P理基因組、生物通路、疾病、藥物和化學(xué)物質(zhì)之間聯(lián)系的集成數(shù)據(jù)庫(K源家族蛋白數(shù)據(jù)庫(COG)和 GO 數(shù)據(jù)庫。差異表達(dá)基因篩選過 FPKM 計算法[82]對 Unigene 表達(dá)量進(jìn)行標(biāo)準(zhǔn)化,然后使用 DEGseq 軟樣本間相同基因的表達(dá)量進(jìn)行差異分析,并以 FDR(false discovery rate2Fold≥1(快速生長組/慢速生長組)為閾值,篩選差異表達(dá)基因(Differesed unigene)。FPKM[82]是目前計算基因表達(dá)量最常用的方法,其是指每中來自某一基因每千堿基長度的 Fragment 數(shù)目。其計算公式為: =106 3公式 (2-1)
LabchipGX檢測膠圖
【學(xué)位單位】:廣東海洋大學(xué)
【學(xué)位級別】:碩士
【學(xué)位年份】:2018
【中圖分類】:S917.4
【部分圖文】:
圖 2-1 RNA 提取步驟Fig.2-1 The procedure RNA extraction構(gòu)建、Illumina 測序和 De novo 組裝格后,2 組樣品各取 1 μg 用于轉(zhuǎn)錄組文庫構(gòu)建。具體方re Beads富集真核生物mRNA,采用高溫和金屬離子作用mRNA為模板,用六堿基隨機引物(random hexamers)合A 為模板合成二鏈 cDNA[78],再用 VAHTSTM DNAC(3)純化的雙鏈 cDNA 先進(jìn)行末端修復(fù)及加 A 尾操作,然TSTMDNA Clean Beads 進(jìn)行片段大小分選[79];(4)最后STM DNA Clean Beads 純化 PCR 產(chǎn)物,得到最終的文gilent 2100 Bioanalyzer 進(jìn)行質(zhì)檢。格后在 Illumina HiSeq 2000 平臺上進(jìn)行 6G 數(shù)據(jù)量測的原始數(shù)據(jù)(Raw reads),刪除接頭和低質(zhì)量 Reads, Trinity 軟件對高質(zhì)量 Clean reads 進(jìn)行 De novo 組裝[8gicl將其去冗余和進(jìn)一步拼接、同源轉(zhuǎn)錄本聚類,得到最
圖 2-2 De novo 組裝流程Fig.2-2 The process of De novo assemblynigene 功能注釋組裝得到的 Unigene 在六大數(shù)據(jù)庫進(jìn)行比對注釋,并分別對注釋到每個注釋上的 Unigene 數(shù)目進(jìn)行統(tǒng)計。六大數(shù)據(jù)庫分別為:非冗余蛋白序列、非冗余核酸序列數(shù)據(jù)庫 (Nt)、高質(zhì)量非冗余的蛋白數(shù)據(jù)庫 (Swiss-P理基因組、生物通路、疾病、藥物和化學(xué)物質(zhì)之間聯(lián)系的集成數(shù)據(jù)庫(K源家族蛋白數(shù)據(jù)庫(COG)和 GO 數(shù)據(jù)庫。差異表達(dá)基因篩選過 FPKM 計算法[82]對 Unigene 表達(dá)量進(jìn)行標(biāo)準(zhǔn)化,然后使用 DEGseq 軟樣本間相同基因的表達(dá)量進(jìn)行差異分析,并以 FDR(false discovery rate2Fold≥1(快速生長組/慢速生長組)為閾值,篩選差異表達(dá)基因(Differesed unigene)。FPKM[82]是目前計算基因表達(dá)量最常用的方法,其是指每中來自某一基因每千堿基長度的 Fragment 數(shù)目。其計算公式為: =106 3公式 (2-1)
LabchipGX檢測膠圖
【參考文獻(xiàn)】
相關(guān)期刊論文 前10條
1 王忠良;丁q
本文編號:2819115
本文鏈接:http://www.sikaile.net/nykjlw/scyylw/2819115.html
最近更新
教材專著