王建忠，蘭俊，陸珍先，楊珺，呂華麗，潘志新，費曉云

?

基于SSR分子標記的桉樹父本鑒定技術研究

王建忠1，蘭俊1＊，陸珍先1，楊珺2，呂華麗1，潘志新1，費曉云1

(1.廣西國有東門林場，廣西 扶綏 532108；2.北京林業(yè)大學生物科學與技術學院，北京 100083)

利用篩選出的11對SSR引物進行SSR分子標記分析，通過比對子代與母本相同引物對應標記位點的等位基因配置，鑒別自由授粉過程中可能由母本自交產生的個體，并借助Cervus軟件應用最大似然值法和遺傳排除法，從154個基因型候選父本中鑒定出D2、D3、D4、D7、D11、D19和D20等7個半同胞子代的遺傳父本分別為同一試驗林內定植的尾巨桉DH18-04、DH18-20、DH32-29和DH241-2。研究對于今后進一步深入探究桉樹自交不親和與近交衰退機制，準確進行子代父本鑒定和雜交親本重建，以及品種選育等工作具有重要指導意義。

桉樹；父本鑒定；SSR

20世紀80年代，中澳技術合作東門桉樹示范林項目引進了大量桉樹基因資源，建立了亞洲最大的桉樹基因庫，并對遺傳材料進行了卓有成效的改良，特別是在桉樹人工雜交育種方面最為顯著。通過制訂雜交育種策略，培育出了一大批優(yōu)良桉樹雜交種，育種技術處于世界先進水平[1]。雜交育種過程中，如果缺乏相應繁育親本間的雜交系譜信息或者遺傳數(shù)據，就會增加近親交配的機率，導致后代遺傳相似度增大，致使遺傳基礎狹窄。因此，在進行桉樹人工雜交育種之前，首先應該明確親本遺傳信息，為親本選擇提供理論基礎。

通常情況下，親本分析主要采用觀察記錄和遺傳標記分析兩種方法開展[2]。由于育種工作中，往往可以通過觀察記錄種子的來源從而直接獲得母本的信息，因此親本分析以父本鑒定為主。對于木本植物而言，通過對每個子代及其候選父本進行遺傳標記分析，采用遺傳排除法或最大似然值法，可鑒定出候選父本。前者是基于孟德爾遺傳定律，即子代的等位基因一定來自雜交親本，根據多位點遺傳數(shù)據直接排除與后代具有不同位點的親本，從而逐次在候選父本中鑒定得到遺傳父本[3]。后者則最早由MEAGHER[4-5]提出，通過估計親本的似然性模型來預測最可能親本來源，可用于分析已知母本，或親本均不明確的子代對應的遺傳親本，是父本鑒定中最常用的方法。

SSR分子標記為共顯性標記且遵從孟德爾遺傳定律，利用Cervus軟件中已知性別群體鑒定模型[6-7]，并結合遺傳排除法比對子代基因型與母本的對應標記點的等位基因配置[8]，可以對自交或雜交子代進行判定[9]，并在大量候選父本中鑒定出雜交子代的遺傳父本[10-11]。CHAIX等[12]曾利用SSR分子標記對巨桉()種子園進行親本分析。PATTERSON等[13]則進一步通過鑒定1 954個藍桉()子代的遺傳父本，對種子園花粉授粉方式進行研究。GRATTAPAGLIA[14]采用SSR分子標記，對桉樹種子園中經自由授粉產生的優(yōu)良子代進行父本鑒定，通過計算親本特殊配合力來評估種子園中親本植株的優(yōu)越性，進而通過剔除對子代優(yōu)良性狀貢獻率低的親本實現(xiàn)對種子園的改建，從而建立一種向前選擇的策略來提高選育優(yōu)良品種的效率[10]。

本研究通過篩選得到的高多態(tài)性SSR引物，結合SSR分子標記分析和遺傳排除法，對前期誘導獲得的尾巨桉()半同胞材料的親本來源進行鑒定，從中鑒別出自交和雜交個體，并進一步對雜交個體的父本來源進行鑒定，從而建立一套桉樹子代父本鑒定及雜交親本重建的技術方法，為進一步開展桉樹遺傳改良工作奠定技術基礎。

