林麝及其近缘物种编码区微卫星分布规律及功能分析
收稿日期: 2021-12-01
录用日期: 2022-07-05
网络出版日期: 2022-12-02
基金资助
国家自然科学基金(31702032);重庆市自然科学基金(cstc2019jcyj-msxmX0410);重庆市基本科研项目(cstc2021jxjl-jbky10008)
Distinct patterns of microsatellite and functional analysis of forest musk deer and its closely related species
Received date: 2021-12-01
Accepted date: 2022-07-05
Online published: 2022-12-02
麝科和鹿科动物均属于偶蹄反刍类动物，具有重要的经济价值。通过系统的微卫星序列 (Simple sequences repeats, SSRs) 从基因组水平揭示物种间的系统进化关系，探索微卫星序列的基因功能及其富集的信号通路，目前仍缺乏相关研究。随着林麝 (Moschus berezovskii)、原麝 (Moschus moschiferus)、小麂 (Muntiacus reevesi)、赤麂 (Muntiacus vaginalis) 和马鹿 (Cervus elaphus) 基因组测序的完成，本文利用生物信息学方法提取了这些动物蛋白质编码区 (coding sequences, CDS) 序列，统计和分析了其CDS区微卫星序列分布规律及其生物学功能，探索了含SSR 基因富集的信号通路及其与疾病的关联性。结果表明，林麝、原麝、小麂、赤麂和马鹿蛋白质编码区含SSR序列的基因所占比例分别为6.96% (1 696个)、7.18% (2 359个)、7.29% (3 005个)、7.36% (1 916个) 和7.48% (1 924个)，并且这5种动物CDS区SSRs分布模式具有相似性，均是三倍体核苷酸 (即三核苷酸和六核苷酸) SSRs最多，分别为96.85%、94.87%、65.44%、64.23%和88.04%。GO功能富集表明，林麝与其他4种动物蛋白质编码区SSR序列在分子功能、细胞组成和生物学过程3个方面具有较多共同显著富集的功能，包括DNA结合、染色质和生长发育等。KEGG 通路分析表明，林麝及其他4种动物蛋白质编码区SSR序列具有7个共同显著富集的KEGG通路，包括遗传信息调控蛋白家族、转录因子、染色体及相关蛋白、剪接体、转录机制和Notch信号通路和成体糖尿病。通过对林麝编码区含SSR关键免疫基因及其相关联的KEGG通路进行分析，发现10个含SSR的关键免疫基因对应的KEGG通路与疾病密切相关。
赵琪 , 张琪 , 李浩玲 , 兰月 , 鄢行安 , 赵贵军 , 戚文华 . 林麝及其近缘物种编码区微卫星分布规律及功能分析[J]. 兽类学报, 2022 , 42(6) : 705 -715 . DOI: 10.16829/j.slxb.150645
Muscidae and Cervidae belong to artiodactyl ruminants and have important economic value. There is still a lack of relevant studies on using simple sequence repeats (SSRs) to reveal the phylogenetic relationship among species at the genomic level and to explore the gene functions of SSRs and their enriched signaling pathways. With the completion of the genome sequence of the forest musk deer (Moschus berezovskii), Siberian musk deer (Moschus moschiferus), the red deer (Cervus elaphus), the Reeves’ muntjac (Muntiacus reevesi), and the red muntjac (Muntiacus vaginalis), we extracted the coding sequences (CDS) of these species by bioinformatics method and analyzed the distribution patterns of coding SSR and their biological functions. The containing-SSR genes were enriched in the signal pathways and their association with disease were explored. The results prove that the proportion of their coding SSRs were 6.96% (1 696), 7.18% (2 359), 7.29% (3 005), 7.36% (1 916), and 7.48% (1 924), respectively. In addition, the distribution pattern of coding SSRs of the five species was similar, triplet repeats SSR (Tri- and Hexa-SSRs) of the five species were the most similar, which were 96.85%, 94.87%, 65.44%, 64.23%, and 88.04%, respectively. The GO functional enrichment showed that the SSR sequences in the CDS of the five species had many same enriched functions in three aspects: molecular function, cellular component, and biological process, including DNA binding, chromatin, growth, and development. The result of the KEGG pathway enrichment analysis showed that there are 7 KEGG pathways that were significantly co?enriched in forest musk deer and other species.These include protein families: genetic information processing, transcription factors, chromosome and associated proteins, spliceosome, transcription machinery, Notch signaling pathway, and maturity?onset diabetes of the young. By analyzing the SSR-containing key immune genes and their associated KEGG pathways in the CDS region of forset musk deer, 10 SSR-containing key immune genes corresponding to the KEGG pathway were found to be closely associated with the disease.
Key words： Forest musk deer; Siberian musk deer; Red deer; Small suede; Red muntjac; CDS; SSR; Functional analysis
| null | Bana N á, Nyiri A, Nagy J, Frank K, Nagy T, Stéger V, Schiller M, Lakatos P, Sugár L, Horn P, Barta E, Orosz L. 2018. The red deer Cervus elaphus genome CerEla1.0: sequencing, annotating, genes, and chromosomes. Molecular Genetics and Genomics, 293 (3): 665-684. |
| null | Chen C J, Chen H, Zhang Y, Thomas H R, Margaret H F, He Y H, Xia R. 2020. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Molecular Plant, 13 (8): 1194-1202. |
| null | Du L, Li Y, Zhang X, Yue B. 2013. A user?friendly program for reporting distribution and building databases of microsatellites from genome sequences. The Journal of Heredity, 104: 154-157. |
| null | Du L M, Zhang C, Liu Q, Zhang X Y, Yue B S. 2018. Krait: an ultrafast tool for genome-wide survey of microsatellites and primer design. Bioinformatics, 34 (4): 681-683. |
| null | Duval A, Hamelin R. 2002. Mutations at coding repeat sequences in mismatch repair-deficient human cancers: toward a new concept of target genes for instability. Cancer Research, 62 (9): 2447-2454. |
| null | Duval A, Gayet J, Zhou X P, Iacopetta B, Thomas G, Hamelin R. 1999. Frequent frameshift mutations of the TCF-4 gene in colorectal cancers with microsatellite instability. Cancer Research, 59 (17): 4213-4215. |
| null | Fan Z X, Li W J, Jin J, Cui K, Yan C C, Peng C J, Jian Z Y, Bu P, Price M, Zhang X Y, Shen Y M, Li J, Qi W H, Yue B S. 2018. The draft genome sequence of forest musk deer (Moschus berezovskii). Giga Science, 7 (4): giy038. |
| null | Jiang X M, Hu T Z, Xiang X S, Qi W H. 2015. Statistics and bioinformatics analysis of microsatellite sequence in poplar whole genomes. Southwest China Journal of Agricultural Sciences, 28 (2): 527-533. (in Chinese) |
| null | Jin P, Warren S T. 2000. Understanding the molecular basis of fragile X syndrome. Human Molecular Ggenetics, 9 (6): 901-908. |
| null | Kashi Y, King D, Soller M. 1997. Simple sequence repeats as a source of quantitative genetic variation. Trends in Genetics, 13 (2): 74-78. |
| null | Li W J, Li Y Z, Du L M, Huang J, Shen Y M, Zhang X Y, Yue B S. 2014. Comparative analysis of microsatellite sequences distribution in the genome of giant panda and polar bear. Sichuan Journal of Zoology, 33 (6): 874-878 (in Chinese) |
| null | Li Y C. 2004. Microsatellites within genes: structure, function, and evolution. Molecular Biology and Evolution, 21 (6): 991-1007. |
| null | Lu T, Wang C, Du C, Liu S, Shen Y M, Zhang X Y, Yue B S. 2017. Distribution regularity of microsatellites in Moschus berezovskii genome. Sichuan Journal of Zoology, 36 (4): 420-424. (in Chinese) |
| null | Metzgar D, Wills C. 2000. Evidence for the adaptive evolution of mutation rates. Cell, 101 (6): 581-584. |
| null | Montgelard C, Catzeflis F, Douzery E. 1997. Phylogenetic relationships among cetartiodactyls andcetaceans as deduced from the comparison of cytochrome b and 12S RNA mitochondrial sequences. Molecular Biology Evolution, 14: 550-559. |
| null | O’Dushlaine C T, Edwards R J, Park S D, et al. 2005. Tandem repeat copy-number variation in protein-coding regions of human genes. Genome Biology, 6 (8): R69. |
| null | Ohtaishi N, Gao Y. 1990. A review of the distribution of all species of deer (Tragulidae, Moschidae and Cervidae) in China. Mammal Review, 20 (2-3): 125-144. |
| null | Qi W H, Jiang X M, Xiao G S, Huang X Y, Du L M. 2013. Seeking and bioinformatics analysis of microsatellite sequence in the genomes of cow and sheep. Acta Veterinaria et Zootechnica Sinica, 44 (11): 1724-1733. (in Chinese) |
| null | Qi W H, Lu T, Zheng C L, Jiang X M, Jie H, Zhang X Y, Yue B S, Zhao G J. 2020.Distribution patterns of microsatellites and development of its marker in different genomic regions of forest musk deer genome based on high throughput sequencing. Aging, 12 (5): 4445-4462. |
| null | Sermon K, Seneca S, Rycke M D, Goossens V, Van de Velde H, De Vos A, Platteau P, Lissens W, Van Steirteghem A, Liebaers I. 2001. PGD in the lab for triplet repeat diseases?myotonic dystrophy, Huntington’s disease and Fragile?X syndrome. Molecular and Cellular Endocrinology, 183: S77-S85. |
| null | Schwartz S, Yamamoto H, Navarro M, Maestro M, Reventós J, Perucho M. 1999. Frameshift mutations at mononucleotide repeats in caspase - 5 and other target genes in endometrial and gastrointestinal cancer of the microsatellite mutator phenotype. Cancer Research, 59 (12): 2995-3002. |
| null | Subramanian S, Mishra ＲK, Singh L. 2003. Genome?wide analysis of microsatellite repeats in humans: their abundance and density in specific genomic regions. Genome Biology, 4 (2): 1-10. |
| null | Timchenko N A, Iakova P, Cai Z J, Smith J R, Timchenko L T. 2001. Molecular basis for impaired muscle differentiation in myotonic dystrophy. Molecular and Cellular Biology, 21 (20): 6927-6938. |
| null | Vassileva V, Millar A, Briollais L, Chapman W, Bapat B. 2002. Genes involved in DNA repair are mutational targets in endometrial cancers with microsatellite instability. Cancer Research, 62 (14): 4095-4099. |
| null | Wang X H, Pan X F, Li H Q, Duan F. 2016. Advances in the studies of the expansion of (CAG)n·(CTG)n trinucleotide repeats and mecha-nisms underlying its related diseases. International Journal of Genetics, 39 (5): 274-281. (in Chinese) |
| null | Wang Z L, Huang J, Du L M, Li W J, Yue B S, Zhang X Y. 2013. Comparison of microsatellites between the genomes of Tetranychus urticae and Ixodes scapularis . Sichuan Journal of Zoology, 32 (4): 481-486. (in Chinese) |
| null | Yi L, Dalai M G, Su R N, Lin W L, Erdenedalai M, Luvsantseren B, Chimedtseren C, Wang Z, Hasi S R. 2020. Whole?genome sequencing of wild Siberian musk deer (Moschus moschiferus) provides insights into its genetic features. BMC Genomics, 21(1): 108. |
| null | Yin Y, Fan H Z, Zhou B T, Hu Y B, Fan G Y, Wang J H, Zhou F, Nie W H, Zhang C Z, Liu L, Zhong Z Y, Zhu W B, Liu G C, Lin Z S, Liu C, Zhou J, Huang G P, Li Z H, Yu J P, Zhang Y L, Yang Y, Zhuo B Z, Zhang B W, Chang J, Qian H Y, Peng Y M, Chen X Q, Chen L, Li Z P, Zhou Q, Wang W, Wei F W. 2021. Molecular mechanisms and topological consequences of drastic chromosomal rearrangements of muntjac deer. Nature Communication, 12 (1): 6858. |
| null | Yuki M, Masumi I, Shujiro O, Yoshizawa A C, Minoru K. 2007. KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Research, 35 (Web Server issue): W182-W185. |
| null | 丁戈, 姚南, 吴琼, 刘恒, 郑国锠. 2008. 着丝粒结构与功能研究的新进展. 植物学通报, 25 (2): 149-160. |
| null | 王希恒, 潘学峰, 李红权, 段斐. 2016. (CAG)n·(CTG)n三核苷酸重复序列扩增及相关疾病机制研究进展. 国际遗传学杂志, 39 (5): 274-281. |
| null | 卢婷, 王晨, 杜超, 刘姝, 沈咏梅, 张修月, 岳碧松. 2017. 林麝全基因组微卫星分布规律研究. 四川动物, 36 (4): 420-424. |
| null | 刘志霄, 盛和林. 2000. 我国麝的生态研究与保护问题概述. 动物学杂志, 35 (3): 54-57. |
| null | 李午佼, 李玉芝, 杜联明, 黄杰, 沈咏梅, 张修月, 岳碧松. 2014. 大熊猫和北极熊基因组微卫星分布特征比较分析. 四川动物, 33 (6): 874-878. |
| null | 杨述林, 王志刚, 樊斌, 刘榜, 李奎. 2003. 微卫星及其功能研究. 湖北农业科学, 2003 (6): 91-93. |
| null | 汪自立, 黄杰, 杜联明, 李午佼, 岳碧松, 张修月. 2013. 二斑叶螨和肩突硬蜱基因组微卫星分布规律研究. 四川动物, 32 (4): 481-486. |
| null | 戚文华, 蒋雪梅，肖国生, 黄小云, 杜联明. 2013. 牛和绵羊全基因组微卫星序列的搜索及其生物信息学分析. 畜牧兽医学报, 44 (11): 1724-1733. |
| null | 蒋雪梅, 胡廷章, 向兴胜, 戚文华. 2015. 杨树全基因组微卫星序列的统计及其生物信息学分析. 西南农业学报, 28 (2): 527-533. |
/
| 〈 |
|
〉 |
青公网安备 63010402000199号 青ICP备05000010号-2
