1、提取方法

1)植物葉片：（天根 DP441（廠(chǎng)家：天根，型號：DP441）

2)植物根、莖、種子：艾德萊 RN40（廠(chǎng)家：艾德萊，型號：RN40））試劑盒

3)動(dòng)物常規組織：TRIzol

4)皮膚、毛發(fā)：凱捷 217044（廠(chǎng)家：凱捷。型號：217004）

5)全血PAX管送樣：PAXgene Blood miRNA Kit(50)（廠(chǎng)家：GENENODE。型號：2145A）試劑盒

6)EDTA抗凝血：TRIzol

2、檢測方法

2.1檢測方法

使用 Nanodrop2000 （廠(chǎng)家：賽默飛，型號： Nanodrop2000）對于提取的核酸進(jìn)行濃度純度檢測，并使用LabChip GX（廠(chǎng)家：鉑金埃爾默，型號：LabChip GX Touch HT）對完整性進(jìn)行檢測

2.2判定標準

1)總量≥1（ug），滿(mǎn)足 2 次建庫

2)濃度≥20（ng/ul）

3)體積≥10(ul)

4)OD260/280 在 1.7-2.5 之間，OD260/230 在 0.5-2.5 之間，260nm 吸收峰顯示正常

5)RIN 值≥4.0（非植物），RIN 值≥4.0（植物），同時(shí)若 GX 檢測 6.0-4.0 的需結合基線(xiàn)判斷，28S/18S≥1 且基線(xiàn)無(wú)上抬或輕微上抬

6)樣本狀態(tài)正常，樣本粘稠、顏色異常、有渾濁或不溶物均不合格

3、建庫方法

3.1rRNA去除

1) rRNA與探針結合

A.將對應物種的去rRNA探針混合物加入到準備好的總 RNA 樣品中，混勻離心，68℃孵育確保RNA變性

B.孵育結束后，將PCR管置于室溫，孵育2 min

2)rRNA探針復合物的捕獲與去除

A.吸取 RRB（rRNA Removal Beads）至新1.5 mL Nuclease-free離心管中

B.將上步中的RNA/Probe混合物轉入RRB中并立即迅速混勻，然后再渦旋混勻，室溫孵育5 min

C.孵育結束后再中速渦旋，然后置于磁力架上2 min至上清液澄清

D.轉移上清液到新離心管中，補加Nuclease-free Water，置于磁力架上至上清液澄清，轉移上清液到新的離心管中

3) rRNA-depleted RNA的純化及打斷

A.向上述反應產(chǎn)物加入相應體積的（1.8X）已混勻的磁珠，混勻，室溫孵育而后置于磁力架上，待磁珠吸附到磁力架一側后移棄上清液

B.向反應管中加入新鮮配制的80% ethanol，磁力架上靜置30s，移棄上清液，重復該步驟1次

C.打開(kāi)管蓋，空氣干燥至磁珠表面無(wú)光澤，有細微裂紋，加入Frag/Prime Buffer洗脫，吹吸混勻，室溫靜置，磁力架上靜置，轉移至新的Nuclease-free 0.2 mL離心管中，PCR儀高溫打斷，打斷完成后立即放在冰上

3.2合成 cDNA 第一條鏈

向上步反應的離心管中加入反轉錄合成所需的試劑1st Strand Buffer 2、1st Strand Enzyme Mix 2及Actinomycin D (5mg/ml)，混勻離心后放入PCR儀中設定好程序進(jìn)行cDNA第一條鏈的合成。

3.3合成 cDNA 第二條鏈（含末端修復及接頭連接）

1)向合成的 cDNA 第一條鏈產(chǎn)物中依次加入二鏈合成反應的試劑2nd Strand Buffer 2（with dUTP）、2nd Strand Enzyme Super Mix，混勻離心后，放在PCR儀中孵育完成二鏈的合成

2)運行結束后加入末端修復及接頭連接的試劑：Rapid Ligation Buffer 3、Rapid DNA Ligase 2及IDT接頭，混勻離心后，放在PCR儀中運行程序

3.4末端修復及 3’-末端加A

向上述純化產(chǎn)物中加入末修試劑，放入PCR儀中進(jìn)行反應，反應結束后，立即進(jìn)行接頭連接

3.5連接接頭及產(chǎn)物純化

