A comprehensive analysis of the genetic diversity and environmental adaptability in worldwide Merino and Merino-derived sheep breeds.

BACKGROUND: To enhance and extend the knowledge about the global historical and phylogenetic relationships between Merino and Merino-derived breeds, 19 populations were genotyped with the OvineSNP50 BeadChip specifically for this study, while an additional 23 populations from the publicly available genotypes were retrieved. Three complementary statistical tests, Rsb (extended haplotype homozygosity between-populations), XP-EHH (cross-population extended haplotype homozygosity), and runs of homozygosity (ROH) islands were applied to identify genomic variants with potential impact on the adaptability of Merino genetic type in two contrasting climate zones. RESULTS: The results indicate that a large part of the Merino's genetic relatedness and admixture patterns are explained by their genetic background and/or geographic origin, followed by local admixture. Multi-dimensional scaling, Neighbor-Net, Admixture, and TREEMIX analyses consistently provided evidence of the role of Australian, Rambouillet and German strains in the extensive gene introgression into the other Merino and Merino-derived breeds. The close relationship between Iberian Merinos and other South-western European breeds is consistent with the Iberian origin of the Merino genetic type, with traces from previous contributions of other Mediterranean stocks. Using Rsb and XP-EHH approaches, signatures of selection were detected spanning four genomic regions located on Ovis aries chromosomes (OAR) 1, 6 and 16, whereas two genomic regions on OAR6, that partially overlapped with the previous ones, were highlighted by ROH islands. Overall, the three approaches identified 106 candidate genes putatively under selection. Among them, genes related to immune response were identified via the gene interaction network. In addition, several candidate genes were found, such as LEKR1, LCORL, GHR, RBPJ, BMPR1B, PPARGC1A, and PRKAA1, related to morphological, growth and reproductive traits, adaptive thermogenesis, and hypoxia responses. CONCLUSIONS: To the best of our knowledge, this is the first comprehensive dataset that includes most of the Merino and Merino-derived sheep breeds raised in different regions of the world. The results provide an in-depth picture of the genetic makeup of the current Merino and Merino-derived breeds, highlighting the possible selection pressures associated with the combined effect of anthropic and environmental factors. The study underlines the importance of Merino genetic types as invaluable resources of possible adaptive diversity in the context of the occurring climate changes.

Geographical contrasts of Y-chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions.

By their paternal transmission, Y-chromosomal haplotypes are sensitive markers of population history and male-mediated introgression. Previous studies identified biallelic single-nucleotide variants in the SRY, ZFY and DDX3Y genes, which in domestic goats identified four major Y-chromosomal haplotypes, Y1A, Y1B, Y2A and Y2B, with a marked geographical partitioning. Here, we extracted goat Y-chromosomal variants from whole-genome sequences of 386 domestic goats (75 breeds) and seven wild goat species, which were generated by the VarGoats goat genome project. Phylogenetic analyses indicated domestic haplogroups corresponding to Y1B, Y2A and Y2B, respectively, whereas Y1A is split into Y1AA and Y1AB. All five haplogroups were detected in 26 ancient DNA samples from southeast Europe or Asia. Haplotypes from present-day bezoars are not shared with domestic goats and are attached to deep nodes of the trees and networks. Haplogroup distributions for 186 domestic breeds indicate ancient paternal population bottlenecks and expansions during migrations into northern Europe, eastern and southern Asia, and Africa south of the Sahara. In addition, sharing of haplogroups indicates male-mediated introgressions, most notably an early gene flow from Asian goats into Madagascar and the crossbreeding that in the 19th century resulted in the popular Boer and Anglo-Nubian breeds. More recent introgressions are those from European goats into the native Korean goat population and from Boer goat into Uganda, Kenya, Tanzania, Malawi and Zimbabwe. This study illustrates the power of the Y-chromosomal variants for reconstructing the history of domestic species with a wide geographical range.

Exogenous Jaagsiekte Sheep Retrovirus type 2 (exJSRV2) related to ovine pulmonary adenocarcinoma (OPA) in Romania: prevalence, anatomical forms, pathological description, immunophenotyping and virus identification.

