中國南方β地貧41-42（-CTTT）移碼突變進化歷程的基因轉(zhuǎn)換證據(jù).pdf

1、地中海貧血(簡稱地貧)是一種由于α和β珠蛋白合成比例失衡所導(dǎo)致的嚴重的致死性遺傳病，是世界上最常見的人類常染色體單基因遺傳病之一。地貧在全球的熱帶和亞熱帶地區(qū)有著廣泛的分布和相當高的攜帶率，最新的研究數(shù)據(jù)表明，該種疾病的基因攜帶者占全球人口總數(shù)的4.83％，而我國該病的高發(fā)區(qū)主要分布于長江以南的地處熱帶和亞熱帶的各省區(qū)。地貧可以根據(jù)珠蛋白生成障礙的種類分為α地貧和β地貧兩大類，分別是由α珠蛋白和β珠蛋白合成障礙而引起。在我國β地貧攜帶者

2、當中，β珠蛋白編碼區(qū)密碼子(Codon，CD)41/42位(CD41/42)的4堿基缺失(-CTTT)是最常見的致病突變，該突變在中國南方人群中攜帶率約為1％左右，占所有中國南方β地貧致病突變的40％左右。 地貧之所以在熱帶和亞熱帶有如此高的攜帶率，長期以來被認為是與瘧疾的自然選擇相關(guān)，早在20世紀40年代，英國遺傳學家JBSHaldane就提出了地貧是因為瘧疾的選擇壓力導(dǎo)致其流行的假說，類似的假說和觀點的產(chǎn)生不僅是因為地貧高發(fā)

3、區(qū)域的地理分布與瘧疾流行區(qū)域的地理分布有著很大相似性，而且地貧的攜帶者在流行病學和分子病理學水平被證明可以抵抗瘧疾。進入20世紀九十年代以來，隨著高通量DNA序列測定技術(shù)的發(fā)展，以及測序成本的降低，針對人群的基因組水平足夠大范圍的測序已成為現(xiàn)實，在群體遺傳學水平上，科研人員對葡萄糖-6-磷酸脫氫酶缺陷癥(G6PDDeficiency)、鐮刀型貧血癥(HbS)、血紅蛋白E(HbE)、血紅蛋白C(HbC)等血紅蛋白病與瘧疾的進化關(guān)系都展開了

4、初步的探索，并得到了一些關(guān)于瘧疾的自然選擇壓力在攜帶這些致病突變的染色體上所留下的基因組水平和群體遺傳學水平的證據(jù)，進一步讓我們確信前人關(guān)于這些疾病受到瘧疾的自然選擇的假說是正確的，但是在β地貧方面，還沒有群體遺傳學水平相關(guān)的研究報道。與此同時，亦沒有針對β地貧的分子進化領(lǐng)域的相關(guān)報道。 本研究預(yù)計通過針對CD41/42β地貧相關(guān)樣本的群體遺傳學研究，分析得到瘧疾的自然選擇壓力作用于CD41/42β地貧，從而在基因組水平上留下痕

5、跡的證據(jù)，達到證明CD41/42β地貧和瘧疾的自然選擇關(guān)系的目的。同時根據(jù)研究過程中得到的群體遺傳學數(shù)據(jù)來對CD41/42β地貧等位基因和單倍型在中國南方人群里的發(fā)生、發(fā)展和變化的進化歷程做一些初步的探索，以填補該遺傳疾病在分子進化領(lǐng)域的研究空白，并可以給今后對其他遺傳疾病的類似研究提供參考。 材料與方法： 我們選擇了以中國南方人群β地貧中攜帶率最高的一種突變--CD41/42作為研究靶點，以β珠蛋白基因簇作為基因組上的

6、研究區(qū)域，采集了24個正常人和中國南方九省的共16個CD41/42β地貧攜帶者作為研究樣本，抽提外周血基因組DNA，然后對該樣本每個個體β珠蛋白基因簇的ψβ-globin、δ-globin、RecombinationHotspot和β-globin區(qū)域分別進行長片段PCR擴增，并對經(jīng)純化的PCR產(chǎn)物進行雙向序列測定。上述四個片段總長共15.933kb，覆蓋20.693kb的基因組范圍。 測序結(jié)果經(jīng)DNAStarLasergene

7、軟件比對得到49個單核苷酸多態(tài)性位點(SingleNucleotidePolymorphism，SNP)，對所有位點用HWsim軟件進行了Harding-Weinberg平衡檢驗之后使用Phase2.1.1軟件對這49個SNP位點以及CD41/42位點進行單倍型推斷和重組熱點區(qū)的預(yù)測。在推斷得到的單倍型的基礎(chǔ)上，進行了如下的群體遺傳學分析：首先運用Arlequin3.01軟件對整個樣本進行了Tajima’sD中性檢驗，該結(jié)果是從群體遺傳