1)向上述反應產(chǎn)物加入0.9X已混勻的磁珠，混勻，室溫孵育而后置于磁力架上，待磁珠吸附到磁力架一側后移棄上清液

2)離心管置于磁力架上，加新鮮配制的 80% ethanol，靜置 30s，移棄上清液

3)重復步驟 2 一次，第二次清洗后離心用 10ul 槍頭將上清液去除干凈

4)離心管置于磁力架上，打開(kāi)管蓋，空氣干燥

5)取下離心管，加入Nuclease-free Water，混勻，室溫孵育，而后置于磁力架上靜置，吸取上清液轉入新的 離心管中

3.6PCR 擴增

向上述離心管中加入 PCR 反應所需的各種試劑，混勻后離心，根據 PCR 反應條件，放入 PCR 儀中進(jìn)行 PCR 擴增

注意：primer中包含Index序列，使用時(shí)需嚴格按照上機調度安排的序列添加

3.7PCR 產(chǎn)物純化

1)瞬時(shí)離心 PCR 反應管，加入相應體積的磁珠，充分混勻

2)室溫孵育，然后置于磁力架上靜置，小心移除上清液

3)向反應管中加入新鮮配制的 80% ethanol ，磁力架上靜置 30s ，移棄上清液，重復該步驟 1 次

4)打開(kāi)管蓋，磁珠干燥后加入 Nuclease-free Water 洗脫 DNA ，吹吸混勻，室溫靜置，磁力架上靜置，小心 吸取上清至新的離心管中，并分裝 1.5 μL 于新的 1.5 mL 離心管中用于片段檢測，-20℃保存

3.8文庫質(zhì)檢

文庫通過(guò) Qsep-400 方法進(jìn)行質(zhì)檢，滿(mǎn)足如下指標即可上機檢測：

3.9引物接頭序列

adpter3=”AGATCGGAAGAGCACACGTCTGAACTCCAGTCAC”;

adpter5=”AGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGT”;

3.10建庫試劑耗材

Ribo-off rRNA Depletion Kit(Human/Mouse/Rat)，貨號N406-01

Ribo-off rRNA Depletion Kit(Plant)，貨號N409-02

1)文庫構建試劑盒

Ribo-off rRNA Depletion Kit(Human/Mouse/Rat)，貨號N406-01

Ribo-off rRNA Depletion Kit(Plant)，貨號N409-02

2)產(chǎn)物純化磁珠

VAHTSTM DNA Clean Beads，貨號N411-03

VAHTSTM RNA Clean Beads，貨號N412-02

3)質(zhì)檢使用試劑：Qsep400 標準DNA 卡夾

4)定量使用試劑：Qubit TM dsDNA HS Assay Ki

5)1.5ml離心管、0.2ml PCR管、10ul，200ul槍頭：Axygen

3.11建庫設備

4、測序方法

測序設備：NovaSeq X plus

測序試劑盒：NovaSeqTM X? Series25B Rgt Kit

英文名:Long Non-coding RNA Derived from lncRNA– mRNA Co-expression Networks Modulates the Locust Phase Change

期刊：Genomics Proteomics Bioinformatics

IF：7.051（1區）

通訊作者單位：中科院、河北大學(xué)

研究背景

長(cháng)鏈非編碼RNA（lncRNA）調節基因表達、動(dòng)物行為等各種生物學(xué)過(guò)程。盡管蛋白質(zhì)編碼基因、microRNA和神經(jīng)肽在亞洲飛蝗的表型可塑性調節中起著(zhù)重要作用，但有關(guān)lncRNA在此過(guò)程中功能研究較少。本文應用高通量RNA-seq來(lái)比較蝗蟲(chóng)型變時(shí)程中lncRNA和mRNA的表達模式。結果顯示lncRNA在型變的早期階段反應更快。功能注釋表明，早期改變的lncRNA在分散和群居階段采用了不同的途徑來(lái)應對種群密度的變化。篩選了分散和群居階段的網(wǎng)絡(luò )中兩個(gè)重疊的中樞lncRNA基因座進(jìn)行功能驗證。本文進(jìn)一步證明LNC1010057為潛在的蝗蟲(chóng)型變因子。這項工作為深入了解蝗蟲(chóng)型變的分子機制并擴大lncRNA在動(dòng)物行為中的作用范圍提供了重要的數據。

