Integration of selective sweeps across the sheep genome: understanding the relationship between production and adaptation traits.

BACKGROUND: Livestock populations are under constant selective pressure for higher productivity levels for different selective purposes. This pressure results in the selection of animals with unique adaptive and production traits. The study of genomic regions associated with these unique characteristics has the potential to improve biological knowledge regarding the adaptive process and how it is connected to production levels and resilience, which is the ability of an animal to adapt to stress or an imbalance in homeostasis. Sheep is a species that has been subjected to several natural and artificial selective pressures during its history, resulting in a highly specialized species for production and adaptation to challenging environments. Here, the data from multiple studies that aim at mapping selective sweeps across the sheep genome associated with production and adaptation traits were integrated to identify confirmed selective sweeps (CSS). RESULTS: In total, 37 studies were used to identify 518 CSS across the sheep genome, which were classified as production (147 prodCSS) and adaptation (219 adapCSS) CSS based on the frequency of each type of associated study. The genes within the CSS were associated with relevant biological processes for adaptation and production. For example, for adapCSS, the associated genes were related to the control of seasonality, circadian rhythm, and thermoregulation. On the other hand, genes associated with prodCSS were related to the control of feeding behaviour, reproduction, and cellular differentiation. In addition, genes harbouring both prodCSS and adapCSS showed an interesting association with lipid metabolism, suggesting a potential role of this process in the regulation of pleiotropic effects between these classes of traits. CONCLUSIONS: The findings of this study contribute to a deeper understanding of the genetic link between productivity and adaptability in sheep breeds. This information may provide insights into the genetic mechanisms that underlie undesirable genetic correlations between these two groups of traits and pave the way for a better understanding of resilience as a positive ability to respond to environmental stressors, where the negative effects on production level are minimized.

A meta-analysis approach for annotation and identification of lncRNAs controlling perirenal fat deposition in suckling lambs.

Long non-coding RNAs (lncRNAs) are being studied in farm animals due to their association with traits of economic interest, such as fat deposition. Based on the analysis of perirenal fat transcriptomes, this research explored the relevance of these regulatory elements to fat deposition in suckling lambs. To that end, meta-analysis techniques have been implemented to efficiently characterize and detect differentially expressed transcripts from two different RNA-seq datasets, one including samples of two sheep breeds that differ in fat deposition features, Churra and Assaf (n = 14), and one generated from Assaf suckling lambs with different fat deposition levels (n = 8). The joint analysis of the 22 perirenal fat RNA-seq samples with the FEELnc software allowed the detection of 3953 novel lncRNAs. After the meta-analysis, 251 differentially expressed genes were identified, 21 of which were novel lncRNAs. Additionally, a co-expression analysis revealed that, in suckling lambs, lncRNAs may play a role in controlling angiogenesis and thermogenesis, processes highlighted in relation to high and low fat deposition levels, respectively. Overall, while providing information that could be applied for the improvement of suckling lamb carcass traits, this study offers insights into the biology of perirenal fat deposition regulation in mammals.

Effect of genotyping strategies on the sustained benefit of single-step genomic BLUP over multiple generations.

