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.

A genome-wide analysis of copy number variation in Murciano-Granadina goats.

BACKGROUND: In this work, our aim was to generate a map of the copy number variations (CNV) segregating in a population of Murciano-Granadina goats, the most important dairy breed in Spain, and to ascertain the main biological functions of the genes that map to copy number variable regions. RESULTS: Using a dataset that comprised 1036 Murciano-Granadina goats genotyped with the Goat SNP50 BeadChip, we were able to detect 4617 and 7750 autosomal CNV with the PennCNV and QuantiSNP software, respectively. By applying the EnsembleCNV algorithm, these CNV were assembled into 1461 CNV regions (CNVR), of which 486 (33.3% of the total CNVR count) were consistently called by PennCNV and QuantiSNP and used in subsequent analyses. In this set of 486 CNVR, we identified 78 gain, 353 loss and 55 gain/loss events. The total length of all the CNVR (95.69 Mb) represented 3.9% of the goat autosomal genome (2466.19 Mb), whereas their size ranged from 2.0 kb to 11.1 Mb, with an average size of 196.89 kb. Functional annotation of the genes that overlapped with the CNVR revealed an enrichment of pathways related with olfactory transduction (fold-enrichment = 2.33, q-value = 1.61 × 10-10), ABC transporters (fold-enrichment = 5.27, q-value = 4.27 × 10-04) and bile secretion (fold-enrichment = 3.90, q-value = 5.70 × 10-03). CONCLUSIONS: A previous study reported that the average number of CNVR per goat breed was ~ 20 (978 CNVR/50 breeds), which is much smaller than the number we found here (486 CNVR). We attribute this difference to the fact that the previous study included multiple caprine breeds that were represented by small to moderate numbers of individuals. Given the low frequencies of CNV (in our study, the average frequency of CNV is 1.44%), such a design would probably underestimate the levels of the diversity of CNV at the within-breed level. We also observed that functions related with sensory perception, metabolism and embryo development are overrepresented in the set of genes that overlapped with CNV, and that these loci often belong to large multigene families with tens, hundreds or thousands of paralogous members, a feature that could favor the occurrence of duplications or deletions by non-allelic homologous recombination.

A genome-wide association analysis for body, udder, and leg conformation traits recorded in Murciano-Granadina goats.

Morphological traits are of great importance to dairy goat production given their effect on phenotypes of economic interest. However, their underlying genomic architecture has not yet been extensively characterized. Herein, we aimed to identify genomic regions associated with body, udder, and leg conformation traits recorded in 825 Murciano-Granadina goats. We genotyped this resource population using the GoatSNP50 BeadChip (Illumina Inc., San Diego, CA) and performed genome-wide association analyses using the GEMMA software. We found 2 genome-wide significant associations between markers rs268273468 [Capra hircus (CHI) 16:69617700] and rs268249346 (CHI 28:18321523) and medial suspensory ligament. In contrast, we did not detect any genome-wide significant associations for body and leg traits. Moreover, we found 12, 19, and 7 chromosome-wide significant associations for udder, body, and leg traits, respectively. Comparison of our data with previous studies revealed a low level of positional concordance between regions associated with morphological traits. In addition to technical factors, this lack of concordance could be due to a substantial level of genetic heterogeneity among breeds or to the strong polygenic background of morphological traits, which makes it difficult to detect genetic factors that have small phenotypic effects.

Non-parametric association analysis of additive and dominance effects of casein complex SNPs on milk content and quality in Murciano-Granadina goats.

Goat milk casein proteins (αS1, αS2, ß and κ) are encoded by four loci (CSN1S1, CSN1S2, CSN2 and CSN3, respectively) clustered within 250 kb in chromosome 6. In this study, 159 Murciano-Granadina goats were genotyped for 48 SNPs within the entire casein region. Phenotypes on milk yield and components were obtained from 2,594 dairy registries. Additive and dominance effects on milk composition and quality were studied using non-parametric tests and principal component analysis to prevent SNPs multicollinearity. Two deletions in exon 4 (CSN1S1 and CSN3), one in exon 7 (CSN2) and one in exon 15 (CSN1S2) have been found at frequencies ranging from 0.12 to 0.50. Bonferroni-corrected significant SNP additive and dominance effects were found for milk yield, fat, protein, dry matter and lactose, and somatic cells. Exons 15 and 7 were significantly associated with milk yield and components except for lactose and somatic cells, while exon 4 was significantly associated with milk yield and components except for protein and dry matter. SNPs' associations with somatic cells were less frequent and weaker than those with milk yield and components. As caseins increase, somatic cells decrease, reducing milk enzymatic activity and consumption suitability. Hence, including molecular information in breeding schemes may promote production efficiency, as selecting against undesirable alleles could prevent the compromises derived from their dominance effects.

