Selection signatures in two oldest Russian native cattle breeds revealed using high-density single nucleotide polymorphism analysis

Tracing the Dynamical Genetic Diversity Changes of Russian Livni Pigs during the Last 50 Years with the Museum, Old, and Modern Samples

The pig industry is usually considered an intensive livestock industry, mainly supported by hybrid breeding between commercial pig breeds. However, people's pursuit of a more natural environment and higher meat quality has led to an increasing demand for eco-friendly and diverse pig feeding systems. Therefore, the importance of rearing and conserving local pig breeds is increasing. The Livni pig is a local breed with good adaptability to the environmental and fodder conditions in central Russia. In this study, we aimed to analyze the genetic diversity and population structure of Livni pigs using whole-genome single nucleotide polymorphism (SNP) data. We utilized the Porcine GGP HD BeadChip on genotype samples from old (n = 32, 2004) and modern (n = 32, 2019) populations of Livni pigs. For the museum samples of Livni pigs (n = 3), we extracted DNA from their teeth, performed genomic sequencing, and obtained SNP genotypes from the whole-genome sequences. SNP genotypes of Landrace (n = 32) and Large White (n = 32) pigs were included for comparative analysis. We observed that the allelic richness of Livni pigs was higher than those of Landrace and Large White pigs (AR = 1.775-1.798 vs. 1.703 and 1.668, respectively). The effective population size estimates (NE5 = 108 for Livni pigs, NE5 = 59 for Landrace and Large White pigs) confirmed their genetic diversity tendency. This was further supported by the length and number of runs of homozygosity, as well as the genomic inbreeding coefficient (almost twofold lower in Livni pigs compared to Landrace and Large White pigs). These findings suggest that the Livni pig population exhibits higher genetic diversity and experiences lower selection pressure compared to commercial pig populations. Furthermore, both principal component and network tree analyses demonstrated a clear differentiation between Livni pigs and transboundary commercial pigs. The TreeMix results indicated gene flow from Landrace ancestors to Livni pigs (2019) and from Large White ancestors to Livni pigs (2004), which was consistent with their respective historical breeding backgrounds. The comparative analysis of museum, old, and modern Livni pigs indicated that the modern Livni pig populations have preserved their historical genomic components, suggesting their potential suitability for future design selection programs.

Genetic legacy and adaptive signatures: investigating the history, diversity, and selection signatures in Rendena cattle resilient to eighteenth century rinderpest epidemics

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.

Characterization of runs of homozygosity islands in American mink using whole-genome sequencing data

The genome-wide analysis of runs of homozygosity (ROH) islands can be an effective strategy for identifying shared variants within a population and uncovering important genomic regions related to complex traits. The current study performed ROH analysis to characterize the genome-wide patterns of homozygosity, identify ROH islands and annotated genes within these candidate regions using whole-genome sequencing data from 100 American mink (Neogale vison). After sequence processing, variants were called using GATK and Samtools pipelines. Subsequent to quality control, 8,373,854 bi-allelic variants identified by both pipelines remained for further analysis. A total of 34,652 ROH segments were identified in all individuals, among which shorter segments (0.3-1 Mb) were abundant throughout the genome, approximately accounting for 84.39% of all ROH. Within these segments, we identified 63 ROH islands housing 156 annotated genes. The genes located in ROH islands were associated with fur quality (EDNRA, FGF2, FOXA2 and SLC24A4), body size/weight (MYLK4, PRIM2, FABP2, EYS and PHF3), immune capacity (IL2, IL21, PTP4A1, SEMA4C, JAK2, CCNA2 and TNIP3) and reproduction (ADAD1, KHDRBS2, INSL6, PGRMC2 and HSPA4L). Furthermore, Gene Ontology and KEGG pathway enrichment analyses revealed 56 and 9 significant terms (FDR-corrected p-value < 0.05), respectively, among which cGMP-PKG signalling pathway, regulation of actin cytoskeleton, and calcium signalling pathway were highlighted due to their functional roles in growth and fur characteristics. This is the first study to present ROH islands in American mink. The candidate genes from ROH islands and functional enrichment analysis suggest possible signatures of selection in response to the mink breeding targets, such as increased body length, reproductive performance and fur quality. These findings contribute to our understanding of genetic characteristics, and provide complementary information to assist with implementation of breeding strategies for genetic improvement in American mink.

