一種用於大豆遺傳完整性分析的繁殖更新方法與流程

2023-09-10 08:42:20

本發明屬於作物種質資源學領域，涉及一種用於大豆遺傳完整性分析的繁殖更新方法。

背景技術：

大豆(Glycine max(L.)Merr.)屬於豆科(Leguminosae)，蝶形花亞科(Papilionatae)，大豆屬(Glycine)，Soja亞屬，在我國人民日常生活和國民經濟體系中具有舉足輕重的地位，是人類重要的營養和保健食品。大豆籽粒的蛋白質含量為40％左右，約為其他植物的1～6倍，除蛋氨酸和胱氨酸含量較低外，其他人體必需的胺基酸含量都較高。因此，大豆蛋白對穀物蛋白有較好的增補效果，是僅次於動物蛋白的理想食品。大豆種質資源是大豆新品種選育和生物技術研究的重要物質基礎，也是大豆生產持續發展的基本保證。我國是世界保存大豆種質資源數量最多的國家，約35000餘份，這些材料是我國大豆遺傳研究和作物育種的寶貴財富。

遺傳完整性(genetic integrity)是指群體的遺傳結構得到完全的保持，包括基因型頻率分布和等位基因頻率分布和其原始群體一樣，保持不變。維持種質遺傳完整性就是種質在貯藏及繁殖更新過程中保持子代與親代具有最大程度的遺傳相似性，這是種質保存工作的核心。影響種質遺傳完整性的因素甚多，其中與繁殖更新相關的因素較多，例如，繁殖地點、繁殖群體量、授粉和收穫方式等。因此，必須對作物的繁殖更新技術方法加以規範，以確保最大程度地維持其遺傳完整性。

目前，國內外對於大豆種質遺傳完整性研究報導較少，且其繁殖更新的程序方法也不規範。大豆種質在繁殖更新時因其生活力下降、繁殖地點、繁殖群體量、授粉和收穫方式等不同，而使其遺傳完整性喪失。因此，針對大豆在繁殖更新時容易發生遺傳完整性變化的實際情況，制定出一套科學合理的用於大豆種質遺傳完整性分析時的繁殖更新方法已迫在眉睫。

技術實現要素：

本發明針對現有技術的不足，提供一種用於大豆遺傳完整性分析的繁殖更新方法，該方法完整可行且系統規範，可應用此方法對庫存的高生活力大豆進行繁殖更新以維持其遺傳完整性。

為此，本發明提供一種用於大豆遺傳完整性分析的繁殖更新方法，包括親代繁殖更新和子代繁殖更新。

所述親代繁殖更新：將親代種子經過人工老化處理後進行生活力檢測，按照發芽率分組，選取不同發芽率梯度的親代處理；各處理種子先進行育苗，然後移栽至大田中；單株收穫。

所述子代繁殖更新：將各個子代處理，按照與親代處理同樣的步驟進行育苗、移栽、形態性狀調查、單株收穫。

所述子代單株與其親代是一一對應關係，即子代的各個處理每個單株的編號與其親代單株是一一對應的。

親代和子代繁殖更新期間，進行形態標記分析及SSR分子標記分析以檢測其遺傳完整性。形態標記分析時，從親代和子代的每個處理中按照編號選出30個單株，對其進行11個形態性狀調查。進行性狀調查的30個單株，親代與子代的編號一一對應。同樣，進行SSR標記分析時親代及子代均取96個單株(包括形態性狀調查的30個單株)，子代採集葉片的單株與其親代編號也是一一對應關係。

具體地，本發明用於大豆遺傳完整性分析的繁殖更新方法，步驟如下：

(一)親代繁殖更新

1、繁殖地點選擇：

選擇種質原產地或與原產地生態環境條件相似地區的地勢平坦、地力均勻、形狀規則、排灌方便且最少兩年未種植大豆的田塊；

2、種子老化及生活力檢測

(1)種子人工老化：將播種、收穫整齊度一致的同一份材料的大豆種子，放入人工氣候箱中進行種子含水量平衡，平衡條件為45％RH，25℃，15天，使處理種子的含水量處於7.4％±0.1％；然後用鋁箔袋真空密封分裝，平均分成若干份，放置於人工老化箱內40℃恆溫老化；放置順序採用倒序法，即人工老化時間最長的處理先放入進行老化，人工老化時間最短的最後放入進行老化，試驗結束時一起取出；人工老化結束後在25℃下，將所有處理密封平衡2天。以未經過人工老化的處理為對照。

(2)生活力檢測：參考農作物種子檢驗規程GB/T3543.4中的標準條件，將對照與經過不同老化時間的處理進行發芽試驗，並統計發芽勢、發芽率、發芽指數和活力指數。