1 材料與方法

1.1 材料來源

研究材料來源于2015年收集的尾巨桉優(yōu)樹經自然授粉結實得到的種子，經播種育苗得到的半同胞子代。隨機選擇其中20個基因型作為半同胞子代父本鑒定的實驗材料。

上述半同胞子代來源于共同的母本M(尾巨桉無性系DH18-19)，其定植于廣西國有東門林場基因庫試驗林。該試驗林建于1999年5月，共收集桉樹無性系(家系)200余個，包含東門雜交無性系/家系(DH系列編號)、尾葉桉()家系(U系列編號)、佛羅里達雜交桉家系(F系列編號)、巴西雜交桉家系(E系列編號)等，各無性系(家系)單行5株，株行距為2 m × 3 m(圖1)。因來源于自然授粉結實，所以這些子代基因型的父本未知，候選父本可能來源于母本M定植的試驗林內部，也可能來源于試驗林周邊行道樹或農戶種植的雜交桉無性系商品林。

圖1 廣西國有東門林場基因庫采樣示意圖

注：M為母本(尾巨桉無性系DH18-19)，空位為因無性系植株死亡導致的缺失

1.2 實驗材料采集、DNA提取與SSR分析

上述半同胞子代，以及包含有它們共同母本M在內的基因庫試驗林內的所有無性系(家系)均是用于父本鑒定研究材料的采集對象。該試驗林保護行所用品種尾巨桉DH32-29，試驗林周邊商品林造林品種尾巨桉DH33-9、DH32-26、DH32-28等，以及周邊道路行所用品種雷林1號(No.1)等候選父本均是材料樣本的采集對象。除子代和母本M以外，共采集到154個基因型的候選父本。采集時，選擇幼嫩新葉作為樣本材料，每個基因型的葉片樣本經編號后置入液氮內速凍后放入足量的干冰中保存、運輸，并盡快轉入-80℃超低溫冰箱進行長期保存。

葉片樣本在液氮內低溫研磨，使用植物基因組DNA提取試劑盒(天根生化科技(北京)有限公司，中國)進行總DNA提取，并經瓊脂糖凝膠電泳檢驗提取質量。共使用3種引物，具體包括5’端包含有M13序列(5’TGTAAAACGACGGCCAGT3’)的上游引物，普通下游引物，以及標有FAM(6-carboxyf-luorescein)、HEX(hexachloro-6-carboxyfluorescein)、ROX(6-carboxy-x-rhodamine)和TAMRA(tetrachloro-6-carboxyrhodamine)熒光標記的M13引物(Schuelke，2000)。PCR總反應體系為20 μL(2 μL DNA模板，10 μL PCR Master Mixes(北京擎科新業(yè)生物技術有限公司，中國)，7.2 μL ddH2O，0.08 μL正向引物，0.32 μL反向引物，0.4 μL熒光引物)，程序設計為：94℃持續(xù)4 min；94℃持續(xù)30 s，退火溫度設置為56℃持續(xù)30 s，72℃持續(xù)30 s，共計25次循環(huán)；94℃持續(xù)30 s，退火溫度設置為53℃持續(xù)30 s，72℃持續(xù)30 s，共計25次循環(huán)；72℃持續(xù)10 min；4℃避光保存。使用ABI-3730XL基因分析儀(Appli-ed Biosystems公司，美國)對PCR產物進行測定，測定結果使用GeneMarker(version 2.2.0)軟件(Soft-Genetics，USA)進行讀取和分析[15]。

1.3 多態(tài)性引物篩選

以已有文獻報道的240對桉樹SSR引物作為篩選對象[16]。針對上述引物，以采集到的所有桉樹基因型DNA為模板依次分別進行TP-M13-SSR PCR，篩選在不同基因型之間有良好多態(tài)性的SSR引物，并通過POPGENE (version 1.32)軟件(University of Alberta，Canada)計算所有多態(tài)性引物的等位基因位點數(shù)(NA)、有效等位基因數(shù)(NE)、觀察雜合度(HO)、期望雜合度(HE)等遺傳參數(shù)[17]。

1.4 自交個體的排除