BACKGROUND: Ovine pulmonary adenocarcinoma (OPA) is a neoplastic disease caused by exogenous Jaagsiekte Sheep Retrovirus (exJSRV). The prevalence of JSRV-related OPA in Eastern European countries, including Romania is unknown. We aimed to investigate: the prevalence and morphological features of OPA (classical and atypical forms) in the Transylvania region (Romania), the immunophenotype of the pulmonary tumors and their relationships with exJSRV infection. A total of 2693 adult ewes slaughtered between 2017 and 2019 in two private slaughterhouses from Transylvania region (Romania) was evaluated. Lung tumors were subsequently assessed by cytology, histology, immunocytochemistry, immunohistochemistry, electron microscopy and DNA testing. RESULTS: Out of 2693 examined sheep, 34 had OPA (1.26% prevalence). The diaphragmatic lobes were the most affected. Grossly, the classical OPA was identified in 88.24% of investigated cases and the atypical OPA in 11.76% that included solitary myxomatous nodules. Histopathology results confirmed the presence of OPA in all suspected cases, which were classified into acinar and papillary types. Myxoid growths (MGs) were diagnosed in 6 classical OPA cases and in 2 cases of atypical form. Lung adenocarcinoma was positive for MCK and TTF-1, and MGs showed immunoreaction for Vimentin, Desmin and SMA; Ki67 expression of classical OPA was higher than atypical OPA and MGs. JSRV-MA was identified by IHC (94.11%) in both epithelial and mesenchymal cells of OPA. Immunocytochemistry and electron microscopy also confirmed the JSRV within the neoplastic cells. ExJSRV was identified by PCR in 97.05% of analyzed samples. Phylogenetic analysis revealed the presence of the exJSRV type 2 (MT809678.1) in Romanian sheep affected by lung cancer and showed a high similarity with the UK strain (AF105220.1). CONCLUSIONS: In this study, we confirmed for the first time in Romania the presence of exJSRV in naturally occurring OPA in sheep. Additionally, we described the first report of atypical OPA in Romania, and to the best of our knowledge, in Eastern Europe. Finally, we showed that MGs have a myofibroblastic origin.

A genome-wide perspective about the diversity and demographic history of seven Spanish goat breeds.

BACKGROUND: The main goal of the current work was to infer the demographic history of seven Spanish goat breeds (Malagueña, Murciano-Granadina, Florida, Palmera, Mallorquina, Bermeya and Blanca de Rasquera) based on genome-wide diversity data generated with the Illumina Goat SNP50 BeadChip (population size, N = 176). Five additional populations from Europe (Saanen and Carpathian) and Africa (Tunisian, Djallonké and Sahel) were also included in this analysis (N = 80) for comparative purposes. RESULTS: Our results show that the genetic background of Spanish goats traces back mainly to European breeds although signs of North African admixture were detected in two Andalusian breeds (Malagueña and Murciano-Granadina). In general, observed and expected heterozygosities were quite similar across the seven Spanish goat breeds under analysis irrespective of their population size and conservation status. For the Mallorquina and Blanca de Rasquera breeds, which have suffered strong population declines during the past decades, we observed increased frequencies of large-sized (ROH), a finding that is consistent with recent inbreeding. In contrast, a substantial part of the genome of the Palmera goat breed comprised short ROH, which suggests a strong and ancient founder effect. CONCLUSIONS: Admixture with African goats, genetic drift and inbreeding have had different effects across the seven Spanish goat breeds analysed in the current work. This has generated distinct patterns of genome-wide diversity that provide new clues about the demographic history of these populations.

Identification of an intronic regulatory mutation at the buffalo αS1-casein gene that triggers the skipping of exon 6.

The characterization of casein polymorphism is an essential step in order to understand the genetic basis of milk quality in dairy ruminants. In this work, we report the identification of a regulatory mutation at the buffalo αs1-casein (CSN1S1) gene that alters the normal processing of the primary transcript. Sequencing of CSN1S1 cDNA from individuals harbouring this new variant revealed that its most distinctive feature is the loss of exon 6 that encodes eight amino acids between positions 35-42 of mature protein. In an effort to map the causal mutation, we sequenced a genomic region spanning exons 5-7 of the buffalo CSN1S1 gene. This experiment allowed us to establish that exon 6-skipping is produced by a G to C substitution at the first position of intron 6 that inactivates the donor splice site. This mutation can be typed by PCR-RFLP by using either TaaI or Bpu10I diagnostic restriction enzymes, and it has a frequency of 0.18 in Romanian buffaloes. This exon skipping phenomenon is the first one described in buffalo CSN1S1 locus.

Pharmacological Effects of Methotrexate and Infliximab in a Rats Model of Diet-Induced Dyslipidemia and Beta-3 Overexpression on Endothelial Cells.

