Livestock Research for Rural Development 28 (3) 2016 Guide for preparation of papers LRRD Newsletter

Citation of this paper

Molecular genetic characterization of Kivircik sheep breed raised in Western Anatolia

Onur Yilmaz, İbrahim Cemal, Orhan Karaca and Nezih Ata

Adnan Menderes University Faculty of Agriculture, Department of Animal Science, Aydin, Turkey.
o-yilmaz@live.com

Abstract

Genetic characterization and diversity of Kıvırcık sheep breed (n=246) reared in three different locations (Aydın, Bandırma and Uşak) were investigated by 21 microsatellite markers recommended by FAO.

The microsatellite markers used in this study showed high levels of polymorphism. A total of 461 alleles were detected in Kıvırcık populations. The mean values of polymorphic information content (PIC=0.84), observed heterozygosity (Ho=0.81) and expected heterozygosity (He=0.85) proved that Kıvırcık populations possess remarkable genetic variability. The mean number of alleles ranged from 11.67 (Bandırma) to 17.05 (Eşme). The results demonstrated that the microsatellite markers that were used were adequately polymorphic, and these markers can be successfully used to investigate genetic diversity in these three populations. Our results showed that within-breed diversity was higher than the between breed diversity. This situation can be seen as a chance in terms of the breeding programs and genetic conservation programs for these breeds.

Key words: genetic diversity, microsatellite, population structure


Introduction

The past century has seen the widespread application of statistical methods to the study of genetic inheritance, especially in farm animals. Recent developments in molecular biology and statistics have opened the possibility of identifying and using genomic variation and major genes for the genetic improvement of livestock. Molecular genetic techniques to identify the genetic structure and diversity in farm animals have shown rapid development in recent years and began to be widely used. Various molecular genetics techniques have been developed to achieve genetic structure and diversity. Microsatellites, which are highly polymorphic and abundant, often found in non-coding regions of the genes, are widely used for paternity and genetic diversity studies.

Microsatellites, which are valuable genetic markers due to their dense distribution in the genome, high variation, co-dominant inheritance, easy genotyping and scored, are used rather widely for the definition of intra and interbreed genetic diversity (Togan et al 2005; Sancristobal et al 2003; Schlötterer 2004; Agaviezor et al 2012; Alvarez et al 2012; Arora et al 2011; Cemal et al 2013; Jyotsana et al 2010; Kusza et al 2011; Lasagna et al 2011; Yilmaz and Karaca 2012; Yilmaz et al 2013, Yilmaz et al 2014).

Turkey hosts a high number of breeds in terms of farm animal species and a high level of diversity thanks to its geographical location as a passage among continents and as a center of domestication of domestic animals such as cattle, sheep and buffalo. Kıvırcık sheep breed is the most important genetic resource for the meat production in Turkey. In addition, this breed has contributed to the formation of many types of sheep such as Türkgeldi, Karacabey Merino and Tahirova sheep, are raised in southern and western provinces in the Marmara and Aegean region as well as Thrace.

This study aims to determine the genetic similarities, differences and the genetic structure of the Kıvırcık sheep breed raised different region in Turkey using microsatellites.


Materials and Methods

Animal material of the study consisted of a total of 246 animals raised at Aydın, Bandırma and Usak (Esme) province.

DNA was extracted from 50 ng of whole blood using the ABM™ Genomic DNA Kit (Applied Biological Materials Inc., Canada) according to the manufacturer’s recommendations. DNA quantity and quality were assessed using a NanoDrop 2000 spectrophotometer (Thermo Scientific, USA).

Three multiplex groups formed according to fragment size of 21 microsatellite loci selected from the FAO recommended list (2011). Genomic DNA was amplified with these multiplex groups by the Polymerase Chain Reaction (PCR) in accordance with the touchdown PCR technique (Table 1).

Table 1. Thermal cycling conditions according to Touchdown PCR

Locus
(Florescent Dye)

Multiplex Group

First Denaturation

Denaturation

Annealing

Extension

Cycle

Final Extension

BM1818 (D4)

1

95 ºC

(5 min)

