基于轉錄組和WGCNA篩選絲瓜果長和果徑發育調控相關基因

陳敏氡，王彬，劉建汀，李永平，白昌輝，葉新如，裘波音，溫慶放，朱海生

基于轉錄組和WGCNA篩選絲瓜果長和果徑發育調控相關基因

陳敏氡，王彬，劉建汀，李永平，白昌輝，葉新如，裘波音，溫慶放，朱海生

福建省蔬菜遺傳育種重點實驗室/福建省蔬菜工程技術研究中心/福建省農業科學院作物研究所，福州 350013

【目的】鑒定絲瓜果長和果徑發育基因共表達模塊并篩選關鍵調控基因，為后續絲瓜果形控制的分子機制研究提供理論依據?！痉椒ā恳越z瓜9個發育階段（開花前2 d以及花后0、2、4、6、8、10、15和20 d）果實為研究材料，測定各階段的果長和果徑，利用WGCNA方法聯合分析轉錄組與果長和果徑數據，鑒定與果長和果徑發育相關的共表達基因模塊，篩選關鍵調控基因?！窘Y果】利用WGCNA鑒定出14個共表達模塊，與果長和果徑顯著相關（相關系數的絕對值=0.9）的共表達模塊有2個，顯著正相關的模塊為Turquoise模塊，顯著負相關的模塊為Lightpink4模塊。KEGG富集分析發現，Turquoise模塊顯著富集在內吞作用和苯丙烷生物合成通路，與果實膨大、伸長調控密切相關，可作為研究絲瓜果長和果徑的關鍵基因模塊。根據Turquoise模塊內基因的連接度以及功能注釋，篩選出10個關鍵調控基因，包括木葡聚糖內轉葡糖基酶/水解酶基因、肌動蛋白解聚因子基因、伴侶蛋白基因、擴展蛋白基因（、和）、驅動蛋白基因、生長素反應基因和。RT-qPCR結果顯示，10個調控基因的表達量均在果實進入快速生長期（花后8 d）后顯著升高，增加倍數約為2—50倍。通過構建基因互作網絡，發現部分調控基因與WRKY、bHLH和HSF轉錄因子家族存在互作關系?！窘Y論】獲得了絲瓜果長和果徑基因共表達重要模塊Turquoise模塊，篩選到10個調控基因可作為絲瓜果形控制的潛在候選基因，并發現絲瓜果長和果徑的發育調控主要涉及細胞壁的重構、細胞的發育與分化、生長素的調控等過程。

轉錄組；WGCNA；絲瓜；果長；果徑；調控基因

0 引言

【研究意義】絲瓜（L.）為葫蘆科絲瓜屬作物，是我國夏秋供應市場的主要瓜類蔬菜之一。絲瓜果肉清香軟滑、營養豐富且具有較高的藥用功能，可用于治療筋骨酸痛、月經不調或產后乳汁不通等癥狀，因此深受消費者喜愛[1-2]。不同地區對絲瓜的消費習慣并不一致，南方地區比較傾向于棒形的絲瓜，而北方地區以圓筒形絲瓜和皺皮絲瓜為主。隨著現代家庭結構的改變，單瓜重約500 g的短棒、厚肉型絲瓜也逐漸受到消費者的追捧[3]?！Ｑ?號’[4]、‘雙豐1號’[5]、‘早優3號’[6]等短肉型絲瓜新品種相繼問世，現已在我國福建、浙江、廣東等沿海地區大面積種植，具有廣闊的市場前景和巨大的經濟潛力。由此可見，果長和果徑已成為絲瓜育種的重要考量指標，解析絲瓜果形控制的分子調控機制顯得尤為必要?！厩叭搜芯窟M展】果長和果徑是絲瓜重要的商品性狀，反映絲瓜果實的發育動態，與品質、產量性狀緊密相關。前人對于絲瓜果長、果徑的遺傳研究已有少量文獻報道，高軍等[7]認為，普通絲瓜的果長遺傳由一對加性主基因控制，同時存在多基因的修飾作用。林明寶等[8]研究發現，有棱絲瓜中控制果長性狀的最少基因數目為4對。蘇小俊等[9]研究顯示，普通絲瓜果長、果徑的主基因遺傳力均較大，宜在分離早代采取表型選擇法進行選擇。崔竣杰等[10-11]在絲瓜連鎖群scaf1上定位到1個貢獻率為28.90%的果長QTL位點，表明絲瓜果長受環境影響較大，認為需要在高世代中進行選擇。近年來，隨著絲瓜全基因組的公布，少數果實發育基因被研究。曾美娟等[12]發現6個生長素酰胺合成酶家族基因（gretchen hagen 3，）在普通絲瓜長、短果品種中差異表達，推測這些基因與絲瓜果長發育相關。隨后又發現SUN-domain蛋白家族基因在絲瓜花后12 d短果中的表達量極顯著高于長果，可能是育種潛力基因[13]。然而，目前已獲得的這些基因依然是對前人研究結果的驗證，調控絲瓜果長和果徑的主效基因及功能尚未可知。權重基因共表達網絡分析（WGCNA）方法能夠有效地整合轉錄組、代謝組和蛋白質組等多組學數據，正在成為發掘新基因功能和探索基因網絡與特定性狀間關聯關系的重要方法[14]。近幾年，在臍橙[15]、黃瓜[16]、西瓜[17]等果蔬形態發生研究中均可見利用WGCNA方法的報道。馮貴芝[15]發現介導柑橘果實ABA信號調節過程，在果實膨大期間起主要作用。李情[16]在黃瓜果實形態發生過程共鑒定出235個細胞分裂相關、1 020個細胞伸長相關的差異表達轉錄本，并發現、在其中起著至關重要的作用。孫蕾[17]發現胞壁酶和激素相關的基因，如聚半乳糖醛酸酶基因、果膠酯酶基因、生長素反應蛋白基因、脫落酸受體基因、乙烯反應轉錄因子、糖基轉移酶基因和擴展蛋白基因等參與調控了西瓜果實的膨大、成熟軟化過程。【本研究切入點】WGCNA為果蔬形態調控基因的挖掘提供了新途徑，目前尚未見采用WGCNA方法篩選絲瓜果長和果徑發育調控基因的報道。【擬解決的關鍵問題】本研究以普通絲瓜9個發育階段的果實為研究材料，分析果實在各個發育階段果長和果徑的變化動態，通過WGCNA方法聯合分析轉錄組與果長和果徑數據劃分基因共表達模塊，挖掘與果長和果徑相關聯的關鍵模塊及調控基因，并構建這些基因的互作調控網絡，為后續絲瓜果形控制的遺傳轉化研究提供潛在候選基因。

1 材料與方法

1.1 材料

試驗于2022年在福建省農業科學院作物研究所福清東張基地和福州科研綜合研究中心進行，選擇絲瓜雜交種‘福研3號’為試驗材料。

1.2 方法

1.2.1 果實果長、果徑測定 采集開花前2 d（記為花后-2 d）以及花后0、2、4、6、8、10、15和20 d長勢均勻的果實各3個，共計27個樣品，用游標卡尺測量其果長和果徑。