桉樹雌雄同花，可能存在自交?；诿系聽栠z傳定律理論，通過比對子代與母本的等位基因配置，若在子代的等位基因配置中存在至少一個在母本的相應標記位點上未被發(fā)現(xiàn)的等位基因，則表明這一子代特有等位基因可能來源于其他親本，即可初步判定這一子代基因型并非來源于自交。

1.5 桉樹半同胞子代父本鑒定

本研究采用SSR多態(tài)性引物對已知母本的尾巨桉半同胞家系子代進行父本鑒定。使用Cervus(version 3.0.7)軟件(Field Genetics，UK)，利用已知性別的群體鑒定模型進行分析[6-7]，對半同胞家系中雜交子代的父本進行分析，通過似然比值(LOD)來篩選最似父本。其中，LOD值表示與隨機候選父本相比，疑似父本傳遞給子代基因的可能性的大小。LOD值越大，表明鑒定結果的準確性越高。當LOD值大于0時，表明所鑒定出的親本為遺傳親本的可能性大于假設親本。若鑒定過程中出現(xiàn)相同LOD值的候選父本時，則進一步采用△-統(tǒng)計對候選父本LOD值之間的差距進行分析，通過在85% ~ 95%的置信區(qū)間內運行10 000次模擬來確定△-臨界值[6]。同時，使用Cervus軟件計算親本排除率EP1和EP2，其中EP1表示在雙親未知情況下，在候選父本群體中隨機取樣，被排除為疑似父本的概率，EP2表示在母本已知情況下，在候選父本群體中隨機取樣，隨機父本被排除為疑似親本的概率[3,18-19]。

2 結果與分析

2.1 尾巨桉半同胞子代中自交個體的排除

經過對已有文獻報道的240對SSR引物的初步篩選[16]，其中具有良好多態(tài)性的96對SSR引物被用來鑒定尾巨桉半同胞子代是否源于自交。結果顯示(圖2)，在7個半同胞子代D2、D3、D4、D7、D11、D19和D20中存在至少一個在母本M中的相應標記位點上未發(fā)現(xiàn)的等位基因，表明這些子代特有等位基因來源于母本以外的其他親本，即這些子代基因型并非來源于母本M的自交。

圖2 尾巨桉半同胞子代與母本SSR等位基因配置示意圖(部分示意)

2.2 尾巨桉半同胞子代父本鑒定引物的篩選

以半同胞子代、試驗林及周邊采集到的共163個基因型DNA樣本為模板，篩選得到表1中的11對具有較高多態(tài)性的引物：EUCeSSR0037、EUCeSSR0166、EUCeSSR0176、EUCeSSR0224、EUCeSSR0465、EUCeSSR0539、EUCeSSR0592、EUCeSSR0606、EUCeSSR0703、EUCeSSR1001、EUCeSSR1063。上述引物被用于7個半同胞子代的父本鑒定，并通過POPGENE軟件對每個多態(tài)性位點的等位基因數(shù)(NA)、有效等位基因數(shù)(NE)、觀察雜合度(HO)、預期雜合度(HE)和多態(tài)性信息含量(PIC)進行估算，結果顯示：在163個基因型中每個位點檢測到3至16個等位基因；觀察和預期雜合度(HO和HE)分別為0.119至0.748(平均值為0.527)和0.311至0.814(平均值為0.618)；多態(tài)性信息含量為0.281至0.784(平均值為0.572)。通過Cervus軟件計算雙親未知情況下在候選父本群體中隨機取樣的排除率EP1，和已知母本情況下在候選父本群體中隨機取樣的排除率EP2，結果顯示11對引物對父本的排除率EP2大于99.8%，表明可以有效排除候選父本中的非遺傳父本。

表1 廣西國有東門林場基因庫內無性系群體SSR位點遺傳變異參數(shù)

注：HO為觀察雜合度，HE為預期雜合度；PIC為多態(tài)性信息含量，EP1為在候選父本群體中隨機取樣，被排除為疑似父本的概率；EP2為當母本已知時，在候選父本群體中隨機取樣，被排除為疑似父本的概率。

2.3 尾巨桉半同胞子代父本鑒定