3、育苗及移栽：

按照遺傳完整性分析的試驗要求，選取包括對照在內的4～5個適合的發芽率梯度的大豆及野生大豆親代處理進行育苗。

(1)製備育苗培養基質：將砂壤土、草炭土、蛭石和珍珠巖按照體積比為4：3：2：1比例混合，向混合基質中加入多菌靈，比例為每kg基質加入1g多菌靈；加入複合肥(N：P2O5：K2O＝15：15：15)，比例為每kg基質加入5g複合肥，加水攪拌混勻後，製得育苗培養基質。

(2)播種定盤：即將種子播種於育苗盤的過程。

(3)溫室育苗：將步驟(2)中的育苗盤放入溫室中，其上覆蓋地膜保溼，溫度控制在24.5～25.5℃，相對溼度為74.5～75.5％，種子出苗後將地膜掀掉。待幼苗長至三葉期後即可進行移栽。

(4)移栽大田：畝施農家肥1000～2000kg，磷酸二銨15kg。大豆每個處理田間定植200株，行距為50cm，株距為30cm，挖移栽穴，深10cm，將育苗盤中的幼苗移栽到大田中，及時填土澆水。

4、田間管理

(1)掛牌編號：待幼苗長至50cm高時，將定植於大田中每個處理的植株進行單株掛牌編號。

(2)栽培措施：常規管理即可。

5、形態性狀調查及去雜

繁殖更新期間，從每個處理中按照編號選出30個單株，對其進行11個形態性狀調查，這些性狀包括：花色、粒色、粒形、臍色、茸毛色、莢色、結莢習性、株高、單株粒數、單株莢數和每莢粒數。並除去雜株。

6、收穫、脫粒、乾燥和保存

每個處理按照編號進行單株收穫，脫粒，乾燥至種子含水量8.0％以下，清選後保存。

(二)子代繁殖更新

1、繁殖地點選擇：將大豆不同生活力的親代各個處理繁殖更新所收穫的種子次年原地進行繁殖更新。子代處理的繁殖的試驗地及配套條件與其親代一致。

2、生活力檢測：子代處理繁殖更新時不需要經過種子老化。參考農作物種子檢驗規程GB/T3543.4中的標準條件，將各個子代處理進行發芽試驗，並統計發芽勢、發芽率、發芽指數和活力指數。

3、育苗及移栽：為保證與親代各個處理的生長環境條件及種植密度一致，對子代處理進行育苗及移栽。培養基質、播種定盤、溫室育苗和移栽大田過程與親代一致。子代各個處理在育苗時，每格播種的種子是其親代處理每個單株收穫的種子，即育苗時每格長出的幼苗編號與其親代單株是一一對應關係。每格播種2粒種子以保證成苗率。移栽時每格優選健壯的幼苗進行移栽。

4、田間管理：子代各個處理田間管理過程與其親代一致。

5、形態性狀調查及去雜：子代各個處理調查的形態性狀與親代完全一致。子代各個處理調查的每個單株與其親代的每個單株也是一一對應關係，同樣也是調查30個單株。去雜過程與親代一致。

6、收穫、脫粒、乾燥和保存：子代處理的收穫、脫粒、乾燥和保存等程序與其親代一致。

本發明提出了一整套完整可行且系統規範的用於大豆遺傳完整性分析時的繁殖更新方法，特別是提出將不同生活力的親代處理進行育苗移栽，對親代處理的單株進行編號，並進行單株形態性狀調查與收穫，子代繁殖更新及形態性狀調查時與其親代處理是一一對應關係等技術創新。綜合採用形態標記與SSR分子標記相結合的分析方法對提出的繁殖更新方法與常規的繁殖更新方法進行了遺傳完整性對比分析，從而進一步驗證了本發明提出的繁殖更新方法的可靠性和準確性。

本發明關鍵技術點及有益效果如下：

1、本發明提出的繁殖更新方法重點之一在於繁殖更新時子代單株與其親代單株的一一對應關係，即子代育苗時的各個處理每個單株的編號與其親代單株是一一對應的；進行形態性狀調查時，子代各個處理調查的單株與其親代單株也是一一對應關係，其目的是為了在進行大豆及野生大豆遺傳完整性分析時明確子代與其親代的遺傳對應關係，使得出的試驗數據具有更高的可靠性和準確性，得出的試驗結論具有更強的說服力。

2、將不同生活力的親代處理進行育苗移栽：由於大豆種子經過人工老化後，其種子活力較差，若是直播種子，田間出苗率很低，不能保證每個處理的種植密度一致，因此本發明採用先育苗，然後移栽的方法，每個處理定植200株，以保證種植密度一致，使得形態性狀調查數據可靠準確，避免了因種植密度不一致而造成形態性狀調查數據失真。

