山羊LF基因序列分析及SNPs研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
本试验根据GenBank上已发表的牛的乳铁蛋白(Lactoferrin, LF)基因DNA序列和山羊mRNA序列进行引物设计,扩增山羊LF基因的部分5’UTR至外显子17序列片断,测序总长度为30306bp,已提交到GenBank,序列号为FJ609300。为比较反刍动物LF基因遗传变异与分化,同时测定了绵羊LF基因相应序列30340bp,序列号为:FJ541507。
     根据测序结果检测到山羊LF基因序列上共30个突变位点,其中处于外显子的有1个,内含子的29个。依据序列号FJ609300上的位置来命名这些突变位点,处于外显子上的1个点突变为7605C>T,引起氨基酸Arg136Trp突变。处于内含子上的29个突变位点分别为2428C>G、2761C>G、6053A>G、6538A>C、7363A>G、8642C>T、8655C>T、9284C>T、13719C>T、15783A>G、15820A>G、16820A>G、17774A>T、17775A>G、17793A>G、17808A>G、17855C>G、18266C>G、18428A>G、18930A>G、19764C>T、19997C>T、20740A>G、20989C>T、21146C>T、21281C>G、23146C>T、24300C>T、24720A>G、25249A>G。同时检测到绵羊LF基因序列上16个突变位点,都是内含子变异,依据序列号FJ541507上的位置来命名这些突变位点分别为:2432C>T、2764C>G、4197C>T、4274A>T、4633C>T、11038C>T、13045A>G、13075A>G、13482A>G、18072A>G、20660C>T、21022C>T、21625T>G、21706T>G、22583T>G、22594C>T。在山羊DNA序列中发现21个短散布元件和14个长散布元件。在绵羊DNA序列中发现19个短散布元件和14个长散布元件。另外,在山羊和绵羊DNA序列中各发现5个简单重复序列(SSR)。
     7605C>T位点上检测到C和T两个等位基因。所研究的10个中国地方山羊品种/品系(承德无角山羊、成都麻羊、济宁青山羊、雷州山羊、辽宁绒山羊、南江黄羊黑色系、南江黄羊高繁系、唐山奶山羊、内蒙古绒山羊、武安山羊)共300个个体。T等位基因频率(50.12%)稍高于C等位基因(49.88%)。从基因型分布看,10个中国地方山羊群体的基因型频率在各群体间差异较大。遗传多样性分析表明,10个中国地方山羊在7605C>T位点具有较丰富的遗传多样性。群体遗传分化分析表明,10个中国地方山羊群体间的平均基因流为2.9433,平均遗传分化指数为0.0783。
     本试验获得了山羊和绵羊LF基因完整编码区(CDS)序列,另外,从GenBank获得了60条不同物种LF基因完整CDS序列,比对分析LF基因在物种间和物种内的遗传变异和分化。大部分物种终止密码子使用TAA,(1个人变异个体使用TAG),而小鼠使用TGA( 1个变异体使用TAA)。物种间长度变异较大(2055bp~2190bp)。根据美国国家生物技术信息中心(NCBI)分类法,相近物种的LF基因的CDS序列长度相同或相差较小。挪威鼠最长为2190bp,是由于终止密码子突变(TAA>GAA)导致的66bp的延长。一个黑猩猩变异体(XM_001150238)LF基因CDS长度最短为2055bp,是由于81bp缺失片段造成的。物种内也存在长度变异,狗长度变异最复杂,其它物种LF基因CDS的变异较小。遗传多样性分析发现多数物种LF基因具有较丰富的遗传多样性,群体间的遗传多样性大于群体内。通过物种内的LF氨基酸序列比对,发现了新的氨基酸变异位点(人8个、小鼠6个、山羊2个、牛10个、猪20个、绵羊3个)。物种间的遗传分化明显,系统发生树所显示的物种的聚类关系与NCBI物种的分类学一致。
