Comprehensive genomic and phenotypic analyses reveal the genetic basis of fruit quality in litchi

研究背景

荔枝是東南亞一種重要的經(jīng)濟(jì)水果作物，以其獨(dú)特的風(fēng)味和營(yíng)養(yǎng)價(jià)值而聞名。荔枝起源于中國(guó)云南省，具有2000多年的栽培史，由于長(zhǎng)期的自然選擇和在不同氣候及種植條件下的人工育種，孕育了豐富的種質(zhì)資源，可用于性狀持續(xù)的遺傳改良。至今已有許多工作是從表型性狀來研究和評(píng)估荔枝的多樣性，然而，相應(yīng)的遺傳變異未被充分的挖掘利用，尤其是能夠代表大量遺傳信息的核心種質(zhì)，群體遺傳結(jié)構(gòu)不清晰，許多重要育種性狀如果實(shí)大小、風(fēng)味、種子大小等調(diào)控果實(shí)品質(zhì)的遺傳基礎(chǔ)仍有待深入探索。

研究?jī)?nèi)容

該研究中對(duì)來自國(guó)家荔枝種質(zhì)資源圃（廣州）從全球收集的276份具有代表性荔枝種質(zhì)資源進(jìn)行了重測(cè)序，測(cè)序深度20X，并鑒定出3450萬(wàn)個(gè)高質(zhì)量的雙等位基因SNP，構(gòu)成了迄今為止荔枝（乃至整個(gè)無(wú)患子科）最全面的遺傳變異圖譜。

隨后，分析了群體結(jié)構(gòu)，在基因組水平上揭示了群體之間的遺傳變異水平和親緣關(guān)系，分析表明該群體內(nèi)存在四個(gè)主要亞群：YNG（云南亞群，來自云南野生、廣西大興的種質(zhì)）、HNG（海南亞群，來自海南和廣地博白的種質(zhì)）、FGG1, 和FGG2亞群（包含來自福建和部分廣東、廣西的種質(zhì)）。YNG的核苷酸多樣性低于HNG、FGG1和FGG2的，LD遺傳距離的大于其他三個(gè)亞群(圖1)。

圖1-荔枝群體結(jié)構(gòu)和遺傳多樣性

通過對(duì)21個(gè)果實(shí)品質(zhì)相關(guān)性狀進(jìn)行表型分析，發(fā)現(xiàn)果實(shí)性狀在許多方面表現(xiàn)出極大的多樣性，例如果實(shí)大小、種子大小、糖酸含量等方面。進(jìn)一步并結(jié)合全基因組關(guān)聯(lián)研究，鑒定了總計(jì)460個(gè)與果實(shí)性狀顯著相關(guān)的SNP和1,807個(gè)候選基因，其中包含大量負(fù)責(zé)調(diào)控種子和果實(shí)品質(zhì)性狀的候選基因和基因組位點(diǎn)，并篩選出了幾個(gè)與種子早期發(fā)育相關(guān)的候選基因。特別是，利用高效液相色譜（HPLC）方法測(cè)定了190個(gè)荔枝品種果實(shí)的糖組分和含量，隨后，以蔗糖和葡萄糖的含量作為性狀進(jìn)行GWAS分析，在7號(hào)染色體上檢測(cè)到一個(gè)強(qiáng)關(guān)聯(lián)峰這個(gè)關(guān)聯(lián)峰位于基因LITCHI009111（LcSAI）的基因體區(qū)域內(nèi)，該基因編碼β-果糖呋喃苷糖酶，也稱為轉(zhuǎn)化酶（invertase），該酶催化蔗糖分解為葡萄糖和果糖，并通過原核表達(dá)和轉(zhuǎn)基因等方法驗(yàn)證了該基因的功能。通過比較LcSAI蔗糖和還原糖種質(zhì)中的差異，開發(fā)了一個(gè)與糖組分和含量緊密連鎖的分子標(biāo)記，可用于荔枝育種的早期篩選，縮短育種周期。

圖2-荔枝果實(shí)性狀表型的多樣性

圖3-LcSAI 在荔枝果實(shí)的糖分積累中的作用

研究總結(jié)

綜上所述，該研究揭示了荔枝的遺傳多樣性和群體結(jié)構(gòu)，通過全面的基因組分析和表型分析，闡明了果實(shí)品質(zhì)的遺傳基礎(chǔ)。這些發(fā)現(xiàn)為理解荔枝果實(shí)品質(zhì)的遺傳基礎(chǔ)提供了寶貴資源和知識(shí)，為進(jìn)一步的遺傳改良貢獻(xiàn)了理論基礎(chǔ)。

以上內(nèi)容來源于網(wǎng)絡(luò)，侵刪

辣椒（Capsicum annuum）是全球最早且廣泛栽培的蔬菜作物之一，其產(chǎn)生的包括辣椒紅素和辣椒素在內(nèi)的多種次生代謝產(chǎn)物在食品，醫(yī)療和軍工業(yè)都具有重要作用。一直以來，對(duì)于辣椒的改良主要依賴于傳統(tǒng)育種方法。相較于同為茄科作物的番茄，辣椒具有龐大而復(fù)雜的基因組，且其中超過80%為轉(zhuǎn)座子或重復(fù)序列。研究表明，轉(zhuǎn)座子和重復(fù)序列是轉(zhuǎn)錄調(diào)控元件的重要來源，由此暗示了辣椒基因組中含有豐富的順勢(shì)調(diào)控元件以及多樣的轉(zhuǎn)錄調(diào)控方式，而調(diào)控元件的變化是動(dòng)植物進(jìn)化和馴化的主要驅(qū)動(dòng)力之一。

該研究以樟樹港辣椒為模型，系統(tǒng)揭示了辣椒在全基因組范圍內(nèi)的表觀特征，并在此基礎(chǔ)上推測(cè)包括啟動(dòng)子和增強(qiáng)子在內(nèi)的瞬時(shí)調(diào)控元件。研究發(fā)現(xiàn)，辣椒啟動(dòng)子主要由H3K4me3,?H3K27ac和開放染色質(zhì)共同修飾，而增強(qiáng)子的特征為開放染色質(zhì)。在此基礎(chǔ)上，該研究對(duì)生長(zhǎng)發(fā)育以及脅迫響應(yīng)的調(diào)控元件進(jìn)行了挖掘，并對(duì)部分元件進(jìn)行實(shí)驗(yàn)驗(yàn)證。此外，該研究報(bào)道了辣椒基因組基因間區(qū)H3K4me1和H3K27me3共同修飾的區(qū)域，推測(cè)為暫時(shí)未表達(dá)（poised）的調(diào)控元件，其具體分子功能有待驗(yàn)證。

綜上，該研究系統(tǒng)挖掘辣椒基因組中重要的順勢(shì)調(diào)控元件，以及其調(diào)控的靶標(biāo)基因，不僅為基因調(diào)控研究提供了基礎(chǔ)，也為通過靶向遺傳和表觀遺傳操作改良辣椒品種提供了理論支持。

以上內(nèi)容來源于北京大學(xué)現(xiàn)代農(nóng)業(yè)研究院，侵刪

百邁客生物為該研究提供了SLAF-seq測(cè)序分析服務(wù)。

研究背景

油茶（Camellia oleifera Abel.）是我國(guó)重要的木本油料作物之一，以其高含油量與豐富的不飽和脂肪酸而享有“東方橄欖油”的美譽(yù)。我國(guó)已有超過2000年的油茶栽培歷史，并具有約4500萬(wàn)畝的種植面積。油茶不僅具有經(jīng)濟(jì)價(jià)值，更具重要生態(tài)功能。其廣泛分布于我國(guó)南方多個(gè)氣候生態(tài)帶，是研究森林生態(tài)適應(yīng)性與基因資源保護(hù)的重要對(duì)象。然而，當(dāng)前油茶品種識(shí)別主要依賴形態(tài)學(xué)方法，容易造成品種混淆、資源管理混亂和知識(shí)產(chǎn)權(quán)糾紛等一系列問題，難以滿足現(xiàn)代林業(yè)種質(zhì)資源保護(hù)與良種選育的需求。為推動(dòng)油茶良種選育，產(chǎn)業(yè)升級(jí)和響應(yīng)國(guó)家森林四庫(kù)相關(guān)文件的精神，需加快構(gòu)建高效、廉價(jià)且準(zhǔn)確的DNA指紋圖譜。

研究結(jié)果

1.用于SLAF測(cè)序的最優(yōu)限制性內(nèi)切酶組合