1.2.2 絲瓜果肉RNA的提取、純化、質檢和文庫構建及其數據分析 選取花后-2 d至15 d共24個樣品進行轉錄組測序。切下樣品中部約5 cm的果肉，液氮中速凍30 min，80℃保存，利用Plant RNA Purification Reagent（Invitrogen）提取總RNA，瓊脂糖凝膠電泳分析RNA的完整性及是否存在DNA污染，采用NanoDropOne（Thermo Fisher Scientific，USA）對RNA質量進行初步檢測（RIN≥6.5，OD260/280=1.8—2.2，OD260/230≥2.0）。按照Illumina TruSeqTMRNA sample preparation Kit方法構建轉錄組文庫，庫檢合格后采用Illumina MiSeq PE300（Illumina，USA）進行測序，由北京奧維森基因科技有限公司完成。高通量測序儀得到的圖像數據經CASAVA堿基識別轉化為序列數據（raw read）。對原始數據進行去除帶接頭、poly-N及低質量read后獲得高質量clean read。并對Q20、Q30、GC含量及重復序列水平進行統計；將clean read用HISAT2軟件與絲瓜參考基因組（https://peanutbase. org/data/public/Arachis_hypogaea/Tifrunner.gnm1.KYV3/）進行基因比對。

1.2.3 加權基因共表達網絡分析 利用RESM軟件對基因表達數據進行背景校正和標準化，過濾表達量過低和變異系數較小的基因[18]。利用R軟件WGCNA包構建基因共表達網絡分析（https://horvath.genetics. ucla.edu/html/CoexpressionNetwork/Rpackages/WGCNA/index.html）[19]。為了確保符合無尺度網絡分布，利用WGCNA包中的pickSoftThreshold函數計算合適的加權系數，的選取標準即滿足相關系數的平方接近0.8，同時還需要保證一定的基因連接度[20]。采用動態切割法識別共表達模塊[21]。利用自動網絡構建函數blockwiseModules構建網絡，其中最小模塊大小為50，合并相似性為0.75的模塊[20]。

1.2.4 特異性共表達模塊的鑒定、KEGG分析和基因表達趨勢分析 定義相關系數絕對值在0.5以上且顯著性值達到0.05水平的模塊為顯著模塊[22]。利用R語言的clusterProfiler包對顯著模塊內基因進行功能富集（KEGG）分析，列出顯著性前20個代謝通路做氣泡圖[23]。利用STEM軟件對特異性共表達模塊基因進行趨勢分析，設定模塊數量為20，為0.05，以獲得顯著性和非顯著性基因表達模型[24]。

1.2.5 特異性模塊核心基因的篩選及互作網絡的構建 利用softConnectivity函數計算基因的連接度，篩選出模塊內連接度前10的基因，結合已報道的與葫蘆科果長和果徑相關的基因（基因的連通性在模塊排名的前10%），定義這些基因為模塊的調控基因[25]。計算調控基因與模塊內全部基因的皮爾遜相關性系數，分別挑選與調控基因正負相關性前10的基因繪制網絡圖，利用Cytoscape v3.6.1進行可視化展示[23]。

1.2.6 RT-qPCR分析 采用通用植物總RNA提取試劑盒（北京百泰克生物技術有限公司）提取絲瓜花后-2 d至15 d共8個發育階段果實總RNA，利用HiScript III 1st Strand cDNA Synthesis Kit（+gDNA wiper）反轉錄試劑盒（南京諾唯贊生物科技股份有限公司）合成cDNA。利用Primer Premier 5軟件設計定量引物，具體引物序列見表1，以絲瓜為內參基因[26]。按照SYBR Premix Ex TaqTM Ⅱ熒光定量試劑盒（寶生物工程大連有限公司）說明書在ABI7500熒光定量PCR儀（Applied Biosystems，USA）上進行分析。擴增程序為95 ℃預變性30 s；95 ℃變性5 s，60 ℃退火34 s，40個循環。每個處理3個重復，并使用2-ΔΔCT方法計算基因的相對表達量[27]，采用Excel軟件作圖。

表1 本研究所用實時熒光定量PCR引物

2 結果

2.1 絲瓜果實發育過程中果長和果徑的生長動態分析

由圖1、2可知，絲瓜果實發育過程可以大致分為3個時期。花后-2 d至2 d為果實發育前期，這期間絲瓜果長和果徑生長緩慢；花后2—15 d為果實發育中期，這期間果長和果徑明顯增長、增大，其中花后8—15 d為果實的快速生長期和成熟期，果實膨大和伸長極顯著（<0.001）；花后15—20 d果實進入發育后期，這期間絲瓜果長和果徑基本不發生變化。

2.2 加權基因共表達網絡的構建

對絲瓜8個發育階段（花后-2 d至15 d）的24個絲瓜果實樣品進行RNA-seq數據分析，對表達矩陣中變異度較低的基因進行過濾后，選取26 620個基因進行WGCNA分析。當最優軟閾值=6時，2＞0.8，平均連通性趨于0（附圖1），表明該軟閾值適合用來構建無尺度網絡。根據基因的表達量進行聚類分析，聚類程度高的基因將被分到一個模塊（附圖2）。本研究共獲得14個基因共表達模塊，各模塊中所包含的基因數量差異明顯，其中Turquoise模塊中的基因數量最多（6 128個），Honeydew1模塊中的基因數量最少（86個）（圖3）。

2.3 特異性共表達模塊的鑒定

將14個基因共表達模塊分別與果長和果徑進行相關性分析，結果發現與絲瓜果長和果徑發育顯著相關的共表達基因模塊基本一致。其中，Turquoise模塊（=0.9，=3×10-9和=0.9，=1×10-9）和Brown4模塊（=0.64，=7×10-4和=0.65，=5×10-4）與果長和果徑顯著正相關，Lightpink4模塊（=-0.9，=2×10-9和=-0.9，=2×10-9）和Ivory模塊（=-0.67，=4×10-4和=-0.65，=6×10-4）與果長和果徑顯著負相關（圖4）。

圖3 各模塊的基因數量

2.4 基因共表達模塊的KEGG分析及表達趨勢分析

選取相關性最高的Lightpink4和Turquoise基因共表達模塊進行KEGG分析。結果顯示，Lightpink4模塊的基因主要顯著富集在氨基酸合成通路、剪接體通路、C5-二元酸代謝通路和丁酸代謝通路（圖5-A）。在生物體中，氨基酸合成通路、C5-二元酸代謝通路和丁酸代謝通路是產生能量、脂肪酸、氨基酸的重要途徑。剪接體負責pre-mRNA剪接，是基因表達與調控的重要環節之一[28]。2016年，哈佛公衛學院Heintz等[29]研究團隊首次將剪接體與生物體的壽命相關聯，并提出在未來剪接體可能作為衰老的一種生物標志物或是早期標記。本研究中Lightpink4模塊的基因與絲瓜果長和果徑顯著負相關（=-0.9，=2×10-9和=-0.9，=2×10-9），表明這些基因可能更多地參與到絲瓜果實發育后期至衰老過程。Turquoise模塊的基因主要顯著富集在內吞作用通路、氨基糖和核苷酸糖代謝通路、植物-病原相互作用通路和苯丙烷生物合成通路（圖5-B）。其中內吞作用和苯丙烷生物合成通路內顯著富集的基因數量最大（分別為52和41個）。內吞作用是植物細胞吸收水分、礦物質等的常見方式之一。德國弗萊堡大學DüNSER等[30]課題組發現內吞運輸能夠同步細胞和液泡表面的增大，推測這可能是植物細胞在細胞擴張速度方面表現出色的原因。苯丙烷生物合成的重要代謝物木質素能夠為植物提供機械支撐，同時參與導管的形成以及水分和礦質元素的運輸[31]。可見，這兩個通路均與果實膨大、伸長過程緊密相關。本研究中Turquoise模塊的基因與絲瓜果長和果徑顯著正相關（=0.9，=3×10-9和=0.9，=1×10-9），表明該模塊可作為研究絲瓜果長和果徑調控相關基因的重要模塊。因此，本研究重點關注Turquoise模塊。對Turquoise模塊中6 128個基因進行表達趨勢分析（圖6），Turquoise模塊中的基因可以歸納為20種表達趨勢模型，其中Profile 0、1、16、19為4個顯著基因表達模型（<0.05），Profile 0為顯著下調表達模型，包含349個基因；Profile 19為顯著上調表達模型，包含2 607個基因。