95 ºC

(40sec)

63-54 ºC

(40 Sec)

72 ºC

(60 sec)

40

72 ºC

(10 min)

D5S2 (D4)

INRA0132 (D4)

INRA0023 (D3)

OarAE0129 (D2)

OarCP34 (D4)

OarFCB193 (D3)

OarFCB20 (D2)

OarFCB304 (D3)

BM8125 (D3)

2

95 ºC

(5 min)

95 ºC

(40sec)

60-50 ºC

(40 sec)

72 ºC

(60 sec)

34

72 ºC

(10 min)

CSRD0247 (D3)

HSC (D2)

BM1329 (D2)

MAF214 (D4)

McM0527 (D3)

OarFCB128 (D2)

OarJMP29 (D4)

INRA005 (D2)

3

95 ºC

(5 min)

95 ºC

(40sec)

63-50 ºC

(40 sec)

72 ºC

(60 sec)

42

72 ºC

(10 min)

OARFCB0011 (D3)

DYMS1 (D3)

MAF0065 (D4)

Each 25 μl PCR reaction mixture contained dNTP (0.2 mM), MgCl2 (2.0 mM), primers (0.10 μM), 5X PCR buffer, Taq DNA polymerase (1U), 50 ng genomic DNA and nuclease free water. PCR fragments were separated by capillary electrophoresis in the

Analyses of fragments were performed using an automated DNA sequencer (GenomeLab™ GeXP Genetic Analysis System, Beckman Coulter, Inc., USA) and Beckman Coulter GeXP Fragment Analysis Software (Beckman Coulter, Inc., USA)

The polymorphism statistics, Hardy-Weinberg equilibrium and F statistics (Weir and Cockerham 1984) were estimated using GenAlEx (Peakall and Smouse 2006), Fstat ver.2.9.3 (Goudet 2001), POPGENE (Yeh et al 1997) and MEGA 4 (Tamura et al 2007). A bootstrap-supported (1,000 replications) dendrogram was constructed according to Nei's (1972) minimum distances.

The population structures were analyzed by cluster techniques based on the Bayesian approach, using the STRUCTURE 2.1 software ( Pritchard et al 2000). The burn-in and MCMC (Markov Chain Monte Carlo) lengths were 20,000 and 100,000, respectively. The analyses were realized at different K values (2-4). The most appropriate cluster number (cluster-K) was detected using the method (ΔK = m|L''(K)|/s[L(K)]) reported by Evanno et al (2005). The true K was determined using Structure Harvester Web version 0.6.93 (Earl 2012).


Results and Discussion

A total of 461 alleles were detected across the 21 loci investigated (Table 2).

Table 2. Polymorphism statistics of microsatellites and Wright’s F-statistics

Locus

N

Na

MNa

Ne

PIC

Ho

He

FIS

FIT

FST

DST

GST

HT

HWE

OARFCB304

239

24

18.0

6.0

0.82

0.78

0.83

0.067

0.097

0.032

0.023

0.027

0.847

ns

OARFCB193

236

21

17.7

7.5

0.86

0.95

0.87

-0.104

-0.077

0.024

0.018

0.020

0.876

***

INRA0023

213

20

13.0

4.1

0.74

0.45

0.76

0.320

0.401

0.119

0.087

0.113

0.772

***

OARFCB20

207

37

17.7

12.3

0.91

0.83

0.92

0.073

0.147

0.080

0.069

0.074

0.929

ns

OARAE0129

222

22

13.0

4.5

0.75

0.6

0.78

0.221

0.245

0.031

0.019

0.025

0.768

***

BM1818

229

27

19.0

9.1

0.88

0.95

0.89

-0.127

-0.066

0.054

0.045

0.050

0.899

***

INRA0132

238

22

17.0

11.4

0.91

0.81

0.91

0.083

0.123

0.043

0.035

0.039

0.909

ns

OARCP034

238

14

12.0

6.1

0.82

0.9

0.84

-0.151

-0.088

0.055

0.043

0.051

0.840

ns

D5S2

237

16

11.7

5.3

0.79

0.64