Association of three single nucleotide polymorphisms in the LPIN1 gene with milk production traits in cows of the Yaroslavl breed

Lipin-1 is a member of the evolutionarily conserved family of proteins and is expressed predominantly in adipose tissue and skeletal muscle. On the one hand, lipin-1 is an enzyme that catalyzes the dephosphorylation of phosphatidic acid to diacylglycerol (DAG) and thus participates in the metabolic pathways of biosynthesis of storage lipids in the cell, membrane phospholipids, and intracellular signaling molecules. On the other hand, lipin-1 is able to be transported from the cytoplasm to the nucleus and is a coactivator of lipid metabolism gene transcription. It was shown, using the analysis of single nucleotide polymorphism (SNP) associations, that the lipin-1 coding gene (LPIN1) is a promising candidate gene for milk production traits in Holstein and Brown Swiss cows. However, it is unclear how much of its effect depends on the breed. The Yaroslavl dairy cattle breed was created in the 18-19 centuries in Russia by breeding northern Great Russian cattle, which were short and poor productive, but well adapted to local climatic conditions and bad food base. It was shown by whole genome genotyping and sequencing that the Yaroslavl breed has unique genetics compared to Russian and other cattle breeds. The aim of the study was to assess the frequency of alleles and genotypes of three SNPs in the LPIN1 gene and to study the association of these SNPs with milk production traits in Yaroslavl cows. Blood samples from 142 cows of the Yaroslavl breed were obtained from two farms in the Yaroslavl region. Genotyping of SNPs was carried out by polymerase chain reaction-restriction fragment length polymorphism method. Associations of SNPs with 305-day milk yield, fat yield, fat percentages, protein yield, and protein percentages were studied from the first to the fourth lactation. Statistical tests were carried out using a mixed linear model, taking into account the relationship between individuals. We identified three SNPs - rs110871255, rs207681322 and rs109039955 with a frequency of a rare allele of 0.042-0.261 in Yaroslavl cows. SNP rs110871255 was associated with fat yield during the third and fourth lactations. SNP rs207681322 was associated with milk yield for the second, third and fourth lactations, as well as protein yield for the third lactation. Thus, we identified significant associations of SNPs rs207681322 and rs110871255 in the LPIN1 gene with a number of milk production traits during several lactations in Yaroslavl cows.

Genetic diversity and positive signatures of selection in indigenous cattle breeds of Iran

Animal domestication, climate changes over time, and artificial selection have played significant roles in shaping the genome structure of various animal species, including cattle. These processes have led to the emergence of several indigenous cattle breeds with distinct genetic characteristics. This study focused on unraveling the genetic diversity and identifying candidate genomic regions in eight indigenous cattle breeds of Iran. The data consisted of ∼777 962 single nucleotide polymorphisms (SNPs) of 89 animals from Iranian indigenous cattle scattered throughout the country. We employed various methods, including integrated haplotype score, FST, and cross-population composite likelihood ratio, to conduct a genome scan for detecting selection signals within and between cattle populations. Average observed heterozygosity across the populations was 0.36, with a range of 0.32-0.40. In addition, negative and low rates of inbreeding (FIS) in the populations were observed. The genome-wide analysis revealed several genomic regions that harbored candidate genes associated with production traits (e.g., MFSD1, TYW5, ADRB2, BLK, and CRTC3), adaptation to local environmental constraints (CACNA2D1, CXCL3, and GRO1), and coat color (DYM). Finally, the study of the reported quantitative trait loci (QTL) regions in the cattle genome demonstrated that the identified regions were associated with QTL related to important traits such as milk composition, body weight, daily gain, feed conversion, and residual feed intake. Overall, this study contributes to a better understanding of the genetic diversity and potential candidate genes underlying important traits in Iranian indigenous cattle breeds, which can inform future breeding and conservation efforts.

Selective Sweeps in Cattle Genomes in Response to the Influence of Urbanization and Environmental Contamination

A genomic study was conducted to identify the effects of urbanization and environmental contaminants with heavy metals on selection footprints in dairy cattle populations reared in the megacity of Bengaluru, South India. Dairy cattle reared along the rural-urban interface of Bengaluru with/without access to roughage from public lakeshores were selected. The genotyped animals were subjected to the cross-population-extended haplotype homozygosity (XP-EHH) methodology to infer selection sweeps caused by urbanization (rural, mixed, and urban) and environmental contamination with cadmium and lead. We postulated that social-ecological challenges contribute to mechanisms of natural selection. A number of selection sweeps were identified when comparing the genomes of cattle located in rural, mixed, or urban regions. The largest effects were identified on BTA21, displaying pronounced peaks for selection sweeps for all three urbanization levels (urban_vs_rural, urban_vs_mixed and rural_vs_mixed). Selection sweeps are located in chromosomal segments in close proximity to the genes lrand rab interactor 3 (RIN3), solute carrier family 24 member 4 (SLC24A4), tetraspanin 3 (TSPAN3), and proline-serine-threonine phosphatase interacting protein 1 (PSTPIP1). Functional enrichment analyses of the selection sweeps for all three comparisons revealed a number of gene ontology (GO) and KEGG terms, which were associated with reproduction, metabolism, and cell signaling-related functional mechanisms. Likewise, a number of the chromosomal segments under selection were observed when creating cattle groups according to cadmium and lead contaminations. Stronger and more intense positive selection sweeps were observed for the cadmium contaminated group, i.e., signals of selection on BTA 16 and BTA19 in close proximity to genes regulating the somatotropic axis (growth factor receptor bound protein 2 (GRB2) and cell ion exchange (chloride voltage-gated channel 6 (CLCN6)). A few novel, so far uncharacterized genes, mostly with effects on immune physiology, were identified. The lead contaminated group revealed sweeps which were annotated with genes involved in carcass traits (TNNC2, SLC12A5, and GABRA4), milk yield (HTR1D, SLCO3A1, TEK, and OPCML), reproduction (GABRA4), hypoxia/stress response (OPRD1 and KDR), cell adhesion (PCDHGC3), inflammatory response (ADORA2A), and immune defense mechanism (ALCAM). Thus, the findings from this study provide a deeper insight into the genomic regions under selection under the effects of urbanization and environmental contamination.

Identification of homozygosity-rich regions in the Holstein genome

In this study, 371 Holstein cows from six herds and 26 Holstein bulls, which were used in these herds, were genotyped by the Illumina BovineSNP50 array. For runs of homozygosity (ROH) identification, consecutive and sliding runs were performed by the detectRUNS and Plink software. The missing calls did not significantly affect the ROH data. The mean number of ROH identified by consecutive runs was 95.4 ± 2.7, and that by sliding runs was 86.0 ± 2.6 in cows, while this number for Holstein bulls was lower 58.9 ± 1.9. The length of the ROH segments varied from 1 Mb to over 16 Mb, with the largest number of ROH having a length of 1-2 Mb. Of the 29 chromosomes, BTA 14, BTA 16, and BTA 7 were the most covered by ROH. The mean coefficient of inbreeding across the herds was 0.111 ± 0.003 and 0.104 ± 0.004 based on consecutive and sliding runs, respectively, and 0.078 ± 0.005 for bulls based on consecutive runs. These values do not exceed those for Holstein cattle in North America. The results of this study confirmed the more accurate identification of ROH by consecutive runs, and also that the number of allowed heterozygous SNPs may have a significant effect on ROH data.

Genetic Structure Analysis of 155 Transboundary and Local Populations of Cattle (Bos taurus, Bos indicus and Bos grunniens) Based on STR Markers

Every week, 1-2 breeds of farm animals, including local cattle, disappear in the world. As the keepers of rare allelic variants, native breeds potentially expand the range of genetic solutions to possible problems of the future, which means that the study of the genetic structure of these breeds is an urgent task. Providing nomadic herders with valuable resources necessary for life, domestic yaks have also become an important object of study. In order to determine the population genetic characteristics, and clarify the phylogenetic relationships of modern representatives of 155 cattle populations from different regions of the world, we collected a large set of STR data (10,250 individuals), including unique native cattle, 12 yak populations from Russia, Mongolia and Kyrgyzstan, as well as zebu breeds. Estimation of main population genetic parameters, phylogenetic analysis, principal component analysis and Bayesian cluster analysis allowed us to refine genetic structure and provided insights in relationships of native populations, transboundary breeds and populations of domestic yak. Our results can find practical application in conservation programs of endangered breeds, as well as become the basis for future fundamental research.

Genome-Wide Association Study for Body Conformation Traits and Fitness in Czech Holsteins

The aim of this study was a genome-wide association study (GWAS) on conformation traits using 25,486 genotyped Czech Holsteins, with 35,227 common SNPs for each genotype. Linear trait records were collected between 1995 and 2020. The Interbull information from Multiple Across Country Evaluation (MACE) was included for bulls that mostly had daughter records in a foreign country. When using the Bonferroni correction, the number of SNPs that were either significant or approached the significance threshold was low-dairy capacity composite on BTA4, feet and legs composite BTA21, total score BTA10, stature BTA24, body depth BTA6, angularity BTA20, fore udder attachment BTA10. Without the Bonferroni correction, the total number of significant or near of significance SNPs was 32. The SNPs were localized on BTA1,2,4,5,6,7,8,18,22,25,26,28 for dairy capacity composite, BTA15,21 for feet and legs composite, BTA10 for total score, BTA24 stature, BTA6,23 body depth, BTA20 angularity, BTA2 rump angle, BTA9,10 rear legs rear view, BTA2,19 rear legs side view, BTA10 fore udder attachment, BTA2 udder depth, BTA10 rear udder height, BTA12 central alignment, BTA24 rear teat placement, BTA8,29 rear udder width. The results provide biological information for the improvement of body conformation and fitness in the Holstein population.

Selection and Drift: A Comparison between Historic and Recent Dutch Friesian Cattle and Recent Holstein Friesian Using WGS Data

Simple Summary

Over the last century, genetic diversity in the cattle species has been affected by the replacement of many local, dual-purpose breeds with a few specialized, high-output dairy breeds. This replacement caused a sharp decline in the population size of local breeds. In the Netherlands, the local Dutch Friesian breed has gradually been replaced by the Holstein Friesian. This resulted in a rapid decrease in numbers of the Dutch Friesian breed with an associated risk of loss of genetic diversity due to drift. The objective of this study is to investigate genomewide genetic diversity between a group of historic and recent Dutch Friesian bulls and a group of recently used Holstein Friesian bulls. Our findings showed that a large amount of diversity is shared between the three groups, but each of them has some unique genetic identity (12% of the single nucleotide polymorphism were group-specific). The genetic diversity of the Dutch Friesians reduced over time, but this did not lead to higher inbreeding levels—especially, inbreeding due to recent ancestors has not increased. Genetically, the recent Dutch Friesians were slightly more different from Holstein Friesians than the historic Dutch Friesians. Our results also highlighted the presence of several genomic regions that differentiated between the groups.

Abstract

Over the last century, genetic diversity in many cattle breeds has been affected by the replacement of traditional local breeds with just a few milk-producing breeds. In the Netherlands, the local Dutch Friesian breed (DF) has gradually been replaced by the Holstein Friesian breed (HF). The objective of this study is to investigate genomewide genetic diversity between a group of historically and recently used DF bulls and a group of recently used HF bulls. Genetic material of 12 historic (hDF), 12 recent DF bulls (rDF), and 12 recent HF bulls (rHF) in the Netherlands was sequenced. Based on the genomic information, different parameters—e.g., allele frequencies, inbreeding coefficient, and runs of homozygosity (ROH)—were calculated. Our findings showed that a large amount of diversity is shared between the three groups, but each of them has a unique genetic identity (12% of the single nucleotide polymorphisms were group-specific). The rDF is slightly more diverged from rHF than hDF. The inbreeding coefficient based on runs of homozygosity (Froh) was higher for rDF (0.24) than for hDF (0.17) or rHF (0.13). Our results also displayed the presence of several genomic regions that differentiated between the groups. In addition, thirteen, forty-five, and six ROH islands were identified in hDF, rDF, and rHF, respectively. The genetic diversity of the DF breed reduced over time, but this did not lead to higher inbreeding levels—especially, inbreeding due to recent ancestors was not increased.

Genomic Diversity Profiling and Breed-Specific Evolutionary Signatures of Selection in Arunachali Yak

Arunachali yak, the only registered yak breed of India, is crucial for the economic sustainability of pastoralist Monpa community. This study intended to determine the genomic diversity and to identify signatures of selection in the breed. Previously available double digest restriction-site associated DNA (ddRAD) sequencing data of Arunachali yak animals was processed and 99,919 SNPs were considered for further analysis. The genomic diversity profiled based on nucleotide diversity, π (π = 0.041 in 200 bp windows), effective population size, Ne (Ne = 83) and Runs of homozygosity (ROH) (predominance of shorter length ROHs) was found to be optimum. Subsequently, 207 regions were identified to be under selective sweeps through de-correlated composite of multiple signals (DCMS) statistic which combined three individual test statistics viz. π, Tajima’s D and |iHS| in non-overlapping 100 kb windows. Mapping of these regions revealed 611 protein-coding genes including KIT, KITLG, CDH12, FGG, FGA, FGB, PDGFRA, PEAR1, STXBP3, olfactory receptor genes (OR5K3, OR5H6 and OR1E1) and taste receptor genes (TAS2R1, TAS2R3 and TAS2R4). Functional annotation highlighted that biological processes like platelet aggregation and sensory perception were the most overrepresented and the associated regions could be considered as breed-specific signatures of selection in Arunachali yak. These findings point towards evolutionary role of natural selection in environmental adaptation of Arunachali yak population and provide useful insights for pursuing genome-wide association studies in future.

Unveiling Comparative Genomic Trajectories of Selection and Key Candidate Genes in Egg-Type Russian White and Meat-Type White Cornish Chickens

Comparison of genomic footprints in chicken breeds with different selection history is a powerful tool in elucidating genomic regions that have been targeted by recent and more ancient selection. In the present work, we aimed at examining and comparing the trajectories of artificial selection in the genomes of the native egg-type Russian White (RW) and meat-type White Cornish (WC) breeds. Combining three different statistics (top 0.1% SNP by F_ST value at pairwise breed comparison, hapFLK analysis, and identification of ROH island shared by more than 50% of individuals), we detected 45 genomic regions under putative selection including 11 selective sweep regions, which were detected by at least two different methods. Four of such regions were breed-specific for each of RW breed (on GGA1, GGA5, GGA8, and GGA9) and WC breed (on GGA1, GGA5, GGA8, and GGA28), while three remaining regions on GGA2 (two sweeps) and GGA3 were common for both breeds. Most of identified genomic regions overlapped with known QTLs and/or candidate genes including those for body temperatures, egg productivity, and feed intake in RW chickens and those for growth, meat and carcass traits, and feed efficiency in WC chickens. These findings were concordant with the breed origin and history of their artificial selection. We determined a set of 188 prioritized candidate genes retrieved from the 11 overlapped regions of putative selection and reviewed their functions relative to phenotypic traits of interest in the two breeds. One of the RW-specific sweep regions harbored the known domestication gene, TSHR. Gene ontology and functional annotation analysis provided additional insight into a functional coherence of genes in the sweep regions. We also showed a greater candidate gene richness on microchromosomes relative to macrochromosomes in these genomic areas. Our results on the selection history of RW and WC chickens and their key candidate genes under selection serve as a profound information for further conservation of their genomic diversity and efficient breeding.

Biodiversity of Russian Local Sheep Breeds Based on Pattern of Runs of Homozygosity

Russian sheep breeds traditionally raised in specific environments are valuable parts of sociocultural heritage and economic component of the regions. However, the import of commercial breeds negatively influences the population sizes of local sheep populations and might lead to biodiversity loss. Estimation of the runs of homozygosity (ROH) in local sheep genomes is an informative tool to address their current genetic state. In this work, we aimed to address the ROH distribution and to estimate genome inbreeding based on SNP data to evaluate genetic diversity in Russian local sheep breeds. Materials for this study included SNP-genotypes from twenty-seven Russian local sheep breeds which were generated using the Illumina OvineSNP50 BeadChip (n = 391) or the Illumina Ovine Infinium HD BeadChip (n = 315). A consecutive runs method was used to calculate ROH which were estimated for each animal and then categorized in the ROH length classes. The ROH were found in all breeds. The mean ROH length varied from 86 to 280 Mb, while the ROH number ranged from 37 to 123. The genomic inbreeding coefficient varied from 0.033 to 0.106. Our findings provide evidence of low to moderate genomic inbreeding in major local sheep populations.

Genome-wide SNP analysis clearly distinguished the Belarusian Red cattle from other European cattle breeds

Local breeds can serve as an important source of genetic variability in domestic animal species. This study aimed to assess the genetic diversity and population structure of Belarusian Red cattle and their differentiation from other European cattle populations based on genome-wide SNP genotypes. Twenty pedigree-recorded non-closely related cows of Belarusian Red cattle were genotyped using the Illumina BovineHD BeadChip. Genotypes of 22 other European cattle breeds were included in the study for comparison. A total of 28 562 SNPs passed through the quality control checks and were selected for analysis. The Belarusian Red cattle displayed a moderate level of genetic variability (_U H_E  = 0.341, H_O  = 0.368), and the highest heterozygote excess (_U F_IS  = -0.066), among the studied breeds; this reflects the contribution of multiple breeds to their formation. The principal component analysis, F_ST -based Neighbor-Net tree and Admixture clustering, clearly distinguished the Belarusian Red cattle from the other European cattle breeds. Moreover, the presence of ancestral genomic components of Danish Red and Brown Swiss breeds were clearly visible, which agrees with the breed's history and its recent development. Our study highlights the importance of maintaining the specific genomic components, which makes a significant contribution to the global genetic diversity in the modern population of Belarusian Red cattle, allowing us to consider them a valuable national genetic resource. Our research results will be useful for the development of conservation programs for this local cattle breed.

Genome-Wide Analysis Reveals Selection Signatures Involved in Meat Traits and Local Adaptation in Semi-Feral Maremmana Cattle

The Maremmana cattle is an ancient Podolian-derived Italian breed raised in semi-wild conditions with distinctive morphological and adaptive traits. The aim of this study was to detect potential selection signatures in Maremmana using medium-density single nucleotide polymorphism array. Putative selection signatures were investigated combining three statistical approaches designed to quantify the excess of haplotype homozygosity either within (integrated haplotype score, iHS) or among pairs of populations (Rsb and XP-EHH), and contrasting the Maremmana with a single reference population composed of a pool of seven Podolian-derived Italian breeds. Overall, the three haplotype-based analyses revealed selection signatures distributed over 19 genomic regions. Of these, six relevant candidate regions were identified by at least two approaches. We found genomic signatures of selective sweeps spanning genes related to mitochondrial function, muscle development, growth, and meat traits (SCIN, THSD7A, ETV1, UCHL1, and MYOD1), which reflects the different breeding schemes between Maremmana (semi-wild conditions) and the other Podolian-derived Italian breeds (semi-extensive). We also identified several genes linked to Maremmana adaptation to the environment of the western-central part of Italy, known to be hyperendemic for malaria and other tick-borne diseases. These include several chemokine (C-C motif) ligand genes crucially involved in both innate and adaptive immune responses to intracellular parasite infections and other genes playing key roles in pulmonary disease (HEATR9, MMP28, and ASIC2) or strongly associated with malaria resistance/susceptibility (AP2B1). Our results provide a glimpse into diverse selection signatures in Maremmana cattle and can be used to enhance our understanding of the genomic basis of environmental adaptation in cattle.