趙　瑩， 趙麗亞， 晁天柱， 張　蓉， 邢正弘， 陳國(guó)強(qiáng)， 肖君華

（1. 上海西普爾-必凱實(shí)驗(yàn)動(dòng)物有限公司， 上海 201203；2. 上海實(shí)驗(yàn)動(dòng)物研究中心， 上海 201203； 3. 東華大學(xué)生物科學(xué)與技術(shù)研究所， 上海 201620）

染色體替換系的方案構(gòu)建及評(píng)價(jià)

趙瑩1，2， 趙麗亞1，2， 晁天柱3， 張蓉1，2， 邢正弘1，2， 陳國(guó)強(qiáng)1，2， 肖君華3

（1. 上海西普爾-必凱實(shí)驗(yàn)動(dòng)物有限公司， 上海 201203；2. 上海實(shí)驗(yàn)動(dòng)物研究中心， 上海 201203； 3. 東華大學(xué)生物科學(xué)與技術(shù)研究所， 上海 201620）

目的 為縮短復(fù)雜性狀相關(guān)基因精細(xì)定位時(shí)間，構(gòu)建染色體替換系（CSSs），并對(duì)構(gòu)建方案進(jìn)行評(píng)價(jià)。 方法 以供體小鼠作為母本、受體C57BL/6小鼠作為父本進(jìn)行雜交， 獲得N1代小鼠。將N1代雄性小鼠與C57BL/6雌性小鼠雜交， 獲得N2代小鼠。選擇N2代1號(hào)染色體非重組的雄性小鼠與C57BL/6雌性小鼠回交， 獲得N3代回交小鼠； 通過(guò)短串聯(lián)重復(fù)序列（STR）與單核苷酸多態(tài)性（SNP）對(duì)每代一定數(shù)量的子代進(jìn)行遺傳分析， 依次回交10代， 選擇N10代1號(hào)染色體未重組的樣本， 雌雄交配， 挑選1號(hào)染色體為純合且來(lái)自供體品系的子代小鼠， 用于繁殖純合后代。通過(guò)基因芯片掃描鑒定替換系的純合度，并與期望值比較。 結(jié)果 根據(jù)構(gòu)建的CSSs方案， 得到2個(gè)不同的CSSs， 通過(guò)基因芯片掃描得到的替換系純合度與期望值接近。 結(jié)論 構(gòu)建的CSSs方案可完成（部分）近交系小鼠品系目標(biāo)染色體替換， 供體品系背景基因組中99.59%以上已被受體C57BL/6小鼠替換，證明構(gòu)建的CSSs方案可行。

染色體替換系（CSSs）； 方案構(gòu)建； 評(píng)價(jià)； 非重組率

