Next Article in Journal
Assessing the Potential of Diverse Forage Mixtures to Reduce Enteric Methane Emissions In Vitro
Previous Article in Journal
The Addition of Nature Identical Flavorings Accelerated the Virucidal Effect of Pure Benzoic Acid against African Swine Fever Viral Contamination of Complete Feed
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Diversity and Population Genetic Structure of a Guzerá (Bos indicus) Meta-Population

by
Maria Gabriela C. D. Peixoto
1,*,
Maria Raquel S. Carvalho
2,
Andrea A. Egito
3,
Raphael S. Steinberg
2,
Frank Ângelo T. Bruneli
1,
Marco Antônio Machado
1,
Fernanda C. Santos
2,4,
Izinara C. Rosse
2,5 and
Pablo Augusto S. Fonseca
2,6
1
Embrapa Gado de Leite, Juiz de Fora 36038-330, Brazil
2
Departamento de Biologia Geral, Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais, Belo Horizonte 31270-901, Brazil
3
Embrapa Gado de Corte, Campo Grande 79106-550, Brazil
4
Centre for Addiction and Mental Health (CAMH), Campbell Family Mental Health Research Institute, Toronto, ON M5T 1R8, Canada
5
Departamento de Farmácia, Escola de Farmácia, Universidade Federal de Ouro Preto, Ouro Preto 35400-000, Brazil
6
Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON N1G 2W1, Canada
*
Author to whom correspondence should be addressed.
Animals 2021, 11(4), 1125; https://doi.org/10.3390/ani11041125
Submission received: 23 February 2021 / Revised: 25 March 2021 / Accepted: 2 April 2021 / Published: 14 April 2021
(This article belongs to the Section Ecology and Conservation)

Abstract

:

Simple Summary

The Guzerá breed is one of the most relevant animal genetic resources for tropical and extreme environments, that is, low and high temperatures, and dry or humid environments. These animals were selected for beef, milk, or dual-purpose (beef and milk), and were extensively used to produce crossbred animals. Consequently, subjecting the breed to intense genetic bottlenecks in Brazil. The local scientific community and breeders have pursued a sustainable management and conservation program over the last 50 years. To evaluate the efficacy of these efforts, we characterized genetic diversity and structure in a Guzerá meta-population. DNA samples of 744 animals from one dairy, nine dual-purpose, and five beef herds were genotyped for 21 microsatellite loci. The genetic diversity estimates suggested a low fixation index, revealing a large genetic diversity in Guzerá herds. The dual-purpose herds/selection lines are the most uniform subpopulation, while the beef one preserved larger amounts of genetic diversity, representing a genetic diversity reservoir for the breed. In addition, the dairy herd showed to be genetically distant from other herds. Taken together, these results suggest that this Guzerá meta-population has higher genetic diversity, a lower degree of population subdivision, and a lower inbreeding level.

Abstract

The Brazilian Guzerá population originated from a few founders introduced from India. These animals adapted well to the harsh environments in Brazil, were selected for beef, milk, or dual-purpose (beef and milk), and were extensively used to produce crossbred animals. Here, the impact of these historical events with regard to the population structure and genetic diversity in a Guzerá meta-population was evaluated. DNA samples of 744 animals (one dairy, nine dual-purpose, and five beef herds) were genotyped for 21 microsatellite loci. Ho, He, PIC, Fis, Fit, and Fst estimates were obtained considering either farms or lineages as subpopulations. Mean Ho (0.73) and PIC (0.75) suggest that genetic diversity was efficiently conserved. Fit, Fis and Fst values (95% CI) pointed to a low fixation index, and large genetic diversity: Fit (Farms = 0.021–0.100; lineages = 0.021–0.100), Fis (Farms = –0.007–0.076; lineages = −0.014–0.070), and Fst (Farms = 0.0237–0.032; lineages = 0.029–0.038). The dual-purpose herds/selection lines are the most uniform subpopulation, while the beef one preserved larger amounts of genetic diversity among herds. In addition, the dairy herd showed to be genetically distant from other herds. Taken together, these results suggest that this Guzerá meta-population has high genetic diversity, a low degree of population subdivision, and a low inbreeding level.

1. Introduction

Maximum conservation of genetic diversity is the best guarantee of long-time survival of a breed or a population under artificial selection. Significant loss of genetic variation due to Bulmer effect, inbreeding, and, especially in small populations, genetic drift in the populations under directional selection could imply in low response to genetic selection, decreasing genetic progress rate as well as adaptive performance [1,2,3]. On the other hand, the rescue of genetic variation by the strategy of introducing animals from other populations or other breeds is an attempt that could reduce the future genetic progress and disrupt the population phenotypic uniformity and adaptive performance. Many strategies were proposed to deal with this problem and the use of genomic tools has been presented as an opportunity to deal with sustainable genetic selection [4,5]. Therefore, monitoring genetic diversity is fundamental in any artificial selection process. The genetic conservation of local and adapted breeds is an important concern due to its contribution to current or future scenarios for genetic improvements and trait selection [5,6].
Guzerá (Bos indicus) cattle are originally from India, particularly from the dry lands of the Kankrej region. A relatively small population was brought from India at the end of the 19th century and expanded to Brazil, where it became adaptable to various, especially harsh, environmental conditions [7,8]. In Brazil, this breed has been selected for beef, milk, and, recently, dual-purpose (beef and milk) production. This breed remained prevalent in bovine livestock until the 1930s. By the middle of the last century, however, the Guzerá population decreased due to its extensive use in crossbreeding schemes [7,9]. The herdbook of the Guzerá breed was initiated in 1936 and closed in 1971. Subsequently, animals with an unknown pedigree but with the typical morphology of the Guzerá breed were registered in an “open herdbook”. A total of 211,107 animals were registered, with 84% in the closed herdbook. In 2013, 11% of males and 89% of females were registered, out of 17,295 animals born.
Analyzing genealogical data, Faria et al. (2009) [9] obtained an effective population size of 104 for all registered Guzerá individuals. A similar value (98) was calculated by Peixoto et al. (2010) [7] for Guzerá herds selected for milk purposes, suggesting genetic variability losses in this subpopulation.
In 1993, FAO included the Guzerá in the World List for Domestic Animals Diversity. This list includes breeds of economic, scientific, cultural, and agricultural interest to be conserved through management. Motivated by the importance of this breed as a genetic resource, two focused, interconnected breeding programs were initiated: one for milk traits (National Program for the Improvement of Guzerá for Milk Traits, PNMGuL) other for beef traits [10].
Nevertheless, because of the selection process based on genetic evaluation, further reduction in the number of sires used by breeders, due to the selection process based on genetic evaluation ranks, has become a concern. In the medium—and long-term, this could result in the reduced genetic variability of the population and increased levels of inbreeding [7,9]. From this perspective, better knowledge of the population structure and genetic variability within and among Guzerá subpopulations could be useful in future definitions of selection and mating strategies based on marker-assisted conservation [11,12,13,14].
Genetic diversity and population structure have been studied in bovine populations using microsatellite markers [15,16,17,18] or single nucleotide polymorphism (SNP) markers [19,20,21,22]. The genome of the Guzerá breed was published in 2016 and a large number of variants with putative functional impacts were described [23]. Moreover, the impact of the recent bottlenecks in the current genetic diversity estimates for the Guzerá were described using a genome-wide SNP genotyping approach [24]. However, the genetic composition within lineages, farms, and selection lines, as well as the genetic variability across the breed are still unclear, making the identification of genetic reservoirs and the implementation of genetic management harder. Microsatellite markers, specifically, have been largely used to estimate the magnitude of genetic variability in animal populations following FAO recommendations as well as population genetic structure [25].
The main objective of this study was to ascertain genetic diversity and population structure in a meta-population of the Guzerá breed. A meta-population was defined in the current study as the fraction of the population to which we had access, that is, the subpopulation made up of herds that take part in the breeding programs and among which some degree of gene flow occurs. From the results of this study, we expect to effectively contribute to the sustainability of this breed by minimizing the loss of diversity and allowing long-term genetic progress.