A genomic map of climate adaptation in Mediterranean cattle breeds.

Domestic species such as cattle (Bos taurus taurus and B. t. indicus) represent attractive biological models to characterize the genetic basis of short-term evolutionary response to climate pressure induced by their post-domestication history. Here, using newly generated dense SNP genotyping data, we assessed the structuring of genetic diversity of 21 autochtonous cattle breeds from the whole Mediterranean basin and performed genome-wide association analyses with covariables discriminating the different Mediterranean climate subtypes. This provided insights into both the demographic and adaptive histories of Mediterranean cattle. In particular, a detailed functional annotation of genes surrounding variants associated with climate variations highlighted several biological functions involved in Mediterranean climate adaptation such as thermotolerance, UV protection, pathogen resistance or metabolism with strong candidate genes identified (e.g., NDUFB3, FBN1, METTL3, LEF1, ANTXR2 and TCF7). Accordingly, our results suggest that main selective pressures affecting cattle in Mediterranean area may have been related to variation in heat and UV exposure, in food resources availability and in exposure to pathogens, such as anthrax bacteria (Bacillus anthracis). Furthermore, the observed contribution of the three main bovine ancestries (indicine, European and African taurine) in these different populations suggested that adaptation to local climate conditions may have either relied on standing genomic variation of taurine origin, or adaptive introgression from indicine origin, depending on the local breed origins. Taken together, our results highlight the genetic uniqueness of local Mediterranean cattle breeds and strongly support conservation of these populations.

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.

The Southwestern fringe of Europe as an important reservoir of caprine biodiversity.

BACKGROUND: Portugal and Spain, with six and 22 officially recognized caprine breeds, encompass 25 % of the European Union goat census. Many of these populations have suffered strong demographic declines because of competition with exotic breeds and the phasing-out of low income rural activities. In this study, we have investigated the consequences of these and other demographic processes on the genetic diversity, population structure and inbreeding levels of Iberian and Atlantic goats. METHODS: A sample of 975 individuals representing 25 officially recognized breeds from Portugal and Spain, two small populations not officially recognized (Formentera and Ajuí goats) and two ecotypes of the Tinerfeña and Blanca Celtibérica breeds were genotyped with a panel of 20 microsatellite markers. A wide array of population genetics methods was applied to make inferences about the genetic relationships and demography of these caprine populations. RESULTS: Genetic differentiation among Portuguese and Spanish breeds was weak but significant (FST = 0.07; P < 0.001), which is probably the consequence of their short splitting times and extensive gene flow due to transhumance. In contrast, Canarian goats were strongly differentiated because of prolonged geographic isolation. Most populations displayed considerable levels of diversity (mean He = 0.65). CONCLUSIONS: High diversity levels and weak population structures are distinctive features of Portuguese and Spanish breeds. In general, these local breeds have a reduced census, but are still important reservoirs of genetic diversity. These findings reinforce the need for the implementation of management and breeding programs based on genetic data in order to minimize inbreeding, maintain overall genetic and allelic diversities and breed identities, while at the same time taking into account the within-breed genetic structure.

Analysis of ddRAD-seq data provides new insights into the genomic structure and patterns of diversity in Italian donkey populations.

