水產(chǎn)養(yǎng)殖環(huán)境脅迫對(duì)魚類表觀遺傳的影響研究進(jìn)展

楊震飛， 劉波、2， 戈賢平、2， 宋長(zhǎng)友， 張慧敏， 單凡

(1.南京農(nóng)業(yè)大學(xué)無錫漁業(yè)學(xué)院，江蘇無錫214081；2.中國(guó)水產(chǎn)科學(xué)研究院淡水漁業(yè)研究中心，江蘇無錫214081)

在一定的環(huán)境條件下，由于外界環(huán)境因素的負(fù)面影響，使動(dòng)物無法維持正常的生理狀態(tài)，稱為環(huán)境脅迫 (environmental stress)，也叫環(huán)境應(yīng)激。生物總會(huì)受到生存環(huán)境中各種環(huán)境因子的刺激，且刺激與生物有機(jī)體總是處于相互作用的動(dòng)態(tài)過程中。近年來，隨著工廠化水產(chǎn)養(yǎng)殖模式迅速發(fā)展，其集約化程度越來越高。在集約化養(yǎng)殖中，養(yǎng)殖密度的提高會(huì)產(chǎn)生一系列不利于魚類生長(zhǎng)的環(huán)境脅迫，引起水質(zhì)惡化，影響魚類的成活率、生長(zhǎng)、健康、福利和產(chǎn)量[1]。

經(jīng)典遺傳學(xué)認(rèn)為，遺傳的分子基礎(chǔ)是核酸，核酸的堿基序列上包含生物體的所有遺傳信息，幾乎所有的生命活動(dòng)均由基因調(diào)控。然而，生物體是一個(gè)復(fù)雜的系統(tǒng)，不是一種因素就能決定整個(gè)生命活動(dòng)。隨著遺傳學(xué)的發(fā)展，人們發(fā)現(xiàn)基因表達(dá)模式的改變是由一些DNA或染色質(zhì)水平的修飾引起的。表觀遺傳 (epigenetics)是指DNA序列不發(fā)生變化，但基因表達(dá)卻發(fā)生了可遺傳的改變。這種改變是細(xì)胞內(nèi)除遺傳信息外其他可遺傳物質(zhì)發(fā)生的改變，且這種改變?cè)诎l(fā)育和細(xì)胞增殖過程中能穩(wěn)定傳遞[2]。它是一種全新的遺傳機(jī)制，因?yàn)樗簧婕癉NA序列的改變，不符合孟德爾的遺傳方式。在遺傳學(xué)研究中，人們發(fā)現(xiàn)了表觀遺傳有多種修飾方式，包括DNA甲基化作用、組蛋白修飾作用、非編碼RNA作用等，而且探究了其錯(cuò)綜復(fù)雜的生物學(xué)作用。研究發(fā)現(xiàn)，表觀遺傳與外界環(huán)境刺激有密切的關(guān)系[3]，由于外界環(huán)境對(duì)機(jī)體刺激的影響，任何一種調(diào)節(jié)機(jī)制發(fā)生紊亂都能導(dǎo)致細(xì)胞活動(dòng)發(fā)生異常，從而引起各種疾病[4]。已有研究表明，環(huán)境脅迫能調(diào)控水生生物基因的表達(dá)，影響水生生物的表型性狀，并且在一定程度上改變遺傳信息并傳遞到下一代，進(jìn)而影響下一代的表型性狀。本研究中，就水產(chǎn)養(yǎng)殖中環(huán)境脅迫與表觀遺傳調(diào)控的研究進(jìn)展做一綜述，以期為進(jìn)一步研究環(huán)境因素與基因互作關(guān)系提供理論參考。

1 表觀遺傳學(xué)的基本要素

表觀遺傳是指基因的表達(dá)發(fā)生可遺傳的變異現(xiàn)象，而基因序列不發(fā)生改變[5]，是生物體內(nèi)遺傳物質(zhì)和環(huán)境因素共同作用的結(jié)果，主要包括4個(gè)方面:DNA甲基化、組蛋白修飾、染色質(zhì)重塑和非編碼RNA調(diào)控等[6]。

1.1 DNA甲基化修飾

1.1.1 DNA甲基化 DNA甲基化是指在DNA甲基轉(zhuǎn)移酶 (DNA methyltransferase，DNMT)的催化作用下，通過S－腺苷－L－甲硫氨酸 (S－adenymethionine，SAH)提供甲基供體，將甲基轉(zhuǎn)移到脫氧胞嘧啶環(huán)第5位碳原子上形成甲基化脫氧胞嘧啶共價(jià)修飾[7]，DNA甲基化過程由DNA轉(zhuǎn)移酶家族(DNMTS)中的維持甲基化轉(zhuǎn)移酶 (DNMT1)和重新甲基化轉(zhuǎn)移酶 (DNMT3A和DNMT3B)參與完成，其中DNMT1主要負(fù)責(zé)DNA復(fù)制后的甲基化，以保證子鏈和親鏈具有相同的甲基化模式；DNMT3A和DNMT3B則主要負(fù)責(zé)對(duì)未甲基化兩條鏈的甲基化過程，以及發(fā)育過程中DNA的重新甲基化，也參與了異常甲基化過程[8－10]。

DNA甲基化主要發(fā)生在CPG島上，CPG島通常存在許多重要的基因啟動(dòng)子區(qū)。DNA的高度甲基化會(huì)影響基因的結(jié)構(gòu)，進(jìn)而影響RNA轉(zhuǎn)錄，引起基因沉默。正常的DNA甲基化是調(diào)控生物體基因表達(dá)的重要方式之一，參與較多生理過程，如維持染色體結(jié)構(gòu)、胚胎發(fā)育、細(xì)胞分化、X染色體失活、肌生成等。而異常的DNA甲基化則會(huì)引發(fā)疾病，近年來研究表明，遺傳過程中存在類似DNA甲基化作用的其他表觀遺傳修飾 (如羥基化)，這些作用可能與魚類的生長(zhǎng)發(fā)育存在密切的關(guān)系[11]。1.1.2 DNA去甲基化 DNA的去甲基化同DNA的甲基化一樣，在基因的表達(dá)過程中扮演著非常重要的作用。DNA去甲基化的兩個(gè)主要過程是主動(dòng)去甲基化 (Active DNA demethylation)和被動(dòng)去甲基化 (Passive DNA demethylation)。在 DNA去甲基酶參與過程下形成的去甲基化為主動(dòng)去甲基化，而在DNA的復(fù)制過程中，由于DNA的甲基化維持機(jī)制受損導(dǎo)致的去甲基化為被動(dòng)去甲基化。

DNA的去甲基化主要是由10－11轉(zhuǎn)置蛋白(Ten－Eleven Translocation，TET)介導(dǎo)的去甲基化過程[12]。 TET1可催化 5－MEC形成 5－HMC(5－h(huán)ydroxymethylcytosine)。 5－HMC 中間體在 5－MEC變成非甲基化C中起著重要作用[13]。沉默抑制因子 (Repressor of Silencing 1，ROSl)是擬南芥Arabidopsis thaliana體內(nèi)的一種DNA去甲基化酶，具有裂解酶和糖基化酶的活性，同時(shí)含有核酸內(nèi)切酶結(jié)合域，從而介導(dǎo)DNA去甲基化。ROS1突變會(huì)導(dǎo)致一些內(nèi)源位點(diǎn)和轉(zhuǎn)基因位點(diǎn)的高度甲基化，同時(shí)也會(huì)對(duì)DNA損傷試劑敏感[14]。ROS4具有乙酰轉(zhuǎn)移酶的活性，是催化組蛋白H3乙酰化的主要活性物質(zhì)，ROS4傾向于結(jié)合不含有修飾的核小體的特異甲基化DNA，組蛋白的乙?；癁槿旧w創(chuàng)造一個(gè)相對(duì)松散的環(huán)境，利于ROS1發(fā)揮作用，同時(shí)DNA甲基化減少蛋白 (Decrease in DNA methylation1，DDM1) 的功能也需要 ROS1的參與[15－16]。ROS5具有Α晶體狀結(jié)構(gòu)規(guī)域，編碼一個(gè)熱激蛋白，可以同ROS4和ROS1互化，同時(shí)也影響DNA去甲基化和組蛋白的修飾[17－18]。近些年相關(guān)研究發(fā)現(xiàn)，AP核酸內(nèi)切酶和堿基切除修復(fù)途徑X線修復(fù)交叉互補(bǔ)基因 (X－ray repair cross－complementing1 XRCC1)也參與 DNA去甲基化途徑[19－20]。同時(shí)也有研究表明，甲基化胞嘧啶結(jié)合蛋白(Methylcytosine－Binding Protein， MBD4) 可結(jié)合含有DNA甲基化的位點(diǎn)，指導(dǎo)DNA的去甲基化[21－22]。

1.2 組蛋白修飾