On the basis of cattle DNA and goat mRNA lactoferrin (LF) gene sequences from GenBank, partial 5’UTR to exon 17 of goat LF gene fragments was amplified. The total length of determined sequence was 30306bp, which had been handed out to GenBank, the accession number was FJ609300. In order to compareing the genetic diversity and differentiation among ruminant, sheep LF gene was sequenced (30340bp), and the accession number was FJ541507.
     30 point mutations were found in goat LF gene sequence, of which 1 was in exon, and 29 was in intron. The mutation in exon was 7605C>T resulting in Arg136Trp. The 29 mutation in intron were 2428C>G, 2761C>G, 6053A>G, 6538A>C, 7363A>G, 8642C>T, 8655C>T, 9284C>T, 13719C>T, 15783A>G, 15820A>G, 16820A>G, 17774A>T, 17775A>G, 17793A>G, 17808A>G, 17855C>G, 8266C>G, 18428A>G, 18930A>G, 19764C>T, 19997C>T, 20740A>G, 20989C>T, 21146C>T, 21281C>G, 23146C>T, 24300C>T, 24720A>G and 25249A>G. Meanwhile, 16 point mutations were detected in LF gene sequence of sheep, which all were in intron. The 16 mutation in intron were 2432C>T, 2764C>G, 4197C>T, 4274A>T, 4633C>T, 11038C>T, 13045A>G, 13075A>G, 13482A>G, 18072A>G, 20660C>T, 21022C>T, 21625T>G, 21706T>G, 22583T>G and 22594C>T. 21 short interspersed nuclear elements and 14 long interspersed nuclear elements were detected in LF gene of goat, and 19 short interspersed nuclear elements and 14 long interspersed nuclear elements were detected in LF gene of sheep. Moreover, 5 simple sequence repeats (SSR) was discovered in LF gene of goat and sheep, respectively.
     We determined a missense mutation C→T at position 7605. Meanwhile the PCR-RFLP of this mutation was performed in 300 individuals from 10 Chinese indigenous goat breeds/strain(sChengde polled goat, Chengdu Ma goat, Jining gray goat, Leizhou goat, Liaoning cashmere goat, two strains of Nanjiang brown goat, Tangshan dairy goat, Neimenggu cashmere goat and Wuan goat). The results showed that the gene frequency of allele T (50.12%) was little higher and that of allele C (49.88%) and the difference of genotypes was great among goat populations. The 7605C>T site performed higher genetic diversity in goat populations. The genetic differentiation was 0.0783 and gene flow was 2.9433 among 10 goat populations.
     The coding regions (CDS) sequence of LF gene of goat was obtained by splicing in this study, and 60 sequences from 11 species were studied to investigate its evolution and differentiation within and among species. Most species use TAA as stop codon for the LF gene, with only one human variation (DQ892855) using TAG. Mus musculus uses TGA, except for one variation (CT010339) with TAG. It was shown that the length of the LF gene with the complete CDS varies greatly among species, ranging from 2055 to 2190bp. Species with close relationships according to the taxonomy in the National Center for Biotechnology Information (NCBI) have a similar length of the LF gene. The longest length (2190bp) of the LF gene, in Rattus norvegicus, is due to the stop codon mutation of TAA to GAA resulting in the 66 bp elongation. And the shortest length (2055bp) is one variant type of Pan roglodytes (XM_001150238) due to an 81bp deletion. There was length variation within species. Intricate length variations were in Canis familiaris with three variant regions, and simple length variations was in other species. Most species had abundant genetic diversity of LF gene, and genetic diversity was higher among species than within species. Novel amino acid variation sites were detected within several species (8 in Homo sapiens, 6 in Mus musculus, 2 in Capra hircus, 3 in Ovis aries, 10 in Bos Taurus and 20 in Sus scrofa). Differentiation of the LF gene was obvious among species, and the clustering result was consistent with the taxonomy in NCBI.
