作為全球重要的油料和蛋白來源，大豆是一種典型的光周期敏感性短日照作物。在高緯度地區，例如中國的東北，夏季日照時間長，這給大豆帶來了一個不小的挑戰——開花延遲甚至無法正常成熟。開花時間是決定大豆適應性、成熟期和最終產量的關鍵生態性狀。長期以來，傳統育種已經成功培育出一些對光周期不敏感的大豆品種，使得大豆種植從東亞擴展到了歐洲和北美。然而，每個特定的大豆品種仍然只適應于一個有限的緯度范圍。要培育能夠適應不同光周期條件的高產大豆品種，特別是針對中國東北這樣具有長日照環境的地區，必須更深入地闡明開花時間的遺傳構成。

以往的研究已經鑒定出多個調控大豆開花和成熟的主效基因和位點，例如E1-E10、Tof11和Tof12系列基因。盡管這些基因在大豆光周期敏感性調控中扮演著重要角色，但對于大豆如何完全適應長期長日照條件的完整遺傳和分子機制，我們仍知之甚少。傳統的數量性狀位點(QTL)定位方法由于分辨率較低，限制了其在育種中的應用。相比之下，全基因組關聯分析(GWAS)因其更高的精度和分辨率，成為基因定位更有效的方法。然而，針對中國東北地區大豆種質開花時間遺傳基礎的研究仍然有限，這阻礙了該地區大豆育種和種植范圍的進一步擴展。

為了填補這一知識空白，一項發表在《The Crop Journal》上的研究，系統性地探索了中國東北大豆種質中開花時間（DTF）的遺傳基礎。這項研究的目標是鑒定穩定的基因組位點和基因，為培育適應多樣化光周期的高產大豆品種提供分子工具。

為開展此項研究，研究人員運用了幾項關鍵技術方法。研究材料來源于中國東北不同地區的200份大豆種質資源。首先，在多個環境點進行了多年份的表型數據采集和分析，包括聯合環境(CE)數據的計算和廣義遺傳力估計。其次，利用Illumina HiSeq平臺對全部材料進行了全基因組重測序，獲得了超過271萬個高質量單核苷酸多態性(SNP)標記用于后續分析。核心的分析策略是整合多種方法：利用GAPIT3軟件包中的六個不同GWAS模型（包括MLM、CMLM、MLMM、BLINK、SUPER和FarmCPU）進行全基因組關聯分析，以鑒定與DTF顯著相關的SNP。基于連鎖不平衡(LD)衰減距離，將穩定的顯著SNP周圍±202.4 kb的基因組區間界定為QTL。最后，結合基因本體(GO)富集分析、表達譜數據、實時定量PCR(qRT-PCR)以及單倍型分析，從QTL區間內篩選和驗證調控開花時間的候選基因。

研究對200份大豆種質在四個獨立環境下的開花時間進行了評估。結果表明，基因型(G)、環境(E)以及基因型與環境互作(G × E)對開花時間均有極顯著影響。開花時間的廣義遺傳力估計為0.43，表明該性狀具有中等的遺傳力，且在不同環境下存在顯著的變異，這反映了開花時間遺傳的復雜性。

利用271萬個高質量SNP進行的群體結構分析表明，該關聯分析群體呈連續分布，沒有明顯的亞群結構。連鎖不平衡(LD)衰減分析顯示，LD衰減到一半（r²= 0.405）時的遺傳距離為202.4 kb，此距離被用于界定候選QTL的物理區間。

通過六個GWAS模型在四個個體環境及聯合環境中的分析，共檢測到15個與DTF顯著相關的SNP，分布在10條大豆染色體上。其中，位于Chr.08、Chr.12、Chr.15和Chr.19上的5個SNP在多個環境和模型中均被檢測到，被視為穩定的標記-性狀關聯。

以上述5個穩定SNP為中心，劃定±202.4 kb的基因組區間，共鑒定出5個穩定的QTL，分別命名為qDTF8、qDTF12、qDTF15.1、qDTF15.2和qDTF19。其中，qDTF8、qDTF12、qDTF15.1和qDTF15.2是本研究首次發現的新的QTL，而qDTF19的區間與之前報道的“First flower 22–3”QTL重疊，但本研究將其物理區間大大縮小。