1.2.1 組蛋白修飾 組蛋白修飾是另一個(gè)主要的表觀遺傳過程 (圖1)。組蛋白修飾包括一系列的轉(zhuǎn)錄后修飾，如甲基化、乙酰化、生物?；?、磷酸化、泛素化等[23－24]。乙酰化組蛋白與活化基因相關(guān)，而去乙酰化通常與轉(zhuǎn)錄沉默相關(guān)[25]。組蛋白甲基化對(duì)基因活性有著不同的影響，組蛋白H3第4位賴氨酸 (K4)甲基化通常和增加基因活性相關(guān)，而組蛋白H3第9位賴氨酸 (K9)甲基化可能導(dǎo)致轉(zhuǎn)錄抑制。更重要的可能是DNA甲基化與組蛋白修飾通常共同起作用[26]，比如，胞嘧啶的甲基化可能有助于H3K9的甲基化。同時(shí)，H3K9甲基化可能促進(jìn)胞嘧啶的甲基化進(jìn)而導(dǎo)致基因下調(diào)。通常DNA甲基化和組蛋白甲基化修飾相互依賴的起作用進(jìn)而調(diào)控整體染色體的狀態(tài)。早期研究發(fā)現(xiàn)，高乙?；M蛋白和轉(zhuǎn)錄活化基因相關(guān)，組蛋白乙?；兄诨虻霓D(zhuǎn)錄[27－28]。實(shí)際上，組蛋白乙?；幕瘜W(xué)性質(zhì)表明了其可能促進(jìn)基因表達(dá)的一個(gè)機(jī)制。組蛋白甲基化被當(dāng)作調(diào)控元件，在許多真核生物啟動(dòng)子區(qū)域富集H3K4me3，但H3K4甲基轉(zhuǎn)移酶較少[29]。在體外，H3K4me3通過果蠅Drosophila melanogaster、人類多梳抑制性復(fù)合物 2(Polycomb RepressiveComplex2， PRC2) 抑制H3K27 三甲基化[30]。此外，H3K9me3 和H3K27me3也與異染色質(zhì)的形成以及多梳沉默相關(guān)。組蛋白磷酸化影響其修飾位點(diǎn)殘基的正電荷，表明組蛋白磷酸化與乙?；邢嗨频恼{(diào)節(jié)核小體動(dòng)力學(xué)機(jī)制[31]。

1.2.2 組蛋白去甲基化去乙?；{(diào)節(jié)機(jī)制 組蛋白甲基化也是動(dòng)態(tài)的，賴氨酸特異性去甲基化酶1(Lysine Specific Demethylase1，LSD1)是一種組蛋白去甲基化酶，也是一種胺氧化酶，可在輔基存在的條件下進(jìn)行胺氧化反應(yīng)催化去除H3K4結(jié)構(gòu)位點(diǎn)的二甲基化和單甲基化，進(jìn)一步研究發(fā)現(xiàn)，LSD1在雄激素受體引導(dǎo)下也可去除H3K4的二甲基化[32]。另一類組蛋白去甲基化酶是含有JmjC結(jié)構(gòu)域的蛋白 (JmjC Domain Containing Protein，JMJ)類蛋白，也是最大的一類去甲基化酶，可去除多甲基化修飾包括無法去除的三甲基化蛋白，并通過活性中心結(jié)構(gòu)域?qū)谆鶊F(tuán)進(jìn)行氧化反應(yīng)，再脫去水釋放甲酸而將甲基去除[33]。類蛋白通過去組蛋白甲基化以及與其他染色質(zhì)修飾相互協(xié)作等方式參與了眾多的基因表達(dá)和染色質(zhì)活動(dòng)的調(diào)控，涉及到了代謝、發(fā)育和環(huán)境適應(yīng)等生物學(xué)過程。

組蛋白的乙?；腿ヒ阴；l(fā)生在不同組蛋白的不同位點(diǎn)。DNA與組蛋白間、核小體與核小體間相互作用的強(qiáng)度可通過組蛋白賴氨酸殘基加上乙?；揎椇笾泻妥陨淼恼姾蓙砀淖僛34]。每種生物均含有多種組蛋白乙酰轉(zhuǎn)移酶和去乙酰化酶，組蛋白乙酰轉(zhuǎn)移酶和去乙?；缸R(shí)別的位點(diǎn)一般不止一個(gè)，這與組蛋白甲基轉(zhuǎn)移酶有較大不同?；虻霓D(zhuǎn)錄激活同樣離不開組蛋白的去乙?；揎棥OS2是Ⅰ型組蛋白去乙?；?，關(guān)于HOS2對(duì)基因的轉(zhuǎn)錄調(diào)控機(jī)制的觀點(diǎn)認(rèn)為:一般表達(dá)量高的基因均需要組蛋白的多種修飾以及多個(gè)轉(zhuǎn)錄激活因子的激活作用，但這些修飾可能成為后期轉(zhuǎn)錄的障礙[35]。因此，對(duì)于需要多輪轉(zhuǎn)錄的基因來說，組蛋白去乙?；窰OS2負(fù)責(zé)恢復(fù)轉(zhuǎn)錄后染色質(zhì)的結(jié)構(gòu)，以方便下一輪轉(zhuǎn)錄。

1.3 染色質(zhì)重塑

染色質(zhì)重塑是指在整個(gè)細(xì)胞分裂周期中染色質(zhì)的位置和結(jié)構(gòu)發(fā)生改變[36]。染色質(zhì)重塑發(fā)生時(shí)，在核小體連接處染色質(zhì)狀態(tài)變得相對(duì)疏松，這使得啟動(dòng)子區(qū)中的順式作用元件暴露，為順式作用元件與轉(zhuǎn)錄因子結(jié)合提供盡可能多的接觸空間。染色質(zhì)重塑與部分依賴能量供給的組蛋白修飾密切聯(lián)系，會(huì)誘導(dǎo)許多蛋白－DNA及蛋白－蛋白間的關(guān)聯(lián)受到破壞。2種蛋白復(fù)合體 (ATP依賴型核小體重構(gòu)復(fù)合體和組蛋白修飾復(fù)合體)會(huì)通過水解作用誘導(dǎo)核小體的構(gòu)型發(fā)生變化[37]，后者可以催化核心組蛋白N端尾部的共價(jià)修飾 (賴氨酸的乙?；胺核鼗?、精氨酸和賴氨酸的甲基化、谷氨酸的多聚ADP核糖基化等)[38]。通過組蛋白修飾酶的作用破壞基因組DNA與組蛋白尾部及核小體之間的相互作用也可以引起染色質(zhì)重塑[39－40]；另外，染色質(zhì)重塑也可為基因表達(dá)相關(guān)的蛋白及高級(jí)染色質(zhì)結(jié)構(gòu)的組織者提供識(shí)別位點(diǎn)，經(jīng)過修飾的組蛋白也可作為染色質(zhì)特異位點(diǎn)的標(biāo)志[23，41]。ATPase亞基存在于所有ATP依賴型核小體重構(gòu)復(fù)合體中，根據(jù)亞基的同源性可分為 ISWI、SWI/SNF、CHD、INO80、RAD54亞家族[42]。酵母交換型轉(zhuǎn)換/蔗糖不發(fā)酵復(fù)合物 (SWI/SNF)是第一個(gè)在釀酒酵母Saccharomyces cerevisiae中發(fā)現(xiàn)并被確認(rèn)的ATP依賴的重塑復(fù)合體[43]。

1.4 非編碼RNA調(diào)控

非編碼RNA是各種不翻譯成蛋白質(zhì)的RNA的統(tǒng)稱，非編碼RNA在基因的表達(dá)過程中同樣扮演重要的角色[44]，包括參與基因轉(zhuǎn)錄、翻譯等調(diào)控過程，催化RNA前體中內(nèi)含子的剪切、DNA復(fù)制、調(diào)控細(xì)胞發(fā)育和分化等[45－48](圖2)。

1.4.1 小分子RNA對(duì)基因轉(zhuǎn)錄的調(diào)控 小分子RNA主要有microRNA、siRNA、piRNA和啟動(dòng)子附近產(chǎn)生的一類RNA，micRNA和siRNA通過在細(xì)胞質(zhì)中組建RISC復(fù)合物誘導(dǎo)Argonaute蛋白與特定區(qū)域的靶標(biāo)RNA結(jié)合實(shí)現(xiàn)對(duì)基因的沉默，然而Argonaute蛋白在 RNA干擾 (RNA interference，RNAi)通路及在細(xì)胞核內(nèi)的存在形式在表觀遺傳調(diào)控上的作用為這套系統(tǒng)賦予了多層面、更復(fù)雜的功能。miR－320由基因POLR3D啟動(dòng)子區(qū)域反義方向轉(zhuǎn)錄產(chǎn)生，miR－320可以介導(dǎo) Argonaute－1(AGO1)蛋白與RRC2復(fù)合物組分EZH2之間的相互作用，也可以將啟動(dòng)子區(qū)域組蛋白H3上第27位賴氨酸的三甲基化 (H3K27me3)，從而抑制POLR3D基因的轉(zhuǎn)錄[49]。

