Association between PPAR-δ gene SNPs, haplotypes and resistances to SGIV and RGNNV in the orange-spotted grouper, Epinephelus coioides
-
摘要:目的
获得石斑鱼Epinephelus spp.抗病分子标记，选育石斑鱼抗病品系以解决石斑鱼病害频发的问题。
方法基于免疫基因PPAR-δ的基因组DNA序列开展单核苷酸多态性(Single nucleotide polymorphism，SNP)位点筛选，并对这些位点分别进行石斑鱼虹彩病毒(Singapore grouper iridovirus, SGIV)和神经坏死病毒(Red-spotted grouper nervous necrosis virus, RGNNV)抗性的关联分析。
结果在SGIV感染的易感组和抗感组样品中共筛查到9个SNP位点，均位于内含子中；多态信息含量(Polymorphism information content, PIC)的范围为0.177~0.375，其中，SNP-S1(g.940T>A)属于低度多态(PIC<0.25)，其余SNP位点属于中度多态(0.25≤PIC<0.50)；关联分析结果显示，SNP-S7(g.4595T>A)基因型频率在SGIV易感组和抗感组中存在显著差异分布(P<0.05)，SNP-S7的TT和AA这2种纯合子基因型与SGIV抗感性状相关，而AT杂合子基因型与SGIV易感性相关。在RGNNV感染的易感组和抗感组样品中共筛查到8个SNP位点，其中，SNP-N1位于外显子中，属于同义突变，其余SNP均位于内含子中；PIC的范围为0.106~0.317，其中，SNP-N1(g.324G>A)和SNP- N2(g.883A>G)属于低度多态(PIC<0.25)，其余SNP属于中度多态(0.25≤PIC<0.50)；关联分析表明，SNP-N5(g.2510C>T)基因型频率在RGNNV易感组和抗感组中存在显著差异分布(P<0.05)，SNP-N5的CT基因型与RGNNV抗感性状相关, CC基因型与RGNNV易感性状相关。
结论本研究在PPAR-δ的基因组DNA序列中分别筛选到与SGIV和RGNNV抗性相关的SNP标记各1个，可以为石斑鱼抗病育种提供技术支持和理论依据。
Abstract:ObjectiveTo obtain the disease-resistant molecular markers in grouper (Epinephelus spp.), and serve for the selective breeding program of disease-resistant grouper strains, so as to solve the problem of frequent occurrence of grouper disease.
MethodSingle nucleotide polymorphisms (SNPs) were screened based on PPAR-δ genomic DNA sequence, and association analysis of resistance to Singapore grouper iridovirus (SGIV) and red-spotted grouper nervous necrosis virus (RGNNV) was performed on these SNPs.
ResultA total of nine SNPs were detected in the susceptible and resistant groups against SGIV infection, all of which were located in the introns, with the polymorphism information contents (PICs) ranging from 0.177−0.375. Among the SNPs, SNP-S1(g.940T>A) showed low degree polymorphism (PIC<0.25), while the rest SNPs showed moderate degree polymorphism(0.25≤PIC<0.50). The association analysis showed that the genotype frequencies of SNP-S7 (g.4595T>A) were significantly different between SGIV susceptible and resistant groups (P<0.05), the TT and AA homozygous genotypes of SNP-S7 were correlated with SGIV resistance traits, while the AT heterozygous genotypes were correlated with SGIV susceptibility traits. In addition, a total of eight SNPs were detected in the susceptible and resistant groups agasinst RGNNV infection, among which SNP-N1 was located in the exon, with a synonymous mutation, and the rest SNPs were located in the introns, with the PICs ranging from 0.106−0.317. Among the SNPs, SNP-N1 (g.324G>A) and SNP-N2 (g.883A>G) showed low degree polymorphism (PIC<0.25), while the rest SNPs showed moderate degree polymorphism (0.25≤PIC<0.50). The association analysis showed that SNP-N5 (g.2510C>T) genotype frequencies were significantly different between RGNNV susceptible and resistant groups (P<0.05), the CT genotype of SNP-N5 was correlated with RGNNV resistance traits and the CC genotype was correlated with RGNNV susceptibility traits.
ConclusionIn this study, we successfully screened one SNP marker related to SGIV resistance and one SNP marker related to RGNNV resistance from PPAR-δ genomic DNA sequence. This finding can offer a technical support and a theoretical basis for resistance breeding of grouper.
-
外泌体(Exosome, EX)是指包含了复杂RNA和蛋白质，直径为30~150 nm的膜性小囊泡。外泌体存在于体液中，是多泡体(Multivesicular body，MVB)与细胞质膜融合并且以胞外方式降解释放的小泡[1]，参与调控精子体内的生理活动[2]。由于细胞类型不同，外泌体通过旁分泌和自分泌的调节方式，与靶细胞融合或传递其内载物，参与免疫应答等生理和病理过程[3-4]。外泌体携带细胞来源的多种蛋白质、脂类、DNA、mRNA、miRNA等，与细胞间通信、细胞迁移、血管新生和免疫调节等过程有关。精液中富含的外泌体称为精浆外泌体(Seminal plasma exosome，spEX)。研究表明，精浆和细胞外囊泡(外泌体和微囊泡)不仅是精子的载体，在其他生殖生理活动中也起着至关重要的作用[5]。spEX具有改善精液质量、调节精子活力、影响精子受精能力的功能[6]。有研究表明，精浆外泌体产生的胞外三磷酸腺苷(Extracelluar adenosine triphosphate，exATP)可以通过孵育新鲜猪精子达到提高精子活力、调节线粒体代谢并控制精子运动的目的[7]。miRNA(microRNA)是一类由内源基因编码、长度约为20~25个核苷酸的非编码单链RNA分子[8]，广泛存在于生物体各种组织或器官中，调控机体生长发育和细胞凋亡。哺乳动物中，miRNA与靶mRNA不完全互补，可在翻译水平上抑制靶基因表达，调控生命活动[9]。有研究表明，miRNA通过识别并结合目标mRNA的部分互补位点，介导精子和精浆两者细胞间的接触，形成诱导的沉默复合物的一部分，miRNA通过翻译抑制或降解靶mRNA而沉默基因[1]，具有调节基因表达的作用[10]。miRNA可以潜在地调控精子发生和成熟的许多方面，包括增殖、分化发育[11]和凋亡[12]。但目前对spEX与miRNA之间的调控关系研究较少。
随着人工受精技术的普及，种公猪的地位显著提高，公猪的精液品质也愈发受人关注。中国是养猪大国，在繁育过程中，良好的猪精液质量可以带来更大的经济效益，深入了解精子活力的影响因素以及寻找决定精子活力的关键因子成为当前亟需解决的技术难题之一[13]。本研究通过对大白公猪精浆外泌体进行分离提取和鉴定[14]，测序分析spEX miRNA表达谱在精子受精中的调控作用，可以更好地了解精子获能与精浆外泌体之间的表达情况，这些数据可为改善精液质量、提高精子受精能力等方面的研究提供一定参考价值。
1. 材料与方法
1.1 材料
选取5头25~32月龄的优秀大白公猪作为研究对象，所有公猪均来自同一猪场(广西某猪场)，处于相同饲养管理条件下，按时接种疫苗。采精间隔5~7 d。室温手握法进行采精，初次射出的5~10 mL的精液不收集。所有样品均采用计算机辅助精液分析仪(CASA)检测精子的活力和形态[2]。样本在17 ℃、2840 r/min离心10 min使精子和精浆初分离，精浆保存于−40 ℃，用干冰运回实验室。
1.2 方法
1.2.1 外泌体分离
37 ℃条件下解冻精浆，取10 mL精浆于4 ℃、5620 r/min离心30 min收集上清液，继续在4 ℃、11240 r/min离心70 min分离出细胞外大囊泡，得到的每份样品液依次使用0.45 μm、0.22 μm针孔滤器过滤，将过滤后的液体在4 ℃、33200 r/min离心80 min，去除上清液后加入PBS溶液反复吹打，白色胶絮状物质即为外泌体，进一步在4 ℃、33200 r/min离心80 min纯化外泌体，加400 μL PBS溶液于−80 ℃存放备用。
1.2.2 透射电镜分析
取10 μL外泌体样品置于碳涂层铜网格上吸附20 min，干燥后避光，用15 μL醋酸双氧铀负染1 min。在Hitachi H-7700透射电子显微镜下拍照，检测外泌体的形态。
1.2.3 纳米颗粒跟踪分析
取5 μL外泌体样品稀释至30 μL，用ZetaView Nanoparticle Tracker测量外泌体的平均直径和浓度，分析外泌体的粒径大小分布情况。
1.2.4 Western blot鉴定分析
取外泌体样品加入适量的SDS-PAGE缓冲液，99 ℃变性10 min。用SDS-PAGE凝胶快速配置试剂盒(上海碧云天生物有限公司)，进行制胶、电泳、转膜、封闭、孵育一抗Calnexin (Abcam， ab133615)、Hsp70 (Abcam， ab181606)、β-Actin (Beyotime，ab8226)和Tsg101(Beyotime，ab125011)，一抗均已以1:1000的体积比稀释，然后根据说明书配制显色液，通过成像系统对条带进行曝光拍照。
1.2.5 miRNA的提取
用miRNeasy®Mini试剂盒(Qiagen)从分离鉴定取得的样品中提取和纯化miRNA。提取后的miRNA用Agilent Bioanalyzer 2100系统(Agilent Technologies)及配套的RNA Nano 6000 Assay Kit进行质控。
1.2.6 cDNA文库构建及测序
用QIAseq miRNA(Library Kit，Qiagen)构建文库。5个单独测序样品在属性序列中额外加入索引代码，最终经Bioanalyzer 2100(Agilent)和qPCR评估miRNA文库构建的质量，并用Illumina Hiseq2000 (Illumina)进行miRNA测序。
1.2.7 miRNA数据处理及分析
用Bowtie软件，将测序得到的miRNA分别与Silva、GTRNAdb、Rfam和Repbase数据库进行序列比对，筛选出miRNA并计算miRNA的表达量(Transcript per million，TPM)，公式为TPM=Readcount ×106/MappedReads(Readcount表示比对到某一miRNA的Read数目；MappedReads表示比对到所有miRNA上的Read数目)。
1.2.8 miRNA靶基因的预测、GO分类及KEGG富集分析
根据鉴定得到的高表达量的miRNA，基于3′UTR作为靶向序列，用miRanda和RNAhybrid软件预测miRNA的靶基因，用GO( http://www.geneonUology.org)和KEGG( http://www.genome.jp/kegg/)对预测的靶基因进行功能注释和通路分析。
2. 结果与分析
2.1 精浆外泌体的鉴定
CASA系统对大白公猪原始精液的精子活力检测(表1)，结果显示精子活力良好，符合国家制定的种猪常温精液质量标准[15]。透射电子显微镜观察精浆中所分离的细胞外囊泡，可见明显的立体膜结构散在分布，呈圆形或椭圆形，大小较为均一，其形态具有典型的茶杯托样特征(图1)；粒径分析显示囊泡粒子平均粒径为96.9 nm，粒径主峰占比为94.9%(图2)；Western blot结果显示囊泡包含Tsg101、Hsp70等外泌体标志性膜蛋白，Tsg101和Hsp70蛋白表达为阳性(图3)。这些数据与前人的研究一致，表明精浆中所分离的细胞外囊泡为外泌体。
表 1 精液质量基本信息Table 1. Basic information of semen quality公猪编号 Boar number 精子活力/% Sperm vitality 精子畸形率/% Semen malformation rate 精子密度/ mL−1Semen density 1 95 3 4.18×108 2 97 1 3.68×108 3 96 2 4.74×108 4 97 2 3.87×108 5 95 5 5.22×108 ![]()
图 1 精浆外泌体透射电镜检测结果Figure 1. Transmission electron microscopic detection of seminal plasma exosomes![]()
图 2 精浆外泌体粒径分析精浆外泌体的平均粒径为96.9 nm；粒径主峰为85.2 nm；主峰占检测所有峰的比例为94.9%；颗粒含量为3.4×107 mL−1Figure 2. Particle size analysis of seminal plasma exosomesThe average particle diameter of seminal plasma exosomes was 96.9 nm; The particle diameter of main peak was 85.2 nm; The main peak accounted for 94.9% of all the detected peaks; The particle content was 3.4×107 mL−1![]()
图 3 精浆外泌体标志性蛋白的Western blot检测2个条带为5个样品中随机挑选的2个样品Figure 3. Western blot assay of marker proteins in seminal plasma exosomesThe two strips were two samples randomly selected out of five samples2.2 猪精浆外泌体中miRNA测序和鉴定
对5头大白公猪spEX样品miRNA测序，检测到Raw read数量为108849665 ，Clean read数量为60882083，碱基质量值(Q30)范围为96.89%~97.40%，共鉴定出329个已知miRNA。miRNA片段长度为19~24 bp，其中22 bp的片段占比超过45%。根据计算miRNA的表达量TPM，筛选前20个代表性miRNA(表2)。
表 2 精浆外泌体中标志性miRNATable 2. Representative miRNAs in seminal plasma exosomesmiRNA 长度/bp Length 序列(5′→3′) Sequence TPM ssc-let-7c 22 UGAGGUAGUAGGUUGUAUGGUU 851579.7906 ssc-let-7a 22 UGAGGUAGUAGGUUGUAUAGUU 655084.4766 ssc-let-7f-5p 22 UGAGGUAGUAGAUUGUAUAGUU 513356.5287 ssc-let-7e 22 UGAGGUAGGAGGUUGUAUAGUU 426156.8525 ssc-miR-10a-5p 22 UACCCUGUAGAUCCGAAUUUGU 181651.3508 ssc-miR-10b 22 UACCCUGUAGAACCGAAUUUGU 181595.6635 ssc-miR-141 21 UAACACUGUCUGGUAAAGAUG 163842.7352 ssc-miR-21-5p 22 UAGCUUAUCAGACUGAUGUUGA 125987.7562 ssc-miR-125b 22 UCCCUGAGACCCUAACUUGUGA 118310.9973 ssc-miR-200b 23 UAAUACUGCCUGGUAAUGAUGAC 99945.7705 ssc-miR-16 22 UAGCAGCACGUAAAUAUUGGCG 96309.6121 ssc-miR-30a-5p 22 UGUAAACAUCCUCGACUGGAAG 93885.9728 ssc-miR-30d 24 UGUAAACAUCCCCGACUGGAAGCU 74152.9510 ssc-miR-148a-3p 22 UCAGUGCACUACAGAACUUUGU 73369.8562 ssc-miR-191 23 CAACGGAAUCCCAAAAGCAGCUG 65102.6570 ssc-miR-26a 22 UUCAAGUAAUCCAGGAUAGGCU 59547.6681 ssc-miR-26b-5p 22 UUCAAGUAAUUCAGGAUAGGUU 53532.4383 ssc-miR-125a 23 UCCCUGAGACCCUUUAACCUGUG 51427.2644 ssc-let-7i-5p 19 UGAGGUAGUAGUUUGUGCU 50378.3688 ssc-miR-29c 22 UAGCACCAUUUGAAAUCGGUUA 44726.2430 2.3 miRNA靶基因预测
用miRanda和RNAhybrid软件预测和筛选得到miR-10家族(miR-10b和miR-10a-5p)、let-7家族(let-7c、let-7a、let-7f-5p和let-7e)等前10个高表达miRNA的靶基因。利用Cytoscape软件对miRNA与靶基因互作分析进行可视化(图4)。富集结果显示，miRNA中高表达的miR-let-7c-5p潜在调控染色质重塑蛋白MORC家族CW型锌指结构蛋白2(MORC family CW-type zinc finger 2，MORC2)、人转化生长因子β受体3(Transforming growth factor beta receptor 3，TGFBR3)、DEAD-box RNA解旋酶3(DEAD-box RNA helicase 3 X-linked，DDX3X)、核糖核苷酸还原酶M2 (Ribonucleotide reductase regulatory subunit M2，RRM2)等多个靶基因；MiR-125a-5p潜在调控镁离子依赖的蛋白磷酸酶(Protein phosphatase, Mg2+/Mn2+dependent 1D，PPM1D)和双特异性酪氨酸磷酸化调节激酶−2(Dual specificity tyrosine phosphorylation regulated kinase 2，DYRK2)，且和miR-let-7c-5p联合调控靶基因鼠双微体4(Murine double minute 4，MDM4)、锌指基质蛋白型蛋白3(Zinc finger matrin-type protein 3，ZMAT3)和胰岛素样生长因子1(Insulin like growth factor 1，IGF1)；miR-16-5p潜在调控血小板衍生生长因子(Platelet derived growth factor subunit B，PDGFB)、核糖体蛋白S6激酶多肽3(Ribosomal protein S6 kinase A3，RPS6KA3)、神经母细胞瘤鼠肉瘤同系物(NRAS proto-oncogene, GTPase，NRAS)、胰岛素样生长因子1受体(Insulin like growth factor 1 receptor，IGF1R)等多个靶基因；miR-let-7a-5p潜在调控活化白细胞黏附分子(Activated leukocyte cell adhesion molecule，ALCAM)、精子相关抗原9(Sperm associated antigen 9，SPAG9)、肝细胞生长因子(Hepatocyte growth factor，HGF)、类剑蛋白p60亚基A1(Katanin catalytic subunit A1 like 1，KATNAL1)等多个靶基因；miR-let-7i-5p潜在调控靶基因Th2型炎症因子(Interleukin 13，IL13)、α−心脏肌动蛋白(Actin alpha cardiac muscle 1，ACTC1)、脂联素2型受体(Adiponectin receptor 2，ADIPOR2)和B淋巴细胞瘤−2(BCL2 apoptosis regulator，BCL2)；miR-125b-5p潜在调控E2F家族转录因子2(E2F transcription factor 2，E2F2)、受体亚家族2C成员2(Nuclear receptor subfamily 2 group C member 2，NR2C2)、碱性神经酰胺酶2(Alkaline ceramidase 2，ACER2)、Ⅲ型纤维连接蛋白域蛋白3A(Fibronectin type III domain containing 3A，FNDC3A)等多个靶基因；其中泛素特异性蛋白酶47(Ubiquitin specific peptidase 47，USP47)、B细胞淋巴瘤−2(BCL2 apoptosis regulator，BCL2)、细胞周期蛋白依赖激酶抑制剂1A(Cyclin dependent kinase inhibitor 1A，CDKN1A)被多个miRNA联合调控。
![]()
图 4 精浆外泌体miRNA与靶基因互作分析橘色方框代表miRNA，蓝色圆圈代表miRNA靶向的靶基因Figure 4. Interaction analysis between miRNAs in seminal plasma exosomes and target genesThe orange box represents miRNA, and the blue circle represents the target gene of the miRNA2.4 猪精浆外泌体中miRNA功能分析
GO和KEGG富集分析结果显示，spEX中miRNA的调控功能主要表现在射精、P53信号通路、前列腺癌、细胞对DNA损伤刺激的反应、负调控凋亡过程、跨膜运输、质膜外侧、顶体膜结合、受精等通路中发挥潜在调控作用(图5)。其中，射精通路中，miR-16-5p的靶基因包括ACVR2A和SLC6A4；P53信号通路中，let-7c-5p的靶基因包括BCL2L1、IGF1、CCND1和CDKN1A；细胞对DNA损伤刺激反应通路中，miR-125a-5p的靶基因包括MORC2、BTG2和TRAF6；精子受精通路中，let-7c-5p的靶基因包括BCL2L1和FNDC3A；质膜外侧的通路中，miR-16-5p的靶基因是TGFBR3，miR-200b-5p的靶基因是IL12RB2，miR-148a-3p的靶基因是ALCAM，let-7i-5p的靶基因是IL13，let-7a-5p的靶基因是THBS1，miR-125a-5p的靶基因是LIFR，miR-125b-5p的靶基因是CDH5。
![]()
图 5 精浆外泌体miRNA靶基因功能GO和KEGG富集分析Figure 5. GO and KEGG enrichment analyses of miRNA target gene functions of seminal plasma exosomes3. 讨论与结论
外泌体可携带多种核酸、蛋白质、脂质以及小分子物质等重要信息物质，形成细胞间通信，直接影响细胞迁移、免疫反应和细胞分化生长等过程。外泌体具有脂质双层膜结构，其组成成分和功能具有特异性，被其包裹的内容物可稳定存在[16]。外泌体的形成以及被包裹内容物的分拣、释放是一系列精细调控的过程，需要很多蛋白质和能量的参与。细胞通过不断吸收释放外泌体，实现了与外界微环境的信息交互，参与到各种生理病理疾病中。本研究提取的大白猪spEX平均粒径约为96.9 nm，电镜下形态呈茶杯托状结构，Western blot结果显示大白公猪精浆外泌体中Tsg101和Hsp70 蛋白表达为阳性，这与其他报道EX的特征结果一致[17]，表明已成功提取到大白猪精液中的spEX。
精液中的EX主要是睾丸、附睾、副性腺分泌物的混合物，即精浆外泌体。其成分因动物物种而异[18]。在外泌体运载的生物分子中，越来越多miRNA的作用被发掘关注。有研究表明miRNA是外泌体中含量最丰富的细胞间信息交流的主要分子之一[19]，外泌体携带的生物活性调节分子参与精子发生、修饰加工及受精能力的调控[20]。一些miRNA的数量和种类一旦改变会导致配子发育停滞、精卵受精障碍和各种生殖系统疾病的产生，显示出EX携带的miRNA在配子发育和精子加工中起重要作用[20-21]。本研究主要对大白公猪spEX miRNA的表达情况进行检测，并对前10个高丰度的miRNA相关功能进行预测，分析发现高表达的前10个miRNA中let-7家族占了4个(let-7c、let-7a、let-7f-5p、let-7e)，显著高于其他miRNA。其中，let-7家族在哺乳动物、昆虫和植物中高度表达和保守，是基本生物学过程中最重要的miRNA调节因子之一[22]。有研究表明在公猪睾丸中miR-let-7家族的成员(let-7a-5p、let-7f-5p、let-7c-5p)高度表达，这与本研究预测结果基本一致，let-7家族可能与公猪睾丸的发育和睾丸中精子发生有关[23]，在生精过程中，高水平的let-7表达可能会改变正常生精所需的蛋白质的产生，从而影响与精子形态和运动性相关的蛋白质[24]。miR-let-7a在猪睾丸[23]、小鼠附睾精液[25]、人精浆[26]中均有发现，提示其可抑制维甲酸诱导的精原细胞分离并介导胚胎附着[27]。精子获能的过程需要能量，miR-26a与精子能量代谢有关，其可能影响精子活力和凋亡[28]，有研究表明PDHX基因可能是miR-26a的潜在靶点，其参与丙酮酸氧化脱酸连接糖酵解和氧化磷酸化的关键步骤，可以调节糖代谢途径，最终影响公猪精子的活力和存活率[29]。范宇[30]研究表明，miR-26a可以调节PTEN表达，通过靶向调控PTEN减少精子细胞凋亡和促存活作用，这可能是雄性不育的原因。在公猪中，miR-26a被多个环状RNA靶向以调节精子活力，这也表明miR-26a 具有作为精子质量和雄性生育能力的非侵入性生物标志物的潜力[31]。精子的发生包括3个主要发育阶段：精原干细胞自我更新、精原干细胞分化和精子细胞变形，在每个阶段中均检测到miRNA表达，许多miRNA参与调节睾丸的发育阶段和控制精子发生[32-33]。有研究报道血浆外泌体中高表达的miRNA没有发现miR-10家族[34]，这排除了miR-10家族来自血浆的可能性。有研究分析miR-10b在不同年龄的小鼠睾丸中表达，且在睾丸发育过程中表达降低[35]，Li等[36]研究发现miR-10b主要在thy1细胞中表达，并且在出生后第6天的小鼠睾丸中高表达，这些结果表明miR-10b在精原干细胞有丝分裂的调节中起作用，可能与精子发生有关。本研究结果显示spEX也检测到miR-10家族，说明了miR-10家族可能是在附睾或副性腺等生殖器官中特异性表达。
GO和KEGG通路分析表明，高表达miRNA的靶基因在P53信号通路、精子受精、负调控凋亡过程、跨膜运输和能量代谢等方面具有潜在的调控作用。P53是细胞内最重要的蛋白之一，能介导细胞应答诸多类型的应激信号。P53蛋白通过参与促进DNA损伤修复、诱导细胞周期阻断、细胞衰老凋亡等通路来维持基因组的稳定性和完整性[37]。有研究表明，P53信号通路与男性生殖能力有关[23]，这也说明了精浆外泌体中的miRNA可能与公猪繁育有关。有些前列腺小体(前列腺分泌的外泌体)可以特异性地将遗传物质和生物活性物质传递给精子，随后可以提高精液质量和公猪的受精率[38]。有研究发现，前列腺癌患者血清中miR-141的含量明显增高，且与特异性抗原(PSA)相比较，两者具有相关性，表明miRNA在前列腺癌的发生中起着非常重要的作用[39]。因此这些表达失调的miRNA可能可以用于癌症的诊断、治疗以及治疗效果的评价。此外，通路分析表明，高表达miRNA的靶基因参与了氧化磷酸化和泛素介导的蛋白水解等活动。蛋白质磷酸化对精子活力至关重要，研究显示，核糖体蛋白S6激酶A3(RPS6KA3)与磷酸化(V-ATPase)运输调节有关，RPS6KA3磷酸化各种底物，这种蛋白质的活性与控制细胞生长和分化有关[40]。泛素蛋白酶体系统可通过降解有毒蛋白质聚集体和错误折叠的蛋白质来参与调节精子蛋白质的稳态[41]。在附睾运输过程中，附睾起始段的精子蛋白与附睾尾段的精子蛋白不一致，意味着运输过程中精子蛋白选择性丢失[42]。泛素介导的蛋白水解过程可能在精子附睾成熟阶段发挥作用。
综上，本研究报道大白公猪精浆外泌体高表达的前10个miRNA，通过靶基因预测和通路富集分析，预测了这些miRNA可能与精子的发生、P53信号通路、线粒体外膜调节、精细胞发育等过程有关，为大白公猪精浆外泌体中miRNA基因表达调控研究提供基础数据，这将有助于确定与精子的生育能力有关的基因及利用miRNA作为生物标志物预测公猪的生育力。但本研究结果基于生物信息学预测，所涉miRNA、靶基因及其相关通路还需进一步试验验证。
-
![]()
图 1 石斑鱼PPAR-δ基因结构及SNP位点
Figure 1. DNA structure of PPAR-δ gene and its SNP loci of grouper
![]()
图 2 石斑鱼SGIV和RGNNV感染组中筛选出PPAR-δ基因SNP标记的连锁不平衡分析
图片顶部的白色长方形代表PPAR-δ基因组DNA序列，黑色细线对应标签代表各SNP位点；图片下半部分的正方形色块内的数值为连锁不平衡参数r2；正方形色块的颜色代表两两SNP位点间的连锁不平衡参数D的范围，白色表示D=0，红色和粉色表示0 < D < 1，蓝色表示D=1；黑框代表该区域内相邻的SNP位点处于连锁不平衡状态，形成单倍型块，用Block标记，括号内的数值表示单倍型块碱基的长度
Figure 2. Linkage disequilibrium analysis of PPAR-δ SNP markers screened from SGIV and RGNNV infected groupers
The white rectangle at the top of image represents the PPAR-δ genome DNA sequence, and the black thin lines and corresponding tags represent the SNP sites; In the bottom half of the picture, the value in the square color block is the linkage unbalance parameter r2; The color of the square color block represents the range of linkage unbalance parameter D between two SNP sites, white represents D=0, red and pink represent 0 < D < 1, and blue represents D=1; The black box indicates that the adjacent SNP sites in this region are in the linkage unbalanced state and form a haplotype block, which is marked by Block, and the value in brackets indicates the length of haplotype block base
表 1 PPAR-δ基因组DNA序列扩增引物
Table 1 Primers for PPAR-δ genomic DNA sequence amplification
引物名称 Primer 引物序列(5′→3′) Primer sequence θ/℃ 产物长度/bp Product length PPAR-δ-1-F CAGGGCTTTTTCCGACGAAC 59 1325 PPAR-δ-1-R GGTTGACTCATCACTCGCCT 59 PPAR-δ-2-F GACCTCCTGTGGCCCTTAAC 60 1398 PPAR-δ-2-R CTGGGAGCGGCTTAAAGGAA 60 PPAR-δ-3-F AAGCCGCTCCCAGTAGAAAC 60 1168 PPAR-δ-3-R TTCCTAGCAATGGCAGCGTT 60 PPAR-δ-4-F CAAACGCTGCCATTGCTAGG 60 1324 PPAR-δ-4-R ACGTTCATTAGCCCAGGACG 60 
下载: 导出CSV
表 2 9个SNP在石斑鱼SGIV抗感组和易感组中的遗传多态性信息
Table 2 Genetic polymorphisms of nine SNPs in SGIV resistant and susceptible grouper groups
组别 Group SNP 观测杂合度 Observed heterozygosity (Ho) 期望杂合度 Expected heterozygosity (He) 多态信息含量 Polymorphism information content (PIC) 有效等位基因数 Number of effective alleles (Ne) Hardy–Weinberg平衡 Hardy-Weinberg equilibrium (H-W) P1) 抗感组 Resistant group SNP-S1 0.160 0.149 0.136 1.173 0.328 0.567 SNP-S2 0.460 0.398 0.317 1.651 1.244 0.265 SNP-S3 0.400 0.407 0.322 1.676 0.016 0.898 SNP-S4 0.660 0.489 0.367 1.937 6.279 0.012* SNP-S5 0.640 0.492 0.369 1.950 4.612 0.032* SNP-S6 0.640 0.492 0.369 1.950 4.612 0.032* SNP-S7 0.240 0.424 0.332 1.724 9.699 0.002* SNP-S8 0.300 0.495 0.370 1.962 7.932 0.005* SNP-S9 0.280 0.492 0.369 1.950 9.676 0.002* 易感组 Susceptible group SNP-S1 0.220 0.198 0.177 1.244 0.688 0.407 SNP-S2 0.460 0.481 0.363 1.908 0.094 0.759 SNP-S3 0.280 0.347 0.284 1.523 1.923 0.166 SNP-S4 0.500 0.500 0.372 1.980 0.000 1.000 SNP-S5 0.500 0.500 0.372 1.980 0.000 1.000 SNP-S6 0.520 0.502 0.373 1.987 0.007 0.796 SNP-S7 0.520 0.453 0.348 1.814 1.109 0.292 SNP-S8 0.360 0.492 0.369 1.950 3.681 0.055 SNP-S9 0.360 0.504 0.375 1.997 4.220 0.040* 1)“*”表示该位点理论基因型频率和观测基因型频率之间的差异显著(P<0.05，χ2检验) 1)“*” indicates significant difference between the theoretical genotype frequency and the observed genotype frequency of the site (P<0.05,χ2 test)
下载: 导出CSV
表 3 8个SNP在石斑鱼RGNNV抗感组和易感组中的遗传多态性信息
Table 3 Genetic polymorphisms of eight SNPs in RGNNV resistant and susceptible grouper groups
组别 Group SNP 观测杂合度 Observed heterozygosity (Ho) 期望杂合度 Expected heterozygosity (He) 多态信息含量 polymorphism information content (PIC) 有效等位基因数 Number of effective alleles (Ne) Hardy–Weinberg平衡 Hardy-Weinberg Equilibrium (H-W) P1) 抗感组 Resistant group SNP-N1 0.180 0.164 0.150 1.200 0.431 0.512 SNP-N2 0.160 0.147 0.136 1.170 0.328 0.567 SNP-N3 0.640 0.435 0.341 1.770 10.700 0.001* SNP-N4 0.540 0.466 0.358 1.870 1.100 0.294 SNP-N5 0.620 0.442 0.344 1.800 7.720 0.005* SNP-N6 0.680 0.448 0.349 1.810 12.800 0.000* SNP-N7 0.660 0.442 0.344 1.790 11.700 0.001* SNP-N8 0.600 0.420 0.332 1.720 8.830 0.003* 易感组 Susceptible group SNP-N1 0.160 0.149 0.136 1.170 0.328 0.567 SNP-N2 0.120 0.114 0.106 1.130 0.168 0.682 SNP-N3 0.480 0.369 0.298 1.570 4.750 0.029* SNP-N4 0.380 0.379 0.304 1.600 0.001 0.982 SNP-N5 0.380 0.358 0.291 1.550 0.200 0.655 SNP-N6 0.440 0.389 0.311 1.630 0.900 0.343 SNP-N7 0.420 0.398 0.317 1.650 0.155 0.694 SNP-N8 0.420 0.335 0.277 1.500 3.340 0.068 1)“*”表示该位点理论基因型频率和观测基因型频率之间的差异显著(P<0.05，χ2检验) 1)“*” indicates significant difference between the theoretical genotype frequency and the observed genotype frequency of the site (P<0.05,χ2 test)
下载: 导出CSV
表 4 石斑鱼PPAR-δ基因在2种病毒抗感组和易感组的SNP统计分析
Table 4 PPAR-δ gene SNPs analysis in resistant and susceptible grouper groups infected by two viruses
SNP (SGIV) 基因型1) Genotype 基因型频率/% Genotype frequency χ2 SNP (RGNNV) 基因型1) Genotype 基因型频率/% Genotype frequency χ2 抗感组 Resistant group 易感组 Susceptible group 抗感组 Resistant group 易感组 Susceptible group SNP-S1 TT(81) 84 78 0.585 SNP-N1 GA(17) 18.8 16.3 0.099 AT(19) 16 22 GG(80) 81.3 83.7 SNP-S2 CC(44) 50 38 4.418 SNP-N2 AA(83) 83.3 87.8 0.384 GC(46) 46 46 GA(14) 16.7 612.2 GG(10) 4 16 SNP-N3 AT(54) 62.5 49.0 1.796 SNP-S3 AC(34) 40 28 1.680 TT(43) 37.5 51.0 AA(58) 52 64 SNP-N4 AA(45) 37.5 55.1 2.520 CC(8) 8 8 GA(44) 52.1 38.8 SNP-S4 CC(28) 26 3 3.818 GG(8) 10.4 6.1 CT(58) 66 50 SNP-N5 CC(46) 37.5 57.1 4.260* TT(14) 8 2 CT(48) 60.4 38.8 SNP-S5 AA(28) 26 30 2.669 TT(3) 2.1 4.1 AT(57) 64 50 SNP-N6 CC(2) 0 4.1 3.820 TT(15) 10 20 GC(54) 66.7 44.9 SNP-S6 GA(58) 64 52 2.324 GG(41) 33.3 51.0 GG(27) 26 28 SNP-N7 AA(3) 0 6.1 3.410 AA(15) 10 20 GA(52) 64.6 42.9 SNP-S7 TT(50) 60 40 8.491* GG(42) 35.4 51.0 AT(38) 24 52 SNP-N8 CT(49) 58.3 42.9 2.320 AA(12) 16 8 TT(48) 41.7 57.1 SNP-S8 TT(26) 28 24 0.451 AT(33) 15 18 AA(41) 42 40 SNP-S9 TT(39) 44 34 1.176 AT(32) 28 36 AA(29) 28 30 1)括号中的数字为基因型数量；2)使用多组样品的χ2检验统计，“*”表示该SNP位点不同基因型在抗感组和易感组的分布差异显著，即该位点不同基因型与抗感/易感性状相关联(P<0.05) 1) Numbers in the brackets are numbers of the genotypes; 2) χ2 test statistics of samples from multiple groups are used, “*” indicates the distribution of different genotypes of the SNP locus is significantly different between the resistant group and the susceptible group, that is, different genotypes of this locus are associated with the resistant/susceptible traits (P<0.05)
下载: 导出CSV
表 5 石斑鱼PPAR-δ基因SNP连锁单倍型与2种病毒抗性的关联分析
Table 5 Association analysis between PPAR-δ SNPs linked haplotype and resistance of two viruses
病毒 Virus 单倍块 Haploblock 单倍型 Haplotype 分布频率/% Distribution frequency 抗感组频率/% Frequency in resistant group 易感组频率/% Frequency in susceptible group χ2 P1) SGIV 单倍块1 Block 1 TTT 42.50 43.00 42.00 0.020 0.890 AAA 32.00 30.00 34.00 0.370 0.540 TAA 13.00 15.00 14.00 0.040 0.840 TAT 12.50 12.00 10.00 0.200 0.650 单倍块2 Block 2 CAG 56.00 58.00 54.00 0.330 0.570 TTA 43.00 41.00 45.00 0.000 0.990 CAA 1.00 1.00 1.00 0.000 1.000 RGNNV 单倍块1 Block 1 TACGGT 23.20 20.00 27.40 2.950 0.086 AGTCAC 12.70 13.30 11.90 0.177 0.674 TGCGGT 12.70 13.80 11.30 0.528 0.468 TACCAT 12.50 12.40 12.50 0 0.987 AGTCAT 13.50 14.20 12.50 0.245 0.621 AACGGT 12.00 12.00 11.90 0.001 0.977 TGTCAT 13.50 14.20 12.50 0.245 0.621 1) P>0.05表示该位点的单倍型分布频率在抗感组和易感组之间差异不显著(χ2检验) 1) P>0.05 indicates the distribution frequency of haplotype of the site is not significantly different between the resistant group and the susceptible group (χ2 test)
下载: 导出CSV
-
[1] 罗鸣, 陈傅晓, 刘龙龙, 等. 我国石斑鱼养殖疾病的研究进展[J]. 水产科学, 2013, 32(9): 549-554. doi: 10.3969/j.issn.1003-1111.2013.09.010 [2] QIN Q W, CHANG S F, NGOH-LIM G H, et al. Characterization of a novel ranavirus isolated from grouper Epinephelus tauvina[J]. Diseases of Aquatic Organisms, 2003, 53(1): 1-9.
[3] COSTA J Z, THOMPSON K D. Understanding the interaction between Betanodavirus and its host for the development of prophylactic measures for viral encephalopathy and retinopathy[J]. Fish and Shellfish Immunology, 2016, 53: 35-49. doi: 10.1016/j.fsi.2016.03.033
[4] 黄剑南, 郭志勋, 冯娟, 等. 鱼类诺达病毒及其所导致的疾病[J]. 水产学报, 2006, 30(6): 831-836. [5] VIGNAL A, MILAN D, SANCRISTOBAL M, et al. A review on SNP and other types of molecular markers and their use in animal genetics[J]. Genetics Selection Evolution, 2002, 34(3): 275-305. doi: 10.1186/1297-9686-34-3-275
[6] 高焕, 赖晓芳, 孟宪红, 等. 中国明对虾抗菌肽基因应答WSSV侵染的表达及其SNP分析[J]. 中国水产科学, 2011, 18(3): 646-653. [7] BAO Y, LI L, ZHANG G. Polymorphism of the superoxide dismutase gene family in the bay scallop (Argopecten irradians) and its association with resistance/susceptibility to Vibrio anguillarum[J]. Developmental and Comparative Immunology, 2010, 34(5): 553-561. doi: 10.1016/j.dci.2009.12.016
[8] FU G H, BAI Z Y, XIA J H, et al. Characterization of the LECT2 gene and its associations with resistance to the big belly disease in Asian seabass[J]. Fish and Shellfish Immunology, 2014, 37(1): 131-138. doi: 10.1016/j.fsi.2014.01.019
[9] DUNN S E, BHAT R, STRAUS D S, et al. Peroxisome proliferator-activated receptor delta limits the expansion of pathogenic Th cells during central nervous system autoimmunity[J]. Journal of Experimental Medicine, 2010, 207(8): 1599-1608. doi: 10.1084/jem.20091663
[10] SUN L, SHI Y, WANG G, et al. PPAR-delta modulates membrane cholesterol and cytokine signaling in malignant B cells[J]. Leukemia, 2018, 32(1): 184-193. doi: 10.1038/leu.2017.162
[11] NEEL J G, GRIMALDI P A. Physiological functions of peroxisome proliferator-activated receptor beta[J]. Physiological Reviews, 2014, 94(3): 795-858. doi: 10.1152/physrev.00027.2013
[12] JI J D, KIM H J, RHO Y H, et al. Inhibition of IL-10-induced STAT3 activation by 15-deoxy-Delta12, 14-prostaglandin J2[J]. Rheumatology, 2005, 44(8): 983-988. doi: 10.1093/rheumatology/keh657
[13] WANG Y, YU Y, WANG Q, et al. PPAR-δ of orange-spotted grouper exerts antiviral activity against fish virus and regulates interferon signaling and inflammatory factors[J]. Fish & Shellfish Immunology, 2019, 94: 38-49.
[14] HUANG X, HAUNG Y, SUN J, et al. Characterization of two grouper Epinephelus akaara cell lines: Application to studies of Singapore grouper iridovirus (SGIV) propagation and virus-host interaction[J]. Aquaculture, 2009, 292(3/4): 172-179.
[15] YANG M, WEI J, LI P, et al. MHC polymorphism and disease resistance to Singapore grouper iridovirus (SGIV) in the orange-spotted grouper, Epinephelus coioides[J]. Science Bulletin, 2016, 61(9): 693-699. doi: 10.1007/s11434-016-1055-5
[16] YANG M, WANG Q, CHEN J, et al. Identification of candidate SNPs and genes associated with anti-RGNNV using GWAS in the red-spotted grouper, Epinephelus akaara[J]. Fish & Shellfish Immunology, 2021, 112: 31-37.
[17] NORTON N, WILLIAMS N M, WILLIAMS H J, et al. Universal, robust, highly quantitative SNP allele frequency measurement in DNA pools[J]. Human Genetics, 2002, 110(5): 471-478. doi: 10.1007/s00439-002-0706-6
[18] PARENTEAU J, ABOU ELELA S. Introns: Good day junk is bad day treasure[J]. Trends in Genetics, 2019, 35(12): 923-934. doi: 10.1016/j.tig.2019.09.010
[19] SHAUL O. How introns enhance gene expression[J]. International Journal of Biochemistry & Cell Biology, 2017, 91: 145-155.
[20] BOTSTEIN D, WHITE R L, SKOLNICK M, et al. Construction of a genetic linkage map in man using restriction fragment length polymorphisms[J]. American Journal of Human Genetics, 1980, 32(3): 314-331.
[21] HUBERT S, BUSSEY J T, HIGGINS B, et al. Development of single nucleotide polymorphism markers for Atlantic cod (Gadus morhua) using expressed sequences[J]. Aquaculture, 2009, 296(1/2): 7-14.
[22] 周欣, 高风英, 卢迈新. 鱼类抗病育种研究进展[J]. 大连海洋大学学报, 2021, 36(3): 510-523. doi: 10.16535/j.cnki.dlhyxb.2020-156 [23] WAN Q Y, SU J G, CHEN X H, et al. Genomic sequence comparison, promoter activity, SNP detection of RIG-I gene and association with resistance / susceptibility to grass carp reovirus in grass carp (Ctenopharyngodon idella)[J]. Developmental & Comparative Immunology, 2013, 39(4): 333-342.
[24] HENG J F, SU J G, HUANG T, et al. The polymorphism and haplotype of TLR3 gene in grass carp (Ctenopharyngodon idella) and their associations with susceptibility / resistance to grass carp reovirus[J]. Fish & Shellfish Immunology, 2011, 30(1): 45-50.
[25] SHEN Y B, ZHANG J B, XU X Y, et al. A new haplotype variability in complement C6 is marginally associated with resistance to Aeromonas hydrophila in grass carp[J]. Fish & Shellfish Immunology, 2013, 34(5): 1360-1365.
[26] 高风英, 卢迈新, 曹建萌, 等. 尼罗罗非鱼 NOD1 基因 SNP 位点和单倍型与抗无乳链球菌感染的关联分析[J]. 农业生物技术学报, 2018, 26(11): 1949-1961. [27] FU G H, LIU F, XIA J H, et al. The LBP gene and its association with resistance to Aeromonas hydrophila in tilapia[J]. International Journal of Molecular Sciences, 2014, 15(12): 22028-22041. doi: 10.3390/ijms151222028
[28] FU G H, WAN Z Y, XIA J H, et al. The MCP-8 gene and its possible association with resistance to Streptococcus agalactiae in tilapia[J]. Fish & Shellfish Immunology, 2014, 40(1): 331-336.
[29] LATRUFFE N, VAMECQ J. Peroxisome proliferators and peroxisome proliferator activated receptors (PPARs) as regulators of lipid metabolism[J]. Biochimie, 1997, 79(2/3): 81-94.
[30] 纪永吉. 牦牛PPARα、PPARγ、PPARδ基因的多态性与生长性状相关性研究[D]. 兰州: 甘肃农业大学, 2016. [31] 丁克越. 不同单倍型块定义在单倍型块结构推断与htSNPs选择中的效应[D]. 北京: 中国协和医科大学, 2004.
计量
- 文章访问数: 84
- HTML全文浏览量: 10
- PDF下载量: 20
