Identification of genome diversity in different breeds of Iranian native sheep using the whole genome sequencing method
Purpose: Iran is considered to be one of the oldest centers of domestication and breeding of livestock and poultry species in the world. Currently, different ecotypes of indigenous sheep are kept in different geographical regions of the country, which have obvious differences from each other in term...
Main Authors: | , , , |
---|---|
Format: | Article |
Language: | fas |
Published: |
Shahid Bahonar University of Kerman
2023-11-01
|
Series: | مجله بیوتکنولوژی کشاورزی |
Subjects: | |
Online Access: | https://jab.uk.ac.ir/article_4028_5862a5a5cc5156bc77e562811f9afd5a.pdf |
_version_ | 1797388985892339712 |
---|---|
author | Zeinab Amiri Ghanatsaman Masoud Asadi Fouzi Hojjat Asadollahpour nanaei Mohammad hossein Banabazi |
author_facet | Zeinab Amiri Ghanatsaman Masoud Asadi Fouzi Hojjat Asadollahpour nanaei Mohammad hossein Banabazi |
author_sort | Zeinab Amiri Ghanatsaman |
collection | DOAJ |
description | Purpose: Iran is considered to be one of the oldest centers of domestication and breeding of livestock and poultry species in the world. Currently, different ecotypes of indigenous sheep are kept in different geographical regions of the country, which have obvious differences from each other in terms of appearance and production characteristics. So far, there has not been a comprehensive study on the whole genome level to identify the genetic diversity of native Iranian sheep. Therefore, the aim of this study is to identify the genomic characteristics of these native reserves in order to organize appropriate programs for their exploitation and protectionMaterials and methodsIn this study, the whole genome sequences of 29 native Iranian sheep were downloaded from the NCBI database and analyzed. Whole genome sequencing of the studied data has been done by Hiseq2000 and Hiseq X Ten sequencer devices. Quality control of raw data sequences was done by FastQC program. To align the sequence data with the sheep reference genome (Oar v.4.0, https://www.ncbi.nlm.nih.gov/assembly/GCF_000298735.2), the BWA-MEM algorithm used in the BWA software package was used. Picard program was used to remove PCR duplicates from mapping outputs. The alignment outputs with the reference genome were processed in two stages, including re-alignment of deletions and small insertions and recalibration of the base quality score using the GATK program. The average coverage depth and alignment percentage for alignment output with the reference genome were calculated using depth and flagstat commands used in samtools software. Single nucleotide polymorphisms (SNPs) were identified by the UnifiedGenotyper tool used in the GATK program. Nucleotide diversity values and genomic inbreeding coefficient were calculated based on homozygous SNPs for each individual using the het command used in the VCFtools program.ResultsThe average coverage depth of the used data in this study was 18.39 X. The average percentage of alignment of short sequences with the sheep reference genome was 99.89%. The values of genomic inbreeding coefficient in Iranian native sheep breeds ranged from 0.01 to 0.12. The lowest value of genomic inbreeding coefficient was observed in the genome of Mughani sheep (0.01) and the highest value of inbreeding coefficient was observed in the genome of Afshari sheep. Also, the average values of observed and expected percentage of heterozygosity calculated for single nucleotide polymorphisms in the genome of Iranian native sheep ecotypes ranged from 20.67 to 23.06 and 32.415 to 32.421. |
first_indexed | 2024-03-08T22:48:46Z |
format | Article |
id | doaj.art-01318d6c614040b0a116325b3416e29e |
institution | Directory Open Access Journal |
issn | 2228-6705 2228-6500 |
language | fas |
last_indexed | 2024-03-08T22:48:46Z |
publishDate | 2023-11-01 |
publisher | Shahid Bahonar University of Kerman |
record_format | Article |
series | مجله بیوتکنولوژی کشاورزی |
spelling | doaj.art-01318d6c614040b0a116325b3416e29e2023-12-16T17:42:52ZfasShahid Bahonar University of Kermanمجله بیوتکنولوژی کشاورزی2228-67052228-65002023-11-0115414516010.22103/jab.2023.22018.15024028Identification of genome diversity in different breeds of Iranian native sheep using the whole genome sequencing methodZeinab Amiri Ghanatsaman0Masoud Asadi Fouzi1Hojjat Asadollahpour nanaei2Mohammad hossein Banabazi3Animal Science Research Department, Fars Agricultural and Natural Resources, Research and Education Center, Agricultural Research, Education and Extension Organization (AREEO), Shiraz, IranDepartment of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, Kerman, Iran.Postdoctoral Researcher, Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, Kerman, IranAgricultural Research, Education and Extension Organization (AREEO), Animal Science Research Institute of Iran, Karaj, IranPurpose: Iran is considered to be one of the oldest centers of domestication and breeding of livestock and poultry species in the world. Currently, different ecotypes of indigenous sheep are kept in different geographical regions of the country, which have obvious differences from each other in terms of appearance and production characteristics. So far, there has not been a comprehensive study on the whole genome level to identify the genetic diversity of native Iranian sheep. Therefore, the aim of this study is to identify the genomic characteristics of these native reserves in order to organize appropriate programs for their exploitation and protectionMaterials and methodsIn this study, the whole genome sequences of 29 native Iranian sheep were downloaded from the NCBI database and analyzed. Whole genome sequencing of the studied data has been done by Hiseq2000 and Hiseq X Ten sequencer devices. Quality control of raw data sequences was done by FastQC program. To align the sequence data with the sheep reference genome (Oar v.4.0, https://www.ncbi.nlm.nih.gov/assembly/GCF_000298735.2), the BWA-MEM algorithm used in the BWA software package was used. Picard program was used to remove PCR duplicates from mapping outputs. The alignment outputs with the reference genome were processed in two stages, including re-alignment of deletions and small insertions and recalibration of the base quality score using the GATK program. The average coverage depth and alignment percentage for alignment output with the reference genome were calculated using depth and flagstat commands used in samtools software. Single nucleotide polymorphisms (SNPs) were identified by the UnifiedGenotyper tool used in the GATK program. Nucleotide diversity values and genomic inbreeding coefficient were calculated based on homozygous SNPs for each individual using the het command used in the VCFtools program.ResultsThe average coverage depth of the used data in this study was 18.39 X. The average percentage of alignment of short sequences with the sheep reference genome was 99.89%. The values of genomic inbreeding coefficient in Iranian native sheep breeds ranged from 0.01 to 0.12. The lowest value of genomic inbreeding coefficient was observed in the genome of Mughani sheep (0.01) and the highest value of inbreeding coefficient was observed in the genome of Afshari sheep. Also, the average values of observed and expected percentage of heterozygosity calculated for single nucleotide polymorphisms in the genome of Iranian native sheep ecotypes ranged from 20.67 to 23.06 and 32.415 to 32.421.https://jab.uk.ac.ir/article_4028_5862a5a5cc5156bc77e562811f9afd5a.pdfiranian native sheepwhole genome sequencingsingle nucleotide polymorphisms |
spellingShingle | Zeinab Amiri Ghanatsaman Masoud Asadi Fouzi Hojjat Asadollahpour nanaei Mohammad hossein Banabazi Identification of genome diversity in different breeds of Iranian native sheep using the whole genome sequencing method مجله بیوتکنولوژی کشاورزی iranian native sheep whole genome sequencing single nucleotide polymorphisms |
title | Identification of genome diversity in different breeds of Iranian native sheep using the whole genome sequencing method |
title_full | Identification of genome diversity in different breeds of Iranian native sheep using the whole genome sequencing method |
title_fullStr | Identification of genome diversity in different breeds of Iranian native sheep using the whole genome sequencing method |
title_full_unstemmed | Identification of genome diversity in different breeds of Iranian native sheep using the whole genome sequencing method |
title_short | Identification of genome diversity in different breeds of Iranian native sheep using the whole genome sequencing method |
title_sort | identification of genome diversity in different breeds of iranian native sheep using the whole genome sequencing method |
topic | iranian native sheep whole genome sequencing single nucleotide polymorphisms |
url | https://jab.uk.ac.ir/article_4028_5862a5a5cc5156bc77e562811f9afd5a.pdf |
work_keys_str_mv | AT zeinabamirighanatsaman identificationofgenomediversityindifferentbreedsofiraniannativesheepusingthewholegenomesequencingmethod AT masoudasadifouzi identificationofgenomediversityindifferentbreedsofiraniannativesheepusingthewholegenomesequencingmethod AT hojjatasadollahpournanaei identificationofgenomediversityindifferentbreedsofiraniannativesheepusingthewholegenomesequencingmethod AT mohammadhosseinbanabazi identificationofgenomediversityindifferentbreedsofiraniannativesheepusingthewholegenomesequencingmethod |