材料方法

實(shí)驗材料：散居化處理為將群居型蝗蟲(chóng)單獨飼養0h、4h、8h和16h后取出手機蝗蟲(chóng)大腦進(jìn)行測序；

群居化處理為將10只散居型蝗蟲(chóng)與20只群居型蝗蟲(chóng)飼養于一個(gè)小籠子（10厘米×10厘米×10厘米）中，于0h、4h、8h和16h后取出，收集蝗蟲(chóng)大腦進(jìn)行測序。在同一時(shí)間點(diǎn)采集樣品的三個(gè)生物學(xué)重復，提取RNA，建立cDNA文庫并進(jìn)行RNA-seq。對lncRNA和mRNA進(jìn)行差異表達分析及qRT-PCR驗證，對不同時(shí)間點(diǎn)的樣本進(jìn)行STEM(ShortTime-seriesExpressionMiner)分析，并構建lncRNA–mRNA共表達網(wǎng)絡(luò )的構建。對目標中心節點(diǎn)lncRNA進(jìn)行RNAi，采集其行為錄像數據進(jìn)行行為數據分析。

實(shí)驗結果

1、蝗蟲(chóng)lncRNA的外顯子比mRNA少但更長(cháng)

研究發(fā)現大約78％的lncRNA包含2個(gè)外顯子，而mRNA中包含的外顯子數量為1至120（圖1C）。因此，lncRNA的外顯子明顯長(cháng)于mRNA的外顯子（平均長(cháng)度1142bpvs.263bp，圖1D，左）。同時(shí)，lncRNA的內含子明顯短于mRNA的內含子（平均長(cháng)度：10086bpvs.12442bp，圖1D，右）。表達水平分析表明，lncRNA的總體表達水平明顯低于mRNA的表達水平（平均值為0.6vs.1.9；圖1E）。但是，表達特異性分析表明，lncRNA的表達受時(shí)間限制的程度要高于mRNA（平均值為0.672對0.436圖1F）?；谙鄬τ趍RNA的lncRNA基因組位置，蝗蟲(chóng)lncRNA被分類(lèi)為基因間、重疊區、有義內含子、反義內含子、有義外顯子和反義外顯子?；认x(chóng)77％以上的lncRNA是長(cháng)基因間的ncRNA（lincRNA，圖1G）。這些結果表明，蝗蟲(chóng)lncRNA與mRNA在結構和表達上有很大不同?；认x(chóng)lncRNAs更長(cháng)，具有更少但更長(cháng)的外顯子和更短的內含子。此外，lncRNA的表達模式顯示出比mRNA更高的時(shí)間特異性。

圖1lncRNA和mRNA之間的不同結構和表達

2、群居型蝗蟲(chóng)中特異性表達的lncRNA比在散居型蝗蟲(chóng)中表達的更多

在散居型和群居型蝗蟲(chóng)大腦中共表達9722個(gè)lncRNA（73.9％），而散居型蝗蟲(chóng)中特異性表達962個(gè)lncRNA（7.3％），群居型蝗蟲(chóng)中特定表達2469個(gè)lncRNA（18.8％）（圖2A，頂部）。分別在散居型和群居型蝗蟲(chóng)中特異性表達了479和1090個(gè)mRNA（3.1％和7.1％）（圖2A，底部）。這些結果表明，與mRNA相比，在兩個(gè)蝗蟲(chóng)相型特異性表達的lncRNA的百分比更高，而在群居蝗蟲(chóng)中比散居型蝗蟲(chóng)中表達的lncRNAs更多。與散居型蝗蟲(chóng)相比，群居型蝗蟲(chóng)中335個(gè)lncRNA和779個(gè)mRNA的表達水平下調，而313個(gè)lncRNA和261個(gè)mRNA的表達上調[倍數變化（FC）>2和P<0.05；圖2B]。lincRNA在下調和上調的lncRNA中所占的比例高（分別為81.2％和87.8％），其次分別是反義外顯子和有義內含子lncRNA（圖2B）。在表達中按FC排名的前10個(gè)lncRNA基因顯示在圖2C中。這表明散居型和群居型蝗蟲(chóng)大腦中表達的lncRNA的數量及其表達水平明顯不同。

圖2lncRNA在蝗蟲(chóng)型變中顯示出不同的表達變化模式

