Comparative transcriptome analysis in livers of female and male marine medaka (Oryzias melastigma)
-
摘要: 海水青鳉 (Oryzias melastigma) 作为雌雄异体的模式动物，在研究外来污染物毒性效应性别差异上具有优势。运用转录组学技术系统研究了雌、雄海水青鳉肝脏组织中的差异表达基因，结果显示雌、雄青鳉肝脏中共有1 351个显著差异表达基因，其中683个在雌鱼肝脏中高表达，668个在雄鱼肝脏中高表达。雌鱼肝脏中高表达的差异基因主要涉及生殖和性激素合成相关通路，如卵黄蛋白原和雌激素受体。雄鱼肝脏中高表达的差异基因主要涉及能量代谢、细胞骨架和肌肉收缩等相关生物过程，如丙酮酸激酶、肌酸激酶、肌球蛋白和肌钙蛋白等。实时荧光定量PCR验证结果显示，除DNA错配修复蛋白基因以外，其他17个基因差异表达倍数与RNA-seq对应基因差异表达趋势基本一致，表明转录组分析数据基本可靠。研究表明，雌、雄海水青鳉肝脏中基因表达具有差异调控模式，研究获得的差异基因和调控通路将为海水青鳉对外来污染物性别差异响应分子机制研究提供一定的理论基础。Abstract: As a gonochoristic model animal, marine medaka (Oryzias melastigma) is good for studying the sex-specific responses of organisms to xenobiotic pollutants. We used comparative transcriptomics technology to systematically investigate the differentially expressed genes (DEGs) between the liver tissues of female and male medaka. We identified 683 significantly up-regulated DEGs in the females, and 668 DEGs in the males. The high expressed DEGs in the females were involved in the reproductive and sex hormone synthesis pathways, such as vitellogenin and estrogen receptor. The top twenty DEGs in the males were involved in energy metabolism, cytoskeleton and muscle contraction, such as pyruvate kinase, creatine kinase, myosin and troponin. Except for the DNA mismatch repair protein, all the 17 DEGs had similar magnitude and expression trends by both qRT-PCR and RNA-seq analyses, which confirms the reliability of the RNA-seq data. The results demonstrate that the gene expression patterns are different in the livers of female and male medaka, and the DEGs provide a theoretical basis for promoting the molecular mechanism of sex-specific responses of medaka to xenobiotic pollutants.
-
Keywords:
- Oryzias melastigma /
- Liver /
- Sex /
- Transcriptome /
- Differential expression genes
-
鲟是鲟形目鱼类的总称，为世界上最古老的鱼类之一[1]，具有较高的经济价值[2]。西伯利亚鲟 (Acipenser baerii) 和施氏鲟 (A. schrenckii) 均属于鲟形目、鲟科、鲟属。西伯利亚鲟俗称“尖吻鲟”，在国内仅分布于新疆额尔齐斯河流域[3]；施氏鲟俗称“七粒浮子”，是我国最早开发利用的鲟类之一[4]，具有生长速度快、抗病力强等优点。目前，养殖纯种鲟普遍存在生长速度慢、病害频发、鱼子酱品质较低等缺点[5]。研究表明，鲟属之间的杂交所获得的杂交种在生长和繁殖方面通常显示出一定的优势[6]。因此，鲟鱼杂交种越来越受到生产者的青睐。本试验采用的杂交鲟 (西伯利亚鲟♀×施氏鲟♂) 具有生长速度快、养殖病害少、适温范围广、运输应激小等优点[7]，近年来已成为我国主要的人工养殖淡水鲟类，市场前景十分广阔。
研究鱼类早期发育对鱼类人工繁殖和苗种培育具有重要意义[8]，并且是构建鱼类人工繁育体系的基础。掌握鱼类早期发育的特点可以有效提高苗种成活率，如徐秋艳[9]在养殖西杂鲟时，分别从亲鱼选择与培育、人工催产与授精、苗种培育以及亲鱼产后护理等方面阐述了提高苗种成活率的具体措施，对于扩大种群资源具有重要意义[10]。
目前，众多学者已对西伯利亚鲟[11]、中华鲟 (A. sinensis)[12]、施氏鲟[13]、匙吻鲟 (Polyodon spathula)[14]、小体鲟 (A. ruthenus)[15]、施杂鲟 (施氏鲟♀×西伯利亚鲟♂)[16] 等多种鲟鱼的胚胎发育进行了研究，但对于杂交鲟 (西伯利亚鲟♀×施氏鲟♂) 的胚胎发育过程以及早期仔鱼发育情况，尚缺乏较为详尽的研究报道。本研究对杂交鲟胚胎发育及早期仔鱼生长进行了系统观察，并详细记录了各阶段的形态特征，旨在丰富杂交鲟繁殖生物学基础研究，为进一步发掘鲟鱼的繁殖性能提供基础数据和技术支撑，以期推动鲟鱼养殖业的可持续发展。
1. 材料与方法
1.1 亲鱼来源
试验所用亲鱼选自四川润兆渔业有限公司。其中母本西伯利亚鲟80尾，10龄以上，体质量为 (47.55±5.85) kg，全长为 (130.00±7.67) cm；父本施氏鲟11尾，8龄以上，体质量为 (12.35±0.95) kg，全长为 (103.00±2.00) cm。所选亲鱼体质健壮，外表光滑无伤，生长发育状况良好。人工催产期间水温为 (16.60±0.60) ℃。
1.2 方法
1.2.1 亲鱼选择与成熟度鉴定
生殖期的雌鱼体表黏液较多，生殖孔红肿外翻 (图1-a)，腹部膨大且按压会回弹；而生殖期的雄鱼从外观上没有明显的特征 (图1-b)，外表难以鉴别。目前普遍采用穿刺法判断卵细胞的成熟度：雌鱼成熟卵卵径 >2.8 mm，灰黑色、圆粒形、有弹性且两级分化明显，动物极一端出现色素环，中央有白色极斑 (图1-c)。卵母细胞的成熟度主要根据极化指数 (Polarization Index，PI值) 判断，PI值是指卵核到动物极卵膜内侧 (即图1-d中AB) 与卵长径 (动物极到植物极的卵径，即图1-d中AC) 的比值，PI值越小，表示卵细胞的成熟度越好[17]。将卵子用开水煮5 min，固定后用刀片沿着动物极中心轴切开，正放于载玻片上，在体视显微镜 (Olympus, SZX16) 下观察并测量细胞核的极化程度 (图1-d)，选择PI值 <0.1的雌鱼作为催产亲鱼。雄鱼成熟度的鉴别同样采用穿刺法，以采集到白色块状物判断其成熟度。
![]()
图 1 亲鱼雌雄与成熟度鉴定注：a为生殖期的雌鱼生殖孔；b为生殖期的雄鱼生殖孔；c为成熟卵细胞 (动物极观)；d中AC为卵长径，AB为卵核到动物极卵膜内侧长度，PI=AB/AC。Figure 1. Identification of male and female parent fish and maturityNote: a. The genital orifice of a female fish; b. The genital orifice of a male fish; c. The mature egg cell (animal view) ; in figure d, AC represents the egg length diameter, while AB represents the length from the egg nucleus to the inner side of the animal's polar oval membrane, PI=AB/AC.1.2.2 催产及人工授精
采用促黄体素释放激素A2 (LHRH-A2, 4~10 μg·kg−1) 和马来酸地欧酮 (DOM, 1~5 mg·kg−1) 的混合药物作为催产剂，所有雌鱼在注射时使用胸鳍基部2次注射法，第1次注射全部剂量的10%~20%，第2次注射剩余剂量；雄鱼采用1次注射法，在雌鱼注射第2针时同时注射，注射剂量为雌鱼的一半，注射部位与雌鱼注射部位一致。雄鱼的效应时间为6~8 h，雌鱼效应的时间为8~20 h。
挤出精、卵后采用干法授精，每4万~7万粒卵加入1~2 mL精液，加入清水激活精子，搅拌均匀后静置5 min，漂洗2~3次以去除黏液和多余的杂质。再利用滑石粉持续脱黏1 h左右后放入孵化筛中进行孵化[18]，孵化期间水温为16.3~17.1 ℃，平均溶解氧 (DO) 质量浓度为9.0 mg·L−1，pH 为7~8，每天早晚用塑料吸管挑出死卵。
1.2.3 胚胎发育和仔鱼生长观察
胚胎发育观察：每次随机选取同批次的受精卵30粒，在体视显微镜下观察胚胎发育各个时期的形态特征变化，并记录发育时间、温度和发育特征。受精卵阶段每隔10~20 min取样观察，卵裂阶段和囊胚阶段每隔1~2 h取样观察，原肠胚阶段到出膜阶段每隔2~4 h取样观察。参考已有的方法[13-14,19-20] 并根据显微镜下观察到的特征，对胚胎发育过程进行分期，以50%以上个体出现同一新特征时作为发育到一个新时期的起始时间。
早期仔鱼观察：以同批次孵化出膜的仔鱼为观察对象，每日随机选取10尾，观察记录仔鱼形态变化和器官发育情况。按照公式计算卵黄囊体积：卵黄囊体积=卵黄囊长径×(卵黄囊短径)2×π/6[21]。
1.3 繁殖效率评价指标
催产率、受精率和孵化率的计算公式分别为：催产率=产卵鱼数/催情鱼总数×100%，受精率=原肠中期活卵数/产卵总数×100%，孵化率=水花鱼苗/受精卵总数×100%。
数据经Excel 2021、Graphpad Prism 9.1和Adobe Photoshop CC 2023等软件整理并作图，结果以“平均值±标准误 ($\overline { x}\pm s_{\overline { x}} $)”表示。
2. 结果
2.1 人工繁殖情况
2022年3—4月，共催产11批杂交鲟 (母本80尾，父本11尾)，平均催产率为100%，受精率为 (80.08±8.57)%，孵化率为 (63.69±10.71)%。其中，4月30日催产的第11批受精率最高，达到90.33%，孵化率最高的是3月30日催产的第8批，达到77.53%，3月31日催产的第9批受精率和孵化率最低，分别为65.58%和26.47% (表1)。
表 1 2022年3—4月杂交鲟 (西伯利亚鲟♀×施氏鲟♂) 人工繁殖情况Table 1. Artificial breeding of hybrid sturgeon (A. baerii ♀×A. schrenckii ♂) during March–April 2022批次
Batch催产日期
Spawning date雌鱼
Number of female/尾雄鱼
Number of male/尾催产水温
Spawning water temperature/℃催产率
Spawning rate/%第1批 1st batch 2022-03-01 12 1 16.00 100.00 第2批 2nd batch 2022-03-08 6 1 16.20 100.00 第3批 3rd batch 2022-03-09 8 1 16.10 100.00 第4批 4th batch 2022-03-11 9 1 16.20 100.00 第5批 5th batch 2022-03-12 9 1 16.70 100.00 第6批 6th batch 2022-03-23 7 1 16.90 100.00 第7批 7th batch 2022-03-24 10 1 16.70 100.00 第8批 8th batch 2022-03-30 5 1 16.60 100.00 第9批 9th batch 2022-03-31 7 1 16.80 100.00 第10批 10th batch 2022-04-20 4 1 17.00 100.00 第11批 11th batch 2022-04-30 3 1 17.20 100.00 平均值 Mean 7 1 16.58 100.00 批次
Batch产卵数
Number of eggs/
万粒受精数
Number of
fertilized eggs/万粒受精率
Fertilization
rate/%出苗数
Number of
fry hatched/万粒孵化率
Hatching
rate/%第1批 1st batch 176.00 140.80 80.00 78.00 55.40 第2批 2nd batch 90.80 61.49 67.72 33.00 53.67 第3批 3rd batch 90.00 81.17 90.19 58.20 71.70 第4批 4th batch 96.00 66.22 68.98 26.85 50.54 第5批 5th batch 118.05 99.28 84.10 73.13 73.66 第6批 6th batch 109.20 97.90 89.65 60.00 61.29 第7批 7th batch 128.00 103.56 80.91 75.70 73.09 第8批 8th batch 112.80 93.04 82.48 72.13 77.53 第9批 9th batch 118.80 77.91 65.58 20.62 46.47 第10批 10th batch 77.60 62.84 80.98 11.92 70.84 第11批 11th batch 39.60 35.77 90.33 24.60 66.41 平均值 Mean 105.20 83.63 80.08 48.56 63.69 2.2 胚胎发育时序及特征
参照其他工作者的研究结果，根据实际观察到的胚胎发育顺序及形态特点，将杂交鲟胚胎发育分为7个阶段，共34个时期 (表2、图2)。在水温 16.3~17.1 ℃的条件下，经124~134 h孵化出膜，发育总积温为 2 083~2 251 h·℃。
表 2 杂交鲟 (西伯利亚鲟♀×施氏鲟♂) 胚胎发育Table 2. Embryonic development of hybrid sturgeon (A. baerii ♀×A. schrenckii ♂)胚胎发育阶段
Embryonic
development发育时期
Development
stage受精后时间
Time after
fertilization发育时间
Development
time水温
Water
temperature/℃积温
Accumulative
temperature/(h·℃)受精卵阶段 Fertilized egg 两极转动期 Pole conversion period 0 0 17.00 0 卵周隙形成期 Perivitelline space formation 30 min 30 min 17.10 8.55 胚盘隆起期 Blastodisc stage 1 h 30 min 1 h 17.00 17.00 卵裂阶段 Cleavage stage 2细胞期 2-cell stage 3 h 30 min 2 h 17.00 34.00 4细胞期 4-cell stage 4 h 30 min 1 h 17.10 17.10 8细胞期 8-cell stage 5 h 30 min 1 h 16.90 16.90 16细胞期 16-cell stage 6 h 30 min 1 h 16.90 16.90 32细胞期 32-cell stage 7 h 30 min 1 h 16.90 16.90 多细胞期 Multicellular stage 9 h 30 min 2 h 16.80 33.60 囊胚阶段 Blastula stage 囊胚初期 Initial blastula stage 11 h 30 min 2 h 16.70 33.40 囊胚中期 Mid-blastula stage 12 h 30 min 1 h 16.70 16.70 囊胚晚期 Late blastula stage 14 h 30 min 2 h 16.70 33.40 原肠胚阶段 Gastrula stage 原肠初期 Initial gastrula stage 18 h 3 h 30 min 16.70 58.45 原肠早期 Early gastrula stage 22 h 4 h 16.90 67.60 原肠中期 Midgastrula stage 26 h 4 h 16.70 66.80 大卵黄栓期 Big yolk plug stage 29 h 3 h 16.90 50.70 小卵黄栓期 Small yolk plug stage 31 h 2 h 16.60 33.20 隙状胚孔期 Lyriform blastopore formation 36 h 30 min 5 h 30 min 16.30 89.65 神经胚阶段 Neural stage 神经胚早期 Warly neurula stage 43 h 6 h 30 min 16.50 107.25 宽神经板期 Wide neural plate formation 44 h 30 min 1 h 30 min 16.50 24.75 神经褶靠拢期 Neural fold closing stage 49 h 4 h 30 min 16.90 76.05 神经胚晚期 Late neurula stage 51 h 30 min 2 h 30 min 16.80 42.00 神经管闭合期 Neural tube closing stage 52 h 30 min 1 h 16.80 16.80 器官形成阶段 Organogenesis stage 眼囊形成期 Eye sac formation 54 h 1 h 30 min 16.70 25.05 尾芽形成期 Caual bud appearance 59 h 30 min 5 h 30 min 16.70 91.85 尾芽分离期 Caual bud separating stage 66 h 30 min 7 h 16.70 116.90 短管心脏期 Shorttubular heart formation 69 h 30 min 3 h 16.40 49.20 长管心脏期 Longtubular heart formation 71 h 30 min 2 h 16.50 33.00 听板形成期 Otic placode appearance 76 h 4 h 30 min 16.50 74.25 肌肉效应期 Stage of muscular effect 79 h 3 h 16.40 49.20 心跳期 Stage of heart beating 95 h 30 min 16 h 30 min 16.40 270.60 尾达头部期 The tail touches the head 104 h 30 min 9 h 17.00 153.00 出膜阶段 Hatching stage 出膜前期 Hatching prophase 118 h 13 h 30 min 16.80 226.80 出膜期 Hatching stage 124~134 h 6 ~16 h 16.80 100.8~268.8 ![]()
图 2 杂交鲟 (西伯利亚鲟♀×施氏鲟♂) 胚胎发育过程注：A-1、A-2. 两极转动期；B. 卵周隙形成期；C. 胚盘隆起期；D、E、F. 2细胞期、4细胞期、8细胞期；G-1、G-2. 16细胞期 (动物极观)、16细胞期 (植物极观)；H-1、H-2. 32细胞期 (动物极观)、32细胞期 (植物极观)；I. 多细胞期；J. 囊胚初期；L. 囊胚中期；M. 囊胚晚期；N. 原肠初期；O. 原肠早期；P. 原肠中期；Q. 大卵黄栓期；K. 小卵黄栓期；S. 隙状胚孔期；T-1、T-2. 神经胚早期；U-1、U-2. 宽神经板期；V. 神经褶靠拢期；W. 神经胚晚期；X. 神经管闭合期；Y. 眼囊形成期；Z. 尾芽形成期；AA. 尾芽分离期；AB-1、AB-2、AB-3. 短管心脏期；AC-1、AC-2. 长管心脏期；AD-1、AD-2. 听板形成期；AE-1、AE-2. 肌肉效应期；AF. 心跳期；AG. 尾达头部期；AH. 出膜前期；AI. 出膜期；AJ-1. 尾部先出膜方式；AJ-2. 头部先出膜方式。Figure 2. Embryonic development of hybrid sturgeon (A. baerii ♀×A. schrenckii ♂)Note: A-1, A-2. Bipolar rotation stage; B. Perioval gap formation stage; C. Placental uplift; D, E, F. 2-cell stage, 4-cell stage, 8-cell stage; G-1,G-2. 16 cell stage (Animal polar stage), 16-cell stage (Plant polar stage); H-1, H-2. 32 cell stage (Animal polar stage), 32 cell stage (Plant polar stage); I. multicellular stage; J. early blastocyst stage; L. mid-blastocyst stage; M. late blastocyst stage; N. Primary gastrul stage; O. Early gastrul stage; P. Mid-gastrul stage; Q. Large yolk suppository stage; K. Small yolk suppository stage; S. Gap embryonic stage; T-1, T-2. Early neuroembryonic stage; U-1, U-2. Wide neural plate stage; V. Neural fold convergence stage; W. Late neuroembryonic stage; X. Neural tube closure stage; Y. Eye capsule formation stage; Z. Tail bud formation stage; AA. tail bud isolation stage; AB-1, AB-2, AB-3. short tube cardiac stage; AC-1, AC-2. long tube heart stage; AD-1, AD-2. Listening plate formation stage; AE-1, AE-2. Muscular effector phase; AF. Heartbeat phase; AG. Tail reaching head stage; AH. Premembrane exit; AI. Membrane release phase; AJ-1. Tail first membrane exit method; AJ-2. Head first membrane exit method.2.2.1 受精卵阶段
杂交鲟的受精卵为沉性卵，卵膜3层，卵径 (2.50±0.56) mm，灰黑色，不透明，卵黄丰富。受精卵的动物极始终朝上，植物极朝下，动物极中央出现明亮的极性斑 (图2-A-1、2-A-2)；受精后30 min，中央极性斑消失，卵周隙扩大 (图2-B)。受精后1 h，原生质向动物极集中，隆起形成胚盘 (图2-C)。
2.2.2 卵裂阶段
受精后3 h 30 min，胚盘开始分裂，动物极被分裂形成2个大小相等的细胞 (图2-D)；1 h后，再次分裂为4个等大的分裂球 (图2-E)，前两次卵裂的分裂面互相垂直，受精后5 h 30 min~7 h 30 min，依次进行第3、第4、第5次卵裂 (图2-F—2-H-2)，最终形成32个排列不规则且大小不一的细胞，每次形成的分裂面都与上一次的分裂面垂直，细胞间的轮廓比较清晰。2 h后动物极被分裂成多个不规则细胞，并出现色素，植物极被完全分裂，为多细胞期 (图2-I)。
2.2.3 囊胚阶段
随着卵裂的继续进行，分裂球变得更小。受精后11 h 30 min，动物极中央色素深、外沿色素浅，细胞变小，且细胞核分裂同步，细胞界限清晰 (图2-J)；1 h后，动物极细胞继续变小，细胞间隙逐渐明显 (图2-L)；2 h后，细胞越来越小且界限模糊，形成囊胚腔 (图2-M)。
2.2.4 原肠胚阶段
囊胚层细胞继续沿卵黄囊外周向植物极端下包，受精后18 h，赤道处出现一条深色素带 (图2-N)，为原肠初期。受精后22 h，色素带下延超过赤道产生一个短而浅的狭缝状胚孔 (图2-O)；4 h后，外胚层覆盖胚胎表面的2/3 (图2-P)，为原肠中期；受精后29 h，植物极的卵黄栓像一个大塞子嵌在胚环内 (图2-Q)，为大卵黄栓期；2 h后，植物极的卵黄栓似小塞子嵌在胚环内 (图2-R)，为小卵黄栓期。
2.2.5 神经胚阶段
受精后36 h 30 min，卵黄栓消失，胚孔靠近呈现隙状 (图2-S)；受精后43 h，可见背部形成神经板，隙状胚孔处形成神经沟，脑部出现神经褶 (图2-T、2-U)，进入神经胚早期；1 h 30 min后，神经板周围可以看到“马蹄”形神经褶，神经板中央出现神经沟 (图2-U-1、2-U-2)，为宽神经板期；受精后49 h，神经褶靠拢，神经板下陷变窄，出现排泄系统原基 (图2-V)，即为神经褶靠拢期；受精后51 h 30 min，头部神经褶靠近，躯干部出现排泄系统原基 (图2-W)，即为神经胚晚期；1 h后，神经板闭合形成神经管，头部分化出前、中、后3个脑泡，头和尾部明显，躯干部可见直线肌节 (图2-X)，为神经管闭合期。
2.2.6 器官形成阶段
受精后54 h，前脑出现2个眼突，即眼原基，中脑出现第1对咽弧原基 (图2-Y)，为眼囊形成期；受精后59 h 30 min，脑室稍隆起，侧板达到头部前端，排泄系统分化出肾小管原基，尾部扁平 (图2-Z)，为尾芽形成期；受精后66 h 30 min，侧板融合一侧出现心脏原基，头部可见3对月牙型咽弧原基，尾芽突出呈棒状，尾尖略游离卵黄囊 (图2-AA)，为尾芽分离期；受精后69 h 30 min，心脏呈短管状，较短，略超过头部 (图2-AB-1、2-AB-2、2-AB-3)，为短管心脏期；2 h后心脏呈长管状，略弯曲，视泡明显 (图2-AC-1、2-AC-2)，为长管心脏期；受精后76 h，心脏呈小“C”型，椭圆形听板出现 (图2-AD-1、2-AD-2)，称为听板形成期；受精后79 h，头尾开始摆动，肌肉颤动，尾的末端接近心脏，有微弱搏动，出现嗅板和听泡 (图2-AE-1、2-AE-2)，为肌肉效应期；受精后95 h，心脏呈大“C”形，心脏搏动有节律性，血液呈乳白色，肌肉能收缩，尾端接触心脏，尾端游离胚胎1/3，可大幅摆动 (图2-AF)，称为心跳期；9 h后尾端接触头部，鳍褶变宽，出现肛门原基，胚胎在膜内经常转动，出现黑色素颗粒，眼呈新月形 (图2-AG)，为尾达头部期。
2.2.7 出膜阶段
受精后118 h，尾端略过头部，卵黄囊上血管明显，心跳加快，胚胎剧烈扭动，此时卵膜极易破裂 (图2-AH)，为出膜前期；受精后124~134 h，尾的末端达到间脑，鳃盖形成，血液循环到达尾部，仔鱼陆续出膜 (图2-AI)，为出膜期。出膜方式有尾部先出膜 (图2-AJ-1) 和头部先出膜 (图2-AJ-2) 2种。
2.3 早期仔鱼生长发育
杂交鲟早期仔鱼发育可分为3个阶段[18]：内源营养期，即从仔鱼刚出膜 (0 d) 到仔鱼开口摄食 (7 d)；混合营养期，即从具有开口摄食能力 (7 d) 到卵黄囊完全被吸收 (10 d)；10 d后完全摄食外源营养，进入外源营养期。
如图3所示，初孵仔鱼全长 (10.61±0.58) mm，卵黄囊约占全长的2/3。头部比较小，紧贴于卵黄囊。卵黄囊分为前后两部分，前半部分为圆球状，后半部分为短棒状，此时眼部黑色素颗粒较少。1日龄仔鱼的头部上抬与身体保持水平，口腔开始分化，卵黄囊体积逐渐变小，眼部黑色素颗粒增多；2日龄仔鱼的口腔和消化系统进一步分化，卵黄囊体积明显变小，眼部黑色素明显增多，体色素增多，可见肠道以及肠内出现少量黑色素栓；3日龄仔鱼的卵黄囊体积持续变小，尾部出现黑色素颗粒，游泳能力逐渐增强，但无法控制平游，只能做上下垂直运动；4日龄仔鱼的卵黄囊呈条状，全身分布不规则的黑色素颗粒，可见清晰的红色血液循环路径；5日龄仔鱼的头部变宽，脑部变大，背鳍和尾鳍开始分开，胸鳍出现；6日龄仔鱼的背鳍和尾鳍进一步分开，体色加深，体内黑色栓排出，池底可见大量黑色素栓，大部分仔鱼聚集在池底或沿壁活动，此时能平稳水平游动；7日龄仔鱼的头部前伸，背鳍和尾鳍完全分开，脊索向上弯曲，卵黄囊仅剩一小部分，此时仔鱼全长增长到 (17.17±0.64) mm，且可以自主摄取外界食物，但效果较差，容易死亡。
![]()
图 3 杂交鲟 (西伯利亚鲟♀×施氏鲟♂) 卵黄囊期仔鱼发育特征注：0—7分别表示出膜日龄；拍摄方位：头朝前，左侧观。Figure 3. Developmental characteristics of hybrid sturgeon (A. baerii ♀×A. schrenckii ♂) larvae at yolk sac stageNote: 0–7 represent age at film emergence respectively; Location: Head forward, left view.如图4所示，杂交鲟卵黄囊期仔鱼从初孵到卵黄完全吸收共历时10 d，卵黄囊吸收曲线V=12.53e−0.562t (r2=0.967 4)，呈指数衰减型函数。仔鱼消耗卵黄的能力随生长发育阶段不断加强，在第7天开口并进入混合营养期，到第10天时卵黄被完全吸收，进入外源营养期，混合营养期仅持续3 d。
![]()
图 4 杂交鲟 (西伯利亚鲟♀×施氏鲟♂) 卵黄囊期仔鱼卵黄吸收过程Figure 4. Process of yolk sac absorption in hybrid sturgeon (A. baerii ♀×A. schrenckii ♂)通过拟合生长曲线，杂交鲟卵黄囊期仔鱼全长与日龄的关系符合指数增长型函数 (图5)，生长方程为L=10.964e0.065 6t (r2=0.978 9, 1≤t≤8)。杂交鲟卵黄囊期仔鱼体质量与日龄的生长方程为M= −0.000 2t2+0.002 9t+0.012 9 (r2=0.983 6, 1≤t≤8)。
![]()
图 5 杂交鲟 (西伯利亚鲟♀×施氏鲟♂) 卵黄囊期仔鱼全长、体质量与日龄的关系Figure 5. Relationship between total length, body mass and age of day of hybrid sturgeon (A. baerii ♀×A. schrenckii ♂) at yolk sac stage3. 讨论
3.1 影响杂交鲟人工繁殖效率的因素
研究鲟鱼人工繁殖技术，掌握影响繁殖效率的重要因素，对于解决我国鲟鱼苗种供应问题、拓展我国鲟鱼产业的深度、保护野生鲟鱼种质资源等方面具有重要意义[22]。人工繁殖的效率受诸多因素的影响，如催产药物、水温、DO以及亲鱼质量等[23]。目前鲟类人工繁殖常用的催产剂有LHRH-A2、鲤鱼脑垂体、DOM等[19,24]。
本试验人工催产杂交鲟所采用的试剂及方法与朱华等[25] 总结出的人工催产俄罗斯鲟 (A. gueldenstaedtii) 的最佳剂量与注射方式基本一致，而罗江等[26] 总结出人工催产中华鲟的剂量与本研究相比较低，这可能与亲鱼的体质量以及成熟度有关。DO是影响人工繁殖效率的重要因素之一，鲟鱼对DO要求较高，水中DO下降会阻碍其生长、繁殖等[9]。本试验中，杂交鲟孵化期间水中的平均DO质量浓度为9.00 mg·L−1，高于肖鑫鑫等[27] 在催产杂交鲟 (达氏鳇Huso dauricus ♀×施氏鲟A. schrenckii ♂) 时的水平，表明DO充足。除了这两个因素之外，鱼类在产卵期间对水温也极为敏感，对于冷水鱼，水温突然升高会导致产卵停止或流产[28]，进而影响繁殖。本试验催产及孵化期间的水温保持在 (16.60±0.60) ℃，符合鲟鱼的最佳繁殖条件。
本试验催产时间处于鲟鱼的最佳繁殖季节[29]，水温适宜，且DO充足，理论上其繁殖效率应该偏高。但结果显示杂交鲟的平均受精率为80.08%，平均孵化率为63.69%；且在温差不超过0.3 ℃、催产时间间隔1 d的情况下，2022年3月30日催产的第8批和2022年3月31日催产的第9批的受精率和孵化率相差较大。而肖鑫鑫等[27]在对杂交鲟进行人工繁殖的试验中，在水温 16.00~17.20 ℃的条件下，其平均受精率达到82.90%，平均孵化率达到70.44%，高于本试验平均值，并且所催产的五批杂交鲟中没有出现受精率和孵化率相差较大的情况。有研究指出，曾经产过卵的亲鱼，其性腺发育与未产过卵的亲鱼存在差异，若同时催产会导致受精率和孵化率较低[30]，因此本试验繁殖效率较低的原因可能在于部分亲鱼在筛选时已产过卵，但这需要进一步的研究验证。
3.2 杂交鲟与其他鲟类胚胎发育的异同
如表3所示，杂交鲟与母本西伯利亚鲟胚胎发育的过程基本相似，呈“偏母遗传”特性，这与施氏鲟♀×西伯利亚鲟♂[20]、达氏鳇♀×施氏鲟♂[27]、鲇Silurus asotus ♀×南方鲇S. meridionalis ♂[31] 等杂交种的胚胎发育与母本相似的研究结果一致。与西伯利亚鲟的不同之处表现在：杂交鲟的平均卵径较大[32]；出膜前期杂交鲟尾的末端到达延脑，而西伯利亚鲟尾的末端到达间脑靠后。整体来看杂交鲟比西伯利亚鲟发育更快，表现出杂种优势。
表 3 杂交鲟 (西伯利亚鲟♀×施氏鲟♂) 与其父母本、反交组合胚胎发育特征的比较Table 3. Comparison of embryonic developmental characteristics of hybrid sturgeon (A. baerii ♀×A. schrenckii ♂) with its parents and inbred combinations项目
Item本试验杂交鲟
A. baeri ♀×A. schrenckii ♂西伯利亚鲟[11]
A. baeri施氏鲟[13]
A. schrenckii反交组合[20]
A. schrenckii ♀×A. baeri ♂卵径
Egg diameter/mm3.0~3.3 2.9~3.2 3.2~3.8 3.4~4.5 原肠中期
Mid-gastrula stage动物极覆盖胚胎表面的2/3 动物极覆盖胚胎表面的2/3 动物极覆盖胚胎表面的3/4 动物极覆盖胚胎表面的2/3 心跳期
Stage of heart beating心脏呈“C”形，心脏有节律搏动时尾末端到达心脏 心脏呈“C”形，心脏有节律搏动时尾末端到达心脏 心脏呈“S”形，心脏有节律搏动时尾末端未到达心脏 心脏膨大呈“S”形，心脏有节律搏动时尾末端未到达心脏 出膜期
Hatching stage尾的末端达到延脑 尾的末端达到间脑 尾的末端达到或略过头部 尾的末端超过头部 初孵仔鱼规格
Size of first hatchlings/mm9.2~10.7 8.9~9.5 9.4~10.0 发育水温
Developmental water temperature/℃16.3~17.1 15.5~18.0 16.0~17.0 18.5 发育时间
Developmental time/h124~134 133~145 114~146 120~130 积温
Temperature accumulation/h·℃2 083~2 251 2 173~2 369 1 856~2 108 2 054~2 405 与父本施氏鲟相比：受精卵均为沉性，但发育到原肠中期时，杂交鲟胚胎的动物极覆盖表面2/3，而施氏鲟胚胎的动物极覆盖表面3/4；发育到心跳期时，杂交鲟心脏呈“C”形，施氏鲟心脏呈“S”形；心脏有节律搏动时杂交鲟尾末端到达心脏，而施氏鲟尾末端未到达心脏；相同条件下，施氏鲟胚胎发育所需的时间较短，总积温较低[13]；出膜前期，杂交鲟尾的末端到达延脑，而施氏鲟尾的末端仅达到或略过头部；刚孵出的杂交鲟比刚孵出的施氏鲟的个体稍大，这与郭长江[33]的研究结果一致：鲟科鱼类出膜前期，其胚胎尾部末端所到达的位置与出膜仔鱼的体长成正比。
与反交组合 (施氏鲟♀×西伯利亚鲟♂) 相比，杂交鲟卵径较小，其他不同之处与父本的不同一致[30]。
3.3 杂交鲟卵黄囊期仔鱼的混合营养期
鱼苗从具有开口摄食能力到完全依赖外源性营养的时期称为混合营养期[34]。在此阶段，仔鱼可以开口摄食，但卵黄囊并没有被完全吸收，主要完成口、眼、消化道、血液循环系统和鳍的初步发育，并且基本可以实现平游[35]。不同鱼类的混合营养期长短不同。研究表明，达氏鳇9日龄开口，第13天时卵黄囊完全消失，混合营养期持续4 d[33]；匙吻鲟卵黄囊期仔鱼在6日龄左右开口，10日龄左右卵黄完全消失，混合营养期持续4 d[34]；施杂鲟出膜第7天开口，第9天卵黄囊吸收完毕，混合营养期持续2 d[16]；黄条鰤 (Seriola aureovittata) 3日龄开口，10日龄卵黄和油球消耗完毕，混合营养期持续7 d[36]；黑龙江三角鲂 (Megalobrama terminalis) 出膜第4天开口，出膜第6天卵黄被完全吸收，混合营养期持续2 d[37]。本研究中杂交鲟在出膜第7天开口，第10天卵黄被完全吸收，混合营养期仅持续3 d，与施杂鲟的情况相似。有学者指出，当仔鱼的卵黄囊被吸收殆尽时，必须马上建立主动摄食能力以完成内源性营养向外源性营养的转换，来维持早期正常的新陈代谢和生长发育[38]。本研究中仔鱼混合营养期与具有开口摄食能力的时间基本一致，说明杂交鲟在消耗卵黄的同时也在建立摄食外源性营养的能力。研究表明仔鱼从开口摄食到饥饿不可逆点的时间越长，适应开口饵料的能力越强[39]。
4. 结论
综上所述，本试验杂交鲟的平均催产率为100%，受精率为 (80.08±8.57)%，孵化率为 (60.96±15.12)%。在水温16.30~17.10℃条件下，杂交鲟由受精卵到孵化出仔鱼需要124~134 h，发育总积温为2 083~2 251 h·℃。通过观察杂交鲟胚胎发育各阶段的特征变化，将杂交鲟胚胎发育过程分为7个阶段、34个时期；卵黄囊期仔鱼在出膜第7天开口，第10天卵黄被完全吸收，混合营养期持续3 d；将杂交鲟与其父母本以及反交组合的胚胎发育过程比较发现，杂交鲟的发育特征与母本西伯利亚鲟相似，具有“偏母遗传”特性。本研究系统记录并完善了杂交鲟胚胎发育及早期仔鱼生长各个阶段的变化特征，为其高质量育种提供了科学依据。
-
![]()
图 1 雌、雄海水青鳉肝脏基因差异表达分析火山图
Figure 1. Volcano plot of deferential expressed genes in livers of female and male O. melastigma
![]()
图 2 qRT-PCR验证结果
Figure 2. Comparison of relative expression levels between RNA-Seq and qRT-PCR results
表 1 靶基因qRT-PCR引物序列
Table 1 Primer sequences of tested genes used in qRT-PCR
Gene 编号
Gene ID注释结果
Gene description引物序列 (5'—3')
Primer sequence (5'–3')ENSOMEG00000014272 丙酮酸激酶 pyruvate kinase GATGCAGGTTCTTCCGTTAT
GTTTCAGCGTGGTATTCGTGENSOMEG00000023209 肌酸激酶 creatine kinase M-type TACAAGCCCACCGACAAGCA
GCCAGAGCCTCAATGGACAGENSOMEG00000008770 ATP依赖6-磷酸果糖激酶 ATP-dependent 6-phosphofructokinase, muscle type TCCTCCAGGGACATCAGA
CCAGGAAAGCCTCAAAGCENSOMEG00000013229 异柠檬酸脱氢酶 isocitrate dehydrogenase [NADP], mitochondrial ACATTCCTCGGCTTGTTC
ACATTCCTCGGCTTGTTCENSOMEG00000016412 葡萄糖激酶 glucokinase GCGATTTCTTGGCTTTGG
GGAGTACATTTGGTTCGTENSOMEG00000013940 肌球蛋白重链 myosin heavy chain AAGGCTAACAGTGAGGTGG
CTCCACATCAATCATAAGGTCENSOMEG00000022055 肌钙蛋白C troponin C ACTCACCCAAACGGACCCAG
GCCCAGCATCCTCATCACCTENSOMEG00000017631 卵黄蛋白原1 vitellogenin-1 ACCCTCTACTCTGTCAACGA
TATCTTCTGGCACTCCTCACENSOMEG00000017938 卵黄蛋白原2 vitellogenin-2 TGAAAGATGTACCGAGTGCG
TCAATGGGTGTTTGGAGGAGENSOMEG00000015512 3-β-羟基类固醇异构酶 3-beta-hydroxysteroid-Delta (8), Delta (7)-isomerase TGGAGTTCCGTGTTTGAGTT
GGTCATTAGCCACATAGTTAGGENSOMEG00000010141 17-β-羟基类固醇脱氢酶 17-beta-hydroxysteroid dehydrogenase 7 GTGCTGCCAAAGAAATCAAA
ACATAGAAGAGGCTCCAACTENSOMEG00000009991 雌激素受体 estrogen receptor ATCAGCCCAGCCTCCTCA
GGGTCCGTTCGTCTCCATENSOMEG00000011913 脂肪酸脱羧酶2 fatty acid desaturase 2 CGGCACTGCTGGCAACTT
GGCTGAACGGCTCCTAAAENSOMEG00000009600 细胞色素P450 cytochrome P4502K1 ACGAGCCAACGAGACAATAC
TAGCAGTCCACAAATCCCCTENSOMEG00000011275 DNA复制许可因子 DNA replication licensing factor MCM6 CCACGGAAACGACGAGGTA
TTTGGCGGTACTCGGGTCTENSOMEG00000001336 DNA错配修复蛋白 DNA mismatch repair protein GGCGTCGCATCATAGTAGCT
CCCTTCTTCCTCCTCCTCTTENSOMEG00000001928 DNA修复蛋白 DNA repair protein RAD51 homolog 1 CGAGTTTGGTGTTGCCGTAG
AGGTCAGATTTCCCAGCGTCENSOMEG00000014847 DNA聚合酶 DNA polymerase epsilon subunit 2 GTGCCCAGATACATTTACAACG
GTGACCATCCCGAGTACGATTA18S核糖体RNA 18S ribosomal RNA GACAAATCGCTCCACCAACT
CCTGCGGCTTAATTTGACCC
下载: 导出CSV
表 2 海水青鳉肝脏转录组测序结果统计
Table 2 Statistics of liver transcriptome of O. melastigma
样品名称
Sample原始数据
Raw reads净序列
Clean reads错误率
Errorrate/%质量值大于30
Q30/%比对到基因组上的序列
Total mapped reads to genome/%CF1 46 032 406 42 172 460 0.02 94.30 91.89 CF2 46 965 052 43 062 498 0.03 93.76 91.71 CF3 45 711 652 41 474 580 0.03 93.90 91.92 CM1 44 214 558 41 630 522 0.03 93.70 90.06 CM2 47 037 770 44 496 970 0.03 93.99 86.77 CM3 47 890 610 44 474 994 0.03 94.00 87.79
下载: 导出CSV
表 3 雌雄海水青鳉肝脏基因表达变化倍数最高的前20个基因
Table 3 Top twenty significantly up-regulated genes in livers of female and male O. melastigma
基因ID
Gene ID变化倍数
log2 fold change基因注释
Gene description雌鱼 Female ENSOMEG00000001115 10.20 溶质载体家族41成员 solute carrier family 41 member 1 ENSOMEG00000003787 10.05 WAP型 (乳清酸性蛋白质) 四二硫核心 WAP-type (Whey Acidic Protein) 'four-disulfide core' ENSOMEG00000009600 9.90 细胞色素 P450 cytochrome P450 2K1 novel.1314 9.90 卵黄蛋白原1 vitellogenin-1 ENSOMEG00000017631 9.80 卵黄蛋白原1 vitellogenin-1 ENSOMEG00000003805 9.69 WAP型 (乳清酸性蛋白质) 四二硫核心 WAP-type (Whey Acidic Protein) 'four-disulfide core' ENSOMEG00000014709 9.65 卵黄蛋白原1 vitellogenin-1 ENSOMEG00000017938 9.63 卵黄蛋白原2 vitellogenin-2 ENSOMEG00000016745 9.35 脑特异性血管生成抑制因子1 brain-specific angiogenesis inhibitor 1 ENSOMEG00000014656 9.31 卵黄蛋白原1 vitellogenin-1 ENSOMEG00000012102 9.20 蛋白赖氨酸6-氧化酶 protein-lysine 6-oxidase ENSOMEG00000016916 9.11 神经元正五肽受体 neuronal pentraxin receptor ENSOMEG00000013237 8.76 cornichon同源蛋白3 protein cornichon homolog 3 ENSOMEG00000017881 8.73 ras相关蛋白 ras-related protein Rab-39B ENSOMEG00000014740 8.73 卵黄蛋白原1 vitellogenin-1 ENSOMEG00000016835 8.67 细胞外丝氨酸/苏氨酸蛋白激酶 extracellular serine/threonine protein kinase FAM20C ENSOMEG00000016456 8.17 卵黄蛋白原1 vitellogenin-1 ENSOMEG00000017903 7.67 ras相关蛋白 ras-related protein Rab-38 ENSOMEG00000008961 6.89 透明带精子结合蛋白4 zona pellucida sperm-binding protein 4 ENSOMEG00000009991 6.45 雌激素受体 estrogen receptor 雄鱼 Male ENSOMEG00000014272 13.42 丙酮酸激酶 pyruvate kinase PKM ENSOMEG00000015762 11.86 肌酸激酶 creatine kinase M-type ENSOMEG00000018541 11.74 热休克蛋白71 heat shock cognate 71 kDa protein ENSOMEG00000013940 11.72 肌球蛋白重链 myosin heavy chain, fast skeletal muscle ENSOMEG00000021803 11.30 肌球蛋白调节轻链2 myosin regulatory light chain 2, skeletal muscle isoform ENSOMEG00000021829 11.00 磷酸甘油酸突变酶2 phosphoglycerate mutase 2 ENSOMEG00000005322 10.74 肌球蛋白重链 myosin heavy chain, fast skeletal muscle ENSOMEG00000015087 10.60 肥大细胞蛋白酶8 mast cell protease 8 ENSOMEG00000000307 10.56 肌钙蛋白C troponin C, skeletal muscle ENSOMEG00000012484 10.55 蛋白酪氨酸磷酸酶 protein tyrosine phosphatase type IVA 3 ENSOMEG00000007501 10.50 ADP/ATP转位酶 ADP/ATP translocase 1 ENSOMEG00000023209 10.48 肌酸激酶 creatine kinase M-type ENSOMEG00000008770 10.47 ATP依赖6-磷酸果糖激酶 ATP-dependent 6-phosphofructokinase, muscle type ENSOMEG00000002691 10.25 辛肌动蛋白结合重复序列蛋白 Xin actin-binding repeat-containing protein 2 ENSOMEG00000006520 10.24 小清蛋白 parvalbumin beta ENSOMEG00000008672 10.21 LIM结构域结合蛋白 LIM domain-binding protein 3 ENSOMEG00000018391 10.15 肌球蛋白结合蛋白C myosin-binding protein C, fast-type ENSOMEG00000020650 9.98 ryanodine受体 ryanodine receptor 1 ENSOMEG00000001483 9.96 结蛋白 desmin ENSOMEG00000015807 9.96 AMP脱氨酶 AMP deaminase 1
下载: 导出CSV
表 4 差异表达基因显著富集KEGG通路
Table 4 Significantly enriched KEGG pathways
KEGG 通路
KEGG pathway雌鱼
Female雄鱼
MaleP 通路 ID
Pathway IDDNA复制与修复过程 DNA replication and repair process DNA复制 DNA replication 20 0 1.91×10−11 ola03030 错配修复 Mismatch repair 11 0 5.03×10−6 ola03430 嘧啶代谢 Pyrimidine metabolism 12 3 3.31×10−3 ola00240 嘌呤代谢 Purine metabolism 10 13 2.52×10−2 ola00230 核苷酸切除修复 Nucleotide excision repair 9 0 3.18×10−2 ola03420 碱基切除修复 Base excision repair 7 1 4.95×10−2 ola03410 脂类代谢过程 Lipid metabolism process 甾体生物合成 Steroid biosynthesis 13 0 2.64×10−7 ola00100 PPAR信号通路 PPAR signaling pathway 8 8 6.39×10−3 ola03320 脂肪酸代谢 Fatty acid metabolism 12 2 1.70×10−2 ola01212 不饱和脂肪酸生物合成 Biosynthesis of unsaturated fatty acids 7 2 2.52×10−2 ola01040 脂肪酸生物合成 Fatty acid biosynthesis 6 0 3.78×10−2 ola00061 萜骨架生物合成 Terpenoid backbone biosynthesis 6 0 3.78×10−2 ola00900 氨基酸与蛋白生物合成 Amino acid and protein biosynthesis 真核生物核糖体生物发生 Ribosome biogenesis in eukaryotes 21 1 2.64×10−7 ola03008 内质网蛋白加工 Protein processing in endoplasmic reticulum 27 3 3.31×10−3 ola04141 精氨酸和脯氨酸代谢 Arginine and proline metabolism 4 8 1.22×10−2 ola00330 氨基酸生物合成 Biosynthesis of amino acids 9 7 1.30×10−2 ola01230 氨酰tRNA生物合成 Aminoacyl-tRNA biosynthesis 10 0 1.31×10−2 ola00970 蛋白质输出 Protein export 7 0 1.81×10−2 ola03060 其他 Others N-聚糖生物合成 N-Glycan biosynthesis 11 0 2.52×10−2 ola00510
下载: 导出CSV
-
[1] CHEN X P, LI L, CHENG J P, et al. Molecular staging of marine medaka: a model organism for marine ecotoxicity study[J]. Mar Pollut Bull, 2011, 63(5/6/7/8/9/10/11/12): 309-317.
[2] KIM B, KIM J, CHOI I Y, et al. Omics of the marine medaka (Oryzias melastigma) and its relevance to marine environmental research[J]. Mar Environ Res, 2016, 113: 141-152. doi: 10.1016/j.marenvres.2015.12.004
[3] YIN X H, LIU Y, ZEB R, et al. The intergenerational toxic effects on offspring of medaka fish Oryzias melastigma from parental benzo[a]pyrene exposure via interference of the circadian rhythm[J]. Environ Pollut, 2020, 267: 115437. doi: 10.1016/j.envpol.2020.115437
[4] LU W J, LONG L, ZHAO P Q, et al. Perfluorinated compounds disrupted osmoregulation in Oryzias melastigma during acclimation to hypoosmotic environment[J]. Ecotoxicol Environ Saf, 2021, 223: 112613. doi: 10.1016/j.ecoenv.2021.112613
[5] FONG C C, SHI Y F, YU W K, et al. iTRAQ-based proteomic profiling of the marine medaka (Oryzias melastigma) gonad exposed to BDE-47[J]. Mar Pollut Bull, 2014, 85(2): 471-478. doi: 10.1016/j.marpolbul.2014.04.024
[6] WANG R F, ZHU L M, ZHANG J, et al. Developmental toxicity of copper in marine medaka (Oryzias melastigma) embryos and larvae[J]. Chemosphere, 2020, 247: 125923. doi: 10.1016/j.chemosphere.2020.125923
[7] ZHANG Y B, WANG J, LU L, et al. Genotoxic biomarkers and histological changes in marine medaka (Oryzias melastigma) exposed to 17α-ethynylestradiol and 17β-trenbolone[J]. Mar Pollut Bull, 2020, 150: 110601. doi: 10.1016/j.marpolbul.2019.110601
[8] CONG Y, JIN F, WANG J Y, et al. The embryotoxicity of ZnO nanoparticles to marine medaka, Oryzias melastigma[J]. Aquat Toxicol, 2017, 185: 11-18. doi: 10.1016/j.aquatox.2017.01.006
[9] ZHANG Y T, CHEN M Y, HE S Q, et al. Microplastics decrease the toxicity of triphenyl phosphate (TPhP) in the marine medaka (Oryzias melastigma) larvae[J]. Sci Total Environ, 2020, 763: 143040.
[10] CONG Y, JIN F, TIAN M, et al. Ingestion, egestion and post-exposure effects of polystyrene microspheres on marine medaka (Oryzias melastigma)[J]. Chemosphere, 2019, 228: 93-100. doi: 10.1016/j.chemosphere.2019.04.098
[11] KANG H M, BYEON E, JEONG H, et al. Different effects of nano- and microplastics on oxidative status and gut microbiota in the marine medaka Oryzias melastigma[J]. J Hazard Mater, 2020, 405: 124207.
[12] 靳非, 田淼, 穆景利, 等. 聚苯乙烯微塑料长期暴露对海水青鳉 (Oryzias melastigma) 亲代生长、繁殖及子代发育的影响[J]. 生态毒理学报, 2021, 16(4): 216-223. [13] WANG J, LI Y J, LU L, et al. Polystyrene microplastics cause tissue damages, sex-specific reproductive disruption and transgenerational effects in marine medaka (Oryzias melastigma)[J]. Environ Pollut, 2019, 254: 113024. doi: 10.1016/j.envpol.2019.113024
[14] HE S W, YU D D, LI P, et al. A new perspective on endocrine disrupting effects of triphenyltin on marine medaka: from brain transcriptome, gut content metabolome and behavior[J]. Chemosphere, 2022, 307: 136190. doi: 10.1016/j.chemosphere.2022.136190
[15] LIANG P P, SAQIB H S A, NI X M, et al. Long-read sequencing and de novo genome assembly of marine medaka (Oryzias melastigma)[J]. BMC Genomics, 2020, 21(1): 640. doi: 10.1186/s12864-020-07042-7
[16] LAI K P, TAM N, WANG S Y, et al. Hypoxia causes sex-specific hepatic toxicity at the transcriptome level in marine medaka (Oryzias melastigma)[J]. Aquat Toxicol, 2020, 224: 105520. doi: 10.1016/j.aquatox.2020.105520
[17] LIANG P P, SAPIB H S A, LIN Z Y, et al. RNA-seq analyses of marine medaka (Oryzias melastigma) reveals salinity responsive transcriptomes in the gills and livers[J]. Aquat Toxicol, 2021, 240: 105970. doi: 10.1016/j.aquatox.2021.105970
[18] YE R R, LEI E N Y, LAM M H W, et al. Gender-specific modulation of immune system complement gene expression in marine medaka Oryzias melastigma following dietary exposure of BDE-47[J]. Environ Sci Pollut Res Int, 2011, 19(7): 2477-2487.
[19] YU W K, SHI Y F, FONG C C, et al. Gender-specific transcriptional profiling of marine medaka (Oryzias melastigma) liver upon BDE-47 exposure[J]. Comp Biochem Physiol D, 2013, 8(3): 255-262.
[20] YIN X H, ZEB R, WEI H, et al. Acute exposure of di (2-ethylhexyl) phthalate (DEHP) induces immune signal regulation and ferroptosis in Oryzias melastigma[J]. Chemosphere, 2021, 265: 129053. doi: 10.1016/j.chemosphere.2020.129053
[21] QIAO Q, MANACH S L, SOTTON B, et al. Deep sexual dimorphism in adult medaka fish liver highlighted by multi-omic approach[J]. Sci Rep, 2016, 26(6): 32459.
[22] ZHENG W L, XU H Y, LAM S H, et al. Transcriptomic analyses of sexual dimorphism of the zebrafish liver and the effect of sex hormones[J]. PLoS One, 2013, 8(1): e53562. doi: 10.1371/journal.pone.0053562
[23] MORTAZAVI A, WILLIAMS B A, MCCUE K, et al. Mapping and quantifying mammalian transcriptomes by RNA-Seq[J]. Nature methods, 2008, 5(7): 621-628. doi: 10.1038/nmeth.1226
[24] LOVE M I, HUBER W, ANDERS S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2[J]. Genome Biol, 2014, 15(12): 550. doi: 10.1186/s13059-014-0550-8
[25] YOUNG M D, WAKEFIELD M J, SMYTH G K, et al. Gene ontology analysis for RNA-seq: accounting for selection bias[J]. Genome Biol, 2010, 11(2): R14. doi: 10.1186/gb-2010-11-2-r14
[26] KANEHISA M, ARAKI M, GOTO S. KEGG for linking genomes to life and the environment[J]. Nucleic Acids Res, 2008, 36: 480-484.
[27] CHEN L G, ZHANG W P, YE R, et al. Chronic exposure of marine medaka (Oryzias melastigma) to 4, 5-Dichloro-2-n-octyl-4-isothiazolin-3-one (DCOIT) reveals its mechanism of action in endocrine disruption via the hypothalamus-pituitary-gonadal-liver (HPGL) axis[J]. Environ Sci Technol, 2016, 50(8): 4492-4501. doi: 10.1021/acs.est.6b01137
[28] LIVAK K J, SCHMITTGEN T D. Analysis of relative gene expression data using realtime quantitative PCR and the 2(-Delta Delta C(T)) method[J]. Methods, 2001, 25: 402-408. doi: 10.1006/meth.2001.1262
[29] WU L Y, CHEN H G, RU H Y, et al. Sex-specific effects of triphenyltin chloride (TPT) on thyroid disruption and metabolizing enzymes in adult zebrafish (Danio rerio)[J]. Toxicol Lett, 2020, 331: 143-151. doi: 10.1016/j.toxlet.2020.06.004
[30] WANG G D, WANG T Y, ZhANG X L, et al. Sex-specific effects of fluoride and lead exposures on histology, antioxidant physiology, and immune system in the liver of zebrafish (Danio rerio)[J]. Ecotoxicology, 2022, 31(3): 396-414. doi: 10.1007/s10646-022-02519-5
[31] BAO S P, TANG W, FANG T. Sex-dependent and organ-specific toxicity of silver nanoparticles in livers and intestines of adult zebrafish[J]. Chemosphere, 2020, 249: 126172. doi: 10.1016/j.chemosphere.2020.126172
[32] 董忠典, 黎学友, 廖健, 等. 雌、雄弓背青鳉 (Oryzias curvinotus) 肝脏转录组比较分析[J]. 海洋与湖沼, 2020, 51(5): 1203-1213. [33] ARUKWE A, KULLMAN S W, HINTON D E. Differential biomarker gene and protein expressions in nonylphenol and estradiol-17β treated juvenile rainbow trout (Oncorhynchus mykiss)[J]. Comp Biochem Physiol C, 2001, 129(1): 1-10.
[34] CHEN X P, LI V W T, YU R M K, et al. Choriogenin mRNA as a sensitive molecular biomarker for estrogenic chemicals in developing brackish medaka (Oryzias melastigma)[J]. Ecotoxicol Environ Saf, 2008, 71(1): 200-208. doi: 10.1016/j.ecoenv.2007.10.005
[35] YE T, KANG M, HUANG Q S, et al. Exposure to DEHP and MEHP from hatching to adulthood causes reproductive dysfunction and endocrine disruption in marine medaka (Oryzias melastigma)[J]. Aquat Toxicol, 2014, 146: 115-126. doi: 10.1016/j.aquatox.2013.10.025
[36] WANG X F, YANG Y J, ZHANG L P, et al. Endocrine disruption by di-(2-ethylhexyl) phthvalate in Chinese rare minnow (Gobiocypris rarus)[J]. Environ Toxicol Chem, 2013, 32(8): 1846-1854. doi: 10.1002/etc.2261
[37] UREN-WEBSTER T M, LEWIS C, FILBY A L, et al. Mechanisms of toxicity of di(2-ethylhexyl) phthalate on the reproductive health of male zebrafish[J]. Aquat Toxicol, 2010, 99(3): 360-369. doi: 10.1016/j.aquatox.2010.05.015
[38] 罗志嘉, 王佩, 彭娜, 等. 水产动物生长性别差异研究进展[J]. 水产学杂志, 2017, 30(6): 56-60. doi: 10.3969/j.issn.1005-3832.2017.06.011 [39] ALVES A M, EUVERINK G J, BIBB M J, et al. Identification of ATP-dependent phosphofructokinase as a regulatory step in the glycolytic pathway of the actinomycete Streptomyces coelicolor A3(2)[J]. Appl Environ Microbiol, 1997, 63(3): 956-961. doi: 10.1128/aem.63.3.956-961.1997
[40] SCHORMANN N, HAYDEN K L, LEE P, et al. An overview of structure, function, and regulation of pyruvate kinases[J]. Protein Sci, 2019, 28(10): 1771-1784. doi: 10.1002/pro.3691
[41] DUMINIL P, DAVANTURE M, OURY C, et al. Arabidopsis thaliana 2,3-bisphosphoglycerate-independent phosphoglycerate mutase 2 activity requires serine 82 phosphorylation[J]. Plant J, 2021, 107(5): 1478-1489. doi: 10.1111/tpj.15395
[42] MCLEISH M J, KENYON G L. Relating structure to mechanism in creatine kinase[J]. Crit Rev Biochem Mol Biol, 2005, 40(1): 1-20. doi: 10.1080/10409230590918577
[43] KÜLTZ D. Molecular and evolutionary basis of the cellular stress response[J]. Annu Rev Physiol, 2005(67): 225-257.
[44] MINTON A P. Confinement as a determinant of macromolecular structure and reactivity[J]. Cell, 1992, 63(3): 1090-1100.
[45] 吴聪颖. 微丝的基本性质与细胞核肌动蛋白[J]. 中国细胞生物学学报, 2019, 41(3): 381-386. [46] 陈剑清, 张耀洲. EF手图像超家族成员-肌钙蛋白C的研究进展[J]. 生物工程学报, 2007, 23(3): 375-380. doi: 10.3321/j.issn:1000-3061.2007.03.003 [47] 吉成龙. 典型溴系阻燃剂对紫贻贝毒理效应的组学研究[D]. 烟台: 中国科学院烟台海岸带研究所, 2014: 107-109. [48] ROSE E, FLANAGAN S P, JONES A G. The effects of synthetic estrogen exposure on the sexually dimorphic liver transcriptome of the sex-role-reversed Gulf pipefish[J]. PLoS One, 2015, 10: e0139401. doi: 10.1371/journal.pone.0139401
[49] CHEN H P, JIANG D N, LI Z Y, et al. Comparative physiological and transcriptomic profiling offers insight into the sexual dimorphism of hepatic metabolism in size-dimorphic spotted scat (Scatophagus argus)[J]. Life (Basel), 2021, 11(6): 589.
-
期刊类型引用(0)
其他类型引用(1)
计量
- 文章访问数: 419
- HTML全文浏览量: 121
- PDF下载量: 50
- 被引次数: 1
粤公网安备 44010502001741号