8、學的角度判斷研究區(qū)域是否受到自然選擇和受到何種自然選擇的依據(jù)：運用DnaSP4.0軟件對樣本進行四配子檢驗，用以測定在研究區(qū)域是否有重組發(fā)生和重組在何處發(fā)生；運用Network4.1.1.2軟件對樣本的單倍型進行了中介網(wǎng)絡(luò)圖的構(gòu)建，用以分析各種單倍型之間的親緣關(guān)系；同時用.Arlequin3.01軟件進行連鎖不平衡(LinkagcDisequilibrium，LD)計算并繪制連鎖不平衡圖譜，可以用來給之前各種分析提供參考。 把上

9、述各種分析得出的結(jié)果進行綜合總結(jié)和討論，可以對CD41/42β地貧等位基因的進化歷程和各種單倍型之間的進化關(guān)系得出一個較為可信的推測。 因為我們的樣本并非按照CD41/42的自然攜帶率進行采集，而是人工構(gòu)建的樣本(artificialsample)，所以我們在中性檢驗和LD計算的過程中都充分考慮了人工樣本給結(jié)果帶來的偏差，并在文中做出了討論。 結(jié)果與討論： 用DNAStar軟件對PCR產(chǎn)物雙向測序結(jié)果比對分析，共

10、找到49個SNP位點，MonteCarloPermutationtest的結(jié)果顯示，除去一個位點之外，其余位點均符合Harding-Weinberg平衡。Phase軟件預(yù)測出研究區(qū)域含有一個重組熱點區(qū)，這與前人的研究結(jié)果和我們之后的LD分析結(jié)果是一致的，為了避免重組熱點區(qū)帶來的不確定因素的干擾，我們僅僅對重組熱點區(qū)下游的包含CD41/42位點的β-globin區(qū)域共16個多態(tài)位點進行了單倍型推斷，共得到20種單倍型，其中包括16種正常單

11、倍型和4種突變單倍型。 Tajima'sD中性檢驗在所研究樣本的β-globin區(qū)域得到的是負值，且結(jié)果具有統(tǒng)計學意義，這個結(jié)果表明，該區(qū)域受到平衡選擇(balancingselection)的作用，而已知的能對該區(qū)域帶來這種作用的候選對象只有瘧疾。ψβ-globin、δ-globin、RecombinationHotspot三個區(qū)域的中性檢驗結(jié)果各不相同，且結(jié)果均不具有統(tǒng)計學意義，這也說明了重組熱點區(qū)確實會給我們的分析帶來不穩(wěn)

12、定的干擾?？紤]到我們的樣本是人工構(gòu)建的非隨機樣本，我們還把樣本中的突變?nèi)旧w條數(shù)按不同組合刪減到符合實際攜帶率的水平，并對由此產(chǎn)生的多個樣本也進行了Tajima'sD中性檢驗，而得出的結(jié)果與原始樣本的結(jié)果是一致的，這也說明我們的樣本在一定程度上是可以代表隨機樣本的。 四配子檢驗和LD分析的結(jié)果提示在β-globin區(qū)域內(nèi)部也發(fā)生過一定程度的重組，那么這些重組事件會否對CD41/42等位基因和單倍型的播散帶來影響?根據(jù)單倍型構(gòu)建的

13、中介網(wǎng)絡(luò)圖顯示中國南方人群在β-globin區(qū)域的單倍型分為三個大支，每個支都有其主要單倍型，而CD41/42β地貧單倍型(βCD41/42haplotype)只分布在其中兩個較為年輕的支上，在無限位點模型下，討論推斷CD41/42突變最早發(fā)生于兩支之一的主要單倍型上，然后通過與另一支的主要單倍型發(fā)生基因轉(zhuǎn)換(geneconversion)從而把CD41/42等位基因傳播到另一支上，從而產(chǎn)生新的遺傳背景完全不同的βCD41/42hapl

14、otype，同時我們還可以觀測到在我們這20種單倍型中，有多種單倍型都很可能是由于另外兩種單倍型之間發(fā)生基因轉(zhuǎn)換而產(chǎn)生的。這樣的研究結(jié)果表明基因轉(zhuǎn)換對于中國南方人群β-globin區(qū)域單倍型的多樣性和致病突變單倍型--βCD41/42單倍型的播散起著主要作用。 本研究得到的上述成果在群體遺傳學水平上證明了瘧疾與β地貧的自然選擇關(guān)系，并且通過對所研究樣本β-globin區(qū)域單倍型的分析得知基因轉(zhuǎn)換在βCD41/42單倍型的發(fā)展和播