2. Materials and Methods

2.1. Sampling

This study used blood samples from 744 animals, representing around 10% of the animals in each of the 15 major Guzerá herds located in Minas Gerais state, in southeastern Brazil. The majority of herds are concentrated in the Minas Gerais State, which provides proven genetic material from Guzerá nationwide. Currently, some herds exploit Guzerá animals for dual-purpose (dairy and beef), using proven sires from both breeding programs. Other herds use Guzerá strictly for beef purposes (Supplementary Table S1). Currently, the breed is under a breeding scheme, and the use of AI, as well as the availability of sires with accurate genetic merit, has increased in the years since the release of the first sire summary for the breed [7,9].
A herd was defined as the group of animals present in a unique breeding location or farm. Many farms take part in the traditional breeding programs, owning herds composed of their own lineages or a mix composed of their own and external lineages from the whole Brazilian Guzerá livestock. Furthermore, the sampled herds have different selection objectives (selection lines), which were defined here as the group of animals selected for beef, dairy or dual-purpose. Two approaches were applied to understand the structure of this meta-population. In the first approach, each farm was considered a subpopulation; in the second approach, each lineage was considered a subpopulation. The comparison of genetic diversity among different lineages or farms was used to infer the population structure and gene flow in subpopulations.

2.2. DNA Extraction and Genotyping

Genomic DNA extraction, PCR amplification, and capillary electrophoresis were developed as described previously [26]. The microsatellites used in this study were chosen based on a set of loci, internationally recommended by both FAO and ISAG (International Society of Animal Genetics), to investigate the genetic diversity of animal species [27]. These microsatellites were identified from the bovine genetic linkage map available at MARC/USDA (Meat Animal Research Center/United States Department of Agriculture). A panel of 21 microsatellites (ILSTS093, MNB-208, BM1237, BMS2614, BM7169, RM150, BMS2252, JAB8, DIK5382, DIK4383, NRDIKM004, NLBCMK13, DIK4593, MNB-88, MNS-20, DIK5307, DIK5183, DIK4513, DIK2279, DIK5300, and DIK1143) was defined as the most informative, specifically for the Guzerá animals, based on their polymorphism information content (PIC).

2.3. Polymorphic Information Content and Average Exclusion Probabilities

Allele frequencies were extrapolated from direct counting. The number of alleles present in the population, numbers of individuals genotyped, expected (He) and observed heterozygosities (Ho), Polymorphism Information Content (PIC), and non-exclusion probabilities were calculated for each locus using the Cervus 3.0 [28] software and averaged for the 21 loci.

2.4. Global Hardy–Weinberg Tests

The Hardy–Weinberg equilibrium (HWE) was tested to support the assessment of the degree of allelic and genotypic diversity in the population. Departures from the HWE could indicate the impact of artificial selection (due to non-random mating) or reduced population size. Heterozygote deficiency was accessed using a multi-sample version of the Score test or U test with the GenePop 4.0 software [29]. The multi-sample score test of Rousset and Raymond (1995) is a global test across loci and/or samples performed using a Markov chain (MC) algorithm. The exact P-value of this test [30] was estimated through a 10,000 dememorization step, followed by 20 batches of 5000 iterations each. Two results were provided for each test by the MC algorithm: the estimated P-value associated with the null hypothesis of HWE, and the standard error (SE) of this estimate. The Hardy-Weinberg equilibrium was also tested using the Cervus 3.0 [28] software using a chi-square goodness-of-fit test.

2.5. F-Statistics and Genetic Distances

To estimate genetic differentiation and relationships among subpopulations, Wright (1978) [31] allele frequency-based correlations (individual fixation index (Fis), fixation index (Fst), overall fixation index (Fit)) were calculated following the variance-based method of Weir and Cockerham (1984) [32], using the FSTAT (version 2.9.3.2) program [33]. Differentiation was estimated by averaging the 21 locus single estimates of those F statistics. Estimates were obtained using two approaches. In the first approach, the 21 loci were analyzed by grouping the subpopulations according to the sample collection farms, comprising 744 animals in 15 farms. This approach was designed to assess the influence of the geographical proximity of farms on the exchange of animals or the migration between them. In the second approach, the subpopulations were grouped according to the lineage or origin of animals. In this approach, lineages with a small number of individuals were removed from the analysis, and the final sample comprised 664 animals from 15 lineages. Fis and pairwise Fst values were calculated for each subpopulation using the Molkin v3.0 software [34]. The AMOVA algorithm, available in the ARLEQUIN computational package, was used to evaluate the partitioning of the genetic variance between different sources of variation and to test their significance using a non-parametric permutation approach [35].
Pair-wise standard genetic distances (D) [36] were calculated from allele frequencies and were used to build a dendrogram, summarizing genetic relationships using a neighbor-joining method [36] implemented in Poptree2 [2]. Bootstrap re-sampling (n = 1000) was performed to evaluate the robustness of the clusters. Additionally, a Neighbor-Net graph (Bryant and Moulton, 2004), based on the Fst estimated for lineages as subpopulations, was constructed using the SplitsTree4 [37] program.
The program structure was used to obtain the cluster assignment of lineages through a Bayesian approach, implementing the Markov Chain method with increasing numbers of inferred populations [38]. The tests were conducted assuming the admixture model based on the correlations between allelic frequencies and applying burn-in periods of 50,000 and 300,000 iterations for each dataset. Two to 17 inferred clusters were performed with three independent runs each. The determination of the real number of selection lines (K) was carried out by implementing the Evanno method [39] to assign lineages to each selection line, supposing some differentiation between them. The graphic display was made using the DISTRUCT software [40].

2.6. Quantification of the Contribution to Diversity

The contribution of each lineage to meta-population genetic diversity was assessed using the method proposed by Caballero and Toro (2002) [41] with the use of the Molkin v3.0 software [26]. With this method, the criterion was the maintenance of the maximum overall Nei (1987) [36] genetic diversity (GD) in the set of lineages.