BACKGROUND: Hyperlipidemia and inflammation are critical components in the pathophysiology of endothelial disorder, which can lead to vascular complications. Our study aimed to evaluate the effects of immunomodulatory therapy (methotrexate and infliximab) in a diet-induced hyperlipidemia rat model. METHODS: Sprague-Dawley (wild type (WT), male, n = 32) rats were divided into four groups: one group fed with standard diet (SD), one group fed with high lipid diet (HLD), and two groups that received HLD and drug treatment (methotrexate (Mtx) or infliximab (Ifx)). In order to evaluate if modifications to the endothelial cells may influence the risk of vascular complications following hyperlipidemia or treatment reactivity, each group was doubled by a rats group that overexpressed beta-3 receptors on the endothelial cells (transgenic (TG-beta 3), male, n = 32). Serum lipid profile, liver enzymes, oxidative stress, and inflammation markers were determined. Histopathologic analysis of the liver and aorta was performed. RESULTS: After 9 weeks of HLD, rats exhibited significant pathologic serum lipid profiles, elevated oxidative stress, and pro-inflammatory markers. Additionally, the aortic histopathological analysis revealed aorta media-intima thickening (p < 0.05) in the transgenic group. Methotrexate and infliximab significantly decreased inflammation and oxidative stress parameters, but presented opposing effects on lipid profiles (methotrexate decreased, whereas infliximab increased the atherosclerosis index). Drug treatment decreased the aorta media-intima thickness (p < 0.05) only in transgenic rats. CONCLUSIONS: HLD was associated with hyperlipidemia, inflammation and oxidative stress. The overexpression of beta-3 receptors on endothelial cells increased aortic thickening in response to the HLD. Methotrexate and infliximab reduced oxidative stress and inflammation in all groups, but led to favorable histopathologic vascular results only in the transgenic groups.

Author Correction: Hotspots of recent hybridization between pigs and wild boars in Europe.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Geographical origin and species differentiation of Transylvanian cheese. Comparative study of isotopic and elemental profiling vs. DNA results.

The present study proposed the cheese differentiation, according to geographical production area and with respect to species (cow, sheep) of two traditional cheese specialties, (salty and ripened), produced in Transylvania, Romania. For this purpose, the elemental profile and carbon isotopic ratios (13C/12C) of cheese and extracted casein were corroborated through statistic supervised techniques to get the best discrimination markers. The manganese content, along with Rare Earth Elements (REE) concentrations, proved to be very powerful predictors, for the traditional salted cheese mainly, due to the direct influence of the local salted water. Despite that proposed techniques are not acknowledged methods for species differentiation, this approach allowed a successful discrimination of the animal species that produced the raw materials for cheese manufacturing (milk). The results generated by the developed chemometric model, for species differentiation, were compared with those obtained using Isoelectric focusing (IEF) and DNA tests. The proposed association of isotopic and elemental markers allowed a differentiation better than 92% for geographical provenance, of each investigated cheese specialties while, for species discrimination (cow vs. sheep) a percentage of 100% was obtained.

Hotspots of recent hybridization between pigs and wild boars in Europe.

After a strong demographic decline before World War II, wild boar populations are expanding and the species is now the second-most abundant ungulate in Europe. This increase raises concerns due to wild boar impact on crops and natural ecosystems and as potential vector of diseases. Additionally, wild boar can hybridize with domestic pigs, which could increase health risks and alter wild boar adaptive potential. We analysed 47,148 Single Nucleotide Polymorphisms in wild boar from Europe (292) and the Near East (16), and commercial (44) and local (255) pig breeds, to discern patterns of hybridization across Europe. We identified 33 wild boars with more than 10% domestic ancestry in their genome, mostly concentrated in Austria, Bosnia and Herzegovina, Bulgaria and Serbia. This difference is probably due to contrasting practices, with free-ranging vs. industrial farming but more samples would be needed to investigate larger geographic patterns. Our results suggest hybridization has occurred over a long period and is still ongoing, as we observed recent hybrids. Although wild and domestic populations have maintained their genetic distinctiveness, potential health threats raise concerns and require implementation of management actions and farming practices aimed at reducing contact between wild and domestic pigs.

Distribution and Functionality of Copy Number Variation across European Cattle Populations.