首先，研究共采集25個(gè)油茶栽培品種，共100個(gè)樣本,利用高通量SLAF-seq建庫(kù)策略進(jìn)行數(shù)據(jù)獲取。在建庫(kù)酶切方案篩選方面，團(tuán)隊(duì)預(yù)先開展酶切模擬試驗(yàn)比較篩選最優(yōu)組合，F(xiàn)igure 1 顯示，RsaI + HaeIII組合在標(biāo)簽分布均勻性、多態(tài)性與文庫(kù)復(fù)雜度上表現(xiàn)最佳。共獲得362,152個(gè)高質(zhì)量SLAF 標(biāo)簽，其中254,357個(gè)為多態(tài)性標(biāo)簽，平均GC含量41.30%，平均Q30高達(dá)93.44%，顯示測(cè)序數(shù)據(jù)具備極高的可靠性。RsaI + HaeIII雙酶切組合可實(shí)現(xiàn)高效的特異性位點(diǎn)擴(kuò)增與標(biāo)簽標(biāo)準(zhǔn)化，為后續(xù)大規(guī)模SNP檢測(cè)奠定了數(shù)據(jù)基礎(chǔ)。

2.基因組變異特征分析

經(jīng)測(cè)序，每個(gè)樣本平均測(cè)序深度為13.95×，共鑒定出470,397個(gè)多態(tài)性SLAF標(biāo)簽，并初步檢測(cè)到36,317,329個(gè)SNP位點(diǎn)，經(jīng)過質(zhì)控過濾后，最終保留了5,685,673個(gè)高質(zhì)量SNP供后續(xù)分析。這些SNP廣泛分布于油茶基因組的各個(gè)染色體區(qū)域，具有較強(qiáng)的代表性和應(yīng)用潛力，為油茶的遺傳多樣性研究和分子育種提供了寶貴的遺傳資源。 Figure 2展示了SNP位點(diǎn)在各染色體上的分布密度、覆蓋深度和變異質(zhì)量參數(shù)，為后續(xù)的遺傳結(jié)構(gòu)分析及品種識(shí)別提供了堅(jiān)實(shí)的基礎(chǔ)數(shù)據(jù)支持。

3.各油茶品種的系統(tǒng)發(fā)育關(guān)系

為了系統(tǒng)闡明25個(gè)油茶品種之間的親緣關(guān)系與遺傳結(jié)構(gòu)，研究結(jié)合了鄰接法（NJ樹）、主成分分析（PCA）與群體結(jié)構(gòu)分析（ADMIXTURE）三種方法，進(jìn)行了系統(tǒng)發(fā)育與群體結(jié)構(gòu)分析（Figure 3）。NJ樹的結(jié)果顯示，大多數(shù)樣本能夠根據(jù)其品種進(jìn)行合理聚類，呈現(xiàn)出明顯的系統(tǒng)發(fā)育規(guī)律。然而，部分品種則呈現(xiàn)出混合的聚類模式，具體包括華碩HS與華鑫HX、湘林XL-3與湘林XL-210、以及湘林XL-1與長(zhǎng)林CL-4。而達(dá)林DL-22品種的個(gè)體聚類較為分散且無(wú)序，表現(xiàn)出較為復(fù)雜的遺傳結(jié)構(gòu)。在PCA分析中，前兩個(gè)主成分（PC1和PC2）解釋了大部分遺傳變異（PC1和PC2分別占11.20%和8.89%），并顯示出明顯的群體結(jié)構(gòu)分布，這與NJ樹的聚類結(jié)果高度一致。進(jìn)一步的群體結(jié)構(gòu)分析表明，當(dāng)K值為3時(shí)，分析結(jié)果最佳，將品種劃分為三大類群：第一類群包括達(dá)林DL-1、長(zhǎng)林CL-4和湘林XL-1；第二類群包括川榮CR-153、湘林XL-3和湘林XL-210；第三類群則包括其余所有品種。需要注意的是，DL-22品種的個(gè)體表現(xiàn)出較高的遺傳混雜性，這可能與其基因漸滲、引種或雜交歷史相關(guān)。這些分析結(jié)果為進(jìn)一步理解油茶品種的遺傳多樣性、進(jìn)化關(guān)系以及品種改良提供了重要的遺傳信息。

4.遺傳多樣性揭示品種差異

為了全面評(píng)估油茶種質(zhì)資源的遺傳多樣性，研究統(tǒng)計(jì)了包括雜合度（Heterozygosity）、核苷酸多樣性指數(shù)（π）、單態(tài)位點(diǎn)數(shù)量（Singleton）以及Watterson’s θ值等關(guān)鍵遺傳參數(shù)（Figure 4）。結(jié)果表明，川榮（CR）系列品種（如CR-153、CR-55）在各項(xiàng)遺傳變異指標(biāo)上均表現(xiàn)出較高水平，尤其是核苷酸多樣性指數(shù)（π）顯著高于其他品種。具體而言，singleton在長(zhǎng)林（CL）系列品種（如CL-3、CL-4、CL-40）、翠屏（CP）系列（如CP-16）和江安（JA）系列品種（如JA-54）中較少，而在川榮CR-153和CR-55中則為最多（Figure 4B）。核苷酸多樣性（π）值的差異也較為顯著，川榮品種CR-50、CR-156、CR-55和CR-153的π值明顯高于其他品種，表明川榮系列品種具有較高的遺傳多樣性，可能為油茶品種改良和分子育種提供重要的遺傳資源（Figure 4C）。類似地，Wattersonθ分析進(jìn)一步揭示了品種間分離位點(diǎn)的變化，川榮CR-153、湘林XL-3和CR-55表現(xiàn)出較高的θ值，與singleton的趨勢(shì)一致（Figure 4D）。近交系數(shù)（F）分析結(jié)果與遺傳多樣性模式一致，顯示川榮CR-156、CR-55和CR-447的F值最低，表明其遺傳多樣性較高，而長(zhǎng)林CL系列品種（如CL-3、CL-40、CL-4）的F值較高，反映其遺傳多樣性較低。總體而言，遺傳多樣性分析結(jié)果表明，長(zhǎng)林（CL）系列品種由于長(zhǎng)期的人工選育和馴化，遺傳背景較為單一，表現(xiàn)出較低的遺傳多樣性；而川榮系列品種則保持較高的遺傳多樣性。該研究結(jié)果為篩選核心種質(zhì)資源和制定育種策略提供了有力的遺傳依據(jù)。

5.遺傳分化分析明確親緣結(jié)構(gòu)

通過計(jì)算品種兩兩之間的Fst值與構(gòu)建遺傳距離矩陣，研究進(jìn)一步明確了油茶品種間的遺傳分化程度和親緣關(guān)系（Figure 5）。Fst值的分析結(jié)果顯示，長(zhǎng)林CL-4、達(dá)林DL-1與湘林XL-1三者之間的Fst值極低（0.0106–0.0110），這表明這些品種可能源自相同或相似的遺傳背景。類似地，川榮CR-153與湘林XL-3、湘林XL-210之間的Fst值也表現(xiàn)出較高的遺傳相似性（0.0036–0.0442），提示這些品種可能具有共同的祖源關(guān)系。相比之下，達(dá)林DL-22與川榮CR-50（0.0816）、川榮CR-55（0.0870）和達(dá)林DL-1（0.0759）之間的Fst值較高，顯示出較為顯著的遺傳分化，反映其與這些品種的親緣關(guān)系較為緊密，但仍保持一定的遺傳獨(dú)特性。其它品種的Fst值范圍從0.1013到0.2908，表明這些品種的遺傳起源更加多樣，可能涉及不同的地理區(qū)域或品種間的較大差異（Figure 5A）。研究人員進(jìn)一步使用GCTA軟件計(jì)算的G矩陣分析，結(jié)果顯示大多數(shù)油茶品種之間的親緣關(guān)系較為明顯，遺傳分化較大（Figure 5B）。然而，一些品種群體的親緣關(guān)系值相似，例如XL-1、CL-4與DL-1；CR-153、XL-3與XL-210；華碩HS與華鑫HX，表明這些品種可能源自相同或相近的祖先背景。與此不同，達(dá)林DL-22品種在其群體內(nèi)部未顯示出顯著的親緣關(guān)系，反映出其群體內(nèi)存在較高的遺傳分化。總體而言，這些分析為理解油茶品種的遺傳背景和親緣結(jié)構(gòu)提供了重要依據(jù)，揭示了不同品種之間的親緣關(guān)系及其遺傳多樣性。

6.構(gòu)建核心SNP指紋圖譜實(shí)現(xiàn)品種識(shí)別