3. Results and Discussion

3.1. Performance of the Microsatellite Panel

The microsatellite panel was evaluated in terms of the number of alleles, He, Ho, PIC, and non-exclusion probabilities, presented in Supplementary Table S2. About 74% (15) of microsatellite loci were highly polymorphic. The number of alleles present in each locus ranged from 8 to 27. Five loci had a high allelic richness (≥18) and the other 10 showed several alleles ranging between 11 and 15. This provides evidence for the existence of substantial genetic variability in the Guzerá population. Combined results, including the 21 loci, are presented in Table 1.
The mean observed heterozygosity ranged from 0.47 to 0.85, diverging numerically from the expected values for some loci, mainly for DIK4513. The average Ho was 0.73, slightly below expected (0.77). Such differences are further discussed in light of the HWE tests. Polymorphic information content (PIC) ranged from 0.57 (BM1237) to 0.88 (DIK5183). Combined with the very small non-exclusion probabilities (Table 1), these results indicate the suitability of this set of markers for studying genetic diversity in the Guzerá breed.

3.2. Global Hardy–Weinberg Equilibrium Tests

Considering the meta-population, significant (p < 0.05) results were obtained from the HWE test for each locus (Supplementary Table S2). Heterozygote deficiency was observed in nine loci (BM7169, DIK5382, DIK4383, NLBCMK13, DIK4593, DIK5307, DIK5183, DIK4513, and DIK2279). Among the possible causes are non-random mating, such as inbreeding or unidirectional gene flow caused by migration, mainly through artificial insemination (AI); or the linkage disequilibrium between the microsatellite and a locus under selection [42,43].

3.3. F-Statistics and Genetic Distances

Table 2 presents an average of 21 loci F-statistics [32] for both approaches, with the subpopulations considered according to farm or lineage. Fit values suggest a low fixation index, revealing a large amount of genetic diversity in the Guzerá herds sampled. The Fst and Fis values obtained in this study were lower than those reported previously in the pedigree data analysis of Guzerá animals registered in the 1994–1998 period: Fst = 1.75 and Fis = 1.37 [9]. These differences could be attributed to differences in methodology between the studies. However, despite the Fis and Fst values found in this study being below 0.05, it must be highlighted that there are situations in which Fst estimates lower than 0.05 do not necessarily imply negligible genetic differentiation [31]. These values support the hypothesis that these subpopulations have been going through a process of genetic differentiation in the approaches used. One hypothesis to explain the differences in the Fst would be that, before AI, most farmers tended to use their own or few lineages preferentially, which resulted in lower genetic diversity. This process was reverted by the introduction of AI as well as by the availability of sires of accurate genetic merit, from other herds or lineages, since the release of the first sire summary with the consequent higher gene flow (migration) between herds [7,9].
The approach considering lineages as subpopulations showed the impact of the distinct lineages in the population structure. The Fis was slightly lower (−0.007) and the Fst was slightly higher (+0.007) than those for farms. The Fst for the lineage approach was 0.034, compared to 0.027 for subpopulations grouped according to farm (Table 2). Table 3 presents the Fis values for each farm subpopulation and pairwise Fst, as well as the genetic distances between farms and lineages. These values are displayed graphically in Figure 1.
Estimated values for Fis using the farm approach ranged from −0.068 to 0.077. This indicates that genetic variability within subpopulations varied among farms, and these differences allow the reconstruction of historical aspects of the Guzerá breed in Brazil. Some farms developed their own lineages and kept them closed over time, avoiding the use of animals from other lineages in their breeding. FARM6, for example, developed LINE7. It is the only farm selecting exclusively for milk production and has been completely closed since the beginning of the 20th century. Nonetheless, LINE7, being one of the oldest, contributed to all the dual-purpose farms and even to some of the beef farms. Out of the 15 lineages, 5 were selected only for beef traits, each one coming from a specific farm: FARM1 (LINE4), FARM2 (LINE6), FARM3 (LINE10), FARM4 (LINE1), and FARM5 (LINE12).
Farms selecting for dual-purpose were originally beef farms or newly formed. For example, FARM11 (LINE9) was also closed for many years and, only from the end of the 20th century, started to use bulls of different lineages. On the contrary, FARM7 and FARM15 were recently formed by animals of different lineages. In the dual-purpose selection group, the correspondence between farms and lineages was far less strict. These historical aspects affected both the genetic diversity, Fis, Fst, and genetic distances.
For example, some farms (FARM6, FARM10, and FARM14) showed lower levels of heterozygosis, probably resulting from the intensive or exclusive use of their own sires or a small number of sires from few origins. At the same time, other farms (FARM1, FARM7, FARM11, and FARM15) showed high heterozygosis, attributable to the broad use of planned mating and/or semen from animals of different origins or lineages.
Estimated values for Fis using the lineage approach ranged from −0.109 to 0.061. This indicates that genetic diversity also varied among lineages. Some lineages were developed as part of the selection strategy of farms. In those cases, higher Fis were observed for both the farms and the lineages developed within them. For example, some lineages (LINE3, LINE6, and LINE7) showed higher Fis and lower heterozygosity, probably resulting again from the intensive or exclusive use of their own sires or a small number of sires from few origins.
Other lineages (LINE8, LINE9, and LINE14) showed high heterozygosity. LINE8 originated from LINE7, hence the high genetic similarity between them, and it has also been composed of other lineages resulting in its high heterozygosity. LINE14 came from FARM15, which contains a herd recently formed from a large genetic basis, hence its high degree of heterozygosity.
The pairwise Fst estimates are shown in Table 3. Pairwise Fst ranged from 0.0029 (FARM15 × FARM1) to 0.048 (FARM6 × FARM3) using the farm approach. They ranged more widely, from 0.007 (LINE8 × LINE11) to 0.087 (LINE7 × LINE10), using the lineage approach. In general, higher average Fst values were estimated for beef lineages. This means that less gene flow has occurred (pairwise Fst), both in the meta-population and in the whole (average Fst). This contrasts with the higher gene flow among dual-purpose lineages expressed by the lower average Fst estimates. One exception can be observed regarding LINE6 (FARM2), which has the lowest pairwise Fst average among the beef lineages (0.018 in the lineage and 0.015 in the farm approaches). Guzerá beef cattle are more extensively bred, generally using little artificial insemination and thus less exchange of genetic material between herds or lineages. However, considering the results of this study, beef herds kept a larger amount of the genetic diversity in Guzerá breeds.
In general, the lower average Fst estimates for the dual-purpose selection lines reveal high genetic variability. This is probably a result of the wide genetic basis for these subpopulations. One clear exception was observed among the dual-purpose selection lines: the average Fst estimated for LINE7 was 0.052, the highest value among all the subpopulations considering dairy, dual-purpose, and beef selection lines. This subpopulation was closed for many generations, and this Fst value reinforces the available information on the breeding policy for this lineage [44]. It should be emphasized that this lineage was selected from one of the first Guzerá animals introduced into Brazil in the early 20th century. Fixation index values were less variable among dual-purpose herds, which indicates that gene flow was greater among farms than among lineages. Thus, some fraction of these herds is most likely kept closed.
Nei’s standard genetic distances ranged from 0.05 (FARM8 ×FARM6, or FARM11 × FARM7) to 0.031 (FARM9 × FARM3, or FARM6 × FARM3) for the farm grouping approach (Figure 2a). In the lineage approach (Figure 2b), these distances ranged from 0.05 (LINE8 × LINE13) to 0.53 (LINE7 × LINE10). On average, LINE7 was placed on a separate branch of the dendrogram (Figure 2b). It is the most genetically distant subpopulation, despite having contributed to all dual-purpose subpopulations. It was especially distant from LINE4 and LINE10, both of which were selected for beef purposes. The values obtained for the genetic distances were generally in agreement with the F-statistics results, validating the conclusions about population structure and gene flow observed among dual-purpose herds (Figure 2).
In the dendrogram, LINE4, LINE6, and LINE11 were clustered to the same branch. Indeed, LINE11 was derived mainly from LINE4, LINE8, and LINE6 (Figure 1b). In addition, LINE4 and LINE6 were derived from geographically close farms that exchange some genetic material. The same has occurred with LINE9, LINE1, and LINE14. LINE9 contributed to LINE1, LINE14, LINE12, and LINE5. In addition, LINE12 also contributed to LINE5 and they both received contributions from LINE7. LINE3 introduced genetic material from LINE15. LINE10 and LINE2, which are geographically close, exchanged genetic material more frequently. Therefore, the dendrogram shown in Figure 2 reflects the historical register of the Guzerá lineage formation.
Table 4 shows the results of some studies that evaluated the genetic diversity of bovine breeds. These studies were selected based on two criteria: publication after 2000 and sample size greater than 100. There is an interesting relationship between the sample size and the number of microsatellite markers used in the present study. Despite this comparison, it is important to highlight the differences in experimental design, analysis, and result presentation/interpretation. Consequently, the comparison among the results must be interpreted carefully. The 744 animal sample size is larger than the samples of all studies except one [45]. However, the number of microsatellite markers used in that study was smaller than the number used in the present study. The Fst values obtained using the farm (0.025) and lineage (0.034) subpopulations were lower than almost all the results shown in Table 4. Comparing heterozygosity, the values obtained for the Guzerá breed (Ho = 0.73 and He = 0.77) in the present study were generally greater than those obtained for all breeds used in the studies shown in Table 4. These results suggest that the Guzerá breed has a greater genetic diversity, a lower degree of population subdivision, and a lower inbreeding level than the breeds analyzed in these other studies. Interestingly, the results obtained in the current study are similar to those obtained using microsatellites and single nucleotide polymorphisms in other Brazilian breeds [15,22].