附圖說明

圖1為本發明大豆親代繁殖更新與子代繁殖更新時單株的一一對應關係示意圖。

其中A表示進行人工老化；B表示田間繁殖過程(子代單株與親代單株是一一對應關係)。

具體實施方式

下面結合具體試驗方法和附圖對本發明的技術方案及其所產生的技術效果進行進一步的闡述，下述說明僅是為了解釋本發明，但不以任何方式對本發明加以限制，基於本發明教導所作的任何變換或替換，均屬於本發明的保護範圍。

本發明中所使用的方法如無特殊說明，均為本領域常規方法。下述實施例中所用的試驗材料、試劑等，如無特殊說明，均可從商業途徑得到。

實施例一應用本發明方法和常規方法對大豆進行繁殖更新並進行遺傳完整性分析

(一)試驗材料：以大豆地方品種大黃豆種子為試驗材料。

(二)試驗方法：

1、方法一：利用本發明提供的方法對大豆地方品種大黃豆進行繁殖更新。步驟如下：

(1)大黃豆親代處理繁殖更新

①繁殖地點選擇：

a、繁殖地區：選擇種質原產地。由於地方品種大黃豆原產地是山東省濟寧市梁山縣，因此在當地找到一農戶，租賃其土地進行繁殖更新，能夠滿足繁殖更新材料的生長發育及其性狀的表達。

b、試驗地：地勢平坦、地力均勻、形狀規則、排灌方便且兩年內未種植大豆；遠離汙染源、無人畜侵擾、附近無高大建築；避開了病蟲害多發區、重發區和檢疫對象發生區。試驗地塊為每個小區寬1.5m，長30m，一個小區內進行一個處理的繁殖更新。

c、配套條件：具備了育苗、收穫、晾曬、貯藏等試驗條件和設施。

②種子老化及生活力檢測

a、種子人工老化：將播種、收穫整齊度一致的大豆地方品種大黃豆親代種子10000粒，放入人工氣候箱中進行種子含水量平衡，平衡條件為45％RH，25℃，15天，使處理種子的含水量處於7.4％±0.1％；然後用鋁箔袋真空密封分裝，平均分成若干份，每份1000粒，放置於人工老化箱內40℃恆溫老化；放置順序採用倒序法，即人工老化時間最長的處理先放入進行老化，人工老化時間最短的最後放入進行老化，試驗結束時一起取出；人工老化結束後在25℃下，將所有處理密封平衡2天。以未經過人工老化的處理為對照。

b、生活力檢測：參考農作物種子檢驗規程GB/T3543.4中的標準條件，將對照與經過不同老化時間的處理進行發芽試驗，並統計發芽勢、發芽率、發芽指數和活力指數(表1)。

③育苗及移栽：

按照遺傳完整性分析的試驗要求，選取包括對照在內的4個適合的發芽率梯度的親代處理，即處理GP-CK1、GP-I1、GP-II1和GP-III1進行育苗。由於大豆種子經過人工老化後，其種子活力較差，若是直播種子，田間出苗率很低，不能保證每個處理的種植密度一致，因此需要預先進行育苗，然後移栽到大田中，每個處理定植200株，以保證種植密度一致，使得形態性狀調查數據可靠準確，否則會因種植密度不一致而造成形態性狀調查數據失真。育苗步驟如下：

a、製備育苗培養基質：將砂壤土、草炭土、蛭石和珍珠巖按照體積比為4：3：2：1比例混合，向混合基質中加入多菌靈，比例為每kg基質加入1g多菌靈；加入史丹利複合肥(N：P2O5：K2O＝15：15：15)，比例為每kg基質加入5g史丹利複合肥，加水混合後攪拌混勻後，製得育苗培養基質。

b、播種定盤：即將種子播種於育苗盤的過程。利用5×10格育苗盤育苗，將步驟a製備的培養基質填充於育苗盤中，填實抹平後，用噴壺均勻灑水，1～2小時後，用圓頭小木棍或鑷子等器具在每個育苗盤格子上打孔，直徑1cm，深度2cm。將不同生活力處理的大豆親代種子播種於小孔中。發芽率在80％以上的處理，每格播種1粒；發芽率在50～80％的處理，每格播種2粒；發芽率在30～50％的處理，每格播種3粒；發芽率在30％以下的處理，每格播種4粒。播種後用培養基質覆蓋壓實，再次灑水澆透即可。共育苗4盤。

c、溫室育苗：將步驟b中的育苗盤放入溫室中，其上覆蓋地膜保溼，溫度控制在24.5～25.5℃，相對溼度為74.5～75.5％，待種子出苗後將地膜掀掉。帶幼苗長至三葉期後即可進行移栽。

