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.

Unravelling transmission ratio distortion across the bovine genome: identification of candidate regions for reproduction defects.

BACKGROUND: Biological mechanisms affecting gametogenesis, embryo development and postnatal viability have the potential to alter Mendelian inheritance expectations resulting in observable transmission ratio distortion (TRD). Although the discovery of TRD cases have been around for a long time, the current widespread and growing use of DNA technologies in the livestock industry provides a valuable resource of large genomic data with parent-offspring genotyped trios, enabling the implementation of TRD approach. In this research, the objective is to investigate TRD using SNP-by-SNP and sliding windows approaches on 441,802 genotyped Holstein cattle and 132,991 (or 47,910 phased) autosomal SNPs. RESULTS: The TRD was characterized using allelic and genotypic parameterizations. Across the whole genome a total of 604 chromosomal regions showed strong significant TRD. Most (85%) of the regions presented an allelic TRD pattern with an under-representation (reduced viability) of carrier (heterozygous) offspring or with the complete or quasi-complete absence (lethality) for homozygous individuals. On the other hand, the remaining regions with genotypic TRD patterns exhibited the classical recessive inheritance or either an excess or deficiency of heterozygote offspring. Among them, the number of most relevant novel regions with strong allelic and recessive TRD patterns were 10 and 5, respectively. In addition, functional analyses revealed candidate genes regulating key biological processes associated with embryonic development and survival, DNA repair and meiotic processes, among others, providing additional biological evidence of TRD findings. CONCLUSIONS: Our results revealed the importance of implementing different TRD parameterizations to capture all types of distortions and to determine the corresponding inheritance pattern. Novel candidate genomic regions containing lethal alleles and genes with functional and biological consequences on fertility and pre- and post-natal viability were also identified, providing opportunities for improving breeding success in cattle.

Genome-wide identification and characterization of Fusarium circinatum-responsive lncRNAs in Pinus radiata.

BACKGROUND: One of the most promising strategies of Pine Pitch Canker (PPC) management is the use of reproductive plant material resistant to the disease. Understanding the complexity of plant transcriptome that underlies the defence to the causal agent Fusarium circinatum, would greatly facilitate the development of an accurate breeding program. Long non-coding RNAs (lncRNAs) are emerging as important transcriptional regulators under biotic stresses in plants. However, to date, characterization of lncRNAs in conifer trees has not been reported. In this study, transcriptomic identification of lncRNAs was carried out using strand-specific paired-end RNA sequencing, from Pinus radiata samples inoculated with F. circinatum at an early stage of infection. RESULTS: Overall, 13,312 lncRNAs were predicted through a bioinformatics approach, including long intergenic non-coding RNAs (92.3%), antisense lncRNAs (3.3%) and intronic lncRNAs (2.9%). Compared with protein-coding RNAs, pine lncRNAs are shorter, have lower expression, lower GC content and harbour fewer and shorter exons. A total of 164 differentially expressed (DE) lncRNAs were identified in response to F. circinatum infection in the inoculated versus mock-inoculated P. radiata seedlings. The predicted cis-regulated target genes of these pathogen-responsive lncRNAs were related to defence mechanisms such as kinase activity, phytohormone regulation, and cell wall reinforcement. Co-expression network analysis of DE lncRNAs, DE protein-coding RNAs and lncRNA target genes also indicated a potential network regulating pectinesterase activity and cell wall remodelling. CONCLUSIONS: This study presents the first comprehensive genome-wide analysis of P. radiata lncRNAs and provides the basis for future functional characterizations of lncRNAs in relation to pine defence responses against F. circinatum.

Unrevealing functional candidate genes for bovine fertility through RNA sequencing meta-analysis and regulatory elements networks of co-expressed genes and lncRNAs.