3.4. Quantifying Contributions to Diversity

The contribution of each subpopulation to the overall genetic diversity is presented in Table 5. Some lineages contributed more to the increase of overall diversity because of their internal diversity. This group includes LINE15, LINE11, LINE6, and LINE8. On the other hand, some lineages had very low internal diversity, decreasing the overall genetic diversity. LINE7 is an example of a subpopulation with very low internal diversity, presenting the largest loss of genetic diversity (Loss-gain ratio: 0.1056). This herd was established using a small number of animals and, since then, has remained closed. Additionally, it is the only herd that did not have genetic flow from other farms and/or subpopulations. The mean distance between subpopulations is another feature contributing to the overall genetic diversity. Thus, LINE7 and LINE10 could be contributing to the increase in the mean distance between subpopulations, since they are the most closed lineages.

3.5. Numbers of Clusters and Lineages

An interesting result was obtained with respect to meta-population structure (Figure 3). If lineages were used as subpopulations, K = 3 was the most probable one, separating the beef (mainly LINE4 and LINE10) from the dual-purpose (mainly LINE3, LINE5, LINE8, LINE9, LINE13, LINE14, and LINE15) and dairy (LINE7) lineages. LINE7 is the only exclusively dairy lineage and it was used to generate all the dual-purpose cattle. The dual-purpose lineages, shown in Figure 3, have contributions from the green component, which corresponds to the dairy lineages, reinforcing the historical contribution of LINE7 to the formation of the dual-purpose. Beef cattle were separated into two subpopulations: one, consisting of lineages from closed herds with little gene flow; the other, consisting of lineages from less closed herds that strongly contributed to the dual-purpose herds (LINE4 and LINE10, respectively). The inclusion of LINE8 and LINE11 in one of the dual-purpose cattle clusters was also interesting. LINE11 was sequentially derived from LINE8, which was derived from LINE7. LINE11 was among the dual-purpose cattle for a period after which selection pressure was imposed for milk trait improvement but only for beef purposes, or better, for improving the developmental and weight trait performance. LINE11, LINE5, and LINE14 were all grouped in both beef and dual-purpose clusters, signaling that beef traits are also very important for meeting their selection objectives. Moreover, some dual-purpose lineages showed a strong contribution to the dairy lineage (LINE7) genetic diversity (LINE3, LINE13, LINE14, and LINE15). Therefore, indicating that these lineages can be used as genetic reservoirs for dairy genetic improvement purposes.

4. Conclusions

Despite the historical events that the Guzerá breed went through in the last decades, such as a strong founder effect and several population bottlenecks, the breed nevertheless preserved a good amount of genetic diversity. In the cluster analysis, it was possible to recover both the origin of the lineages and the selection purposes. Reasonably, in the best model (K = 3), dual-purpose selection lines clustered with the beef selection lines from which they originated. Also, gene flow between the different lineages was detected in the clusters formed. Modeling the genetic structure of the Guzerá meta-population helped to understand the impact of planned crosses on the genetic diversity and the distribution of the genetic diversity among lineages selected for different purposes (beef, dairy, and dual-purpose). The separation of animals for a new selection purpose, if not implying a stronger selection pressure, decreases the risk of loss of diversity. Taken together, these results suggest that, despite the population structure produced by selection strategies based originally on relatively closed farms/lineages, the amount of genetic diversity conserved in the Guzerá breed is high. In particular, beef farms/lineages emerge as a reservoir for both genetic selection and conservation programs.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ani11041125/s1, Table S1: Sample distribution according to selection purpose and farms; Table S2: Microsatellite loci description.

Author Contributions