d、移栽大田：幼苗移栽前，需造墒，精細整地，重施底肥，一般畝施農家肥1000～2000kg，磷酸二銨15kg。大豆每個處理田間定植200株。行距為50cm，株距為30cm，挖移栽穴，深10cm，將育苗盤中的幼苗移栽到大田中，及時填土澆水。移栽最好在下午進行，提高幼苗成活率。

④田間管理

a、掛牌編號：待幼苗長至50cm高時，將定植於大田中每個處理的植株進行單株掛牌編號。

b、栽培措施：封壟前及時進行中耕除草、培土；加強花莢期管理，爭取多花多莢，防花莢脫落，花莢期施肥可根據長勢酌情處理，如長勢弱，花期可適當追肥，壯苗不追肥防止徒長，花莢期和鼓粒期可要適當澆水，防止受旱，以免影響百粒重及產量。需要注意地是每個處理的澆水施肥量必須相同，以免產生系統誤差。適時除草，適期防病，科學治蟲。播後芽前進行土壤封閉除草；盛花至結莢鼓粒期造橋蟲、大豆卷葉螟等害蟲容易發生，如田間有低齡幼蟲啃食的網狀和鋸齒狀葉片出現，要及時用藥物防治。

⑤形態性狀調查及去雜

a、性狀調查：在地方品種大黃豆親代處理繁殖更新期間，從每個處理中按照編號選出30個單株，對其進行11個形態性狀調查，這些性狀包括：花色、粒色、粒形、臍色、茸毛色、莢色、結莢習性、株高、單株粒數、單株莢數和每莢粒數。

b、去雜：在形態性狀調查期間，核對每個單株的生育期、花色、茸毛色、葉形、生長習性、結莢習性和株高等主要性狀與主體類型明顯不一致的個體，當作雜株拔除。

⑥收穫、脫粒、乾燥和保存

a、收穫：適時收穫。對每個處理按照編號進行單株收穫，每個單株的種子收穫後進行編號，以便繁殖其子代。

b、脫粒：每個處理每個單株收穫的種子在脫粒前，必須打掃乾淨脫粒場地、機械、用具等，嚴防混雜；將每個單株的種子單粒脫，單袋裝；種子袋標籤編號須與每個處理的每個單株編號一致，袋內外各附標籤，避免寫(掛)錯標籤。

c、乾燥：脫粒裝袋後及時晾曬，防止發熱黴變。乾燥至種子含水量8.0％以下，可停止乾燥。

d、清選：去除破損粒、病蟲粒和泥沙等雜質。

e、保存：將清選過的種子用鋁箔袋密封包裝後保存於-4℃中期庫內。

(2)大黃豆子代處理繁殖更新

①繁殖地點選擇：將大豆地方品種大黃豆的親代各個處理繁殖更新所收穫的種子次年原地進行繁殖更新。子代處理的繁殖的試驗地及配套條件與其親代一致。

②生活力檢測：子代處理繁殖更新時不需要經過種子老化。子代處理的繁殖更新程序與其親代大體一致，但不需要經過種子老化。參考農作物種子檢驗規程GB/T3543.4中的標準條件，將各個子代處理，即處理GF1-CK1、GF1-I1、GF1-II1和GF1-III1進行發芽試驗，並統計發芽勢、發芽率、發芽指數和活力指數(表1)。

表1大豆地方品種大黃豆發芽數據統計(方法一)

③育苗及移栽：為保證與親代各個處理的生長環境條件及種植密度一致，對子代處理進行育苗及移栽。培養基質、播種定盤、溫室育苗和移栽大田過程與親代一致。子代各個處理在育苗時，每格播種的種子是其親代處理每個單株收穫的種子，即育苗時每格長出的幼苗編號與其親代單株是一一對應關係。每格播種2粒種子以保證成苗率。移栽時每格優選健壯的幼苗進行移栽。

④田間管理：子代各個處理田間管理過程與其親代一致。

⑤形態性狀調查及去雜：子代各個處理調查的形態性狀與親代完全一致。子代各個處理調查的每個單株與其親代的每個單株也是一一對應關係，同樣也是調查30個單株。去雜過程與親代一致。

⑥收穫、脫粒、乾燥和保存：子代處理的收穫、脫粒、乾燥和保存等程序與其親代一致。

本發明中，子代繁殖更新時單株與其親代是一一對應關係，即子代育苗時的各個處理每個單株的編號與其親代單株是一一對應的，進行形態性狀調查時，子代各個處理調查的單株與其親代單株也是一一對應關係，其目的是為了在進行大豆地方品種大黃豆遺傳完整性分析時明確子代與其親代的遺傳對應關係，使得出的試驗數據具有更高的可靠性和準確性，得出的試驗結論具有更強的說服力。親代與子代各個處理關係如圖1所示。

