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

Citation of this paper

Genetic diversity of Senepol cattle in Colombia using ten multiplexed microsatellites

Natalia Londoño Díaz, Mario Fernando Cerón Muñoz1 and Iván Darío Soto Calderón2

Laboratorio de Genética Animal, Grupo de Investigación GaMMA, Universidad de Antioquia, Calle 70 Nº 52 - 21, Medellín, Antioquia, Colombia
ivan.soto@udea.edu.co
1 Facultad de Ciencias Agrarias, Universidad de Antioquia, Calle 70 Nº 52 - 21, Medellín, Antioquia, Colombia
2 Instituto de Biología, Universidad de Antioquia, Calle 70 Nº 52 - 21, Medellín, Antioquia, Colombia

Abstract

An automated multiplex PCR assay of 10 polymorphic microsatellite loci was optimized here in Senepol bovine cattle. Allele frequencies were estimated for Colombian herds. The power of this set of markers for parentage and kinship analyses was also tested. Population genetic parameters, combined probability of identity (CPI) and combined probabilities of exclusion (CPE) were estimated.

The average number of alleles per locus was 14, the effective number of alleles was 5, the observed and expected heterozygosities were 0.635 and 0.776, respectively, the Shannon Index was 2.70 and the polymorphic information content was 0.751. These values evidenced elevated levels of polymorphism and an underestimation of genetic diversity through the use of manual genotyping methods. Inbreeding coefficient (Fis=0.206) in Colombian Senepol herds was moderate. Altogether, data from the used set of microsatellite loci yielded a CPI between two random samples of 1.07E-12. The CPE in parentage testing were >0.998 for the first parent and >0.999 for both the second parent and parent pairs. This set of polymorphic markers can be genotyped in a relatively easy, effective and reliable way for relatedness and forensic analyses. Finally, we recommend our approach for population genetic analyses of other Senepol populations and future genetic monitoring of inbreeding in Senepol herds.

Key words: Bos taurus, genetic relatedness, genetic variation, molecular marker, parentage



Diversidad genética de ganado Senepol en Colombia empleando diez microsatélites combinados

Resumen

Se optimizó en ganado bovino Senepol una prueba automatizada de PCR multiplex de 10 loci microsatélite polimórficos. Adicionalmente, se estimaron las frecuencias alélicas en hatos colombianos y se evaluó el poder de este grupo de marcadores en pruebas de paternidad y parentesco. Se estimaron parámetros genético-poblacionales, probabilidad de identidad combinada (CPI) y probabilidades de exclusión combinadas (CPE).

El número promedio de alelos por locus fue 14, el número efectivo de alelos fue 5, las heterocigocidades observada y esperada fueron 0.635 y 0.776, respectivamente, el índice de Shannon fue 2. 70 y el contenido de información polimórfica fue 0.751. Estos valores evidenciaron altos niveles de polimorfismo y una subestimación de la diversidad genética cuando se usan métodos manuales de genotipificación. El coeficiente de endogamia en hatos colombianos de Senepol fue moderado (Fis=0.206). En general, los datos del grupo utilizado de microsatélites arrojaron un CPI entre dos muestras aleatorias de 1.07E-12. El CPE in las pruebas de parentesco fue >0.998 para el primer parental y >0.999 tanto para el segundo parental como para el par de parentales. Este grupo de marcadores microsatélite es relativamente fácil de genotipificar en una forma efectiva y confiable para análisis forenses y de parentesco. Finalmente, recomendamos nuestra aproximación para análisis genéticos-poblacionales de otras poblaciones de Senepol y futuros monitoreos genéticos de endogamia en hatos de Senepol.

Palabras clave: Bos taurus, marcador molecular, parentesco genético, paternidad, variación genética


Introduction