Conceptualization, M.G.C.D.P., M.R.S.C., F.Â.T.B., A.A.E.; Methodology, P.A.S.F., F.C.S., I.C.R., R.S.S., M.A.M.; Software, P.A.S.F., F.C.S., I.C.R.; Validation, P.A.S.F., M.G.C.D.P.; Formal Analysis, P.A.S.F., M.G.C.D.P., M.R.S.C., F.C.S.; Investigation, M.G.C.D.P., P.A.S.F., M.R.S.C.; Data Curation, M.G.C.D.P., P.A.S.F.; Writing—Original Draft Preparation, Formal Analysis, P.A.S.F., F.C.S., M.G.C.D.P., M.R.S.C.; A.A.E.; Writing—Review and Editing, M.G.C.D.P., M.R.S.C., P.A.S.F.; Visualization, P.A.S.F.; Supervision, M.G.C.D.P., M.R.S.C., P.A.S.F.; Project Administration, M.G.C.D.P.; Funding Acquisition, M.G.C.D.P., M.R.S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work funded by Empresa Brasileira de Pesquisa Agropecuária (EMBRAPA) with the following two projects: 03.08.00.080.00.00 and 02.13.05.010.00.00. MRSC was awarded a fellowship from the Brazilian National Research Council (Conselho Nacional de Pesquisa, CNPq, grant number 312068/2015-8 and 312405/2018-9) and is supported by grants from the Fundação de Amparo à Pesquisa de Minas Gerais (FAPEMIG) APQ-01093-15, APQ-02003-15, APQ-01377-17; by the Instituto Nacional de Ciência e Tecnologia de Ciência Animal; Pró-Reitoria de Pesquisa de Universidade Federal de Minas Gerais (PRPq/UFMG).

Institutional Review Board Statement

Ethical review and approval were waived for this study, due to the fact that only DNA samples stored in the genetic biobank of the research groups involved in the current study were used.

Acknowledgments