Copy number variation (CNV), which is characterized by large-scale losses or gains of DNA fragments, contributes significantly to genetic and phenotypic variation. Assessing CNV across different European cattle populations might reveal genetic changes responsible for phenotypic differences, which have accumulated throughout the domestication history of cattle as consequences of evolutionary forces that act upon them. To explore pattern of CNVs across European cattle, we genotyped 149 individuals, that represent different European regions, using the Illumina Bovine HD Genotyping array. A total of 9,944 autosomal CNVs were identified in 149 samples using a Hidden Markov Model (HMM) as employed in PennCNV. Animals originating from several breeds of British Isles, and Balkan and Italian regions, on average, displayed higher abundance of CNV counts than Dutch or Alpine animals. A total of 923 CNV regions (CNVRs) were identified by aggregating CNVs overlapping in at least two animals. The hierarchical clustering of CNVRs indicated low differentiation and sharing of high-frequency CNVRs between European cattle populations. Various CNVRs identified in the present study overlapped with olfactory receptor genes and genes related to immune system. In addition, we also detected a CNV overlapping the Kit gene in English longhorn cattle which has previously been associated with color-sidedness. To conclude, we provide a comprehensive overview of CNV distribution in genome of European cattle. Our results indicate an important role of purifying selection and genomic drift in shaping CNV diversity that exists between different European cattle populations.

Differential distribution of Y-chromosome haplotypes in Swiss and Southern European goat breeds.

The analysis of Y-chromosome variation has provided valuable clues about the paternal history of domestic animal populations. The main goal of the current work was to characterize Y-chromosome diversity in 31 goat populations from Central Eastern (Switzerland and Romania) and Southern Europe (Spain and Italy) as well as in reference populations from Africa and the Near East. Towards this end, we have genotyped seven single nucleotide polymorphisms (SNPs), mapping to the SRY, ZFY, AMELY and DDX3Y Y-linked loci, in 275 bucks from 31 populations. We have observed a low level of variability in the goat Y-chromosome, with just five haplotypes segregating in the whole set of populations. We have also found that Swiss bucks carry exclusively Y1 haplotypes (Y1A: 24%, Y1B1: 15%, Y1B2: 43% and Y1C: 18%), while in Italian and Spanish bucks Y2A is the most abundant haplotype (77%). Interestingly, in Carpathian goats from Romania the Y2A haplotype is also frequent (42%). The high Y-chromosome differentiation between Swiss and Italian/Spanish breeds might be due to the post-domestication spread of two different Near Eastern genetic stocks through the Danubian and Mediterranean corridors. Historical gene flow between Southern European and Northern African goats might have also contributed to generate such pattern of genetic differentiation.

A signature protein-based method to distinguish Mediterranean water buffalo and foreign breed milk.

A novel genetic variant at the αs1-casein locus of water buffalo (WB), 8-residue shorter than its wild-type has been found and sequenced. The internal deletion of the peptide E(35)KVNELsT(42) was confirmed by the isolation of the junction peptide. The 8-residue deletion mutant has a molecular weight that is 919 Da less than that of the wild-type. The novel isoform with a unique f35-42 deletion could be the result of the skipping of exon 6, generating an exon 6-deleted variant of αs1-casein. The wild-type and its shortened αs1-casein forms were found to co-exist in many individual milk samples. In contrast, the 8-residue, internally deleted αs1-casein variant did not occur in water buffaloes of the Mediterranean breed reared in Italy. Wild-type αs1-casein has 6 to 8 phosphate groups (P) while the internally deleted form 6 and 7P per molecule.

A high throughput genotyping approach reveals distinctive autosomal genetic signatures for European and Near Eastern wild boar.

The lack of a Near Eastern genetic signature in modern European porcine breeds indicates that, although domestic pigs from the Fertile Crescent entered Europe during the Neolithic, they were completely replaced by their European counterparts in a short window of time. Whilst the absence of such genetic signature has been convincingly demonstrated at the mitochondrial level, variation at the autosomal genomes of European and Near Eastern Sus scrofa has not been compared yet. Herewith, we have explored the genetic relationships among 43 wild boar from Europe (N = 21), Near East (N = 19) and Korea (N = 3), and 40 Iberian (N = 16), Canarian (N = 4) and Mangalitza (N = 20) pigs by using a high throughput SNP genotyping platform. After data filtering, 37,167 autosomal SNPs were used to perform population genetics analyses. A multidimensional scaling plot based on genome-wide identity-by-state pairwise distances inferred with PLINK showed that Near Eastern and European wild boar populations are genetically differentiated. Maximum likelihood trees built with TreeMix supported this conclusion i.e. an early population split between Near Eastern and European Sus scrofa was observed. Moreover, analysis of the data with Structure evidenced that the sampled Iberian, Canarian and Mangalitza pigs did not carry any autosomal signature compatible with a Near Eastern ancestry, a finding that agrees well with previous mitochondrial studies.