ÓÄÊ 636.2:575.174:577.212

THE ANALYSIS OF 30 MICROSATELLITE MARKERS IN LOCAL CATTLE POPULATIONS

T.Y. Kiselyova1, B.Y. Podoba2, Y.Y. Zabludovskiy2, V.P. Terletskiy1, N.I. Vorobyev1, J. Kantanen3

An artificial insemination used in commercial cattle breeds has resulted in a loss of genetic diversity. Local cattle breeds may possess unique alleles and their combinations as a result of adaptation to different environmental challenges. Microsatellite markers are widely used for population genetic studies of cattle. 30 microsatellite markers (BM1824, BM2113, BM1818, CSRM60, CSSM66, ETH3, ETH10, ETH152, ETH185, ETH225, HAUT24, HAUT27, HEL1, HEL5, HEL9, HEL13, INRA023, INRA032, INRA035 INRA037, INRA063, ILSTS005, ILSTS006, MM122, INRA005, SPS115, TGLA227, TGLA126, TGLA122, TGLA53) in 238 individuals from 6 local cattle populations (Suksun, Istoben, Yaroslavl, Kholmogory, Grey Ukrainian and Pechora type of Kholmogory) were analyzed. In total, 268 alleles were detected from the 30 loci surveyed. Each locus was characterized according to its length, number of alleles, expected heterozygosity, Fis criterion. Genetic estimation of each population was carried out according to total and mean number of alleles, number of private and effective alleles, observed and expected (unbiased) heterozygosity, Fis criteria and interval confidence (IC). As a result wide polymorphism was revealed both in the loci and the populations studied. All populations had stable genetic structure.

Key words: microsatellite markers, alleles, local breeds, cattle.

 

Intraspecific biodiversity of cattle breeds is the subject of interest in terms of population genetics and preservation of the gene pool, and this is also a source of valuable genotypes for practical breeding. The data about genetic structure of species are important in studying the evolution of farm animals. It should be noted that local populations of cattle retain higher individual variability than commercial breeds. An extensive use of artificial insemination has resulted in loss of several valuable alleles. In experiments with Argentine Creole cattle (a local breed from South America) and American Holstein, it has been shown a higher degree of genetic diversity for loci affecting milk production in local cattle (1). Evaluation of Nordic commercial breeds showed the absence of some alleles of the locus for transferrin (2). In Russia, the ongoing tendency for crowding out of local cattle by commercial breeds can lead to irreversible loss of valuable domestic gene pool formed during a long-term adaptation to local environmental conditions.

Local breeds can carry unique combinations of alleles, which is the important reserve of genetic variability for commercial breeds of domestic animals. These alleles must be detected, conserved and used in selection practice (3-5). In this regard, it is extremely actual to identify genetically unique breeds.

Modern studies of population structure use the most common method for genotyping of alleles of microsatellite loci. Their location is important for searching the candidate genes responsible for commercially valuable quantitative traits (QTLs) and the genes providing manifestation of hereditary diseases (which is used in development of methods for prevention and treatment), because microsatellite alleles can be inherited linked to alleles of such genes. Identification of loci responsible for economically valuable traits is based on an the analysis of length of microsatellite DNA fragments using appropriate statistical programs containing the data about dairy and meat productivity (6, 7). The large number of alleles provides accurate identification of any individual since microsatellites are found almost in all living organisms. Currently, evaluation of animal’s origin upon microsatellite markers is more precise than upon blood groups. That’s why the international practice of sale abroad of breeding animals obligates to confirm their origin by the analysis of microsatellite DNA. There are recommendations by FAO about types of microsatellites used for genetic evaluation of cattle (8).

According to INRA (French National Institute for Agricultural research) database, there are 2402 microsatellites identified in cattle today, including 2,244 mapped ones (microsatellites are present on all 30 pairs of chromosomes). Genetic analysis of  populations involves using packages of computer programs - Arlequin, GDA, Gen-AlEx, GENEPOP, GeneClass, Genetic Studio, MicroSat, PowerMarker, which are widely used in agricultural genetics.

The purpose of this work was to study the polymorphism of 30 microsatellite markers in six local populations of cattle.

Technique. The analysis was performed upon DNA extracted from blood samples of 238 animals from populations the following local breeds: Kholmogory (pedigree farm “Arkhangel’skii”, Kholmogorsky district, Arkhangelsk province), Pechora type of Kholmogory cattle (pedigree farm “Ukhta-97”, Ukhta, Komi Republic), Gray Ukrainian (subsidiary farm of Kiev-Pechersk Lavra, Boryspil district, Ukraine), Istoben (pedigree farm “Istobenskii”, Orichevsky district, Kirov province), Suksun (the farm “Savinskoe”, Suksun district, Perm province) and Yaroslavl (the farm “Savinskoe” and the experimental farm “Tutaevo”, Tutaev district, Yaroslavl province).