Authors thank the Brazilian Guzerá breeders and the Centro Brasileiro de Melhoramento Genético do Guzerá (CBMG2) for their logistic support and for providing biological materials, FAPEMIG for financial support, as well as the technicians of the Laboratory of Molecular Genetics of the Embrapa Dairy Cattle for DNA extraction and genotyping, and the creators of the free software used for the analysis performed in the current study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. de Rochambeau, H.; Fournet-Hanocq, F.; Vu Tien Khang, J. Measuring and managing genetic variability in small populations. Ann. Zootech. 2000, 49, 77–93. [Google Scholar] [CrossRef]
  2. Taberlet, P.; Valentini, A.; Rezaei, H.R.; Naderi, S.; Pompanon, F.; Negrini, R.; Ajmone-Marsan, P. Are cattle, sheep, and goats endangered species? Mol. Ecol. 2008, 17, 275–284. [Google Scholar] [CrossRef] [PubMed]
  3. Kristensen, T.N.; Hoffmann, A.A.; Pertoldi, C.; Stronen, A.V. What can livestock breeders learn from conservation genetics and vice versa? Front. Genet. 2015, 5, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Boichard, D.; Ducrocq, V.; Fritz, S. Sustainable dairy cattle selection in the genomic era. J. Anim. Breed. Genet. 2015, 132, 135–143. [Google Scholar] [CrossRef]
  5. Toro, M.A.; Fernández, J.; Caballero, A. Molecular characterization of breeds and its use in conservation. Livest. Sci. 2009, 120, 174–195. [Google Scholar] [CrossRef]
  6. In Vivo Conservation of Animal Genetic Resources; FAO: Rome, Italy, 2013; ISBN 0306-5251.
  7. Peixoto, M.G.C.D.; Poggian, C.F.; Verneque, R.S.; Egito, A.A.; Carvalho, M.R.S.; Penna, V.M.; Bergmann, J.A.G.; Viccini, L.F.; Machado, M.A. Genetic basis and inbreeding in the Brazilian Guzerat (Bos indicus) subpopulation selected for milk production. Livest. Sci. 2010, 131, 168–174. [Google Scholar] [CrossRef]
  8. Santana, M.L.; Pereira, R.J.; Bignardi, A.B.; El Faro, L.; Pires, M.F.Á.; Andrade, R.G.; Perez, B.C.; Bruneli, F.A.T.; Peixoto, M.G.C.D. Dual-purpose Guzerá cattle exhibit high dairy performance under heat stress. J. Anim. Breed. Genet. 2020, 137, 486–494. [Google Scholar] [CrossRef] [PubMed]
  9. Faria, F.J.C.; Filho, A.E.V.; Madalena, F.E.; Josahkian, L.A. Pedigree analysis in the Brazilian Zebu breeds. J. Anim. Breed. Genet. 2009, 126, 148–153. [Google Scholar] [CrossRef] [PubMed]
  10. Peixoto, M.G.C.D.; Verneque, R.S.; Teodoro, R.L.; Penna, V.M.; Martinez, M.L. Genetic trend for milk yield in Guzerat herds participating in progeny testing and MOET nucleus schemes. Genet. Mol. Res. 2006, 5, 454–465. [Google Scholar] [PubMed]
  11. Lenstra, J.A. Marker-assisted conservation of European cattle breeds: An evaluation. Anim. Genet. 2006, 37, 475–481. [Google Scholar]
  12. Shivaji, S.; Kholkute, S.D.; Verma, S.K.; Gaur, A.; Umapathy, G.; Singh, A.; Sontakke, S.; Shailaja, K.; Reddy, A.; Monika, S.; et al. Conservation of wild animals by assisted reproduction and molecular marker technology. Indian J. Exp. Biol. 2003, 41, 710–723. [Google Scholar] [PubMed]
  13. Allendorf, F.W.; Hohenlohe, P.A.; Luikart, G. Genomics and the future of conservation genetics. Nat. Rev. Genet. 2010, 11, 697–709. [Google Scholar] [CrossRef] [PubMed]
  14. Eusebi, P.G.; Martinez, A.; Cortes, O. Genomic tools for effective conservation of livestock breed diversity. Diversity 2020, 12, 8. [Google Scholar] [CrossRef] [Green Version]
  15. Egito, A.A.; Paiva, S.R.; Albuquerque, M.D.S.M.; Mariante, A.S.; Almeida, L.D.; Castro, S.R.; Grattapaglia, D. Microsatellite based genetic diversity and relationships among ten Creole and commercial cattle breeds raised in Brazil. BMC Genet. 2007, 8, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Dadi, H.; Tibbo, M.; Takahashi, Y.; Nomura, K.; Hanada, H.; Amano, T. Microsatellite analysis reveals high genetic diversity but low genetic structure in Ethiopian indigenous cattle populations. Anim. Genet. 2008, 39, 425–431. [Google Scholar] [CrossRef]
  17. Dalvit, C.; De Marchi, M.; Dal Zotto, R.; Zanetti, E.; Meuwissen, T.; Cassandro, M. Genetic characterization of the Burlina cattle breed using microsatellites markers. J. Anim. Breed. Genet. 2008, 125, 137–144. [Google Scholar] [CrossRef] [PubMed]
  18. Groeneveld, L.F.; Lenstra, J.A.; Eding, H.; Toro, M.A.; Scherf, B.; Pilling, D.; Negrini, R.; Finlay, E.K.; Jianlin, H.; Groeneveld, E.; et al. Genetic diversity in farm animals—A review. Anim. Genet. 2010, 41, 6–31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. McKay, S.D.; Schnabel, R.D.; Murdoch, B.M.; Matukumalli, L.K.; Aerts, J.; Coppieters, W.; Crews, D.; Neto, E.D.; Gill, C.A.; Gao, C.; et al. An assessment of population structure in eight breeds of cattle using a whole genome SNP panel. BMC Genet. 2008, 9, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Flury, C.; Tapio, M.; Sonstegard, T.; Drögemüller, C.; Leeb, T.; Simianer, H.; Hanotte, O.; Rieder, S. Effective population size of an indigenous Swiss cattle breed estimated from linkage disequilibrium. J. Anim. Breed. Genet. 2010, 127, 339–347. [Google Scholar] [CrossRef]
  21. Sanz, A.; Martin-Burriel, I.; Cons, C.; Reta, M.; Poblador, A.; Rodellar, C.; Zaragoza, P. Genetic diversity, structure and individual assignment of Casta Navarra cattle: A well-differentiated fighting bull population. J. Anim. Breed. Genet. 2014, 131, 11–18. [Google Scholar] [CrossRef] [PubMed]
  22. Campos, B.M.; do Carmo, A.S.; do Egito, A.A.; da Mariante, A.S.; do Albuquerque, M.S.M.; de Gouveia, J.J.S.; Malhado, C.H.M.; Verardo, L.L.; da Silva, M.V.G.B.; Carneiro, P.L.S. Genetic diversity, population structure, and correlations between locally adapted zebu and taurine breeds in Brazil using SNP markers. Trop. Anim. Health Prod. 2017, 49, 1677–1684. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Rosse, I.C.; Assis, J.G.; Oliveira, F.S.; Leite, L.R.; Araujo, F.; Zerlotini, A.; Volpini, A.; Dominitini, A.J.; Lopes, B.C.; Arbex, W.A.; et al. Whole genome sequencing of Guzerá cattle reveals genetic variants in candidate genes for production, disease resistance, and heat tolerance. Mamm. Genome 2017, 28, 66–80. [Google Scholar] [CrossRef]
  24. de Souza Fonseca, P.A.; dos Santos, F.C.; Rosse, I.C.; Ventura, R.V.; Brunelli, F.Â.T.; Penna, V.M.; da Silva Verneque, R.; Machado, M.A.; da Silva, M.V.G.B.; Carvalho, M.R.S.; et al. Retelling the recent evolution of genetic diversity for Guzerá: Inferences from LD decay, runs of homozygosity and Ne over the generations. Livest. Sci. 2016, 193, 110–117. [Google Scholar] [CrossRef]
  25. Thomson, K. World Agriculture: Towards 2015/2030: An FAO Perspective; FAO: Rome, Italy, 2003. [Google Scholar]
  26. Machado, M.A.; Schuster, I.; Martinez, M.L.; Campos, A.L. Genetic Diversity of Four Cattle Breeds Using Microsatellite Markers. Rev. Bras. Zootec. 2003, 32, 93–98. [Google Scholar] [CrossRef] [Green Version]
  27. Ruttan, V.W.; Alexandratos, N. World Agriculture: Towards 2010: An FAO Study; FAO: Rome, Italy, 1996. [Google Scholar]
  28. Kalinowski, S.T.; Taper, M.L.; Marshall, T.C. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol. Ecol. 2007, 16, 1099–1106. [Google Scholar] [CrossRef] [PubMed]
  29. Rousset, F. GENEPOP’007: A complete re-implementation of the GENEPOP software for Windows and Linux. Mol. Ecol. Resour. 2008, 8, 103–106. [Google Scholar] [CrossRef]
  30. Guo, S.W.; Thompson, E.A. Performing the Exact Test of Hardy-Weinberg Proportion for Multiple Alleles. Biometrics 1992, 48, 361–372. [Google Scholar] [CrossRef] [PubMed]
  31. Wright, S. Variability within and among natural populations. In Evolution and the Genetics of Populations; University of Chicago Press: Chicago, IL, USA, 1978. [Google Scholar]
  32. Weir, B.S.; Cockerham, C.C. Estimating F-statistics for the analysis of population structure. Evolution 1984, 36, 1358–1370. [Google Scholar]
  33. Goudet, J. FSTAT (Version 1.2): A Computer Program to Calculate F-statistics. J. Hered. 1995, 86, 485–486. [Google Scholar] [CrossRef]
  34. Gutiérrez, J.P.; Royo, L.J.; Álvarez, I.; Goyache, F. MolKin v2.0: A computer program for genetic analysis of populations using molecular coancestry information. J. Hered. 2005, 96, 718–721. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Excoffier, L.; Laval, G.; Schneider, S. Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol. Bioinform. 2005, 1, 47–50. [Google Scholar] [CrossRef] [Green Version]
  36. Saitou, N.; Nei, M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 1987, 4, 406–425. [Google Scholar] [PubMed]
  37. Huson, D.H.; Bryant, D. Estimating Phylogenetic Trees and Networks Using SplitsTree4. 2014. Available online: www.splitstree.org (accessed on 23 February 2021).
  38. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar] [PubMed]
  39. Evanno, G.; Regnaut, S.; Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol. Ecol. 2005, 14, 2611–2620. [Google Scholar] [CrossRef] [Green Version]
  40. Rosenberg, N.A. DISTRUCT: A program for the graphical display of population structure. Mol. Ecol. Notes 2004, 4, 137–138. [Google Scholar] [CrossRef]
  41. Caballero, A.; Toro, M.A. Analysis of genetic diversity for the management of conserved subdivided populations. Conserv. Genet. 2002, 3, 289–299. [Google Scholar] [CrossRef]
  42. Chen, B.; Cole, J.W.; Grond-Ginsbach, C. Departure from Hardy Weinberg Equilibrium and genotyping error. Front. Genet. 2017. [Google Scholar] [CrossRef] [PubMed]
  43. Garnier-Géré, P.; Chikhi, L. Population Subdivision, Hardy-Weinberg Equilibrium and the Wahlund Effect. In eLS; John Wiley & Sons: Hoboken, NJ, USA, 2013. [Google Scholar]
  44. Panetto, J.C.C.; Gutiérrez, J.P.; Ferraz, J.B.S.; Cunha, D.G.; Golden, B.L. Assessment of inbreeding depression in a Guzerat dairy herd: Effects of individual increase in inbreeding coefficients on production and reproduction. J. Dairy Sci. 2010, 93, 4902–4912. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Delgado, J.V.; Martínez, A.M.; Acosta, A.; Álvarez, L.A.; Armstrong, E.; Camacho, E.; Cañõn, J.; Cortés, O.; Dunner, S.; Landi, V.; et al. Genetic characterization of Latin-American Creole cattle using microsatellite markers. Anim. Genet. 2012, 43, 2–10. [Google Scholar] [CrossRef]
  46. Kim, K.S.; Yeo, J.S.; Choi, C.B. Genetic diversity of north-east Asian cattle based on microsatellite data. Anim. Genet. 2002, 33, 201–204. [Google Scholar] [CrossRef] [PubMed]
  47. Mukesh, M.; Sodhi, M.; Bhatia, S.; Mishra, B.P. Genetic diversity of Indian native cattle breeds as analysed with 20 microsatellite loci. J. Anim. Breed. Genet. 2004, 121, 416–424. [Google Scholar] [CrossRef]
  48. Sasazaki, S.; Honda, T.; Fukushima, M.; Oyama, K.; Mannen, H.; Mukai, F.; Tsuji, S. Genealogical relationship between pedigree and microsatellite information and analysis of genetic structure of a highly inbred Japanese black cattle strain. Asian Australas. J. Anim. Sci. 2004, 17, 1355–1359. [Google Scholar] [CrossRef]
  49. Chikhi, L.; Goossens, B.; Treanor, A.; Bruford, M.W. Population genetic structure of and inbreeding in an insular cattle breed, the Jersey, and its implications for genetic resource management. Heredity 2004, 92, 396–401. [Google Scholar] [CrossRef]
  50. Maretto, F.; Cassandro, M. Temporal variation in genetic diversity and population structure of Burlina cattle breed. Ital. J. Anim. Sci. 2014. [Google Scholar] [CrossRef] [Green Version]
  51. Rendo, F.; Iriondo, M.; Jugo, B.M.; Aguirre, A.; Mazón, L.I.; Vicario, A.; Gómez, M.; Estonba, A. Analysis of the genetic structure of endangered bovine breeds from the Western Pyrenees using DNA microsatellite markers. Biochem. Genet. 2004, 42, 99–108. [Google Scholar] [CrossRef] [PubMed]
  52. Pham, L.D.; Do, D.N.; Binh, N.T.; Nam, L.Q.; Van Ba, N.; Thuy, T.T.T.; Hoan, T.X.; Cuong, V.C.; Kadarmideen, H.N. Assessment of genetic diversity and population structure of Vietnamese indigenous cattle populations by microsatellites. Livest. Sci. 2013, 155, 17–22. [Google Scholar] [CrossRef]