3、lncRNA對種群密度變化的快速應答

為了進(jìn)一步分析lncRNA和mRNA的表達模式，使用了短時(shí)間序列表達挖掘器（ShortTime-seriesExpressionMiner，STEM）對基于表達的轉錄本進(jìn)行聚類(lèi)。在CS期間，將214個(gè)lncRNA和242個(gè)mRNA分別分為8個(gè)和9個(gè)重要的表達譜（圖3D）。其中lncRNA表達譜15以及mRNA表達譜15、13和16直到群居化處理后8或16小時(shí)才改變。這些配置文件稱(chēng)為后期更改的配置文件（模式d）。與后期更改的配置文件相反，早期更改的配置文件以4小時(shí)時(shí)間點(diǎn)開(kāi)始的轉錄表達變化為特征。根據4小時(shí)后的表達變化，將早期變化的輪廓細分為模式a（早期變化）、模式b（早期中變化）和模式c（可持續變化）。在IG過(guò)程中，269個(gè)lncRNA和639個(gè)mRNA分別聚集成6個(gè)和7個(gè)顯著(zhù)表達譜。所有lncRNA配置文件均為早期變化（圖3E）。在mRNA表達譜中，有6個(gè)是早期改變的，而譜12是晚期改變的。在CS和IG期間，lncRNA的聚類(lèi)分布圖的數量與mRNA的分布無(wú)明顯差異。但是，通過(guò)計算譜圖中的轉錄本數量，發(fā)現早期改變的lncRNA的百分比比CS中的mRNA的百分比高61.1％（88.3％vs.54.8％，圖3F）。同樣，IG中早期改變的lncRNA的百分比高于IG中的mRNA（100％vs.81.3％，圖3F）。因此，lncRNA比mRNA對散居化和群居化處理的反應更快。在CS中，早期改變的lncRNA的表達變化速率快于285個(gè)mRNA的表達變化速率（0.55vs.0.24）和IG（0.77vs.0.39）（圖3G）。因此，早期改變的lncRNA的比例和表達變化率高于早期改變的mRNA。因此，蝗蟲(chóng)lncRNA比mRNA對種群密度的變化更敏感。

圖3lncRNA對散居化和群居化處理的快速反應

4、lncRNA參與CS和IG早期變化的不同途徑

在CS中，自噬調節、剪接體和肌醇磷酸代謝途徑在上位和至少兩個(gè)模塊中均過(guò)代表（圖4C）。因此，CS中早期改變的lncRNA可能在調節自噬、RNA剪接和信號轉導途徑中發(fā)揮重要作用。lncRNA對自噬的調控參與了群居化的初期和中期，而對剪接體和磷酸肌醇代謝的調控則參與了整個(gè)過(guò)程。與神經(jīng)遞質(zhì)釋放有關(guān)的突觸小泡循環(huán)途徑的lncRNA僅在模式a（早期改變）模塊中被過(guò)代表（圖4C）。該結果表明，早期改變的lncRNA可能在CS期間調節神經(jīng)系統中具有重要作用。此外，一些lncRNA與已知的型變相關(guān)基因密切相關(guān)。例如，早期改變的lncRNALNC531328.2、LNC1425451.2、LNC1088763.7和LNC492755.1與基因NPYR、NPF1a和Vat1相關(guān)（圖4A）。持續變化的lncRNALNC494161.1和LNC1065048.14與多巴胺途徑中的基因Ebony和Vat1相關(guān)。在CS網(wǎng)絡(luò )中，程度值前5％的lncRNA被視為hublncRNA（圖4A）。計算了lncRNA基因座的程度值，在CS中，基于度數排名前5％的10個(gè)lncRNA基因座被鑒定為hublncRNA。其中，四個(gè)基因表達上調，其他六個(gè)基因表達下調（圖4D）。在IG期間，包括谷氨酸能突觸、多巴胺能突觸和膽堿能突觸途徑在內的與突觸有關(guān)的途徑得以豐富。這些途徑中的大多數都富含早期改變的模塊。同時(shí)，多巴胺能突觸和膽堿能突觸途徑僅參與模式a（早期改變）的模塊（圖4E）。IG中也豐富了信號轉導途徑，例如鈣信號傳導、Ras信號傳導、胰島素信號傳導和MAPK信號傳導途徑。功能注釋的結果表明，與CS相比，IG中早期變化的lncRNA參與的突觸相關(guān)和信號處理途徑更多。