在5個穩定QTL的基因組區間內共存在84個基因模型。通過文獻調研、生物信息學分析、基因注釋和RNA-seq數據分析，初步篩選出20個可能的候選基因。這些基因的功能涉及分生組織從營養狀態向生殖狀態的轉變、與花起始和發育相關的激素信號（如生長素、細胞分裂素、乙烯和油菜素甾醇）、細胞分裂、花發育、蛋白質泛素化、磷酸化等過程。

進一步的qRT-PCR分析顯示，有6個基因（Glyma.08?g278400、Glyma.19?g200200、Glyma.19?g199300、Glyma.19?g200700、Glyma.19?g240700和Glyma.19?g197600）在早花和晚花品種的開花前(V4)和開花后(R2)階段表達存在顯著差異。同時，對20個候選基因進行的單倍型分析發現，其中10個基因的不同單倍型等位基因與DTF的表型變異顯著相關。綜合qRT-PCR和單倍型分析的結果，最終確定了14個基因是調控大豆DTF的潛在候選基因。其中，Glyma.19?g200700和Glyma.19?g197600在兩個分析中均得到驗證，被視為最有潛力的關鍵基因。

深入研究單倍型與基因表達的關系發現，僅Glyma.19?g197600基因的不同單倍型等位基因與其表達水平及DTF表型之間存在顯著關聯。例如，單倍型Hap_3與晚開花和較高的基因表達水平相關，而Hap_1則與早開花和較低的基因表達水平相關。這表明Glyma.19?g197600是早開花的負調控因子。此外，該基因的單倍型在中國不同緯度地區的大豆種質中呈現差異分布，暗示其在適應不同光周期環境中的重要作用。

研究還對已知的開花時間基因（如E1、E2、E3等）進行了序列變異和單倍型分析。結果顯示，在當前的東北大豆種質中，E1、E2、E3、E8、E10和Dt2等基因的單倍型對開花時間有顯著影響。然而，與之前研究不同的是，Tof12基因的單倍型在當前群體中差異不顯著，這可能是因為適應東北地區的早開花單倍型已在群體中趨于固定。這些結果共同表明，除了已知基因外，一套新的基因集合在調控適應中國東北環境的大豆開花時間中發揮著關鍵作用。

在討論部分，作者強調了開花時間對于大豆適應性和產量的極端重要性。本研究采用的整合策略——結合GWAS、QTL定位、單倍型和候選基因分析——成功揭示了東北大豆開花時間的復雜遺傳架構。與以往使用低通量標記和雙親群體的研究相比，本研究的GWAS方法顯著提高了QTL定位的分辨率，將qDTF19的區間大大縮小，并發現了4個全新的QTL。

研究鑒定出的候選基因功能豐富，涉及多個與開花調控密切相關的生物學過程和通路，如激素信號傳導、蛋白質翻譯后修飾等。特別是Glyma.19?g197600被證實為早開花的負調控因子，這為了解大豆在長日照條件下延遲開花的分子機制提供了新線索。該基因單倍型在中國不同生態區的差異分布，進一步印證了其在生態適應中的關鍵角色。

綜上所述，本研究系統解析了中國東北大豆開花時間的遺傳基礎，取得了多項重要結論：首先，鑒定出15個與開花時間顯著相關的SNP和5個穩定QTL，其中4個為首次報道。其次，通過綜合生物信息學、表達和遺傳分析，篩選出14個潛在的候選基因，并重點驗證了Glyma.19?g200700和Glyma.19?g197600的功能潛力，首次提出Glyma.19?g197600是早開花的負調控因子。最后，研究揭示了一套新的基因集合與已知開花基因共同調控著東北大豆的開花時間，凸顯了該性狀遺傳調控的復雜性。

這項研究的意義在于，它不僅為理解大豆在長日照環境下開花的遺傳機制提供了新的見解，更重要的是為分子設計育種提供了寶貴的基因資源和精確的分子靶點。所發現的位點和基因可用于開發分子標記，輔助選擇具有理想開花時間的大豆品種，從而培育出能夠適應短生長季、長光周期環境的高產、穩產品種，對于保障糧食安全、擴大大豆種植區域具有重要的理論和應用價值。