小鼠與人類共患許多受到遺傳調(diào)控的常見(jiàn)疾病，其中既有單基因參與調(diào)控的疾病，也有多基因參與調(diào)控的疾病，而多基因參與調(diào)控的疾病屬于復(fù)雜性狀。為了定位復(fù)雜性狀相關(guān)基因， 目前主要有三種方法： 第一，傳統(tǒng)方法，兩個(gè)品系通過(guò)回交N2代或雜交N2代初步鑒定相應(yīng)的數(shù)量性狀基因座（QTL）［1］， 然后通過(guò)構(gòu)建同類系鑒定相應(yīng)QTL［2］。最后， 通過(guò)精細(xì)定位鑒定QTL。連鎖定位和同類系構(gòu)建需要10～15代， 對(duì)于小鼠來(lái)說(shuō)大概需要3～5年， 對(duì)于大鼠來(lái)說(shuō)則需更長(zhǎng)時(shí)間［3］。第二， 構(gòu)建超大規(guī)模重組近交系（Recombinant inbred lines）， 將C57BL/6J、A/J、CAST/Ei、NOD/LtJ、NZO、WSB/Ei、PWK/Ph 與129S1/SvImJ等8個(gè)品系互相雜交20代以上［4］，形成上千個(gè)重組近交系，形成的重組近交系包含有來(lái)自8個(gè)近交系的基因。第三，染色體替換系（chromosome substitution strain， CSSs）［3］，CSSs主要是由一個(gè)受體品系與一個(gè)供體品系構(gòu)建，CSSs的一條染色體來(lái)自供體品系，其余都是以受體品系為背景［3，5］，通過(guò)CSSs與受體品系有限的回交或雜交，比較CSSs與背景品系在感興趣表型的差異，從而鑒定相應(yīng)染色體上的基因座。CSSs已經(jīng)用于鑒定幾個(gè)復(fù)雜性狀的QTL，例如，明暗場(chǎng)選擇、開(kāi)放區(qū)域選擇、血清中甾醇類和血清中部分氨基酸的濃度都與1號(hào)染色體有關(guān)［6］； 肥胖和動(dòng)脈粥樣硬化是由環(huán)境和遺傳因素引起的，發(fā)現(xiàn)與動(dòng)脈粥樣硬化相關(guān)的QTL有兩個(gè)，與體質(zhì)量有關(guān)的QTL有數(shù)十個(gè)，并找到三個(gè)與能量消耗有關(guān)的QTL和一個(gè)與食物吸收相關(guān)的QTL［7］； 睪丸癌是目前影響年輕男性健康的一種常見(jiàn)癌癥，已發(fā)現(xiàn)與睪丸癌有關(guān)的QTL定位于19號(hào)染色體［8］； 與焦慮有關(guān)QTL定位于1號(hào)、4號(hào)、6號(hào)、15號(hào)、17號(hào)染色體［9］； 與性發(fā)育有關(guān)的QTL定位于6號(hào)、13號(hào)染色體［10］； 與骨質(zhì)疏松有關(guān)的QTL定位于1號(hào)染色體［11］。不同數(shù)量性狀相關(guān)疾病的QTL位于不同染色體上。

小鼠1號(hào)染色體是目前復(fù)雜性狀研究成果最集中的染色體， 在已知的4 000多個(gè)QTL中， 有400多個(gè)被定位在1號(hào)染色體（informatics.jax.org）， 占總數(shù)的10%，雖然有如此多QTL被成功定位，但明確鑒定導(dǎo)致復(fù)雜性狀差異的基因卻只有20個(gè)左右［12］。其根本原因是小鼠遺傳資源的積累與研究導(dǎo)致了復(fù)雜性狀初步定位的蓬勃發(fā)展； 而另一方面，初步定位到QTL后對(duì)基因的精細(xì)定位卻舉步維艱［13］，因此構(gòu)建染色體替換系非常重要。

1　材料與方法

1.1動(dòng)物飼育

選擇C57BL/6作為受體品系，C3H/He、FVB/ N為供體品系。所有動(dòng)物均為SPF級(jí)，購(gòu)自上海西普爾-必凱實(shí)驗(yàn)動(dòng)物有限公司［SCXK（滬）2008-0016］。所有動(dòng)物飼養(yǎng)于上海西普爾-必凱實(shí)驗(yàn)動(dòng)物有限公司［SYXK（滬）2008-0058］，動(dòng)物實(shí)驗(yàn)遵守1988年《實(shí)驗(yàn)動(dòng)物管理?xiàng)l例》。

1.2構(gòu)建染色體替換系

以供體小鼠作為母本、受體C57BL/6品系小鼠作為父本進(jìn)行雜交，獲得N1代小鼠。將N1代雄性小鼠與C57BL/6品系雌性小鼠雜交，獲得N2代小鼠。通過(guò)構(gòu)建的STR與SNP分型方案鑒定N2代小鼠1號(hào)染色體的重組情況［14-16］，以雄性小鼠為分型對(duì)象，選擇1號(hào)染色體上所有位點(diǎn)都未重組的樣本； 將選擇到的雄性小鼠與C57BL/6品系雌性小鼠回交，獲得N3代回交小鼠； 依次回交10代［17］； 從而獲得2個(gè)不同的染色體替換系。對(duì)N10代染色體替換系雌雄小鼠分型，挑選1號(hào)染色體未重組的樣本，雌雄交配，鑒定其后代1號(hào)染色體的重組情況，選擇1號(hào)染色體為純合且來(lái)自供體品系的子代小鼠，即為得到的小鼠染色體替換系（圖1）。

圖1　構(gòu)建染色體替換系方案Figure 1 The strategy for constructing CSSs

1.3基因組背景純度鑒定