Figure 1. Dendrogram summarizing genetic relationships among Guzerá farms (a) and selection lines (b), using lineages Fst (F distances) generated by Poptree2 for the Guzerá meta-population in Minas Gerais state, Brazil.
Figure 1. Dendrogram summarizing genetic relationships among Guzerá farms (a) and selection lines (b), using lineages Fst (F distances) generated by Poptree2 for the Guzerá meta-population in Minas Gerais state, Brazil.
Animals 11 01125 g001
Figure 2. Dendrogram summarizing genetic relationships among Guzerá farms (a) and lineages (b), using Nei’s standard genetic distances generated by Poptree2 for the Guzerá meta-population in Minas Gerais state, Brazil.
Figure 2. Dendrogram summarizing genetic relationships among Guzerá farms (a) and lineages (b), using Nei’s standard genetic distances generated by Poptree2 for the Guzerá meta-population in Minas Gerais state, Brazil.
Animals 11 01125 g002
Figure 3. Histogram showing the population structure for the 15 lineages of the Guzerá cattle meta-population in Minas Gerais state, Brazil, generated using the model-based STRUCTURE software. Each animal is represented by a single vertical line divided into K colors, where K is the number of clusters assumed and the colors show the estimated individual proportions of cluster membership. K = 3: cluster 1—red (beef cattle); cluster 2—dark blue (beef cattle); and cluster 3—green (dual-purpose).
Figure 3. Histogram showing the population structure for the 15 lineages of the Guzerá cattle meta-population in Minas Gerais state, Brazil, generated using the model-based STRUCTURE software. Each animal is represented by a single vertical line divided into K colors, where K is the number of clusters assumed and the colors show the estimated individual proportions of cluster membership. K = 3: cluster 1—red (beef cattle); cluster 2—dark blue (beef cattle); and cluster 3—green (dual-purpose).
Animals 11 01125 g003
Table 1. Combined analysis of 21 microsatellite loci of the Guzerá cattle at Minas Gerais, Brazil.
Table 1. Combined analysis of 21 microsatellite loci of the Guzerá cattle at Minas Gerais, Brazil.
Number of individuals:723
Number of loci:21
Mean number of alleles per locus:13.67
Mean proportion of individuals typed:0.9187
Mean expected heterozygosity:0.7746
Mean polymorphic information content (PIC):0.7473
Combined non-exclusion probability (first parent):1.17 × 10−5
Combined non-exclusion probability (second parent):5.64 × 10−9
Combined non-exclusion probability (parent pair):8.41 × 10−15
Combined non-exclusion probability (identity):2.16 × 10−24
Combined non-exclusion probability (sib identity):1.63 × 10−9
Table 2. F-statistics (Fis, Fst, Fit) means, standard errors (SE), and 95% confidence intervals obtained for the Guzerá dairy metapopulation, with the subpopulations defined as farm sampled (on the left) or as the lineage (on the right).
Table 2. F-statistics (Fis, Fst, Fit) means, standard errors (SE), and 95% confidence intervals obtained for the Guzerá dairy metapopulation, with the subpopulations defined as farm sampled (on the left) or as the lineage (on the right).
F-Statistics
Farms as SubpopulationsLineages as Subpopulations
LimitsFisFstFitFisFstFit
Value (SE)-0.029
(0.022)
0.027
(0.002)
0.055
(0.021)
0.022
(0.022)
0.034
(0.002)
0.055
(0.021)
95% confidence intervalMin:−0.0070.0230.021−0.0140.0290.021
Max:0.0760.0310.1000.0700.0380.100
SE: Standard error. Min: minimum and Max: Maximum.
Table 3. Lineage Fis, farm Fis, pair-wise lineages Fst (above diagonal), and pair-wise farms Fst (below diagonal). Average values of Fst and D’s for each farm are presented below the matrix (italic), while average values of Fst and D’s for each lineage (LINE) are presented above the matrix.
Table 3. Lineage Fis, farm Fis, pair-wise lineages Fst (above diagonal), and pair-wise farms Fst (below diagonal). Average values of Fst and D’s for each farm are presented below the matrix (italic), while average values of Fst and D’s for each lineage (LINE) are presented above the matrix.
LineageLINE1LINE2LINE3LINE4LINE5LINE6LINE7LINE8LINE9LINE10LINE11LINE12LINE13LINE14LINE15
FARMS
FARM10.018
0.003
0.0190.0310.0170.0100.0130.0580.0100.0080.0310.0100.0350.0130.0240.008LlNE1
FARM20.0130.008
0.033
0.0190.0360.0200.0170.0580.0240.0220.0270.0250.0240.0190.0220.019LINE2
FARM30.0400.0250.061
0.020
0.0350.0110.0190.0600.0140.0180.0390.0200.0300.0170.0280.009LINE3
FARM40.0320.0150.0250.003
0.022
0.0240.0140.0660.0180.0220.0390.0090.0320.0300.0230.027LINE4
FARM50.0320.0160.0420.0230.027
0.031
0.0130.0250.0130.0130.0270.0120.0130.0140.0110.013LINE5
FARM60.0360.0220.0480.0270.0230.037
0.048
0.0520.0110.0120.0250.0090.0200.0130.0180.014LINE6
FARM70.0250.0130.0310.0170.0140.0130.038
0.006
0.0300.0410.0870.0440.0620.0470.0630.029LINE7
FARM80.0250.0130.0280.0190.0140.0090.011−0.030
0.033
0.0140.0340.0070.0160.0180.0130.014LINE8
FARM90.0170.0110.0310.0170.0350.0160.0080.010−0.019
0.018
0.0400.0110.0210.0120.0130.009LINE9
FARM100.0290.0140.0300.0170.0190.0210.0080.0140.0120.020
0.046
0.0340.0420.0290.0360.034LINE10
FARM110.0170.0090.0370.0160.0290.0190.0080.0130.0160.0130.028
−0.068
0.0190.0190.0140.016LINE11
FARM120.0240.0140.0400.0200.0230.0140.0080.0110.0090.0110.0090.031
0.024
0.0200.0390.014LINE12
FARM130.0100.0060.0240.0160.0100.0120.0110.0090.0050.0130.0060.0090.032
0.017
0.0160.011LINE13
FARM140.0030.0140.0390.0180.0260.0140.0060.0090.0230.0140.0180.0150.008−0.109
0.077
0.011LINE14
FARM150.00290.0180.0430.0240.0350.0160.0100.0130.0210.0170.0200.0150.0090.0170.005
−0.043
LINE15
FARM1FARM2FARM3FARM4FARM5FARM6FARM7FARM8FARM9FARM10FARM11FARM12FARM13FARM14FARM15
Average Fst0.0210.0250.0250.0280.0160.0180.0520.0170.0180.0370.0180.0280.0200.0240.016Lineage
0.0240.0150.0350.0200.0240.0210.0130.0140.0170.0170.0170.0160.0110.0160.021Farm
Average D’s0.200.190.210.200.120.130.330.140.140.250.140.220.140.170.13Lineage
0.1710.1080.2410.1510.1900.1440.1030.1080.1710.1180.1260.1100.1120.1340.142Farm
Table 4. Results of the evaluation of genetic diversity in other breeds.
Table 4. Results of the evaluation of genetic diversity in other breeds.
BreedNumber of IndividualsNumber of MarkersMean Expected Heterozygosity (He)Mean Observed Heterozygosity (Ho)F-StatisticsReference
Ethiopian Cattle351300.7260.674Fst = 0.013[16]
Creole Cattle857190.7380.718Fis = 0.028[45]
Korean Cattle
Chinese Cattle
Japanese Black Cattle
European Holstein
200130.728
0.744
0.471
0.693
0.721
0.745
0.515
0.753
Fst = 0.109[46]
Sahiwal
Hariana
Deoni
136200.61
0.66
0.70
0.42
0.53
0.59
Fst = 0.067[47]
Japanese Black Cattle252200.6180.623Fst = 0.151[48]
Jersey223120.6430.6355Fst = 0.013
Fst = 0.035 ††
[49]
Burlina Cattle279240.690.63Fst = 0.036[50]
Pirenaica
Betizu
Terreña
Monchina
302110.688
0.715
0.747
0.762
0.682
0.651
0.737
0.757
Fst = 0.041[51]
Caracu; Criolo Lageano; Curraleiro; Mocho Nacional and Pantaneiro; Holstein and Jersey623220.7930.695Fst = 0.061[15]
Nellore; Gyr and Guzerá292220.7480. 645Fst = 0.040[15]
Vietnamese indigenous cattle410270.7600.680Fst = 0.04[52]
Between parishes; †† Between Farms.
Table 5. Genetic diversity (GD), within-subpopulation contribution to GD (Internal_Diversity), between-subpopulations contribution to GD (Mean_Distance), and total contribution to GD (Loss/Gain) after hypothetically removing each subpopulation.
Table 5. Genetic diversity (GD), within-subpopulation contribution to GD (Internal_Diversity), between-subpopulations contribution to GD (Mean_Distance), and total contribution to GD (Loss/Gain) after hypothetically removing each subpopulation.
LineageGDInternal DiversityMean DistanceLoss/Gain
LINE10.77490.0046−0.0358−0.0312
LINE20.7731−0.1360−0.1267−0.2627
LINE30.77500.0484−0.0671−0.0186
LINE40.77590.2753−0.17400.1013
LINE50.7764−0.08830.24670.1585
LINE60.7745−0.21570.1370−0.0787
LINE70.77590.5530*−0.4474*0.1056
LINE80.7751−0.19290.1903−0.0026
LINE90.7747−0.17120.1201−0.0511
LINE100.77390.2400−0.3954−0.1554
LINE110.7746−0.23970.1690−0.0707
LINE120.77460.0260−0.0925−0.0665
LINE130.77590.03620.06090.0971
LINE140.77590.1356−0.03590.0997
LINE150.7739−0.34570.1868−0.1589
* The higher values found for Internal Diversity and Mean distance are indicated in bold.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Peixoto, M.G.C.D.; Carvalho, M.R.S.; Egito, A.A.; Steinberg, R.S.; Bruneli, F.Â.T.; Machado, M.A.; Santos, F.C.; Rosse, I.C.; Fonseca, P.A.S. Genetic Diversity and Population Genetic Structure of a Guzerá (Bos indicus) Meta-Population. Animals 2021, 11, 1125. https://doi.org/10.3390/ani11041125