為建立油茶品種的快速鑒定體系，研究人員進(jìn)一步構(gòu)建了25個(gè)主栽油茶品種的DNA指紋圖譜（Figure 6）。在初步獲得的234個(gè)候選SNP中，有196個(gè)位點(diǎn)被注釋為功能明確的編碼區(qū)域突變，富集于DNA連接、RNA加工、果糖代謝等重要通路（Figure 6A）。為提升區(qū)分效率，進(jìn)一步篩除掉存在多變異干擾或功能不明的位點(diǎn)，最終篩選得到15個(gè)核心SNP位點(diǎn)，可有效區(qū)分全部25個(gè)油茶品種，實(shí)現(xiàn)穩(wěn)定的遺傳身份標(biāo)識(shí)。為確保指紋圖譜的準(zhǔn)確性與應(yīng)用價(jià)值，研究通過Sanger測(cè)序?qū)?5個(gè)候選SNP位點(diǎn)進(jìn)行了引物驗(yàn)證，最終篩選出8個(gè)多態(tài)性強(qiáng)、峰型清晰的核心SNP位點(diǎn)用于構(gòu)建指紋圖譜（Figure 6B）。這些位點(diǎn)主要分布于功能基因區(qū)，涉及次生代謝、胚胎發(fā)育、自噬調(diào)控等關(guān)鍵生物過程，如編碼α-萜品醇合酶的LG12_G02320和控制胚胎發(fā)育與芽形成的FLA8蛋白基因LG03_G02445等，展現(xiàn)出良好的應(yīng)用潛力。為進(jìn)一步提升指紋信息的實(shí)用性，研究人員還開發(fā)了包含8個(gè)位點(diǎn)基因型信息的二維碼系統(tǒng)，其中涵蓋了油脂品質(zhì)、農(nóng)藝性狀及適宜種植區(qū)域等數(shù)據(jù)，為油茶品種的數(shù)字化管理與選育決策提供便捷、高效的技術(shù)支持（Figure 6C）。二維碼系統(tǒng)為油茶品種的數(shù)字化管理與選育決策提供了便捷、高效的技術(shù)支持。通過將分子標(biāo)記信息與油茶生產(chǎn)應(yīng)用相結(jié)合，研究成果不僅實(shí)現(xiàn)了分類識(shí)別功能，還為分子育種提供了重要參考價(jià)值。

總結(jié)

綜上，該DNA指紋圖譜的構(gòu)建為油茶種質(zhì)資源的快速鑒定、精準(zhǔn)管理及產(chǎn)業(yè)化育種提供了高效工具和理論依據(jù)。其開發(fā)的二維碼系統(tǒng)進(jìn)一步推動(dòng)了油茶品種管理的智能化和數(shù)據(jù)化，為我國(guó)油茶資源保護(hù)、優(yōu)質(zhì)品種選育及區(qū)域性產(chǎn)業(yè)布局提供了科學(xué)支持。

內(nèi)容來源于林木科學(xué)評(píng)論 侵刪

百邁客生物為該研究提供群體、轉(zhuǎn)錄組測(cè)序服務(wù)。

研究背景

高粱是世界第五大禾谷類作物，養(yǎng)活了世界干旱或半干旱地區(qū)5億多人口。高粱起源于非洲薩赫勒地區(qū)，對(duì)非生物脅迫，比如干旱、高溫、貧瘠、鹽堿等有良好的耐受性，是研究植物響應(yīng)逆境脅迫機(jī)制的理想模式物種。然而，有限的遺傳和基因組資源以及遺傳轉(zhuǎn)化技術(shù)的瓶頸限制了高粱功能基因的研究。

研究?jī)?nèi)容

研究團(tuán)隊(duì)利用ONT Ultra-long、PacBio HiFi和Hi-C技術(shù)，成功組裝了高粱E048 的T2T無(wú)缺口高質(zhì)量基因組圖譜。該基因組大小為729.5 Mb，預(yù)測(cè)蛋白編碼基因35,549個(gè)，并完整解析了全部端粒和著絲粒序列信息。通過與BTx623、Ji2055、HYZ等另外3個(gè)高粱T2T基因組的比較分析，研究團(tuán)隊(duì)鑒定出E048特異存在的2.9 Mb基因組區(qū)域，包含187個(gè)基因。這些基因主要參與信號(hào)轉(zhuǎn)導(dǎo)、免疫響應(yīng)和代謝調(diào)控等生物學(xué)過程，為解析E048高抗病、抗逆、抗倒伏和高含糖量等優(yōu)異性狀的分子機(jī)制提供了重要線索。

為了深入解析高粱E048的轉(zhuǎn)錄調(diào)控過程，研究團(tuán)隊(duì)測(cè)定了E048在13個(gè)不同發(fā)育階段多種組織（包括幼苗、葉片、根、莖、花序和種子等）的轉(zhuǎn)錄組數(shù)據(jù)，構(gòu)建了全面的基因表達(dá)參考圖譜。分析表明，E048基因組中93.39%（33,200個(gè)）的基因在至少一種組織中表達(dá)，20,475個(gè)基因（57.59%）在所有13種組織中均檢測(cè)到表達(dá)。其中在開花前的花序組織中特異表達(dá)的基因數(shù)量最多，GO富集分析發(fā)現(xiàn)些特異表達(dá)基因主要參與生殖相關(guān)的生物學(xué)過程，反映了高粱從營(yíng)養(yǎng)生長(zhǎng)階段向生殖生長(zhǎng)階段過渡時(shí)的關(guān)鍵轉(zhuǎn)錄特征。

研究團(tuán)隊(duì)進(jìn)一步開發(fā)了E048大規(guī)模飽和EMS突變體庫(kù)，在M1代收獲了13,226個(gè)單株并對(duì)179個(gè)M2植株進(jìn)行了全基因組重測(cè)序。對(duì)變異位點(diǎn)進(jìn)行分析發(fā)現(xiàn)，突變體包含了2,291,074個(gè)SNPs或InDels變異，突變位點(diǎn)幾乎均勻分布在每條染色體上，涵蓋了97.54%的基因，這些突變體為高粱功能基因解析與遺傳改良提供了重要種質(zhì)資源。研究團(tuán)隊(duì)利用E048基因組，通過MutMap和MapMap+的方法，快速定位了黃化突變體和多葉突變體的目的基因，驗(yàn)證了利用E048基因組和突變體資源進(jìn)行高粱功能基因研究的高效性。

為促進(jìn)高粱功能基因組研究以及突變體庫(kù)的廣泛應(yīng)用，研究團(tuán)隊(duì)還建立了高粱基因組和突變體庫(kù)數(shù)據(jù)庫(kù)SGMD（https://sorghum.genetics.ac.cn/SGMD）。該數(shù)據(jù)庫(kù)整合了高粱E048基因組信息、基因表達(dá)圖譜、突變體庫(kù)部分表型變異和基因變異信息等，是高粱功能基因研究的“超級(jí)工具包”。

內(nèi)容來源于MPlant植物科學(xué)，侵刪

以下小編列舉三篇成功案例，解析群體在QTL定位中的應(yīng)用。

當(dāng)我們針對(duì)一群體，關(guān)注性狀較多時(shí)，可參照案例一，整篇文章通過群體圖譜構(gòu)建，QTL定位后幾乎沒有驗(yàn)證工作，可以幫我們拿到一個(gè)群體QTL的基礎(chǔ)數(shù)據(jù)。

案例一

發(fā)表期刊：Plant Physiology and Biochemistry

影響因子：6.1

合作單位：江蘇省徐淮區(qū)徐州農(nóng)業(yè)科學(xué)研究所農(nóng)業(yè)農(nóng)村部甘薯生物學(xué)與遺傳改良重點(diǎn)實(shí)驗(yàn)室

實(shí)驗(yàn)方法：Xin24×Yushu10雜交中選擇212個(gè)F1材料: 特異性位點(diǎn)擴(kuò)增片段測(cè)序（SLAF-seq）；遺傳圖譜構(gòu)建；數(shù)量性狀位點(diǎn)（QTL）定位；GWAS關(guān)聯(lián)（F1群體）

表型檢測(cè)：多年多表型收集

百邁客生物為該研究提供了群體測(cè)序及部分?jǐn)?shù)據(jù)分析服務(wù)。

甘薯作為一個(gè)全球重要的作物，其含有豐富的營(yíng)養(yǎng)成分以及對(duì)不同環(huán)境的適應(yīng)能力。然而，由于其自交不親和，高雜合度等特性，使其遺傳特征的研究相對(duì)較少。作者從Xin24×Yushu10雜交中選擇212個(gè)F1材料，SLAF-seq測(cè)序，獲得親本26.73×，子代52.25×的測(cè)序數(shù)據(jù)，依據(jù)SNP及百邁客生物自主研發(fā)的HighMap構(gòu)圖軟件，生成一個(gè)長(zhǎng)度為2441.56 cM、平均圖距為0.51 cM的遺傳圖譜。

基于連鎖圖譜，鑒定出26個(gè)QTL，解釋了6.3-10%的表型變異，包括6個(gè)最長(zhǎng)藤蔓長(zhǎng)度 QTL、6個(gè)單株產(chǎn)量 QTL、10個(gè)干物質(zhì)含量QTL、1個(gè)淀粉含量 QTL、一個(gè)可溶性糖含量QTL和2個(gè)類胡蘿卜素含量QTL。該研究結(jié)果對(duì)甘薯的標(biāo)記輔助育種和基因克隆具有重要意義。

圖1-表型檢測(cè)