圖4早期改變的lncRNA參與CS和IG的不同途徑

5、LNC1010057可能調節蝗蟲(chóng)型變

為為了驗證LNC1010057和LNC992414在蝗蟲(chóng)型變中的功能，進(jìn)行了一系列分子生物學(xué)實(shí)驗。首先，克隆了在LNC1010057和LNC992414基因座中鑒定出的長(cháng)的轉錄本。LNC1010057由重復的A、B和C元素組成，這些元素依次分布并重復4次，但C元素重復三次（圖5A）。序列比對證明LNC992414在5’端與LNC1010057共享相似的重復序列，但是它們是從不同的基因組基因座轉錄而來(lái)的（圖5A）。盡管LNC992414的定量表達水平可以通過(guò)特異性引物檢測，但LNC1010057的表達水平為重復元件的總表達水平。如預期的那樣，LNC1010057和LNC992414的表達模式極為相似。在CS期間，兩種lncRNA的表達持續增加，并且在16h幾乎增加了四倍（圖5B）。在IG期間，4h后表達水平顯著(zhù)下降，之后保持相對穩定（圖5B）。LNC1010057和LNC992414之間的序列結構和表達模式的相似性表明它們是同源的lncRNA。但是，實(shí)時(shí)qPCR分析表明，在散居型和群居型蝗蟲(chóng)的大腦中，LNC1010057的表達水平比LNC992414的表達水平高約1000倍（圖5C）。其次，為了測試LNC1010057和LNC992414是否參與蝗蟲(chóng)型變調節，在通過(guò)RNAi抑制蝗蟲(chóng)大腦中的表達后進(jìn)行了行為分析。顯著(zhù)降低LNC1010057的表達水平（7470對3764610）和LNC992414（1.9vs.1.0）（圖5D）后行為分析表明，群居型蝗蟲(chóng)其行為顯著(zhù)改變?yōu)樯⒕訝顟B(tài)（圖5E）。此外，多個(gè)與型有關(guān)的行為參數發(fā)生了變化?？傄苿?dòng)距離（TDM）和總移動(dòng)持續時(shí)間（TDMV）顯著(zhù)降低（圖5F），但是，移動(dòng)速度沒(méi)有區別。同時(shí)群居蝗蟲(chóng)的特定喜好行為顯著(zhù)下降（58.3對-13.5圖5F）。但LNC992414表達降低并未引起從群居狀態(tài)到散居狀態(tài)的轉變（圖5G和H）。與對照相比，與型相關(guān)的行為參數（包括運動(dòng)和特定物種的優(yōu)選行為）沒(méi)有差異（圖5I）。LNC992414的RNAi實(shí)驗不會(huì )引起行為變化，因此排除了LNC992414調節蝗蟲(chóng)型變的可能性。這些結果表明是LNC1010057而不是LNC992414潛在地調節了蝗蟲(chóng)的型變。

圖5LNC1010057潛在地調節了蝗蟲(chóng)的型變

總結

lncRNA被證實(shí)是生物過(guò)程的關(guān)鍵調控因子，調節包括mRNA轉錄、穩定性、翻譯和翻譯后修飾等多種生物途徑。之前研究表明昆蟲(chóng)lncRNA參與了例如殺蟲(chóng)劑抗性、繁殖力和腺體凋亡等生物過(guò)程，但是尚未有證據證明lncRNA可以調節非模型昆蟲(chóng)的行為?；认x(chóng)是世界范圍內的一種農業(yè)害蟲(chóng)，表現出顯著(zhù)的表型可塑性。為了鑒定與其型變相關(guān)的lncRNA，本文系統地分析了蝗蟲(chóng)lncRNA的表達并注釋其功能，證實(shí)與mRNAs相比lncRNAs顯示對型變更敏感的響應，并證實(shí)了其中一個(gè)lncRNA可調節型變相關(guān)行為。本研究揭示了lncRNA在蝗蟲(chóng)型變中的重要作用以及表型中蛋白編碼基因和lncRNA之間的相互作用。深入解析蝗蟲(chóng)散居型和群居型轉變的分子機制，為可持續治理蝗蟲(chóng)的新策略和新方法的開(kāi)發(fā)提供了基礎。