The high genetic heterogeneity and environmental effects of subfertility in livestock species make the elucidation of the genetic mechanisms associated with reproductive efficiency a difficult task. Network and co-expression network meta-analyses were applied alongside genetic variant calling and long non-coding RNA (lncRNA) characterization to identify functionally relevant target genes and regulatory subnetworks associated with fertility in dairy cattle. In total, 505 lncRNAs (441 previously annotated in the bovine reference genome ARS-UCD 1.2 and 64 novel lncRNAs) were identified. Seven differentially expressed genes between high-fertile (HF) and sub-fertile (SF) Holstein cows were identified in the network meta-analysis (CA5A, ENSBTAG00000051149, ENSBTAG00000003272, DEFB7, DIO2, TRPV3, and COL4A4). Additionally, seven functional candidate differentially co-expressed (DcoExp) modules with a differential regulatory pattern (|z-score|>2) were identified between HF and SF cows. The functional candidate genes and DcoExp modules identified were associated with fertility relevant processes such as the regulation of embryonic implantation and proliferation, interaction and molecule transfer between the fetus and the cow, and the immune system. These results help to better understand the genetic mechanisms associated with reproductive efficiency in dairy cattle through the identification of potential biomarkers and genetic variants associated with differentially expressed regulatory gene and lncRNAs regulatory element networks.

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.

Generalized severe junctional epidermolysis bullosa with congenital absence of skin in churra lambs.

BACKGROUND: Up to 0.5% of churra lambs from two genetically related flocks showed congenital skin lesions of variable severity, jeopardizing the life of the lambs in the most severe cases. HYPOTHESIS/OBJECTIVES: The primary objective of this study was to classify the type of congenital epithelial disease suffered by these animals, based on the description of the macroscopic skin defects, the histological and ultrastructural changes and the hereditary nature of the condition. ANIMALS: Thirty affected newborn lambs from two genetically related flocks were studied. Three additional lambs acquired from two other flocks, which had no grossly apparent skin lesions and had died of infectious diseases, were studied as unaffected control animals. METHODS: Histological and ultrastructural examinations of skin and oral mucosa samples were performed. Pedigree analyses were used to investigate genealogical relationships. RESULTS: Generalized severe junctional epidermolysis bullosa with congenital absence of skin was described in all lambs studied and an autosomal recessive mode of inheritance was identified. CONCLUSIONS AND CLINICAL IMPORTANCE: The pathological findings and mode of inheritance in these lambs are similar to an inherited epidermolysis bullosa subtype of humans, which has not been reported previously in veterinary medicine.

Hereditary lissencephaly and cerebellar hypoplasia in Churra lambs.

BACKGROUND: Lissencephaly is a rare developmental brain disorder in veterinary and human medicine associated with defects in neuronal migration leading to a characteristic marked reduction or absence of the convolutional pattern of the cerebral hemispheres. In many human cases the disease has a genetic basis. In sheep, brain malformations, mainly cerebellar hypoplasia and forms of hydrocephalus, are frequently due to in utero viral infections. Although breed-related malformations of the brain have been described in sheep, breed-related lissencephaly has not been previously recorded in a peer reviewed publication. RESULTS: Here we report neuropathological findings in 42 newborn lambs from a pure Churra breed flock, with clinical signs of weakness, inability to walk, difficulty in sucking and muscular rigidity observed immediately after birth. All the lambs showed near-total agyria with only a rudimentary formation of few sulci and gyri, and a severe cerebellar hypoplasia. On coronal section, the cerebral grey matter was markedly thicker than that of age-matched unaffected lambs and the ventricular system was moderately dilated. Histologically, the normal layers of the cerebral cortex were disorganized and, using an immunohistochemical technique against neurofilaments, three layers were identified instead of the six present in normal brains. The hippocampus was also markedly disorganised and the number and size of lobules were reduced in the cerebellum. Heterotopic neurons were present in different areas of the white matter. The remainder of the brain structures appeared normal. The pathological features reported are consistent with the type LCH-b (lissencephaly with cerebellar hypoplasia group b) defined in human medicine. No involvement of pestivirus or bluetongue virus was detected by immunohistochemistry. An analysis of pedigree data was consistent with a monogenic autosomal recessive pattern inheritance. CONCLUSIONS: The study describes the clinical and pathological findings of lissencephaly with cerebellar hypoplasia in Churra lambs for which an autosomal recessive inheritance was the most likely cause. Histopathological features observed in the cerebral cortex and hippocampus are consistent with a possible failure in neuronal migration during brain development. This report suggests that lissencephaly should be considered in the differential diagnosis of congenital neurological disease in newborn lambs showing weakness, inability to walk and difficulty sucking.