選擇含有1 449個(gè)SNP位點(diǎn)的芯片， 通過(guò)基質(zhì)輔助激光吸電離飛行時(shí)間質(zhì)譜（MALDI-TOF MS）技術(shù)， 對(duì)2個(gè)CSSs N6代和N10代進(jìn)行基因組背景純度鑒定，每個(gè)品系每個(gè)代數(shù)選擇一個(gè)樣本，共4個(gè)樣本，根據(jù)所得數(shù)據(jù)對(duì)該染色體替換方案進(jìn)行評(píng)估，該步實(shí)驗(yàn)由晶能生物技術(shù)有限公司完成。

1.4染色體替換系N2代的非重組率及替換系子代種群大小

在構(gòu)建CSSs 過(guò)程中，每個(gè)品系需要保證有5～8只雄性傳代小鼠，才能在理想時(shí)間內(nèi)獲得數(shù)量足夠大的子代群體。通過(guò)構(gòu)建的STR與SNP基因分型方法鑒定每個(gè)替換系的N2代小鼠1號(hào)染色體的重組情況，分型對(duì)象為雄性小鼠，選擇1號(hào)染色體上所有位點(diǎn)都未重組的樣本，非重組率為每個(gè)品系篩選出的樣本數(shù)占所有用于分型的樣本總數(shù)，替換系子代種群大小等于需要傳代的小鼠數(shù)量乘以2再除以非重組率。

1.5替換系子代小鼠體質(zhì)量稱量

從N1代到N10代， 對(duì)每個(gè)品系每代小鼠10日齡和20日齡子代小鼠稱重，考察體質(zhì)量隨代數(shù)的變化情況。

1.6統(tǒng)計(jì)學(xué)處理

2　結(jié)果

2.1構(gòu)建的2個(gè)染色體替換系及基因組背景純度

借助經(jīng)典遺傳學(xué)， 以供體做母本， 受體C57BL/6品系做父本，雜交得到N1子代，該N1代雄性小鼠的X染色體來(lái)自于供體，Y染色體來(lái)自于受體C57BL/6， 而N1代雌性小鼠的X染色體為雜合子，一條來(lái)自于供體， 一條來(lái)自于受體C57BL/6，因此N1代使用雄性小鼠傳代，即N1代雄性小鼠與受體C57BL/6品系回交，得到的N2代雄性小鼠的線粒體和X染色體都來(lái)源于受體C57BL/6背景，從N3代采用雄性小鼠或者雌性小鼠與受體C57BL/6回交傳代，得到子代的線粒體和性染色體都是來(lái)源于受體C57BL/6背景，最終達(dá)到替換1號(hào)染色體的目的。通過(guò)構(gòu)建的CSSs方案， 目前已成功培育了2個(gè)染色體替換系，2個(gè)CSSs 分別為： C57BL/6.C3H/ He-Chr1、C57BL/6.FVB/N-Chr1（見(jiàn)圖2）。替換系的線粒體和性染色體完全來(lái)自C57BL/6，其1號(hào)染色體已被不同的供體品系替換，同時(shí)背景基因組的純度已達(dá)到期望值。

利用1 449個(gè)高密度SNP基因芯片對(duì)2個(gè)染色體替換系N6代、N10代進(jìn)行基因組純度鑒定，C57BL/6.C3H/He-Chr1背景基因組N6代和N10代純度分別為98.90%和99.65%，C57BL/6.FVB/N-Chr1背景基因組N6代和N10代純度分別為99.38%和99.59%（表1）。

圖2　兩個(gè)染色體替換系示意圖Figure 2 The diagram of two chromosome substitution strains

2.2染色體替換系N2代的非重組率及替換系子代種群大小

通過(guò)構(gòu)建的STR與SNP基因分型方法對(duì)2個(gè)染色體替換系N2代小鼠進(jìn)行分型，篩選出1號(hào)染色體未重組的樣本（表2）。發(fā)現(xiàn)近交系來(lái)源小鼠1號(hào)染色體的平均非重組率為17.7%。其中， C57BL/6. C3H/He-Chr1替換系的非重組率為15.5%， C57BL/6. FVB/N-Chr1替換系的非重組率為23.8%。自N1代起選用雄性小鼠傳代， 每代需5～8只小鼠傳代， 結(jié)合本方案中每個(gè)替換系N2代1號(hào)染色體非重組率的實(shí)際值， C57BL/6.C3H/He-Chr1與C57BL/6.FVB/N-Chr1替換系每代分別需要65～104只和42～68只個(gè)體。同時(shí)，N2代的非重組率及實(shí)際種群數(shù)量也為每個(gè)替換系后續(xù)回交繁殖群體數(shù)量的控制提供了依據(jù)。