With more than 150 recognized breeds, donkeys assume relevant economic importance, especially in developing countries. Even if the estimated number of heads worldwide is 53M, this species received less attention than other livestock species. Italy has traditionally been considered one of the cradles of European donkey breeding, and despite a considerable loss of biodiversity, today still counts nine autochthonous populations. A total of 220 animals belonging to nine different populations were genotyped using the double-digest restriction site associated DNA (ddRAD) sequencing to investigate the pattern of diversity using a multi-technique approach. A total of 418,602,730 reads were generated and successfully demultiplexed to obtain a medium-density SNP genotypes panel with about 27K markers. The diversity indices showed moderate levels of variability. The genetic distances and relationships, largely agree with the breeding history of the donkey populations under investigation. The results highlighted the separation of populations based on their genetic origin or geographical proximity between breeding areas, showed low to moderate levels of admixture, and indicated a clear genetic difference in some cases. For some breeds, the results also validate the success of proper management conservation plans. Identified runs of homozygosity islands, mapped within genomic regions related to immune response and local adaptation, are consistent with the characteristics of the species known for its rusticity and adaptability. This study is the first exhaustive genome-wide analysis of the diversity of Italian donkey populations. The results emphasized the high informativeness of genome-wide markers retrieved through the ddRAD approach. The findings take on great significance in designing and implementing conservation strategies. Standardized genotype arrays for donkey species would make it possible to combine worldwide datasets to provide further insights into the evolution of the genomic structure and origin of this important genetic resource.

Donkeys assume relevant economic importance in several countries worldwide. However, the genetic structure of these populations is less investigated compared to other species. The aim of this study was to investigate the genetic background of nine different Italian donkey populations. A total of 220 animals were genotyped with about 27K markers extracted by the double-digest restriction site associated DNA sequencing. The consistency of the results across different approaches agreed with the demographic history, the origin, and previous results on the nine donkey populations, suggesting that our conclusions are robust. Moreover, the results of the present study highlighted low to moderate levels of admixture and, for some breeds, confirmed the success of proper management conservation plans.

Analysis of conservation priorities of Iberoamerican cattle based on autosomal microsatellite markers.

BACKGROUND: Determining the value of livestock breeds is essential to define conservation priorities, manage genetic diversity and allocate funds. Within- and between-breed genetic diversity need to be assessed to preserve the highest intra-specific variability. Information on genetic diversity and risk status is still lacking for many Creole cattle breeds from the Americas, despite their distinct evolutionary trajectories and adaptation to extreme environmental conditions. METHODS: A comprehensive genetic analysis of 67 Iberoamerican cattle breeds was carried out with 19 FAO-recommended microsatellites to assess conservation priorities. Contributions to global diversity were investigated using alternative methods, with different weights given to the within- and between-breed components of genetic diversity. Information on Iberoamerican plus 15 worldwide cattle breeds was used to investigate the contribution of geographical breed groups to global genetic diversity. RESULTS: Overall, Creole cattle breeds showed a high level of genetic diversity with the highest level found in breeds admixed with zebu cattle, which were clearly differentiated from all other breeds. Within-breed kinships revealed seven highly inbred Creole breeds for which measures are needed to avoid further genetic erosion. However, if contribution to heterozygosity was the only criterion considered, some of these breeds had the lowest priority for conservation decisions. The Weitzman approach prioritized highly differentiated breeds, such as Guabalá, Romosinuano, Cr. Patagonico, Siboney and Caracú, while kinship-based methods prioritized mainly zebu-related breeds. With the combined approaches, breed ranking depended on the weights given to the within- and between-breed components of diversity. Overall, the Creole groups of breeds were generally assigned a higher priority for conservation than the European groups of breeds. CONCLUSIONS: Conservation priorities differed significantly according to the weight given to within- and between-breed genetic diversity. Thus, when establishing conservation programs, it is necessary to also take into account other features. Creole cattle and local isolated breeds retain a high level of genetic diversity. The development of sustainable breeding and crossbreeding programs for Creole breeds, and the added value resulting from their products should be taken into consideration to ensure their long-term survival.

Genetic structure, relationships and admixture with wild relatives in native pig breeds from Iberia and its islands.