2、方法二：利用常規方法對大豆地方品種大黃豆進行繁殖更新。步驟如下：

(1)大黃豆親代處理繁殖更新

①繁殖地點選擇：為保證對比試驗的準確性和可靠性，繁殖地區、試驗地和配套條件與方法一相同。

②種子老化及生活力檢測：種子人工老化的材料、程序、方法與方法一相同，即同一批親代處理種子。生活力檢測的數據，包括發芽勢、發芽率、發芽指數和活力指數與方法一相同(表2)。

③播種：常規繁殖更新方法中沒有育苗及移栽的步驟，而是將親代處理種子直播。根據遺傳完整性分析的試驗要求，選取包括對照在內的4個適合的發芽率梯度的大豆，即處理GP-CK2、GP-I2、GP-II2和GP-III2進行直播。試驗地塊大小與方法一相同。每個小區根據發芽率播種，保證定苗200株。

④田間管理：播種前造墒，精細整地，重施底肥，一般畝施農家肥1000～2000kg，磷酸二銨15kg。出苗前後進行查苗，補苗，及時補種，但由於經過人工老化後，其活力較低，補苗後仍然出苗較差，因此不能保證種植密度一致。其餘栽培措施與方法一相同。需要注意地是每個處理的澆水施肥量必須相同，以免產生系統誤差。與方法一最大不同就是對單株不進行編號。

⑤形態性狀調查及去雜：在地方品種大黃豆親代處理繁殖更新期間，從每個處理中隨機選出30個單株，對其進行11個形態性狀調查，這些性狀包括：花色、粒色、粒形、臍色、茸毛色、莢色、結莢習性、株高、單株粒數、單株莢數和每莢粒數。

⑥收穫、脫粒、乾燥和保存：收穫方式為按照處理進行混合收穫，即以處理為單位進行收穫，以便繁殖其子代。脫粒方式為以處理為單位進行混合脫粒，紗網袋裝，標籤編號與處理編號一致，袋內外各附標籤。乾燥和清選方式與方法一相同。用紗網袋盛裝後保存於-4℃中期庫內。

(2)大黃豆子代處理繁殖更新

將親代各個處理繁殖更新所收穫的子代種子次年原地進行繁殖更新。子代處理的繁殖更新程序與其親代大體一致，但不需要經過種子老化。參考農作物種子檢驗規程GB/T3543.4中的標準條件，將各個子代處理，即處理GF1-CK2、GF1-I2、GF1-II2和GF1-III2進行發芽試驗，並統計發芽勢、發芽率、發芽指數和活力指數(表2)。

表2大豆地方品種大黃豆發芽數據統計(方法二)

需要指出地是利用方法一與方法二進行繁殖更新的不同生活力親代處理是同一批材料，因此生活力檢測數據是一致的。

(3)利用形態標記對經過兩種方法繁殖更新的地方品種大黃豆進行遺傳完整性分析

本發明所統計形態標記性狀如前所述共11個，包括：花色、粒色、粒形、臍色、茸毛色、莢色、結莢習性、株高、單株粒數、單株莢數和每莢粒數。按照方法一和方法二給出的程序進行形態性狀調查，方法一和方法二每個處理的樣本量均為30個單株；樣本量為30在統計學中算作大樣本，兩組採用同樣的樣本量，可以有效地消除系統誤差，使數據更加準確可靠。方法一是按照編號調查且親代處理與其子代處理調查的單株是一一對應關係，方法二則是隨機調查，親代處理與其子代處理沒有一一對應關係。

形態性狀多樣性的計算採用Shannon-Weaver信息指數，即H』＝-∑PilnPi，Pi為某性狀第i個代碼值出現的概率。質量性狀如花色、粒色和粒形等予以賦值(表3)。將數量性狀如株高、單株粒數、單株莢數和每莢粒數進行10級分類，1級<X-2δ，10級≥X+2δ，中間每級間差0.5δ，X為平均值，δ為標準差。11個形態性狀數據按不同的變異類型分別轉換成AA、BB和CC等字母格式，利用生物統計軟體Popgene等計算不同形態性狀變異的Shannon-Weaver指數。然後利用SAS V9.1將兩種方法繁殖更新的各個親代處理的Shannon-Weaver指數與各自對照進行顯著性t檢驗(表4)。

表3大豆地方品種大黃豆形態性狀鑑定項目及標準

表4經過兩種方法繁殖更新的地方品種大黃豆Shannon-Weaver指數t檢驗