圖2-遺傳圖譜構(gòu)建

圖3-QTL定位

前文是對(duì)多個(gè)性狀的連鎖分析，當(dāng)我們關(guān)注單個(gè)性狀時(shí)，BSA無(wú)疑是高性價(jià)比的初定位選擇，當(dāng)然，這就意味著我們得做到基因的精細(xì)定位與克隆，除去傳統(tǒng)圖位克隆的方式，轉(zhuǎn)錄組，蛋白組，自然群體GWAS，基因組都可助力基因克隆，如果想發(fā)高水平的文章，基因的功能探索也是必不可少的。

成功案例二

發(fā)表期刊：Plant Biotechnology Journal

影響因子：10.1

發(fā)表單位：河南農(nóng)業(yè)大學(xué)

實(shí)驗(yàn)方法：莢果大小/重量差異顯著2份Virginia-type花生材料( ND _ L和ND _ S)雜交，產(chǎn)生遺傳群體；F2:3群體BSA-seq（20+20混池）；F6:7 和F6:8群體精細(xì)定位；基因克??；系統(tǒng)進(jìn)化分析；亞細(xì)胞定位；免疫熒光；酵母雙雜交；pull-down；CO-IP；番茄擬南芥轉(zhuǎn)化實(shí)驗(yàn)等。

百邁客生物為該研究提供了群體測(cè)序及部分?jǐn)?shù)據(jù)分析服務(wù)。

花生莢果大小是決定花生產(chǎn)量的關(guān)鍵農(nóng)藝性狀，為了鑒定控制花生莢果大小的基因，該研究對(duì)188份核心種質(zhì)進(jìn)行鑒定，并選取莢果大小/重量差異顯著的2份花生材料構(gòu)建F2群體，BSA-Seq獲得了288.58 Gb的原始數(shù)據(jù)。利用285914個(gè)高質(zhì)量SNPs 和70 759個(gè) InDel，將控制莢果大小的基因定位在07染色體1.17 Mb的區(qū)間內(nèi)。作者從F6：7和F6：8群體中開發(fā)了15個(gè)多態(tài)性標(biāo)記并對(duì)個(gè)體進(jìn)行基因分型，精細(xì)定位QTL到KASP 9 和In Del15之間20.4 kb的區(qū)間內(nèi)。該區(qū)間僅包含1個(gè)預(yù)測(cè)的非同義突變基因和InDels，將該基因命名為PSW1。

PSW1編碼一個(gè)LRR – RLK蛋白激酶，等位基因PSW1HapII賦予了PSW1更高的表達(dá)水平和對(duì)其輔助受體AhBAK1更強(qiáng)的親和力，以上調(diào)PSW1 – based途徑，調(diào)節(jié)花生莢果大小。此外，PSW1HapII的過表達(dá)增加了多種植物的種子/果實(shí)大小。

圖4-表型檢測(cè)

圖5-BSA分析

當(dāng)然，如果想讓我們的文章影響因子再上一臺(tái)階，兼顧群體的“廣度”和“深度”，是更好的選擇。

成功案例三

發(fā)表期刊：Nature?Genetics

影響因子：30.8

發(fā)表單位：浙江省農(nóng)業(yè)科學(xué)院等

實(shí)驗(yàn)方法：菜用豇豆G98，糧用豇豆G323 基因組Denovo；重測(cè)序GWAS：344份全世界收集的豇豆核心種質(zhì)，其中包括342份栽培豇豆（87份糧用豇豆、244份菜用豇豆和11份未知用途豇豆）和2份野生豇豆；Illumina測(cè)序，10x深度；基因單倍型驗(yàn)證：菜用豇豆地方品種 ‘ZN016’ 和菜用豇豆育成品種‘Zhijiang282構(gòu)建的RIL群體（183 lines）G98和G323構(gòu)建的F2群體（165 individuals）

百邁客生物為該研究提供了群體測(cè)序、基因組測(cè)序及部分?jǐn)?shù)據(jù)分析服務(wù)。

豇豆起源于非洲，在世界范圍內(nèi)作為糧食、蔬菜或牲畜飼料種植。該研究結(jié)合PacBio、Hi-C和二代測(cè)序，組裝了糧用豇豆和菜用豇豆的染色體水平基因組。對(duì)包括地方品種、野生品種和育成品種的344個(gè)材料進(jìn)行二代測(cè)序，以闡明豇豆基因組的系統(tǒng)進(jìn)化。

為了研究自然或人工選擇對(duì)豇豆分化的影響，作者通過選擇清除分析比較了三個(gè)豇豆亞群之間的基因組選擇特征，鑒定出239個(gè)與豇豆馴化和改良相關(guān)的基因。此外，通過GWAS，挖掘到裂莢性、莢長(zhǎng)、單莢粒數(shù)、千粒重、可溶性糖、總淀粉和粗蛋白質(zhì)含量相關(guān)基因，并在遺傳群體中驗(yàn)證。同時(shí)揭示了兩個(gè)亞種之間基因組結(jié)構(gòu)變異（SVs）的全圖譜，為豇豆在全基因組選擇下的馴化與改良提供了見解。產(chǎn)量性狀和品質(zhì)性狀的差異基因組選擇將有助于建立糧用豇豆和菜用豇豆雙向改良的遺傳資源。

圖6-群體選擇與GWAS分析

今日分享結(jié)束，期待下期精彩內(nèi)容~~~

發(fā)表期刊：Poultry Science

影響因子：3.8

發(fā)表單位：南京農(nóng)業(yè)大學(xué)，西藏農(nóng)牧學(xué)院等

研究對(duì)象：金陵白鴨

研究方法：全基因組關(guān)聯(lián)分析（GWAS）

百邁客生物為該研究提供了全基因組關(guān)聯(lián)分析（GWAS）測(cè)序及部分?jǐn)?shù)據(jù)分析服務(wù)。

研究背景

鴨子為人類消費(fèi)提供肉、蛋、羽毛等產(chǎn)品，中國(guó)是世界上最大的鴨肉消費(fèi)國(guó)。金陵白鴨做為我國(guó)優(yōu)質(zhì)的選育品種，具有優(yōu)良的生長(zhǎng)速度和肉質(zhì)品質(zhì)。從孵化到上市的整個(gè)時(shí)期，鴨子的體重逐漸增加，這一過程由復(fù)雜的生理和生化機(jī)制調(diào)控，器官和骨骼的發(fā)育是體重增加的關(guān)鍵因素，全基因組關(guān)聯(lián)分析（GWAS）是許多牲畜和家禽物種研究生長(zhǎng)和發(fā)育的遺傳機(jī)制的常用方法，然而，現(xiàn)有的一些關(guān)于鴨生長(zhǎng)發(fā)育性狀的GWAS研究測(cè)序深度較低，發(fā)現(xiàn)的候選基因很少。

材料方法

201只雄鴨全基因組測(cè)序深度為10×；

GWAS分析：表型：出生體重（BWB）、1周體重（BW1）、3周體重（BW3）、5周體重（BW5）和7周體重（BW7）分析模型：FaST-LMM；EMMAX；LMM；LM 閾值線：log10(p)=5如果一個(gè)SNP在2個(gè)或更多的模型中被關(guān)聯(lián)，作者認(rèn)為它是一個(gè)與特定性狀相關(guān)的高可信度SNP。

進(jìn)化分析：進(jìn)化樹構(gòu)建；群體結(jié)構(gòu)分析；PCA分析；LD衰減分析等

研究結(jié)果

1.金陵白鴨群體體重表型統(tǒng)計(jì)

該研究對(duì)金陵白鴨（n = 201）在不同時(shí)期的體重進(jìn)行統(tǒng)計(jì)，BWB、BW1、BW3、BW5和BW7的體重均值分別為46.93g、127.19g、697.17g、1182.43g和1,888.52g。生長(zhǎng)曲線結(jié)果顯示：金陵白鴨第3周-第5周生長(zhǎng)最快，第5周后生長(zhǎng)放緩。體重的分散度隨著年齡和體重的增加而增加。相鄰2周的體重之間存在較強(qiáng)的相關(guān)性，BW3與BW1（r = 0.6）、BW5（r = 0.68）和BW7（r = 0.6）顯著正相關(guān)，BW5與BW7顯著正相關(guān)（r = 0.82）。但隨著時(shí)間間隔的延長(zhǎng)，性狀間的相關(guān)性降低，BW1與BW7的相關(guān)性不顯著（r = 0.21），BWB則與BW3 (r = 0.18)， BW5 (r = 0.058)，BW7 (r = 0.45)均不顯著相關(guān)。

圖1-表型分析

2.金陵白鴨系統(tǒng)發(fā)育、群體遺傳結(jié)構(gòu)解析