篩選得到的11對SSR引物被用于7個半同胞子代的父本鑒定，經過Cervus軟件計算，結果顯示(表2)：在母本已知的情況下進行父本鑒定，7個基因型均未出現(xiàn)與母本對應SSR位點上相異的等位基因；7個半同胞子代的LOD值均為正值，介于1.34至5.09之間；7個基因型均未與最似父本出現(xiàn)位點錯配情況，因而鑒定得到唯一父本。

表2 Cervus軟件父本鑒定分析表　　

3 討論

由生殖生物學特性所決定，桉樹具有較高的可能性發(fā)生自交[20-22]。盡管不同種桉樹均有不同步的花期[22]、特異性的花絲或柱頭長度[23]、胚珠或合子早期敗育[24]等自交不親和機制來減少自交的發(fā)生，但桉樹在自然授粉時仍有較高比率的自交后代產生[25-26]。近年來，由于不需要繁瑣的去雄工作，桉樹SVP(Single Visit Pollination)[27]、OSP(One Stop Pollination)[28]等“一步授粉”技術被廣泛應用于雜交育種工作，該類方法雖然在蒴蓋開裂前提前授粉，但仍有接近5%自交個體的產生[13,29]。由于它們當中的一些個體在8個月后才表現(xiàn)出胸徑生長的劣勢[30]，因此在育種工作中，自交個體對后續(xù)的苗期評價和選擇工作產生著負面影響。SSR分子標記為共顯性標記且遵從孟德爾遺傳定律?；谧哟牡任换蛞欢▉碜云潆p親的定律，通過比對子代基因型與母本的對應標記位點的等位基因配置，可以對自交或異交子代進行判定，準確率可達97% ~ 99%[9-10]。本研究應用此方法，在20個半同胞基因型中，排除了自交個體，鑒定得到7個半同胞雜交子代。

在桉樹中，已被確認并定名的種、亞種或變種超過900個[31]，并且亞屬內各種之間可以自由雜交[32]。在開展桉樹遺傳改良工作時，通過明確現(xiàn)有優(yōu)良個體的親本來源并實現(xiàn)雜交親本重建，可有效提升雜交育種工作效率。SSR分子標記可以有效地分析子代個體某個位點的基因組成，進而快速的檢測其親本來源，并在桉樹、楊樹()等上得到應用[10-11,33]。本研究應用此方法，并結合采用位點逐一排除的遺傳排除法[8]，鑒定出同一母本M的7個半同胞子代分別來源于4個遺傳父本。其中構成半同胞子代D7、D11和D20的雄配子來源于定植于母本M附近的同一遺傳父本DH18-20，這是因為桉樹作為典型蟲媒傳粉植物且蜜源豐富，在實驗中觀察到該母本在散粉季節(jié)吸引了大量螞蟻吸食花蜜，并往返于同一或相鄰植株的不同花枝和蟻穴之間，這可能是在鑒定異交子代時發(fā)現(xiàn)了較多自交個體的原因；其余4個基因型的遺傳父本則與母本相距超過60 m，這是由于借助蜜蜂等飛行昆蟲，桉樹花粉傳播距離最遠可達數(shù)十公里[34-36]。

本研究利用篩選出的11對SSR引物進行SSR分子標記分析，通過比對子代與母本相同引物對應標記位點的等位基因配置，可用來鑒別自由授粉過程中可能由母本自交產生的個體；借助Cervus軟件應用最大似然值法可以對尾巨桉半同胞子代的遺傳父本進行鑒定，此基礎上結合遺傳排除法，鑒定出7個半同胞子代的遺傳父本分別為同一試驗林內定植的尾巨桉DH18-04、DH18-20、DH32-29和DH241-2。上述研究對于未來進一步深入探究桉樹自交不親和和近交衰退機制，準確進行子代父本鑒定和雜交親本重建，以及開展品種選育等工作具有重要指導意義。

[1] RICHARD P,李海佳,李宏偉,等.廣西桉樹人工雜交育種研究[J].廣西林業(yè)科學,2006,35(4):238-242.

[2] 陳小勇.自然植物種群的親本分析及其在生態(tài)學研究中的應用[J].生態(tài)學雜志,1999,18(2):30-35.

