两性生殖动物具有两大细胞系,即生殖细胞系和体细胞系,生殖细胞系负责遗传信息的传递,而体细胞系维持生物个体生长和发育
[1]。原始生殖细胞(primordial germ cells,PGCs)是生殖细胞的祖细胞,为早期胚胎发育过程中最早建立的一群生殖干细胞,其形成后逐渐迁移到性腺原基,经过有丝分裂和减数分裂最终发育为成熟的配子
[1]。
dnd (dead end)基因是脊椎动物生殖细胞的特异性标记基因,编码RNA结合蛋白,对原始生殖细胞的形成和迁移至关重要
[2]。科研人员对多种鱼类的研究结果表明,
dnd基因功能的缺失会导致PGCs无法正常迁移,造成生殖细胞缺失,从而发育为不育的个体
[3-5]。
dnd基因最初是在斑马鱼(
Danio rerio)中发现的,敲除斑马鱼
dnd基因会导致PGCs异常迁移,但体细胞的发育不受影响
[6]。人们对非洲爪蟾(
Xenopus laevis)的研究表明,
dnd基因仅在卵巢中表达,且当
dnd基因的转录和翻译受到抑制时,PGCs无法正常迁移
[7]。在鱼类中,
dnd基因的表达模式相对保守,为性腺特异表达,例如在稀有鮈鲫(
Gobiocypris rarus)
[8]、异育银鲫(
Carassius auratus gibelio)
[9]、褐牙鲆(
Paralichthys olivaceus)
[10]、大菱鲆(
Scophthalmus maximus)
[11]和中华鲟(
Acipenser sinensis)
[12]中,均能在性腺中检测到
dnd基因的表达,且在卵巢的表达量高于精巢。另外,
dnd基因在爬行动物中华鳖(
Pelodiscus sinensis)
[13]中也具有相似的表达模式。上述研究表明,
dnd基因在性腺发育过程中发挥着重要的作用。
大鳍鳠(
Hemibagrus macropterus)俗称石扁头、石胡子、江鼠等,隶属于硬骨鱼纲、鲇形目、鲿科、鳠属,广泛分布于长江和珠江各水系的干流及其支流中,在长江上游江段渔获物中曾占有一定比例
[14]。大鳍鳠肉质细嫩、味道鲜美、营养丰富,具有较高的食用价值。目前科研人员已经开展了大鳍鳠人工繁殖、养殖、胚胎及仔稚鱼发育方面的研究
[15-17],而对其生殖细胞了解甚少。因此,本试验采用RACE方法克隆大鳍鳠
dnd基因(
Hmdnd)的全长cDNA序列,通过实时荧光定量PCR(RT-qPCR)研究了
Hmdnd基因在大鳍鳠不同组织和发育时期性腺中的表达特征,并通过原位杂交技术鉴定了
Hmdnd基因在性腺细胞的定位特征,以期为大鳍鳠生殖细胞的鉴定及发育研究奠定基础。
1 材料与方法
1.1 材 料
试验所用大鳍鳠样本来自于重庆市水产研究所长寿基地,1龄雌鱼和雄鱼体长分别为(19.5±1.2)和(18.7±0.4) cm,体质量分别为(45.5±2.1)和(41.1±2.6) g/尾;2龄雌鱼和雄鱼体长分别为(25.4±0.4)和(27.5±1.8) cm,体质量分别为(123.3±9.8)和(125.6±2.9) g/尾;4龄雌鱼和雄鱼体长分别为(44.6±0.6)和(41.9±0.7) cm,体质量分别为(383.7±21.2)和(381.5±2.5) g/尾。取1,2,4龄鱼各6尾(雌雄各半),断尾放血后分别采集心脏、肝脏、脾脏、肾脏、肠、脑、鳃及性腺组织精巢和卵巢,置于RNA保存液(Vivacell,中国)中4 ℃过夜,然后保存于-20 ℃冰箱,分别用于后续RNA提取。同时,取每个个体的性腺,用波恩氏液(Scientific Phygene,中国)和40 g/L多聚甲醛(paraformaldehyde,PFA)固定,用于后续石蜡切片制备和原位杂交试验。所有试验程序均按照中国水产科学院长江水产研究所制定的《实验动物养护和使用指导原则》实施。
1.2 Hmdnd 基因cDNA全长序列克隆与分析
1.2.1 序列克隆
用FastPure Cell/Tissue Total RNA Isolation Kit V2(Vazyme,中国)试剂盒提取4龄大鳍鳠卵巢组织总RNA,用琼脂糖凝胶电泳和微量分光光度计(Thermo Scientific,美国,NanoDrop One)检测RNA质量和浓度。按照SMARTer® RACE 5′/3′ Kit(TaKaRa,日本)说明书,以大鳍鳠卵巢RNA为模板,分别合成SMART 3′ cDNA和5′cDNA,置于-20 ℃冰箱保存备用。
根据大鳍鳠转录组获得
dnd基因部分序列,用Prime5软件设计RACE引物(
表1),进行巢式PCR扩增。以合成的卵巢SMART cDNA为模板,使用引物Hmdnd-5′RACE-1和Hmdnd-3′RACE-1与RACE通用引物UPM分别进行5′RACE和3′RACE的第1次PCR扩增;将PCR产物稀释10倍,作为第2次PCR的模板,用引物Hmdnd-5′RACE-2和Hmdnd-3′ RACE-2与RACE通用引物NUP分别进行5′RACE和3′ RACE的第2次PCR扩增。反应体系:上下游引物各0.5 μL,mix Vazyme 2×
Taq Plus Master MIXⅡ(Dye Plus)(Vazyme,中国)10 μL,模板cDNA 1 μL,无菌去离子水8 μL。设置PCR反应条件为:95 ℃预变性5 min;95 ℃变性30 s,59 ℃退火30 s,72 ℃延伸10 s,35次循环;72 ℃终延伸10 min。用1.5%的琼脂糖凝胶电泳检测PCR产物,并回收目的片段,然后连接到pMD19-T载体上(16 ℃连接12~16 h),将连接产物转入大肠杆菌感受态细胞DH5α,在超净工作台中用涂布棒将菌液均匀涂布在含有氨苄青霉素的培养皿上,37 ℃培养过夜,次日挑选单个菌落进行PCR鉴定,PCR产物进行1.5%琼脂糖凝胶电泳检测。将PCR鉴定阳性样品送至武汉天一华煜基因公司测序。最后用DNAMAN 9.0软件根据3′和5′端序列所包含的重叠序列进行拼接,得到
Hmdnd基因全长cDNA序列。
1.2.2 序列分析
通过NCBI网站上的ORF finder 在线工具(
https://www.ncbi.nlm.nih.gov/orffinder)获得
Hmdnd基因开放阅读框(open read frame,ORF)。使用DNAMAN 9.0软件将
Hmdnd基因核苷酸序列翻译成氨基酸序列,并通过Vector NTI 11.5.1比对分析核苷酸序列;通过NCBI数据库工具(Protein BLAST:Align two or more sequences using BLAST(nih.gov))计算物种间dnd蛋白序列的一致性;采用Mega 11.0软件中的Clustal W进行比对后,用Neighbor-joining法进行系统发育分析,Bootstrap值设置为1 000;通过MEGA 11.0软件和ESPript 3.0(
https://espript.ibcp.fr/ESPript/cgi-bin/ESPript.cgi)进行氨基酸序列多重比对和绘图。
1.3 Hmdnd基因在不同组织中的表达特征
提取1龄雌、雄大鳍鳠不同组织的总RNA,使用微量分光光度计(Thermo Scientific,美国,NanoDrop One)测定RNA浓度,通过1.5%琼脂糖凝胶电泳检测RNA质量,选择260和280 nm处吸光度比值(A260/A280)在1.8~2.0,且电泳条带完整的样品用于后续cDNA的合成,质检合格的总RNA于-80 ℃冰箱储存备用。使用TaKaRa Prime Script RT reagent Kit with gDNA Eraser Kit对RNA进行反转录,获得cDNA,-20 °C冰箱保存,用于RT-qPCR试验。
根据
Hmdnd基因的基因组和cDNA序列设计出跨内含子的RT-qPCR引物Hmdnd-RT-F/R,选择
β⁃actin作为内参基因并设计引物
β-
actin-RT-F/R(
表1)。通过标准曲线和熔解曲线进行引物筛选,所有定量引物扩增效率为90%~110%,
R2≥0.99。RT-qPCR反应体系:上、下游引物各0.35 μL,SYBR Mix 10 μL,ddH
2O 8.3 μL,模板cDNA 1 μL,共20 μL。反应程序:50 ℃ 2 min,95 ℃ 10 min;95 ℃ 15 s,58 ℃ 15 s,72 ℃ 20 s,循环40次;60 ℃ 1 min,95 ℃ 15 s。最终采用2
-ΔΔCt法计算得出各组织的表达量。
1.4 Hmdnd基因在不同发育时期性腺中的表达特征
对波恩氏液固定的1,2,4龄大鳍鳠性腺样本进行脱水包埋后制备石蜡切片,用苏木精-伊红染液进行染色,封片后在显微镜(DM5000B,Leica)下观察并拍照。确定性腺样本处于不同发育时期后,进行RNA提取和反转录,参照1.3节中引物Hmdnd-RT-F/R和反应体系进行RT-qPCR试验,通过2-ΔΔCt法得到Hmdnd基因在不同发育时期性腺的表达量。
1.5 Hmdnd mRNA在性腺中的细胞定位分析
采用原位杂交试验分析Hmdnd mRNA在精巢和卵巢中的细胞定位。性腺样本经40 g/L多聚甲醛固定过夜后,在4 ℃冰箱中用0.3 kg/L的蔗糖溶液进行渗透,次日用OCT包埋剂(Tissue-Tek)包埋,最终使用Leica低温恒温切片机制备冰冻切片。
探针合成:在大鳍鳠cDNA序列的3′端设计探针引物ISH-dnd(
表1),取1龄大鳍鳠的卵巢cDNA进行PCR扩增,产物经琼脂糖凝胶电泳后回收目的片段,获得制备探针的线性化模板,使用NanoDrop One微量分光光度计和1.5%琼脂糖凝胶电泳分别检测其浓度和质量后进行探针合成。使用Digoxigenin (DIG) RNA Labeling Kit 合成正义和反义探针,经DNase和LiCl纯化后,用1.5%琼脂糖凝胶电泳检测合成探针的质量。
用PBS将冰冻切片复水,经蛋白酶K(10 μg/mL)室温消化10 min后,再用40 g/L 多聚甲醛固定15 min,最后加入含有探针的杂交液65 ℃杂交16~20 h;次日,在65 ℃孵育箱中用wash buffer、2×SSCT和0.2×SSCT充分摇洗,再用MABT室温摇洗后加入杂交封闭阻断液封闭3~4 h,然后加入Fab-AP抗地高抗体反应液4 ℃过夜;用MABT洗干净玻片后加入BCIP/NBT显色液对探针信号进行避光显色,显色过程中每间隔1 h在显微镜下观察一次,待显色达到合适效果后用体积分数50%的甘油封片,在显微镜下拍照。
1.6 数据分析
利用IBM SPSS Statistics 27.0.1(IBM,美国)对试验数据进行单因素方差分析(One-way ANOVA),采用LSD和Duncan法进行多重比较,分析组间的差异显著性,显著水平定为P<0.05。所有结果以“平均值±标准差”表示。使用Graphpad prism 5.0软件(Graphpad,美国)制图。
2 结果与分析
2.1 Hmdnd基因克隆与序列特征分析
经RACE-PCR克隆测序获得了2 035 bp的
Hmdnd基因全长cDNA序列(GenBank登录号为PQ626051),包括5′端起始密码子上游长度为100 bp的非编码区、3′端终止密码子下游长度为804 bp的非编码区和编码376个氨基酸的1 131 bp的ORF(
图1),其中高度保守的RNA识别结构域(RNA recognition motif,RRM)为主要的功能结构域(Glu
54-Val
126)。
对大鳍鳠dnd蛋白与其他物种同源蛋白序列进行多重比对,结果表明,大鳍鳠与其他物种dnd氨基酸序列相似,具有6个保守的结构域,其中包括1个N端结构域(N-terminal region,NR)、1个高度保守的RRM和4个C端结构域(C-terminal region(CR),CR1~CR 4)(
图2)。大鳍鳠与黄颡鱼dnd氨基酸序列的一致性最高(79.24%),与其他硬骨鱼类如大西洋鲑鱼、青鳉、斑马鱼序列的一致性分别为51.45%,50.00%和48.42%;与两栖动物非洲爪蟾序列的一致性为37.71%;与爬行动物中华鳖和哺乳动物人、鼠、鸡序列的一致性比较低,分别为37.23%,37.43%,35.59%和35.81%。
dnd蛋白系统进化分析表明,除外类群黑腹果蝇外,脊椎动物dnd蛋白主要分为2支,一支由硬骨鱼类组成,其中大鳍鳠和黄颡鱼聚在一起,其遗传距离最近;另一支由四足类动物组成(
图3)。
2.2 Hmdnd基因在不同组织和发育时期性腺的表达特征
利用RT-qPCR方法检测1龄大鳍鳠不同组织中
Hmdnd mRNA表达水平,结果(
图4)显示,
Hmdnd仅在精巢和卵巢中表达,且在卵巢的表达量显著高于精巢(
P<0.05)。
通过石蜡切片确定性腺的发育时期,结果表明,1龄大鳍鳠卵巢主要由卵黄发生前期的Ⅰ期和Ⅱ期初级卵母细胞组成(
图5-A);此时,精巢处于有丝分裂时期,主要由精原细胞组成(
图5-D)。2龄大鳍鳠卵巢主要由卵黄发生期的Ⅲ期初级卵母细胞组成(
图5-B);此时,精巢处于减数分裂时期,主要由初级精母细胞、次级精母细胞和少量的精子细胞组成(
图5-E)。4龄大鳍鳠卵巢主要由卵黄成熟期的Ⅳ期初级卵母细胞组成(
图5-C);此时,精巢处于精子生成时期,主要由精子组成(
图5-F)。RT-qPCR检测结果表明,
Hmdnd基因在1龄大鳍鳠卵黄发生前期的卵巢(
图5-G)和有丝分裂时期的精巢(
图5-H)中的表达量最高,随着配子发生的进行,
Hmdnd表达量显著降低(
图5-G和H,
P<0.05)。
2.3 Hmdnd mRNA在性腺中的细胞定位
为了进一步了解
Hmdnd基因在性腺的表达特征,利用RNA原位杂交技术研究其在2龄大鳍鳠卵巢和精巢细胞的定位表达情况。结果显示,
Hmdnd主要在生殖细胞的细胞质中表达。卵巢中,
Hmdnd在Ⅰ期初级卵母细胞有少量表达,随着初级卵母细胞发育到Ⅱ期,其表达量急剧增加(
图6-A);然后,当初级卵母细胞发育至Ⅲ期时,
Hmdnd表达量恢复到Ⅰ期初级卵母细胞的表达水平。精巢中,
Hmdnd mRNA在精原细胞的表达量最高,在初级精母细胞中的表达量次之,而在次级精母细胞和精子细胞中几乎不表达(
图6-B)。
Hmdnd mRNA正义探针在卵巢和精巢中均未检测到任何表达信号(
图6-C,D)。
3 讨 论
本研究采用RACE克隆方法获得
Hmdnd基因全长cDNA序列,通过氨基酸序列比对分析发现,大鳍鳠dnd蛋白具有6个保守结构域,表明这几个结构域在进化过程中非常保守。dnd蛋白系统进化分析结果表明,大鳍鳠与同为鲿科鱼类的黄颡鱼首先聚为一支,亲缘关系最近,随后与其他硬骨鱼类聚为一支,这与其形态学上的分类关系相一致
[18]。结果表明,
Hmdnd基因在进化过程中非常保守。
鱼类
dnd基因的研究表明,
dnd mRNA仅在性腺中表达,且在卵巢的表达量高于精巢
[19]。在银鲫
[9]和中华鲟
[12]中,
dnd基因在卵巢的表达量大约是精巢的2倍。同样地,在斑马鱼
[7]、西里伯斯青鳉(
Oryzias celebensis)
[20-21]、星突江鲽(
Platichthys stellatus)
[22]和黄鳝(
Monopterus albus)
[23]中的研究也发现,
dnd基因在卵巢的表达量显著高于精巢。本研究中,大鳍鳠
Hmdnd基因在卵巢的表达量远高于精巢,且表达量随着配子发育的推进呈现下降的趋势,在成熟性腺的表达量最低。对太平洋蓝鳍金枪鱼(
Thunnus orientalis)
[24]的研究发现,
dnd基因在早期和晚期精巢的表达量较低,在青春期精巢的表达量最高。上述结果表明,虽然鱼类
dnd基因只在性腺中表达,但是不同鱼类间的表达特征存在一定差异。
基因的细胞表达定位对于了解基因的功能至关重要。目前,通过RNA原位杂交研究发现,鱼类
dnd基因主要在生殖细胞的细胞质中表达。对稀有鮈鲫
[8]、异育银鲫
[9]、中华鲟
[12]、西里伯斯青鳉
[20]和大西洋牙鲆(
Paralichthys dentatus)
[25]的研究表明,
dnd mRNA在卵巢卵原细胞和Ⅰ期初级卵母细胞中弱表达,在Ⅱ期初级卵母细胞中强表达,之后随着卵子的发生逐渐减弱。在精巢中,多数鱼类的
dnd mRNA在精原细胞中的表达最强,在初级精母细胞和次级精母细胞中次之,在精子细胞和精子中不表达
[8-9,12,20]。本研究中,大鳍鳠
Hmdnd基因也在性腺的生殖细胞中特异表达,且其动态表达趋势与上述鱼类的研究结果相类似。然而,在中华鳖中,
dnd在初级和次级精母细胞中表达最强,在精原细胞中减弱,在精子细胞中不表达
[13]。这可能是由于中华鳖属于爬行动物,与鱼类的亲缘关系相对较远,在进化过程中,
dnd mRNA的表达特征发生了变化。总体而言,
dnd基因为生殖细胞特异表达基因,其在生殖细胞的发育过程中起着重要作用。
本研究结果表明,大鳍鳠Hmdnd基因在生殖细胞中特异性表达,在配子发生过程中其表达水平呈现动态变化特征,说明其在生殖细胞的发育过程中发挥重要作用,但其具体作用机制有待通过基因编辑技术来进一步阐明。
重庆市科研机构绩效激励引导专项(CSTB2023JXJL-YFX0084)
重庆市财政项目(C20240103)
国家重点研发计划项目(2022YFD2400101)
中央级公益性科研院所基本科研业务费专项资金项目(2023TD23)