表1　基因芯片檢測(cè)結(jié)果Table 1 The result of microarray detection

2.3各品系替換過(guò)程中子代的體質(zhì)量分析結(jié)果

每個(gè)染色體替換系每代10日齡和20日齡小鼠的平均體質(zhì)量和標(biāo)準(zhǔn)差如表3、表4。

C57BL/6.C3H/He-Chr1和C57BL/6.FVB/N-Chr1替換系的N1代10日齡與20日齡體質(zhì)量均高于其它代數(shù)的10日齡與20日齡體質(zhì)量，且有顯著差異（P＜0.05）。從N2代到N10代， C57BL/6.C3H/He-Chr1 和C57BL/6.FVB/N-Chr1替換系每代雄鼠與雌鼠的10日齡與20日齡體質(zhì)量基本處于一條直線上，無(wú)代數(shù)反應(yīng)關(guān)系（P＞0.05）。

表2　各替換系N2代1號(hào)染色體的非重組率Table 2 The rate of non-recombination of N2 mice in each chromosome 1 substitution strain

表3　各品系替換過(guò)程中子代10 d的體質(zhì)量Table 3 The body weight for 10 days progeny mice in chromosome 1 substitution strains

表4　各品系替換過(guò)程中子代20 d的體質(zhì)量Table 4 The body weight for 20 days progeny mice in chromosome 1 substitution strains

2.4CSSs的方案評(píng)價(jià)

染色體替換過(guò)程中各代基因組背景純度的理論值可以由計(jì)算而得［10］，N1代可替換約50%的基因組，即N1代C57BL/6遺傳背景的純度理論值約為50%。以此類推， N2代到N10代C57BL/6遺傳背景的純度理論值分別為75.00%、87.50%、93.75%、96.88%、98.44%、99.22%、99.61%、99.81%、99.90%。N6代純度理論值已接近99%，并趨于穩(wěn)定， 因此選N6代與N10代進(jìn)行純合度檢測(cè)。根據(jù)SNP基因芯片掃描結(jié)果， C57BL/6.C3H/He-Chr1替換系N6代與N10代的純合度分別為98.90%、99.65%，C57BL/6.FVB/N-Chr1替換系N6代與N10代的純合度分別為99.38%、99.59%。2個(gè)CSSs N6代的實(shí)際純合度均高于理論值98.44%，2個(gè)CSSs N10代的實(shí)際純合度與理論值99.90%接近，證明文中構(gòu)建的染色體替換方案可行。

3　討論

在本研究中，主要關(guān)注CSSs的方案構(gòu)建，通過(guò)供體與受體的雜交后代不斷與受體回交，采用STR與SNP遺傳位標(biāo)鑒定子代的1號(hào)染色體重組情況，篩選出1號(hào)染色體未重組的樣本，在替換過(guò)程中首先解決線粒體、X染色體、Y染色體替換的難題，通過(guò)1號(hào)染色體的非重組率估計(jì)每個(gè)替換系子代的種群數(shù)量，并通過(guò)基因芯片掃描進(jìn)一步評(píng)價(jià)該方案，CSSs系的構(gòu)建需要大量的分型工作和飼養(yǎng)工作， 大概需要2～3年，以及大量的費(fèi)用，但是一旦形成，它們對(duì)于公眾的染色體研究就是一個(gè)可持續(xù)和非常有價(jià)值的資源。到目前為止，文獻(xiàn)中報(bào)道［3］已經(jīng)構(gòu)建的CSSs主要有： CSS C57BL/6J-chr （i）A/J，i代表1～21，CSS 129/Sv-chr19MOLF/Ei。