Senepol is a Bos taurus cattle breed endemic to the New World (Hupp 1978), that emerged in 1810’s in St. Croix, Virgin Islands, after crossing N’Dama and Red Poll. It was introduced in Colombia from the United States, Virgin Islands, Venezuela and Brazil in the early 2000’s and it is mainly distributed in the North West in the departments of Antioquia, Córdoba and Cundinamarca. Senepol was subsequently introduced in other departments and the current populations size is estimated at around 2700 individuals (Asosenepol 2015). Two previous population genetic studies of Colombian Senepol have made use of ten polymorphic microsatellite loci recommended by the International Society for Animal Genetics (ISAG) (Montoya et al 2010; Sepúlveda et al 2012). These studies revealed remarkable levels of genetic variability and evidenced the utility of such markers for forensic and relatedness analyses. Such analyses have been based on single amplification of each marker followed by manual genotyping and allele scoring through analysis of polyacrilamede gel electrophoresis (PAGE) profiles. However, this approach is expensive as well as time and effort consuming, but also allele scoring is prone to errors, thus limiting practicity, replicability and reliability.

Therefore, the aim of this study was to optimize a routine methodology of multiplex amplification and automated genotyping of the same ten polymorphic microsatellite loci in a sample of the Senepol Colombian population and assess the potential of these markers for parentage and individual identification analyses.


Materials and methods

Sampling and genotyping

We selected 302 Senepol individuals, representing over 10% of the current population. Genomic DNA of 302 Senepol individuals from 25 farms in seven Colombian departments (Cundinamarca, Antioquia, Boyacá, Casanare, Córdoba, Meta and Tolima) had been previously extracted from blood samples kindly provided by Asociación Colombiana de Criadores de Ganado Senepol y sus Cruces (ASOSENEPOL) (Sepúlveda et al 2012).

A set of ten polymorphic microsatellites previously used in parentage tests and recommended by The International Society of Animal Genetics (ISAG) was selected for this study. PCR reactions were carried out in a total volume of 15µl containing 1X PCR buffer, 0.2 mM dNTPs, between 1.0 and 3.0 mM MgCl 2, 0.2µM of each primer, 0.5U de Taq polymerase (Thermo Fisher Scientific, Waltham, MA, USA) and 30 to 50µg of template DNA (Table 1). Initial denaturation at 94°C for 5 min was followed by 35 cycles of denaturation at 94°C for 30 s, primer annealment for 30 s (between 50 and 65°C) and extension at 72°C for 1 min with a final extension step at 72°C for 10 min.

Primer sets with similar thermal profiles and MgCl2 requirements were pooled in the same tube and products were coamplified. To do this, we first adjusted the amplicon size of four loci by redesigning one or both primers from the reference Bos taurus genome (Genbank Accession number GCA_000003055.3) with the software Primer3 v4.0 (Untergrasser et al 2012). The 5’ end of all forward primers was labeled with one of the commercial dyes FAMTM, VIC®, PET® or NEDTM (Applied Biosystems, Foster City, CA, USA), making sure that coamplified markers had different labels (Table 1).

Amplicons were run in an ABI genetic Analyzer 3130 (Applied Biosystems) using LIZ500® as molecular weight standard (size range between 35-500 bp). Genotype editing and allele calling was performed making use of the software Geneious v6.1.6 (Biomatters Ltd., Auckland, New Zealand).

Data analysis

Deviation from the Hardy-Weinberg equilibrium (HWE) was tested across all loci with a modified version of the Markov chain random walk algorithm (Guo and Thompson, 1992). HWE, observed and expected heterozygosities (Ho and He, respectively) and Inbreeding coefficient (Fis) were estimated in the package Arlequin v3.11 (Excoffier and Lischer 2010).

Table 1. Primer sequences, chromosome location, allele size and multiplexed amplification conditions of ten microsatellite loci in Senepol cattle.

Locus

Primer sequence 5’-3’
(Fluorochrome)
**

Chm
***

Size range
(bp)

[MgCl2]
mM

Temp.
(°C)

Multiplex

BM2113

F:GCTGCCTTCTACCAAATACCC
(FAM)
R:CTTCCTGAGAGAAGCAACACC

2

118-148

2.0

65.0

M1

SPS115

F:AAAGTGACACAACAGCTTCTCCAG
(PET)
R:AACGAGTGTCCTAGTTTGGCTGTG

