Volume: 54 Issue: 2
Year: 2023, Page: 524-531, Doi: https://doi.org/10.51966/jvas.2023.54.2.524-531
Received: Nov. 15, 2022 Accepted: Jan. 31, 2023 Published: June 30, 2023
Ankamali pig is a domesticated native variety of Kerala which is well known for its disease resistance, lean meat and adaptability to hot tropical environments. Recent breakthrough in genome sequencing technologies have created unparalleled prospects to characterize individual genomic landscapes and identifying mutations between and within populations. The current study aims to determine the genetic variations in Ankamali pigs using whole genome sequencing. The GATK HaplotypeCaller was used to identify the variants. There were over 26 million (26,604,589) single nucleotide variants (SNVs), including more than 21 million SNPs and over 5 million indels. In Ankamali pigs, the total genome length obtained was more than 2.5 billion with an average variant rate of one variant in every 94 bases. The significance of different variant rate on 18 chromosomes were analysed using the chi-square statistics. The variant rates in Sus scrofa chromosomes10 and 13 were significantly different (p<0.01%) in Ankamali pigs. The significantly higher variable rate on chromosome 10 was observed with one variant per 64 bases. Whereas, significantly lower variable rate was observed on chromosome 13, with one variant in every 122 bases. The variant rate reported on Sus scrofa chromosome 12 (SSC12) was also significantly higher (p<0.05%), having one variant per 72 bases. The variability of many genes and QTLs associated with several haematological traits and meat quality traits located on these chromosomes may contribute the phenotypic and genetic uniqueness of Ankamali animals
Keywords: GATK Haplotype caller, single nucleotide variants, indels, haematological traits, missense, quantitative trait loci (QTL)
Baes, C.F., Dolezal, M.A., Koltes, J.E., Bapst, B., Fritz-Waters, E., Jansen, S., Flury, C., Signer-Hasler, H., Stricker, C., Fernando, R. and Fries, R. 2014. Evaluation of variant identification methods for whole genome sequencing data in dairy cattle. BMC genomics. 15: 1-18.
Bainbridge, M.N., Wang, M., Wu, Y., Newsham, I., Muzny, D.M., Jefferies, J.L., Albert, T.J., Burgess, D.L. and Gibbs, R.A. 2011. Targeted enrichment beyond the consensus coding DNA sequence exome reveals exons with higher variant densities. Gen. biol. 12: 1-12.
Bartonicek, N., Clark, M.B., Quek, X.C., Torpy, J.R., Pritchard, A.L., Maag, J.L.V., Gloss, B.S., Crawford, J., Taft, R.J., Hayward, N.K. and Montgomery, G.W. 2017. Intergenic disease-associated regions are abundant in novel transcripts. Genome Biol. 18: 1-16.
Behl, R., Sheoran, N., Behl, J. and Vijh, R.K. 2006. Genetic analysis of Ankamali pigs of India using microsatellite markers and their comparison with other domesticated Indian pig types. J. Anim. Breed. Genet. 123: 131-135.
Bickhart, D.M. and Liu, G.E. 2014. The challenges and importance of structural variation detection in livestock. Front. Genet. 5: 37.
Bolger, A.M., Lohse, M. and Usadel, B. 2014. Trimmomatic: a flexible trimmer for Illumina sequence data. J. Bioinform. 30: 2114-2120.
Cingolani, P., Platts, A., Wang, L.L., Coon, M., Nguyen, T., Wang, L., Land, S.J., Lu, X. and Ruden, D.M. 2012. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff SNPs in the genome of Drosophila melanogaster strain. Fly. 6: 80-92.
Cock, P.J.A., Fields, C.J., Goto, N., Heuer, M.L. and Rice, P.M. 2009. The Sanger FASTQ file format for sequences with quality scores, and the Solexa/Illumina FASTQ variants. Nucleic Acids Res. 38: 1767-1771.
D'Alessandro, E., Sapienza, I., Giosa, D., Giuffrè, L. and Zumbo, A. 2019. In silico analysis of meat quality candidate genes among Nero Siciliano, and Italian heavy pigs genomes. Large Anim. Rev. 25: 137-140.
DePristo, M.A., Banks, E., Poplin, R., Garimella, K.V., Maguire, J.R., Hartl, C., Philippakis, A.A., Del Angel, G., Rivas, M.A., Hanna, M. and McKenna, A. 2011. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 43: 491-498.
Guo, Y., Ye, F., Sheng, Q., Clark, T. and Samuels, D.C. 2013. Three-stage quality control strategies for DNA re-sequencing data. Brief. Bioinform. 15(6): 879-89.
Lefouili, M. and Nam, K., 2022. The evaluation of Bcftools mpileup and GATK HaplotypeCaller for variant calling in non-human species. Sci. Rep. 12:1-8.
Li, H. and Durbin, R. 2009. Fast and accurate short read alignment with Burrows–Wheeler transform. J. Bioinform. 25: 1754-1760.
Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., Marth, G., Abecasis, G. and Durbin, R. 2009. The sequence alignment/map format and SAMtools. J. Bioinform. 25: 2078-2079.
Marykutty, T., Radhika, G., Aravindakshan, T.V., Thirupathy, R., Raji, K. and Shynu, M. 2021. Linkage disequilibrium over short physical genomic distances measured using medium density SNP beadchip in native goat breeds of Kerala. J. Vet. Anim. Sci. 52: 14-18.
McKenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., Garimella, K., Altshuler, D., Gabriel, S., Daly, M. and DePristo, M.A. 2010. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20: 1297-1303.
Mei, C., Wang, H., Liao, Q., Wang, L., Cheng, G., Wang, H., Zhao, C., Zhao, S., Song, J., Guang, X. and Liu, G.E. 2018. Genetic architecture and selection of Chinese cattle revealed by whole genome resequencing. Mol. Biol. Evol. 35: 688-699.
Park, S.T. and Kim, J. 2016. Trends in next-generation sequencing and a new era for whole genome sequencing. Int. Neurol. J. 20: S76.
Ren, S., Ahmed, N., Bertels, K. and Al-Ars, Z. 2019. GPU accelerated sequence alignment with traceback for GATK HaplotypeCaller. BMC Genom. 20: 103-116.
Reshma, R.S., Aravindakshan, T.V., Radhika, G., Naicy, T. and Raji, K. 2020. Exome wide variant discovery by next generation DNA sequencing in Vechur cattle of Kerala. J. Vet. Anim. Sci. 51: 201-206.
Rohrer, G.A. 2000. Identification of quantitative trait loci affecting birth characters and accumulation of backfat and weight in a Meishan-White Composite resource population. J. Anim. Sci. 78: 2547-2553.
Snedecor, G.W. and Cochran, W.G. 1989. Statistical Methods. (8th Ed.). Iowa State University Press, Ames, Iowa.
Yu, B.T., Wang, L., Tuggle, C.K. and Rothschild, M.F. 1999. Mapping genes for fatness and growth on pig chromosome 13: a search in the region close to the pig PIT1 gene. J. Anim. Breed. Genet. 116:269-280.
Yu, J., Zhao, P., Zheng, X., Zhou, L., Wang, C. and Liu, J.F. 2020. Genome-wide detection of selection signatures in Duroc revealed candidate genes relating to growth and meat quality. G3: Genes Genom. Genet. 10: 3765-3773
Zhang, W., Yang, M., Zhou, M., Wang, Y., Wu, X., Zhang, X., Ding, Y., Zhao, G., Yin, Z. and Wang, C. 2020. Identification of signatures of selection by whole-genome resequencing of a Chinese native pig. Front. Genet. 11: 566255.
© 2023 Michelle et al. This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Michelle, E.R., Manoj, M., Rojan, P.M., Tina, S., Aravindakshan, T.V., Usha, A.P. and mUnnikrishnan, M.P. 2023. Identification of genetic variants by whole genome sequencing inm Ankamali pigs of Kerala. J. Vet. Anim. Sci. 54(2):524-531
DOI: https://doi.org/10.51966/jvas.2023.54.2.524-531