BACKGROUND: Single-step genomic best linear unbiased prediction (ssGBLUP) allows the inclusion of information from genotyped and ungenotyped individuals in a single analysis. This avoids the need to genotype all candidates with the potential benefit of reducing overall costs. The aim of this study was to assess the effect of genotyping strategies, the proportion of genotyped candidates and the genotyping criterion to rank candidates to be genotyped, when using ssGBLUP evaluation. A simulation study was carried out assuming selection over several discrete generations where a proportion of the candidates were genotyped and evaluation was done using ssGBLUP. The scenarios compared were: (i) three genotyping strategies defined by their protocol for choosing candidates to be genotyped (RANDOM: candidates were chosen at random; TOP: candidates with the best genotyping criterion were genotyped; and EXTREME: candidates with the best and worse criterion were genotyped); (ii) eight proportions of genotyped candidates (p); and (iii) two genotyping criteria to rank candidates to be genotyped (candidates' own phenotype or estimated breeding values). The criteria of the comparison were the cumulated gain and reliability of the genomic estimated breeding values (GEBV). RESULTS: The genotyping strategy with the greatest cumulated gain was TOP followed by RANDOM, with EXTREME behaving as RANDOM at low p and as TOP with high p. However, the reliability of GEBV was higher with RANDOM than with TOP. This disparity between the trend of the gain and the reliability is due to the TOP scheme genotyping the candidates with the greater chances of being selected. The extra gain obtained with TOP increases when the accuracy of the selection criterion to rank candidates to be genotyped increases. CONCLUSIONS: The best strategy to maximise genetic gain when only a proportion of the candidates are to be genotyped is TOP, since it prioritises the genotyping of candidates which are more likely to be selected. However, the strategy with the greatest GEBV reliability does not achieve the largest gain, thus reliability cannot be considered as an absolute and sufficient criterion for determining the scheme which maximises genetic gain.

Differences within Churra breed sheep in the early immune response to the infection by Teladorsagia circumcincta.

This study describes early immunological mechanisms that underlie resistance to Teladorsagia circumcincta infection in adult Churra sheep. After a first experimental infection, 6 animals were classified as resistant (RG) and 6 as susceptible (SG) to T. circumcincta infection based on their cumulative faecal egg count (cFEC) at the end of the infection. RG showed higher IgA levels against somatic antigen of T. circumcincta fourth-larvae stage (L4) in serum at day 3 post-infection (pi) (p < 0.05) and close to significance at day 21 pi (p = 0.06). Moreover, a strong negative correlation between cFEC and specific IgA was only significant in RG at day 3 pi (r = - 0.870; p < 0.05), but absent in SG. At the end of this infection, sheep were treated with moxidectin and infected again 3 weeks later to be slaughtered at day 7 pi. At necropsy, the specific IgA levels in gastric mucosa were similar between groups; the absence differences at day 7 pi could be due to a previous increase in the IgA response, probably around day 3 pi, as described during the first infection. L4 burden, 68% lower in RG than in SG, was influenced by the specific IgA in gastric mucus and the number of γδ T cells. RG group showed a positive correlation between γδ T cells and eosinophils (r = 0.900; p = 0.037); however, this correlation was not found in SG. These results show that these two phenotypes show different early immune response pattern to T. circumcincta infection in Churra sheep.

Microbiota characterization of sheep milk and its association with somatic cell count using 16s rRNA gene sequencing.

This work aimed to use 16S ribosomal RNA sequencing with the Illumina MiSeq platform to describe the milk microbiota from 50 healthy Assaf ewes. The global observed microbial community for clinically healthy milk samples analysed was complex and showed a vast diversity. The core microbiota of the sheep milk includes five genera: Staphylococcus, Lactobacillus, Corynebacterium, Streptococcus and Escherichia/Shigella. Although there are some differences, some of these genera are common with the microbiota core pattern of milk from other species, especially with dairy cows. The microbial composition of the studied samples, based on the definition of amplicon sequence variants, was analysed through a correlation network. A preliminary analysis by grouping the milk samples based on their somatic cell count (SCC), which is considered an indicator of subclinical mastitis (SM), showed certain differences for the core of the samples identified as SM. The differences in the microbiota diversity pattern among samples might also suggest that subclinical mastitis would be associated with the significant increase in some genera that are inhabitants of the mammary gland and a remarkable concomitant reduction in the microbial diversity. Additionally, we have also presented here a preliminary analysis to assess the impact of the sheep milk microbiome on SCC, as an indicator of subclinical mastitis. The results here reported provide a first characterization of the sheep milk microbiota and settle the basis for future studies in this field.

Exploring the mechanisms of resistance to Teladorsagia circumcincta infection in sheep through transcriptome analysis of abomasal mucosa and abomasal lymph nodes.

The present study exploited the RNA-seq technology to analyze the transcriptome of target tissues affected by the Teladorsagia circumcincta infection in two groups of adult ewes showing different statuses against gastrointestinal nematode (GIN) infection with the aim of identifying genes linked to GIN infection resistance in sheep. For this, based on the accumulated faecal egg count of 18 adult Churra ewes subjected to a first experimental infection with T. circumcincta, six ewes were classified as resistant and six others as susceptible to the infection. These 12 animals were dewormed and infected again. After humanitarian sacrifice of these 12 animals at day 7 post-infection, RNA samples were obtained from abomasal mucosa and lymph node tissues and RNA-Seq datasets were generated using an Illumina HiSeq 2000 sequencer. The distribution of the genes based on their expression level were very similar among the two different tissues and conditions. The differential expression analysis performed with two software (DESeq and EdgeR) only identified common differentially expressed genes (DEGs), a total of 106, in the lymph node samples which were considered as GIN-activated. The enrichment analysis performed for these GIN-activated genes identified some pathways related to cytokine-mediated immune response and the PPARG signaling pathway as well as disease terms related to inflammation and gastro-intestinal diseases as enriched. A systematic comparison with the results of previous studies confirmed the involvement of genes such as ITLN2, CLAC1 and galectins, in the immune mechanism activated against T. circumcincta in resistant sheep.

Variant discovery in the sheep milk transcriptome using RNA sequencing.

BACKGROUND: The identification of genetic variation underlying desired phenotypes is one of the main challenges of current livestock genetic research. High-throughput transcriptome sequencing (RNA-Seq) offers new opportunities for the detection of transcriptome variants (SNPs and short indels) in different tissues and species. In this study, we used RNA-Seq on Milk Sheep Somatic Cells (MSCs) with the goal of characterizing the genetic variation within the coding regions of the milk transcriptome in Churra and Assaf sheep, two common dairy sheep breeds farmed in Spain. RESULTS: A total of 216,637 variants were detected in the MSCs transcriptome of the eight ewes analyzed. Among them, a total of 57,795 variants were detected in the regions harboring Quantitative Trait Loci (QTL) for milk yield, protein percentage and fat percentage, of which 21.44% were novel variants. Among the total variants detected, 561 (2.52%) and 1,649 (7.42%) were predicted to produce high or moderate impact changes in the corresponding transcriptional unit, respectively. In the functional enrichment analysis of the genes positioned within selected QTL regions harboring novel relevant functional variants (high and moderate impact), the KEGG pathway with the highest enrichment was "protein processing in endoplasmic reticulum". Additionally, a total of 504 and 1,063 variants were identified in the genes encoding principal milk proteins and molecules involved in the lipid metabolism, respectively. Of these variants, 20 mutations were found to have putative relevant effects on the encoded proteins. CONCLUSIONS: We present herein the first transcriptomic approach aimed at identifying genetic variants of the genes expressed in the lactating mammary gland of sheep. Through the transcriptome analysis of variability within regions harboring QTL for milk yield, protein percentage and fat percentage, we have found several pathways and genes that harbor mutations that could affect dairy production traits. Moreover, remarkable variants were also found in candidate genes coding for major milk proteins and proteins related to milk fat metabolism. Several of the SNPs found in this study could be included as suitable markers in genotyping platforms or custom SNP arrays to perform association analyses in commercial populations and apply genomic selection protocols in the dairy production industry.

High-resolution analysis of selection sweeps identified between fine-wool Merino and coarse-wool Churra sheep breeds.

BACKGROUND: With the aim of identifying selection signals in three Merino sheep lines that are highly specialized for fine wool production (Australian Industry Merino, Australian Merino and Australian Poll Merino) and considering that these lines have been subjected to selection not only for wool traits but also for growth and carcass traits and parasite resistance, we contrasted the OvineSNP50 BeadChip (50 K-chip) pooled genotypes of these Merino lines with the genotypes of a coarse-wool breed, phylogenetically related breed, Spanish Churra dairy sheep. Genome re-sequencing datasets of the two breeds were analyzed to further explore the genetic variation of the regions initially identified as putative selection signals. RESULTS: Based on the 50 K-chip genotypes, we used the overlapping selection signals (SS) identified by four selection sweep mapping analyses (that detect genetic differentiation, reduced heterozygosity and patterns of haplotype diversity) to define 18 convergence candidate regions (CCR), five associated with positive selection in Australian Merino and the remainder indicating positive selection in Churra. Subsequent analysis of whole-genome sequences from 15 Churra and 13 Merino samples identified 142,400 genetic variants (139,745 bi-allelic SNPs and 2655 indels) within the 18 defined CCR. Annotation of 1291 variants that were significantly associated with breed identity between Churra and Merino samples identified 257 intragenic variants that caused 296 functional annotation variants, 275 of which were located across 31 coding genes. Among these, four synonymous and four missense variants (NPR2_His847Arg, NCAPG_Ser585Phe, LCORL_Asp1214Glu and LCORL_Ile1441Leu) were included. CONCLUSIONS: Here, we report the mapping and genetic variation of 18 selection signatures that were identified between Australian Merino and Spanish Churra sheep breeds, which were validated by an additional contrast between Spanish Merino and Churra genotypes. Analysis of whole-genome sequencing datasets allowed us to identify divergent variants that may be viewed as candidates involved in the phenotypic differences for wool, growth and meat production/quality traits between the breeds analyzed. The four missense variants located in the NPR2, NCAPG and LCORL genes may be related to selection sweep regions previously identified and various QTL reported in sheep in relation to growth traits and carcass composition.

Detection and replication of QTL underlying resistance to gastrointestinal nematodes in adult sheep using the ovine 50K SNP array.

BACKGROUND: Persistence of gastrointestinal nematode (GIN) infection and the related control methods have major impacts on the sheep industry worldwide. Based on the information generated with the Illumina OvineSNP50 BeadChip (50 K chip), this study aims at confirming quantitative trait loci (QTL) that were previously identified by microsatellite-based genome scans and identifying new QTL and allelic variants that are associated with indicator traits of parasite resistance in adult sheep. We used a commercial half-sib population of 518 Spanish Churra ewes with available data for fecal egg counts (FEC) and serum levels of immunoglobulin A (IgA) to perform different genome scan QTL mapping analyses based on classical linkage analysis (LA), a combined linkage disequilibrium and linkage analysis (LDLA) and a genome-wide association study (GWAS). RESULTS: For the FEC and IgA traits, we detected a total of three 5 % chromosome-wise significant QTL by LA and 63 significant regions by LDLA, of which 13 reached the 5 % genome-wise significance level. The GWAS also revealed 10 significant SNPs associated with IgAt, although no significant associations were found for LFEC. Some of the significant QTL for LFEC that were detected by LA and LDLA on OAR6 overlapped with a highly significant QTL that was previously detected in a different half-sib population of Churra sheep. In addition, several new QTL and SNP associations were identified, some of which show correspondence with effects that were reported for different populations of young sheep. Other significant associations that did not coincide with previously reported associations could be related to the specific immune response of adult animals. DISCUSSION: Our results replicate a FEC-related QTL located on OAR6 that was previously reported in Churra sheep and provide support for future research on the identification of the allelic variant that underlies this QTL. The small proportion of genetic variance explained by the detected QTL and the large number of functional candidate genes identified here are consistent with the hypothesis that GIN resistance/susceptibility is a complex trait that is not determined by individual genes acting alone but rather by complex multi-gene interactions. Future studies that combine genomic variation analysis and functional genomic information may help elucidate the biology of GIN disease resistance in sheep.

Major Histocompatibility Complex class IIB polymorphism in an ancient Spanish breed.

Genes from the Major Histocompatibility Complex class II region are involved in the presentation of antigens. Therefore, they have the key role in regulating the immune response and in the resistance to infections. We investigated the Major Histocompatibility Complex class IIB genes, DRB and DQB, in Churra sheep, one of the most important indigenous breeds of Spain. These genes are among the most polymorphic in the mammalian genome. Furthermore, often different numbers of class IIB genes per haplotype exist, complicating the genotyping and sequencing of these genes. Especially the DQB region is only partially characterized in sheep and the repertoire of DRB and DQB alleles in Churra sheep, an ancient breed, is unknown. Here, we sequenced the class IIB genes for 15 rams that are the pedigree heads of a selection Nucleus herd. In total, we found 12 DRB and 25 DQB alleles. From these, 3 and 15 were new, respectively. Fourteen haplotypes carrying one or two DQB alleles could be deduced and the evolutionary relationship of these was investigated by phylogenetic trees. Based on the sequences of these most common class II alleles, a more efficient genotyping system for larger numbers of Churra sheep will be developed.

Feed efficiency in dairy sheep: An insight from the milk transcriptome.

Introduction: As higher feed efficiency in dairy ruminants means a higher capability to transform feed nutrients into milk and milk components, differences in feed efficiency are expected to be partly linked to changes in the physiology of the mammary glands. Therefore, this study aimed to determine the biological functions and key regulatory genes associated with feed efficiency in dairy sheep using the milk somatic cell transcriptome. Material and methods: RNA-Seq data from high (H-FE, n = 8) and low (L-FE, n = 8) feed efficiency ewes were compared through differential expression analysis (DEA) and sparse Partial Least Square-Discriminant analysis (sPLS-DA). Results: In the DEA, 79 genes were identified as differentially expressed between both conditions, while the sPLS-DA identified 261 predictive genes [variable importance in projection (VIP) > 2] that discriminated H-FE and L-FE sheep. Discussion: The DEA between sheep with divergent feed efficiency allowed the identification of genes associated with the immune system and stress in L-FE animals. In addition, the sPLS-DA approach revealed the importance of genes involved in cell division (e.g., KIF4A and PRC1) and cellular lipid metabolic process (e.g., LPL, SCD, GPAM, and ACOX3) for the H-FE sheep in the lactating mammary gland transcriptome. A set of discriminant genes, commonly identified by the two statistical approaches, was also detected, including some involved in cell proliferation (e.g., SESN2, KIF20A, or TOP2A) or encoding heat-shock proteins (HSPB1). These results provide novel insights into the biological basis of feed efficiency in dairy sheep, highlighting the informative potential of the mammary gland transcriptome as a target tissue and revealing the usefulness of combining univariate and multivariate analysis approaches to elucidate the molecular mechanisms controlling complex traits.

Intergenerational impact of dietary protein restriction in dairy ewes on epigenetic marks in the perirenal fat of their suckling lambs.

In sheep, nutrition during the prepubertal stage is essential for growth performance and mammary gland development. However, the potential effects of nutrient restriction in a prepuberal stage over the progeny still need to be better understood. Here, the intergenerational effect of maternal protein restriction at prepubertal age (2 months of age) on methylation patterns was evaluated in the perirenal fat of Assaf suckling lambs. In total, 17 lambs from ewes subjected to dietary protein restriction (NPR group, 44% less protein) and 17 lambs from control ewes (C group) were analyzed. These lambs were ranked based on their carcass proportion of perirenal and cavitary fat and classified into HighPCF and LowPCF groups. The perirenal tissue from 4 NPR-LowPCF, 4 NPR-HighPCF, 4 C-LowPCF, and 4 C-HighPCF lambs was subjected to whole-genome bisulfite sequencing and differentially methylated regions (DMRs) were identified. Among other relevant processes, these DMRs were mapped in genes responsible for regulating the transition of brown to white adipose tissue and nonshivering thermoregulation, which might be associated with better adaptation/survival of lambs in the perinatal stage. The current study provides important biological insights about the intergenerational effect on the methylation pattern of an NPR in replacement ewes.

Influence of a temporary restriction of dietary protein in prepubertal ewe lambs on first lactation milk traits and response to a mammary gland inflammatory challenge.

This study evaluated the influence of a temporary nutritional protein restriction (NPR) performed, under commercial conditions, in prepubertal female lambs on first lactation milk production traits and the inflammatory response triggered by an inflammatory challenge of the. From 40 Assaf female lambs, we defined a control group (Cn = 20), which received a standard diet for replacement lambs and the NPR group (n = 20), which received the same diet but without soybean meal between 3 and 5 months of age. About 150 days after lambing, 24 of these ewes (13 NPR, 11C) were subjected to an intramammary infusion of E. coli lipopolysaccharide (LPS). Our dynamic study identified indicator traits of local (SCC) and systemic (rectal Ta, IL-6, CXCL8, IL-10, IL-36RA, VEGF-A) response to the LPS challenge. The NPR did not show significant effects on milk production traits and did not affect the SCC and rectal Ta after the LPS challenge. However, the NPR had a significant influence on 8 of the 14 plasma biomarkers analysed, in all the cases with higher relative values in the C group. The effects observed on VEGF-A (involved in vasculogenesis during mammary gland development and vascular permeability) and IL-10 (a regulatory cytokine classically known by its anti-inflammatory action) are the most remarkable to explain the differences found between groups. Whereas further studies should be undertaken to confirm these results, our findings are of interest considering the current concern about the future world's demand for protein and the need for animal production systems to evolve toward sustainability.

Transcriptome analysis of perirenal fat from Spanish Assaf suckling lamb carcasses showing different levels of kidney knob and channel fat.

Introduction: Suckling lamb meat is highly appreciated in European Mediterranean countries because of its mild flavor and soft texture. In suckling lamb carcasses, perirenal and pelvic fat depots account for a large fraction of carcass fat accumulation, and their proportions are used as an indicator of carcass quality. Material and Methods: This study aimed to characterize the genetic mechanisms that regulate fat deposition in suckling lambs by evaluating the transcriptomic differences between Spanish Assaf lambs with significantly different proportions of kidney knob and channel fat (KKCF) depots in their carcasses (4 High-KKCF lambs vs. 4 Low-KKCF lambs). Results: The analyzed fat tissue showed overall dominant expression of white adipose tissue gene markers, although due to the young age of the animals (17-36 days), the expression of some brown adipose tissue gene markers (e.g., UCP1, CIDEA) was still identified. The transcriptomic comparison between the High-KKCF and Low-KKCF groups revealed a total of 80 differentially expressed genes (DEGs). The enrichment analysis of the 49 DEGs with increased expression levels in the Low-KKCF lambs identified significant terms linked to the biosynthesis of lipids and thermogenesis, which may be related to the higher expression of the UCP1 gene in this group. In contrast, the enrichment analysis of the 31 DEGs with increased expression in the High-KKCF lambs highlighted angiogenesis as a key biological process supported by the higher expression of some genes, such as VEGF-A and THBS1, which encode a major angiogenic factor and a large adhesive extracellular matrix glycoprotein, respectively. Discussion: The increased expression of sestrins, which are negative regulators of the mTOR complex, suggests that the preadipocyte differentiation stage is being inhibited in the High-KKCF group in favor of adipose tissue expansion, in which vasculogenesis is an essential process. All of these results suggest that the fat depots of the High-KKCF animals are in a later stage of development than those of the Low-KKCF lambs. Further genomic studies based on larger sample sizes and complementary analyses, such as the identification of polymorphisms in the DEGs, should be designed to confirm these results and achieve a deeper understanding of the genetic mechanisms underlying fat deposition in suckling lambs.

Linkage disequilibrium and inbreeding estimation in Spanish Churra sheep.

BACKGROUND: Genomic technologies, such as high-throughput genotyping based on SNP arrays, have great potential to decipher the genetic architecture of complex traits and provide background information concerning genome structure in domestic animals, including the extent of linkage disequilibrium (LD) and haplotype blocks. The objective of this study was to estimate LD, the population evolution (past effective population size) and the level of inbreeding in Spanish Churra sheep. RESULTS: A total of 43,784 SNPs distributed in the ovine autosomal genome was analyzed in 1,681 Churra ewes. LD was assessed by measuring r2 between all pairs of loci. For SNPs up to 10 kb apart, the average r2 was 0.329; for SNPs separated by 200-500 kb the average r2 was 0.061. When SNPs are separated by more than 50 Mbp, the average r2 is the same as between non-syntenic SNP pairs (0.003). The effective population size has decreased through time, faster from 1,000 to 100 years ago and slower since the selection scheme started (15-25 generations ago). In the last generation, four years ago, the effective population size was estimated to be 128 animals. Inbreeding coefficients, although differed depending on the estimation approaches, were generally low and showed the same trend, which indicates that since 2003, inbreeding has been slightly increasing in the studied resource population. CONCLUSIONS: The extent of LD in Churra sheep persists over much more limited distances than reported in dairy cattle and seems to be similar to other ovine populations. Churra sheep show a wide genetic base, with a long-term viable effective population size that has been slightly decreasing since selection scheme began in 1986. The genomic dataset analyzed provided useful information for identifying low-level inbreeding in the sample, whereas based on the parameters reported here, a higher marker density than that analyzed here will be needed to successfully conduct accurate mapping of genes underlying production traits and genomic selection prediction in this sheep breed. Although the Ovine Assembly development is still in a draft stage and future refinements will provide a more accurate physical map that will improve LD estimations, this work is a first step towards the understanding of the genetic architecture in sheep.

Elucidating genes and gene networks linked to individual susceptibility to milk fat depression in dairy goats.

Dietary supplementation with marine lipids modulates ruminant milk composition toward a healthier fatty acid profile for consumers, but it also causes milk fat depression (MFD). Because the dairy goat industry is mainly oriented toward cheese manufacturing, MFD can elicit economic losses. There is large individual variation in animal susceptibility with goats more (RESPO+) or less (RESPO-) responsive to diet-induced MFD. Thus, we used RNA-Seq to examine gene expression profiles in mammary cells to elucidate mechanisms underlying MFD in goats and individual variation in the extent of diet-induced MFD. Differentially expression analyses (DEA) and weighted gene co-expression network analysis (WGCNA) of RNA-Seq data were used to study milk somatic cell transcriptome changes in goats consuming a diet supplemented with marine lipids. There were 45 differentially expressed genes (DEGs) between control (no-MFD, before diet-induced MFD) and MFD, and 18 between RESPO+ and RESPO-. Biological processes and pathways such as "RNA transcription" and "Chromatin modifying enzymes" were downregulated in MFD compared with controls. Regarding susceptibility to diet-induced MFD, we identified the "Triglyceride Biosynthesis" pathway upregulated in RESPO- goats. The WGCNA approach identified 9 significant functional modules related to milk fat production and one module to the fat yield decrease in diet-induced MFD. The onset of MFD in dairy goats is influenced by the downregulation of SREBF1, other transcription factors and chromatin-modifying enzymes. A list of DEGs between RESPO+ and RESPO- goats (e.g., DBI and GPD1), and a co-related gene network linked to the decrease in milk fat (ABCD3, FABP3, and PLIN2) was uncovered. Results suggest that alterations in fatty acid transport may play an important role in determining individual variation. These candidate genes should be further investigated.

Integrated analyses of the methylome and transcriptome to unravel sex differences in the perirenal fat from suckling lambs.

In sheep, differences were observed regarding fat accumulation and fatty acid (FA) composition between males and females, which may impact the quality and organoleptic characteristics of the meat. The integration of different omics technologies is a relevant approach for investigating biological and genetic mechanisms associated with complex traits. Here, the perirenal tissue of six male and six female Assaf suckling lambs was evaluated using RNA sequencing and whole-genome bisulfite sequencing (WGBS). A multiomic discriminant analysis using multiblock (s)PLS-DA allowed the identification of 314 genes and 627 differentially methylated regions (within these genes), which perfectly discriminate between males and females. These candidate genes overlapped with previously reported QTLs for carcass fat volume and percentage of different FAs in milk and meat from sheep. Additionally, differentially coexpressed (DcoExp) modules of genes between males (nine) and females (three) were identified that harbour 22 of these selected genes. Interestingly, these DcoExp were significantly correlated with fat percentage in different deposits (renal, pelvic, subcutaneous and intramuscular) and were associated with relevant biological processes for adipogenesis, adipocyte differentiation, fat volume and FA composition. Consequently, these genes may potentially impact adiposity and meat quality traits in a sex-specific manner, such as juiciness, tenderness and flavour.

Microbial community in resistant and susceptible Churra sheep infected by Teladorsagia circumcincta.

Gastrointestinal nematodes (GIN) are a major threat to health and welfare in small ruminants worldwide. Teladorsagia circumcincta is a nematode that inhabits the abomasum of sheep, especially in temperate regions, causing important economic losses. Given that T. circumcincta and microbiome share the same niche, interactions between them and the host are expected. Although it is known that within a sheep breed there are animals that are more resistant than others to infection by GIN, it is not known if the microbiome influences the phenotype of these animals. Under this condition, 12 sheep were classified according to their cumulative faecal egg count (cFEC) at the end of a first experimental infection, 6 as resistant group (RG) and 6 as susceptible group (SG) to T. circumcincta infection. Then, all sheep were experimentally infected with 70,000 L3 of T. circumcincta and at day 7 days post-infection were euthanized. At necropsy, gastric mucosa and gastric content from abomasum were collected to extract bacterial DNA and sequence V3-V4 region from 16S rRNA gene using Ilumina technology. After bioanalysis performed, results showed that α-diversity and ß-diversity remained similar in both groups. However, resistant phenotype sheep showed a higher number of bacteria butyrate-fermenting species as Clostridium sensu stricto 1 (abundance in RG: 1.29% and in SG: 0.069%; p = 0.05), and Turicibacter (abundance in RG: 0.31% and in SG: 0.027%; p = 0.07) in gastric content but also Serratia spp in gastric mucosa (abundance in RG: 0.12% and in SG: 0.041%; p = 0.07). A trend towards a significant negative correlation between cFEC and Clostridium sensu stricto 1 abundance in gastric content was detected (r = - 0.537; p = 0.08). These data suggest that microbiome composition could be another factor associated with the development of the resistant phenotype modifying the interaction with the host and the in last instance affecting the individual risk of infection.

Information content in genome-wide scans: concordance between patterns of genetic differentiation and linkage mapping associations.

BACKGROUND: Scanning the genome with high density SNP markers has become a standard approach for identifying regions of the genome showing substantial between-population genetic differentiation, and thus evidence of diversifying selection. Such regions may contain genes of large phenotypic effect. However, few studies have attempted to address the power or efficacy of such an approach. RESULTS: In this study, the patterns of allele frequency differences between two cattle breeds based on the Bovine HapMap study were compared with statistical evidence for QTL based on a linkage mapping study of an experimental population formed by a cross between the same breeds. Concordance between the two datasets was seen for chromosomes carrying QTL with strong statistical support, such as BTA5 and BTA18, which carry genes associated with coat color. For these chromosomes, there was a correspondence between the strength of the QTL signal along the chromosome and the degree of genetic differentiation between breeds. However, such an association was not seen in a broader comparison that also included chromosomes carrying QTL with lower significance levels. In addition, other chromosomal regions with substantial QTL effects did not include markers showing extreme between-breed genetic differentiation. Furthermore, the overall consistency between the two studies was weak, with low genome-wide correlation between the statistical values obtained in the linkage mapping study and between-breed genetic differentiation from the HapMap study. CONCLUSIONS: These results suggest that genomic diversity scans are capable of detecting regions associated with qualitative traits but may be limited in their power to detect regions associated with quantitative phenotypic differences between populations, which may depend on the marker resolution of the study and the level of LD in the populations under investigation.

Accuracy of Imputation of Microsatellite Markers from a 50K SNP Chip in Spanish Assaf Sheep.

Transitioning from traditional to new genotyping technologies requires the development of bridging methodologies to avoid extra genotyping costs. This study aims to identify the optimum number of single nucleotide polymorphisms (SNPs) necessary to accurately impute microsatellite markers to develop a low-density SNP chip for parentage verification in the Assaf sheep breed. The accuracy of microsatellite marker imputation was assessed with three metrics: genotype concordance (C), genotype dosage (length r2), and allelic dosage (allelic r2), for all imputation scenarios tested (0.5-10 Mb microsatellite flanking SNP windows). The imputation accuracy for the three metrics analyzed for all haplotype lengths tested was higher than 0.90 (C), 0.80 (length r2), and 0.75 (allelic r2), indicating strong genotype concordance. The window with 2 Mb length provides the best accuracy for the imputation procedure and the design of an affordable low-density SNP chip for parentage testing. We additionally evaluated imputation performance under two null models, naive (imputing the most common allele) and random (imputing by randomly selecting the allele), which in comparison showed weak genotype concordances (0.41 and 0.15, respectively). Therefore, we describe a precise methodology in the present article to impute multiallelic microsatellite genotypes from a low-density SNP chip in sheep and solve the problem of parentage verification when different genotyping platforms have been used across generations.