30 microsatellite markers were used: BM1824, BM2113, BM1818, CSRM60, CSSM66, ETH3, ETH10, ETH152, ETH185, ETH225, HAUT24, HAUT27, HEL1, HEL5, HEL9, HEL13, INRA023, INRA032, INRA035, INRA037, INRA063, ILSTS005, ILSTS006, MM122, INRA005, SPS115, TGLA227, TGLA126, TGLA122, TGLA53.

DNA was isolated by standard phenol method (9). The resulting DNA and primers were diluted with deionized water to concentrations 10 ng/ul and 10 pmol/ul. Conditions of polymerase chain reaction (PCR) were set according to the database of  Roslin Institute (UK). Microsatellites with sharply differing lengths of alleles were analyzed in pairs. PCR products were separated in 6% polyacrylamide gel. Lengths of microsatellites were measured on the devices ALF Express automated DNK sequencer (“Pharmacia”, Sweden) and MegaBACE™ 500 DNA sequencer (“Amersham Biosciences”, USA). To identify alleles of microsatellite loci, the program ALF Win Fragment Analyser 1.00.36 was used. True length of allele was determined by comparing it with lengths of alleles of marker animals. On the obtained charts, locations of peaks were assessed to determine the length of a microsatellite (bp).

Each locus was evaluated by its length, number of alleles, observed and expected heterozygosity, Fis factor. The polymorphism within a population was assessed by total number of identified alleles, observed and expected heterozygosity (unbiased estimate), the average number of alleles per locus, the number of private and effective alleles, Fis coefficient with allowance for the confidence interval.

Statistical analysis of data was performed using the programs GE-NEPOP, Genetix, Convert, and own approved computer programs developed on the basis of common mathematical models (10, 11).

Results. The analysis of 30 microsatellites located on  the 21st pair of chromosomes (1st - 3rd, 5th, 7th -12th, 14th, 15th – 23rd and 26th) revealed polymorphism in all loci (Table 1).

                                                                                                  

1. Characteristics of 30 microsatellite loci identified in six local breed populations of cattle

Microsatellite

Chromosome

Length of allele, bp

Maximum number of alleles

Expected
heterozigosity

Fis

BM1824

1st

179-191

7

0,7379

0,028

BM2113

2nd

122-142

10

0,8258

0,065

ETH10

5th

211-223

7

0,6915

0,080

ETH225

9th

139-153

8

0,7749

0,085

ETH3

19th

115-129

8

0,7786

0,196

HEL5

21st

151-167

7

0,7893

0,110

ILSTS005

10th

184-194

3

0,3871

0,075

INRA023

3rd

197-219

10

0,8020

0,056

INRA005

12th

139-145

4

0,5769

0,020

INRA063

18th

175-183

5

0,6135

0,074

BM1818

23rd

256-268

7

0,6664

0,053

CSSM66

14th

179-199

12

0,8578

0,075

ETH152

5th

193-207

7

0,7184

0,004

HEL1

15th

101-113

7

0,7268

0,061

HEL13

11th

178-194

6

0,6125

0,146

HEL9

8th

143-171

13

0,7816

0,087

ILSTS006

7th

285-301

10

0,7932

0,067

INRA032

11th

172-186

7

0,6598

0,076

INRA037

10 th

120-147

13

0,6122

0,109

TGLA227

18 th

79-105

13

0,8840

0,067

TGLA126

20 th

115-129

8

0,7393

0,032

TGLA122

21st

138-182

17

0,7785

0,039

HAUT24

22nd

106-128

10

0,7981

0,066

INRA035

16 th

102-122

6

0,4645

0,509

HAUT27

26 th

128-156

10

0,6915

0,142

CSRM60

10 th

93-107

8

0,7824

0,016

MM12

9 th

113-133

11

0,7633

-0,014

ETH185

17 th

220-244

11

0,8517

0,047

TGLA53

16 th

151-183

17

0,8755

0,085

SPS115

15 th

244-257

6

0,6902

0,094

Note. Fis — inbreeding coefficient reflecting deficit or excess of heterozygotes in population

Different microsatellites had unequal maximum number of alleles – from 3 (ILSTS005) to 17 (TGLA122, TGLA53). The microsatellite ILSTS005 was designed for zebu-like cattle which has five alleles (12). Totally, 268 alleles in 30 microsatellite loci were found. Expected heterozigosity was quite high for all microsatellites (0,6122-0,8840) except INRA035 (0,4645) and ILSTS005 (0,3871). The microsatellite INRA035 demonstrated fairly close inbreeding: Fis coefficient (shows the decrease of observed heterozigosity compared with expected one) for INRA035 equaled to 0,509 as the sum for all populations. Fis was negative for only microsatellite (MM12; Fis = -0,014) indicating a weak outbreeding; for other microsatellites, it was weakly positive.