橫坐標代表不同表型，縱坐標代表不同模塊。圖中每組數據表示模塊與表型的相關性系數R值及顯著性P值（括號內）。紅色代表模塊與表型具有正相關，藍色代表模塊與表型具有負相關

圖5 Lightpink4（A）和Turquoise（B）基因共表達模塊的KEGG富集分析

2.5 絲瓜果長和果徑發育調控相關基因篩選

選取Turquoise模塊內連接度最高的前10個基因作為該模塊的核心基因，并通過與葫蘆科同源基因的比對和注釋來預測核心基因的功能。結果表明，Turquoise模塊內連接度最高的前10個基因中有6個基因有注釋信息，分別編碼木葡聚糖內轉葡糖基酶/水解酶蛋白（xyloglucan endotransglucosylase/ hydrolase protein 23，XTH23）（Maker00007620）、肌動蛋白解聚因子（actin-depolymerizing factor 2，ADF2）（Maker00033581）、伴侶蛋白DnaJ10（chaperone protein dnaJ 10，DnaJ10）（Maker00000726）、磷脂酰肌醇（phosphatidylinositol，PI）（Maker00012467）、同型半胱氨酸S-甲基轉移酶（homocysteine S- methyltransferase 1，HMT1）（Maker00010036）和鈣依賴性蛋白激酶（calcium-dependent protein kinase 1，CDPK1）（Maker00009822）。PI、HMTI和CPK1在植物體中主要參與激活信號傳遞[32]、甲基傳遞[33]和鈣信號轉導[34]等過程。而XTH23、ADF2和DnaJ10在植物細胞壁重組[35]、細胞生長[36]和葉綠體的形成及分化[37]等生命過程起重要作用，與果長和果徑發育密切相關，因此，認為這3個基因可作為絲瓜果長和果徑發育候選基因。繼續在Turquoise模塊中篩選葫蘆科中已報道與果長和果徑發育相關的基因，結果在該模塊中共發現7個已報道過的基因，且基因的連通性均在模塊排名的前10%，包括4個擴展蛋白基因（Maker00000768、Maker00039481、Maker00026031、Maker00009740）、1個驅動蛋白基因（Maker00024997）以及2個生長素反應蛋白基因（Maker00001372）和（Maker00038241）（表2）。因此，本研究將以上10個基因作為絲瓜果長和果徑發育的候選基因。

方框內數字為P值 The number in the box denoting the P-value

表2 絲瓜果長和果徑發育調控的候選基因

2.6 調控基因的表達模式分析

進一步利用RT-qPCR方法對10個調控基因的表達情況進行驗證，其結果與轉錄組數據基本一致（圖7），說明轉錄組數據結果的可信度較高。RT-qPCR結果顯示，10個候選基因大致包括3種表達模式，Maker00000726、Maker00038241、Maker00001372和Maker00024997這4個基因隨著果實的生長發育，表達量逐漸升高，其中花后10—15 d后升高明顯，表達量增加約2—6倍。Maker00007260和Maker00033581的表達量在發育前期（花后-2 d至2 d）基本不變，進入發育中期后共有2次表達量的顯著升高，花后4—8 d表達量增加約2—6倍，花后10—15 d表達量增加顯著，約10—50倍。Maker00026031、Maker00000768、Maker00039481和Maker00009740的表達量在花后8 d前基本不變，花后10—15 d表達量顯著升高，增加約10—40倍。

2.7 調控基因互作網絡的構建

挑選與10個調控基因正負相關性前10的基因繪制網絡圖（附表1），利用Cytoscape v3.6.1進行可視化展示（圖8）。結果顯示與（Maker00007620）負相關性最高的基因編碼UV-B誘導蛋白（Maker00013740）（相關系數=0.980），正相關性最高的基因編碼富亮氨酸重復受體蛋白激酶（Maker00026072）（=0.985）；與2（Maker00033581）負相關性最高的基因編碼谷氨酰-tRNA還原酶1（Maker00033993）（=0.930），正相關性最高的基因編碼DA1相關蛋白2（Maker00009584）（=0.970）；與（Maker00000726）負相關性最高的基因編碼3-脫氫奎尼酸脫水酶/莽草酸脫氫酶（Maker00027020）（=0.943），正相關性最高的基因編碼真核生物翻譯起始因子3亞基J（Maker00008454）（=0.986）；與（Maker00000768）負相關性最高的基因編碼羥?；入赘孰乃饷福∕aker00028749）（= 0.926），正相關性最高的基因編碼受體激酶At5g67200（Maker00012477）（=0.965）；與（Maker00039481）負相關性最高的基因編碼葡萄糖-6-磷酸1-脫氫酶（Maker00001558）（=0.921），正相關性最高的基因編碼生長素外排載體（Maker00029480）（=0.939）；與（Maker00026031）負相關性最高的基因編碼受體蛋白激酶At5g15080（Maker00025844）（=0.934），正相關性最高的基因編碼生長素誘導蛋白22d（Maker00039055）（=0.988）；與（Maker00024997）負相關性最高的基因編碼轉運抑制反應蛋白1（Maker00016411）（=0.907），正相關性最高的基因編碼TIFY 9蛋白（Maker000）（=0.924）；與（Maker00009740）負相關性最高的基因編碼E3泛素蛋白連接酶RING1（Maker00022044）（=0.893），正相關性最高的基因編碼轉錄因子bHLH149（Maker00003629）（= 0.912）；與（Maker00001372）負相關性最高的基因編碼Csa_5G517870蛋白（Maker00029937）（=0.817），正相關性最高的基因編碼LOC101210747亞型X1蛋白（Maker00004186）（=0.858）；與（Maker00038241）負相關性最高的基因編碼非典型聚（A）聚合酶PAPD5（Maker00016990）（=0.968），正相關性最高的基因編碼胚芽水解酶7亞型X1（Maker00006161）（=0.977）。此外，在網絡圖中注釋到3類轉錄因子，分別編碼（Maker00039492）、9（Maker00036518）、（Maker0002332）、（Maker00003629）和（Maker00028475）。、、、分別與、、、存在互作關系。前人已在擬南芥[38]、苜蓿[39]、楊樹[40]中證實WRKY和bHLH類轉錄因子參與木質素生物合成的調控。木質素是植物細胞壁的主要成分，、、、均與植物細胞壁合成相關。由此推測這2類轉錄因子可能也參與到絲瓜果實生長膨大過程中細胞壁的合成、松弛和伸展。調控植物生長發育，在西瓜果實發育[41]和番茄植株矮化[42]過程中均發揮著重要作用。與存在互作關系，推測二者共同參與絲瓜果實細胞的生長、分化。

圖7 10個調控基因的RT-qPCR和RNA-seq結果

圖中連線顏色代表調控基因與關聯基因的相關性類型，紅色代表正相關，藍色代表負相關；線條粗細代表相關性數值大小。調控基因及其關聯的轉錄因子用黃色標注

3 討論