Identifying hepatic genes regulating the ovine response to gastrointestinal nematodes using RNA-Sequencing.

Gastrointestinal nematode (GIN) infections are considered the most important disease of grazing sheep and due to increasing anthelmintic resistance, chemical control alone is inadequate. Resistance to Gastrointestinal nematode infection is a heritable trait, and through natural selection many sheep breeds have higher resistance. Studying the transcriptome from GIN-exposed and GIN-unexposed sheep using RNA-Sequencing technology can provide measurements of transcript levels associated with the host response to Gastrointestinal nematode infection, and these transcripts may harbor genetic markers that can be used in selective breeding programs to enhance disease resistance. The objective of this study was to compare liver transcriptomes of sheep naturally exposed to Gastrointestinal nematode s, with either high or low parasite burdens, to GIN-unexposed control sheep in order to identify key regulator genes and biological processes associated with Gastrointestinal nematode infection. Differential gene expression analysis revealed no significant differentially expressed genes (DEG) between sheep with a high or low parasite burden (p-value ≤0.01; False Discovery Rate (FDR) ≤ 0.05; and Fold-Change (FC) of > ±2). However, when compared to the control group, low parasite burden sheep showed 146 differentially expressed genes (64 upregulated and 82 downregulated in the low parasite burden group relative to the control), and high parasite burden sheep showed 159 differentially expressed genes (57 upregulated and 102 downregulated in the low parasite burden group relative to the control) (p-value ≤0.01; FDR ≤0.05; and FC of > ±2). Among these two lists of significant differentially expressed genes, 86 differentially expressed genes (34 upregulated, 52 downregulated in the parasited group relative to the control) were found in common between the two parasite burden groups compared to the control (GIN-unexposed sheep). Functional analysis of these significant 86 differentially expressed genes found upregulated genes involved in immune response and downregulated genes involved in lipid metabolism. Results of this study offer insight into the liver transcriptome during natural Gastrointestinal nematode exposure that helps provide a better understanding of the key regulator genes involved in Gastrointestinal nematode infection in sheep.

Genetic diversity of two native sheep breeds by genome-wide analysis of single nucleotide polymorphisms.

Wallachian and Sumava sheep are autochthonous breeds that have undergone a significant bottleneck effect and subsequent restoration efforts. The first objective of this study was to evaluate the degree of genetic variability of both breeds and, therefore, the current management of the breeding. The second was to determine whether these two breeds still retain their genetic uniqueness in relation to each other and other breeds, despite regenerative interventions. Our data consisted of 48 individuals of Sumava and 37 individuals of Wallachian sheep. The comparison data contained 25 other breeds (primarily European) from the HapMap dataset generated by the International Sheep Genomics Consortium. When comparing all 27 breeds, the Czech breeds clustered with 15 other breeds and formed a single branch with them according to Nei's distances. At the same time, however, the clusters of both breeds were integral and easily distinguishable from the others when displayed with principal component analysis (PCA). Population substructure analysis did not show any common genetic ancestry of the Czech national breeds and breeds used for regeneration or, eventually, breeds whose ancestral population was used for regeneration. The average values of FST were higher in Wallachian sheep (FST = 0.14) than in Sumava sheep (FST = 0.08). The linkage disequilibrium (LD) extension per autosome was higher in Wallachian than in Sumava sheep. Consequently, the Ne estimates five generations ago were 68 for Sumava versus 34 for Wallachian sheep. Both native Czech breeds exhibit a wide range of inbreeding based on the excess of homozygosity (FHOM) among individuals, from -0.04 to 0.16 in Sumava and from -0.13 to 0.12 in Wallachian. Average inbreeding based on runs of homozygosity was 0.21 in Sumava and 0.27 in Wallachian. Most detected runs of homozygosity (ROH) were less than 5 Mb long for both breeds. ROH segments longer than 15 Mb were absent in Wallachian sheep. Concerning putative selection signatures, a total of 471 candidate genes in Wallachian sheep within 11 hotspots and 653 genes within 13 hotspots in Sumava sheep were identified. Czech breeds appear to be well differentiated from each other and other European breeds. Their genetic diversity is low, especially in the case of the Wallachian breed. Sumava is not so threatened by low diversity but has a larger share of the non-native gene pool.

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.

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.

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.

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.

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.