The analysis of Hardy-Weinberg genetic equilibrium using the exact method and the Bonferroni correction revealed significant deviations of genotype frequencies of the microsatellite INRA035 in four populations: Suksun (P < 10-4), Yaroslavl (P < 10-4), Kholmogory (P = 3×10-4) and in the Pechora type of Kholmogory (P = 10-3), since the critical level of significance of individual tests for 30 loci shall not exceed 0,0017 to be significant at 0,05%. In Istoben and Gray Ukrainian breed populations, P-values were, respectively 0,0051 and 0,0541.

2. Genetic characteristics of six local breed populations of cattle upon the polymorphism of 30 microsatellites

Indi-
cator

Suksun

Istoben

Yaroslavl

Kholmogory

Pechora type of Kholmogory

Gray Ukrainian

1

182

184

194

176

164

152

2

6,07±0,36

6,13±0,34

6,47±0,45

5,87±0,37

5,43±0,32

5,07±0,30

3

7

10

13

6

1

12

4

3,33

3,32

3,31

2,90

2,93

2,86

5

0,700±
0,019

0,685±
0,021

0,698±
0,022

0,656±
0,025

0,659±
0,025

0,651±
0,022

6

0,681±
0,028

0,690±
0,031

0,698±
0,030

0,634±
0,030

0,644±
0,026

0,643±
0,024

7

+0,0125

-0,0104

-0,0163

+0,0180

+0,0088

+0,0040

8

 (-0,02267… +0,04858)

 (-0,04582… +0,01266)

 (-0,04036… +0,01808)

 (-0,02271… +0,04840)

 (-0,03750… +0,04620)

 (-0,05115… +0,03876)

Note. 1 – total number of alleles, 2 – mean number of alleles per locus, 3 – number of unique alleles, 4 – number of effective alleles, 5 – expected heterozigosity, 6 – observed heterozigosity, 7 - Fis, 8 – IC (confidence interval, 95 %).

All populations differed in number of identified alleles - from 152 to 194 (Table 2). The mean number of alleles per locus for 30 microsatellites amounted to 8,93. According to the literature, this indicator in German Simmental equals 6,64, in Swedish Holstein  - 6,52 for 25 microsatellite markers included in the abovementioned list (5). In Black-and-White cattle, this parameter was 9,00 in the analysis of 7 microsatellites with maximum number of alleles while the higher diversity in length of microsatellites: TGLA227 - 75-105, ETH225 - 132-160, ETH10 - 207-231, BM1818 - 248-278, BM2113 - 151-155, BM1824 - 176-196 and SPS115 - 234-258 bp (13). In this study, in cases when different breed populations had equal number of microsatellite alleles, lengths of these alleles were distinct. Thus, in Yaroslavl and Gray Ukrainian cattle, five alleles of the microsatellite BM1824 were found, whose lengths were, respectively 179, 181, 183, 189, 191 and 179, 181, 183, 187, 189 bp. In this case, the alleles of 187 and 191 bp length were specific for these breeds: the first one was detected in Gray Ukrainian and wasn’t found in Yaroslavl cattle, the second allele was present in Yaroslavl breed and absent - in Gray Ukrainian.

Unique alleles discovered in only one population are of the great interest. In Suksun population, 7 such alleles were identified, in Istoben - 5, in Yaroslavl - 13, in Kholmogory - 6, in Gray Ukrainian - 12 and in the Pechora type of Kholmogory cattle - 1. Gray Ukrainian and Yaroslavl breeds had almost similar  number of unique alleles. It is possible that during the evolution of such breed, its long-time farming in small populations with forced inbreeding resulted in development of inner mechanisms preventing the restriction of genetic diversity. This fact was indicated by the analysis of inheritance of alleles for B-system blood groups in the genetic fund of Gray Ukrainian cattle the herd of pedigree farm “Polivanovka” (14). In this population, actual number of heterozygotes was 18% while the expected value of 28% , and it has been found the intra-gene recombinations resulted in emergence of new alleles (15, 16).

The number of effective alleles per locus (provides main contribution into total homo- and heterozygosity) in breed populations of cattle varied insignificantly: from 2,86 (Gray Ukrainian) to 3,33 (Suksun) (Table 2). No differences in observed and expected heterozygosity as well as in average number of alleles per locus were detected. Coefficient of inbreeding (Fis) showed the negligible prevalence of heterozygous  alleles (the sum over all loci) in Istoben (-0,0104) and in Yaroslavl (- 0,0163) populations, while the other breeds demonstrated low positive values indicating a slight predominance of homozygous alleles (Suksun - +0,0125, Kholmogory – +0,0180, Pechora type of Kholmogory - +0,0088 and Gray Ukrainian - +0,0040). In all groups of animals, Fis was within the confidence interval (95%), which confirmed the absence of both inbreeding and outbreeding, therefore, the studied populations had stable genetic structure.