雖然該研究只有金陵白鴨一個(gè)種群，但考慮到繁殖過程中種群分層的潛力，作者進(jìn)行了系統(tǒng)發(fā)育和群體結(jié)構(gòu)分析。系統(tǒng)發(fā)育樹和PCA結(jié)果顯示，金陵白鴨可分為5個(gè)種群，群體結(jié)構(gòu)表明，當(dāng)K=6時(shí)分群結(jié)果最佳，金陵白鴨種群具有豐富的遺傳多樣性，總SNPs的平均R2為0.24，當(dāng)R2=0.2，LD的衰減距離約為30kbp。

圖2-群體進(jìn)化分析

3.GWAS分析

該研究對(duì)基于測(cè)序產(chǎn)生的2,610.50 Gbp的clean data進(jìn)行分析，Q30為95.55%，樣本與參考基因組平均比對(duì)率為99.59%，平均覆蓋深度為10×，基因組覆蓋度為96.51%。作者利用，797,309,337個(gè)SNPs，4種關(guān)聯(lián)模型進(jìn)行后續(xù)GWAS分析，其中95個(gè)SNP與BWB性狀顯著相關(guān)，這些SNP主要分布在1號(hào)染色體和4號(hào)染色體上。通過篩選和注釋，共檢測(cè)到5個(gè)相關(guān)的候選基因：PUS7、FBXO11、FOXN2、MSH6、SLC4A4。針對(duì)BW1性狀，作者發(fā)現(xiàn)了101個(gè)與BW1性狀顯著相關(guān)的SNPs，這些SNP主要分布在7號(hào)染色體上，共檢測(cè)到2個(gè)與BW1性狀相關(guān)的候選基因：RAG2和TMEFF2。針對(duì)BW3性狀，作者發(fā)現(xiàn)了112個(gè)顯著SNPs。這些SNP主要分布在1號(hào)染色體和11號(hào)染色體上。通過篩選和注釋，共檢測(cè)到4個(gè)與BW3性狀相關(guān)的候選基因：?STARD13、Klotho、ZAR1L和TLE3。同時(shí)確定了92個(gè)與BW5性狀相關(guān)SNPs，這些SNP主要分布在1號(hào)染色體和2號(hào)染色體上，通過注釋STARD13、Klotho和ZAR1L與BW5性狀相關(guān)。此外，作者還鑒定了新的候選基因：KAT2B、KCNH8和SATB1。

BW7性狀與33個(gè)SNP關(guān)聯(lián)，這些SNP主要分布在1號(hào)染色體和2號(hào)染色體上。通過篩選和注釋，共檢測(cè)到6個(gè)與BW7性狀相關(guān)的候選基因：PLXNC1、ATP1A1、CD58、FRYL、OCIAD1和OCIAD2。在分析的四個(gè)模型的結(jié)果，LM識(shí)別出了更多的顯著位點(diǎn)，曼哈頓圖顯示出更高的峰值。然而，QQ-plot顯示，LM模型中的大部分站點(diǎn)都位于對(duì)角線上方，可能具有較高的假陽(yáng)性。其余3個(gè)模型的結(jié)果均表現(xiàn)出一致性，而QQ-plot的結(jié)果則優(yōu)于LM模型。

圖3-GWAS分析

研究總結(jié)

金陵白鴨是一種新開發(fā)的品種，因其生長(zhǎng)速度快，肉質(zhì)優(yōu)良的特點(diǎn)，使其具有重要的經(jīng)濟(jì)價(jià)值和研究潛力；然而，人們對(duì)其體重性狀的遺傳基礎(chǔ)尚不太清楚。該研究對(duì)201只金陵白公鴨進(jìn)行了全基因組重測(cè)序，并進(jìn)行了群體基因組分析，表明金陵白鴨種群具有豐富的遺傳多樣性。

作者對(duì)出生體重（BWB）、1周體重（BW1）、3周體重（BW3）、5周體重（BW5）和7周體重（BW7）進(jìn)行了全基因組關(guān)聯(lián)分析，4種統(tǒng)計(jì)模型比較研究表明，F(xiàn)aST-LMM表現(xiàn)出最優(yōu)的效率，產(chǎn)生更多的結(jié)果和最小的假陽(yáng)性。

作者發(fā)現(xiàn)，PUS7、FBXO11、FOXN2、MSH6和SLC4A4均與BWB相關(guān)。RAG2和TMEFF2是BW1的候選基因，STARD13、Klotho、ZAR1L可能是BW3和BW5的候選基因。PLXNC1、ATP1A1、CD58、FRYL、OCIAD1和OCIAD2與BW7相關(guān)。這些研究結(jié)果為金陵白鴨的選擇和育種提供了遺傳參考，同時(shí)也加深研究者們對(duì)白鴨生長(zhǎng)發(fā)育表型的認(rèn)識(shí)。

沈陽(yáng)醫(yī)學(xué)院高兵教授團(tuán)隊(duì)在?Redox Biology?雜志發(fā)表題為“ Meta-data analysis of kidney stone disease highlights?ATP1A1?involvement in renal crystal formation ”?的研究論文，結(jié)合全外顯子測(cè)序與轉(zhuǎn)錄組等組學(xué)數(shù)據(jù)分析，解釋ATP1A1在腎晶體形成中起到的重要作用，表明ATP1A1可能是治療鈣結(jié)石的潛在治療靶點(diǎn)。ong>

英文標(biāo)題：Meta-data analysis of kidney stone disease highlights?ATP1A1?involvement in renal crystal formation

中文標(biāo)題：腎結(jié)石病的Meta數(shù)據(jù)分析揭示ATP1A1參與腎晶體形成

發(fā)表期刊：Redox Biology

合作單位：沈陽(yáng)醫(yī)學(xué)院

影響因子：10.7

研究對(duì)象：CaOx腎結(jié)石患者及對(duì)照組樣本

研究方法：全外顯子測(cè)序（WES）+表達(dá)譜數(shù)據(jù)分析

百邁客生物為該研究提供了外顯子測(cè)序服務(wù)。

研究背景

腎結(jié)石在全球范圍內(nèi)都威脅著人類的生命健康，近年來發(fā)病率和復(fù)發(fā)率都不斷上升，大多數(shù)腎結(jié)石由草酸鈣（Calcium ox）結(jié)石組成，晶體-細(xì)胞粘附被認(rèn)為是晶體保留和結(jié)石形成中的關(guān)鍵步驟，有研究報(bào)道，受損的腎小管上皮細(xì)胞對(duì)晶體附著的親和力增加，而晶體沉積通過產(chǎn)生過多的活性氧（ROS）反過來又誘導(dǎo)細(xì)胞損傷和炎癥反應(yīng)。Na/K-ATPase（NKA）是一種在腎臟中高度表達(dá)的跨膜離子泵，同時(shí)其由ATP1A1編碼的α1亞基也起到信號(hào)轉(zhuǎn)導(dǎo)器的作用，當(dāng)?ATP1A1?構(gòu)象改變或被敲低時(shí)，Src 就會(huì)被磷酸化并激活下游信號(hào)級(jí)聯(lián)，進(jìn)而導(dǎo)致 ROS 的產(chǎn)生，ROS可以激活MAPK級(jí)聯(lián)反應(yīng)和NF-κB炎癥反應(yīng)，并作為ATP1A1的配體啟動(dòng)ATP1A1/Src信號(hào)通路，從而形成NKA/ROS擴(kuò)增環(huán)，誘導(dǎo)氧化應(yīng)激，然而，ATP1A1在腎結(jié)石形成中的作用尚不清楚。本實(shí)驗(yàn)結(jié)合多組學(xué)分析方法，以期提供一種有用的方法，研究腎結(jié)石形成病理。

材料方法

實(shí)驗(yàn)材料：28例中國(guó)漢族CaOx腎結(jié)石患者及對(duì)照組樣本

組學(xué)方法：全外顯子測(cè)序（WES）+表達(dá)譜數(shù)據(jù)分析

研究結(jié)果

首先該研究從GEO數(shù)據(jù)中獲取了腎臟鈣結(jié)石患者的表達(dá)譜數(shù)據(jù)化，其中患者的RP組織（Randall斑塊）命名為P組，患者的正常組織命名為N組，對(duì)其表達(dá)譜篩選基因進(jìn)行WGCNA分析，選擇與樣本相關(guān)性較高的淺橙色模塊（434個(gè)基因）進(jìn)行后續(xù)研究，GO富集結(jié)果表明這些基因在質(zhì)膜和跨膜轉(zhuǎn)運(yùn)中富集。KEGG富集分析顯示，這些基因參與“醛固酮調(diào)節(jié)的鈉重吸收”、“礦物質(zhì)吸收”、“內(nèi)分泌等因子調(diào)節(jié)的鈣重吸收”和“近端小管碳酸氫鹽回收”等通路，進(jìn)一步查看發(fā)現(xiàn)，FXYD2、FXYD4、ATP1A4、ATP1A1和ATP1B1等5個(gè)編碼NKA的基因參與了跨膜轉(zhuǎn)運(yùn)、吸收和重吸收的過程，表明NKA可能與腎結(jié)石密切相關(guān)。