BACKGROUND: Native pig breeds in the Iberian Peninsula are broadly classified as belonging to either the Celtic or the Mediterranean breed groups, but there are other local populations that do not fit into any of these groups. Most of the native pig breeds in Iberia are in danger of extinction, and the assessment of their genetic diversity and population structure, relationships and possible admixture between breeds, and the appraisal of conservation alternatives are crucial to adopt appropriate management strategies. METHODS: A panel of 24 microsatellite markers was used to genotype 844 animals representing the 17 most important native swine breeds and wild populations existing in Portugal and Spain and various statistical tools were applied to analyze the results. RESULTS: Genetic diversity was high in the breeds studied, with an overall mean of 13.6 alleles per locus and an average expected heterozygosity of 0.80. Signs of genetic bottlenecks were observed in breeds with a small census size, and population substructure was present in some of the breeds with larger census sizes. Variability among breeds accounted for about 20% of the total genetic diversity, and was explained mostly by differences among the Celtic, Mediterranean and Basque breed groups, rather than by differences between domestic and wild pigs. Breeds clustered closely according to group, and proximity was detected between wild pigs and the Mediterranean cluster of breeds. Most breeds had their own structure and identity, with very little evidence of admixture, except for the Retinto and Entrepelado varieties of the Mediterranean group, which are very similar. Genetic influence of the identified breed clusters extends beyond the specific geographical areas across borders throughout the Iberian Peninsula, with a very sharp transition from one breed group to another. Analysis of conservation priorities confirms that the ranking of a breed for conservation depends on the emphasis placed on its contribution to the between- and within-breed components of genetic diversity. CONCLUSIONS: Native pig breeds in Iberia reveal high levels of genetic diversity, a solid breed structure and a clear organization in well-defined clusters.

Variations in Casein Genes Are Associated with Milk Protein and Fat Contents in Sarda Goats (Capra hircus), with an Important Role of CSN1S2 for Milk Yield.

This work aimed to assess the variability of casein genes in a population of 153 bucks and 825 lactating does of the Sarda breed, and to perform association analysis between polymorphic sites and milk yield and composition traits. To genotype the casein genes, we chose an SNP panel including 44 SNPs mapping to the four casein genes CSN1S1, CSN2, CSN1S2, and CSN3. Genotyping (made by KASP™ genotyping assay, based on competitive allele-specific PCR) revealed the high variability of the Sarda goat, and haplotype analysis revealed linkage disequilibrium (LD) between CSN1S1 and CSN2 genes, in addition to two LD blocks within the CSN1S2 and two LD blocks within the CSN3 gene, in bucks and does. Association analysis revealed that variability at all four casein genes was associated with milk protein content, total solids, and milk energy. The three Ca-sensitive casein genes were associated with lipid content, and CSN1S2 showed a unique pattern, with intron variants associated with milk yield, in addition to milk pH, NaCl, and SCS (Somatic Cell Score). This information might prove useful in selection schemes and in future investigations aiming to better understand the biology of lactation, and the direct link between genotype and phenotype.

Comprehensive Gene Expression Profiling Analysis of Adipose Tissue in Male Individuals from Fat- and Thin-Tailed Sheep Breeds.

It has been shown that tail fat content varies significantly among sheep breeds and plays a significant role in meat quality. Recently, significant efforts have been made to understand the physiological, biochemical, and genomic regulation of fat deposition in sheep tails in order to unravel the mechanisms underlying energy storage and adipose tissue lipid metabolism. RNA-seq has enabled us to provide a high-resolution snapshot of differential gene expression between fat- and thin-tailed sheep breeds. Therefore, three RNA-seq datasets were meta-analyzed for the current work to elucidate the transcriptome profile differences between them. Specifically, we identified hub genes, performed gene ontology (GO) analysis, carried out enrichment analyses of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and validated hub genes using machine learning algorithms. This approach revealed a total of 136 meta-genes, 39 of which were not significant in any of the individual studies, indicating the higher statistical power of the meta-analysis. Furthermore, the results derived from the use of machine learning revealed POSTN, K35, SETD4, USP29, ANKRD37, RTN2, PRG4, and LRRC4C as substantial genes that were assigned a higher weight (0.7) than other meta-genes. Among the decision tree models, the Random Forest ones surpassed the others in adipose tissue predictive power fat deposition in fat- and thin-tailed breeds (accuracy > 0.85%). In this regard, combining meta-analyses and machine learning approaches allowed for the identification of three important genes (POSTN, K35, SETD4) related to lipid metabolism, and our findings could help animal breeding strategies optimize fat-tailed breeds' tail sizes.

Genetic diversity in 10 populations of domestic Turkeys by using microsatellites markers.