另外，有報(bào)道［18］利用亨廷頓舞蹈病的疾病修飾基因座評(píng)估CSSs， C57BL/6J-chr10A/J/NaJ（CSS10）與HdhQ111基因敲除小鼠C57BL/6J雜交，對(duì)F1子代檢測(cè)， 證實(shí)染色體替換系在體質(zhì)量上有明顯差異，HdhQ111基因敲除小鼠與肝和核內(nèi)包涵體形成中CAG基因穩(wěn)定性有關(guān)。另外Courtney 等［19］利用C57BL/6J-chr10A/J/NaJ（CSS14）染色體替換系來(lái)評(píng)估小鼠的內(nèi)在運(yùn)動(dòng)能力， 將QTL定位在14號(hào)染色體，并對(duì)C57BL/6J-chr10A/J/NaJ（CSS14）與C57BL/6J回交的F2代進(jìn)行連鎖分析， 在14號(hào)染色體上找到多個(gè)相關(guān)QTL， 證明了CSSs可大大縮短基因定位時(shí)間。

在染色體替換過(guò)程中，假定各種染色體交換之間互不干擾，根據(jù)染色體長(zhǎng)度，子代染色體的非重組率期望值公式為1/2e-d/100［3］， d表示染色體長(zhǎng)度， 小鼠1號(hào)染色體長(zhǎng)度約為112.5 cM， 則小鼠1號(hào)染色體的非重組率期望值為16%。但是各種染色體交換之間并不是完全無(wú)關(guān)，且各品系間遺傳背景也有差異，因此小鼠子代染色體的實(shí)際非重組率與期望值略有差異，捷克、美國(guó)及日本等科學(xué)家進(jìn)行亞種間染色體替換的研究顯示，供體為A/J， 受體為C57BL/6的1號(hào)染色體實(shí)際非重組率為16%［2］。本方案中小鼠1號(hào)染色體的平均非重組率為17.7%，與文獻(xiàn)報(bào)道結(jié)果一致。

染色體替換系群由于僅僅對(duì)比一條染色體對(duì)性狀的貢獻(xiàn)，基因組其它基因調(diào)控因素被完全均一化，因此是研究微效基因獨(dú)特資源。因此，該類資源的研發(fā)可加速人類疾病相關(guān)基因的研究步伐，具有重大科學(xué)應(yīng)用價(jià)值與良好社會(huì)效益。

［1］Lander ES， Botstein D. Mapping mendelian factors underlying quantitative traits using RFLP linkage maps［J］. Genetics，1989， 121（1）：185-199.

［2］Snell GD. Methods for the study of histocompatibility genes ［J］. J Genet， 1948， 49（2）：87-108.

［3］Nadeau JH， Singer JB， Matin A， et al. Analysing complex genetic traits with chromosome substitution strains［J］. Nat Genet， 2000， 24（3）：221-225.

［4］Peters LL， Robledo RF， Bult CJ， et al. The mouse as a model for human biology： A resource guide for complex trait analysis ［J］. Nat Rev Genet， 2007， 8（1）：58-69.

［5］Churchill GA， Airey DC， Allayee H， et al. The collaborative cross， a community resource for the genetic analysis of complex traits［J］. Nat Genet， 2004， 36 （11）：1133-1137.

［6］Singer JB， Hill AE， Burrage LC， et al. Genetic dissection of complex traits with chromosome substitution strains of mice ［J］. Science， 2004， 304（5669）：445-448.

［7］Spiezio SH， Amon LM， Mcmillen TS， et al. Genetic determinants of atherosclerosis， obesity， and energy balance in consomic mice［J］. Mamm Genome， 2014， 25（11-12）：549-563.

［8］Youngren KK， Nadeau JH， Matin A. Testicular cancer susceptibility in the 129.MOLF-Chr19 mouse strain： additive effects， gene interactions and epigenetic modifications［J］.Hum Mol Genet， 2003， 12（4）：389-398.

［9］Singer JB， Hill AE， Nadeau JH， et al. Mapping quantitative trait loci for anxiety in chromosome substitution strains of mice［J］. Genetics， 2005， 169（2）：855-862.

［10］ Krewson TD， Supelak PJ， Hill AE， et al. Chromosomes 6 and 13 harbor genes that regulate pubertal timing in mouse chromosome substitution strains［J］. Endocrinol， 2004， 145（10）：4447-4451.

［11］ Beamer WG， Shultz KL， Ackert-Bicknell CL， et al. Genetic dissection of mouse distal chromosome 1 reveals three linked BMD QTLs with sex-dependent regulation of bone phenotypes［J］. J Bone Miner Res， 2007， 22（8）：1187-1196.