引文
[1]Montreuil J, Tonnelat J, Mullet S. Preparation and properties of lactotransferrin of human milk[J]. Biochimica et Biophysica Acta, 1960, 45, 413-21.
    [2]Sujatas, Tejp S, Krishanl B. Preparation and characterization of the N and C monoferric lobes of buffalo lactoferrin produced by proteolysis using proteinase K[J]. Journal of Dairy Reasearch, 1999, 66: 81-90.
    [3]Cohen MS, Britigan BE, French M, et al. Preliminary observations on LF sretion in human vaginal mucus:Variation during the menstrual cycle, evidence of hormonal regulation, and implications for infection with Neisseria gonorrhoeae[J]. J Obstet Gynecol, 1987, 157: 1122-5.
    [4]Huang PM. Three dimensional solution structure of lactoferricin B, an antimicrobialpepride derived from bovine lactoferrin[J]. Biochemistry, 1998, 37: 4288-98.
    [5]Caccavo D, Pellegrino N M, Altamura M, et al. Antimicrobial and immunoregulatory functions of lactoferrin and it s potentialt herapeutic application[J]. Endotoxin Res, 2002, 8: 403-17.
    [6]Strom MB, Hang BE, Rekual Q, et al. Important structural features of 15-residue lactoferrcin derivatives and met hods for improvements of antimicrobial activity[J]. Biochem, 2002, 80: 65-74.
    [7]Hunter HN, Demcoe AR, Jenssen H, et al. Human lactoferricin is partially folded in aqueous solution and is better stabilized in membrane mimetic solvent[J]. Antimicrob Agents Chemother, 2005, 49: 3387-95.
    [8]Vogel HJ, Schibli D, Jing W, et al. Towards a structurefunction analysis of bovine lactoferricin and related tryptophan and arginine containing peptides[J]. Biochem Cell Biol, 2002, 80: 49-63.
    [9]Marchetti. Metal complexes of bovine lactoferrin inhibit in vitro resplication of herpes simplex virus[J]. Biometals, 1998, 11:89-94.
    [10]Mattsby B. Lactoferrin or a fragment there of inhibits the endotoxin-induced interleukin-6 response in human mononcytic cell[J]. Pediatr Res, 1996, 40: 257-62.
    [11]Lee WL. The protective effects of lactoderrin feeding against endotoxin letheal shock in germ-free piglet[J]. Infect Immun, 1998, 66(4): 1421-6.
    [12]Farnaud S, EvansRW. Lactoferrin: the conductor of the immunological system[M]. New York: Nova Science Publishers, 2005.34-44.
    [13]Chapple DS, Hussain R, Joannou CL, et al. Structure and association of human lactoferrin peptides with Escherichia coli[J]. Antimicrob Agents Chemother, 2004, 48: 2190-8.
    [14]Andersen J H, Jenssen H, Sandvik K, et al. Anti-HSV activity of lactoferricin is dependent on the presence of heparin sulphate at t he cell surface[J]. Med Virol, 2004, 74: 262-71.
    [15]Jenssen H, Andersen JH, Uhlin HL, et al. Anti-HSV activity of lactoferricin analogues is only partly related to t heir affinity for heparan sulfate[J]. Antiviral Res, 2004, 61: 101-9.
    [16]Mohamed JA, DuPont HL, Jiang ZD. A novel single-nucleotide polymorphism in the LF gene is associated with susceptibility to diarrhea in North American travelers to Mexico[J]. Clin Infect Dis, 2007, 44 (7): 945-52.
    [17]Faber HR, Bland T, Day CL, et al. Altered domain closure and iron binding in transferrins: thecrystal structure of the Asp60Ser mutant of the amino-terminal half-molecule of human lactoferrin[J]. Mol Biol, 1996, 256(2): 352-63.
    [18]Wu HM, Church CF. Arginine 25 and Arginine 28 of lactoferrin are critical for effective heparin neutralization in blood[J]. Arch Biochem Biophys, 2003, 412(1): 121-5.
    [19]Araki K, Ando Y, Nakamura M, et al. Lactoferrin Glu561Asp facilitates sondary amyloidosis in the cornea[J]. Ophthalmol, 2005, 89(6): 684-8.
    [20]李国华,张沅,孙东晓,等.奶牛乳铁蛋白基因5′侧翼区PCR-SSCP多态性分析[J].遗传, 2004, 26(6): 827-30.
    [21]张利军,蔡亚非,刘庆华,等.奶牛乳铁蛋白基因启动子区PCR-RFLP分析与乳房炎的相关性[J].福建农林大学学报, 2005, 34(1): 87-91.
    [22]周磊. bLF基因部分序列的PCR_SSCP分析及其作为乳房炎抗性分子标记的可行性研究[D].南京:南京农业大学硕士学位论文, 2006.
    [23]Kaminski S, Olenski K, Brym P, et al. Single nucleotide polymorphism in the promoter region of the lactoferrin gene and its associations with milk performance traits in Polish Holstein-Friesian cows[J]. Genetika, 2006, 42(8): 1117-20.
    [24]Liang QW, Richardson T. Expression and characterization of human lactoferrin in yeast[J]. J Agric Food Chem, 1993, 41: 1800-7.
    [25]Ward PP, Piddington CS, Cunningham GA, et al. A system for production of commercial quantities of human lactoferrin[J]. Biotechnology, 1995, 13: 498-503.
    [26]Salmon V, Dominique L, Marie CS, et al. Production of human Lactoferrin in transgenic tobacco plants[J]. Protein Expression and Purification, 1998, 13: 127-35.
    [27]Chong DK, Langridge WH. Expression of full-length bioactive antimicrobial human lactoferrin in potato plants[J].Transgentic Res, 2000, 9: 71-8.
    [28]Toman PD, Pieper F, Sakai N, et al. Production of recombinant human type1 procollagen homotrimer in the mammany gland of transgenic mice[J]. Transgentic Res, 1999, 8(6): 4-15.
    [29]Patrick HC, Welling MM, Geeits M, et al. Large scale production of recombinant human lactoferrin in the milk of transgentic cows[J]. Nature Biotechnology, 2002, 20: 484-95.
    [30]萨姆布鲁克J,弗里奇E F,曼尼阿蒂斯T.分子克隆实验指南(第三版)[M].北京:科学技术出版社, 2002.
    [31]吴冠芸,潘华珍.生物化学与分子生物学实验常用数据手册[M].北京:科学出版社, 1999.
    [32]徐华.中国荷斯坦奶牛催乳素基因和微卫星DNA多态性与产奶性能的相关分析[D].保定:河北农业大学硕士学位论文, 2004.
    [33]Swindell SR, Plasterer TN. Seqman Contig assembly[J]. Methods Mol Biol, 1997, 70: 75-89.
    [34]Thompson JD, Higgins DG, Gibson TJ. Clustal W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice[J]. Nucleic Acids Res. 1994, 22(22): 4673-80.
    [35]张成岗,贺福初.生物信息学方法与实践[M].北京:科学出版社, 2002.
    [36]Nei M. Analysis of gene diversity in subdivided populations[J]. Proc Natl Acad Sci, 1973, 70: 3321-3.
    [37]Shannon CE, Weaver W. The mathematical theory of communication[M]. Urbana: Univ. of Illinois Press. 1949.
    [38]Slatkin M, Barton NH. A comparison of three indirect methods for estimating average levels of gene flow[J]. Evolution, 1989, 43: 1349-68.
    [39]Nei M. Estimation of average heterozygosity and genetic distance from a small number of individuals [J]. Genetics, 1978, 89: 583-90.
    [40]Tajima F. Evolutionary relationship of DNA sequences in finite populations[J]. Genetics, 1983, 105: 437-60.
    [41]蔡佩玲,童英,王晓强,等.人乳铁蛋白(hLF) cDNA的克隆及序列分析[J].四川大学学报(自然科学版), 2005, 42(6): 1229-32.
    [42]刘涛,张耀洲,吴祥甫.人乳铁蛋白cDNA的克隆及在家蚕细胞中的表达[J].蚕业科学, 2005, 31(3): 280-5.
    [43]Kabilan V, Jeffrey B, Kaplan, et al. One of two human LF variants exhibits increased antibacterial and transcriptional activation activities and is associated with localized juvenile periodontitis[J]. Infect immune, 2003, 11: 6141-7.
    [44]Lee TH, Shimazaki K, Yu SL, et al. Polymorphic sequence of Korean Native goat lactoferrin exhibiting greater antibacterial activity[J]. Anim Genet, 1997, 28(5):367-9.
    [45]杨光,庞庆,吴长毅,等.人乳铁蛋白基因的克隆及基因序列多态性分析[J].徐州医学院学报, 2005, 25(5): 390-2.
    [46]Lynch M, Crease TJ. The analysis of population survey data on DNA sequence variation[J]. Mol Biol Evol, 1990, 7:377-94.
    [47]Li WH, Dan G. Expand fundamentals of molecular evolution[M]. Sunderland: Sinauer Associates, 1991.
    [48]方平.乳铁蛋白基因多态性与奶牛乳房炎的相关性研究[D].杨凌,西北农林科技大学硕士论文, 2007.
    [49]Hoek KS, Milne JM, Grieve PA, et al. Antibacterial activity of bovine lactoferrin-derived peptides[J]. Antimicrob, 1997, 41: 54-9.
    [50]Kang JF, Li XL, Zhou RY, et al. Bioinformatics analysis of lactoferrin gene for several species[J]. Biochem Genet, 2008, 46: 312-22.
    [51]Kimura M, Nam M S, Ohkouchi Y, et al. Antimicrobial peptide of korean native goat lactoferrin and identification of the part essential for this activity, in Biochem Biophys[J]. Res Commun, 2000, 3: 33-6.
    [52]狄冉,何晓红,韩建林,等.中国绒山羊遗传多样性现状和系统发生关系的微卫星分析[J].生物多样性, 2007, 15 (5): 470-8.
    [53]Canon J, Garcia D, Garcia-Atance MA, et al. Geographical partitioning of goat diversity in Europe and the Middle East[J]. Animal Genetics, 2006, 37: 327-334.
    [54]Saitbekova N, Gaillard C, Obexer-Ruff G, et al. Genetic diversity in Swiss goat breeds based on microsatellite analysis[J]. Animal Genetics, 1999, 30: 36-41.
    [55]Aggarwal RK, Dixit SP, Verma NK, et al. Population genetics analysis of Mehsana goat based onmicrosatellite markers[J]. Current Science, 2007, 92: 1133-7.
    [56]李祥龙,张亚平,陈圣偶,等.我国主要地方山羊品种随机扩增多态DNA研究[J].畜牧兽医学报,2000, 31 (5): 416-22.
    [57]Li XL, Valentini A. Genetic diversity of Chinese indigenous goat breeds based on microsatellite markers [J]. Journal of Animal Breeding and Genetics, 2004, 121(5): 350-5.
    [58]Tang CJ, Zhou RY, Li XL, et al. Variation of 423G>T in the agouti gene exon 4 in indigenous Chinese goat breeds[J]. Biochem Genet, 2008, 46:770-80.
    [59]王建涛.山羊毛色稀释位点候选基因RAB27A基因研究[D].保定:河北农业大学硕士学位论文, 2008.
    [60]冯付军.山羊毛色稀释位点候选基因MLPH研究[D].保定:河北农业大学硕士学位论文, 2008.
    [61]郑向忠,徐宏发,陆厚基.动物种群遗传异质性研究进展[J].生物多样性, 1997, 5(3): 210-6.
    [62]Wright S. Evolution in Mendelian population [J]. Genetics, 1931, 16: 91-159.
    [63]Li XL, Zheng GR, Zhou RY, et al. Evolution and Differentiation of MSHR Gene in Different Species[J]. Journal of Heredity, 2007, 98(2): 165-8.
    [64]Ghosh T. Studies on codon usage in Entamoeba histolytica[J]. Int J Parasitol, 2000, 30: 715-22.
    [65]Robert FW. Molecular biology[M]. 2001, University of Kansas-Lawrence, Scientific Press.
    [66]Yang Z, Yoder AD. Comparison of likelihood and Bayesian methods for estimating divergence times using multiple gene Loci and calibration points, with application to a radiation of cute-looking mouse lemur species[J]. Syst Biol, 2003, 52(5): 705-16.
    [67]Wildman DE, Uddin M, Liu G, et al. Implications of natural selection in shaping 99.4% nonsynonymous DNA identity between humans and chimpanzees: enlarging genus Homo[J]. Proc Natl Acad Sci, 2003, 100(12): 7181-8.
    [68]Kumar S, Hedges SB. A molecular timescale for vertebrate evolution[J]. Nature, 1998, 392(6679): 917-9.
    [69]Arnason U, Gullberg A, Gretarsdottir S, et al. The mitochondrial genome of the sperm whale and a new molecular reference for estimating eutherian divergence dates[J]. Mol Evol, 2000, 50(6): 569-78.
    [70]Nikaido M, Kawai K, Cao Y, et al. Maximum likelihood analysis of the complete mitochondrial genomes of eutherians and a reevaluation of the phylogeny of bats and instivores[J]. Mol Evol, 2001, 53(4): 508-16.
    [71]杨银凤,唐博,曹贵方.骆驼β-防御素caBD-1cDNA的克隆及序列分析[J].畜牧兽医学报, 2004, 35(4): 357-61.
    [72]唐博,曹贵方,杨银凤,等.骆驼β-防御素caBD-1 cDNA基因全序列RACE扩增方法的建立[J].中国兽医科学, 2006, 36(2): 151-6.
    [73]Delsuc F, Vizca?′SF, Douzery EJ. Influence of tertiary paleoenvironmental changes on the diversification of South American mammals: a relaxed molecular clock study within xenarthrans [J/OL]. BMC Evol Biol, 2004, 4:11
    [74]张蓉,刁其玉.乳铁蛋白在犊牛生产中的应用及其生物信息学分析[J].乳业科学与技术, 2007, 1: 39-42.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700