[3] POLLEGIONI P, WOESTE K, MUGNOZZA G S, et al. Retrospective identification of hybridogenic walnut plants by SSR fingerprinting and parentage analysis[J]. Molecular Breeding, 2009, 24(4): 321-335.

[4] MEAGHER T R. Analysis of paternity within a natural population of. 1. Identification of most-likely male parents[J].The American Naturalist, 1986, 128(2):199-215.

[5] MEAGHER T R, THOMPSON E. The relationship between single parent and parent pair genetic likelihoods in genealogy reconstruction[J].Theoretical Population Biology,1986, 29(1):87-106.

[6] MARSHALL T C, SLATE J, KRUUK L, et al. Statistical confidence for likelihood-based paternity inference in natural populations[J]. Molecular Ecology, 1998, 7(5): 639-655.

[7] JONES A G, SMALL C M, PACZOLT K A, et al. A practical guide to methods of parentage analysis[J]. Molecular Ecology Resources, 2010, 10(1): 6-30.

[8] ROBLEDO-A J J, GIL L. Patterns of pollen dispersal in a small population ofL. revealed by total- exclusion paternity analysis[J]. Heredity, 2005, 94(1): 13-22.

[9] REIZ C, BRETO M P, ASINS M J. A quick methodology to identify sexual seedlings in citrus breeding programs using SSR marker[J]. Euphytica, 2000, 112(1): 89-94.

[10] GRATTAPAGLIA D, RIBEIRO V J, REZENDE G D S P. Retrospective selection of elite parent trees using paternity testing with microsatellite markers: an alternative short term breeding tactic for[J]. Theoretical Applied Genetics, 2004, 109(1): 192-199.

[11] HAN Z, GAO P, GENG X, et al. Identification of the male parent of superior half-sibindividuals based on SSR markers[J].Molecular Breeding,2017, 37(12):155-166.

[12] CHAIX G, GERBER S, RAZAFIMAHARO V, et al. Gene flow estimation with microsatellites in a Malagasy seed orchard of[J]. Theoretical Applied Genetics, 2003, 107(4): 705-712.

[13] PATTERSOM B, GORE P, POTTS B M, et al. Advances in pollination techniques for large-scale seed production in[J]. Australian Journal of Botany, 2004, 52(6): 781-788.

[14] GRATTAPAGLIA D. Molecular breeding ofMolecular biology of woody plants[M].Netherlands: Springer, 2000.

[15] HULCE D, LI X, SNYDER-L T.GeneMarker? genotyping software: Tools to increase the statistical power of DNA fragment analysis[J]. Journal of Biomolecular Techniques:JBT, 2011, 22(Suppl): 35-36.

[16] ZHOU C, HE X, LI F, et al. Development of 240 novel EST-SSRs inL’Herit[J]. Molecular Breeding, 2014, 33(1): 221-225.

[17] YEH F, YANG R, BOYLE T. POPGENE Version 1.32 Microsoft Windows-based freeware for populations genetic analysis[M]. Edmonton: University of Alberta, 1999.

[18] ISAGI Y, KANAZASHI T, SUZUKI W, et al. Highly variable pollination patterns in Magnolia obovata revealed by microsatellite paternity analysis[J]. International Journal of Plant Sciences, 2004, 165(6): 1047-1053.

[19] VAN E A L, WEABER R L, DRAKE D J, et al. DNA-based paternity analysis and genetic evaluation in a large, commercial cattle ranch setting[J]. Journal of Animal Science, 2007, 85(12): 3159-3169.

[20] GRIFFIN A R, HAND F C.Post-anthesis development of flowers ofF. Muell and the timing of artificial pollination[J].Australian Forest Research,1979,9: 9-15.

[21] ELDRIDGE K G, DAVIDSON J, HARDWOOD C, et al. Eucalypt domestication and breeding[M].Oxford: Clarendon Press, 1993.

[22] POTTS B M, GORE P L. Reproductive biology and controlled pollination of–a review[M]. Tasmania: School of Plant Science UTAS, 1995.