15

240-260

TGLA126

F1:CCCAGATTTCTTGACTCATCG *
(VIC)
R1:AGTCCCAGCCTCTGATCTCT *

20

240-318

ETH225

F1:GGGGATTATCTGTGGCATAA *
(FAM)
R:ACATGACAGCCAGCTGCTACT

9

184-212

ETH10

F:GTTCAGGACTGGCCCTGCTAACA
(NED)
R:CCTCCAGCCCACTTTCTCTTCTC

5

200-218

2.0

65.0

M2

TGLA122

F:CCCTCCTCCAGGTAAATCAGC
(VIC)
R:AATCACATGGCAAATAAGTACATAC

21

132-178

INRA32

F:AAACTGTATTCTCTAATAGCTAC
(NED)
R:GCAAGACATATCTCCATTCCTTT

11

158-196

2.5

62.5

M3

INRA64

F1:GCATCAACACAGGGTTGTTTA *
(PET)
R:CTGAAAGCAGAATGAGGTGC

23

218-228

BM1824

F:GAGCAAGGTGTTTTTCCAATC
(PET)
R:CATTCTCCAACTGCTTCCTTG

1

136-192

2.5

59.5

M4

INRA37

F1:GGCGAAGGTTTGGATACAT *
(FAM)
R:AAAATTCCATGGAGAGAGAAAC

10

254-288

* Modified from original primer sequence. ** Fluorochrome colors: Blue (FAM), Green (VIC), red (PET) and yellow (NED). *** Chromosome.

Polymorphism level and genetic diversity were assessed through the estimation of allele frequencies, total number of alleles per locus (Na), effective number of alleles (Ne) and log2-base Shannon Index (SI) in the software GenAlEx v6.5 (Peakall and Smouse 2006). We also calculated Polymorphic Information Content (PIC), Probability of Identity (PI) and Probability of Exclusion (PE) with the software Cervus v3.0.6 (Kalinowsky et al 2007). PE is defined as the probability of excluding a candidate parent or parent pair from parentage of a given offspring in three possible scenarios: testing one parent when the other parental genotype is unavailable (PE1 or first parent); testing one parent when the genotype of the other parent is available (PE2 or second parent) and; when a parent pair is tested (PE3 or parent pair) (Jamieson and Taylor, 1997; Marshall et al 1998).


Results

The ten microsatellite loci were successfully amplified and scored in four multiplex reactions consisting on one set of four (M1) and three sets of two (M2, M3 and M4) pooled markers (Table 1). The average Na was 14.3, which ranged between six (INRA64) and 28 (TGLA126) alleles per locus. Ne ranged from 1.9 to 12.1 with an average of 5.6. Average values for SI and PIC were 2.70 and 0.752, respectively (Tables 2 and 3). All ten tested loci exhibited significant deviation from the Hard-Weinberg expectations (p<0.001) as determined by observed heterozygosity deficit ( Ho vs He) and inbreeding levels that reached 21% (Fis = 0.206) and. The average PI per locus was 0.087 with a range between 0.013 and 0.302. The Combined PI (CPI) of the ten-locus set was 1.07E-12. As far as the probability of exclusion, PE1 ranged between 0.126 and 0.716, PE2 between 0.279 and 0.834 and, PE3 between 0.447 and 0.955. The combined analysis of all ten microsatellite loci yielded a CPE1 >0.998, whereas both CPE2 and CPE3 were >0.999.

Table 2. Genetic diversity indices and power for parentage testing of ten microsatellite loci in Colombian Senepol cattle.

TGLA
122

TGLA
126

INRA
32

ETH
10

BM
2113

ETH
225

INRA
37

BM
1824

INRA
64

SPS
115

Average

N

295

298

297

297

299

298

299

297

298

298

N a

16

28

18

8

14

12

16

15

6

10

14.3

N e

4.69

12.1

7.51

5.46

5.15

3.78

7.03

3.52

1.93

4.43

5.56

Ho

0.732

0.762

0.670

0.801

0.696

0.517