AMA Style

Peixoto MGCD, Carvalho MRS, Egito AA, Steinberg RS, Bruneli FÂT, Machado MA, Santos FC, Rosse IC, Fonseca PAS. Genetic Diversity and Population Genetic Structure of a Guzerá (Bos indicus) Meta-Population. Animals. 2021; 11(4):1125. https://doi.org/10.3390/ani11041125

Chicago/Turabian Style

Peixoto, Maria Gabriela C. D., Maria Raquel S. Carvalho, Andrea A. Egito, Raphael S. Steinberg, Frank Ângelo T. Bruneli, Marco Antônio Machado, Fernanda C. Santos, Izinara C. Rosse, and Pablo Augusto S. Fonseca. 2021. "Genetic Diversity and Population Genetic Structure of a Guzerá (Bos indicus) Meta-Population" Animals 11, no. 4: 1125. https://doi.org/10.3390/ani11041125

APA Style

Peixoto, M. G. C. D., Carvalho, M. R. S., Egito, A. A., Steinberg, R. S., Bruneli, F. Â. T., Machado, M. A., Santos, F. C., Rosse, I. C., & Fonseca, P. A. S. (2021). Genetic Diversity and Population Genetic Structure of a Guzerá (Bos indicus) Meta-Population. Animals, 11(4), 1125. https://doi.org/10.3390/ani11041125

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop