- Research Article
- Open access
- Published:
Genetic legacy and adaptive signatures: investigating the history, diversity, and selection signatures in Rendena cattle resilient to eighteenth century rinderpest epidemics
Genetics Selection Evolution volume 56, Article number: 32 (2024)
Abstract
Background
Rendena is a dual-purpose cattle breed, which is primarily found in the Italian Alps and the eastern areas of the Po valley, and recognized for its longevity, fertility, disease resistance and adaptability to steep Alpine pastures. It is categorized as 'vulnerable to extinction' with only 6057 registered animals in 2022, yet no comprehensive analyses of its molecular diversity have been performed to date. The aim of this study was to analyse the origin, genetic diversity, and genomic signatures of selection in Rendena cattle using data from samples collected in 2000 and 2018, and shed light on the breed's evolution and conservation needs.
Results
Genetic analysis revealed that the Rendena breed shares genetic components with various Alpine and Po valley breeds, with a marked genetic proximity to the Original Braunvieh breed, reflecting historical restocking efforts across the region. The breed shows signatures of selection related to both milk and meat production, environmental adaptation and immune response, the latter being possibly the result of multiple rinderpest epidemics that swept across the Alps in the eighteenth century. An analysis of the Rendena cattle population spanning 18 years showed an increase in the mean level of inbreeding over time, which is confirmed by the mean number of runs of homozygosity per individual, which was larger in the 2018 sample.
Conclusions
The Rendena breed, while sharing a common origin with Brown Swiss, has developed distinct traits that enable it to thrive in the Alpine environment and make it highly valued by local farmers. Preserving these adaptive features is essential, not only for maintaining genetic diversity and enhancing the ability of this traditional animal husbandry to adapt to changing environments, but also for guaranteeing the resilience and sustainability of both this livestock system and the livelihoods within the Rendena valley.
Background
Rendena is an indigenous, dual-purpose cattle breed, which is valued by farmers for both milk and meat quality, and reared principally in the Italian Alps and the eastern areas of the Po valley. The breed is named after the Rendena valley, a UNESCO World Heritage Site nestled between the Brenta Dolomites and the Eastern Alps in the Province of Trento, Italy. This breed is relatively small, with short robust legs, small hooves and good grazing ability, and is well adapted to grazing on the steep Alpine summer pastures. Rendena cattle are also popular locally for their longevity, high fertility and disease resistance [1]. However with only 6057 animals in the national Herd Book (http://www.anare.it/) in 2022 among which 4170 are more than two-year old females, the FAO Domestic Animal Diversity Information System (https://www.fao.org/dad-is/en/) has categorized the Rendena breed as ‘vulnerable to extinction’.
Successive colonization by diverse human populations across the centuries has led to substantial human cultural diversity in the Trentino-Alto Adige/Südtirol region [2]. The Rendena valley is no exception, with a rich history of human settlements dating back to prehistoric times, according to archaeological studies [3, 4]. This historical trajectory has given rise to roughly 30 local cattle breeds across the Alpine arc [5], with the first written evidence of Rendena cattle dating to the beginning of the eighteenth century [6]. During that same period, successive epidemics of rinderpest spread throughout the Alpine cattle populations of northern Italy, including Val Rendena. To restock their herds, farmers imported cattle of Swiss origin which were phenotypically similar to the autochthonous population [6]. When the epidemics ended in the late eighteenth century, the importation of cattle ceased, and crossbreeding with neighbouring breeds was mainly limited by geographical barriers. During the nineteenth century, the Rendena breed expanded its distribution in northern Italy, especially to Veneto and Lombardy, establishing itself for a time as the most wide-spread dairy breed in the area [6]. This expansion was interrupted between the two World Wars (c. 1919–1938) when the Italian government imposed the replacement of local breeds with more productive ones. As a consequence, the number of Rendena cattle declined to a few thousand individuals [6]. In spite of such politically-motivated disincentives, the breed survived thanks to the dedication of local farmers in the Rendena valley, who continued to rear their native cattle. More recently, the Rendena breed was fully recognised and protected with the establishment of the National Association of Breeders of Rendena Cattle (A.N.A.RE) in 1981.
To date, there has been no comprehensive analyses of the molecular diversity of the Rendena breed and information on its genetic origin, molecular diversity and conservation status are still lacking. Here, we analyse a genomic single nucleotide polymorphism (SNP) dataset that includes 28 Rendena individuals sampled in 2000 and 140 individuals sampled in 2018, as well as those from other local and highly selected breeds to (i) assess changes in population genetic diversity and structure of this breed over the last two decades; (ii) establish the origin of the breed, and (iii) identify breed-specific genomic signatures of selection.
Methods
Sample collection, genotyping, and quality control
Fresh blood samples of 140 Rendena cattle were collected from 27 farms located in the Province of Trento, and six farms located in the Province of Padua in the region of Veneto (Italy). To maximise the representativeness of the dataset, samples were selected to minimise the number of paternal half sibs and avoid sampling within the same maternal lineage. The birth years of the sampled animals ranged from 2010 to 2017. Rendena heifers typically calve for the first time at around three years of age. These animals are maintained in the herd for four or five lactations; therefore, the average lifespan of the animals is around eight years. Blood samples were collected following the Italian and European legislation on animal welfare (D.lgs n. 146/2001, Council Directive 98/58/CE) by a licensed veterinarian, according to the European directive 2010/63 during annual screening campaigns. In addition, 28 Rendena samples from the year 2000, collected during a previous European research project (EU RESGEN CT 98–118), were included in our analysis.
Data from 17 local breeds from the Alpine region, including 24 Rendena samples collected in the region of Veneto, were available from an independent investigation on Alpine cattle biodiversity [5]. Moreover, publicly available genotype data of 11 Central European [5, 7, 8] and two unpublished East European cattle breeds were included in the dataset (see Additional file 1: Table S1).
DNA was extracted at the laboratory of Università Cattolica del Sacro Cuore (Italy) from the 168 Rendena samples using the GenElute Mammalian Genomic DNA Miniprep kit (Sigma, St Luis, MO, USA) according to the manufacturer’s instructions. DNA samples were genotyped with the GeneSeek (Neogene, Lincoln, Nebraska) GGP Bovine 100 K Illumina SNP chip (Illumina Inc.). For this study, two datasets were built: the first comprising the data from the 168 new Rendena samples genotyped with the GGP Bovine 100 K Illumina SNP chip, and the second also including data from the 30 Alpine, Central and East European breeds, after updating marker positions of all the samples to the ARS-UCD1.2 bovine reference genome [9], and retaining only consensus markers between different SNP chips. In the merged dataset, Rendena cattle from this investigation were labelled RENgen (after the Rendenagen project), and the subsets of animals sampled in years 2000 and 2018 named Rendena2000 and Rendena2018, respectively. The publicly available data from 24 extra Rendena samples were labelled REN.
Both datasets were quality-controlled using the PLINK v1.9 software [10] to achieve a SNP call rate > 0.98, an individual call rate > 0.98, and a minor allele frequency (MAF) > 0.01. Markers located on sex chromosomes or with an unknown map position were removed. Pruning for linkage disequilibrium (LD) was performed with the PLINK flag –indep-pairwise setting an r2 value of 0.2, a window of 50 and a step size of 5.
Characterization of the Rendena breed within local and European genetic contexts
The merged dataset was used to investigate potential gene flow and admixture that may have occurred over centuries between the Rendena and other breeds. We calculated observed (HO) and expected (HE) heterozygosity for each breed included in this study using PLINK; instead, inbreeding coefficient (FIS) values were calculated using the Arlequin v3.5.2.2 software, with the threshold for statistical significance set at p < 0.05 [11]. To estimate the effective population size (Ne), two different approaches were used: SneP v1.1 [12] with default parameters (to facilitate comparison with results from previous publications on local cattle breeds) and currentNe v1.0 [13]. Runs of homozygosity (ROH)-derived inbreeding coefficients (FROH) were calculated independently for each population, as implemented in the R package detectRUNS v 0.9.6 [14]. ROH values were calculated with PLINK using the following parameters: (i) sliding window of 50 SNPs; (ii) a maximum of one heterozygous genotype and one missing genotype allowed per ROH; (iii) minimum number of SNPs in a ROH calculated according to the formula described by Purfield and colleagues [15]; (iv) a minimum ROH length of 500 kb; (v) a minimum density of one SNP per 100 kb and; (vii) a maximum gap of 1 Mb between consecutive homozygous SNPs.
Changes in population structure were assessed with a principal component analysis (PCA) that was performed using PLINK and plotted with the R software [16]. The software Admixture v1.3.0 implements a maximum-likelihood based approach to infer ancestry proportions by evaluating a K number of theoretical ancestral populations [17]. Here, we tested ancestry solutions for K ranging from 2 to 30, and the built-in cross-validation (CV) error approach was used to identify the best fitting value of K. To alleviate sample-size bias in the analysis, we used the BITE v1.2 R package [18] to subset the most numerous populations to a maximum of 30 individuals, selected to mimic the population structure of the full set (the subset of 30 Rendena samples includes animals from previously published data and from this study).
Relationships among populations were explored by pairwise Reynolds genetic distances (computed with a custom script) and visualized with a Neighbour-net graph produced with the SplitsTree v4.14.6 software [19]. To assess population splits and gene flow events, the Treemix software [20] was run on the whole dataset and on a reduced dataset comprising the following breeds: Rendena, Brown Swiss (BSW), Original Braunvieh (OBV) and Murnau-Werdenfelser (MWF), using Jersey (JER) as the root. We set the block size for jackknifing (-k) to 500 SNPs and tested a number of migration events (m) from 0 to 10 in the whole dataset and from 0 to 5 in the reduced set.
Evaluation of the recent genetic management of the Rendena breed
We used the SNP chip data from 168 Rendena individuals to evaluate how selection has affected the genetic structure and level of inbreeding of the Rendena cattle population over the past 18 years, by comparing data from samples collected in the year 2000 with data from samples collected for this study in the year 2018. ROH were calculated independently for each of the two groups with PLINK according to the parameters given above, and used to evaluate the inbreeding level in Rendena2000 and Rendena2018, using individual and per-group FROH values calculated with the R package detectRUNS. To enhance the reliability of the comparison between populations with different sample sizes (N = 28 in Rendena2000 and N = 140 in Rendena2018), for Rendena2018, ROH mean statistics were computed following 1000 bootstrap iterations of 28 randomly selected animals. Confidence intervals (95%) were computed on the results obtained from the two Rendena populations. Changes in population structure were assessed with PCA performed with PLINK and visualised in R. To further test the genetic divergence between the two populations, a discriminant analysis of principal components (DAPC) was performed and visualised in R using the package Adegenet [21].
Analyses of signatures of selection
ROH, cross-population extended haplotype homozygosity (XP-EHH) and integrated haplotype score (iHS) were used to scan for signatures of selection in the Rendena breed, which may have occurred at different times during the history of the breed.
ROH
ROH values were investigated in the RENgen sample set to detect human and environmental-mediated selection. To identify regions putatively under selection, we used the R package detectRUNS to calculate the frequency with which each SNP falls within a ROH [14]. SNPs with values in the top 1% of the distribution were identified as potential locations under selection [22]. The R package GALLO [23] was used to annotate each significant SNP, considering a flanking region of 17.8 kb, upstream and downstream of the target SNP, corresponding to half the mean distance between adjacent markers.
EHH-based metrics
To evaluate genome-wide selection, we used EHH and estimated the decay of haplotype homozygosity as genetic distance increases. Two EHH-based metrics were calculated, iHS within population and XP-EHH between pairs of populations using the rehh v2.0 package [24]. The first analysis was performed on the RENgen samples only, while the XP-EHH analysis was also performed on the Original Braunvieh and modern Brown Swiss breeds (BSW and BSW_IT), and results were compared to identify shared and divergent signatures of selection.
To polarise the iHS analysis, we used the ancestral and derived allele cattle data from Naji and colleagues [25]. However, as information on the ancestral alleles was not available for all markers in the dataset, we performed the analysis in two stages. First, the subset of markers with available ancestral allele information was analysed, then the analysis was repeated using the entire dataset, assuming the major allele is the ancestral allele. The significance threshold was set at a p-value = 0.01.
The XP-EHH metric identifies selective sweeps where an allele undergoes strong directional selection in one population while remaining polymorphic in the overall population [26]. We applied XP-EHH to compare the Rendena versus both the Original Braunvieh and the Modern Brown Swiss breeds. Markers with an XP-EHH score in the top 1% of the standardized distribution were identified as under putative positive selection. Only signatures of selection relative to the Rendena breed are discussed below. Markers with a value above the threshold were annotated with the GALLO R package considering a flanking region of 17.8 kb, upstream and downstream the target SNP as described above.
Results
Sample collection, genotyping and quality control
Following quality control of the two datasets, 75,157 SNPs and 167 out of 168 genotyped individuals were retained in the dataset that consisted exclusively of the Rendena samples from the Rendenagen project (RENgen), and 16,785 SNPs and 908 individuals for the dataset including Rendenagen samples, previously published Rendena samples and the breeds chosen for biodiversity analyses.
Characterization of the Rendena breed within local and European genetic contexts
The values for HO and HE ranged from 0.30 and 0.29, respectively, for the JER breed, to 0.40 and 0.40 for the RENgen samples (Table 1). P-values associated with the inbreeding coefficient FIS estimates were not statistically significant (significance at p < 0.05; see Additional file 2: Table S2). Genomic inbreeding estimated from ROH gave values ranging from 0.02 (Cika) to 0.12 (Jersey), with RENgen having values of 0.06 and 0.09 for the Rendena2000 and Rendena2018 populations, respectively. Contemporary Ne estimated by the SNeP software were 180 and 237 for the RENgen populations sampled in the years 2000 and 2018, respectively (Table 1). Ne values of 180 and 269 were obtained from currentNe for the Rendena2000 and Rendena2018 populations, respectively, with 90% confidence intervals ranging from 154 to 209 for the Rendena2000 population, and from 236 to 284 for the Rendena2018 population (see Additional file 3: Table S3).
In the PCA, the first and second principal components (PC) explained 2.82 and 1.85% of the total variance, respectively. The first PC discriminated the modern Brown Swiss group of breeds (BSW_IT, BSW, BRV) from all the other breeds included in the analysis (Fig. 1). In the first PC, Rendena clusters with the Original Braunvieh and Murnau-Werdenfelser Alpine breeds, but forms a separate group from the modern Brown Swiss cluster. The second PC separates two groups, one including the Alpine breeds, positioned according to a west-to-east geographic gradient, and a second comprising Podolian breeds. Burlina falls within this latter group, because of its proximity with Holstein, located at the upper extreme of the second PC (Fig. 1, and see Additional file 4: Fig. S1).
ADMIXTURE analysis at K = 2 separates the modern Brown Swiss breeds from the other breeds, with Original Braunvieh and Rendena sharing more of the Brown Swiss-like ancestral component than the rest of the populations in the dataset (Fig. 2). K = 3 showed a genetic component shared by all the alpine breeds, but it was more markedly present in breeds from the western area of the Alps (FLV, SIM, PRO, MON, ABO).
As illustrated in Fig. 2, at K = 8 an ancestry component (brown colour) was detected, primarily shared by Rendena, Original Braunvieh and Murnau-Werdenfelser, but also present in all the alpine breeds as well as in the Varzese breed from the north-western part of the Apennine region. All K values greater than 8 progressively distinguished individual breeds (Fig. 2, see Additional file 5: Fig. S2). The lowest CV error value corresponded to K = 23 (see Additional file 6: Fig. S3) at which all breeds had distinguishable patterns, with the exception of the modern Brown Swiss group, the Fleckvieh and Simmental, and the REN and RENgen samples.
The Neighbour-net graph based on Reynold’s genetic distances showed that the Rendena breed has an ancestral relationship with the Brown Swiss cluster (Fig. 3). However, as seen in the PCA, the Rendena breed appears to be closer to the root of the branch and to the Original Braunvieh than to the modern Brown Swiss population. The Murnau-Werdenfelser and the Varzese breeds were located nearby, confirming the genetic pattern shown by the Admixture analysis (Fig. 2).
The Treemix analysis performed on the whole dataset for m from 2 to 10 did not identify gene flow events involving the Rendena breed. Consequently, focussed Treemix analysis was performed on a subset of populations including Rendena and the breeds that were identified by the Neighbour-net analysis as genetically closer to Rendena (i.e., BSW, BRV, BSW_IT, OBV, MWF, VAR and PMT), using Jersey as outgroup. The analysis identified gene flow events from Original Braunvieh towards the Rendena, Varzese and the Murnau-Werdenfelser (Fig. 4).
Evaluation of recent genetic management of the Rendena breed
Data from Rendena cattle sampled in the year 2000 (Rendena2000, n = 28) and year 2018 (Rendena2018, n = 140) were analysed independently to assess changes in population structure and molecular diversity over the past ~ 20 years. HO and HE values calculated for the two populations were almost identical, respectively 0.41 ± 0.14 and 0.4 ± 0.11 for the Rendena2000, and 0.40 ± 0.13 and 0.39 ± 0.12 for Rendena2018. In contrast, the mean FROH of Rendena2000 (0.06) was significantly lower than that of Rendena2018 (0.09) (p < 0.001) (Table 2). A mean number of ROH per individual of 27.68 for Rendena2000 and 30.15 for Rendena2018 (p < 0.001; Table 2) was recorded. The average ROH length was 5.21 Mb and 7.01 Mb for Rendena2000 and Rendena2018, respectively (p < 0.001; Table 2). Rendena2000 had a higher frequency of short ROH segments, between 0–2 Mb and 2–4 Mb, while Rendena2018 showed a higher occurrence of medium to long segments, with almost 20% of the ROH segments being longer than 16 Mb (Table 2). PCA identified two partially overlapping clusters (Fig. 5), also confirmed by DAPC analysis (see Additional file 7: Fig. S4).
Analyses of signatures of selection
ROH
Signatures of selection were investigated in RENgen dataset assessing the proportion of times each SNP falls inside an ROH (Fig. 6). This analysis led to the identification of 769 candidate SNPs. These SNPs were subsequently annotated, resulting in a total of 114 genes located on ten chromosomes (BTA3, 4, 5, 6, 10, 12, 13, 16, 22, and 25; see Additional file 8: Table S4). The most represented chromosome was BTA10, on which 48 annotated genes were identified, followed by BTA6 and BTA22 with 19 and 17 genes, respectively.
EHH-based metrics
iHS analysis of the Rendena cattle was initially performed on a dataset of 167 animals and only on the 7626 SNPs for which the ancestral allele information was available. A threshold of significance was set to a nominal p-value of 0.01 genome-wide; however, no SNP was significant. The analysis was then performed on the whole set of SNPs, considering the major allele as ancestral (see Additional file 9: Fig. S5). In this case, two SNPs on BTA7 were identified as being under selection, although no gene was annotated in the 17.8 kb region upstream and downstream of these loci.
XP-EHH analysis was used to compare EHH profiles of the Rendena breed with the Original Braunvieh and the modern Brown Swiss breeds, to detect signatures indicating divergent selection (Fig. 7). In the comparison between Rendena and Original Braunvieh, 22 SNPs with values falling within the top 1% of the standardized distribution (equivalent to an XP-EHH value equal or greater than 2.52) were annotated, including SNPs on chromosomes BTA1, 2, 3, 8, 9, 14, 15, 16, and 20. In addition, in the comparison between Rendena and Modern Brown Swiss, 35 SNPs exceeded the threshold of 2.23, corresponding to the top 1% of the standardized distribution, and were distributed across 13 chromosomes (BTA1, 2, 3, 7, 8, 9, 11, 14, 15, 16, 17, 18, and 23). Among identified markers, 17 SNPs on eight different chromosomes scored a XP-EHH value in the top 1% of the standardized distribution in both comparisons (see Additional file 10: Table S5). In the 17.8 kb region upstream and downstream of the identified SNPs, ten genes were annotated in the comparison between Rendena and Original Braunvieh and 12 in the comparison with Modern Brown Swiss. Seven of the annotated genes were in common between the two comparisons (see Additional file 10: Table S5). There was no shared signature of selection between the two methods.
Discussion
The documented history of Rendena cattle can be traced back to the eighteenth century when restocking efforts were initiated following a series of devastating rinderpest epidemics, which required the introduction of cattle from neighbouring Alpine valleys in different countries [6]. Subsequently, between 1920 and 1940, national legislation mandated the replacement of local cattle with more productive breeds. However, in the Rendena valley, livestock farmers strongly opposed and ultimately resisted this mandate, leading to an exemption from the law in 1937 [6].
Consistent with these historical documents, our admixture analysis at K = 8 (Fig. 2), revealed a genetic component shared among a number of breeds from the Alps and the Po valley. This result is further supported by the Treemix results (Fig. 4), which indicate three gene flow events affecting the northern Italian cattle. The source population was always the Original Braunvieh with three receiving breeds: the Varzese from the western Po valley, the Murnau-Werdenfelser and the Rendena. These genetic signatures likely reflect historical restocking efforts across the Alps, extending down into the Po valley, areas significantly affected by rinderpest [27].
The PCA and Neighbour-net analyses underscored the genetic proximity of the Rendena breed to the Original Braunvieh, a breed with roots tracing back to mediaeval Switzerland and subsequently exported to western Austria and the Italian (southern) Alps in the fourteenth century. In addition, the Rendena breed shows genetic affinity to the highly endangered Murnau-Werdenfelser breed, which originated from crossbreeding between the local Bavarian breeds Oberinntaler and Braunvieh (Figs. 1, 2, and 4) [28]. These findings align with historical records indicating that, in the eighteenth century, the Rendena valley was primarily restocked with Swiss-origin cattle that shared similar morphological traits and suitability to mountainous environments. However, in the twentieth century, the breed also experienced an introgression from the transboundary Brown Swiss population imposed by national legislation, making it challenging to distinguish the genetic signatures of these two events.
Our comparative analysis of Rendena animals sampled at an 18-year interval (2000 vs 2018) revealed a notable mean inbreeding level in 2018. This increase was corroborated by the comparison of the mean number of ROH per individual, which was significantly larger in the 2018 sample (Table 2). In addition, the mean length of ROH segments in the 2018 sample exceeded that recorded in 2000, underscoring a recent escalation in inbreeding. However, this level of inbreeding remains comparable to values recorded for other local breeds from the Alps. Notably, Ne was larger for the 2018 sample, which is potentially due to a sampling strategy tailored to select only female samples from diverse maternal lineages and this may not have been the case for the 2000 samples, although a change in demographic dynamics cannot be excluded. The breeding strategy currently used for the Rendena breed is designed to preserve the dual-purpose nature of the breed; however, there is more emphasis on milk characteristics, with dairy accounting for 65% of the total value of the selection index, and the remaining 35% for beef traits (http://www.anare.it/) [29]. In fact, signatures of selection, revealed by the ROH and XP-EHH analyses that compare the Rendena population with the Original Braunvieh and Brown Swiss breeds, identified distinct signatures of selection in genes related to both meat and milk production (Figs. 6 and 7). Interestingly, a signature of selection on chromosome 6, identified in the ROH analysis, encompasses several genes that are linked to milk production, such as ABCG2 [30,31,32], HERC3 [33, 34], and HERC6 [35, 36]. In addition, it includes genes associated with meat-related traits, such as FAM184B [37,38,39,40] and NCAPG [40,41,42]. The LAP3 gene, which plays a role in hormone level regulation and protein maturation, has previously been associated with body conformation and carcass traits, growth, calving ease and milk production traits in cattle [41, 43,44,45,46,47, 49]. Chromosomes 3 and 5 also displayed signatures of selection in genes previously associated with milk production (RNPC3 and ST8SIA1) [48, 49]. The VPS13D gene, identified through XP-EHH analysis, is involved in cholesterol regulation and lactase persistence in humans [50, 51]. Furthermore, a homozygous segment on chromosome 5 was shared among multiple individuals and contained genes related to feed efficiency (C2CD5) [52] and skeletal muscle development (PYROXD1) [53]. The CADM2 gene on chromosome 1 and the SOX6 gene on chromosome 10, which were identified in a signature of selection through XP-EHH analysis, have previously been linked to growth traits and fatty acid composition in multiple studies [54,55,56,57,58,59]. The change in the selection index, which took place in 2003 by incorporating muscularity and mammary traits, has not led to discernible genetic changes between 2000 and 2018, possibly due to the brief interval between the two sampling periods, but also to the lower weighting of these traits in the selection index as compared to milk yield.
Additional signatures of selection were identified near genes associated with the immune system, in particular TNFRSF1B on chromosome 16, which is known for its involvement in immune response and was previously observed to be under selection in Russian native cattle breeds [60]. In addition, the HIVEP1 gene, identified through XP-EHH analysis between Rendena and modern Brown Swiss, has been linked to susceptibility to paratuberculosis infection in Chinese Holstein and US Jersey cattle [61,62,63] and may also play a role in trypanotolerance in Zebu Crossbred cattle in Burkina Faso [64]. Moreover, a consistently homozygous region on chromosome 22 encompassed eight genes involved in adaptive immune responses (Fig. 6; see Additional file 8: Table S4) [65]. The large number of immune response genes under selection may be the result of the Rendena having survived a series of historical epidemics that affected the Alps over several centuries. This hypothesis finds a further confirmation in the signatures of selection related to disease resistance and immunity that were found in other breeds from the Alpine arc [58, 66, 67].
Interestingly, a signature of selection was also identified in a region on chromosome 10 that harbours six genes associated with heat tolerance, lipid and energy metabolism, and adaptation to high-altitude and cold environments [68,69,70,71] (Fig. 6). Lastly, the RCL1 gene on chromosome 8 was also found to be under positive selection in another Alpine breed, the Valdostana Black Pied [58]. These findings are particularly notable because Rendena has been bred for centuries to thrive in Alpine valleys, where animals are typically kept outdoors during the summer months. These animals graze on steep pastures at altitudes that can, in some cases, exceed 2000 m above mean sea level, often facing extreme weather conditions which could potentially impact their productivity and health.
Conclusions
Although sharing a common origin with breeds descending from the Brown Swiss group, Rendena has retained distinctive genetic characteristics. For example, disease resistance and adaptation to local food sources and climatic conditions appear to have evolved in response to environmental challenges, and were preserved thanks to the dedication of local breeders. Retention of these traits is fundamental to maintaining this reservoir of genetic diversity, to maximize the capacity of this traditional livestock system to effectively adapt to a changing environment. Genetic signatures of selection in contemporary Rendena are likely a consequence of human pressure to improve productivity and extrinsic challenges, such as pathogens. With more dense marker data or whole-genome sequences, future research could further explore the functional significance of the genes associated with immune response and environmental adaptation in Rendena cattle, to better elucidate mechanisms underlying disease resistance and climate resilience. Investigating the genomic basis of local adaptation may inform targeted breeding strategies that are aimed at enhancing the sustainability and productivity of indigenous breeds in the face of changing and challenging conditions. In addition, comparing data from the same population sampled over time proved a valuable tool for assessing demographic trajectories and monitoring the risk of an increase in inbreeding, since local breeds that are characterized by small population sizes and regional distributions, are particularly vulnerable to increases in inbreeding and the loss of specific adaptive features. In this context, molecular data is of fundamental importance for planning breeding strategies to maintain diversity in the population, while selecting for production and adaptative traits. Moreover, locally-adapted breeds such as the Rendena breed, represent substantial social and economic value to local communities: this breed is not only a source of cultural legacy and tradition, but also plays a role in sustaining the livelihoods of the human population, contributing to the local economy not only through dairy and meat production, but also through tourism. Given the distinctive genetic heritage of the Rendena breed and its adaptation to the Alpine environment, developed over centuries of human selection, along with its substantial social, economic, and cultural importance, preserving this breed is essential for ensuring the resilience and sustainability of both the livestock system and the economy of the Rendena valley.
Availability of data and materials
The datasets used during the current study are available from the corresponding author on reasonable request.
References
Mantovani R, Gallo L, Carnier, P, Cassandro M, Bittante G. The use of a juvenile selection scheme for genetic improvement of small populations: the example of Rendena breed. In: Proceedings of the 48th Annual Meeting of the European Association for Animal Production: 25–28 August 1997; Vienna. 1997.
Coia V, Boschi I, Trombetta F, Cavulli F, Montinaro F, Destro-Bisol G, et al. Evidence of high genetic variation among linguistically diverse populations on a micro-geographic scale: A case study of the Italian Alps. J Hum Genet. 2012;57:254–60.
Cavada E. Storia del Trentino. Il territorio: popolamento, abitati, necropoli. Vol. II. L'età romana. Bologna: Il Mulino. 2000.
Lanzinger M, Marzatico F, Pedrotti A. Storia del Trentino. Vol. I. La preistoria e la protostoria. Bologna: Il Mulino; 2000.
Senczuk G, Mastrangelo S, Ciani E, Battaglini L, Cendron F, Ciampolini R, et al. The genetic heritage of Alpine local cattle breeds using genomic SNP data. Genet Sel Evol. 2020;52:40.
Bonsembiante M, Bittante G, Ramanzin M, Neri C. Caratteristiche, evoluzione e miglioramento della razza rendena. Padova: Pragmark; 1988.
Barbato M, Hailer F, Upadhyay M, Del Corvo M, Colli L, Negrini R, et al. Adaptive introgression from indicine cattle into white cattle breeds from Central Italy. Sci Rep. 2020;10:1279.
Sempéré G, Moazami-Goudarzi K, Eggen A, Laloë D, Gautier M, Flori L. WIDDE: A Web-Interfaced next generation database for genetic diversity exploration, with a first application in cattle. BMC Genomics. 2015;16:940.
Rosen BD, Bickhart DM, Schnabel RD, Koren S, Elsik CG, Tseng E, et al. De novo assembly of the cattle reference genome with single-molecule sequencing. Gigascience. 2020;9:giaa021.
Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. 2015;4:7.
Excoffier L, Lischer HEL. Arlequin suite ver3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10:564–7.
Barbato M, Orozco-terWengel P, Tapio M, Bruford MW. SNeP: a tool to estimate trends in recent effective population size trajectories using genome-wide SNP data. Front Genet. 2015;6:109.
Santiago E, Caballero A, Köpke C, Novo I. Estimation of the contemporary effective population size from SNP data while accounting for mating structure. Mol Ecol Resour. 2024;24: e13890.
Biscarini F, Cozzi P, Gaspa G, Marras G. detectRUNS: Detect runs of homozygosity and runs of heterozygosity in diploid genomes. R package version 0.9.6. 2019. https://CRAN.R-project.org/package=detectRUNS/. Accessed 17 May 2021.
Purfield DC, Berry DP, McParland S, Bradley DG. Runs of homozygosity and population history in cattle. BMC Genet. 2012;13:70.
R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2023.
Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009;19:1655–64.
Milanesi M, Capomaccio S, Vajana E, Bomba L, Garcia JF, Ajmone-Marsan P, et al. BITE: an R package for biodiversity analyses. bioRxiv. 2017. https://www.biorxiv.org/content/early/2017/08/29/181610
Huson DH, Bryant D. Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006;23:254–67.
Pickrell JK, Pritchard JK. Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet. 2012;8: e1002967.
Jombart T, Ahmed I. adegenet 1.3-1: New tools for the analysis of genome-wide SNP data. Bioinformatics. 2011;27:3070–1.
Selli A, Ventura RV, Fonseca PAS, Buzanskas ME, Andrietta LT, Balieiro JCC, et al. Detection and visualization of heterozygosity-rich regions and runs of homozygosity in worldwide sheep populations. Animals (Basel). 2021;11:2996.
Fonseca PAS, Suárez-Vega A, Marras G, Cánovas Á. GALLO: An R package for genomic annotation and integration of multiple data sources in livestock for positional candidate loci. Gigascience. 2020;9:giaa149.
Gautier M, Klassmann A, Vitalis R. rehh 2.0: a reimplementation of the R package rehh to detect positive selection from haplotype structure. Mol Ecol Resour. 2017;17:78–90.
Naji MM, Utsunomiya YT, Sölkner J, Rosen BD, Mészáros G. Investigation of ancestral alleles in the Bovinae subfamily. BMC Genomics. 2021;22:108.
Sabeti PC, Reich DE, Higgins JM, Levine HZP, Richter DJ, Schaffner SF, et al. Detecting recent positive selection in the human genome from haplotype structure. Nature. 2002;419:832–7.
Taylor WE, Gibbs EPJ, Bandyopadhyay SK, Pastoret PP, Atang P. Rinderpest and its eradication. Rome: FAO/OIE; 2022.
Simčič M, Medugorac I. A genome-assisted Moet design for inbreeding prevention in the endangered Murnau-Werdenfelsed cattle. In: Proceedings of the 24th International Symposium Animal Science Days: 21–23 September 2016; Ptuj. 2016.
Mancin E, Tuliozi B, Pegolo S, Sartori C, Mantovani R. Genome wide association study of beef traits in local alpine breed reveals the diversity of the pathways involved and the role of time stratification. Front Genet. 2022;12: 746665.
Cohen-Zinder M, Seroussi E, Larkin DM, Loor JJ, Everts-Van Der Wind A, Lee JH, et al. Identification of a missense mutation in the bovine ABCG2 gene with a major effect on the QTL on chromosome 6 affecting milk yield and composition in Holstein cattle. Genome Res. 2005;15:936–44.
Olsen HG, Nilsen H, Hayes B, Berg PR, Svendsen M, Lien S, et al. Genetic support for a quantitative trait nucleotide in the ABCG2 gene affecting milk composition of dairy cattle. BMC Genet. 2007;8:32.
Ron M, Cohen-Zinder M, Peter C, Weller JI, Erhardt G. Short communication: a polymorphism in ABCG2 in Bos indicus and Bos taurus cattle breeds. J Dairy Sci. 2006;89:4921–3.
Pedrosa VB, Schenkel FS, Chen SY, Oliveira HR, Casey TM, Melka MG, et al. Genomewide association analyses of lactation persistency and milk production traits in Holstein cattle based on imputed whole-genome sequence data. Genes (Basel). 2021;12:1830.
Sanchez MP, Govignon-Gion A, Croiseau P, Fritz S, Hozé C, Miranda G, et al. Within-breed and multi-breed GWAS on imputed whole-genome sequence variants reveal candidate mutations affecting milk protein composition in dairy cattle. Genet Sel Evol. 2017;49:68.
Do DN, Bissonnette N, Lacasse P, Miglior F, Sargolzaei M, Zhao X, et al. Genome-wide association analysis and pathways enrichment for lactation persistency in Canadian Holstein cattle. J Dairy Sci. 2017;100:1955–70.
Alberto FJ, Boyer F, Orozco-TerWengel P, Streeter I, Servin B, De Villemereuil P, et al. Convergent genomic signatures of domestication in sheep and goats. Nat Commun. 2018;9:813.
Xia J, Fan H, Chang T, Xu L, Zhang W, Song Y, et al. Searching for new loci and candidate genes for economically important traits through gene-based association analysis of Simmental cattle. Sci Rep. 2017;7:42048.
Bovo S, Schiavo G, Kazemi H, Moscatelli G, Ribani A, Ballan M, et al. Exploiting within-breed variability in the autochthonous Reggiana breed identified several candidate genes affecting pigmentation-related traits, stature and udder defects in cattle. Anim Genet. 2021;52:579–97.
Raza SHA, Khan S, Amjadi M, Abdelnour SA, Ohran H, Alanazi KM, et al. Genome-wide association studies reveal novel loci associated with carcass and body measures in beef cattle. Arch Biochem Biophys. 2020;694: 108543.
Ramirez-Diaz J, Cenadelli S, Bornaghi V, Bongioni G, Montedoro SM, Achilli A, et al. Identification of genomic regions associated with total and progressive sperm motility in Italian Holstein bulls. J Dairy Sci. 2023;106:407–20.
Lindholm-Perry AK, Sexten AK, Kuehn LA, Smith TPL, King DA, Shackelford SD, et al. Association, effects and validation of polymorphisms within the NCAPG - LCORL locus located on BTA6 with feed intake, gain, meat and carcass traits in beef cattle. BMC Genet. 2011;12:103.
Niu Q, Zhang T, Xu L, Wang T, Wang Z, Zhu B, et al. Integration of selection signatures and multi-trait GWAS reveals polygenic genetic architecture of carcass traits in beef cattle. Genomics. 2021;11:3325–36.
An B, Xu L, Xia J, Wang X, Miao J, Chang T, et al. Multiple association analysis of loci and candidate genes that regulate body size at three growth stages in Simmental beef cattle. BMC Genet. 2020;21:32.
Bertolini F, Schiavo G, Bovo S, Sardina MT, Mastrangelo S, Dall’Olio S, et al. Comparative selection signature analyses identify genomic footprints in Reggiana cattle, the traditional breed of the Parmigiano-Reggiano cheese production system. Animal. 2020;14:921–32.
Bongiorni S, Mancini G, Chillemi G, Pariset L, Valentini A. Identification of a short region on chromosome 6 affecting direct calving ease in Piedmontese Cattle Breed. PLoS One. 2012;7: e50137.
Xu L, Bickhart DM, Cole JB, Schroeder SG, Song J, Van Tassell CP, et al. Genomic signatures reveal new evidences for selection of important traits in domestic cattle. Mol Biol Evol. 2015;32:711–25.
Zheng X, Ju Z, Wang J, Li Q, Huang J, Zhang A, et al. Single nucleotide polymorphisms, haplotypes and combined genotypes of LAP3 gene in bovine and their association with milk production traits. Mol Biol Rep. 2011;38:4053–61.
Seo M, Lee HJ, Kim K, Caetano-Anolles K, Jeong JY, Park S, et al. Characterizing milk production related genes in Holstein using RNA-seq. Asian-Australas J Anim Sci. 2016;29:343–51.
Nayeri S, Sargolzaei M, Abo-Ismail MK, May N, Miller SP, Schenkel F, et al. Genome-wide association for milk production and female fertility traits in Canadian dairy Holstein cattle. BMC Genet. 2016;17:75.
Wagh K, Bhatia A, Alexe G, Reddy A, Ravikumar V, Seiler M, et al. Lactase persistence and lipid pathway selection in the Maasai. PLoS One. 2012;7: e44751.
Serranito B, Cavalazzi M, Ciani E, Taurisson-Mouret D, Ciani E, Rouvellac E, et al. Local adaptations of Mediterranean sheep and goats through an integrative approach. Sci Rep. 2021;11:21363.
Brunes LC, Baldi F, Lopes FB, Lôbo RB, Espigolan R, Costa MFO, et al. Weighted single-step genome-wide association study and pathway analyses for feed efficiency traits in Nellore cattle. J Anim Breed Genet. 2021;138:23–44.
O’Grady GL, Best HA, Sztal TE, Schartner V, Sanjuan-Vazquez M, Donkervoort S, et al. Variants in the oxidoreductase PYROXD1 cause early-onset myopathy with internalized nuclei and myofibrillar disorganization. Am J Hum Genet. 2016;99:1086–105.
Feitosa FLB, Pereira ASC, Mueller LF, de Fonseca PAS, Braz CU, Amorin S, et al. Genome-wide association study for beef fatty acid profile using haplotypes in Nellore cattle. Livest Sci. 2021;245:104396.
Ge F, Jia C, Chu M, Liang C, Yan P. Copy number variation of the CADM2 gene and its association with growth traits in Yak. Animals (Basel). 2019;9:1008.
Li X, Huang C, Liu M, Dai R, Wu X, Ma X, et al. Copy number variation of the SOX6 gene and its associations with growth traits in Ashidan Yak. Animals (Basel). 2022;12:3074.
Xu Z, Wang X, Zhang Z, An Q, Wen Y, Wang D, et al. Copy number variation of CADM2 gene revealed its association with growth traits across Chinese Capra hircus (goat) populations. Gene. 2020;741: 144519.
Strillacci MG, Vevey M, Blanchet V, Mantovani R, Sartori C, Bagnato A. The genomic variation in the Aosta cattle breeds raised in an extensive alpine farming system. Animals (Basel). 2020;10:2385.
Fan R, Gu Z, Guang X, Marín JC, Varas V, González BA, et al. Genomic analysis of the domestication and post-Spanish conquest evolution of the llama and alpaca. Genome Biol. 2020;21:159.
Zinovieva NA, Dotsev AV, Sermyagin AA, Deniskova TE, Abdelmanova AS, Kharzinova VR, et al. Selection signatures in two oldest Russian native cattle breeds revealed using high-density single nucleotide polymorphism analysis. PLoS One. 2020;15: e0242200.
Gao Y, Jiang J, Yang S, Cao J, Han B, Wang Y, et al. Genome-wide association study of Mycobacterium avium subspecies Paratuberculosis infection in Chinese Holstein. BMC Genomics. 2018;19:972.
Zare Y, Shook GE, Collins MT, Kirkpatrick BW. Genome-wide association analysis and genomic prediction of Mycobacterium avium subspecies paratuberculosis infection in US Jersey cattle. PLoS One. 2014;9: e88380.
Seong J, Yoon H, Kong HS. Identification of microRNA and target gene associated with marbling score in Korean cattle (Hanwoo). Genes Genomics. 2016;38:529–38.
Yougbaré B, Ouédraogo D, Tapsoba ASR, Soudré A, Zoma BL, Orozco-terWengel P, et al. Local ancestry to identify selection in response to trypanosome infection in Baoulé x Zebu crossbred cattle in Burkina Faso. Front Genet. 2021;12: 670390.
Goyache F, Pérez-Pardal L, Fernández I, Traoré A, Menéndez-Arias NA, Álvarez I. Ancient autozygous segments subject to positive selection suggest adaptive immune responses in West African cattle. Gene. 2021;803: 145899.
Moravíková N, Simčič M, Mészáros G, Sölkner J, Kukučková V, Vlček M, et al. Genomic response to natural selection within alpine cattle breeds. Czech J Anim Sci. 2018;63:136–43.
Signer-Hasler H, Casanova L, Barenco A, Maitre B, Bagnato A, Vevey M, et al. Genomic regions underlying positive selection in local, Alpine cattle breeds. Anim Genet. 2023;54:239–53.
Cunningham HC, Cammack KM, Hales KE, Freetly HC, Lindholm-Perry AK. Differential transcript abundance in adipose tissue of mature beef cows during feed restriction and realimentation. PLoS One. 2018;13: e0194104.
Ge Q, Guo Y, Zheng W, Zhao S, Cai Y, Qi X. Molecular mechanisms detected in yak lung tissue via transcriptome-wide analysis provide insights into adaptation to high altitudes. Sci Rep. 2021;11:7786.
Huang H, Cao J, Hanif Q, Wang Y, Yu Y, Zhang S, et al. Genome-wide association study identifies energy metabolism genes for resistance to ketosis in Chinese Holstein cattle. Anim Genet. 2019;50:376–80.
Weldenegodguad M, Pokharel K, Niiranen L, Soppela P, Ammosov I, Honkatukia M, et al. Adipose gene expression profiles reveal insights into the adaptation of northern Eurasian semi-domestic reindeer (Rangifer tarandus). Commun Biol. 2021;4:1170.
Acknowledgements
The authors would like to thank the numerous farmers who actively participated in the project and granted permission for genotyping their animals. They also extend their acknowledgement to Federazione Provinciale Allevatori and ANARE, with their Director, Dr. Dario Tonietto, for their encouragement and advice and for access to pedigree records. The authors also thank the two reviewers for useful suggestions to improve this paper.
Funding
This research was partially funded by the Project “Azioni per l’individuazione, caratterizzazione, divulgazione della biodiversità agricola e alimentare in Provincia di Trento”, Bando Ministero delle politiche agricole alimentari, forestali e del turismo (prot. 31943 del 16/11/2018) and partially by Università Cattolica own funds. NF received a PhD Scholarship from the Fondazione E. Mach and the Università Cattolica del Sacro Cuore, Piacenza.
Author information
Authors and Affiliations
Contributions
EP, NF, MK, AA, HCH, PAM developed the research idea and sampling design and acquired funding; EP, NF, MK were involved in sample collection; ES, MB, LC performed the data analyses; ES, MB, AMCO, LC, RM, FP, AA, HCH and PAM performed the data interpretation. ES wrote the first draft and all authors contributed to editing of the final manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
All samples were analysed with the explicit written consent of the owners and consisted of subsamples of blood taken during scheduled screening regimes; in addition, the genetic characteristics of each herd were summarized and presented in the form of a short report to each livestock farmer.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1:Table S1.
Breed acronyms, extended breed names, SNP chip used and data source.
Additional file 2: Table S2.
Breed acronyms, FIS values and associated p-value.
Additional file 3: Table S3.
Contemporary Ne estimates computed with the currentNe software. The table reports breed acronyms, number of samples (N samples) in the working dataset, contemporary Ne computed with currentNe (S-Ne), lower bounds of 90% confidence interval (CI lower), upper bounds of 90% confidence interval (CI upper).
Additional file 4: Figure S1.
Principal component analysis (PC1 vs. PC3). The variance explained by each component is given as percentage in brackets. See Additional file 1: Table S1 for breed abbreviations.
Additional file 5: Figure S2.
Circular plot of Admixture results for K from 2 to 30 displayed using a colour-blind friendly palette. Breeds are ordered accordi g to K = 2 values. See Additional file 1: Table S1 for breed abbreviations.
Additional file 6: Figure S3.
Cross-validation error values for K values from 2 to 30.
Additional file 7: Figure S4.
Discriminant analysis of principal component. Discriminant function 1 computed for the Rendena2000 and Rendena2018 populations.
Additional file 8: Table S4.
For each bovine chromosome, are listed the genes annotated for the most represented SNPs in the ROH analysis for the RENgen individuals.
Additional file 9: Figure S5.
iHS results for the Rendena breed. In grey the analysis performed with all SNPs, in red the analysis performed using only SNPs for which the information on the ancestral allele was available. The dotted red and grey lines represent the significance threshold for iHS with and without ancestral allele information, respectively.
Additional file 10: Table S5.
Chromosome, marker, genomic position, XP-EHH values in the top 1% scored in the comparison between RENgen and Original Braunvieh data, RENgen and Brown Swiss cluster and annotated genes. In bold are markers scoring a value of XP-EHH above the threshold in both comparisons.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Somenzi, E., Partel, E., Barbato, M. et al. Genetic legacy and adaptive signatures: investigating the history, diversity, and selection signatures in Rendena cattle resilient to eighteenth century rinderpest epidemics. Genet Sel Evol 56, 32 (2024). https://doi.org/10.1186/s12711-024-00900-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12711-024-00900-y