為進(jìn)一步探究腎結(jié)石相關(guān)關(guān)鍵基因，該研究對(duì)28例中國(guó)漢族CaOx結(jié)石患者進(jìn)行了WES分析，共鑒定得到76個(gè)SNPs，其中錯(cuò)義SNPs31個(gè)，編碼同義SNPs30個(gè)，3′UTR7個(gè)SNPs，5′UTR區(qū)8個(gè)SNPs，作為候選SNPs，這些突變位點(diǎn)定位于67個(gè)基因，將這67個(gè)基因與WGCNA分析中的434個(gè)基因做韋恩圖取交集，發(fā)現(xiàn)ATP1A1是兩個(gè)數(shù)據(jù)集交叉點(diǎn)中的唯一基因，一定程度上說明了ATP1A1在鈣結(jié)石形成中的潛在作用。

接下來，該研究從WES數(shù)據(jù)中提取了ATP1A1的SNP，rs11540947 T等位基因在患者組中的突變頻率頻率為16.1%，而對(duì)照組僅為2.4%。為進(jìn)一步證實(shí)rs11540947與鈣結(jié)石形成的相關(guān)性，收集214例鈣結(jié)石患者和232例匹配對(duì)照，通過HRM分析檢測(cè)rs11540947的基因型，結(jié)果顯示，CT+TT攜帶者在患者中比對(duì)照組更常見，且在男性中與鈣結(jié)石有關(guān)，但在女性中無(wú)顯著相關(guān)性，SNPrs11540947 （NM_000701:c.-78C > T）位于?ATP1A1?的 5′UTR 中；在ATP1A1的5′UTR中發(fā)現(xiàn)了Sp1結(jié)合位點(diǎn)和潛在的TATA盒，表明5′UTR具有啟動(dòng)子活性，雙熒光素酶報(bào)告基因結(jié)果也證明了這一點(diǎn)，rs11540947的T等位基因介導(dǎo)了ATP1A1啟動(dòng)子活性的顯著降低，這可能會(huì)降低ATP1A1在轉(zhuǎn)錄水平的表達(dá)。對(duì)草酸鈣一水合物（COM）暴露處理的HK2細(xì)胞（人腎近曲小管上皮細(xì)胞）進(jìn)行指標(biāo)測(cè)定，結(jié)果顯示處理3h后，NKA酶活性和ATP1A1 mRNA水平顯著下降，6h后ATP1A1的蛋白水平升高，隨后降低，同時(shí)COM 暴露后細(xì)胞內(nèi) ROS 水平顯著增加，這可能與ATP1A1/Src信號(hào)通路的激活有關(guān)，而ATP1A1的過表達(dá)抑制了COM誘導(dǎo)的Src、p38、JNK、p65和p50激活；WB結(jié)果也顯示COM激活的半胱天冬酶3被ATP1A1的過表達(dá)以及ATP1A1/Src信號(hào)復(fù)合物的特異性拮抗劑pNaKtide顯著抑制，從而阻斷了細(xì)胞凋亡過程。

晶體與細(xì)胞的粘附是結(jié)石形成的關(guān)鍵過程，晶體沉積會(huì)損傷腎細(xì)胞，因此，作者研究了ATP1A1在晶體-細(xì)胞粘附中的作用。該研究觀察到用 Ad-hATP1A1 感染和用 pNaKtide 處理后，Ca²⁺濃度與對(duì)照組相比顯著降低，表明ATP1A1過表達(dá)和pNaKtide可以阻止晶體與細(xì)胞的粘附，從而保護(hù)細(xì)胞免受CaOx晶體誘導(dǎo)的損傷并減少腎結(jié)石的形成。DNA啟動(dòng)子中CpG島的甲基化是基因沉默的常見機(jī)制，COM暴露后，DNA甲基轉(zhuǎn)移酶（DNMTs）的mRNA和蛋白質(zhì)水平均顯著升高，而用DNA甲基化抑制劑5Aza-2dc預(yù)處理HK2細(xì)胞再進(jìn)行COM暴露處理后，ATP1A1表達(dá)量下調(diào)；ATP1A1的蛋白水平也被逆轉(zhuǎn)。這些結(jié)果表明，改變的DNA甲基化參與了COM誘導(dǎo)的ATP1A1降低。

為了證實(shí)ATP1A1對(duì)體內(nèi)晶體形成的影響，該研究構(gòu)建了CaOx腎病大鼠模型，用 pNaKtide 治療后尿草酸鹽排泄量顯著減少，與體外研究結(jié)果類似，HYP+CaCl2處理降低了 ATP1A1 mRNA和蛋白質(zhì)表達(dá)水平，激活了 Src、p38、JNK、p65 和 p50，降低 Nrf2 表達(dá)，而這些均被pNaKtide逆轉(zhuǎn)，這些結(jié)果強(qiáng)調(diào)了ATP1A1和ATP1A1/Src信號(hào)通路在腎結(jié)石形成中的重要性。

研究總結(jié)

綜上所述，該研究從遺傳水平和環(huán)境因子影響的轉(zhuǎn)錄水平探究腎結(jié)石相關(guān)關(guān)鍵易感基因，將鈣結(jié)石形成者Randall斑塊（RP）組織的基因表達(dá)譜數(shù)據(jù)與CaOx患者的WES數(shù)據(jù)相結(jié)合，鑒定出ATP1A1這一關(guān)鍵基因，還研究了ATP1A1遺傳變異與鈣結(jié)石風(fēng)險(xiǎn)的關(guān)聯(lián)，并進(jìn)一步探討了ATP1A1/Src/ROS信號(hào)在體外和體內(nèi)腎結(jié)石形成中的作用，這項(xiàng)研究的發(fā)現(xiàn)揭示了ATP1A1基因變異及其表達(dá)降低參與腎晶體形成的機(jī)制，可能為CaOx結(jié)石提供潛在的治療靶點(diǎn)。

參考文獻(xiàn)：

Li Y, Lu X, Yu Z, Wang H, Gao B. Meta-data analysis of kidney stone disease highlights?ATP1A1?involvement in renal crystal formation.?Redox Biol. 2023 May;61:102648. doi: 10.1016/j.redox.2023.102648. Epub 2023 Feb 27. PMID:36871182.

該研究以瓢雞和仙居雞為親本建立回交家系，以探索瓢雞無(wú)尾特征的遺傳機(jī)制和分子基礎(chǔ)。通過全基因組關(guān)聯(lián)和連鎖分析，研究人員將無(wú)尾性狀的候選區(qū)域精確定位到798.5 kb的區(qū)間內(nèi)（chr2:86.9-87.7 Mb）。

通過對(duì)包含家系內(nèi)特殊基因型個(gè)體的突變位點(diǎn)進(jìn)行綜合篩選和分析后，一個(gè)4.2 kb的缺失被確定與瓢雞的無(wú)尾表型完全相關(guān)。在對(duì)致因區(qū)間內(nèi)的基因表達(dá)結(jié)果進(jìn)行深入探究后，一個(gè)全新的基因Rum（長(zhǎng)度大于22 kb，無(wú)內(nèi)含子），其表達(dá)呈現(xiàn)出無(wú)尾表達(dá)缺失，雜合子是顯性純合子表達(dá)量的一半的現(xiàn)象。深入研究發(fā)現(xiàn)Rum的表達(dá)具有胚胎特異性，研究人員認(rèn)為Rum的表達(dá)缺失以及雜合子中的單倍劑量不足，使其無(wú)法正常調(diào)控MSGN1基因長(zhǎng)轉(zhuǎn)錄本以及TBX6等基因的正常表達(dá)，進(jìn)而使得尾部骨骼發(fā)育出現(xiàn)異常。