0.696

0.609

0.205

0.661

0.635

He

0.788

0.919

0.868

0.818

0.808

0.737

0.859

0.717

0.484

0.776

0.777

SI

2.83

3.93

3.27

2.55

2.77

2.48

3.06

2.24

1.42

2.44

2.70

PIC

0.762

0.912

0.853

0.790

0.784

0.710

0.842

0.673

0.448

0.741

0.752

PI

0.070

0.013

0.031

0.060

0.060

0.095

0.036

0.123

0.302

0.084

0.087

CPI

Combined Probability of Identity

1.07E-12

PE1

0.432

0.716

0.580

0.455

0.457

0.359

0.502

0.313

0.126

0.309

0.425

CPE1

Combined Probability of Exclusion (first parent)

>0.998

PE2

0.609

0.834

0.735

0.631

0.633

0.545

0.714

0.488

0.279

0.568

0.604

CPE2

Combined Probability of Exclusion (second parent)

>0.999

PE3

0.803

0.955

0.896

0.808

0.820

0.750

0.878

0.679

0.447

0.754

0.779

CPE3

Combined Probability of Exclusion (parent pair)

>0.999

N: Sample size; Na: number of alleles; Ne: effective number of alleles; Ho: observed heterozygosity; He: expected heterozygosity; SI: Shannon Index; PIC: polymorphic information content; PI: probability of identity; PE1: probability of exclusion (first parent); PE2: probability of exclusion (second parent); PE3: probability of exclusion (parent pair).



Table 3. Allele frequencies and allele sizes (in parenthesis) of ten microsatellite loci in Colombian Senepol cattle.

Allele

TGLA
122

TGLA
126

INRA
32

ETH
10

BM
2113

ETH
225

INRA
37

BM
1824

INRA
64

SPS
115

N

295

298

297

297

299

298

299

297

298

298

1

0.008
(132)

0.003
(240)

0.020
(158)

0.003
(200)

0.007
(118)

0.003
(184)

0.003
(254)

0.007
(136)

0.008
(218)

0.003
(240)

2

0.017
(134)

0.003
(242)

0.010
(160)

0.224
(206)

0.012
(120)

0.008
(186)

0.002
(256)

0.002
(140)

0.035
(220)

0.008
(244)

3

0.015
(136)

0.032
(244)

0.005
(162)

0.202
(208)

0.035
(122)

0.029
(188)

0.007
(258)

0.002
(142)

0.176
(222)

0.304
(246)

4

0.017
(138)

0.020
(246)

0.005
(164)

0.039
(210)

0.192
(124)

0.015
(190)

0.012
(260)

0.012
(144)

0.693
(224)

0.138
(248)

5

0.053
(140)

0.035
(248)

0.007
(170)

0.013
(212)

0.050
(126)

0.154
(192)

0.025
(262)

0.003
(152)

0.057
(226)

0.300
(250)

6

0.253
(142)

0.042
(250)

0.066
(172)

0.150
(214)

0.129
(128)

0.060
(194)

0.007
(264)

0.002
(158)

0.030
(228)

0.144
(252)

7

0.025
(144)

0.097
(252)

0.133
(174)

0.177
(216)

0.012
(130)

0.458
(196)

0.055
(266)

0.005
(172)


0.029
(254)

8

0.005
(146)

0.035
(254)

0.098
(176)

0.192
(218)

0.104
(132)

0.092
(198)

0.196
(268)

0.007
(178)


0.022
(256)

9

0.015
(148)

0.029
(256)

0.19
(178)


0.100
(134)

0.015
(206)

0.092
(270)

0.133
(180)


0.045
(258)

10

0.361
(150)

0.022
(258)

0.221
(180)


0.339
(136)

0.128
(208)

0.169
(272)

0.288
(182)


0.007
(260)

11

0.078
(152)

0.027
(260)

0.099
(182)


0.007
(138)

0.034
(210)

0.134
(274)

0.421
(184)


12

0.068
(154)

0.002
(262)

0.069
(184)


0.01
(140)

0.003
(212)