2008年，Langfelder和Horvath[19]提出了加權共表達網絡分析（WGCNA）方法[19]，國內外學者通過該方法尋找到了許多與植物器官形態建成和發育調控的關鍵基因[43-45]。因此，本研究將WGCNA方法引入絲瓜果長和果徑發育調控基因的篩選。

3.1 絲瓜果長和果徑基因共表達重要模塊的鑒定

在共表達網絡分析中共獲得14個基因共表達模塊，并發現與絲瓜果長和果徑發育顯著相關的模塊基本一致，推測絲瓜果長和果徑發育可能由許多相同的基因調控。通過共表達模塊與表型的相關性分析以及模塊內基因的KEGG分析，關注到Turquoise模塊。該模塊中有41個共表達基因顯著富集在苯丙烷生物合成途徑。苯丙烷類生物合成途徑是植物3條主要次生代謝途徑之一，該途徑中的一個重要產物木質素與植物細胞的分化及器官生長過程密切相關。晏巢等[46]研究顯示油茶果皮中木質素積累過程影響了果實發育，木質素積累越早，果實越小，并且種子的形態也與木質素具有直接的關系。朱海英等[47]研究表明，絲瓜果實的發育過程伴隨著大量木質素的累積。課題組前期也發現5種木質素合成蛋白（PAL、4-CL、CAD、POD和CCoAOMT）的表達水平在絲瓜果實發育過程顯著升高，說明絲瓜果實發育過程中苯丙烷生物合成活躍，木質素大量合成[48]。由此表明，Turquoise模塊可作為研究絲瓜果長和果徑基因的重要模塊。未來可以結合細胞生物學和代謝組學手段研究絲瓜果實內木質素的結構類型、積累規律以及活躍區域等，進一步探究木質素與絲瓜果長和果徑發育調控的關系。

3.2 絲瓜果形控制潛在的候選基因及功能預測