Deviations from Mendelian Inheritance on Bovine X-Chromosome Revealing Recombination, Sex-of-Offspring Effects and Fertility-Related Candidate Genes.

Transmission ratio distortion (TRD), or significant deviations from Mendelian inheritance, is a well-studied phenomenon on autosomal chromosomes, but has not yet received attention on sex chromosomes. TRD was analyzed on 3832 heterosomal single nucleotide polymorphisms (SNPs) and 400 pseudoautosomal SNPs spanning the length of the X-chromosome using 436,651 genotyped Holstein cattle. On the pseudoautosomal region, an opposite sire-TRD pattern between male and female offspring was identified for 149 SNPs. This finding revealed unique SNPs linked to a specific-sex (Y- or X-) chromosome and describes the accumulation of recombination events across the pseudoautosomal region. On the heterosomal region, 13 SNPs and 69 haplotype windows were identified with dam-TRD. Functional analyses for TRD regions highlighted relevant biological functions responsible to regulate spermatogenesis, development of Sertoli cells, homeostasis of endometrium tissue and embryonic development. This study uncovered the prevalence of different TRD patterns across both heterosomal and pseudoautosomal regions of the X-chromosome and revealed functional candidate genes for bovine reproduction.

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.

Comparison between methods for measuring fecal egg count and estimating genetic parameters for gastrointestinal parasite resistance traits in sheep.

Fecal egg count (FEC) is an indicative measurement for parasite infection in sheep. Different FEC methods may show inconsistent results. Not accounting for inconsistencies can be problematic when integrating measurements from different FEC methods for genetic evaluation. The objectives of this study were to evaluate the difference in means and variances between two fecal egg counting methods used in sheep-the Modified McMaster (LMMR) and the Triple Chamber McMaster (LTCM); to estimate variance components for the two FEC methods, treating them as two different traits; and to integrate FEC data from the two different methods and estimate genetic parameters for FEC and other gastrointestinal parasite resistance traits. Fecal samples were collected from a commercial Rideau-Arcott sheep farm in Ontario. Fecal egg counting was performed using both LMMR and the LTCM methods. Other parasite resistance trait records were collected from the same farm including eye score (FAMACHA), body condition score (BCS), and body weight (WT). The two FEC methods were highly genetically (0.94) and phenotypically (0.88) correlated. However, the mean and variance between the two FEC methods were significantly different (P < 0.0001). Therefore, re-scaling is required prior to integrating data from the different methods. For the multiple trait analysis, data from the two fecal egg counting methods were integrated (LFEC) by using records for the LMMR when available and replacing missing records with re-standardized LTCM records converted to the same mean and variance of LMMR. Heritability estimates were 0.12 ± 0.04, 0.07 ± 0.05, 0.17 ± 0.06, and 0.24 ± 0.07 for LFEC egg count, FAMACHA, BCS, and WT, respectively. The estimated genetic correlations between FEC and the other parasite resistance traits were low and not significant (P > 0.05) for FAMACHA (r = 0.24 ± 0.32) and WT (r = 0.22 ± 0.19), and essentially zero for BCS (r = -0.03 ± 0.25), suggesting little to no benefit of using such traits as indicators for LFEC.