註：*表示5％水平上顯著差異；**表示1％水平上顯著差異

結果顯示，方法一繁殖更新的親代處理中，發芽率為83.0％的處理GP-I1，其Shannon-Weaver指數與對照GP-CK1相比差異不顯著，發芽率61.0％和30.0％的處理GP-II1和GP-III1，其Shannon-Weaver指數與對照GP-CK1相比差異顯著或極顯著。這表明，較高生活力的處理(發芽率≥83.0％)與對照相比，其遺傳完整性得到了較好的保持，而較低生活力的處理(發芽率≤61.0％)與對照相比，其遺傳完整性發生了顯著變化。子代中各個處理(包括GF1-CK1、GF1-I1、GF1-II1和GF1-III1)的發芽率均為90.0％以上，其Shannon-Weaver指數與對照相比差異均不顯著，特別是處理GF1-CK1與對照GP-CK1相比幾乎沒有任何差異，表明高生活力的子代群體，與對照相比其遺傳完整性得到了較好的保持。這說明形態多樣性與生活力密切相關，原因是形態多樣性主要是數量性狀造成的，而質量形狀的差異很小，高生活力子代群體的數量性狀與對照群體的差異較小，而低生活力子代群體的數量性狀與對照群體的差異較大。

方法二繁殖更新的親代處理中，發芽率同為83.0％的處理GP-I2，其Shannon-Weaver指數與對照GP-CK2相比差異顯著，發芽率61.0％和30.0％的處理GP-II2和GP-III2，其Shannon-Weaver指數與對照GP-CK2相比差異極顯著。這表明由於繁殖方法不當較高生活力的親代處理(發芽率≥83.0％)與對照GP-CK2相比，其遺傳完整性就發生了顯著變化，而較低生活力的處理(發芽率≤61.0％)與對照GP-CK2相比，其遺傳完整性發生了極顯著變化。子代中處理GF2-CK2的發芽率為89.0％，與對照GP-CK2相比差異不顯著，但t值和概率值顯示，與方法一相比，經過方法二繁殖更新的子代處理，雖然其親代生活力很高，為93.0％，但是其子代仍然發生了形態多樣性變化；處理GF1-I2和GF1-II2的發芽率為88.0％和86.0％，與方法一相比，其發芽率有所下降，而Shannon-Weaver指數與對照GP-CK2相比差異雖不顯著，但其概率值已接近0.05，其差異水平快要達到顯著水平；處理GF1-III2的發芽率為85.0％，其Shannon-Weaver指數與對照GP-CK2相比差異顯著，表明其遺傳完整性就發生了顯著變化。

以上數據表明，由於繁殖方法不當，較高生活力的子代處理(發芽率≥85.0％)，其遺傳完整性也發生了顯著變化。推其原因：首先，親代處理中由於種子活力下降而導致田間成苗率下降，造成種植密度不一致，一些單株長得過於高大，其形態性狀，特別是數量性狀發生較大變化；其次，由於活力較差，田間成苗率大幅下降，親代傳遞給子代的遺傳多樣性降低，出現所謂的「奠基者效應」，導致子代的遺傳完整性發生顯著變化。

(4)利用SSR分子標記對經過兩種方法繁殖更新的地方品種大黃豆進行遺傳完整性分析

①分析樣本：

a、方法一：不同生活力親代處理及其子代處理在田間繁殖更新期間，選取96個單株(包括形態性狀調查的30個單株)，採集幼嫩葉片，-80℃下保存備用。子代採集葉片的單株與其親代也是一一對應關係。

b、方法二：不同生活力親代處理及其子代處理在田間繁殖更新期間，隨機選取96個單株(不一定包括形態性狀調查的30個單株)，採集幼嫩葉片，-80℃下保存備用。子代採集葉片的單株與其親代不存在一一對應關係。

②基因組DNA提取：採用SDS法單株提取上述兩種方法繁殖更新的親代及子代各個處理的基因組DNA。

③SSR引物擴增及聚丙烯醯氨凝膠檢測：採用40對SSR核心引物(選取自王棟等，SSR標記分析種子老化及繁殖世代對大豆種質遺傳完整性的影響，植物遺傳資源學報2010，11(2)：192-199。每個連鎖群上分布兩個位點，如表2所示)，對上述兩種方法繁殖更新的上述兩種方法繁殖更新的親代及其子代各個處理基因組DNA進行PCR擴增，採用6％聚丙烯醯氨凝膠電泳分離擴增產物，功率75W，時間為50min，電泳結束後採用常規銀染法檢測。

表5大豆地方品種大黃豆遺傳完整性分析40對SSR核心引物