1.4.2 長(zhǎng)鏈非編碼RNA對(duì)基因轉(zhuǎn)錄的調(diào)控 長(zhǎng)鏈非編碼RNA是一類不具有編譯蛋白質(zhì)能力的RNA，研究發(fā)現(xiàn)，影響DNA甲基化的RNA許多為編碼蛋白質(zhì)的mRNA的反義轉(zhuǎn)錄產(chǎn)物，如LUC7L和P53的反義RNA，WRAP53、HBA1(Haemoglobin α1)的反義RNA。WRAP53是通過與CTCF結(jié)合調(diào)節(jié)P53基因的表達(dá)，而LUC7L則是通過甲基化啟動(dòng)子區(qū)CPG島抑制基因轉(zhuǎn)錄。作為支架招募組蛋白修飾相關(guān)的復(fù)合物是長(zhǎng)鏈非編碼RNA調(diào)控基因轉(zhuǎn)錄的一種重要方式，通過對(duì)染色質(zhì)組蛋白進(jìn)行甲基化或乙?；揎椪{(diào)節(jié)基因的表達(dá)[50]。

2 環(huán)境脅迫對(duì)養(yǎng)殖魚類表觀遺傳的影響

水產(chǎn)養(yǎng)殖環(huán)境脅迫按性質(zhì)可分為:物理、化學(xué)、生物和人為因子等4種[51]。物理脅迫因子是由物理因素對(duì)機(jī)體造成的應(yīng)激，包括高溫、低溫及其溫度的劇變，以及光照、噪音等；化學(xué)脅迫因子源包括有機(jī)物和重金屬污染，亞硝酸鹽及氨氮濃度的過高，以及鹽度、pH變化等；生物脅迫因子包括養(yǎng)殖密度過高導(dǎo)致的擁擠應(yīng)激、營(yíng)養(yǎng)缺乏或過剩、魚類雌雄分化和生活習(xí)性等；人為脅迫因子包括養(yǎng)殖過程中人為的捕捉、分池、運(yùn)輸和病害處理等造成的應(yīng)激。

2.1 主要物理脅迫對(duì)水產(chǎn)養(yǎng)殖動(dòng)物的表觀遺傳效應(yīng)

2.1.1 溫度 魚類屬于變溫動(dòng)物，有適宜的生長(zhǎng)與繁育溫度，存在生命活動(dòng)的低溫與高溫極限。Liu等[52]研究發(fā)現(xiàn)，高溫應(yīng)激可顯著提高團(tuán)頭魴Megalobrama amblycephala血清促腎上腺激素和皮質(zhì)醇水平，并誘導(dǎo)氧化應(yīng)激，促進(jìn) HSP60、HSP70、HSP90基因的表達(dá)。在對(duì)歐洲鱸歐洲鱸Dicentrarchus labrax的研究中發(fā)現(xiàn)，在性腺成熟前溫度可以影響性別比例，在水溫較高情況下，雄魚比例大于雌魚，雄性鱸魚幼魚性腺內(nèi)的CYP19A啟動(dòng)子甲基化水平顯著高于雌性幼魚，在高溫刺激下，雌魚CYP19A基因的啟動(dòng)子甲基化水平明顯上升，說明DNA甲基化參與基因調(diào)控誘導(dǎo)性別分化，并且在該過程中CYP19A基因的表達(dá)與DNA甲基化水平呈現(xiàn)反比的關(guān)系[53]，個(gè)體CYP19A基因的啟動(dòng)子甲基化水平的變化過程發(fā)生在生殖腺而非其他組織中，說明這不是一個(gè)普通的溫度效應(yīng)，而且這些影響在性別未分化和未經(jīng)雌激素處理的魚中同樣可以檢測(cè)到，因此，在溫度誘導(dǎo)雄魚產(chǎn)生的過程中，CYP19A基因啟動(dòng)子甲基化扮演了重要的角色。在體外條件下，溫度會(huì)誘導(dǎo)歐洲鱸魚個(gè)體CYP19A基因啟動(dòng)子發(fā)生甲基化，SF－1和FoxI2基因的表達(dá)會(huì)受到抑制，推測(cè)CYP19A基因啟動(dòng)子甲基化可能是生物體在漫長(zhǎng)的進(jìn)化過程中，形成的一種對(duì)溫度適應(yīng)的機(jī)制[53]。Jabbari等[54]發(fā)現(xiàn)，變溫的魚類和兩棲類動(dòng)物DNA甲基化水平是恒溫的鳥類和哺乳類動(dòng)物的2倍，并且他還發(fā)現(xiàn)，南極魚類DNA甲基化水平高于熱帶魚類。有研究表明，魚類的肌纖維生成過程及魚類肌肉分化相關(guān)的基因表達(dá)與孵化溫度有密切的關(guān)系，斑馬魚Danio rerio胚胎在26℃孵化條件下，其成年個(gè)體肌纖維數(shù)量比低溫 (22℃)和高溫 (31℃)條件孵化的個(gè)體增加了14%～19%[55]。在對(duì)刺參Apostichopus japonicus的研究中發(fā)現(xiàn)，高溫應(yīng)激會(huì)促進(jìn) DNMT1、HDAC3和MLL5基因的表達(dá)[56]。河豚魚胚胎孵化后期，溫度升高會(huì)促進(jìn)FoxK1基因轉(zhuǎn)錄水平[57]。斑馬魚胚胎在23、27、31℃條件下孵化，選擇發(fā)育過程中的6個(gè)時(shí)間點(diǎn)檢測(cè)DNMTS基因的表達(dá)情況，結(jié)果顯示，DNMT3基因表達(dá)在發(fā)育前期有顯著的動(dòng)態(tài)變化趨勢(shì)，尤其是DNMT3A的mRNA的轉(zhuǎn)錄水平在孵化時(shí)期顯著增加，同時(shí)在相同發(fā)育階段，發(fā)現(xiàn)不同孵化溫度組中DNMT3A和DNMT3B基因表達(dá)有顯著差異，由此推測(cè)，魚類發(fā)育前期階段，DNMT3A和DNMT3B基因在溫度的表觀遺傳調(diào)控中發(fā)揮著不同的作用[58]。羅非魚 Oreochromis mossambicus耐寒品系DNA甲基化總體水平相比于對(duì)照組基因組明顯下降，經(jīng)過連續(xù)多代的低溫脅迫可引起尼羅羅非魚DNA甲基化水平及狀態(tài)發(fā)生改變，這說明羅非魚抗寒性能與DNA甲基化密切相關(guān)[59]。

總的來說，溫度是水產(chǎn)養(yǎng)殖過程中的主要脅迫因子，長(zhǎng)期的進(jìn)化過程中魚類已衍化出抵御溫度應(yīng)激的生理調(diào)節(jié)機(jī)制。綜上所述，機(jī)體可能會(huì)通過表觀遺傳修飾方式的變化以適應(yīng)溫度應(yīng)激，通過進(jìn)一步研究溫度脅迫對(duì)魚類表觀遺傳的影響，有助于尋找與環(huán)境適應(yīng)力相關(guān)的重要基因標(biāo)志物，從而解釋表觀遺傳的作用機(jī)理，為表觀遺傳學(xué)在魚類抗逆育種及其抗脅迫中的應(yīng)用提供新的思路和理論參考。2.1.2 光照 在自然界，所有生物的能量均來源于光照，光的變化具有穩(wěn)定性和規(guī)律性，光的周期性變化可影響魚類的繁殖、代謝等生理機(jī)能。已有研究表明，光的周期性變化可通過表觀遺傳影響魚類的表型。干旱期間，低耗能狀態(tài)和極端溫度下穴居的有條紋青蛙骨骼肌基因的轉(zhuǎn)錄沉默，DNA甲基轉(zhuǎn)移酶和DNA甲基化顯著擴(kuò)增。夏眠成年蛙DNA甲基轉(zhuǎn)移酶DNMT1的轉(zhuǎn)錄水平和轉(zhuǎn)錄協(xié)同抑制因子SIN3A顯著高于非成年蛙。與此類似，光周期性變化通過DNA甲基化調(diào)節(jié)肌肉生長(zhǎng)速率的機(jī)制在大西洋鱈魚Gadus morhus中也有報(bào)道，連續(xù)暴露在光照下的大西洋鱈肌肉中DNMT1和DNMT3A基因的表達(dá)量顯著高于自然光周期下飼養(yǎng)的魚。這種現(xiàn)象會(huì)一直維持下去，直到自然光周期等于連續(xù)照明的光周期，這表明，飼養(yǎng)期間連續(xù)的光照對(duì)于DNMT基因的表達(dá)影響可能是持久的[60]。目前，光照對(duì)魚類表觀遺傳具體作用機(jī)制還有待進(jìn)一步研究。

2.2 主要化學(xué)脅迫對(duì)水產(chǎn)動(dòng)物的表觀遺傳效應(yīng)

2.2.1 有機(jī)污染物 在集約化的高密度水產(chǎn)養(yǎng)殖環(huán)境中，環(huán)境污染物的存在是不可避免的。已有研究表明，多數(shù)環(huán)境污染物能通過基因組水平調(diào)控表觀遺傳機(jī)制，如調(diào)控基因啟動(dòng)子甲基化水平和組蛋白修飾等，這在生物生長(zhǎng)、發(fā)育過程中具有重要的作用。這種表觀遺傳狀態(tài)的變化可以通過生殖系細(xì)胞中的一系列調(diào)控從親本傳遞到生殖細(xì)胞，最后引起后代表型的變化。已有大量研究表明，生物生長(zhǎng)后期相關(guān)基因的疾病是由基因組甲基化狀態(tài)的變化引起。苯菌靈 (Benomy1)處理后的斑馬魚Danio rerio胚胎將會(huì)有一半死亡，存活的胚胎也會(huì)出現(xiàn)畸形、心臟伸長(zhǎng)、脊柱彎曲、心包水腫、心率降低等癥狀[61]。用雌二醇 (E2)處理后的雄性三刺魚性腺基因組會(huì)處于高甲基化狀態(tài)[7]。水蚤暴露于烯菌酮中DNA甲基化水平會(huì)發(fā)生持續(xù)穩(wěn)定的降低，且這種變化會(huì)遺傳給后代[62]。用烯菌酮處理F0代小鼠，可在F3代小鼠的生殖細(xì)胞中檢測(cè)到DNA甲基化狀態(tài)的改變[63]。暴露17α－炔雌醇一段時(shí)間后的成年斑馬魚，其肝臟卵黃蛋白原Ⅰ基因的5′側(cè)翼甲基化水平明顯下降，這表明雌激素可通過誘導(dǎo)DNA甲基化水平調(diào)控軟黃蛋白原的表達(dá)[64]。暴露于E2或BPA環(huán)境下的小鼠，其前列腺4型磷酸二酯酶變體4酶基因會(huì)處于低甲基化狀態(tài)，這表明發(fā)育時(shí)期小鼠的前列腺表觀基因組會(huì)受到低濃度雌激素的影響[65]。

野生魚在污染水域中性腺增長(zhǎng)相關(guān)基因的甲基化水平會(huì)異常升高，且其性腺增長(zhǎng)速度顯著低于在正常水環(huán)境中的野生魚[66]。在阿特拉津和毒死蜱農(nóng)藥脅迫下，鯉Cyprinus carpio的性腺組織和腦甲基化水平也會(huì)隨之發(fā)生改變[67]。先前有研究表明，當(dāng)環(huán)境中的苯并 [a]芘濃度達(dá)到1 μg/L時(shí)能夠傳代影響F3代青鳉Oryzias latipes幼魚的骨骼發(fā)育，苯并 [a]芘處理后的子三代雄性成魚通過Ct測(cè)量和組織病理學(xué)分析發(fā)現(xiàn)，骨骼的厚度減少，微裂紋增加，通過比對(duì)青鳉魚對(duì)照組和處理組miRNA的表達(dá)，檢測(cè)到5種方式mRNA/miRNA(Osx/Mir－214， Col2a1b/Mir－29b， Runx2/miR－204， Sox9b/miR－199a－3p， Apc/Mir－27b)[68]， 這表明苯并[a]芘能通過miRNA調(diào)節(jié)表觀遺傳機(jī)制。有研究表明，Cr(VI)可與組蛋白去乙?；?－DNA甲基化轉(zhuǎn)移酶1的復(fù)合物 (HDAC－DNMTL)相結(jié)合，通過交聯(lián)到Cypla1啟動(dòng)子區(qū)域的染色質(zhì)的方式，抑制由AhR調(diào)節(jié)的基因轉(zhuǎn)錄，引起組蛋白H3甲基化、磷酸化等，從而促進(jìn)了B[a]P和DNA形成加合物，導(dǎo)致DNA損傷[69](圖3)。有機(jī)污染物對(duì)魚類表觀遺傳影響的數(shù)據(jù)越來越多，并引起廣大學(xué)者的重視，更多的作用機(jī)制還有待進(jìn)一步研究。2.2.2 重金屬 環(huán)境中的重金屬，如鎘、鎳、鉻和甲基汞等可以引起機(jī)體發(fā)生氧化還原反應(yīng)，產(chǎn)生大量氧自由基，可氧化損傷DNA，干擾DNMTs與DNA的相互作用，從而影響DNA甲基化狀態(tài)[70]。有研究表明，毛蚶Scapharca subcrenata在經(jīng)過Cd2＋和Cu2＋處理之后，鰓組織甲基化水平異常升高[71]。也有研究顯示，在經(jīng)砷暴露數(shù)周后的大鼠Rattus norvegicus肝臟全基因組DNA呈低甲基化狀態(tài)[72?。同樣，當(dāng)金魚Carassius auratus經(jīng)砷暴露數(shù)周后，肝臟全基因組呈高甲基化狀態(tài)[73]。周新文等[74]研究發(fā)現(xiàn)，鯽魚Carassius auratus在經(jīng)0.02 mg/L的Pb2＋染毒48 h后，腎臟和肝臟中DNA總甲基化顯著提高。王丙蓮等[75]研究發(fā)現(xiàn)，泥鰍 Misgurnus anguillicaudatus在經(jīng) Pb2＋染毒 14 d后，肝胰臟DNA甲基化水平明顯升高，但在染毒28 d后，腎臟DNA甲基化水平異常降低。將圍產(chǎn)期大鼠經(jīng)飲水暴露乙酸鉛2個(gè)月后發(fā)現(xiàn)，新生小鼠海馬體全基因組呈高甲基化狀態(tài)[76]。在用鎘對(duì)歐洲鰻魚Anguilla anguilla的慢性脅迫試驗(yàn)中發(fā)現(xiàn)，相關(guān)基因通過甲基化水平的變化對(duì)基因表達(dá)產(chǎn)生負(fù)調(diào)控[77]。通過研究重金屬對(duì)水生動(dòng)物表觀遺傳效應(yīng)的影響，闡明相關(guān)的作用機(jī)制，可為重金屬污染防治提供一些參考依據(jù)。

2.2.3 鹽度 鹽度對(duì)于水生動(dòng)物在淡水湖泊、河流和海洋間的遷移是十分重要的。洄游魚類必須通過調(diào)節(jié)基因的表達(dá)途徑，以適應(yīng)從出生地到進(jìn)入海洋 (或從海洋返回產(chǎn)卵場(chǎng))間不斷變化的鹽度。從孵化場(chǎng)遷往海洋的魚存活率往往較低，這可能與適應(yīng)孵化環(huán)境有關(guān)。在對(duì)鮭魚Oncorhynchus的研究中發(fā)現(xiàn)，不同生長(zhǎng)條件的體質(zhì)下降現(xiàn)象與基因組的甲基化有關(guān)[78]，但是這種變化在淡水孵化場(chǎng)養(yǎng)殖的魚中未被發(fā)現(xiàn)[79]。在對(duì)三疣梭子蟹Portunus trituberculatus的低鹽脅迫試驗(yàn)中發(fā)現(xiàn)，其鰓中有12個(gè)差異表達(dá)的mRNA，這些mRNA的靶基因主要參與甲殼素的代謝及離子轉(zhuǎn)運(yùn)等滲透壓的調(diào)節(jié)過程[80]。滲透脅迫轉(zhuǎn)錄因子1(Osmotic Stress Transcription Factor 1，OSTF1)在水生動(dòng)物的滲透壓調(diào)節(jié)過程中起著重要作用。Yan等[81]研究發(fā)現(xiàn)，尼羅羅非魚miR－429可通過調(diào)控OSTF1的表達(dá)來調(diào)控血漿中的離子濃度。Flynt等[82]研究發(fā)現(xiàn)，在斑馬魚Danio rerio的胚胎發(fā)育過程中，miR－8家族可通過調(diào)節(jié)Nherf1基因的表達(dá)來調(diào)控離子的轉(zhuǎn)運(yùn)，并且在離子細(xì)胞中高度表達(dá)。

圖1 組蛋白修飾[24]Fig.1 Histone tail modifications[24]

圖2 非編碼RNA的調(diào)控機(jī)制[49]Fig.2 Regulation of gene transcription by non－coding RNAs[49]

圖3 B[a]P通過與DNA共價(jià)形成DNA加合物導(dǎo)致基因損傷[71]Fig.3 B[a]P can cause gene damage by covalent binding of DNA adducts with DNA[71]

一些廣鹽性的水生動(dòng)物如凡納濱對(duì)蝦Litopenaeus vannamei、羅非魚、梭魚Mugil soiuy等在淡水和海水中均可進(jìn)行人工養(yǎng)殖，目前對(duì)其滲透壓的調(diào)節(jié)機(jī)制尚不明確，從上述研究結(jié)果來看，鹽度的變化能影響水生動(dòng)物的一些miRNA及其相關(guān)基因的表達(dá)。針對(duì)具體物種的滲透壓調(diào)節(jié)機(jī)制，分析養(yǎng)殖水體環(huán)境的水質(zhì)、離子等成分可在必要時(shí)補(bǔ)充相關(guān)離子，對(duì)水生動(dòng)物的生長(zhǎng)具有重要的意義。因此，研究鹽度對(duì)水產(chǎn)動(dòng)物的表觀遺傳調(diào)控機(jī)制顯得尤為重要，相關(guān)內(nèi)容還待進(jìn)一步地解析。

2.2.4 亞硝酸鹽、氨氮 亞硝酸鹽和氨氮是普遍存在于水環(huán)境中的污染物，亞硝酸鹽由細(xì)菌消化過程和反硝化過程產(chǎn)生，氨氮是水產(chǎn)動(dòng)物蛋白質(zhì)代謝的最終產(chǎn)物。當(dāng)水體中的亞硝酸鹽和氨氮濃度過高時(shí)，會(huì)降低水產(chǎn)動(dòng)物的生長(zhǎng)率、免疫力，并造成其生理功能紊亂，甚至導(dǎo)致死亡。Guo等[83]對(duì)凡納濱對(duì)蝦的研究發(fā)現(xiàn)，亞硝酸鹽暴露顯著上調(diào)caspase－3 mRNA水平，同時(shí)誘導(dǎo)血細(xì)胞凋亡。Sun等[84]對(duì)團(tuán)頭魴的研究發(fā)現(xiàn)，亞硝酸鹽暴露誘導(dǎo)鰓細(xì)胞凋亡，并伴隨著caspase－8和caspase－9 mRNA的上調(diào)。腫瘤抑制因子P53參與多種生物學(xué)過程，包括DNA損傷修復(fù)，在對(duì)暗紋東方鲀Takifugu obscurus的研究中發(fā)現(xiàn)，在氨氮應(yīng)激下，機(jī)體會(huì)通過誘導(dǎo)NO和ROS的產(chǎn)生并上調(diào)P53的表達(dá)量來對(duì)抗應(yīng)激帶來的損傷[85]。目前，關(guān)于亞硝酸鹽、氨氮對(duì)水產(chǎn)動(dòng)物的表觀遺傳機(jī)制的影響研究較少，但根據(jù)表觀遺傳的作用機(jī)制可以分析發(fā)現(xiàn)，表觀遺傳修飾可影響多種生物學(xué)過程，包括DNA的損傷修復(fù)，誘導(dǎo)相關(guān)抗氧化基因的上調(diào)，從而抵抗外界環(huán)境脅迫對(duì)水產(chǎn)動(dòng)物機(jī)體造成的損傷。開展此方面的研究有助于了解亞硝酸鹽、氨氮對(duì)水生動(dòng)物表觀遺傳修飾的作用規(guī)律，為揭示水產(chǎn)動(dòng)物氧化應(yīng)激調(diào)節(jié)機(jī)制具有重要意義。

2.3 主要生物應(yīng)激對(duì)水生動(dòng)物的表觀遺傳效應(yīng)

2.3.1 營(yíng)養(yǎng)素 許多營(yíng)養(yǎng)刺激可以通過改變特定細(xì)胞的表觀遺傳標(biāo)記，影響魚類發(fā)育過程中細(xì)胞的分化，目前，營(yíng)養(yǎng)與表觀遺傳學(xué)的研究也受到廣大學(xué)者的關(guān)注。虹鱒Oncorhynchus mykiss的幼體發(fā)育需要經(jīng)歷內(nèi)源性營(yíng)養(yǎng)、混合營(yíng)養(yǎng)和外源性營(yíng)養(yǎng)3個(gè)階段，有研究發(fā)現(xiàn)，虹鱒肌肉中miR－1/133在這3個(gè)階段中特異表達(dá)，且表達(dá)量逐漸下降，這表明隨著魚類的食性變化，miR－1/133會(huì)通過調(diào)控肌細(xì)胞的形成參與魚類的生長(zhǎng)過程[86]。對(duì)禁食和飽食后的鱖Siniperca chuatsi研究發(fā)現(xiàn)，miR－10c等7個(gè)miRNA在肌肉中差異表達(dá)，這表明食物攝取能通過miRNA調(diào)控骨骼肌的發(fā)育[87]。有研究發(fā)現(xiàn)，飼喂高碳水化合物飼料能使虹鱒患有糖尿病，通過檢測(cè)mRNA轉(zhuǎn)錄水平發(fā)現(xiàn)，H3K9處于低乙?；癄顟B(tài)、G6pcb2處于高甲基化狀態(tài)[88]。缺乏維生素(A、D、C)的鱸魚其生長(zhǎng)因子相關(guān)基因表達(dá)的時(shí)序和協(xié)調(diào)性會(huì)受到破壞，影響成骨細(xì)胞的分化，從而使某些細(xì)胞轉(zhuǎn)化成脂肪細(xì)胞，導(dǎo)致魚體畸形[88]。成骨細(xì)胞轉(zhuǎn)化成脂肪細(xì)胞也可以發(fā)生在海洋魚類早期生命階段膳食中不飽和脂肪酸較高的情況下[89]。鮭魚脂肪前體細(xì)胞具有根據(jù)環(huán)境條件分化免疫譜系細(xì)胞 (最有可能是類巨噬細(xì)胞)的能力[90]。目前，日糧蛋白質(zhì)缺乏對(duì)機(jī)體表觀遺傳影響的報(bào)道相對(duì)較多，Rees等[91]研究發(fā)現(xiàn)，給妊娠期小鼠飼喂低蛋白質(zhì)的日糧會(huì)導(dǎo)致子代肝臟組織DNA高度甲基化，這也是動(dòng)物胚胎時(shí)期營(yíng)養(yǎng)缺乏與子代表觀遺傳修飾變化關(guān)系的較早報(bào)道。有諸多試驗(yàn)表明，在動(dòng)物妊娠期間，日糧蛋白質(zhì)缺乏會(huì)導(dǎo)致子代DNA某些核心區(qū)域甲基化模式發(fā)生改變，且這種變化在動(dòng)物成年后會(huì)變得非常穩(wěn)定。總的來說，現(xiàn)有研究表明，母體蛋白質(zhì)缺乏會(huì)通過表觀遺傳修飾的方式影響子代基因的表達(dá)。Burdge等[92]研究表明，脂肪酸，尤其是多不飽和脂肪酸 (n－3 HUFA's)可以改變表觀基因組。n－3 HUFA是類花生酸類的重要前體物質(zhì)，對(duì)魚類的生長(zhǎng)和發(fā)育至關(guān)重要，其以ATP的形式為魚類代謝提供能量，且參與細(xì)胞膜的形成[93]。雌性鮭魚性腺發(fā)育過程中，日糧中的脂肪源發(fā)生改變，會(huì)影響魚卵的脂肪酸成分，結(jié)果導(dǎo)致胚胎發(fā)育異常[94]。增加真鯛Pagrosomus major親魚日糧中的魚油添加比例，會(huì)影響孵化3 d后仔魚的存活率、魚苗的長(zhǎng)度，卵的直徑，以及蛋白質(zhì)和脂肪含量，并影響脂肪細(xì)胞的代謝和形態(tài)[95－96]。對(duì)人類的研究表明，妊娠期α－亞麻酸的攝取會(huì)誘發(fā)后代肝臟出現(xiàn)表觀遺傳學(xué)變異[97]。Benatti等[98]研究發(fā)現(xiàn)，一些脂肪酸會(huì)作為激素控制轉(zhuǎn)錄因子的活性。Seierstad等[99]研究發(fā)現(xiàn)，飼喂不同脂肪源的大西洋鮭Salmo salar長(zhǎng)至3～5 kg時(shí)，就已經(jīng)出現(xiàn)與典型生活方式相關(guān)的癥狀，如動(dòng)脈硬化、肝臟和心臟等器官中的脂肪異常沉積，導(dǎo)致魚類器官功能受損和各種循環(huán)問題。

日糧營(yíng)養(yǎng)與有機(jī)體的關(guān)系一直是學(xué)者重點(diǎn)研究的方向，有關(guān)日糧營(yíng)養(yǎng)素對(duì)動(dòng)物DNA甲基化的影響已有較多的研究，甲基供體對(duì)DNA甲基化的調(diào)控也受到了較多關(guān)注，如甜菜堿對(duì)肉雞脂蛋白酯酶(LPL)[100]和脂肪合成酶 (FAS)[101]啟動(dòng)子甲基化的影響。研究營(yíng)養(yǎng)素對(duì)水產(chǎn)動(dòng)物表觀遺傳影響的分子機(jī)制，有助于從源頭上闡述營(yíng)養(yǎng)素對(duì)魚類的影響，但目前營(yíng)養(yǎng)素對(duì)魚類表觀遺傳調(diào)控的影響研究剛剛起步，隨著研究的逐步深入，其中的具體調(diào)控機(jī)制將會(huì)逐步明確。

2.3.2 雌、雄個(gè)體 有研究表明，雌魚哺育后代能通過表觀遺傳影響后代，通過對(duì)三刺魚Gasterosteus aculeatus的研究，雄魚哺育后代同樣能增加后代對(duì)早期環(huán)境壓力的抗性，提高后代成活率，而這種變化與腦組織中DNA甲基化密切相關(guān)[102]。歐洲鱸魚的雄魚相比于雌魚芳香化酶啟動(dòng)子表現(xiàn)出較高水平的DNA甲基化[103]。類似的研究在兩棲類動(dòng)物紅耳龜Trachemys scripta和密西西比短吻鱷Alligator mississippiensis中也有報(bào)道[104]。一些種類的海龜 Chelonia mydas、蜥蜴和鱷魚 Crocodylus siamensis的性別會(huì)受到溫度的影響，在對(duì)鱷魚的研究中發(fā)現(xiàn)，雄性性腺的SOX9基因的啟動(dòng)子序列基因相比于雌性表現(xiàn)出較高的DNA甲基化水平[105]。在遺傳育種方面，紅鯽與鯉的異源四倍體雜交過程同樣被證實(shí)受到DNA甲基化調(diào)控[106]。研究表明，CYP11a1基因在卵巢和精巢發(fā)育過程中的表達(dá)規(guī)律不同，在精巢成熟期表達(dá)量最高，CYP11a1基因在一定程度上與斑馬魚的性腺發(fā)育相關(guān)，其可能在促進(jìn)斑馬魚精子成熟中發(fā)揮一定作用[107]。不同種魚類有不同的生活習(xí)性及生長(zhǎng)特性，通過研究其與表觀遺傳的作用機(jī)制，對(duì)于理解表觀遺傳機(jī)制的多功能性和解析水產(chǎn)動(dòng)物生長(zhǎng)特性均具有重要意義。

2.3.3 養(yǎng)殖密度 在集約化水產(chǎn)養(yǎng)殖中，養(yǎng)殖密度的提高可有效增加養(yǎng)殖效益，但隨著養(yǎng)殖密度的提高，魚類也受到了更多的生存壓力。Sadhu等[108]研究發(fā)現(xiàn)，海鱸Lateolabrax japonicus血糖含量在高密度組顯著高于低密度組，Xie等[109]研究發(fā)現(xiàn)，鯉的血液皮質(zhì)醇、葡萄糖、溶菌酶和肝丙二醛含量在擁擠脅迫后顯著升高，肝臟中過氧化氫酶和超氧化物歧化酶活性在擁擠脅迫后減少。Costas等[110]研究發(fā)現(xiàn)，塞內(nèi)加爾舌鰨Solea senegalensis養(yǎng)殖63 d后，高密度組中魚血漿滲透壓顯著提高，這表明過高的養(yǎng)殖密度會(huì)導(dǎo)致機(jī)體滲透壓功能和離子平衡的改變。Caipang等[111]研究表明，大西洋鱈在擁擠脅迫24 h和72 h后，HSP70基因表達(dá)顯著上調(diào)。在對(duì)養(yǎng)殖密度的研究中發(fā)現(xiàn)，高密度組的海鱸和虹鱒HSP70基因表達(dá)顯著高于低密度組[112－113]。而對(duì)塞內(nèi)加爾鰨的研究中，高密度組魚HSP70和HSP90基因的表達(dá)顯著下調(diào)[114]。Gornati等[112]研究發(fā)現(xiàn)，海鱸 MT mRNA水平在高密度養(yǎng)殖條件下顯著上調(diào)，表明MT可作為一種應(yīng)激指標(biāo)來評(píng)估機(jī)體的生理狀態(tài)。CYP1A是混合功能氧化酶，可作生物標(biāo)志物，用來評(píng)估機(jī)體生理狀態(tài)。研究發(fā)現(xiàn)，海鱸肝、脾和腎CYP1A基因在高密度養(yǎng)殖條件下均顯著上調(diào)[112]。促性腺激素釋放激素GnRH1是一種十肽結(jié)構(gòu)的神經(jīng)激素，具有刺激腺垂體釋放促性腺激素 (FSH)和黃體生成素 (LH)釋放的功能，控制著脊椎動(dòng)物的繁殖[115]。在對(duì)非洲慈鯛的研究中發(fā)現(xiàn)，妊娠期雌魚正常生長(zhǎng)情況下和擁擠脅迫條件下腦組織中GnRH1基因甲基化程度在一些CG位點(diǎn)明顯不同，且其后代中CpG島－1146、－1128和－1110位點(diǎn)呈低甲基化狀態(tài)[116]。在對(duì)虹鱒的研究發(fā)現(xiàn)，高密度能顯著提高其鰓和腎臟中PYGMA基因的表達(dá)，使肝臟中HIF1A信號(hào)通路上調(diào)，且與擁擠脅迫下血糖濃度的變化相一致[117]。從以上結(jié)果可以發(fā)現(xiàn)，擁擠脅迫可以影響魚的生長(zhǎng)，以及代謝酶和基因的表達(dá)。雖然某些基因變化的具體調(diào)節(jié)機(jī)制尚不明確，但可以明確的是，基因表達(dá)水平的變化是通過外界環(huán)境刺激誘導(dǎo)表觀遺傳修飾，進(jìn)而影響基因的表達(dá)。目前，密度應(yīng)激與酶及應(yīng)激相關(guān)信號(hào)通路 (如Nrf2－Keap1)的關(guān)系已有報(bào)道，但與表觀遺傳變化的具體調(diào)節(jié)機(jī)制還有待研究。

2.4 人為應(yīng)激的表觀遺傳效應(yīng)

管理脅迫是指在水產(chǎn)養(yǎng)殖過程中由于人為管理操作而對(duì)水產(chǎn)動(dòng)物所造成的脅迫。高強(qiáng)度的管理脅迫會(huì)對(duì)養(yǎng)殖動(dòng)物的生理狀態(tài)、繁殖、免疫力和生長(zhǎng)造成不利影響。如在拉網(wǎng)、分池、捕捉、苗種運(yùn)輸、麻醉、消毒、病害用藥等。生物前期生活經(jīng)歷的環(huán)境效應(yīng)和前期母體效應(yīng)將會(huì)影響各種精神機(jī)能的發(fā)育。對(duì)大鼠的研究表明，母體效應(yīng)的變化信息能夠穩(wěn)定的影響后代DNA甲基化[118]。研究發(fā)現(xiàn)，親代長(zhǎng)期受到電離輻射會(huì)通過生殖細(xì)胞將輻射危害信息傳遞給下一代，增加癌癥發(fā)病率和突變的可能性。研究人員推測(cè)，這可能是環(huán)境中的電離輻射會(huì)影響機(jī)體基因組DNA甲基化水平，且該影響具有相對(duì)穩(wěn)定的遺傳特性[118]。對(duì)表觀遺傳機(jī)制與人為應(yīng)激的研究較少，許多具體機(jī)制有待進(jìn)一步研究。

綜上所述，表觀遺傳修飾主要包括DNA甲基化、組蛋白修飾、染色質(zhì)重塑和非編碼RNA調(diào)控等4個(gè)方面。這4個(gè)方面既相互作用又相互獨(dú)立，共同調(diào)節(jié)著胚胎期及出生后整個(gè)生命周期中的基因表達(dá)。DNA甲基化是表觀遺傳修飾最主要的方式，近些年，在對(duì)哺乳動(dòng)物的研究中已獲得了較多數(shù)據(jù)支持，對(duì)水生生物的研究結(jié)果表明，有較多的應(yīng)激方式均可引起DNA甲基化的變化，如高溫可誘導(dǎo)鱸魚CYP19A啟動(dòng)子和刺參相關(guān)促甲基化酶的表達(dá)，光照、營(yíng)養(yǎng)成分、水環(huán)境變化等均會(huì)影響相關(guān)基因的甲基化程度。污染物B[a]P可引起組蛋白H3甲基化、磷酸化等，從而促進(jìn)了B[a]P和DNA形成加合物，影響基因的表達(dá)。目前，關(guān)于非編碼RNA對(duì)水生動(dòng)物表觀遺傳效應(yīng)的影響研究已有較多的報(bào)道，亞硝酸鹽、氨氮、鹽度的變化等均會(huì)引起相關(guān)非編碼RNA的變化。總之，水生動(dòng)物在生長(zhǎng)過程中會(huì)受到物理、化學(xué)、生物、人為等應(yīng)激因素影響，并通過DNA甲基化、組蛋白修飾、染色質(zhì)重塑和非編碼RNA調(diào)控等的表觀遺傳調(diào)控，最終影響水生動(dòng)物的生長(zhǎng)，以及發(fā)育障礙、疾病等表型性狀 (圖4)。

4 展望

外界環(huán)境刺激會(huì)影響動(dòng)物生長(zhǎng)發(fā)育，引起表現(xiàn)型發(fā)生改變，其未涉及遺傳信息的改變，這屬于表觀遺傳學(xué)的范疇。在宏觀上，生物體通過表現(xiàn)型的變化適應(yīng)生活環(huán)境的改變，在分子水平上，環(huán)境脅迫和表觀遺傳間的關(guān)系還有待深入研究。近年來，隨著蛋白質(zhì)組學(xué)技術(shù)和基因組測(cè)序技術(shù)的快速發(fā)展，研究基因調(diào)控和基因功能與環(huán)境應(yīng)激的關(guān)系已經(jīng)成為可能，如運(yùn)用二代測(cè)序技術(shù)可系統(tǒng)檢測(cè)水生動(dòng)物的基因組及轉(zhuǎn)錄組，但研究表觀基因組需要考慮的因素也許多，如研究處于什么發(fā)展階段的表觀基因組，判斷這些表觀遺傳變化是否遺傳就必須研究親本和子代早期生活階段的表觀基因組。在某些情況下，在魚的發(fā)育過程中進(jìn)行表觀遺傳標(biāo)記是必要的。通過對(duì)表觀遺傳機(jī)理的深入研究，能讓養(yǎng)殖工作者更多地了解環(huán)境因素對(duì)魚類健康的影響，從而針對(duì)不同生長(zhǎng)階段的水產(chǎn)動(dòng)物，設(shè)計(jì)更精細(xì)、更合理的養(yǎng)殖環(huán)境，使魚類的生產(chǎn)性能得到更大程度發(fā)揮的同時(shí)減少魚類患病幾率，提高經(jīng)濟(jì)效益，同時(shí)為育種學(xué)家提供新的思路，如尋找鑒定新的miRNA，并對(duì)靶基因進(jìn)行標(biāo)識(shí)，實(shí)現(xiàn)對(duì)水產(chǎn)動(dòng)物的組織定向分化?；?qū)NA甲基化用作分子遺傳標(biāo)記手段對(duì)目標(biāo)性狀進(jìn)行輔助選擇標(biāo)記，調(diào)控特異性基因?qū)崿F(xiàn)對(duì)性別的控制，從而改善水產(chǎn)動(dòng)物品質(zhì)，加快推動(dòng)集約化高效養(yǎng)殖的現(xiàn)代化漁業(yè)進(jìn)程。

圖4 表觀遺傳調(diào)節(jié)機(jī)制Fig.4 Epigenetic regulation mechanisms

參考文獻(xiàn):

[1] 麥康森，艾慶輝，徐瑋，等.水產(chǎn)養(yǎng)殖中的環(huán)境脅迫及其預(yù)防——營(yíng)養(yǎng)學(xué)途徑[J].中國(guó)海洋大學(xué)學(xué)報(bào)，2004，34(5):767－774.

[2] Deans C，Maggert K A.What do you mean，“epigenetic”?[J].Genetics，2015，199(4):887－896.

[3] Seo M K，Ly N N，Lee C H，et al.Early life stress increases stress vulnerability through BDNF gene epigenetic changes in the rat hippocampus[J].Neuropharmacology，2016，105:388－397.

[4] Weidman J R，Dolinoy D C，Murphy S K，et al.Cancer susceptibility:epigenetic manifestation of environmental exposures[J].Cancer Journal，2007，13(1):9－16.

[5] Goldberg A D，Allis C D，Bernstein E.Epigenetics:a landscape takes shape[J].Cell，2007，128(4):635－638.

[6] Cox B D.On the difficulty in getting out of historical ruts:waddington and an argument for behavioral epigenetics[J].New Ideas in Psychology，2013，31(3):374－389.

[7] Aniagu S O，Williams T D，Allen Y，et al.Global genomic methylation levels in the liver and gonads of the three－spine stickleback(Gasterosteus aculeatus)after exposure to hexabromocyclododecane and 17－β oestradiol[J].Environment International，2008，34(3):310－317.

[8] Nelissen E C M，van Montfoort A P A，Dumoulin J C M.Epigenetics and the placenta[J].Human Reproduction Update，2010，17(3):397－417.

[9] 王杰，徐友信，刁其玉，等.非孟德爾遺傳模式:表觀遺傳學(xué)及其應(yīng)用研究進(jìn)展[J].中國(guó)農(nóng)學(xué)通報(bào)，2016，32(14):37－43.

[10] Kurdyukov S，Bullock M.DNA methylation analysis:choosing the right method[J].Biology，2016，5(1):3.

[11] Si Yufeng，Ding Yuxia，He Feng，et al.DNA methylation level of cyp19a1a and Foxl2 gene related to their expression patterns and reproduction traits during ovary development stages of Japanese flounder(Paralichthys olivaceus)[J].Gene，2016，575(2):321－330.

[12] Rasmussen K D，Helin K.Role of TET enzymes in DNA methylation，development，and cancer[J].Genes ＆ Development，2016，30(7):733－750.

[13] Kohli R M，Zhang Yi.TET enzymes，TDG and the dynamics of DNA demethylation[J].Nature，2013，502(7472):472－479.

[14] Gong Zhizhong，Morales－Ruiz T，Ariza R R，et al.ROS1，a re-pressor of transcriptional gene silencing in Arabidopsis，encodes a DNA glycosylase/lyase[J].Cell，2002，111(6):803－814.

[15] Li Xiaojie，Qian Weiqiang，Zhao Yusheng，et al.Antisilencing role of the RNA－directed DNA methylation pathway and a histone acetyltransferase in Arabidopsis[J].Proceedings of the National Academy of Sciences of the United States of America，2012，109(28):11425－11430.

[16] Qian Weiqiang，Miki D，Zhang Heng，et al.A histone acetyltransferase regulates active DNA demethylation in Arabidopsis[J].Science，2012，336(6087):1445－1448.

[17] Zhao Yusheng，Xie Shaojun，Li Xiaojie，et al.Repressor of silencing5 encodes a member of the small heat shock protein family and is required for DNA demethylation in arabidopsis[J].Plant Cell，2014，26(6):2660－2675.

[18] Qian Weiqiang，Miki D，Lei Mingguang，et al.Regulation of active DNA demethylation by an α－crystallin domain protein in Arabidopsis[J].Molecular Cell，2014，55(3):361－371.

[19] Li Yan，Córdoba－Ca?ero D，Qian Weiqiang，et al.An AP endonuclease functions in active DNA demethylation and gene imprinting in Arabidopsis[J].PLoS Genetics，2015，11(1):e1004905.

[20] Lee J，Jang H，Shin H，et al.AP endonucleases process 5－methylcytosine excision intermediates during active DNA demethylation in Arabidopsis[J].Nucleic Acids Research，2014，42(18):11408－11418.

[21] Lang Zhaobo，Lei Mingguang，Wang Xingang，et al.The Methyl－CpG－binding protein MBD7 facilitates active DNA demethylation to limit DNA hyper－methylation and transcriptional gene silencing[J].Molecular Cell，2015，57(6):971－983.

[22] Li Qi，Wang Xiaokang，Sun Han，et al.Regulation of active DNA demethylation by a methyl－CpG－binding domain protein in Arabidopsis thaliana[J].PLoS Genetics，2015，11(5):e1005210.

[23] Kouzarides T.Chromatin modifications and their function[J].Cell，2007，128(4):693－705.

[24] Jiang C，Pugh B F.Nucleosome positioning and gene regulation:advances through genomics[J].Nature Reviews Genetics，2009，10(3):161－172.

[25] Sun Weijian，Zhou Xiang，Zheng Jihang，et al.Histone acetyltransferases and deacetylases:molecular and clinical implications to gastrointestinal carcinogenesis[J].Acta Biochimica et Biophysica Sinica，2012，44(1):80－91.

[26] Clements E G，Mohammad H P，Leadem B R，et al.DNMT1 modulates gene expression without its catalytic activity partially through its interactions with histone－modifying enzymes[J].Nucleic Acids Research，2012，40(10):4334－4346.

[27] Allfrey V G，F(xiàn)aulkner R，Mirsky A E.Acetylation and methylation of histones and their possible role in the regulation of RNA synthesis[J].Proceedings of the National Academy of Sciences of the United States of America，1964，51(5):786－794.

[28] Pogo B G T，Allfrey V G，Mirsky A E.RNA synthesis and histone acetylation during the course of gene activation in lymphocytes[J].Proceedings of the National Academy of Sciences of the United States of America，1966，55(4):805－812.

[29] Hon G C，Hawkins R D，Ren B.Predictive chromatin signatures in the mammalian genome[J].Human Molecular Genetics，2009，18(R2):R195－R201.

[30] Schmitges F W，Prusty A B，F(xiàn)aty M，et al.Histone methylation by PRC2 is inhibited by active chromatin marks[J].Molecular Cell，2011，42(3):330－341.

[31] Banerjee T，Chakravarti D.A peek into the complex realm of histone phosphorylation[J].Molecular and Cellular Biology，2011，31(24):4858－4873.

[32] He Yingzi，Tang Dongmei，Cai Chengfu，et al.LSD1 is required for hair cell regeneration in zebrafish[J].Molecular Neurobiology，2016，53(4):2421－2434.

[33] Wang Yang，Li Yuejiao，Guo Chen，et al.ISL1 and JMJD3 synergistically control cardiac differentiation of embryonic stem cells[J].Nucleic Acids Research，2016，44(14):6741－6755.

[34] Tse C，Sera T，Wolffe A P，et al.Disruption of higher－order folding by core histone acetylation dramatically enhances transcription of nucleosomal arrays by RNA polymerase III[J].Molecular and Cellular Biology，1998，18(8):4629－4638.

[35] Wang A，Kurdistani S K，Grunstein M.Requirement of Hos2 histone deacetylase for gene activity in yeast[J].Science，2002，298(5597):1412－1414.

[36] Riddihough G，Pennisi E.The evolution of epigenetics[J].Science，2001，293(5532):1063.

[37] Racki L R，Narlikar G J.ATP－dependent chromatin remodeling enzymes:two heads are not better，just different[J].Current Opinion in Genetics ＆ Development，2008，18(2):137－144.

[38] Wang G G，Allis C D，Chi Ping.Chromatin remodeling and cancer，part I:covalent histone modifications[J].Trends in Molecular Medicine，2007，13(9):363－372.

[39] Barrett R M，Wood M A.Beyond transcription factors:The role of chromatin modifying enzymes in regulating transcription required for memory[J].Learning ＆ Memory，2008，15(7):460－467.

[40] Taverna S D，Li Haitao，Ruthenburg A J，et al.How chromatinbinding modules interpret histone modifications:Lessons from professional pocket pickers[J].Nature Structural＆Molecular Biology，2007，14(11):1025－1040.

[41] Li Bing，Carey M，Workman J L.The role of chromatin during transcription[J].Cell，2007，128(4):707－719.

[42] Gangaraju V K，Bartholomew B.Mechanisms of ATP dependent chromatin remodeling[J].Mutation Research/Fundamental and Molecular Mechanisms of Mutagenesis，2007，618(1－2):3－17.

[43] Peterson C L，Herskowitz I.Characterization of the yeast SWI1，SWI2，and SWI3 genes，which encode a global activator of transcription[J].Cell，1992，68(3):573－583.

[44] Mattick J S，Makunin I V.Non－coding RNA.Human Molecular Genetics，2006，15(S1):R17－R29.

[45] Ma Lina，Bajic V B，Zhang Zhang.On the classification of long non－coding RNAs[J].RNA Biology，2013，10(6):924－933.

[46] Quinn J J，Chang H Y.Unique features of long non－coding RNA biogenesis and function[J].Nature Reviews Genetics，2016，17(1):47.

[47] Cech T R，Steitz1 J A.The noncoding RNA revolution—trashing old rules to forge new ones[J].Cell，2014，157(1):77－94.

[48] 劉旭，單革.非編碼RNA對(duì)基因轉(zhuǎn)錄的調(diào)控[J].生命科學(xué)，2016，28(5):569－575.

[49] Kim D H，Saetrom P，Sn?ve Jr O，et al.MicroRNA－directed transcriptional gene silencing in mammalian cells[J].Proceedings of the National Academy of Sciences of the United States of America，2008，105(42):16230－16235.

[50] Salda?a－Meyer R，González－Buendía E，Guerrero G，et al.CTCF regulates the human p53 gene through direct interaction with its natural antisense transcript，Wrap53[J].Genes ＆ Development，2014，28(7):723－734.

[51] Pickering A D.Environmental stress and the survival of brown trout，salmo trutta[J].Freshwater Biology，1989，21(1):47－55.

[52] Liu B，Xu P，Brown P B，et al.The effect of hyperthermia on liver histology，oxidative stress and disease resistance of the Wuchang bream，Megalobrama amblycephala[J].Fish＆ Shellfish Immunology，2016，52:317－324.

[53] Navarro－Martín L，Vi?as J，Ribas L，et al.DNA methylation of the gonadal aromatase(cyp19a)promoter is involved in temper ature－dependent sex ratio shifts in the european sea bass[J].PLoS Genetics，2011，7(12):e1002447.

[54] Jabbari K，Cacciò S，de Barros J P P，et al.Evolutionary changes in CpG and methylation levels in the genome of vertebrates[J].Gene，1998，205(1－2):109－118.

[55] Johnston I A，Lee H T，Macqueen D J，et al.Embryonic temperature affects muscle fibre recruitment in adult zebrafish:genomewide changes in gene and microRNA expression associated with the transition from hyperplastic to hypertrophic growth phenotypes[J].Journal of Experimental Biology，2009，212(Pt 12):1781－1793.

[56] 李尚俊，孫國(guó)華，李雪燕，等.高溫脅迫下仿刺參表觀遺傳調(diào)控相關(guān)基因的表達(dá)特征[J].中國(guó)水產(chǎn)科學(xué)，2017，24(3):470－476.

[57] Fernandes J M，MacKenzie M G，Kinghorn J R，et al.FoxK1 splice variants show developmental stage－specific plasticity of expression with temperature in the tiger pufferfish[J].Journal of Experimental Biology，2007，210(Pt 19):3461－3472.

[58] Campos C，Valente L M P，F(xiàn)ernandes J M O.Molecular evolution of zebrafish dnmt3 genes and thermal plasticity of their expression during embryonic development[J].Gene，2012，500(1):93－100.

[59] 朱華平，盧邁新，黃樟翰，等.低溫對(duì)羅非魚基因組DNA甲基化的影響[J].水產(chǎn)學(xué)報(bào)，2013，37(10):1460－1467.

[60] Giannetto A，Nagasawa K，F(xiàn)asulo S，et al.Influence of photoperiod on expression of DNA(cytosine－5)methyltransferases in Atlantic cod[J].Gene，2013，519(2):222－230.

[61] Kim D J，Seok S H，Baek M W，et al.Benomyl induction of brain aromatase and toxic effects in the zebrafish embryo[J].Journal of Applied Toxicology，2009，29(4):289－294.

[62] Vandegehuchte M B，Lemière F，Vanhaecke L，et al.Direct and transgenerational impact on Daphnia magna of chemicals with a known effect on DNA methylation[J].Comparative Biochemistry and Physiology Part C:Toxicology ＆ Pharmacology，2009，151(3):278－285.

[63] Guerrero－Bosagna C，Settles M，Lucker B，et al.Epigenetic transgenerational actions of vinclozolin on promoter regions of the sperm epigenome[J].PLoS One，2010，5(9):e13100.

[64] Str?mqvist M，Tooke N，Brunstr?m B.DNA methylation levels in the 5′flanking region of the vitellogenin I gene in liver and brain of adult zebrafish(Danio rerio)—Sex and tissue differences and effects of 17α－ethinylestradiol exposure[J].Aquatic Toxicology，2010，98(3):275－281.

[65] Prins G S，Tang W Y，Belmonte J，et al.Perinatal exposure to oestradiol and bisphenol a alters the prostate epigenome and increases susceptibility to carcinogenesis[J].Basic＆Clinical Pharmacology ＆ Toxicology，2008，102(2):134－138.

[66] Pierron F，Colombier S B D，Moffett A，et al.Abnormal ovarian DNA methylation programming during gonad maturation in wild contaminated fish[J].Environmental Science ＆ Technology，2014，48(19):11688－11695.

[67] Xing Houjuan，Wang Chao，Wu Hongda，et al.Effects of atrazine and chlorpyrifos on DNA methylation in the brain and gonad of the common carp[J].Comparative Biochemistry and Physiology Part C:Toxicology ＆ Pharmacology，2015，168:11－19.

[68] Seemann F，Jeong C B，Zhang Ge，et al.Ancestral benzo[a]pyrene exposure affects bone integrity in F3 adult fish(Oryzias latipes)[J].Aquatic Toxicology，2016，183:127－134.

[69] 夏菠.Cr(Ⅵ)和B[a]P對(duì)16HBE細(xì)胞的聯(lián)合毒性效應(yīng)及表觀遺傳改變的研究[D].長(zhǎng)沙:中南大學(xué)，2012.

[70] Baccarelli A，Bollati V.Epigenetics and environmental chemicals[J].Current Opinion in Pediatrics，2009，21(2):243－251.

[71] Lee Y W，Broday L，Costa M.Effects of nickel on DNA methyltransferase activity and genomic DNA methylation levels[J].Mutation Research/Genetic Toxicology and Environmental Mutagenesis，1998，415(3):213－218.

[72] Uthus E O，Davis C.Dietary arsenic affects dimethylhydrazine－induced aberrant crypt formation and hepatic global DNA methylation and DNA methyltransferase activity in rats[J].Biological Trace Element Research，2005，103(2):133－145.

[73] Bagnyukova T V，Luzhna L I，Pogribny I P，et al.Oxidative stress and antioxidant defenses in goldfish liver in response to shortterm exposure to arsenite[J].Environmental and Molecular Mutagenesis，2007，48(8):658－665.

[74] 周新文，朱國(guó)念，Mwalilino J，等.Cu、Zn、Pb、Cd 及其混合重金屬離子對(duì)鯽魚(Carassius auratus)DNA甲基化水平的影響[J].中國(guó)環(huán)境科學(xué)，2001，21(6):546－552.

[75] 王丙蓮，張迎梅，譚玉鳳，等.鎘鉛對(duì)泥鰍DNA甲基化水平的影響[J].毒理學(xué)雜志，2006，20(2):78－80.

[76] Sánchez－Martín F J，Lindquist D M，Landero－Figueroa J，et al.Sex－and tissue－specific methylome changes in brains of mice perinatally exposed to lead[J].NeuroToxicology，2015，46:92－100.

[77] Pierron F，Baillon L，Sow M，et al.Effect of low－dose cadmium exposure on DNA methylation in the endangered European eel[J].Environmental Science ＆ Technology，2014，48(1):797－803.

[78] Chilcote M W C，Goodson K W G，F(xiàn)alcy M R.Reduced recruitment performance in natural populations of anadromous salmonids associated with hatchery－reared fish[J].Canadian Journal of Fisheries and Aquatic Sciences，2011，68(3):511－522.

[79] Blouin M S，Thuillier V，Cooper B，et al.No evidence for large differences in genomic methylation between wild and hatchery steelhead(Oncorhynchus mykiss)[J].Canadian Journal of Fisheries and Aquatic Sciences，2010，67(2):217－224.

[80] Lv Jianjian，Liu Ping，Gao Baoquan，et al.The identification and characteristics of salinity－related microRNAs in gills of Portunus trituberculatus[J].Cell Stress and Chaperones，2016，21(1):63－74.

[81] Yan Biao，Zhao Lihui，Guo Jinliang，et al.miR－429 regulation of osmotic stress transcription factor 1(OSTF1)in tilapia during osmotic stress[J].Biochemical and Biophysical Research Communications，2012，426(3):294－298.

[82] Flynt A S，Thatcher E J，Burkewitz K，et al.miR－8 microRNAs regulate the response to osmotic stress in zebrafish embryos[J].The Journal of Cell Biology，2009，185(1):115－127.

[83] Guo Hui，Xian Jian'an，Li Bin，et al.Gene expression of apo ptosis－related genes， stress protein and antioxidant enzymes in hemocytes of white shrimp Litopenaeus vannamei under nitrite stress[J].Comparative Biochemistry and Physiology Part C:Toxicology ＆ Pharmacology，2013，157(4):366－371.

[84] Sun Shengming，Ge Xiangping，Zhu Jian，et al.Identification and mRNA expression of antioxidant enzyme genes associated with the oxidative stress response in the Wuchang bream(Megalobrama amblycephala Yih)in response to acute nitrite exposure[J].Comparative Biochemistry and Physiology Part C:Toxicology＆Pharmacology，2014，159:69－77.

[85] Cheng Changhong，Yang Fangfang，Ling Renzhi，et al.Effects of ammonia exposure on apoptosis，oxidative stress and immune response in pufferfish(Takifugu obscurus)[J].Aquatic Toxicology，2015，164:61－71.

[86] Mennigen J A，Skiba－Cassy S，Panserat S.Ontogenetic expression of metabolic genes and microRNAs in rainbow trout alevins during the transition from the endogenous to the exogenous feeding period[J].The Journal of Experimental Biology，2013，216(Pt 9):1597－1608.

[87] Zhu Xin，Chen Dunxun，Hu Yi，et al.The microRNA signature in response to nutrient restriction and refeeding in skeletal muscle of Chinese perch(Siniperca chuatsi)[J].Marine Biotechnology，2015，17(2):180－189.

[88] Mazurais D，Darias M J，Gouillou－Coustans M F，et al.Dietary vitamin mix levels influence the ossification process in European sea bass(Dicentrarchus labrax)larvae[J].American Journal of Physiology.Regulatory， Integrative and Comparative Physiology，2007，294(2):R520－R527.

[89] ?kugor S，?kugor A，Todorˇcevic'M，et al.Exposure to lipopolysaccharide induces immune genes in cultured preadipocytes of Atlantic salmon[J].Fish ＆ shellfish immunology，2010，29(5):817－824.

[90] Marandel L，Lepais O，Arbenoits E，et al.Remodelling of the hepatic epigenetic landscape of glucose－intolerant rainbow trout(Oncorhynchus mykiss)by nutritional status and dietary carbohydrates[J].Scientific Reports，2016，6:32187.

[91] Rees W D，Hay S M，Brown D S，et al.Maternal protein deficiency causes hypermethylation of DNA in the livers of rat fetuses[J].Journal of Nutrition，2000，130(7):1821－1826.

[92] Burdge G C，Lillycrop K A.Fatty acids and epigenetics[J].Current Opinion in Clinical Nutrition and Metabolic Care，2014，17(2):156－161.

[93] Sargent J R，Bell M V，Bell J G，et al.Origins and functions of n－3 polyunsaturated fatty acids in marine organisms[M].Champaign:Ceve G and Paltauf F，1995:248－259.

[94] Pickova J，Kiessling A，Pettersson A，et al.Fatty acid and carotenoid composition of eggs from two nonanadromous Atlantic salmon stocks of cultured and wild origin[J].Fish Physiology and Biochemistry，1999，21(2):147－156.

[95] Cruz－Garcia L，Sánchez－Gurmaches J，Bouraoui L，et al.Changes in adipocyte cell size，gene expression of lipid metabolism markers，and lipolytic responses induced by dietary fish oil replacement in gilthead sea bream(Sparus aurata L.)[J].Comparative Biochemistry and Physiology Part A:Molecular＆Integrative Physiology，2011，158(4):391－399.

[96] Zakeri M，Kochanian P，Marammazi J G，et al.Effects of dietary n－3 HUFA concentrations on spawning performance and fatty acids composition of broodstock，eggs and larvae in yellowfin sea bream，Acanthopagrus latus[J].Aquaculture，2011，310(3－4):388－394.

[97] Niculescu M D，Lupu D S，Craciunescu C N.Perinatal manipulation of α－linolenic acid intake induces epigenetic changes in maternal and offspring livers[J].The FASEB Journal，2013，27(1):350－358.

[98] Benatti P，Peluso G，Nicolai R，et al.Polyunsaturated fatty acids:biochemical，nutritional and epigenetic properties[J].Journal of the American College of Nutrition，2004，23(4):281－302.

[99] Seierstad S L，Svindland A，Larsen S，et al.Development of intimal thickening of coronary arteries over the lifetime of Atlantic salmon，Salmo salar L.，fed different lipid sources[J].Journal of Fish Diseases，2008，31(6):401－413.

[100] Xing Jinyi，Kang Li，Jiang Yunliang.Effect of dietary betaine supplementation on lipogenesis gene expression and CpG methylation of lipoprotein lipase gene in broilers[J].Molecular Biology Reports，2011，38(3):1975－1981.

[101] Xing Jinyi，Jiang Yunliang.Effect of dietary betaine supplementation on mRNA level of lipogenesis genes and on promoter CpG methylation of fatty acid synthase(FAS)gene in laying hens[J].African Journal of Biotechnology，2012，11(24):6633－6640.

[102] McGhee K E，Bell A M.Paternal care in a fish:epigenetics and fitness enhancing effects on offspring anxiety[J].Proceedings of the Royal Society B:Biological Sciences，2014，281(1794):20141146.

[103] Mizukami T，Kanai Y，F(xiàn)ujisawa M，et al.Five azacytidine，a DNA methyltransferase inhibitor，specifically inhibits testicular cord formation and Sertoli cell differentiation in vitro[J].Molecular Reproduction ＆ Development，2008，75(6):1002－1010.

[104] Matsumoto Y，Buemio A，Chu R，et al.Epigenetic control of gonadal aromatase(cyp19a1)in temperature－dependent sex determination of red－eared slider turtles[J].PLoS One，2013，8(6):e63599.

[105] Merchant－Larios H，Díaz－Hernández V.Environmental sex determination mechanisms in reptiles[J].Sexual Development，2013，7(1－3):95－103.

[106] Xiao Jun，Song Can，Liu Shaojun，et al.DNA methylation analysis of allotetraploid hybrids of red crucian carp(Carassius auratus red var.)and common carp(Cyprinus carpio L.)[J].PLoS One，2013，8(2):e56409.

[107] 王智誠(chéng)，周曉旭，王昊澤，等.長(zhǎng)鏈非編碼RNA及其在斑馬魚中的研究進(jìn)展[J].大連海洋大學(xué)學(xué)報(bào)，2017，32(2):248－254.

[108] Sadhu N，Sharma S R K，Joseph S，et al.Chronic stress due to high stocking density in open sea cage farming induces variation in biochemical and immunological functions in Asian seabass(Lates calcarifer，Bloch)[J].Fish Physiology and Biochemistry，2014，40(4):1105－1113.

[109] Xie Jun，Liu Bo，Zhou Qunlan，et al.Effects of anthraquinone extract from rhubarb Rheum officinale Bail on the crowding stress response and growth of common carp Cyprinus carpio var.Jian[J].Aquaculture，2008，281(1－4):5－11.

[110] Costas B，Arag?o C，Mancera J M，et al.High stocking density induces crowding stress and affects amino acid metabolism in Senegalese sole Solea senegalensis(Kaup 1858)juveniles[J].Aquaculture Research，2008，39(1):1－9.

[111] Caipang C M A，Brinchmann M F，Berg I，et al.Changes in selected stress and immune－related genes in Atlantic cod，Gadus morhua， following overcrowding[J].Aquaculture Research，2008，39(14):1533－1540.

[112] Gornati R，Papis E，Rimoldi S，et al.Rearing density influences the expression of stress－related genes in sea bass(Dicentrarchus labrax，L.)[J].Gene，2004，341:111－118.

[113] Aksakal E，Ekinci D，Erdoˇgan O，et al.Increasing stocking density causes inhibition of metabolic－antioxidant enzymes and elevates mRNA levels of heat shock protein 70 in rainbow trout[J].Livestock Science，2011，141(1):69－75.

[114] Salas－Leiton E，Anguis V，Martín－Antonio B，et al.Effects of stocking density and feed ration on growth and gene expression in the Senegalese sole(Solea senegalensis):potential effects on the immune response[J].Fish ＆ Shellfish Immunology，2010，28(2):296－302.

[115] Bouligand J，Ghervan C，Tello J A，et al.Isolated familial hypogonadotropic hypogonadism and a GNRH1 mutation[J].New England Journal of Medicine，2009，360(26):2742－2748.

[116] Alvarado S G，Lenkov K，Williams B，et al.Social crowding during development causes changes in GnRH1 DNA methylation[J].PLoS One，2015，10(10):e0142043.

[117] Rebl A，Zebunke M，Borchel A，et al.Microarray－predicted marker genes and molecular pathways indicating crowding stress in rainbow trout(Oncorhynchus mykiss)[J].Aquaculture，2017，473:355－365.

[118] Zhang Tieyuan，Hellstrom I C，Bagot R C，et al.Maternal care and DNA methylation of a glutamic acid decarboxylase 1 promoter in rat hippocampus[J].Journal of Neuroscience，2010，30(39):13130－13137.