瓢雞無(wú)尾表型的定位區(qū)間，與之前報(bào)道的美國(guó)阿勞卡納雞品種無(wú)尾表型區(qū)間非常接近但精細(xì)定位區(qū)間沒有重疊。對(duì)瓢雞、阿勞卡納以及其他具有正常尾巴的多個(gè)種群的群體基因組分析顯示，雖然這兩個(gè)無(wú)尾品種在2號(hào)染色體的相同區(qū)域受到選擇，但選擇具有品種差異性。在阿勞卡納雞2號(hào)染色體上識(shí)別的兩個(gè)候選SNP在無(wú)尾瓢雞中并不存在，并且在瓢雞中證實(shí)的致因缺失突變?cè)诎诳{雞也不存在。即瓢雞和阿勞卡納雞無(wú)尾特征可能源自不同的致病突變。

值得注意的是，阿勞卡納和瓢雞是遺傳上獨(dú)特的品種，兩者沒有共同祖先。瓢雞源自中國(guó)云南，阿勞卡納起源于智利后引入美國(guó)。兩種無(wú)尾雞在外觀表型上也有所不同，包括體型、耳簇和蛋殼顏色等。在當(dāng)前的研究中，確認(rèn)了這兩個(gè)品種特有的致因突變存在。因此，如果這些不同的遺傳背景導(dǎo)致了相同的表型，這意味著自然已經(jīng)進(jìn)化出了替代的遺傳解決方案來實(shí)現(xiàn)所需的表型特征。這種遺傳路徑的冗余匯聚形成了共同的表型，是進(jìn)化多樣性的顯著例證。

研究發(fā)現(xiàn)了瓢雞無(wú)尾的致因突變和關(guān)鍵基因，揭示了關(guān)鍵基因的作用方式。這項(xiàng)研究不僅增進(jìn)了我們對(duì)于雞骨骼發(fā)育調(diào)控的理解。還強(qiáng)調(diào)了遺傳演化和適應(yīng)性的動(dòng)態(tài)本質(zhì)，是鳥類遺傳和進(jìn)化研究領(lǐng)域的又一突破性的發(fā)現(xiàn)。

 

期刊：Nature Genetics

IF：30.8

時(shí)間：2023年 7月

全文：點(diǎn)擊下載

作者通過秈稻品種特青和粳稻品種02428構(gòu)建高級(jí)遺傳群體克隆了產(chǎn)量相關(guān)基因GY3，其啟動(dòng)子區(qū)域的反轉(zhuǎn)座子插入增強(qiáng)了啟動(dòng)子區(qū)域的表觀修飾，降低GY3表達(dá)量，從而減少細(xì)胞分裂素合成前體底物的無(wú)效消耗，提高體內(nèi)活性的細(xì)胞分裂素含量，增加每穗粒數(shù)和谷物產(chǎn)量。GY3可作為后期秈稻高產(chǎn)育種的重要目標(biāo)基因，推動(dòng)產(chǎn)量提升。

野生八倍體草莓單倍型基因組解析

期刊：Nature Plants

IF：18.0

時(shí)間：2023年7月

全文：點(diǎn)擊下載

該研究報(bào)道了兩個(gè)八倍體野生草莓智利草莓 (F. chiloensis)和弗州草莓 (F. virginiana)染色體水平的高質(zhì)量分型基因組組裝。探討了八倍體草莓在馴化過程中的同源偏向性表達(dá)分化，并鑒定到了一些重要的轉(zhuǎn)錄因子在馴化過程中發(fā)生了表達(dá)轉(zhuǎn)變。此研究深度解析了草莓的起源、基因組進(jìn)化和馴化。（百邁客參與測(cè)序工作·）

橡膠樹基因組解析助力馴化研究

期刊：Nature Communications

IF：16.6

時(shí)間：2023年7月

全文：點(diǎn)擊下載

作者利用PacBio CLR＋BioNano＋Hi-C的方法構(gòu)建了橡膠樹參考基因組，結(jié)合127個(gè)栽培和208個(gè)野生品種重測(cè)序數(shù)據(jù)對(duì)其馴化的分子機(jī)制解析，同時(shí)通過208份橡膠樹種質(zhì)材料進(jìn)行GWAS分析，鑒定到一個(gè)調(diào)控乳管數(shù)量的馴化基因。此研究將為后期橡膠樹的高產(chǎn)育種提供了分子基礎(chǔ)。

被子植物6大器官模式圖

1、在與研究目的不沖突的前提下，盡量選取新鮮、幼嫩、生長(zhǎng)狀態(tài)良好的組織部位。植物組織越幼嫩新鮮所含的次生代謝產(chǎn)物就越少，隨著植物組織的逐漸成熟，次生代謝產(chǎn)物的含量會(huì)逐漸增多，而次生代謝產(chǎn)物的存在會(huì)影響核酸的提取效果，次生代謝產(chǎn)物越多，核酸提取越困難，得率也越低。

2、如果植物組織表面污漬較多，在采集之前應(yīng)迅速用預(yù)冷的純凈水或75%乙醇擦拭或沖洗干凈，再用吸水紙吸干組織表面殘留液體。

3、將處理好的組織樣本混合均勻后保存于2 mL或更大體積的旋蓋凍存管中。

4、建議10min中內(nèi)（越快越好）置于液氮中冷凍1h以上，然后轉(zhuǎn)移至-80℃長(zhǎng)期保存，干冰運(yùn)輸。

1）提前準(zhǔn)備好冷凍組織的足量液氮，預(yù)裝樣品的凍存管，并在做好標(biāo)記（盡量不使用漢字命名，命名字符控制在5個(gè)以內(nèi)）

2）活體取下新鮮組織，立即剔除結(jié)締組織等非研究所需的組織類型。對(duì)腫瘤組織的取材，應(yīng)盡可能準(zhǔn)確地判定腫瘤和正常組織，腫瘤組織應(yīng)將周圍的正常組織切除干凈（正常組織也應(yīng)周圍的腫瘤組織切除干凈），腸道組織一定要把內(nèi)容物清洗干凈；

3）迅速用預(yù)冷的PBS溶液（RNase free）或0.9%生理鹽水將組織表面的殘留血液沖洗干凈

4）如果組織體積較大，將組織切成長(zhǎng)寬高均≤0.5 mm 的小塊（即黃豆大?。?/p>

5）將凍處理好的組織樣本混合均勻后保存于2 mL 或更大體積的螺口凍存管（RNase free）中，標(biāo)注編號(hào)。

6）立即（20s內(nèi)）置于液氮中冷凍 3~4 h，然后轉(zhuǎn)移至-80 °C長(zhǎng)期保存

7）干冰運(yùn)輸。

2）震蕩混勻，常溫裂解 5 min 后，使用低溫離心機(jī)12000g離心10min，將上清液轉(zhuǎn)至新的2.0ml離心管中（拍質(zhì)控照片，方便核查），轉(zhuǎn)移至-80℃低溫保存，運(yùn)輸時(shí)選擇干冰寄送。

2）如果組織帶血液或其他體液，需要過夜后更換一次RNAlater??Solution；不要將剛浸入RNAlater Solution中的樣本立即冷凍，需將樣本置于4 °C 保存過夜（使Solution充分浸潤(rùn)組織樣本），然后轉(zhuǎn)移至-20 °C 或-80 °C 長(zhǎng)期保存；

3）RNAlater??Solution 不會(huì)影響組織結(jié)構(gòu)，可以把已保存的組織從 RNAlater??Solution中取出，切下實(shí)驗(yàn)所需的用量，把剩余的組織再放入到原來的保存液中繼續(xù)保存；

4）保存于RNAlater??Solution中樣本，-20 °C 存放時(shí)，樣本不會(huì)結(jié)凍，但可能會(huì)有晶體析出，這并不影響后續(xù)的 RNA 提取工作；-80 °C 存放時(shí)，樣本會(huì)結(jié)凍， 在進(jìn)行 RNA 提取前，需置于冰上融化再進(jìn)行后續(xù)操作，解凍后的樣本可再次放入-80 °C 保存；

5）一般來講，生物樣本保存于RNAlater??Solution中，37 °C可存放1天，25 °C（室溫）可存放1周，4 °C 可存放1個(gè)月，-20 °C 或-80 °C 可長(zhǎng)期保存。但鑒于生物樣本的特殊性及實(shí)驗(yàn)可重復(fù)性，建議所有保存于RNAlater??Solution中樣本，都要置于-20 °C 或-80 °C 長(zhǎng)期保存。

1）顯微鏡下觀察細(xì)菌生長(zhǎng)狀態(tài)，盡量收集生長(zhǎng)期處于對(duì)數(shù)期的細(xì)菌。

2）將適量體積的菌液轉(zhuǎn)移至2 mL 旋蓋尖底離心管（無(wú)菌，無(wú)核酸酶）中，于室溫下14000 ×g 離心1 min。

3）棄掉培養(yǎng)基，將細(xì)菌菌體沉淀迅速置于液氮中冷凍1-3 h以上（凍存時(shí)間視組織量而定，保證樣品凍存充分），然后轉(zhuǎn)移至-80 °C長(zhǎng)期保存。

4）將樣品管置于管架上固定好后埋在干冰箱的中間部位，大體積干冰運(yùn)輸。

1）顯微鏡下觀察酵母菌生長(zhǎng)狀態(tài)，盡量收集生長(zhǎng)期處于對(duì)數(shù)期的酵母菌。

2）將適量體積的酵母菌液轉(zhuǎn)移至將2 mL旋蓋尖底離心管中（無(wú)菌，無(wú)核酸酶），于室溫下14000×g 離心1 min。

3）棄盡培養(yǎng)基，將酵母細(xì)胞沉淀迅速置于液氮中冷凍1-3 h以上（凍存時(shí)間視組織量而定，保證樣品凍存充分），然后轉(zhuǎn)移至-80 °C長(zhǎng)期保存。

4）將樣品管置于管架上固定好后埋在干冰箱的中間部位，大體積干冰運(yùn)輸

一次提取反應(yīng)所需細(xì)胞數(shù)≤1 X 10^7 ?個(gè)，以3 X 10^6-1 X 10^7 ?個(gè)為宜，可以將細(xì)胞經(jīng)裂解液裂解后凍存運(yùn)輸。

1、從培養(yǎng)箱中取出貼壁培養(yǎng)的細(xì)胞，顯微鏡下觀察細(xì)胞，確定生長(zhǎng)狀態(tài)良好（正常細(xì)胞融合度在80 %左右）;

2、棄去培養(yǎng)基，向細(xì)胞培養(yǎng)瓶或培養(yǎng)皿中加入 5mL PBS（RNase free）， 清洗一次后棄盡PBS；

3、加入1mlPBS（RNase free）重懸后將細(xì)胞轉(zhuǎn)移至 1.5 mL 旋蓋尖底離心管（RNase free）中；

4、離心（依據(jù)客戶實(shí)驗(yàn)室細(xì)胞離心步驟，注意細(xì)胞離心要適度，不要使細(xì)胞離心過實(shí)而造成裂解液不能充分滲透，4度離心機(jī)，轉(zhuǎn)速3000g為宜）得到細(xì)胞沉淀，棄去PBS，置于液氮中速凍2h后轉(zhuǎn)移至-80℃低溫保存，送樣時(shí)選擇干冰運(yùn)輸寄送。

離心收集的細(xì)胞迅速溶于 TRIzol 裂解，參考用量為每 5 X 10^6個(gè)細(xì)胞加 1 mL TRIzol；細(xì)胞溶于 TRIzol 后，如出現(xiàn)成團(tuán)，需用吸頭將細(xì)胞團(tuán)吹打散，或者激烈震蕩混勻，使細(xì)胞完全溶于 TRIzol 中充分裂解，室溫靜置5min，之后轉(zhuǎn)移至-80℃低溫保存，送樣時(shí)選擇干冰運(yùn)輸寄送。

液氮速凍法：

離心收集的細(xì)胞直接液氮速凍，速凍2h后轉(zhuǎn)移至-80℃低溫保存，送樣時(shí)選擇干冰運(yùn)輸寄送。

一次提取反應(yīng)所需細(xì)胞數(shù)≤1 X 10^7 ?個(gè)，以3 X 10^6-1 X 10^7 ?個(gè)為宜，可以將細(xì)胞經(jīng)裂解液裂解后凍存運(yùn)輸。

1、確定細(xì)胞生長(zhǎng)狀態(tài)良好；離心得到細(xì)胞沉淀（依據(jù)客戶實(shí)驗(yàn)室細(xì)胞離心步驟，注意細(xì)胞離心要適度，不要使細(xì)胞離心過實(shí)而造成裂解液不能充分滲透，4度離心機(jī)，轉(zhuǎn)速3000g為宜）；

2、棄去培養(yǎng)基，加入 1 mlLPBS（RNase free，室溫），輕輕將細(xì)胞沉淀懸起， 轉(zhuǎn)移至 1.5 mL 旋蓋尖底離心管（RNase free）中；

3、離心（依據(jù)客戶實(shí)驗(yàn)室細(xì)胞離心步驟，注意細(xì)胞離心要適度，不要使細(xì)胞離心過實(shí)而造成裂解液不能充分滲透，4度離心機(jī)，轉(zhuǎn)速3000g為宜）得到細(xì)胞沉淀，棄去PBS，置于液氮中速凍2h后轉(zhuǎn)移至-80℃低溫保存，送樣時(shí)選擇干冰運(yùn)輸寄送。

TRIzol 裂解法（僅限于RNA類產(chǎn)品）：

離心收集的細(xì)胞迅速溶于 TRIzol 裂解，參考用量為每 5 X 10^6個(gè)細(xì)胞加 1 mL TRIzol；細(xì)胞溶于 TRIzol 后，如出現(xiàn)成團(tuán)，需用吸頭將細(xì)胞團(tuán)吹打散，或者激烈震蕩混勻，使細(xì)胞完全溶于 TRIzol 中充分裂解，室溫靜置5min，之后轉(zhuǎn)移至-80℃低溫保存，送樣時(shí)選擇干冰運(yùn)輸寄送。

*注：判斷裂解液加入量是否合適的標(biāo)準(zhǔn)可以根據(jù)細(xì)胞溶解物的黏度來判斷。在細(xì)胞剛?cè)芙鈺r(shí)，可以發(fā)現(xiàn)有絲狀物出現(xiàn)，若裂解液加入量合適，吹打幾次后，絲狀物會(huì)消失，液體黏稠性下降；若裂解液的量過少，絲狀物往往一直存在，液體黏稠性大，應(yīng)繼續(xù)補(bǔ)加裂解液。裂解液加入量過少，會(huì)導(dǎo)致抽提的 RNA 降解。