0.81

0.184

0.241

0.071

0.053

0.065

0.814

***

CSRD247

242

25

18.0

9.1

0.88

0.91

0.89

-0.100

-0.034

0.060

0.050

0.057

0.883

ns

MCM527

242

16

12.0

7.7

0.86

0.78

0.87

0.077

0.128

0.055

0.044

0.051

0.872

ns

BM8125

242

20

11.7

6.9

0.84

0.84

0.86

-0.020

0.005

0.025

0.018

0.021

0.850

**

HSC

207

21

14.3

6.2

0.82

0.72

0.84

0.101

0.128

0.029

0.020

0.024

0.839

***

BM1329

209

23

11.3

6.9

0.84

0.8

0.86

-0.041

0.058

0.095

0.076

0.090

0.841

ns

OARFCB128

238

19

12.7

7.7

0.86

0.82

0.87

0.031

0.052

0.021

0.014

0.017

0.868

ns

OARJMP29

243

30

20.0

9.6

0.89

0.82

0.9

0.043

0.072

0.030

0.023

0.025

0.896

ns

MAF214

242

30

16.3

5.4

0.80

0.89

0.82

-0.147

-0.112

0.030

0.021

0.026

0.810

***

INRA005

233

25

15.3

6.4

0.83

0.86

0.85

-0.072

-0.051

0.020

0.013

0.016

0.836

*

OARFCB11

238

15

12.3

7.9

0.86

0.87

0.87

-0.007

0.015

0.021

0.015

0.017

0.875

ns

DYMS1

238

22

13.0

8.4

0.87

0.93

0.88

-0.127

-0.065

0.054

0.045

0.050

0.886

***

MAF0065

239

12

9.3

5.3

0.79

0.94

0.81

-0.170

-0.158

0.010

0.005

0.006

0.807

***

General

14.5

7.3

0.84

0.81

0.85

0.003

0.048

0.046

0.035

0.041

0.853

Na=number of alleles, MNa=mean number of alleles, Ne=effective number of alleles, PIC= polymorphic information content, Ho= observed heterozygosity, He= expected heterozygosity, FIT, FIS, FST =Wright’s F-statistics, DST = the diversity between breeds, G ST =coefficient of gene differentiation for each locus, HT =Nei’s gene diversity, HWE= Hardy-Weinberg equilibrium, ns: non-significant, ***: P<0.001, **:P<0.01, *: P<0.05

The polymorphism statistics (Na, MNa, Ne, PIC, Ho, and He) obtained from this study were higher than previous research (Grigaliunaite et al 2003; Arora et al 2011; Koban 2004; Stahlberger-Saitbekova et al 2001; Baumung et al 2006; Jyotsana et al 2010; Agaveizor et al 2012; Yılmaz and Karaca 2012; Cemal et al 2013; Yılmaz et al 2013; Pariset et al 2003; Koban 2004; Cerit et al 2004; Mukesh et al 2006; Acar 2010). This situation is an indication that studied microsatellites have quite high polymorphic and genetic variability. It was reported that required He value was range between 0.30 and 0.80 to be used in the genetic diversity studies (Takezaki and Nei 1996). It can be said that used loci enables power of separation even though individuals are close to each other.

It is noteworthy that there is loss heterozygosity in eleven microsatellites when examined obtained FIS values. Similar findings have been reported in the literature (Paiva et al 2005; Alvarez et al 2012). Although, obtained general mean of FST value was higher than those of earlier studies (Cemal et al 2013; Hoda and Marsan, 2012), it was lower than some of these (Mukesh et al 2006; Arora et al 2011).