0.002
(278)

0.057
(186)



13

0.007
(156)

0.003
(264)

0.019
(186)


0.002
(142)


0.003
(282)

0.005
(188)



14

0.017
(158)

0.002
(266)

0.005
(188)


0.002
(148)


0.005
(284)

0.056
(190)



15

0.059
(160)

0.003
(268)

0.040
(190)




0.100
(286)

0.002
(192)



16

0.002
(178)

0.002
(276)

0.003
(192)




0.189
(288)




17


0.002
(280)

0.003
(194)








18


0.002
(284)

0.007
(196)








19


0.002
(286)









20


0.002
(290)









21


0.002
(304)









22


0.077
(306)









23


0.084
(308)









24


0.136
(310)









25


0.087
(312)









26


0.143
(314)









27


0.067
(316)









28


0.035
(318)










Discussion

The microsatellite data showed a large number of alleles and high degree of polymorphism that nonetheless contrast with moderate levels of inbreeding and heterozygosity deficit that lead to departure from the Hardy-Weinberg equilibrium. In particular, this is expected to occur as a consequence of finite size of individual herds, strong trait selection and limited genetic exchange among herds. In case this trend remains unchanged over time, inbreeding levels could substantially increase and facilitate loss of genetic diversity.

Even though only a few previous studies have assessed levels of genetic variation in Senepol, these seem nonetheless to be relatively high in Colombian herds as compared with other populations. Brenneman et al (2007) estimated an average Na of 6.6 from 26 microsatellite loci in a herd in Florida (USA) that was initially established from Caribbean founders. Also, Silva Filho et al (2013) assessed the genetic diversity of a herd in Minas Gerais (Brazil) making use of ten polymorphic microsatellites that yielded an average Na of 5.3 and Ne of 3.1. Unlike the present study in Colombia, which showed a Na of 14.3 and Ne of 5.6, both analyses were focused on single herds and evidenced genetic diversity levels below the values in Colombian herds. Whether the smaller diversity diversity observed in previous studies is determined by limited sample sizes or historical factors of those populations remains to be determined in future studies.

Overall, combined multilocus genotypes provide high power for both individual identification and parentage testing in the Colombian Senepol population (Radko 2010; Souza et al 2012; Stevanović et al 2009). Unlike our two previous genetic surveys of the Colombian Senepol population, which were based on PAGE genotyping (Montoya et al 2010; Sepúlveda et al 2012), the present study implemented an automated approach of multiplex genotyping and allele scoring. Interestingly, the number of alleles per locus increased nearly 100% in the present study, which unveils an underestimation of genetic diversity using manual methods. The mismatch between these two genotyping assays is reasonable given the limited resolution of PAGE-based methods. This limitation is especially evident in the detection of allelic variants in microsatellites with dinucleotide repeat units, which are a common type of repeat in many species and the case of all ten loci implemented in the present study (Schug et al 1998; Li et al 2002).

In summary, the protocol presented here is cost and time effective since it provides a routine alternative to amplify ten microsatellite loci in four multiplex PCR reactions. This marker set revealed substantial levels of polymorphism and evidence of inbreeding in Colombian herds of Senepol cattle. These markers yielded low probability of identity, which is suitable for forensic analysis as well as individual recognition, but they also showed high probability of exclusion under three different scenarios in paternity tests (first parent, second parent and parent pair), which will be useful for future population genetic analysis of Senepol in Colombia and elsewhere.


Acknowledgements

To the University of Antioquia (CODI-Sostenibilidad UdeA 2016 and CODI-Menor 2009) and the Colombian Association of Senepol Cattle Breeders (ASOSENEPOL) for financial support and supply of biological samples.


References

Asosenepol 2015 Origen de la raza Senepol y su ingreso a Colombia. Departamento Técnico, Santafé de Bogotá, Revisado en: http://asosenepolcolombia.com/portal/images/stories/docs/Historia_Senepol_y_en_Colombia.pdf

Brenneman R A, Chase C C, Olson T A, Riley D G and Coleman S W 2007 Genetic diversity among Angus, American Brahman, Senepol and Romosinuano cattle breeds. Animal Genetics 38(1):50–53.