收集的細(xì)胞直接液氮速凍后，轉(zhuǎn)移至-80℃低溫保存，送樣時(shí)選擇干冰運(yùn)輸寄送。

1、采血管選擇

采集血液進(jìn)行檢測(cè)面臨的一個(gè)主要問題就是細(xì)胞內(nèi) RNA 的不穩(wěn)定性，RNA 在采血后數(shù)小時(shí)內(nèi)便會(huì)迅速降解。此外，某些種屬的 RNA 在采血后會(huì)通過基因誘導(dǎo)在體外增加。體外RNA 降解和基因誘導(dǎo)均可導(dǎo)致體內(nèi)相關(guān)基因的轉(zhuǎn)錄本數(shù)目低估或高估。

BD 公司的 PAXgene 血液 RNA 管內(nèi)含能夠穩(wěn)定體內(nèi)基因轉(zhuǎn)錄性狀的添加劑，它能夠減少體外 RNA 的降解，并將基因誘導(dǎo)的水平降低到最小程度，適用于人和靈長(zhǎng)類動(dòng)物全血樣品制備，對(duì)于全血樣品，我們只推薦使用此管送樣。

BD PAXgene 血液 RNA 管（貨號(hào) 762165）

2、采集血液

使用PAXgene血液RNA管前請(qǐng)仔細(xì)閱讀產(chǎn)品說明書，以便進(jìn)行正確的操作。如果PAXgene血液 RNA 管為唯一的抽血管，則將血液抽入 PAXgene 血液 RNA 管之前應(yīng)先抽血入“廢棄管”內(nèi)，使抽血過程中所用的采血器內(nèi)能被血液預(yù)灌注；否則 PAXgene 血液 RNA 管應(yīng)為抽血程序中的最后一只試管。確保血液已停止流入試管再?gòu)某轴樒魃先∠略嚬埽≒AXgene 血液 RNA管內(nèi)的負(fù)壓設(shè)計(jì)可向管內(nèi)抽 2.5ml 的血）。

采血后應(yīng)立即將 PAXgene 血液 RNA 管輕輕地顛倒 8-10 次，以確保管內(nèi)保護(hù)劑和血液充分混合均勻。

3、保存、運(yùn)輸

將 PAXgene 血液 RNA 管豎直置于室溫（18-25℃）靜置 2-24 小時(shí)，之后可貯藏于-20℃或更低溫度；若試管要貯藏于低于-20℃的溫度，先在-20℃冷凍 24 小時(shí)，然后再將其轉(zhuǎn)移到-70℃或-80℃，可保存至少 50 個(gè)月。大體積干冰運(yùn)輸。

2、上下輕輕顛倒混勻十次后，分裝到常規(guī)EP管中，轉(zhuǎn)移至-20或-80 °C長(zhǎng)期保存（實(shí)際根據(jù)采血管的說明要求操作）

3、干冰運(yùn)輸

注：采集完血液之后需將血液轉(zhuǎn)移至EP管中，玻璃材質(zhì)采血管冷凍之后易碎

制備常見問題匯總

所有組織樣品務(wù)必根據(jù)質(zhì)量多少選擇使用下圖的螺紋管保存，嚴(yán)禁直接裝入密封袋或錫箔紙（低溫凍脆易破裂，導(dǎo)致樣品泄露，交叉污染）保存

組織送樣量要求

建議送樣量適用于95%的物種組織，少部分組織因胞內(nèi)核酸量較低，需適當(dāng)增加送樣量，比如植物果實(shí)、根等。送樣量過少或過多均不利于提取實(shí)驗(yàn)，下限為建議量的一半，上限為建議量的3倍。隨著技術(shù)發(fā)展，當(dāng)前大部分產(chǎn)品建庫(kù)起始量較低，不需要太多組織，采樣量過多反而會(huì)造成速凍不徹底、提取取樣困難、凍融降解等危險(xiǎn)。

1、RNA類產(chǎn)品建議組織送樣量

2、二代DNA類產(chǎn)品建議組織送樣量

3、三代DNA類產(chǎn)品建議組織送樣量