Thus, the six studied local breed populations of cattle were analyzed upon 30 microsatellite markers, which has revealed a wide  polymorphism of alleles’ lengths. In all of these populations unique microsatellite alleles have been found. The established levels of  observed and expected heterozygosity, as well as Fis-coefficient indicated a stable genetic structure (no inbreeding or outbreeding were reliably detected). Polymorphic microsatellite markers are widely used in breeding work. In the world practice of cattle farming, there are available commercial kits (DNA tests) for identification the individuals with best tenderness and marbling of meat (http://www.bovigen.com, http://www.idenity.com, http://www.nbcec.org).

REFERENCES

1. Golidjow C.D., Giovambattista G., Ripoli M.V., Dulout F.N. and Lojo M.M., Genetic Variability and Population Structure in Loci Related to Milk Production Traits in Native Argentine Creole and Commercial Argentine Holstain Cattle, Genet. Mol. Biol., 1999, vol. 22, no. 3, pp. 395-398.
2. Kantanen J., Olsaker I., Holm L.-E. et al., Genetic Diversity and Population Structure of 20 North European Cattle Breeds, J. Heredity, 2000, vol. 91, pp. 446-457.
3. Hall S.J.G. and Bradley D.G., Conserving Livestock Breed Biodiversity, Trends Ecol. & Evol., 1995, vol. 10, no. 7, pp. 267-270.
4. Notter D.R., The Importance of Genetic Diversity in Livestock Population of the Future, J. Animal Sci., 1999, vol. 77, pp. 61-69.
5. Grzybowski G. and Prusak B., Genetic Variation in Nine European Cattle Breeds as Determined on the Basis of Microsatellite Markers. III. Genetic Integrity of the Polish Red Cattle Included in the Breeds Preservation Program, Animal Sci. Papers and Reports, 2004, vol. 22, no. 1, pp.  45-46.
6. Velmala R.J., Vilkki H.J. and Elo K.T., A Search for Quantitative Trait Loci for Milk Production Traits on Chromosome 6 in Finnish Ayrshire Cattle, Animal Genet., 1999, vol. 30, no. 2, pp. 136-143.
7. Gutierrez-Gil B., Wiener P., Nute G.R. et al., Detection of Quantitative Trait Loci for Meat Quality Traits in Cattle, Animal Genet., 2008, vol. 39, pp. 51-61.
8. Bradley D.G., Genetic Diversity in Cattle,  Proc. ISAG meeting, Tours, France, 1996 (http://www.projects.roslin.ac.uk/cdiv/markers.html).
9. Sambrook J., Fritsch E.F. and Maniatis T., Molecular Cloning: a Laboratory Manual, Cold Spring Harbor: Cold Spring Harbor Laboratory Press, N.Y., 1989.
10. Khedrik F., Genetika populyatsij (Population Genetics), Moscow, 2003.
11. Weir B.S., Genetic Data Analysis II, Sunderland, USA, 1996.
12. Pandey A.K., Sharma R., Singh Y. et al., Genetic Diversity Studies of Kherigarh Cattle Based on Microsatellite Markers, J. Genet., 2006, vol. 85, no. 2, pp. 117-122.
13. Zhivotovskii L.A., Populatsionnaya biometriya (Population Biometry), Moscow, 1991.
14. Milofeeva L.A., Development of the System for Multiloci Gene Analysis of Microsatellites and its Use at Description of the Allele Fund in the Black-and-White Cattle Breed, Extended Abstract of Cand. Sci. Dissertation, Dubrovitsy, 2009.
15. Eisner F.F., Vlasov V.I. and Podoba B.E., Genetic Problems in Cattle Breeding, Vestn. s.-kh. nauki, 1985, no. 3, pp. 73-81.
16. Podoba B.E. and Misostova N.V., Recombinations in the B Loci for Blood Groups of Cattle, Nauch.-tekhn. bull. NII zhivotnovodstva Lesostepi i Poles’ya USSR, Kyiv, 1974, no. 11, pp. 43-46.

1All-Russia Research and Development Institute of Farm Animal Genetics and Breeding, RAAS, 
St.-Petersburg – Pushkin 196600, Russia
,
e-mail: ktyu@yahoo.com, valeriterl@mail.ru,
nikolai_vorobyov@yahoo.com;
2Institute of Animal Breeding and Genetics, UAAS, Kyiv province, Borispil region, Chubinsky settlement 08321, Ukraine,  
e-mail: mba@kbs.kiev.ua;
3Agricultural Research Centre of Finland
Animal Production Research,

FIN-31600 Jokioinen, Finland,
e-mail: juha.kantanen@mtt.fi

Ïîñòóïèëà â ðåäàêöèþ
6 àïðåëÿ 2010 ãîäà