A. 實(shí)驗(yàn)設(shè)計(jì)示意圖。
B. 堆疊條形圖顯示60日齡和180日齡豬睪丸組織中三種不同甲基化類型(mCG、mCHG和mCHH;H=A、C或T)的平均比例。甲基化類型分別用綠色、橙色和藍(lán)色表示。
C. 小提琴圖描繪了三種不同序列環(huán)境(CG、CHG和CHH;H=A、C或T)中的甲基化密度。每一行表示一個(gè)不同的環(huán)境。Y軸表示甲基化密度,X軸代表不同的樣本(60日齡和180日齡)。
D. 小提琴圖展示三種不同序列環(huán)境(CG、CHG和CHH;H = A、C或T)中甲基化水平分布。
E. 60日齡與180日齡睪丸組織中不同基因組元件的甲基化水平。每一行表示一個(gè)不同的環(huán)境。X軸表示不同的基因組元件,Y軸表示甲基化水平。左側(cè)和右側(cè)面板分別顯示60日齡和180日齡睪丸組織的甲基化水平分布。綠色、紫色和黃色線條分別代表CG、CHG和CHH甲基化水平。
F. 60日齡與180日齡睪丸組織中上下游±2K區(qū)域的甲基化水平分布。
(2)DMR鑒定
圖2:豬睪丸組織中差異甲基化區(qū)域(DMRs)的鑒定。
A. CG、CHG和CHHDMR相關(guān)基因維恩圖。
B. CG環(huán)境中DMRs長(zhǎng)度分布。X軸表示DMRs長(zhǎng)度;Y軸表示每個(gè)長(zhǎng)度的密度,擬合曲線的DMRs長(zhǎng)度分布用紅色標(biāo)出。
C. CG環(huán)境中甲基化水平分布小提琴圖。X軸代表不同年齡組(60日齡和180日齡),Y軸代表平均甲基化水平的值。
D. CG環(huán)境中DMR錨定區(qū)域。X軸和Y軸分別表示不同基因區(qū)域和DMRs數(shù)量。
(3)DMG功能富集分析
圖3:差異甲基化基因(DMGs)的GO和KEGG分析。
A. 條形圖顯示CG中GO富集分析結(jié)果。Y軸描述不同的GO分析,X軸顯示DMGs數(shù)量。BP、CC和MF分別代表生物過(guò)程、細(xì)胞組分和分子功能。
B. KEGG通路富集散點(diǎn)圖示例。Y軸和X軸分別顯示通路名稱和基因比例。點(diǎn)大小表示DMGs數(shù)量,顏色對(duì)應(yīng)于不同的q值。
(4)DNA甲基化與差異表達(dá)基因(DEG)的關(guān)聯(lián)分析
圖4:DNA甲基化與基因表達(dá)之間的相關(guān)性分析。
A. DEGs和DMGs(CG、CHG和CHH環(huán)境)之間的重疊分析維恩圖。
B. CG環(huán)境中的DMGs與DEGs維恩圖,區(qū)分了與高甲基化(hyper)和低甲基化(hypo)相關(guān)基因,以及表達(dá)水平升高(up)和降低(down)的基因。
(5)與DEGs負(fù)相關(guān)DMG關(guān)鍵通路的富集
圖5:與DEGs負(fù)相關(guān)的DMGs的功能富集分析。
A. 條形圖顯示了與DEGs負(fù)相關(guān)的DMGs的GO富集。
B. 散點(diǎn)圖展示了與DEGs負(fù)相關(guān)的DMGs的KEGG通路。
表1:從與DEGs負(fù)相關(guān)的DMGs中篩選睪丸發(fā)育候選基因。
(6)鑒定參與睪丸發(fā)育的基因
圖6:DMGs及其表達(dá)模式的驗(yàn)證。
A. 基因的RT-qPCR表達(dá)水平與RNA-Seq數(shù)據(jù)一致;高甲基化抑制基因表達(dá),而低甲基化刺激基因表達(dá)。
B. BSP檢測(cè)到的DNA甲基化與WGBS數(shù)據(jù)一致。進(jìn)行三次生物學(xué)重復(fù)實(shí)驗(yàn),β-actin為對(duì)照基因。*p<0.05;**p<0.01;***p<0.001。
參考文獻(xiàn):
Feng Y, Zhang Y, Wu J, Qiao M, Zhou J, Xu Z, Li Z, Sun H, Peng X, Mei S. Comprehensive Analysis of Methylome and Transcriptome to Identify Potential Genes Regulating Porcine Testis Development. Int J Mol Sci. 2024 Aug 22;25(16) pii: ijms25169105. doi: 10.3390/ijms25169105. PubMed PMID: 39201790.