a、PCR擴增反應體系為：Premix TaqTM(TaKaRa TaqTM Version2.0)10.0μL，5.0μmol/L Primer pairs(1.0+1.0)μL，20.0ng/μL DNA 5.0μL，ddH2O 3.0μL，總體系為20.0μL。

b、PCR擴增反應程序為：94℃預變性4min；95℃變性45s，47℃退火45s，72℃延伸45s，35個循環；72℃延伸10min。溫度降至4℃時取出，4℃冰箱內保存備用。

⑦結果統計及分析：根據檢測結果，利用POPGENE version 1.31和Powermarker V3.25等軟體對大豆地方品種大黃豆經過上述兩種方法繁殖更新的親代及其子代各個處理進行遺傳結構分析，計算出各個處理的各位點等位基因頻率、每位點等位基因數、遺傳多樣性指數和香農指數，並利用SAS 9.1對各處理群體與對照群體之間的各位點等位基因頻率進行χ2檢驗，對每位點等位基因數、遺傳多樣性指數和香農指數進行顯著性t檢驗。

a、等位基因頻率差異分析

表6大豆地方品種大黃豆處理間等位基因頻率差異分析(方法一)

從表6可以看出，方法一繁殖更新的親代處理GP-I1、GP-II1和GP-III1與對照GP-CK1的等位基因頻率顯著差異和極顯著差異的位點個數以及佔總位點個數百分比，隨著生活力的下降而增加，其中發芽率為30.0％的處理GP-III1的顯著和極顯著差異位點個數最多，分別為8個和4個，佔總位點的個數的百分比分別為20.00％和10.00％，發芽率為61.0％的處理GP-II1次之，分別為5和2個，佔總位點的個數的百分比分別為12.50％和5.00％，發芽率為83.0％的處理GP-I1最少，分別為1個和0個，佔總位點的個數的百分比分別為2.50％和0，這表明生活力下降顯著地影響了大豆種質材料處理內的等位基因頻率分布。由對照GP-CK1繁殖的處理GF1-CK1與對照GP-CK1相比，無顯著差異位點，由處理GP-1I繁殖的處理GF1-I1與對照GP-CK1相比，顯著或極顯著差異的位點個數與處理GP-I1持平，同為為1個和0個，佔總位點的個數的百分比分別為2.50％和0，這表明發芽率為93.0％和83.0％的大豆經過方法一繁殖更新後，各位點的等位基因頻率與對照相比幾乎沒有差異。由處理GP-II1和GP-III1繁殖的處理GF1-II1和GF1-III1，這2個處理與對照GP-CK1顯著或極顯著差異的位點個數較多，分別為6個和3個，9個和5個，佔總位點的個數的百分比分別為15.00％和7.50％，22.50％和12.50％，且子代處理比相應親代處理的顯著或極顯著差異位點數要高，表明發芽率為61.0％和30.0％的大豆種質材料子代處理各位點的等位基因頻率與對照相比差異顯著，且生活力水平愈低差異愈大。方法一結果表明，生活力下降對地方品種大黃豆的等位基因頻率分布影響甚大。

表7大豆地方品種大黃豆處理間等位基因頻率差異分析(方法二)

從表7可以看出，與方法一相比，方法二繁殖更新的各個處理與對照GP-CK2的等位基因頻率顯著差異和極顯著差異的位點個數以及佔總位點個數百分比普遍增加。親代處理GP-I2、GP-II2和GP-III2與對照GP-CK2的顯著和極顯著差異位點個數，同樣也隨著生活力的下降而增加，與方法一不同的是，發芽率為83.0％的處理GP-I2的顯著和極顯著差異位點個數分別為4個和2個，佔總位點的個數的百分比分別為10.00％和5.00％，而方法一相應的顯著和極顯著差異位點個數為1個和0個，佔總位點的個數的百分比分別為2.50％和0，表明經過方法二繁殖更新的較高生活力群體內的等位基因頻率分布發生了顯著變化。發芽率為61.0％和30.0％的處理GP-II2和GP-III2與對照GP-CK2相比，其顯著和極顯著差異位點個數分別為7個和5個，12個和8個，佔總位點的個數的百分比分別為17.50％和12.50％，30.00％和20.00％，都明顯高於方法一相對應的數據。方法二繁殖更新的處理GF1-CK2與對照相比，其顯著和極顯著差異位點個數分別為4個和3個，佔總位點的個數的百分比分別為10.00％和7.50％，而方法一相應的顯著和極顯著差異位點個數分別為0個，佔總位點的個數的百分比也是0，這表明發芽率為93.0％的處理經過方法二繁殖更新後，其等位基因頻率分布與對照相比也產生了顯著差異。由處理GP-I2、GP-II2和GP-III2繁殖的處理GF1-I2、GF1-II2和GF1-III2，這3個處理與對照GP-CK1顯著或極顯著差異的位點個數較多，分別為6個和4個，9個和7個以及14個和10個，佔總位點的個數的百分比分別為15.00％和10.00％，22.50％和17.50％以及35.00％和25.00％，與方法一的相應數據相比，都明顯高出甚多。以上數據說明方法二繁殖更新的各個處理(包括高生活力處理)的等位基因頻率分布與對照之間的差異明顯高於方法一。

