獭兔,典型的皮用兔,因被毛独特而倍受人们喜爱。近年来,在生产中发现,某些獭兔会在腹部和四肢以及臀部出现皱褶,称其为皱襞型獭兔。皱襞型獭兔的个体比非皱襞型稍大,同时具有正常獭兔相同的皮毛品质,这对于以皮张面积和皮毛品质为重要经济效益因素的獭兔来说,意义重大,提供了新的獭兔育种方向。本研究采用RNA-Seq技术,对皱襞型(试验组)和非皱襞型(对照组)白色獭兔的皮肤样本进行测序分析。分析测序结果,筛选出与皮肤生长发育相关的基因。将样本的有效reads合并,通过Trinity软件,使用两端重复进行de novo拼接。两组皮肤样本测序总共获得377,618条Unigenes, Unigenes的平均长度是680 nt,最长的序列长达20,293 nt。将Unigenes与NR、Swiss-Prot、KEGG和KOG等数据库比对(E-Value<10"5),根据比对的注释,分析已知Unigene的功能。由于家兔的全基因组信息还不完善,85.5931%序列无法比对到任何数据库上。根据NR注释及Blast2GO,总共有42,895条Unigene被映射到不同的GO term上,并按照分子功能、细胞组分及生物过程将其分类；KOG的结果显示60,593条Unigenes被分为25类。KEGG数据库比对显示,13,431条Unigene分别被注释到399个KEGG pathway中,其中注释到PI3K-Akt信号通路(ko04151)的基因最多(686,5.108%),该信号通路在皮肤的生长发育过程中具有重要的意义。两种不同表型的样本中,共有308个显著的差异表达基因(P-Value≤0.05且Iog2|foldchange|> 1),其中118个在皱襞型样品中为上调基因,190个为下调基因。这些基因中有99个被注释到不同的GO Term上,其中只有29个基因被注释到KEGG不同的代谢途径中。除此之外的209个基因被称为新基因,在这些基因中,33个只在非皱襞型个体中表达,56个基因只在皱襞型个体中表达。通过qRT-PCR验证测序结果。选择7个候选基因,使用oligo 7.0设计特异性引物,采用△△Ct法计算皱襞型和非皱襞型样品的相对表达量。结果显示,qRT-PCR的结果和转录组测序的结果相符。说明测序结果是正确的,据此初步确定与家兔皮肤发育相关的基因：MYL12B、LAMBS、PLOW、EIF3E、CEBPA、FAD2, ENPP6等。候选基因MYL12B是细胞骨架的重要组分,作为肌球蛋白的调节性轻链,可通过调节轻链的磷酸化与去磷酸化,调节细胞的许多动力学过程。本研究克隆了MYL12B基因的CDS序列,并与人、黑猩猩、小鼠等12个物种的CDS序列和氨基酸序列进行比较。发现,除马外,该基因与哺乳动物CDS相似性高达94%,与鸟类和两栖类的同源性也较高。生物信息学分析显示MYL12B编码的蛋白是一个不稳定的、具有一定亲水性的结构蛋白,同时富含磷酸化和糖基化位点。同时该基因的CDS序列在皱襞型和非皱襞型的个体中没有序列的差异,表达量的差异可能是由肌球蛋白作用机制决定的。本研究采用PCR-SSCP技术,设计候选基因LAMB3第5、6、7、8、9和10等外显子的特异性引物,通过聚丙烯酰胺凝胶电泳检测,将不同带型的PCR产物送生物工程公司测序,结果发现LAMB3第8外显子存在1个潜在的突变位点,通过统计分析,这些突变位点在皱襞型和非皱襞型白色獭兔的皮肤组织样品中具有显著性差异。本研究首次将转录组测序技术应用于白獭兔皱襞性状的探索中,初步确定与皱襞性状相关的基因,并通过该基因的功能注释分析,了解该基因在皮肤的生长发育过程中的作用方式和途径,为进一步研究该基因在分子和蛋白质层面的功能和作用机制提供了依据。