共表達模塊內具有高連通性的基因被視為核心基因，這些基因通常起重要調控作用[23,49]。本研究對Turquoise模塊內連通度前10的基因進行功能注釋，篩選出3個與果長和果徑發育相關的調控基因，分別編碼木葡聚糖內轉葡糖基酶/水解酶XTH23、肌動蛋白解聚因子ADF2和伴侶蛋白DnaJ10。XTH是植物細胞壁重構過程中的關鍵酶，不僅能夠松弛、降解細胞壁，而且能參與細胞壁的強化和合成[50]。有研究表明，XTH在調控黃瓜[51]、番茄[52]、柿子[53]、甜櫻桃[54]等果實伸長膨大、成熟及軟化過程中發揮著重要作用。本研究中，（Maker00007620）是Turquoise模塊內連通性最高的基因。RT-qPCR顯示絲瓜果實進入快速生長期（花后8 d）后，的表達量顯著升高，花后10—15 d基因表達量較發育前期（花后-2 d至2 d）增加約33—48倍?？梢姡谡{控絲瓜果長和果徑發育過程中起著重要調控作用。ADF是植物細胞的結構、生長、發育和分化中必不可少的因素。水稻的和擬南芥的與花粉的形成和花粉管的伸長有關[55-56]。楊梅的ADF蛋白對果實的質地變化起重要作用[57]。Qiao等[58]發現玉米的異源過表達提高了轉基因擬南芥種子的大小，推測該基因可能在玉米籽粒發育中發揮重要作用。本研究中，（Maker00033581）同樣在絲瓜果實的快速生長期大量表達，花后10—15 d基因表達量較發育前期增加約7—11倍，提示該基因可能也與絲瓜果實的伸長及膨大相關。DnaJ是多種作物葉綠體的形成和分化過程中必不可少的重要因子[36]。果實葉綠體的豐度和功能直接影響著果實發育和品質。本研究中，（Maker00000726）的表達量在絲瓜果實的發育中期顯著升高，推測該基因可能在該過程中調控果實葉綠素形成，而葉綠素作為同化器官合成大量同化物用于絲瓜果實發育。

此外，本研究還在Turquoise模塊中篩選到7個已在葫蘆科中報道過與果長和果徑發育相關的基因，它們分別編碼4個擴展蛋白Expansin、1個驅動蛋白Kinesin以及2個生長素反應蛋白SAUR和Aux/IAA。Expansin與XTH的功能相似，也在植物細胞壁的松弛與重構中發揮重要作用[35]。郝西等[59]鑒定了黃瓜擴展蛋白基因家族，發現該家族基因在黃瓜果實發育過程中表達豐度較高。劉昌[60]在甜瓜中克隆到，發現敲除該基因的T1代果實與野生型相比出現果實變小的表型，表明能促進甜瓜果實成長期的膨大生長。本研究篩選到4個（Maker00026031、Maker00000768、Maker00039481和Maker00009740），發現這4個基因的表達量在果實進入快速生長期（花后8 d）前基本不變，隨后表達量均顯著升高，表明它們在果實果長、果徑的發育過程發揮重要作用。其中Maker00039481和Maker00009740編碼的和在果實成熟期（花后10—15 d）出現第二次的高表達，推測它們在果實成熟過程也起著重要作用。研究發現植物具有明顯的時空效應，有些的表達與果實發育、膨大相關，而有些是果實成熟軟化的特異基因[61]。Kinesin參與細胞分裂、生長和形態發生等植物生長發育過程。有研究報道，在黃瓜[62]、西瓜[63]等果實發育早期的細胞分裂和細胞膨大過程中起重要作用。本研究中，（Maker00024997）的表達量隨著絲瓜果實發育逐漸升高，并且在果實快速生長期顯著升高，可見很可能參與絲瓜果實整個生長發育期的細胞快速分裂過程。生長素是植物最重要的激素之一。/和作為兩類主要的生長素早期誘導基因在葫蘆科果實發育調控中發揮關鍵作用[64]。袁靜賢[65]從西瓜中鑒定到4個/與果實發育早期細胞分裂呈正相關，7個/與果實膨大呈正相關。于越等[66]前期通過EMS誘變方法獲得黃瓜短果突變體，對其進行轉錄組分析發現一個關鍵的上調表達差異基因，該基因屬于Aux/IAA家族。LUO等[67]在冬瓜果實中發現一個特異性高表達的，發現在擬南芥中過表達該基因會導致角果變長，由此推測該基因在冬瓜果長調控中起重要作用。本研究中，（Maker00038241）和（Maker00001372）的表達量均隨著果實發育逐漸升高，并在果實快速生長期顯著高表達，可見生長素及其響應基因對絲瓜果實發育過程具有十分重要的調節作用。通過對以上10個調控基因的功能分析，發現絲瓜果長和果徑的發育調控主要涉及細胞壁的重構、細胞的發育與分化、生長素的調控等過程，同時可能也與葉綠素的形成和分化相關。與此同時，研究還發現絕大多數的調控基因與多種作物果實的成熟軟化存在著密不可分的關系，因此，推測這些基因可能也參與絲瓜果實的成熟軟化過程，但其具體調控方式還需要進一步深入研究。

4 結論

本研究基于轉錄組和WGCNA構建了絲瓜果長和果徑發育相關基因共表達網絡，共鑒定出14個共表達模塊，獲得了關鍵基因調控模塊Turquoise模塊，篩選到、、、、、、、和等10個調控基因可作為絲瓜果形控制潛在的候選基因。

[1] 李佳欣, 馮玉. 絲瓜不同部位藥理作用研究進展. 食品工業科技, 2021, 42(10): 355-361.

LI J X, FENG Y. Research progress on pharmacological action of different parts of(L.) Roem. Science and Technology of Food Industry, 2021, 42(10): 355-361. (in Chinese)

[2] OBOH I O, ALUYOR E O.-an emerging cash crop. African Journal of Agricultural Research, 2009, 4(8): 684-688.

[3] 王青青, 王天文, 高安輝. 絲瓜種質資源與育種研究進展. 現代園藝, 2019(21): 33-35.

WANG Q Q, WANG T W, GAO A H. Research progress on germplasm resources and breeding of luffa. Xiandai Horticulture, 2019(21): 33-35. (in Chinese)

[4] 李大忠, 朱海生, 康玉妹, 溫慶放, 李永平, 康建坂. 普通絲瓜新品種‘福研1號’選育. 亞熱帶植物科學, 2016, 45(1): 87-89.

LI D Z, ZHU H S, KANG Y M, WEN Q F, LI Y P, KANG J B. Breeding of a new sponge gourd cultivar ‘Fuyan 1’. Subtropical Plant Science, 2016, 45(1): 87-89. (in Chinese)

[5] 林奕韓, 陳捷凱, 林愛惜. 絲瓜新品種雙豐1號肉絲瓜的選育. 廣東農業科學, 2010, 37(7): 46-47.

LIN Y H, CHEN J K, LIN A X. Breeding of a new luffa variety Shuangfeng No.1 meat luffa. Guangdong Agricultural Sciences, 2010, 37(7): 46-47. (in Chinese)

[6] 王安樂, 朱海泉, 鄧穩橋, 張望隆. 絲瓜新品種早優3號的選育. 中國蔬菜, 2007(11): 27-28.

WANG A L, ZHU H Q, DENG W Q, ZHANG W L. A new sponge gourd F1hybrid - ‘Zaoyou No. 3’. China Vegetables, 2007(11): 27-28. (in Chinese)

[7] 高軍, 徐海, 蘇小俊, 袁希漢. 普通絲瓜果長遺傳規律分析. 江蘇農業科學, 2007, 35(5): 123-125.

GAO J, XU H, SU X J, YUAN X H. Analysis on genetic rule of fruit length in luffa. Jiangsu Agricultural Sciences, 2007, 35(5): 123-125. (in Chinese)

[8] 林明寶, 林師森. 有棱絲瓜果長遺傳效應的初步研究. 華南農業大學學報, 2000, 21(2): 8-9, 24.

LIN M B, LIN S S. Prelimary study on the heredity of fruit length of angular sponge gourd. Journal of South China Agricultural University, 2000, 21(2): 8-9, 24. (in Chinese)

[9] 蘇小俊, 徐海, 高軍, 鄭子松, 宋波, 陳龍正, 袁希漢, 陳勁楓. 普通絲瓜果實性狀的遺傳分析. 江蘇農業學報, 2009, 25(5)1112-1118.

SU X J, XU H, GAO J, ZHENG Z S, SONG B, CHEN L Z, YUAN X H, CHEN J F. Genetic analysis on fruit characters of luffa. Jiangsu Journal of Agricultural Sciences, 2009, 25(5)1112-1118. (in Chinese)

[10] 崔竣杰, 程蛟文, 譚澍, 李衛鵬, 汪國平, 胡開林. 絲瓜果長及果柄長的遺傳規律分析. 廣東農業科學, 2014, 41(8): 52-56.

CUI J J, CHENG J W, TAN S, LI W P, WANG G P, HU K L. Genetic analysis of fruit length and fruit stalk length in luffa. Guangdong Agricultural Sciences, 2014, 41(8): 52-56. (in Chinese)

[11] 崔竣杰, 呂振, 楊天文, 王靜, 洪宇, 曹毅. 絲瓜高密度bin標記遺傳圖譜構建與果長QTL定位. 浙江農業學報, 2023, 35(3): 590-597.

CUI J J, Lü Z, YANG T W, WANG J, HONG Y, CAO Y. Construction of high-density bin marker genetic map and QTL mapping for fruit length inspp. Acta Agriculturae Zhejiangensis, 2023, 35(3): 590-597. (in Chinese)

[12] 曾美娟, 劉建汀, 李祖亮, 陳敏氡, 葉新如, 王彬, 朱海生, 溫慶放. 普通絲瓜基因家族全基因組鑒定及表達分析. 江蘇農業科學, 2022, 50(18): 82-87.

ZENG M J, LIU J T, LI Z L, CHEN M D, YE X R, WANG B, ZHU H S, WEN Q F. Genome-wide identification and expression analysis ofgene family in. Jiangsu Agricultural Sciences, 2022, 50(18): 82-87. (in Chinese)

[13] 曾美娟, 陳敏氡, 劉建汀, 葉新如, 王彬, 林錦輝, 朱海生, 溫慶放. 普通絲瓜基因的克隆及表達分析. 中國細胞生物學學報, 2022, 44(6): 1031-1038.

ZENG M J, CHEN M D, LIU J T, YE X R, WANG B, LIN J H, ZHU H S, WEN Q F. Cloning and expression analysis ofgene from. Chinese Journal of Cell Biology, 2022, 44(6): 1031-1038. (in Chinese)

[14] ZHANG B, HORVATH S. A general framework for weighted gene co-expression network analysis. Statistical Applications in Genetics and Molecular Biology, 2005, 4: 17.

[15] 馮貴芝. 基于時空轉錄組的臍橙果實發育和成熟調控機理的研究[D]. 武漢: 華中農業大學, 2021.

FENG G Z. Study on the regulation mechanism of navel orange fruit development and ripening based on spatio-temporal transcriptome [D]. Wuhan: Huazhong Agricultural University, 2021. (in Chinese)

[16] 李情. 黃瓜果實早期發育中基因的時空表達與調控[D]. 北京: 中國農業科學院, 2019.

LI Q. Temporal and spatial expression and regulation of genes in early development of cucumber fruit [D]. Beijing: Chinese Academy of Agricultural Sciences, 2019. (in Chinese)

[17] 孫蕾. 西瓜中心果肉硬度主效QTL及候選基因分析[D]. 哈爾濱: 東北農業大學, 2021.

SUN L. Analysis of major QTL and candidate genes for firmness of watermelon central pulp [D]. Harbin: Northeast Agricultural University, 2021. (in Chinese)

[18] TRAPNELL C, ROBERTS A, GOFF L, PERTEA G, KIM D, KELLEY D R, PIMENTEL H, SALZBERG S L, RINN J L, PACHTER L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nature Protocols, 2012, 7(3): 562-578.

[19] LANGFELDER P, HORVATH S. WGCNA: An R package for weighted correlation network analysis. BMC Bioinformatics, 2008, 9: 559.

[20] 馬娟, 曹言勇, 王利鋒, 李晶晶, 王浩, 范艷萍, 李會勇. 利用WGCNA鑒定玉米株高和穗位高基因共表達模塊. 作物學報, 2020, 46(3): 385-394.

MA J, CAO Y Y, WANG L F, LI J J, WANG H, FAN Y P, LI H Y. Identification of gene co-expression modules of maize plant height and ear height by WGCNA. Acta Agronomica Sinica, 2020, 46(3): 385-394. (in Chinese)

[21] 周雨青, 楊永飛, 葛常偉, 沈倩, 張思平, 劉紹東, 馬慧娟, 陳靜, 劉瑞華, 李士叢, 趙新華, 李存東, 龐朝友. 基于WGCNA的棉花子葉抗冷相關共表達模塊鑒定. 中國農業科技導報, 2022, 24(4): 52-62.

ZHOU Y Q, YANG Y F, GE C W, SHEN Q, ZHANG S P, LIU S D, MA H J, CHEN J, LIU R H, LI S C, ZHAO X H, LI C D, PANG C Y. Identification of cold-related co-expression modules in cotton cotyledon by WGCNA. Journal of Agricultural Science and Technology, 2022, 24(4): 52-62. (in Chinese)

[22] 王澤涵, 于文濤, 王鵬杰, 劉財國, 樊曉靜, 谷夢雅, 蔡春平, 王攀, 葉乃興. 茶樹禿房與茸房種質花器官差異表達基因的WGCNA分析. 園藝學報, 2023, 50(3): 620-634.

WANG Z H, YU W T, WANG P J, LIU C G, FAN X J, GU M Y, CAI C P, WANG P, YE N X. WGCNA analysis of differentially expressed genes in floral organs of tea germplasms with ovary-glabrous and ovary-trichome. Acta Horticulturae Sinica, 2023, 50(3): 620-634. (in Chinese)

[23] 郭永春, 王鵬杰, 金珊, 侯炳豪, 王淑燕, 趙峰, 葉乃興. 基于WGCNA鑒定茶樹響應草甘膦相關的基因共表達模塊. 中國農業科學, 2022, 55(1): 152-166. doi: 10.3864/j.issn.0578-1752.2022. 01.013.

GUO Y C, WANG P J, JIN S, HOU B H, WANG S Y, ZHAO F, YE N X. Identification of co-expression gene related to tea plant response to glyphosate based on WGCNA. Scientia Agricultura Sinica, 2022, 55(1): 152-166. doi: 10.3864/j.issn.0578-1752.2022.01.013. (in Chinese)

[24] ERNST J, BAR-JOSEPH Z. STEM: A tool for the analysis of short time series gene expression data. BMC Bioinformatics, 2006, 7(1): 191.

[25] 薦紅舉, 張梅花, 尚麗娜, 王季春, 胡柏耿, Vadim Khassanov, 呂典秋. 利用WGCNA篩選馬鈴薯塊莖發育候選基因. 作物學報, 2022, 48(7): 1658-1668.

JIAN H J, ZHANG M H, SHANG L N, WANG J C, HU B G, KHASSANOV V, Lü D Q. Screening candidate genes involved in potato tuber development using WGCNA. Acta Agronomica Sinica, 2022, 48(7): 1658-1668. (in Chinese)

[26] 朱海生, 陳敏氡, 溫慶放, 藍新隆, 李永平, 王彬, 張前榮, 吳衛東. 絲瓜基因克隆及其作為內參基因的應用. 核農學報, 2016, 30(1): 35-41.

ZHU H S, CHEN M D, WEN Q F, LAN X L, LI Y P, WANG B, ZHANG Q R, WU W D. Cloning ofgene fromand its application as an internal standard. Journal of Nuclear Agricultural Sciences, 2016, 30(1): 35-41. (in Chinese)

[27] LIVAK K J, SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2?ΔΔCTmethod. Methods, 2001, 25(4): 402-408.

[28] 施一公. 剪接體的結構與分子機理研究. 科學通報, 2020, 65(27): 2952-2953.

SHI Y G. Structural and mechanistic investigations of the spliceosome. Chinese Science Bulletin, 2020, 65(27): 2952-2953. (in Chinese)

[29] HEINTZ C, DOKTOR T K, LANJUIN A, ESCOUBAS C, ZHANG Y, WEIR H J, DUTTA S, SILVA-GARCíA C G, BRUUN G H, MORANTTE I, HOXHAJ G, MANNING B D, ANDRESEN B S, MAIR W B. Splicing factor 1 modulates dietary restriction and TORC1 pathway longevity in. Nature, 2017, 541(7635): 102-106.

[30] DüNSER K, SCH?LLER M, R??LING A K, L?FKE C, XIAO N N, PA?íZKOVá B, MELNIK S, RODRIGUEZ-FRANCO M, ST?GER E, NOVáK O, KLEINE-VEHN J. Endocytic trafficking promotes vacuolar enlargements for fast cell expansion rates in plants. eLife, 2022, 11: e75945.

[31] 于明革, 楊洪強, 翟衡. 植物木質素及其生理學功能. 山東農業大學學報(自然科學版), 2003, 34(1): 124-128.

YU M G, YANG H Q, ZHAI H. Lignin and physiological function in plant. Journal of Shandong Agricultural University (Natural Science Edition), 2003, 34(1): 124-128. (in Chinese)

[32] 肖玫, 趙旭琴, 許迎科, 李漢兵. 磷脂酰肌醇信號在GLUT4囊泡轉運中的調控作用. 中國細胞生物學學報, 2018, 40(1): 139-145.

XIAO M, ZHAO X Q, XU Y K, LI H B. The regulatory role of phosphatidylinositol signaling in GLUT4 vesicle trafficking. Chinese Journal of Cell Biology, 2018, 40(1): 139-145. (in Chinese)

[33] 張載勵, 丁晨, 王爽, 梁明朗, 李濤, 蘇來曼·買買提, 胡康棣. 番茄、擬南芥和葡萄中同型半胱氨酸S-甲基轉移酶基因家族功能研究. 現代農業科技, 2023(3): 66-71.

ZHANG Z L, DING C, WANG S, LIANG M L, LI T, SULAIMAN MAIAITI, HU K D. Study on function of homocysteine S-methyltransferase gene family in tomato,and grape. Modern Agricultural Science and Technology, 2023(3): 66-71. (in Chinese)

[34] 何樂平, 張蕾. 植物鈣依賴性蛋白激酶及其相關蛋白激酶(CDPKs/CRKs)的研究進展. 安徽農業科學, 2020, 48(18): 26-30, 51.

HE L P, ZHANG L. Research progress of plant calcium-dependent protein kinases(CDPKs) and CDPK-related kinases(CRKs). Journal of Anhui Agricultural Sciences, 2020, 48(18): 26-30, 51. (in Chinese)

[35] 賈鑫磊, 何貝軒, 郭丹丹, 郭美麗. 膨脹素和木葡聚糖內轉葡糖基酶/水解酶基因的功能研究進展. 植物生理學報, 2018, 54(11): 1659-1668.

JIA X L, HE B X, GUO D D, GUO M L. Research progress in the function of expansins and xyloglucan endotransglucosylase/hydrolase. China Industrial Economics, 2018, 54(11): 1659-1668. (in Chinese)

[36] 王楠, 呂香玲, 李亮, 吳金鳳, 尚佳微, 郝轉芳, 李鳳海. 植物肌動蛋白解聚因子ADF的研究進展. 西北植物學報, 2015, 35(11): 2349-2354.

WANG N, Lü X L, LI L, WU J F, SHANG J W, HAO Z F, LI F H. Research progress of ADF/cofilin in plant. Acta Botanica Boreali-Occidentalia Sinica, 2015, 35(11): 2349-2354. (in Chinese)

[37] 樊芳菲, 楊暹, 康云艷, 柴喜榮, 蒙林平. 植物DnaJ蛋白的研究進展. 分子植物育種, 2018, 16(6): 2028-2034.

FAN F F, YANG X, KANG Y Y, CHAI X R, MENG L P. Research progress on DnaJ proteins in plants. Molecular Plant Breeding, 2018, 16(6): 2028-2034. (in Chinese)

[38] YAN L, XU C H, KANG Y L, GU T W, WANG D X, ZHAO S Y, XIA G M. The heterologous expression inof sorghum transcription factor SbbHLH1 downregulates lignin synthesis. Journal of Experimental Botany, 2013, 64(10): 3021-3032.

[39] WANG H Z, AVCI U, NAKASHIMA J, HAHN M G, CHEN F, DIXON R A. Mutation of WRKY transcription factors initiates pith secondary wall formation and increases stem biomass in dicotyledonous plants. Proceedings of the National Academy of Sciences of the United States of America, 2010, 107(51): 22338-22343.

[40] 王依依. 楊樹、葡萄、擬南芥和水稻WRKY Ⅲ基因家族比較及表達分析[D]. 合肥: 安徽農業大學, 2016.

WANG Y Y. Comparison and expression analysis of WRKY Ⅲ gene family in poplar, grape,and rice [D]. Hefei: Anhui Agricultural University, 2016. (in Chinese)

[41] HAMEED SABA. 西瓜熱激轉錄因子家族基因的鑒定及其在果實發育過程中的表達分析[D]. 武漢: 華中農業大學, 2019.

HAMEED S. Identification of heat shock transcription factor family genes in watermelon and their expression analysis during fruit development [D]. Wuhan: Huazhong Agricultural University, 2019. (in Chinese)

[42] 王華容. 番茄熱激轉錄因子SlHSFC1過表達導致轉基因植株矮化的機理研究[D]. 重慶: 西南大學, 2023.

WANG H R. Study on the mechanism of dwarfing transgenic plants caused by overexpression of tomato heat shock transcription factor SlHSFC1 [D]. Chongqing: Southwest University, 2023. (in Chinese)

[43] WANG R C, SHU P, ZHANG C, ZHANG J L, CHEN Y, ZHANG Y X, DU K, XIE Y, LI M Z, MA T, ZHANG Y, LI Z G, GRIERSON D, PIRRELLO J, CHEN K S, BOUZAYEN M, ZHANG B, LIU M C. Integrative analyses of metabolome and genome-wide transcriptome reveal the regulatory network governing flavor formation in kiwifruit (). The New Phytologist, 2022, 233(1): 373-389.

[44] CHEN L, TIAN N, HU M Q, SANDHU D, JIN Q F, GU M Y, ZHANG X Q, PENG Y, ZHANG J L, CHEN Z Y, LIU G Z, HUANG M D, HUANG J N, LIU Z H, LIU S Q. Comparative transcriptome analysis reveals key pathways and genes involved in trichome development in tea plant (). Frontiers in Plant Science, 2022, 13: 997778.

[45] CHEN Z Y, LIAO M L, YANG Z J, CHEN W Y, WEI S H, ZOU J, PENG Z S. Co-expression network analysis of genes and networks associated with wheat pistillody. PeerJ, 2022, 10: e13902.

[46] 晏巢, 姚小華, 殷恒福, 王開良, 尹潤富, 滕建華. 油茶果實生長發育動態和木質素積累規律. 江西農業大學學報, 2020, 42(4): 788-801.

YAN C, YAO X H, YIN H F, WANG K L, YIN R F, TENG J H. Fruit development dynamics and lignin accumulation law of oil tea. Acta Agriculturae Universitatis Jiangxiensis (Natural Sciences Edition), 2020, 42(4): 788-801. (in Chinese)

[47] 朱海英, 李人圭, 王隆華, 黃祥輝, 顏季瓊. 絲瓜果實發育中木質素代謝及有關導管分化的生理生化研究. 華東師范大學學報(自然科學版), 1997(1): 87-94.

ZHU H Y, LI R G, WANG L H, HUANG X H, YAN J Q. Physiological and biochcmical studies on lignin mctabolism and its treachary element differentiation of fruit of. Journal of East China Normal University (Natural Science), 1997(1): 87-94. (in Chinese)

[48] 陳敏氡, 王彬, 朱海生, 李永平, 劉建汀, 溫慶放. 絲瓜果實發育過程的蛋白質差異表達分析. 中國細胞生物學學報, 2018, 40(11): 1806-1818.

CHEN M D, WANG B, ZHU H S, LI Y P, LIU J T, WEN Q F. Proteomic analysis of differentially expressed proteins induring fruit development. China Industrial Economics, 2018, 40(11): 1806-1818. (in Chinese)

[49] 彭佳偉, 張葉, 寇單單, 楊麗, 劉曉飛, 張學英, 陳海江, 田義. ‘倉方早生’桃及其早熟芽變不同發育時期果實的轉錄組分析. 中國農業科學, 2023, 56(5): 964-980. doi:10.3864/j.issn.0578-1752.2023.05. 012.

PENG J W, ZHANG Y, KOU D D, YANG L, LIU X F, ZHANG X Y, CHEN H J, TIAN Y. Transcriptome analysis of peach fruits at different developmental stages in peach kurakato wase and early-ripening mutant. Scientia Agricultura Sinica, 2023, 56(5): 964-980. doi:10.3864/j.issn.0578-1752.2023.05.012. (in Chinese)

[50] 杜麗萍, 沈昕, 陳少良, 胡贊民. 細胞壁重構關鍵酶木葡聚糖內轉糖苷酶/水解酶(XTH)的研究進展. 農業生物技術學報, 2010, 18(3): 604-609.

DU L P, SHEN X, CHEN S L, HU Z M. Research advances on a key cell wall remodeling enzyme xyloglucan endotransglucosylase/ hydrolase(XTH). Journal of Agricultural Biotechnology, 2010, 18(3): 604-609. (in Chinese)

[51] 傅豐慶. 黃瓜果實發育調控的分子機理研究[D]. 杭州: 浙江大學, 2008.

FU F Q. Study on molecular mechanism of cucumber fruit development regulation [D]. Hangzhou: Zhejiang University, 2008. (in Chinese)

[52] TSUCHIYA M, SATOH S, IWAI H. Distribution of XTH, expansin, and secondary-wall-related CesA in floral and fruit abscission zones during fruit development in tomato (). Frontiers in Plant Science, 2015, 6: 323.

[53] HAN Y, ZHU Q G, ZHANG Z K, MENG K, HOU Y L, BAN Q Y, SUO J T, RAO J P. Analysis of xyloglucan endotransglycosylase/ hydrolase () genes and diverse roles of isoenzymes during persimmon fruit development and postharvest softening. PLoS One, 2015, 10(4): e0123668.

[54] ZHAI Z F, FENG C, WANG Y Y, SUN Y T, PENG X, XIAO Y Q, ZHANG X, ZHOU X, JIAO J L, WANG W L, DU B Y, WANG C, LIU Y, LI T H. Genome-wide identification of the xyloglucan endotransglucosylase/hydrolase (XTH) and polygalacturonase (PG) genes and characterization of their role in fruit softening of sweet cherry. International Journal of Molecular Sciences, 2021, 22(22): 12331.

[55] LI L, LI Y X, SONG S F, DENG H F, LI N, FU X Q, CHEN G H, YUAN L P. An anther development F-box (ADF) protein regulated by tapetum degeneration retardation (TDR) controls rice anther development.Planta, 2015, 241(1): 157-166.

[56] 彭世清, 黃冬芬. 擬南芥肌動蛋白解聚因子4基因()在煙草中表達引起植株形態變化. 植物生理與分子生物學學報, 2006, 32(1): 52-56.

PENG S Q, HUANG D F. Expression of anactin- depolymerizing factor 4 gene () in tobacco causes morphological change of plants. Journal of Plant Physiology and Molecular Biology, 2006, 32(1): 52-56. (in Chinese)

[57] 陳義勇. 不同發育時期楊梅果肉差異表達蛋白及其功能研究[D]. 福州: 福建農林大學, 2014.

CHEN Y Y. Study on differentially expressed proteins inpulp at different developmental stages and their functions[D]. Fuzhou: Fujian Agriculture and Forestry University, 2014. (in Chinese)

[58] QIAO D H, DONG Y B, ZHANG L, ZHOU Q, HU C H, REN Y L, LI Y L. Ectopic expression of the maizeZmADFgene inrevealing its functions in kernel development. Plant Cell, Tissue and Organ Culture, 2016, 126(2): 239-253.

[59] 郝西, 理向陽, 臘貴曉, 代丹丹, 楊鐵鋼. 黃瓜擴展蛋白基因家族的鑒定與生物信息學分析. 分子植物育種, 2015, 13(10): 2280-2289.

HAO X, LI X Y, LA G X, DAI D D, YANG T G. Identification and bioinformatic analysis of the expansin gene family in cucumber. Molecular Plant Breeding, 2015, 13(10): 2280-2289. (in Chinese)

[60] 劉昌. 甜瓜擴展蛋白基因家族成員的鑒定及CmEXPA8和CmEXPA12基因在果實發育中功能的初步分析[D]. 呼和浩特: 內蒙古大學, 2018.

LIU C. Identification of members of melon expansin gene family and preliminary analysis of the functions of CmEXPA8 and CmEXPA12 genes in fruit development [D]. Hohhot: Inner Mongolia University, 2018. (in Chinese)

[61] 張安, 曹清河, 周志林, 趙冬蘭, 李昂, 李元元, 唐君. 植物細胞壁松弛因子: Expansin研究進展. 江蘇農業科學, 2013, 41(6): 11-13.

ZHANG A, CAO Q H, ZHOU Z L, ZHAO D L, LI A, LI Y Y, TANG J. Research progress of plant cell wall relaxation factor-Expansin. Jiangsu Agricultural Sciences, 2013, 41(6): 11-13. (in Chinese)

[62] YANG X Y, WANG Y, JIANG W J, LIU X L, ZHANG X M, YU H J, HUANG S W, LIU G Q. Characterization and expression profiling of cucumber kinesin genes during early fruit development: Revealing the roles of kinesins in exponential cell production and enlargement in cucumber fruit. Journal of Experimental Botany, 2013, 64(14): 4541-4557.

[63] TIAN S J, JIANG J, XU G Q, WANG T, LIU Q Y, CHEN X E, LIU M, YUAN L. Genome wide analysis of kinesin gene family inreveals an essential role in early fruit development. BMC Plant Biology, 2021, 21(1): 210.

[64] 李艷艷, 齊艷華. 植物Aux/IAA基因家族生物學功能研究進展. 植物學報, 2022, 57(1): 30-41.

LI Y Y, QI Y H. Advances in biological functions of Aux/IAA gene family in plants. Bulletin of Botany, 2022, 57(1): 30-41. (in Chinese)

[65] 袁靜賢. 生長素相關轉錄因子家族基因在西瓜果實膨大過程中的調控作用分析[D]. 武漢: 華中農業大學, 2015.

YUAN J X. Regulation of auxin-related transcription factor family genes in watermelon fruit expansion [D]. Wuhan: Huazhong Agricultural University, 2015. (in Chinese)

[66] 于越, 陳海燕, 周龍, 郭冬雪, 宮思宇, 郝寧, 杜亞琳, 武濤. 黃瓜果實發育相關基因CsIAA29的克隆與表達分析. 華北農學報, 2020, 35(S1): 11-17.

YU Y, CHEN H Y, ZHOU L, GUO D X, GONG S Y, HAO N, DU Y L, WU T. Cloning and expression analysis of CsIAA29 gene related to fruit development in cucumber. Acta Agriculturae Boreali-Sinica, 2020, 35(S1): 11-17. (in Chinese)

[67] LUO C, YAN J Q, HE C X, LIU W R, XIE D S, JIANG B. Genome-wide identification of the SAUR gene family in wax gourd () and functional characterization of BhSAUR60 during fruit development. International Journal of Molecular Sciences, 2022, 23(22): 14021.

Screening Regulatory Genes Related to Luffa Fruit Length and Diameter Development Based on Transcriptome and WGCNA

CHEN MinDong, WANG Bin, LIU JianTing, LI YongPing, BAI ChangHui, YE XinRu, QIU BoYin, WEN QingFang, ZHU HaiSheng

Fujian Key Laboratory of Vegetable Genetics and Breeding/Fujian Engineering Research Center for Vegetables/Crops Research Institute, Fujian Academy of Agricultural Sciences, Fuzhou 350013

【Objective】 The aim of this study was to identify the co expression modules of luffa fruit length and diameter development and to screen key regulatory genes, so as to provide the theoretical basis for subsequent research on the molecular mechanism of fruit shape control in luffa. 【Method】 The luffa fruits in 9 fruit development stages (2 days before anthesis, and 0, 2, 4, 6, 8, 10, 15, and 20 days after anthesis) were applied as research materials. The fruit length and diameter of each stage were measured. The WGCNA method was used to jointly analyze transcriptome and fruit length and diameter data, to identify co-expressed gene modules of fruit length and diameter development, and to screen out key regulatory genes.【Result】A total of 14 co expression modules were identified by WGCNA, among which two modules (Turquoise and Lightpink4) were significantly correlated with fruit length and diameter (absolute value of correlation coefficient=0.9); Turquoise module was significantly positively correlated, while Lightpink4 module was significantly negatively correlated. KEGG enrichment analysis found that the Turquoise module was significantly enriched in endocytosis and phenylpropanoid biosynthesis pathways, which were closely related to fruit enlargement and growth regulation, and could be used as a key gene module for studying fruit length and diameter in luffa. According to the connectivity and functional annotation of genes in Turquoise module, ten key regulatory genes were screened, including xyloglucan endotransglucosylase/hydrolase gene, actin-depolymerizing factor gene, chaperone protein gene, expansin gene (,and), kinesin gene, auxin response genes, and. The RT-qPCR results showed that the expression levels of ten regulatory genes significantly increased after the fruit entered the rapid growth period (8 day after anthesis), with an increase of 2-50 times approximately. Through constructing a gene interaction network, it was found that some candidate genes interacted with the,, andtranscription factor families.【Conclusion】The Turquoise module, an important co expression module of luffa fruit length and diameter was obtained, and ten potential candidate genes for luffa fruit shape control were screened. It was found that luffa fruit length and diameter development regulation mainly involved the processes of cell wall reconstruction, cell development and differentiation, and auxin regulation.

transcriptome; WGCNA;L.; fruit length; fruit diameter; regulatory genes

10.3864/j.issn.0578-1752.2023.22.012

2023-04-10；

2023-09-05

福建省自然科學基金（2021J01494）、福建省農業科學院英才項目（YC2021004）、福建省農業科學院科技創新團隊項目（CXTD2021003-1）、國家大宗蔬菜產業技術體系福州綜合試驗站項目（CARS-23-G51）、福建省種業創新與產業化工程項目（zycxny2021009）、福建省農業高質量發展超越“5511”協同創新工程項目（XTCXGC2021003）

陳敏氡，E-mail：1053870789@qq.com。通信作者朱海生，E-mail：zhs0246@163.com。通信作者溫慶放，E-mail：fjvrc@163.com

（責任編輯 趙伶俐）