The domestic turkey is a native breed in danger of extinction due to the introduction of new breeds specializing in meat production and yield. Turkeys have lost some prominence in urban areas, and only certain breeds of turkeys are preserved in rural areas. Wild and domestic turkeys are different; rural or indigenous turkeys, with black plumage, were domesticated from Mexican turkeys and have been reproduced throughout Latin America. Some of them were taken to Europe in the 16th century and later arrived in North America, where they crossed with another wild species, from which the bronze turkey emerged: the ancestor of all commercial turkeys. The objective of the present work was to evaluate the genetic diversity in 10 populations of domestic turkeys worldwide by using breeds from Europe: Spain and Italy; America: Mexico, United States and Brazil; and the Near East: Iran and Egypt. A total of 522 blood samples of both sexes were collected from domestic turkey populations. Thirty-four microsatellites were used to obtain genetic parameters, and genetic diversity was evaluated. All microsatellites used were polymorphic, and a total of 427 alleles were detected across the 34 markers investigated. In this study, a mean number of 13.44 alleles was found. The four most diverse breeds were from the Andalusia, Mexico, United States, and wild populations, which had the highest mean heterozygosity expected (0.619, 0.612, 0.650, and 0.773) and heterozygosity observed (0.422, 0.521, 0.429, and 0.627), respectively. The MNT348 marker deviated from the HWE in all populations. Our study has shown that the populations close to the species origin are more diverse than those resulting from posterior expansions. Mexican birds were the most diverse, followed by the Spanish populations because Spain imported a large number of turkeys coming from America. Such information can be complementary to other genotypic data required to validate the evolutionary relationships among turkey populations.

Seven Shades of Grey: A Follow-Up Study on the Molecular Basis of Coat Colour in Indicine Grey Cattle Using Genome-Wide SNP Data.

Shades of grey and brown are a dominant component in mammal coat colours, representing a fundamental trait involved in a great number of processes including cryptism, sexual selection and signalling. The genetic mechanisms of the grey colouration in mammals are very complex and controlled by hundreds of genes whose effects and interactions are still largely unclear. In this study, we adopted a robust multi-cohort Fst outlier approach based on pairwise contrasts between seven grey indicine cattle breeds and both taurine and indicine non-grey cattle breeds in order to find genomic regions potentially related to the grey colouration. On the basis of three main drawn settings, built in order to control both the effect of the sample size and the genetic structure, we have identified some signals common to those obtained in a previous work employing only taurine cattle. In particular, using the top 1% Fst approach, we detected a candidate region (22.6-23.8 megabases) on chromosome 14 in which genes related to pigmentation have been already documented. In addition, when we constructed a phylogenetic tree using the significant markers identified in this study and including also the genotyping data at these loci of both the grey taurine and the extinct wild auroch, we found a topological repartition consistent with breed colour pattern rather than with the known bovine evolutionary history. Thus, on the basis of this evidence, together with the geographical distribution of the current taurine grey cattle, an ancestral indicine origin for the grey phenotype would seem to be a conceivable interpretation. In this context, a higher thermo-tolerance and less UV-induced damage of the grey phenotype might have favoured the retention of advantageous genes into the taurine genome during the post-Neolithic human-mediated cattle expansions.

New Insight into the Genome-Wide Diversity and Admixture of Six Colombian Sheep Populations.

Creole sheep represent a strategic genetic resource for populations living in marginal areas under financial restrictions on the American continent. Six Colombian sheep breeds (two wool (BCL-Boyacá and NCL-Nariño, 12 and 14 samples) and four hair (OPCE-Ethiopian, 54 samples; OPCS-Sudan, 74 samples; OPCP-Pelibeuy, 59 samples; OPCW-Wayúu, 24 samples) were genotyped using the Illumina Ovine SNP50 BeadChip. Data was also included from international 44 breeds from International Sheep Genomics Consortium (ISGC) and from data published in previous a previous work on the Caribbean and African breeds. Although geographically separated, wool (NCL, BCL) and hair types (OPCE, OPCS, OPCW) presented little genetic differentiation (FST 0.05) at a global level but several groups of animals separated suggesting local clustering due to geographical isolation. The OPCP underwent a recent crossing with Mexican Pelibuey, explaining its differentiation. Findings in this work such as the proximity to West African Djallonké (WAD) and Barbados Black Belly (BBB), suggest different introductions of African type animals from the Caribbean region on a pre-existing genetic basis formed by animals deriving from the first importations coming from Europe in colonial times. As expected, Colombian wool breeds showed, in particular in Admixture software results, a greater genomic component in common with European breeds and in particular with Iberian ones (Churra). This study provides a basis for future research into the genetic diversity within and between the Colombian sheep breeds analysed, and scientific data for policy decisions on Farm Animal Genetic Resources (FAnGR).

Population Structure of Curraleiro Pé-Duro Cattle and its Relationship With the Serological Profile Against Pathogens of Economic and Zoonotic Interest.

Curraleiro Pé-Duro (CPD) cattle, a Brazilian local breed, are animals that are highly resistant to infectious and parasitic diseases. Strategies for the conservation of the breed and the genetic resistance to diseases should also consider the characteristics of the breed. The objective of this study was to analyze the diversity and population structure of the CPD breed using microsatellite markers and to correlate the serological profiles for causative agents of brucellosis, leptospirosis, neosporosis, leukosis, infectious bovine rhinotracheitis, and bovine viral diarrhea. DNA samples of 608 bovines were amplified and genotyped using 28 microsatellite markers for breed characterization. The genotypes were assigned to three clusters, indicating a substructure population related to the geographic distance. The observed heterozygosity was lower than that expected in most loci, and fixation index (Fst) in the subpopulation values ranged from 0.03073 (ETH3) to 0.402 (INRA63) on cluster 1, from 0.00 (INRA32) to 0.39359 (INRA63) on cluster 2, and from 0.00 (INRA32) to 0.28483 (TGLA53) on cluster 3. The Brucella spp., Neospora caninum, and Bovine Leukemia Virus frequencies were significantly different (p < 0.05) between clusters. The differences between the occurrences of diseases were not sufficient to indicate a genetic subpopulation with increased resistance to infections.

Genomic patterns of homozygosity and inbreeding depression in Murciano-Granadina goats.

BACKGROUND: Inbreeding depression can adversely affect traits related to fitness, reproduction and productive performance. Although current research suggests that inbreeding levels are generally low in most goat breeds, the impact of inbreeding depression on phenotypes of economic interest has only been investigated in a few studies based on genealogical data. RESULTS: We genotyped 1040 goats with the Goat SNP50 BeadChip. This information was used to estimate different molecular inbreeding coefficients and characterise runs of homozygosity and homozygosity patterns. We detected 38 genomic regions with increased homozygosity as well as 8 ROH hotspots mapping to chromosomes 1, 2, 4, 6, 14, 16 and 17. Eight hundred seventeen goats with available records for dairy traits were analysed to evaluate the potential consequences of inbreeding depression on milk phenotypes. Four regions on chromosomes 8 and 25 were significantly associated with inbreeding depression for the natural logarithm of the somatic cell count. Notably, these regions contain several genes related with immunity, such as SYK, IL27, CCL19 and CCL21. Moreover, one region on chromosome 2 was significantly associated with inbreeding depression for milk yield. CONCLUSIONS: Although genomic inbreeding levels are low in Murciano-Granadina goats, significant evidence of inbreeding depression for the logarithm of the somatic cell count, a phenotype closely associated with udder health and milk yield, have been detected in this population. Minimising inbreeding would be expected to augment economic gain by increasing milk yield and reducing the incidence of mastitis, which is one of the main causes of dairy goat culling.

Refining the Camelus dromedarius Myostatin Gene Polymorphism through Worldwide Whole-Genome Sequencing.

Myostatin (MSTN) is a highly conserved negative regulator of skeletal muscle in mammals. Inactivating mutations results in a hyper-muscularity phenotype known as "double muscling" in several livestock and model species. In Camelus dromedarius, the gene structure organization and the sequence polymorphisms have been previously investigated, using Sanger and Next-Generation Sequencing technologies on a limited number of animals. Here, we carried out a follow-up study with the aim to further expand our knowledge about the sequence polymorphisms at the myostatin locus, through the whole-genome sequencing data of 183 samples representative of the geographical distribution range for this species. We focused our polymorphism analysis on the ±5 kb upstream and downstream region of the MSTN gene. A total of 99 variants (77 Single Nucleotide Polymorphisms and 22 indels) were observed. These were mainly located in intergenic and intronic regions, with only six synonymous Single Nucleotide Polymorphisms in exons. A sequence comparative analysis among the three species within the Camelus genus confirmed the expected higher genetic distance of C. dromedarius from the wild and domestic two-humped camels compared to the genetic distance between C. bactrianus and C. ferus. In silico functional prediction highlighted: (i) 213 differential putative transcription factor-binding sites, out of which 41 relative to transcription factors, with known literature evidence supporting their involvement in muscle metabolism and/or muscle development; and (ii) a number of variants potentially disrupting the canonical MSTN splicing elements, out of which two are discussed here for their potential ability to generate a prematurely truncated (inactive) form of the protein. The distribution of the considered variants in the studied cohort is discussed in light of the peculiar evolutionary history of this species and the hypothesis that extremely high muscularity, associated with a homozygous condition for mutated (inactivating) alleles at the myostatin locus, may represent, in arid desert conditions, a clear metabolic disadvantage, emphasizing the thermoregulatory and water availability challenges typical of these habitats.

An inter-laboratory study of DNA-based identity, parentage and species testing in animal forensic genetics.

The probative value of animal forensic genetic evidence relies on laboratory accuracy and reliability. Inter-laboratory comparisons allow laboratories to evaluate their performance on specific tests and analyses and to continue to monitor their output. The International Society for Animal Genetics (ISAG) administered animal forensic comparison tests (AFCTs) in 2016 and 2018 to assess the limitations and capabilities of laboratories offering forensic identification, parentage and species determination services. The AFCTs revealed that analyses of low DNA template concentrations (≤300 pg/µL) constitute a significant challenge that has prevented many laboratories from reporting correct identification and parentage results. Moreover, a lack of familiarity with species testing protocols, interpretation guidelines and representative databases prevented over a quarter of the participating laboratories from submitting correct species determination results. Several laboratories showed improvement in their genotyping accuracy over time. However, the use of forensically validated standards, such as a standard forensic short tandem repeat (STR) kit, preferably with an allelic ladder, and stricter guidelines for STR typing, may have prevented some common issues from occurring, such as genotyping inaccuracies, missing data, elevated stutter products and loading errors. The AFCTs underscore the importance of conducting routine forensic comparison tests to allow laboratories to compare results from each other. Laboratories should keep improving their scientific and technical capabilities and continuously evaluate their personnel's proficiency in critical techniques such as low copy number (LCN) analysis and species testing. Although this is the first time that the ISAG has conducted comparison tests for forensic testing, findings from these AFCTs may serve as the foundation for continuous improvements of the overall quality of animal forensic genetic testing.

Evaluation of Different Test-Day Milk Recording Protocols by Wood's Model Application for the Estimation of Dairy Goat Milk and Milk Constituent Yield.

Goats have important social and economic roles in many countries because of their ability to survive and be productive in marginal areas. The overarching aim of this study was to compare the application of Wood's model to different test-day milk recording protocols for estimation of total milk, fat, and protein yield in dairy goats. A total of 465 goats were used (Garganica, 78; Girgentana, 81; Jonica, 76; Maltese, 77; Red Mediterranean, 76; Saanen, 77). Milk yield was recorded every 15 days throughout lactation of 210 days, for a total of 14 collection days, during both morning and afternoon milking sessions. Milk samples were collected and analyzed for protein and fat. The fat-corrected milk was standardized at 35g fat/kg of milk. Wood models showed high R2 values, and thus good fitting, in all the considered breeds. Wood model applied to first, second, fourth, and sixth month recordings (C) and ICAR estimation showed total milk yield very close to Wood's model applied to all 14 recordings (A) (p > 0.38). Differently, Wood's model applied to the first, second, third, and fourth month recording (B) estimation showed great differences (p < 0.01). This could be applied for farms that had the necessity to synchronize flock groups for kidding in order to produce kid meat. In farms that apply the estrus induction and/or synchronization for kidding, it would be possible to perform only four test-day milk recordings and to apply the Wood's model on them in order to obtain the estimation of total milk, fat, and protein yield during lactation for animals inscribed, or to be inscribed, to the genealogical book.