Molecular genetic analysis of a cattle population to reconstitute the extinct Algarvia breed
© Ginja et al. 2010
Received: 13 January 2010
Accepted: 11 June 2010
Published: 11 June 2010
Skip to main content
© Ginja et al. 2010
Received: 13 January 2010
Accepted: 11 June 2010
Published: 11 June 2010
Decisions to initiate conservation programmes need to account for extant variability, diversity loss and cultural and economic aspects. Molecular markers were used to investigate if putative Algarvia animals could be identified for use as progenitors in a breeding programme to recover this nearly extinct breed.
46 individuals phenotypically representative of Algarvia cattle were genotyped for 27 microsatellite loci and compared with 11 Portuguese autochthonous and three imported breeds. Genetic distances and factorial correspondence analyses (FCA) were performed to investigate the relationship among Algarvia and related breeds. Assignment tests were done to identify representative individuals of the breed. Y chromosome and mtDNA analyses were used to further characterize Algarvia animals. Gene- and allelic-based conservation analyses were used to determine breed contributions to overall genetic diversity.
Genetic distance and FCA results confirmed the close relationship between Algarvia and southern Portuguese breeds. Assignment tests without breed information classified 17 Algarvia animals in this cluster with a high probability (q > 0.95). With breed information, 30 cows and three bulls were identified (q > 0.95) that could be used to reconstitute the Algarvia breed. Molecular and morphological results were concordant. These animals showed intermediate levels of genetic diversity (MNA = 6.0 ± 1.6, Rt = 5.7 ± 1.4, Ho = 0.63 ± 0.19 and He = 0.69 ± 0.10) relative to other Portuguese breeds. Evidence of inbreeding was also detected (Fis = 0.083, P < 0.001). The four Algarvia bulls had Y-haplotypes H6Y2 and H11Y2, common in Portuguese cattle. The mtDNA composition showed prevalence of T3 matrilines and presence of the African-derived T1a haplogroup. This analysis confirmed the genetic proximity of Algarvia and Garvonesa breeds (Fst = 0.028, P > 0.05). Algarvia cattle provide an intermediate contribution (CB = 6.18, CW = -0.06 and D1 = 0.50) to the overall gene diversity of Portuguese cattle. Algarvia and seven other autochthonous breeds made no contribution to the overall allelic diversity.
Molecular analyses complemented previous morphological findings to identify 33 animals that can be considered remnants of the Algarvia breed. Results of genetic diversity and conservation analyses provide objective information to establish a management program to reconstitute the Algarvia breed.
Breeding practices designed to alleviate production constraints are prejudicial to the survival of traditional domestic animal breeds, and tend to lead to impoverishment of the gene pool . The Food and Agriculture Organization of the United Nations has encouraged a series of conservation measures designed to help prevent irreversible loss of diversity of domesticated animal species . A heightened awareness of the cultural, historical and social heritage represented by traditional breeds has led to increased interest in their preservation . Despite its small geographic area, Portugal hosts a wide variety of domestic breeds , with as many as 13 autochthonous cattle breeds recognized [4, 5]. Analysis of genetic diversity of some of these breeds has used blood protein polymorphisms , microsatellite variation [7–12], mitochondrial DNA (mtDNA) [13, 14] and Y chromosome sequence variation . These studies have shown that, among European cattle, the Mirandesa breed is one of the most important targets for preservation based on its contribution to diversity . In general, southern European cattle appear to represent particularly important reservoirs of genetic diversity .
Algarvia cattle are native to the Algarve region of southern Portugal and were first represented in 1868  and later described with more details [18–23]. Based on its morphology, the breed has been classified in the Aquitanian (or Red Convex) group, together with Alentejana, Mertolenga, Garvonesa and Minhota [23, 24]. Except for Minhota, these breeds are distributed throughout the southern river Tagus valley. Algarvia cattle were used predominantly for meat production and/or draft. Although they were never formally registered as an independent breed, 25,000 - 29,000 animals were officially catalogued between 1940 and 1970 . Since then, the population has declined rapidly and the breed was considered to be effectively extinct by the 1980s , an event that subjected the Regional Agricultural Authorities to criticism for failure to implement proper conservation measures.
The morphology of eight autochthonous Portuguese cattle breeds has been analyzed by numerical taxonomic methods , and these analyses have been extended recently to include animals thought to derive from the Algarvia breed . This study identified a uniform group of animals that shared many of the phenotypic characteristics of that breed. These remnants were preserved mainly for cultural reasons by traditional breeders in different areas of Algarve and southern Alentejo regions.
In the present study, we have characterized these putative Algarvia animals using nuclear microsatellites, Y chromosome markers and mtDNA sequences to determine to which extent they can be genetically distinguished from related autochthonous breeds, and to identify individuals that would represent suitable progenitors for a breeding programme to reconstitute the Algarvia breed.
Forty-six animals (42 females and four males) were chosen from ten independent herds and consisted of individuals thought to derive from the Algarvia cattle on the basis of phenotypic similarity, as judged by traditional farmers with experience on the breed characteristics .
Additionally, samples from Garvonesa (29) and Preta (47) breeds were collected because they are present throughout the southern region of the country and admixture with Algarvia cannot be excluded. The herds in which putative Algarvia individuals were located also had animals from the Alentejana breed. A 9 ml whole blood sample was collected from each individual by jugular venipuncture in tubes containing EDTA-K3 as anticoagulant. Genomic DNA was extracted from leukocytes using the Puregene DNA Isolation Kit (Gentra Systems, Minneapolis, USA).
A set of 27 microsatellite markers (BM1818, BM1824, BM203, BM2113, BM2613, BRRIBO, CSSM36, CYP21, ETH10, ETH152, ETH225, ETH3, HEL11, HEL13, HEL9, ILSTS035, ILSTS065, INRA023, MGTG4B, RM006, RM067, SPS115, TGLA122, TGLA126, TGLA227, TGLA345 and TGLA53) was used to genotype each animal. Multiplex PCR were performed using fluorescence-labelled primers, as described by Mateus et al. . PCR products were separated by capillary electrophoresis on ABI PRISM® 310 instruments (Applied Biosystems, Foster City, CA) and fragment size analysis was done with STRand software . Genotypes from nine Portuguese autochthonous and three imported breeds were obtained from Mateus et al. [12, 29], thus the complete dataset included, in addition to the above breeds, Alentejana (50 individuals), Arouquesa (50), Barrosã (50), Brava de Lide (40), Minhota (50), Marinhoa (51), Maronesa (47), Mertolenga (50) and Mirandesa (50), Charolais (45), Friesian (35) and Limousin (48). Reference samples were included in all PCR assays to standardize allele sizing across datasets.
Allele frequencies were determined with GENALEX version 6 . The software GENEPOP version 3.4  was used to perform global and per locus/per population Hardy-Weinberg Equilibrium (HWE) tests, and to test for genotypic linkage disequilibrium (LD). Exact probability tests were done for loci with four or fewer alleles; otherwise, a Markov chain method was employed  with 10,000 dememorization steps, 500 batches and 5,000 iterations. GENETIX version 4.02  was used to estimate within-population observed (Ho) and unbiased expected (He) heterozygosities , the mean number of alleles (MNA), and inbreeding coefficients (Fis) . The statistical significance of Fis > 0 was obtained based on 1,000 permutations. Pairwise population Fst values were calculated with FSTAT version 2.9.3  and P-values obtained based on 1,000 randomizations. This software was also used to estimate allelic richness (Rt) per locus and population. To investigate breed relationships, neighbour-joining (N-J) dendrograms  were constructed from DA genetic distances  using POPULATIONS version 1.2.28 . Bootstrap values were obtained with 1,000 replicates over loci. A dendrogram based on allele-sharing distances between individuals was also constructed using this software. TREEVIEW version 1.6.6  was used to visualize and edit the dendrograms. A factorial correspondence analysis (FCA) was done with GENETIX to investigate relationships among individuals. Assignment tests were performed to identify individuals most representative of the Algarvia breed and to detect admixture. STRUCTURE version 2.2  was used to estimate the most probable number of population clusters (K). The analysis was done without prior information on populations, assuming correlated allele frequencies and admixture . Ten independent runs with 100,000 Markov Chain Monte Carlo (MCMC) iterations and 10,000 burn-in were performed at each K (1 ≤ K ≤ 9) to calculate ΔK as in Evanno et al. . A longer run (1,000,000 iterations and 100,000 burn-in) was done for the most probable K to determine the number of individuals within each cluster. An assignment test with these settings but including prior breed information was also performed. The partially Bayesian simulation-exclusion procedures of GENECLASS version 2.0  were used for assignment tests with 10,000 Monte Carlo resamplings of individuals [44–46].
The contribution of each population to the overall genetic diversity was analysed considering the within- (CW) and between-breed (CB) diversity components, and aggregate genetic diversity (D1 = Fst*CB + (1-Fst)*CW) as described by Ollivier and Foulley . METAPOP software version 1.0.2.  was used to account for allelic diversity and estimate the contribution of each population (ci) to a pool of maximal genetic diversity . Equal weights were given to within- and between-breed coancestries (λ = l). The average molecular coancestry (f m) of each population was also obtained with this software.
The four putative Algarvia males were genotyped for one SNP (UTY intron 19 AY936543: g.423C > A), one indel (ZFY intron 10 AF241271: g.697_8indelGT), and five microsatellites (DDX3Y_1, BM861, INRA189, UMN0103 and UMN0307) located in the male-specific region of the bovine Y chromosome. Analyses were done as described by Ginja et al.  for a comparison with previously identified patrilines in Portuguese autochthonous breeds.
A 919 bp PCR fragment containing the complete mtDNA control region was obtained and sequenced for Algarvia animals. The analysis was done as described by Ginja et al.  and sequences were aligned with the taurine reference sequence [GenBank: V00654, ] using the Multalin interface http://bioinfo.genotoul.fr/multalin/multalin.html. Haplotypes were identified with GENALEX  and ARLEQUIN version 2.0  was used to calculate haplotype diversity (H), nucleotide diversity (π), and the mean number of pairwise nucleotide differences (MNPD) accounting for heterogeneity of substitution rates per site . mtDNA sequences of the southern Portuguese autochthonous breeds and the imported Limousin were obtained from the GenBank database [accession numbers: FJ815445-59, FJ815525-40, FJ815573-88, FJ815620-35 and FJ815880-95] and used in ARLEQUIN to estimate pairwise-population Fst values (5% significance level obtained with 10,000 permutations). A Median-Joining (MJ) network of haplotypes was constructed with NETWORK version 4.5.10  software to investigate breed relationships.
Results of the Bayesian assignment tests done with STRUCTURE
Without prior information
With prior information
Q ± SD
% Corr. assign.
Q ± SD
% Corr. assign.
0.89 ± 0.14
0.97 ± 0.05
0.66 ± 0.40
0.85 ± 0.27
0.94 ± 0.07
0.97 ± 0.07
0.95 ± 0.05
0.98 ± 0.04
0.90 ± 0.15
0.97 ± 0.07
0.89 ± 0.18
0.94 ± 0.19
0.87 ± 0.22
0.95 ± 0.15
Results of GENECLASS assignments are summarized in Additional file 5 Table S2. Animals were correctly assigned if the genotype probability was higher than the threshold exclusively in their source population. For Algarvia animals, ~22 to 50% of the individuals were correctly assigned with accuracies > 0.94 (ratio between the number of correctly assigned individuals and the sum of correctly and incorrectly assigned). Depending on the threshold considered, between 2 and 20% of the Algarvia animals were excluded (genotype probabilities lower than the threshold in all populations) and 30 to 76% were assigned to multiple populations (genotype probabilities greater than the threshold in at least two populations).
The genetic diversity (MNA = 6.0 ± 1.6, Rt = 5.7 ± 1.4, Ho = 0.63 ± 0.19 and He = 0.69 ± 0.10) of the 33 Algarvia animals identified with STRUCTURE was identical to that found for the related Alentejana breed, and slightly lower than the average estimates across all breeds included in this study (Additional file 3 Table S1, MNA = 6.8 ± 0.6, Rt = 6.0 ± 1.8, Ho = 0.67 ± 0.05 and He = 0.70 ± 0.03). Ten population-specific alleles were found in the Algarvia population but none had a frequency greater than 0.05. Deviations from HWE were significant (P < 0.05) due to heterozygote deficit at ten loci (BM203, BM1824, BM2113, BRRIBO, ETH152, ILSTS035, SPS115, TGLA53, TGLA122, and TGLA345). Evidence of inbreeding within the Algarvia group could be inferred from the Fis estimate (0.083) which was significantly (P < 0.001) greater than zero. The molecular coancestry of Algarvia animals (f m = 0.322) was slightly higher than the overall value of 0.310. Pairwise population Fst estimates showed that Algarvia animals are genetically closer to the Alentejana breed (Fst = 0.045, P <0.05) than to the other southern breeds (results not shown).
Among the four Algarvia bulls, two (ALG3M and ALG4M) had the H11Y2 haplotype which is fixed in the Alentejana breed but also common in other Portuguese breeds, whereas the two other animals (ALG1M and ALG2M) had the H6Y2 haplotype which is fixed in the Garvonesa breed but also found in other Portuguese breeds as well as in Charolais and Limousin breeds .
Complete mtDNA control region sequences (909 bp) were obtained for the 33 Algarvia animals (four bulls and 29 females) assigned to this group with STRUCTURE (sequence quality of AG19 was low and thus was discarded). Sequence alignment is shown in Additional file 6 Figure S4 with polymorphic positions represented. A total of 12 distinct haplotypes was identified [GenBank: G086285-G086317] based on 21 variable sites, of which 11 were phylogenetically informative, nine were singletons, and one was an indel. The European T3 haplogroup (T at nt16255) was the most frequent (29 animals) but the African-derived T1a type (T at nt16050, C at nt16113 and C at nt16255) was also detected in three Algarvia individuals (AG16, AG40 and AG43). Genetic diversity estimates in the Algarvia population were H = 0.81 ± 0.05, π = 0.003 ± 0.002 and MNPD = 3.03 ± 1.62. Pairwise Fst values showed that Algarvia is significantly differentiated from all breeds except from Garvonesa (Fst = 0.028, P > 0.05).
Breed contributions to overall genetic diversity
ci, λ = 1
Brava de Lide
Decisions to initiate costly conservation programmes need to take into account assessment of extant variability and diversity loss, as well as cultural and economic aspects [57–59]. Algarvia cattle were adapted to the climatic and geographical conditions of Algarve, a region that is highly susceptible to suffer from climate change . In addition to its cultural and historical relevance, reconstituting this breed could also contribute to reinforce sustainable agriculture, a non negligible component of the economic activity of the region. The use of molecular data to assess the genetic structure of domestic species in conservation programmes has been described [59, 61–67]. Analysis of morphological and molecular information provides a more solid basis to define the characteristics of a breed than the use of morphological traits alone. DNA markers were not available at the time when the Algarvia breed became nearly extinct, and no biological samples have been preserved. Based on morphology descriptions , putative Algarvia descendants were identified . As preliminary work towards reconstituting this breed, we used molecular markers to investigate to which extent these animals could be distinguished from other Portuguese cattle. In agreement with their morphological classification [20, 23, 27, 68], genetic distances and factorial correspondence analyses showed a close relationship between Algarvia and other breeds of the Aquitanian (Red Convex) group. Algarvia and Garvonesa breeds are considered to be descendants of the Alentejana breed , and our study shows that Algarvia is closely related to these breeds. This genetic proximity probably also reflects recent admixture, because most of the Algarvia animals were found in herds predominantly composed of Alentejana animals. Genetic erosion due to crossbreeding is expected to have occurred since the breed began to decline about 35 years ago.
We used assignment tests to determine whether a cluster of putative Algarvia animals could be distinguished from the related southern breeds and to identify potential candidates to reconstitute the breed. STRUCTURE results confirmed that the most probable partition of the data agreed with the number of populations tested. Although the Algarvia group was the last to emerge as an independent cluster, it was possible to identify 17 cows that belonged to this group with q values > 0.95 without using prior information on populations. With prior information on sample origin, 33 animals (30 cows and three bulls) were classified as Algarvia. The partially Bayesian method of GENECLASS resulted in relatively low percentages of animals classified in each breed but, because of the high accuracy, it was useful to confirm reference animals in each population. Even though GENECLASS methods are conservative , the results were consistent with those of STRUCTURE with 23 animals assigned to Algarvia of which only two (AG34 and AG1M) were not among those selected by STRUCTURE.
Based on their genotypes, three bulls (AG2M, AG3M and AG4M) identified as Algarvia can be used to reconstitute the breed. The remaining bull (AG1M) was possibly admixed with Alentejana, according to the STRUCTURE analysis, although with GENECLASS it was assigned to the source population but with a low probability. Overall, a certain amount of convergence of results from independent approaches is noted, since 19 of these animals (16 cows and three bulls) were also identified as belonging to the core group to reconstitute the Algarvia breed based on numerical taxonomy analyses of morphological characters . Admixture was detected in several Algarvia females, three of which (AG33, AG35 and AG41) did not represent descendants of this breed because they were misclassified. Based on morphology, these animals clustered within a group also containing Alentejana, Garvonesa and Mertolenga animals . The heterogeneous composition of the putative Algarvia population is reflected by the lower average genotype membership coefficients when compared to those of breeds having herdbook registries. This result was not surprising given the expected dilution of Algarvia through crossbreeding.
The Y-haplotypes of the putative Algarvia bulls provided additional evidence of the genetic proximity with Alentejana and Garvonesa breeds, but did not exclude possible admixture with imported breeds such as Limousin or Charolais. mtDNA analysis corroborated the close genetic relationship between the core group of 33 Algarvia animals and the Garvonesa breed through their sharing of haplotypes. Common Iberian matrilines (European T3 and African T1a) were found in the Algarvia population, as well as a distinct haplotype (AG16) possibly related to the more ancestral African-derived AA haplogroup found in Creole cattle . The use of genotypic data, together with morphological analysis, facilitated the definition of a group of animals that could be used to reconstitute the Algarvia breed. How conservative the inclusion criteria should be relative to the acceptable degree of admixture will depend on breeding strategies yet to be defined.
The genetic diversity of the core set of 33 Algarvia animals was slightly lower than that found across Portuguese breeds. A significant heterozygote deficit was detected, possibly due to inbreeding, which is consistent with the strong genetic erosion. These results were not unexpected considering the extremely reduced number of extant Algarvia descendants. Another possible explanation for HWE deviations could be the sampling of Algarvia animals from independent herds, which could generate population subdivision and an increased frequency of homozygotes (e.g. Wahlund effect).
Gene diversity-based estimates indicate that the Algarvia population makes an intermediate contribution to the overall genetic diversity of Portuguese cattle. In contrast, the allelic diversity-based estimates suggest that the variation found in Algarvia is represented in the genetic pool of four other autochthonous breeds. The decision concerning which of these measures should be used in management programs to evaluate breed contributions to overall genetic diversity is not consensual . In the case of Portuguese autochthonous breeds that are considered endangered , particularly for the Algarvia breed, immediate conservation measures should aim at maximizing gene diversity rather than allelic diversity. This approach would maintain allelic diversity and guarantee a more effective response to selection while controlling inbreeding [61, 64].
Although nearly three decades have elapsed since the Algarvia breed was declared effectively extinct, we were able to identify a small group of cows and bulls with phenotypic characteristics of this breed. Analyses of autosomal, maternal and paternal markers have helped refine previous morphological findings to identify 33 animals that can be considered remnants from the Algarvia breed. For cultural and economic reasons, reconstituting the Algarvia breed is relevant to maintain the distinct regional identity of the Algarve. Molecular analyses have characterized the genetic diversity of the core set of animals which, together with information from conservation analyses, can be used to establish a management program to reconstitute the Algarvia breed.
This research was supported by DGV - PIDDAC (Direcção Geral de Veterinária - Programa de Investimentos e Despesas de Desenvolvimento da Administração Central) and the Fundação para a Ciência e a Tecnologia. PRAXIX XXI contract N. 3/3.2/CA/2005/95. C. Ginja was supported by a PhD fellowship (SFRH/BD/13502/2003) from the Fundação para a Ciência e a Tecnologia. We thank Angel Del Valle from the Veterinary Genetics Laboratory, University of California, Davis for assistance with the mtDNA sequencing.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.