Position of Hungarian Merino among other Merinos, within-breed genetic similarity network and markers associated with daily weight gain
Article information
Abstract
Objective
In this study, we aimed to position the Hungarian Merino among other Merino-derived sheep breeds, explore the characteristics of our sampled animals' genetic similarity network within the breed, and highlight single nucleotide polymorphisms (SNPs) associated with daily weight-gain.
Methods
Hungarian Merino (n = 138) was genotyped on Ovine SNP50 Bead Chip (Illumina, San Diego, CA, USA) and positioned among 30 Merino and Merino-derived breeds (n = 555). Population characteristics were obtained via PLINK, SVS, Admixture, and Treemix software, within-breed network was analysed with python networkx 2.3 library. Daily weight gain of Hungarian Merino was standardised to 60 days and was collected from the database of the Association of Hungarian Sheep and Goat Breeders. For the identification of loci associated with daily weight gain, a multi-locus mixed-model was used.
Results
Supporting the breed's written history, the closest breeds to Hungarian Merino were Estremadura and Rambouillet (pairwise FST values are 0.035 and 0.036, respectively). Among Hungarian Merino, a highly centralised connectedness has been revealed by network analysis of pairwise values of identity-by-state, where the animal in the central node had a betweenness centrality value equal to 0.936. Probing of daily weight gain against the SNP data of Hungarian Merinos revealed five associated loci. Two of them, OAR8_17854216.1 and s42441.1 on chromosome 8 and 9 (−log10P>22, false discovery rate<5.5e–20) and one locus on chromosome 20, s28948.1 (−log10P = 13.46, false discovery rate = 4.1e–11), were close to the markers reported in other breeds concerning daily weight gain, six-month weight, and post-weaning gain.
Conclusion
The position of Hungarian Merino among other Merino breeds has been determined. We have described the similarity network of the individuals to be applied in breeding practices and highlighted several markers useful for elevating the daily weight gain of Hungarian Merino.
INTRODUCTION
Merino is a globally distributed fine-wool sheep breed representing about a quarter of the global sheep population [1]. The appearance of the first European fine-wool sheep in the southern region of Spain and Italy dates back to the 1st century BCE. These animals were cross-bred with local and imported Arabian sheep [2], likely resulting in a breed that can be considered an early ancestor of the Merino. Ciani et al [3] lean towards the opinion that the Merino breed has an Iberian origin with other stock contributions from the Mediterranean region. Until 1760, the export sale of Merinos from Spain was prohibited. Thenceforward, small flocks were exported to several European countries [4]. The first flock arrived in Hungary in 1765 due to a decree emitted by Empress Maria Theresa of Austria [5]. This importation was followed by two others in 1773 and 1775 [6]. The last Merino export from Spain to Hungary was recorded in 1802.
The Industrial Revolution in Britain led to considerable improvements in weaving efficiency. The development of the wool textile industry resulted in an increasing global demand for wool. The selection process for wool quantity and quality resulted in many types of Merino sheep worldwide [7]. The Hungarian Combed Merino breed was developed from imported Merinos, the local Racka breed, and German domestic sheep [8]. Importation of Rambouillet sheep from France (1816) and North-German Combed sheep has also contributed to the development of the Hungarian breed [6]. The first Hungarian herd book was introduced in 1859, in this period, Merino and Merino-derived sheep reached 31.5% of the population among Hungarian sheep [4].
At the end of the 19th century, in spite of the Western European meat-breeding trend, wool was still the main production objective in Hungary [9]. After World Wars I and II, crossbreeding resulted in a fine and short-wool population, which became more homogeneous by breeding Précoce and Rambouillet rams [6,9]. In the 1950s, the breed was reshaped by fine-wool variants from the former Soviet Union (Askanian, Caucasian, Stavropol, and Grozny Merinos). This decision led to considerable increases in the body mass and the wool production capabilities of the breed [10]. In the 1960s, French Précoce and German Mutton Merinos were imported to Hungary to improve the breed's meat-production potential. In the 1970s through 1980s, several attempts were made to improve the production parameters and prolificacy of Merino flocks. For that very reason, Kent, Corriedale, and Australian Booroola Merinos were introduced into Hungary [11]. These crosses also breed better muscularity and wool quality [1,10].
In the 1980s, wool production was pushed into the background, and Hungarian Merino was considered a dual-purpose (meat and milk producing) breed [12]. A drastic change in agricultural subsidies followed the democratic transformation of the former socialist regime in 1989. During this period, exported lamb became the main product of sheep husbandry, and the average daily weight gain of lambs was about 250 g/d [13].
In the 2010s, due to considerable declines in breeding activity, Hungarian Merinos became classified as endangered. Subsequently, the Hungarian government began to provide subsidies for the conservation of genetic resources as part of the coordinated breeding program. In 2020, the Hungarian Merino population amounted to 8,167 ewes. At present, only 11 Hungarian genealogical lines are noted [1] out of the 39 recorded in 2004.
We had three aims in this study: i) Since the genetic position of the Hungarian Merino is not known based on single nucleotide polymorphism (SNP) markers, our samples were compared to the database of Ciani et al [3] to determine the genetic position of Hungarian Merinos among other Merino-derived sheep breeds. ii) However, there are publications considering the gene network [14] or genetic connections among the breeds [15], we wanted to explore the characteristics of our sampled animals' genetic similarity network within the breed to increase the ability of sheep breeders (including those keeping other sheep or non-sheep populations in the world) to easily adopt decisions on an individual level. iii) Since nowadays the most important parameter of the Hungarian sheep-breeding industry is the weight-gain of lambs, we decided to deliver markers suitable to improving the daily weight-gain of Hungarian Merino lambs. Our investigation overlaps other activities like investigating growth and meat production traits in several other sheep breeds [16–19] and highlights regions in the sheep genome that are usable not only for Hungarian Merino but for other breeds as well.
MATERIALS AND METHODS
Approval from the animal research ethics committee was not necessary as the ear tag application is part of the routine breeding procedure.
Animals, genotyping, and quality control
The sampling locations covered Hungary, which has a quite uniform, continental climate. One hundred thirty-eight Hungarian Merino ewes were included in this study. They were born in different years spanning 10 years. Their daily weight gain was standardised to 60 days as follows: 60 BM30–80/d, where BM30–80 is the body mass of the lamb measured between 30 and 80 days of age, d is the actual number of days after birth on the day of the measurement.
Ear cartilage samples were collected by the TypiFix (Agrobiogen GmBH, Hilgertshausen, Germany) tissue sampling system [20] during routine ear tag application. For comparative studies, we used the database provided by Ciani et al [3]. SNP typing of Hungarian Merinos was performed on Ovine SNP50 Bead Chip (Illumina, San Diego, USA) by Neogen Corporation (Ayr, UK).
Quality control of data included filtering for mapped and polymorphic markers (minor allele frequency>0.05) with a call rate higher than 0.95. Duplicated samples were not detected (identical by state value, identical by state (IBS)>0.95). Call rates of the genotypes were over 0.95 for all Hungarian Merino samples. The matched database, which included both the Sheep Consortium [3] and Hungarian data, contained a total of 31 breeds, 693 animals (Table 1), and 22,265 SNPs.
In case of within-breed network and association analyses 46,906 SNPs were kept after data filtering.
Population genetics analyses
Determination of observed and expected heterozygosity (Ho, He), subsequent calculation of inbreeding coefficient (FIS), genetic distance (FST), and principal component analysis (PCA) were performed by SNP & Variation Suite 8.1.1 program (GoldenHelix, Bozeman, MT, USA). The inbreeding coefficient is the (observed - expected number of homozygous markers)/(number of genotyped markers - expected number of homozygous markers). Expected number of homozygous markers
Admixture software v.1.3 has been used with --cv option for determining the most probable cluster number (K) from the value of cross-validation error in each Ks [23]. Before analysis, genotypes were transformed by PLINK using --recode12 switch.
Treemix 1.13 [24] software has been applied to calculate a dendrogram; block size was 500, the number of repetitions was 1,000. Bootstrap values were determined by the Phylip 3.697 program package [25].
Network analysis concerning the relation of the Hungarian Merinos was based on the identity-by-state value of pairs calculated by SVS program (GoldenHelix, USA).
Identity-by-state pairwise distance between any two individuals is ([number of markers sharing two alleles+0.5× number of markers sharing one allele]/number of markers).
Within-breed network and association analyses
Based on IBS values, the net properties of the Hungarian Merino samples, such as diameter and betweenness centrality, were calculated and visualised by Python 3.6 software using the libraries networkx 2.3 and matplotlib 3.1.1. Diameter is the largest distance observed between the pairs of nodes or the distance between the two furthest nodes [26]. Betweenness centrality of a given animal/node
For the identification of loci associated with daily weight gain, a multi-locus mixed-model was used [27]. The model was y = Xβ+Zu+e, where y is a continuous variable, the daily weight gain, X is the matrix of fixed effects composed of 46,906 SNPs and covariates (date and place of birth), Z is the incidence matrix relating random effects to each individual, which becomes an identity matrix, because one phenotypic measurement is performed per individual [28]. e stands for the residual effects, and β and u are vectors representing coefficients of fixed and random effects. For further details on mixed linear models, see the link [29].
RESULTS
Diversity and structure of populations
Observed heterozygosity ranged from 0.320 to 0.392, slightly below the expected value of 0.397. Based on the observed and expected heterozygosity values, the inbreeding coefficient of breeds ranged from 0.013 to 0.193 (Table 1). Cordoba Spanish Merino was an outlier with its highest FIS value (0.341), whereas, in Hungarian Merino, this value was 0.054. Pairwise values of FST calculation fell in the range of 0.004 to 0.313 (Supplementary Table S1). FST values of Hungarian Merino relative to the other breeds were between 0.035 and 0.184 (Table 2).
Principal component analysis based on pairwise identity-by-state of the animals (Figure 1A–D) placed the Hungarian Merino population on the outer rim of the investigated populations, close to Rambouillet and Chinese Merino.
The most probable cluster number of the 31 breeds produced by the admixture analysis was nineteen (K = 19, Figure 2). When K was equal to the number of analysed breeds, most of the Hungarian Merinos displayed admixture with Rambouillet (Supplementary Figure S1). On the dendrogram of the breeds, Hungarian Merino was placed on a common branch with Rambouillet, Chinese Merino, Merinizzata, and Estremadura Spanish Merino (Figure 3; Supplementary Figure S2).
Identity-by-state network and markers associated with daily weight gain
For visualisation of the genetic net of the individuals, IBS values (edges) above 0.73 have been maintained. In the analysed population, this value provided an interconnected network of all nodes (individuals). The number of edges was 341, the average degree was 5.014, net diameter was 4 (Figure 4A). The five highest betweenness centrality values were 0.936, 0.035, 0.034, 0.029, and 0.023. When the animal having the highest betweenness centrality (0.936) was removed, the interconnected network contained 108 animals and 199 edges (Figure 4B). In this case, average degree became 3.682, and the diameter of the net expanded to 10. The highest five betweenness centrality values were 0.744, 0.261, 0.235, 0.185, and 0.141.
The values of daily weight gain ranged from 261.6 to 494.7 grams. In multi-locus mixed-model analysis, several loci were identified to be associated with the daily weight gain (Table 3). The top hits were five in number (−log10P>10) on chromosome 2, 3, 8, 9, and 20 with false discovery rates lower than 4.1e–11.
DISCUSSION
FIS value calculated from observed and expected heterozygosity of Hungarian Merinos is below 0.055. Ten out of the thirty-one studied breeds have small FIS values, indicating that there are not considerably more inbred animals than expected.
We have found that the calculated genetic distances and the position of Hungarian Merinos in PCA and on the phylogenetic tree support the contribution of Rambouillet imported twice [6], to our breed. Similar positions of the Chinese Merino and Rambouillet on the PCA plot and the dendrogram indicate a close relationship between the two breeds. As for the Spanish breeds presented in the database, Estremadura Merino is the closest to Hungarian Merino, presenting its strongest influence on the breed (Table 2). Pairwise FST value of the Spanish Cordoba and Hungarian Merino breed was the highest (FST = 0.184), ruling out the Cordoba breed as an origin breed of Hungarian Merino. The phylogenetic study placed Estremadura Merino the fourth closest to Hungarian Merino, while the Rambouillet breed remained the closest (Supplementary Figure S2).
In Admixture analysis, the most probable cluster number (K) was 19 (Figure 2; Supplementary Figure S1); however, the Hungarian Merino’s admixture with Rambouillet was seen at K = 29 and 31 (Supplementary Figure S1).
Network analysis of the individuals revealed (Figure 4A) a spoked wheel structure [26], which is probably induced by the extensive usage of preferable individuals. When the animal with the highest betweenness centrality in the genetic similarity net was removed, the wheel structure disappeared, and more details were revealed. In Figure 4B, we can identify other animals which have strong influences on the breed composition relative to the others. Based on network data, breeders can select suitable animals for optimal mating schemes. Visualisation of edge strengths (a.k.a. pairwise IBS values between the animals) promotes easy recognition of closely related animals (Figure 4). In this way, deterioration of heterozygosity parameters, inbreeding depression, and declination of economic performance can be avoided. Improvement of the visualisation of genetic similarities can further enhance the applicability of this analysis. Improvements could be achieved by pop-up windows carrying cross-references to the breeders, farms, age, ancestors, or by time-lapse interpretation of the connections as members of the populations are changing via newborn and culled animals.
Concerning daily weight gain, we have found five highly correlated loci (Table 3) using 138 animals. However, the sample size is lower than that used by Zhang et al [16] or Lu et al [18], comparing our highly correlated loci with markers revealed by other teams on similar data hints that the main players in a multigenic trait could be exposed with 100+ animals as well. Two of our five markers (OAR8_17854216.1 on chromosome 8 and s42441.1 on chromosome 9) were in close vicinity, −1.3 and +4.8 million bases, respectively, to the top findings of Zhang et al [16]. The reported markers were associated with daily weight gain, six-month weight, and post-weaning weight. A third marker (s28948.1) we highlighted was 11.7, and 1.5 million bases further away from a post-weaning gain [16] and a weaning weight [18] associated marker, respectively. The close vicinity of the above mentioned three markers to those identified by other teams [16,18] might be explained by the high similarity of Hungarian and Chinese Merino revealed by principal components (Figure 1), pairwise genetic distance (Table 2; Supplementary Table S1), and dendrogram (Figure 3; Supplementary Figure S2). Our markers on chromosome 2 and 3, OAR2_137539806.1 and OAR3_188812728.1, or their neighbouring markers closer than twelve million base pairs, were not reported elsewhere related to weight-gain properties.
This study gives insight into the overall genetic state of the Hungarian Merino in contrast to other breeds. The results might be useful for other research conducted on a wider set of Merino breeds.
The within-breed network analysis of the animals based on their genetic similarity and its visualisation could provide an intuitive tool to the breeders of any species to easily have a glimpse of the genetic structure of the population and the most important animals. Later on, when more within-breed structures could be analysed and compared, we could dissect different breeding practices among different breeders or countries, giving further insight into the influence of the breeding policy on the structure of a given breed.
Some of the loci associated with a standardized daily weight gain were very close to the markers presented in other studies on 2.4 and 3.4 higher number of animals (Table 3). Application of these markers could improve the daily weight-gain of lambs, which is important as lambs production provides a considerable income for the breeders.
ACKNOWLEDGMENTS
This research has received a financial grant entitled TJUGEN from the Hungarian Ministry of Agriculture. The authors express their gratitude to the Association of Hungarian Sheep and Goat Breeders for providing the samples and for the supportive conversations regarding our research. Many thanks to the research of Ciani et al [3] for providing the basis of the comparison of Hungarian Merino to other sheep breeds. Thanks to Anneliese Kleinschmidt for proofreading the manuscript. Thanks to our reviewers for providing their valuable notices, comments, and questions, and many thanks to the editorial staff as well.
Notes
CONFLICT OF INTEREST
We certify that there is no conflict of interest with any financial organization regarding the material discussed in the manuscript.
FUNDING
The authors received no financial support for this article.
SUPPLEMENTARY MATERIAL
Supplementary file is available from: https://doi.org/10.5713/ab.21.0459
Supplementary Table S1. Pairwise values of FST and their heat map
ab-21-0459-suppl1.pdfSupplementary Figure S1. Admixture analysis of 31 Merino-derived breeds and for a range of K-values (K = 2 to 31) and cross-validation error at different K values. Each individual is represented by a single column divided into K coloured segments, where K is the number of assumed clusters.
ab-21-0459-suppl2.pdfSupplementary Figure S2. UPGMA tree of Nei’s genetic distance sheep breeds. Bootstrap values are indicated on the nodes.
ab-21-0459-suppl3.pdf