GST value indicates that 4.1% of total genetic variation resulted from the differences between populations; while 95.9% can be explained by the difference between the individuals. General mean of GST was lower than the values reported by Agaveizor et al (2012). Obtained DST values indicate a low genetic diversity among populations. These findings are an expecting result because this research performed in the same breed raised in a different location. General mean of gene diversity between populations (DST) obtained in the present study was lower than the previous literature (Hoda and Marsan 2012). Nei average gene diversity (HT) for loci ranged from 0.768 to 0.929, with an average value of 0.853. Obtained results from this study in terms of Nei average gene diversity confirmed that these markers were appropriate for measuring genetic variation (Takezaki and Nei 1996). Obtained HT values were higher than values reported by Ligda et al (2009) and Hoda and Marsan (2012). It is an important indication that the relatively high genetic diversity in studied populations. Eleven loci revealed significant (P < 0.05) deviations from Hardy–Weinberg equilibrium. These results are encountered as a natural result of the selection programs that have been carried out for several years in these populations.

The results of population based evaluation are summarized in Table 3.

Table 3. Polymorphism statistics, FIS value, number of loci in the Hardy- Weinberg disequilibrium (P < 0.05) and number of private alleles of Kıvırcık populations

Population MNa

Heterozygosity

FIS

Private allele

Ho±SE

He±SE

HWE

NPA
(>%5)

NPA
(<%5)

TNPA

KIV-E

17.1

0.87±0.208

0.85±0.080

-0.027

16.00

19

86

105

KIV-B

11.7

0.76±0.241

0.79±0.072

0.060

17.00

-

24

24

KIV-A

14.0

0.76±0.190

0.80±0.095

0.051

19.00

4

32

36

MNa: Mean number of allele, Ho= observed heterozygosity, He= expected heterozygosity NPA: Number of private allele, TNPA:Total number of private allele

Obtained results from the studied populations raised in three different location in terms of number of allele is an indicator of intra-breed diversity. The mean number of allele and He values were higher than the numerous studies implemented in different sheep breeds (Alvarez et al 2012; Qurioz et al 2008; Baumung et al 2006; Grigaliunaite et al 2013; Paiva et al 2005; Pariset et al 2003; Yilmaz et al 2013).

FIS values, indicated the loss of heterozygosity and that is an important parameter in describing the population characteristics, changed between -0.027 (KIV-E) and 0.060 (KIV-A) in the present study. FIS values was positive in KIV-B and KIV-A population that indicates general risk of inbreeding. It was found that 16 loci in the KIV-E; 17 loci in the KIV-B and 19 loci in the KIV-A population showed significant deviations from Hardy-Weinberg equilibrium. A total of 165 private alleles were observed in three population. However, only 23 of them in KIV-E and KIV-A population have a frequency higher than 5%. The determined 23 (>5%) private alleles can be said to have the property of being able to determine the populations.

When the dendrogram (Figure 1) was examined, it was found that the KIV-E population was in a separate group from KIV-A and KIV-B population.

Figure 1. Dendrogram based on Nei’s minimum genetic distances of three different Kıvırcık populations
(bootstrap resampling methodology (1000 replicates)) (KIV-E: Esme, KIV-B: Bandırma, KIV-A: Aydın).

The results of STRUCTURE analysis containing different numbers of clustering are given in Figure 2.

Figure 2. Estimation of the population structure with different K (K=2, 3 and 4)
values (KIV-E: Eşme, KIV-B: Bandırma, KIV-A: Aydın).

The results obtained from structure analysis have been similar to the dendrogram given in Figure 1. STRUCTURE analysis showed that KIV-E population presents quite complicated population structure for loci investigated. At the same time the results obtained by STRUCTURE analysis concordance with the dendrograms genetic similarities between Kıvırcık populations.

For the purpose of presenting the suitable cluster number (K) in structure analysis, the statistical analysis results received with the method which has been stated by Evanno et al (2005) was given in Table 4.

Table 4. The estimation of posterior probabilities [Ln Pr(X|K)] for different numbers of inferred clusters (K) and ΔK

K

[Ln Pr(X|K)]

ΔK

1

-22851.2

-

2

-21554.3

263.3973

3

-20933.3

108.1885

4

-20752.7

0.951034

ΔK value, which has been taken from 3 studied Kıvırcık population, shows that the most suitable group number is 2 (K=2). It is noteworthy that the mixed population structure in Kıvırcık according to investigated all loci when the Table 5 examined. This situation are an expecting result because this research performed in the same breed raised in a different location.


Conclusions


Acknowledgements

We acknowledge TÜBİTAK, who supplied animal materials for the study, (Project No: KAMAG-109G017), General Directorate for Agricultural Research and Policies (GDAR), Bandırma Sheep Research Station and Adnan Menderes University Agricultural Biotechnology and Food Safety Application and Research Center (ADÜ-TARBİYOMER) providing laboratory facilities for molecular genetics analyses.


References

Agaviezor B O, Peters S O, Adefenwa M A, Yakubu A, Adebambo O A, Ozoje M O, Ikeobi C O, Wheto M, Ajayi O O, Amusan S A, Ekundayo O J, Sanni T M, Okpeku M, Onasanya G O, De Donato M, Ilori B M, Kizilkaya K and Imumorin I G 2012. Morphological and microsatellite DNA diversity of Nigerian indigenous sheep. Journal of Animal Science and Biotechnology 3:38.

Alvarez I, Capote J, Traore A, Fonseca N, Pérez K, Cuervo M, Fernandez I and Goyache F 2012. Genetic relationships of the Cuban hair sheep inferred from microsatellite polymorphism. Small Ruminant Research 104:89– 93.

Arora R J, Bhatia S, Mishra B P, Jain A and Prakash B 2011. Diversity analysis of sheep breeds from Southern peninsular and Eastern regions of India. Tropical Animal Health and Production 43:401–408.

Baumung R, Cubric-Curik V, Schwend K, Achmann R and Solkner J 2006. Genetic characterisation and breed assignment in Austrian sheep breeds using microsatellite marker information. Journal of Animal Breeding and Genetics 123: 265–271.

Cemal İ, Yilmaz O, Karaca O, Binbaş P and Ata N 2013. Analysis of genetic diversity in indigenous Çine Çaparı sheep under conservation by microsatellite markers. Journal of the Faculty of Veterinary Medicine, Kafkas University 19 (3):383-390.

Cerit H, Altınel A, Elmaz Ö and Avanus K 2004. Polymorphism evaluation of various genomic loci in the Kıvırcık sheep breed of Turkey. Turkish Journal of Veterinary and Animal Science 28: 415-425.

Evanno G, Regnaut S and Goudet J 2005. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Molecular Ecology 14: 2611-2620.

Earl-Dent A and vonHoldt Bridgett M 2012. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conservation Genetics Resources vol. 4 (2) pp. 359-361 doi: 10.1007/s12686-011-9548-7

FAO 2011. Molecular genetic characterization of animal genetic resources. FAO Animal Production and Health Guidelines. No: 9, Rome

Goudet J 2001. FSTAT, A Program to estimate and test gene diversities and fixation indices Version 2.9.3. http://www.unil.ch/izea/softwares/fstat.html. (Accessed:15.03.2011)

Grigaliunaite I, Tapio M, Viinalass H, Grislis Z, Kantanen J and Miceikiene I 2003. Microsatellite variation in the Baltic sheep breeds. Veterinarija Ir Zootechnika 21 (43):66-73.

Hoda A and Marsan P A 2012. Genetic characterization of Albanian sheep breeds by microsatellite markers. In: M Caliskan (Ed.), Analysis of Genetic Variation in Animals, InTech, DOI: 10.5772/34554. http://www.intechopen.com/books/analysis-of-genetic-variation-in-animals/genetic-characterization-of-albanian-sheep-breeds-by-microsatellite-markers (Accessed:20.08.2013)

Jyotsana, B, Jakhesara S, Prakash V, Rank D N and Vataliya P H 2010. Genetic features of Patanwadi, Marwari and Dumba sheep breeds (India) inferred by microsatellite markers. Small Ruminant Research 93:57–60.

Koban E 2004. Genetic diversity of native and crossbreed sheep breeds in Anatolia. Department of Biology. PhD thesis, Middle East Technical University, Ankara, p. 125.

Kusza Sz, Ivankovic A, Ramljak J, Nagy I, Javor A and Kukovics S 2011. Genetic structure of Tsigai, Ruda, Pramenka and other local sheep in Southern and Eastern Europe. Small Ruminant Research 99(2-3):130-4.