b、群體遺傳結構分析

表8大豆地方品種大黃豆群體遺傳結構(方法一)

註：*表示5％水平上顯著差異；**表示1％水平上顯著差異

從表8中可以看出，親代處理GP-I1、GP-II1、GP-III1及其子代處理GF1-I1、GF1-II1和GF1-III1，其各項遺傳多樣性參數均低於對照GP-CK1，而且隨老化時間的延長，生活力水平愈低，其各項遺傳多樣性參數也愈低。t檢驗結果顯示，由對照GP-CK1繁殖的處理GF1-CK1，其各項遺傳多樣性參數與對照GP-CK1相比幾乎沒有差異，表明發芽率為93.0％的處理經過方法一繁殖更新後，其群體遺傳結構得到了較好的保持。處理GP-I1的A、H和I等3項遺傳多樣性參數與對照相比無顯著差異，而其子代處理GF1-I1的各項遺傳多樣性參數與對照相比雖有所下降但沒有達到顯著差異水平。處理GP-II1和GP-III1的A、H和I等3項遺傳多樣性參數與對照GP-CK1相比差異顯著或極顯著，其子代處理GF1-II和GF1-III的A、H和I等3個遺傳多樣性參數與對照相比差異極顯著，表明更新發芽率分別為61.0％和30.0％的處理因生活力水平下降，其自身及子代群體內的遺傳多樣性低於對照處理遺傳多樣性。方法一結果表明，生活力下降對大豆地方品種大黃豆的群體遺傳結構影響顯著。

表9大豆地方品種大黃豆群體遺傳結構(方法二)

註：*表示5％水平上顯著差異；**表示1％水平上顯著差異

從表9中可以看出，與方法一相比，方法二繁殖更新的各個處理包括對照的各項遺傳多樣性參數均有所下降。親代處理GP-I2、GP-II2、GP-III2及子代處理GF1-CK2、GF1-I2、GF1-II2和GF1-III2，其各項遺傳多樣性參數均低於對照GP-CK2，而且生活力水平愈低，其各項遺傳多樣性參數也愈低。t檢驗結果顯示，與方法一相比，由對照GP-CK2繁殖的處理GF1-CK2，其各項遺傳多樣性參數與對照GP-CK2相比差異雖未達到顯著水平，但下降的較多，表明發芽率為93.0％的處理經過方法二繁殖更新後，其群體遺傳結構有所改變，但未達到顯著水平。處理GP-I2及其子代處理GF1-I2的A、H和I等3項遺傳多樣性參數與對照GP-CK2相比差異顯著，這點與方法一不同，表明發芽率為83.0％的處理經過方法二繁殖更新後，其群體遺傳結構發生了顯著變化。處理GP-II2和GP-III2及其子代處理GF1-II2和GF1-III2的A、H和I等3項遺傳多樣性參數與對照GP-CK2相比差異極顯著，表明更新發芽率分別為61.0％和30.0％的處理因生活力水平下降，其自身及子代群體內的遺傳多樣性低於對照的遺傳多樣性，這方面兩種方法是一致的。以上結果表明，方法二繁殖更新的各個處理(包括高生活力處理)的群體遺傳結構與對照之間的差異明顯高於方法一，其群體內的遺傳多樣性降低較多。

綜合大豆地方品種大黃豆形態標記分析和SSR標記分析的實驗結果可得出以下結論：方法一繁殖更新的大豆種質材料，當更新發芽率高於83.0％時，親代及其子代在形態性狀水平及DNA分子水平上，其遺傳完整性都得到了較好的保持；更新發芽率低於61.0％時，親代及其子代在DNA分子水平上的遺傳完整性明顯改變。以此推薦其更新發芽率最低不得低於61.0％。方法二繁殖更新的大豆種質材料，即使更新發芽率為83.0％時，親代及其子代在形態性狀水平及DNA分子水平上，其遺傳完整性都發生了顯著變化，且得出的群體內遺傳多樣性低於方法一。因此，方法一較方法二更有利於大豆遺傳完整性的保持，建議採用本發明提供方法對大豆種質遺傳完整性分析時進行繁殖更新。