［12］ Korstanje R， Paigen B. From QTL to gene： the harvest begins ［J］. Nat Genet， 2002， 31（3）：235-236.

［13］ Laurie CC， Nickerson DA， Anderson AD， et al. Linkage disequilibrium in wild mice［J］. PLoS Genet， 2007， 3（8）：e144.

［14］ Xu J， Zhao Q， Du P， et al. Developing high throughput genotyped chromosome segment substitution lines based on population whole-genome re-sequencing in rice［J］. BMC Genomics， 2010， 11（1）：656-669.

［15］ 晁天柱， 陳國(guó)強(qiáng)， 趙瑩， 等. 用于“野生小家鼠來(lái)源一號(hào)染色體替換系”構(gòu)建的PCR- LDR分型系統(tǒng)［J］. 中國(guó)實(shí)驗(yàn)動(dòng)物學(xué)報(bào)， 2011， 19（5）：372-376.

［16］ 趙麗亞， 趙瑩， 邢正弘， 等. 采用STR和SNP對(duì)3個(gè)1號(hào)染色體替換系篩選比較［J］. 實(shí)驗(yàn)動(dòng)物與比較醫(yī)學(xué)， 2012， 32 （5）：396-400.

［17］ Rogner UC， Avner P. Congenic mice： cutting tools for complex immune disorders［J］. Nat Rev Immunol， 2003， 3（3）：243-352.

［18］ Ramos EM， Kovalenko M， Guide JR， et al. Chromosome substitution strain assessment of a Huntington’s disease modifier locus［J］. Mamm Genome， 2015， 26（3-4）：119-130.

［19］ Courtney SM， Massett MP. Effect of chromosome substitution on intrinsic exercise capacity in mice［J］. Version 2. Fl000Res， 2014， 3（9）：1-15.

Construction and Evaluation of Chromosome Substitution Strains

ZHAO Ying1，2， ZHAO Li-ya1，2， CHAO Tian-zhu3， ZHANG Rong1，2，XING Zheng-hong1，2， CHEN Guo-qiang1，2， XIAO Jun-hua3
（1. Sino-British SIPPR/BK Laboratory Animal LTD.， CO， Shanghai 201203， China；2. Shanghai Laboratory Animal Research Center， Shanghai 201203， China；3. Institute of Biological Sciences and Biotechnology， Donghua University， Shanghai 201620， China）

Objective To shorten the fine-structure mapping time， construct and evaluate the programme of chromosome substitution strains （CSSs） . Methods The rst step required making hybrids between C57BL/6 male mice and donor female mice. N1 hybrids male mice were backcrossed to host C57BL/6 female mice. Progeny with a non-recombinant chromosome derived from donor strain， in this case chromosome 1， were identified in this and subsequent backcrosses. N2 Male mice were backcrossed to host C57BL/6 female mice at each generation. A certain number of progeny per generation was screened with short tandem repeat （STR） and single nucleotide polymorphism （SNP）. These mice were backcrossed to host C57BL/6 at each generation until tenth generation. At the tenth backcross generation males and females with the non-recombinant chromosome derived from donor mice were intercrossed. Progeny of this intercross that were homozygous （homosomic） for chromosome 1 were used to propagate the homosomic strain. The homozygosity of CSSs was identified through scanning the whole genomic by microarray to compare with expected proportion. Results We have constructed two CSSs by the programme that had been constructed. The actual homozygosity of progeny with microarray was approaching the expected. Conclusion The constructed CSSs confirmed the replacement of target chromosome and over 99.59% chromosome segments from host C57BL/6. This breeding program to generate a CSSs panel proved to be feasible.

Chromosome substitution strains （CSSs）； Construct the programme； Discuss；Rate of non-recombination

Q95-33

A

1674-5817（2015）06-0478-06

10.3969/j.issn.1674-5817.2015.06.010

2015-06-05

上海市科學(xué)技術(shù)委員會(huì)科研計(jì)劃項(xiàng)目（12140900400、13140900302）

趙瑩（1982-）， 女， 助研， 碩士， 從事群體遺傳學(xué)及藥物毒性安全評(píng)價(jià)研究。E-mail： 12414863@qq.com