[23] LARCOMBE M J, COSTA J, TILYARD P, et al. On the persistence of reproductive barriers in Eucalyptus: The bridging of mechanical barriers to zygote formation by F1 hybrids is counteracted by intrinsic post-zygotic incompatibilities[J]. Annals of Botany, 2016, 118(3): 431-444.

[24] POUND L M, WALLWORK M A B, POTTS B M, et al.. Early ovule development following self- and cross-pollinations inLabill. ssp. globules[J]. Annals of Botany, 2002, 89(5): 613-620.

[25] HARDNER C M, POTTS B M. Postdispersal selection following mixed mating in[J]. Evolution, 1997, 51(1): 103-111.

[26] COSTA E S J, HARDNER C, TILYARD P, et al. The effects of age and environment on the expression of inbreeding depression in[J]. Heredity, 2011, 107(1): 50-60.

[27] WILLIAMS D R, POTTS B M, BLACK P G. Testing single visit pollination procedures forand[J]. Australian Forestry, 1999, 62(4): 346-352.

[28] HARBARD J L, GRIFFIN A R, ESPEJO J. Mass controlled pollination of Eucalyptus globulus: a practical reality[J].Canadian Journal Forestry Research,1999, 29(10): 1457-1463.

[29] HORSLEY T N, JOHNSON S D, MYBURG A A. Comparison of different control-pollination techniques for small-flowered eucalypts[J]. New Forests, 2010, 39(1): 75-88.

[30] HARDNER C M, Potts B M. Inbreeding depression and changes in variation after selfing inssp. globules[J]. Silvae Genetica, 1995, 44: 46-54.

[31] DOUGHTY R W. The Eucalyptus. A natural and commercial history of the gum tree[M]. Baltimore: John Hopkins Univ Press, 2000.

[32] WILCOX M D. A catalogue of the eucalypts[M]. Auckland: Groome Poyry, 1997.

[33] DONG C B, SUO Y J, KANG X Y. Assessment of the genetic composition of triploid hybridusing SSR markers with low recombination frequencies[J]. Canadian Journal Forestry Research, 2014, 44(7): 692-699.

[34] BARBOUR R C, POTTS B M, VAILLANCOUNT R E. Gene flow between introduced and native Eucalyptus species: crossability of native Tasmanian species with exotic E. nitens[J]. Australian Journal of Botany, 2005, 53(5) :465-477.

[35] BARBOUR R C, POTTS B M, Vaillancourt RE. Gene flow between introduced and nativespecies: Early-age selection limits invasive capacity of exotic×F1 hybrids[J]. Forest Ecology and Management, 2006, 228: 206–214.

[36] BARBOUR R C, POTTS B M, VAILLANCOUNT R E. Gene flow between introduced and native Eucalyptus species: Morphological analysis of tri-species and backcross hybrids involving[J]. Silvae Genetica, 2007, 56(3/4) :127-133.

Study on Using SSR Markers for Paternal Identification in

WANG Jianzhong1, LAN Jun1, LU Zhenxian1, YANG Jun2, LV Huali1, PAN Zhixin1, FEI Xiaoyun1

(1.,,,; 2.,)

A technique for male parent identification using SSR markers inhalf-sib progenies was developed in the study reported here. Eleven selected SSR loci with high polymorphism were used to identify the selfed individuals in an open pollinated population, based on the comparison of marker loci with the same primers between progenies and female parents. According to maximum likelihood scores obtained from analyses by Cervus software,clones DH18-04, DH18-20, DH32-29 and DH241-2 were identified as the male parents of the outcrossed progenies D2, D3, D4, D7, D11 and D20, respectively. This research has important significance for guiding further examination of self-incompatibility, paternal identification and hybrid parental reconstruction in these eucalypt taxa.

; paternal identification; SSR

Q943

A

廣西創(chuàng)新驅動發(fā)展基金項目(桂科AA17204087-3)

王建忠(1985― )，男，碩士，工程師，主要從事桉樹遺傳改良及無性系開發(fā)研究，E-mail: 379760245@qq.com

蘭俊(1981― )，男，碩士，高級工程師，主要從事桉樹遺傳改良及無性系開發(fā)研究，E-mail: lanjun107@aliyun.com