Excoffier L and Lischer H E 2010 Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources 10(3): 564–567.

Guo S W and Thompson E A 1992 Performing the exact test of Hardy-Weinberg proportion for multiple alleles. Biometrics 48(2):361–372.

Hupp H D 1978 History and development of Senepol cattle (No. 11). Agricultural Experiment Station, College of the Virgin Islands. 12 p.

Jamieson A and Taylor C S 1997 Comparisons of three probability formulae for parentage exclusion. Animal Genetics 28(6):397–400.

Kalinowsky S T, Taper M L and Marshall T C 2007 Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Molecular Ecology 16(5):1099–1106.

Li Y C, Korol A B, Fahima T, Beiles A and Nevo E 2002 Microsatellites: genomic distribution, putative functions and mutational mechanisms: a review. Molecular Ecology 11(12):2453–2465.

Marshall T C, Slate J B K E, Kruuk L E B and Pemberton J M 1998 Statistical confidence for likelihood-based paternity inference in natural populations. Molecular Ecology 7(5):639–655.

Montoya A E, Cerón M F, Moreno M A, Martínez E, Corrales J D, Tirado J F and Calvo S J 2010 Genetic characterization of the Hartón del Valle, Angus, Brangus, Holstein, and Senepol cattle breeds in Colombia, using ten microsatellite markers. Revista Colombiana de Ciencias Pecuarias 23(3):283–291. http://www.scielo.org.co/scielo.php?pid=S0120-06902010000300003&script=sci_arttext&tlng=es

Peakall R O D and Smouse P 2006 GenAlEx 6.5: Genetic analysis in Excel. Population software for teaching and research. Molecular Ecology Notes 6(1):288–295.

Radko A 2010 Application of a complementary set of 10 microsatellite DNA markers for parentage verification in Polish Red cattle. Annals of Animal Science 10(1):9–15. http://www.izoo.krakow.pl/czasopisma/annals/2010/1/art2.pdf

Schug M D, Wetterstrand K A, Gaudette M S, Lim R H, Hutter C M and Aquadro C F 1998 The distribution and frequency of microsatellite loci in Drosophila melanogaster. Molecular Ecology 7(1):57–70.

Sepúlveda J C, Ángel-Marín P A, Toro A, Corrales J D, Moreno M A and Cerón-Muñoz M F 2012 Genetic variability of Senepol cattle in Colombia using molecular markers. Revista Colombiana de Ciencias Pecuarias 25(2):183-190.   http://rccp.udea.edu.co/index.php/ojs/article/view/774/837

Silva Filho E, Dos Santos Y M, Margarido Y M M, Mendes T J, De Lima Pinheiro L M, De Souza P H and Campelo J E G 2013 Genetic variability of Guzerat and Senepol bovine breeds by microsatellite markers. Journal of Biotechnology and Biodiversity 3(3):276–282. http://www.revista.uft.edu.br/index.php/JBB/article/viewArticle/592

Souza C A, Paiva SR, McManus C M, Azevedo H C, Mariante A S and Grattapaglia D 2012 Genetic diversity and assessment of 23 microsatellite markers for parentage testing of Santa Inês hair sheep in Brazil. Genetics and Molecular Research 11(2):1217–1229.

Stevanović J, Stanimirović Z, Dimitrijević V, Stojić V, Fratrić N and Lazarević M 2009 Microsatellite DNA Polymorphism and its usefulness for pedigree verification in Simmental cattle from Serbia. Acta Veterinaria 59(5-6):621–631. http://www.actaveterinaria.rs/volume/issue/7/39/301

Untergrasser A, Cutcutache I, Koressaar T, Ye J, Faircloth B C, Remm M and Rozen S G 2012 Primer3 - New capabilities and interfaces. Nucleic Acids Res 40(15):e115. http://nar.oxfordjournals.org/content/40/15/e115.short


Received 7 June 2016; Accepted 7 July 2016; Published 1 August 2016

Go to top