Lasagna E, Bianchi M, Ceccobelli S, Landi V, Martínez A M, Pla J L V, Panella F, Bermejo J V D and Sarti F M 2011 . Genetic relationships and population structure in three Italian Merino-derived sheep breeds. Small Ruminant Research 96(2-3):111-119.

Ligda C, Altarayrah J, Georgoudis A and the ECONOGENE Consortium 2009. Genetic analysis of Greek sheep breeds using microsatellite markers for setting conservation priorities. Small Ruminant Research 83: 42-48.

Mukesh, M, Sodhi M and Bhatia S 2006. Microsatellite-based diversity analysis and genetic relationships of three Indian sheep breeds. Journal of Animal Breeding and Genetics 123: 258–264.

Nei M 1972. Genetic distance between populations. The American Naturalist 106: 283-292.

Paiva S R, Faria D A, Silverio V C, McManus C, Egito A A, Dergam J A, Guimaraes S E F, Castro S R, Albuquerque M S M and Mariante A S 2005. Genetic variability among Brazilian sheep using microsatellites. The Role of Biotechnology, 5-7 March, Villa Gualino, Turin, Italy, pp.195-196.

Pariset, L, Savarese M C, Cappuccio I and Valentini A 2003. Use of microsatellites for genetic variation and inbreeding analysis in Sarda sheep flocks of central Italy. Journal of Animal Breeding and Genetics 120: 425–432.

Peakall R and Smouse P E 2006. GENALEX 6: Genetic analysis in Excel. Population genetic software for teaching and research. Molecular Ecology Notes 6: 288-295.

Pritchard J K, Stephens M and Donnelly P 2000. Inference of population structure using multilocus genotype data. Genetics 155: 945-959.

Quiroz J, Martinez A M, Zaragoza L, Perezgrovas R, Vega-Pla J L and Delgado J V 2008. Genetic characterization of the autochthonous sheep populations from Chiapas, Mexico. Livestock Science 116:156–161.

SanCristobal M, Chevalet C, Foulley J L and Ollivier L 2003. Some methods for analysing genetic marker data in biodiversity setting example of the Pigbiodiv data. Archivos de Zootecnia 52:173-183.

Schlötterer C 2004. The evolution of molecular markers - just a matter of fashion? Nature Reviews Genetics 5(1):63-9.

Stahlberger-Saitbekova N, Schlaepfer J, Dolf G and C Gaillard 2001 . Genetic relationships in Swiss sheep breeds based on microsatellite analysis. Journal of Animal Breeding and Genetics 118: 379-387.

Takezaki N and Nei M 1996 Genetic distances and reconstruction of phylogenetic trees from microsatellite DNA. Genetics 144: 389-399.

Tamura K, Dudley J, Nei M and Kumar S 2007. MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Molecular Biology and Evolution 24:1596-1599.

Togan İ, Soysal İ, Berkman C C and Koban E 2005 . Irkların korunmasında moleküler işaretler. Tekirdağ Ziraat Fakültesi Dergisi 2(1): 44-9.

Yeh F C, Yang R C, Boyle T B J, Ye Z H and Mao J X 1997 . POPGENE the user-friendly shareware for population genetic analysis. University of Alberta, Canada.

Yilmaz O, Cemal I, Karaca O, Ata N, Sevim S and Ozturk M 2013. Genetic diversity of Karya and Çine Çapari sheep. Scientific Papers Series D Animal Science LVI:31-35.

Yilmaz O and Karaca O 2012. Paternity analysis with microsatellite markers in Karya sheep. Journal of the Faculty of Veterinary Medicine, Kafkas University 18 (5):807-813.

Yilmaz O, Cemal İ and Karaca O 2014. Genetic diversity in nine native Turkish sheep breeds based on microsatellite analysis. Animal Genetics 45 (4): 604–608.

Weir B S and Cockerham C C 1984. Estimating F-statistics for the analysis of population, structure. Evolution 38:1358-1370.


Received 13 November 2015; Accepted 29 December 2015; Published 1 March 2016

Go to top