A dimensional reduction approach to modulate the core ruminal microbiome associated with methane emissions via selective breeding.

The rumen is a complex microbial system of substantial importance in terms of greenhouse gas emissions and feed efficiency. This study proposes combining metagenomic and host genomic data for selective breeding of the cow hologenome toward reduced methane emissions. We analyzed nanopore long reads from the rumen metagenome of 437 Holstein cows from 14 commercial herds in 4 northern regions in Spain. After filtering, data were treated as compositional. The large complexity of the rumen microbiota was aggregated, through principal component analysis (PCA), into few principal components (PC) that were used as proxies of the core metagenome. The PCA allowed us to condense the huge and fuzzy taxonomical and functional information from the metagenome into a few PC. Bivariate animal models were applied using these PC and methane production as phenotypes. The variability condensed in these PC is controlled by the cow genome, with heritability estimates for the first PC of ~0.30 at all taxonomic levels, with a large probability (>83%) of the posterior distribution being >0.20 and with the 95% highest posterior density interval (95%HPD) not containing zero. Most genetic correlation estimates between PC1 and methane were large (≥0.70), with most of the posterior distribution (>82%) being >0.50 and with its 95%HPD not containing zero. Enteric methane production was positively associated with relative abundance of eukaryotes (protozoa and fungi) through the first component of the PCA at phylum, class, order, family, and genus. Nanopore long reads allowed the characterization of the core rumen metagenome using whole-metagenome sequencing, and the purposed aggregated variables could be used in animal breeding programs to reduce methane emissions in future generations.

Genome-wide expression profiling in muscle and subcutaneous fat of lambs in response to the intake of concentrate supplemented with vitamin E.

BACKGROUND: The objective of this study was to acquire a broader, more comprehensive picture of the transcriptional changes in the L. Thoracis muscle (LT) and subcutaneous fat (SF) of lambs supplemented with vitamin E. Furthermore, we aimed to identify novel genes involved in the metabolism of vitamin E that might also be involved in meat quality. In the first treatment, seven lambs were fed a basal concentrate from weaning to slaughter (CON). In the second treatment, seven lambs received basal concentrate from weaning to 4.71 ± 2.62 days and thereafter concentrate supplemented with 500 mg dl-α-tocopheryl acetate/kg (VE) during the last 33.28 ± 1.07 days before slaughter. RESULTS: The addition of vitamin E to the diet increased the α-tocopherol muscle content and drastically diminished the lipid oxidation of meat. Gene expression profiles for treatments VE and CON were clearly separated from each other in the LT and SF. Vitamin E supplementation had a dramatic effect on subcutaneous fat gene expression, showing general up-regulation of significant genes, compared to CON treatment. In LT, vitamin E supplementation caused down-regulation of genes related to intracellular signaling cascade. Functional analysis of SF showed that vitamin E supplementation caused up-regulation of the lipid biosynthesis process, cholesterol, and sterol and steroid biosynthesis, and it down-regulated genes related to the stress response. CONCLUSIONS: Different gene expression patterns were found between the SF and LT, suggesting tissue specific responses to vitamin E supplementation. Our study enabled us to identify novel genes and metabolic pathways related to vitamin E metabolism that might be implicated in meat quality. Further exploration of these genes and vitamin E could lead to a better understanding of how vitamin E affects the oxidative process that occurs in manufactured meat products.

Interpreting mismatches between linguistic and genetic patterns among speakers of Tanimuka (Eastern Tukanoan) and Yukuna (Arawakan).

Northwestern Amazonia is home to a great degree of linguistic diversity, and the human societies in that region are part of complex networks of interaction that predate the arrival of Europeans. This study investigates the population and language contact dynamics between two languages found within this region, Yukuna and Tanimuka, which belong to the Arawakan and Tukanoan language families, respectively. We use evidence from linguistics, ethnohistory, ethnography and population genetics to provide new insights into the contact dynamics between these and other human groups in NWA. Our results show that the interaction between these groups intensified in the last 500 years, to the point that it is difficult to differentiate between them genetically. However, this close interaction has led to more substantial contact-induced language changes in Tanimuka than in Yukuna, consistent with a scenario of language shift and asymmetrical power relations.

Exploring the ovine sperm transcriptome by RNAseq techniques. I Effect of seasonal conditions on transcripts abundance.

Understanding the cell molecular changes occurring as a results of climatic circumstances is crucial in the current days in which climate change and global warming are one of the most serious challenges that living organisms have to face. Sperm are one of the mammals' cells most sensitive to heat, therefore evaluating the impact of seasonal changes in terms of its transcriptional activity can contribute to elucidate how these cells cope with heat stress events. We sequenced the total sperm RNA from 64 ejaculates, 28 collected in summer and 36 collected in autumn, from 40 Manchega rams. A highly rich transcriptome (11,896 different transcripts) with 90 protein coding genes that exceed an average number of 5000 counts were found. Comparing transcriptome in the summer and autumn ejaculates, 236 significant differential abundance genes were assessed, most of them (228) downregulated. The main functions that these genes are related to sexual reproduction and negative regulation of protein metabolic processes and kinase activity. Sperm response to heat stress supposes a drastic decrease of the transcriptional activity, and the upregulation of only a few genes related with the basic functions to maintain the organisms' homeostasis and surviving. Rams' spermatozoids carry remnant mRNAs which are retrospectively indicators of events occurring along the spermatogenesis process, including abiotic factors such as environmental temperature.

Fungal and ciliate protozoa are the main rumen microbes associated with methane emissions in dairy cattle.

BACKGROUND: Mitigating the effects of global warming has become the main challenge for humanity in recent decades. Livestock farming contributes to greenhouse gas emissions, with an important output of methane from enteric fermentation processes, mostly in ruminants. Because ruminal microbiota is directly involved in digestive fermentation processes and methane biosynthesis, understanding the ecological relationships between rumen microorganisms and their active metabolic pathways is essential for reducing emissions. This study analysed whole rumen metagenome using long reads and considering its compositional nature in order to disentangle the role of rumen microbes in methane emissions. RESULTS: The ß-diversity analyses suggested a subtle association between methane production and overall microbiota composition (0.01 < R2 < 0.02). Differential abundance analysis identified 36 genera and 279 KEGGs as significantly associated with methane production (Padj < 0.05). Those genera associated with high methane production were Eukaryota from Alveolata and Fungi clades, while Bacteria were associated with low methane emissions. The genus-level association network showed 2 clusters grouping Eukaryota and Bacteria, respectively. Regarding microbial gene functions, 41 KEGGs were found to be differentially abundant between low- and high-emission animals and were mainly involved in metabolic pathways. No KEGGs included in the methane metabolism pathway (ko00680) were detected as associated with high methane emissions. The KEGG network showed 3 clusters grouping KEGGs associated with high emissions, low emissions, and not differentially abundant in either. A deeper analysis of the differentially abundant KEGGs revealed that genes related with anaerobic respiration through nitrate degradation were more abundant in low-emission animals. CONCLUSIONS: Methane emissions are largely associated with the relative abundance of ciliates and fungi. The role of nitrate electron acceptors can be particularly important because this respiration mechanism directly competes with methanogenesis. Whole metagenome sequencing is necessary to jointly consider the relative abundance of Bacteria, Archaea, and Eukaryota in the statistical analyses. Nutritional and genetic strategies to reduce CH4 emissions should focus on reducing the relative abundance of Alveolata and Fungi in the rumen. This experiment has generated the largest ONT ruminal metagenomic dataset currently available.

Use of Maximum Likelihood-Mixed Models to select stable reference genes: a case of heat stress response in sheep.

BACKGROUND: Reference genes with stable expression are required to normalize expression differences of target genes in qPCR experiments. Several procedures and companion software have been proposed to find the most stable genes. Model based procedures are attractive because they provide a solid statistical framework. NormFinder, a widely used software, uses a model based method. The pairwise comparison procedure implemented in GeNorm is a simpler procedure but one of the most extensively used. In the present work a statistical approach based in Maximum Likelihood estimation under mixed models was tested and compared with NormFinder and geNorm softwares. Sixteen candidate genes were tested in whole blood samples from control and heat stressed sheep. RESULTS: A model including gene and treatment as fixed effects, sample (animal), gene by treatment, gene by sample and treatment by sample interactions as random effects with heteroskedastic residual variance in gene by treatment levels was selected using goodness of fit and predictive ability criteria among a variety of models. Mean Square Error obtained under the selected model was used as indicator of gene expression stability. Genes top and bottom ranked by the three approaches were similar; however, notable differences for the best pair of genes selected for each method and the remaining genes of the rankings were shown. Differences among the expression values of normalized targets for each statistical approach were also found. CONCLUSIONS: Optimal statistical properties of Maximum Likelihood estimation joined to mixed model flexibility allow for more accurate estimation of expression stability of genes under many different situations. Accurate selection of reference genes has a direct impact over the normalized expression values of a given target gene. This may be critical when the aim of the study is to compare expression rate differences among samples under different environmental conditions, tissues, cell types or genotypes. To select reference genes not only statistical but also functional and biological criteria should be considered. Under the method here proposed SDHA/MDH1 have arisen as the best set of reference genes to be used in qPCR assays to study heat shock in ovine blood samples.

Changes in HSP gene and protein expression in natural scrapie with brain damage.

Heat shock proteins (Hsp) perform cytoprotective functions such as apoptosis regulation and inflammatory response control. These proteins can also be secreted to the extracellular medium, acting as inflammatory mediators, and their chaperone activity permits correct folding of proteins and avoids the aggregation of anomalous isoforms. Several studies have proposed the implication of Hsp in prion diseases. We analysed the gene expression and protein distribution of different members of the Hsp27, Hsp70, and Hsp90 families in the central nervous system of sheep naturally infected with scrapie. Different expression profiles were observed in the areas analysed. Whereas changes in transcript levels were not observed in the cerebellum or medulla oblongata, a significant decrease in HSP27 and HSP90 was detected in the prefrontal cortex. In contrast, HSP73 was over-expressed in diencephalons of scrapie animals. Western blotting did not reveal significant differences in Hsp90 and Hsp70 protein expression between scrapie and control animals. Expression rates identified by real-time RT-PCR and western blotting were compared with the extent of classical scrapie lesions using stepwise regression. Changes in Hsp gene and protein expression were associated with prion protein deposition, gliosis and spongiosis rather than with apoptosis. Finally, immunohistochemistry revealed intense Hsp70 and Hsp90 immunolabelling in Purkinje cells of scrapie sheep. In contrast, controls displayed little or no staining in these cells. The observed differences in gene expression and protein distribution suggest that the heat shock proteins analysed play a role in the natural form of the disease.

Microsatellite based genetic diversity and population structure of the endangered Spanish Guadarrama goat breed.

BACKGROUND: Assessing genetic biodiversity and population structure of minor breeds through the information provided by neutral molecular markers, allows determination of their extinction risk and to design strategies for their management and conservation. Analysis of microsatellite loci is known to be highly informative in the reconstruction of the historical processes underlying the evolution and differentiation of animal populations. Guadarrama goat is a threatened Spanish breed which actual census (2008) consists of 3057 females and 203 males distributed in 22 populations more or less isolated. The aim of this work is to study the genetic status of this breed through the analysis of molecular data from 10 microsatellites typed in historic and actual live animals. RESULTS: The mean expected heterozygosity across loci within populations ranged from 0.62 to 0.77. Genetic differentiation measures were moderate, with a mean FST of 0.074, GST of 0.081 and RST of 0.085. Percentages of variation among and within populations were 7.5 and 92.5, respectively. Bayesian clustering analyses pointed out a population subdivision in 16 clusters, however, no correlation between geographical distances and genetic differences was found. Management factors such as the limited exchange of animals between farmers (estimated gene flow Nm = 3.08) mostly due to sanitary and social constraints could be the major causes affecting Guadarrama goat population subdivision. CONCLUSION: Genetic diversity measures revealed a good status of biodiversity in the Guadarrama goat breed. Since diseases are the first cause affecting the census in this breed, population subdivision would be an advantage for its conservation. However, to maintain private alleles present at low frequencies in such small populations minimizing the inbreeding rate, it would necessitate some mating designs of animals carrying such alleles among populations. The systematic use of molecular markers will facilitate the comprehensive management of these populations, which in combination with the actual breeding program to increase milk yield, will constitute a good strategy to preserve the breed.

Gene Expression and Fatty Acid Profiling in Longissimus thoracis Muscle, Subcutaneous Fat, and Liver of Light Lambs in Response to Concentrate or Alfalfa Grazing.

A better understanding of gene expression and metabolic pathways in response to a feeding system is critical for identifying key physiological processes and genes associated with polyunsaturated fatty acid (PUFA) content in lamb meat. The main objective of this study was to investigate transcriptional changes in L. thoracis (LT) muscle, liver, and subcutaneous fat (SF) of lambs that grazed alfalfa (ALF) and concentrate-fed (CON) slaughtered at 23 kg and using the Affymetrix Ovine Gene 1.1 ST whole-genome array. The study also evaluated the relationship between meat traits in LT muscle, including color, pigments and lipid oxidation during 7 days of display, α-tocopherol content, intramuscular fat (IMF) content and the fatty acid (FA) profile. Lambs that grazed on alfalfa had a greater α-tocopherol concentration in plasma than CON lambs (P < 0.05). The treatment did not affect the IMF content, meat color or pigments (P > 0.05). Grazing increased the α-tocopherol content (P < 0.001) and decreased lipid oxidation on day 7 of display (P < 0.05) in LT muscle. The ALF group contained a greater amount of conjugated linoleic acid (CLA), C18:3 n-3, C20:5 n-3, C22:5 n-3, and C22:6 n-3 than did the CON group (P < 0.05). We identified 41, 96 and four genes differentially expressed in LT muscle, liver, and subcutaneous fat, respectively. The most enriched biological processes in LT muscle were skeletal muscle tissue development, being the genes related to catabolic and lipid processes downregulated, except for CPT1B, which was upregulated in the ALF lambs. Animals grazing alfalfa had lower expression of desaturase enzymes in the liver (FADS1 and FADS2), which regulate unsaturation of fatty acids and are directly involved in the metabolism of n-3 PUFA series. The results found in the current study showed that ingesting diets richer in n-3 PUFA might have negative effects on the de novo synthesis of n-3 PUFA by downregulating the FADS1 and FADS2 expression. However, feeding diets poorer in n-3 PUFA can promote fatty acid desaturation, which makes these two genes attractive candidates for altering the content of PUFAs in meat.

Genome-wide association studies for reproductive seasonality traits in Rasa Aragonesa sheep breed.

Sheep breeds from Mediterranean area show reproductive seasonal patterns of oestrous behaviour and ovulatory activity, mainly regulated by variation in the photoperiod. Maximal reproductive activity is associated with short days from August to March. The aim of this study therefore was, to identify new SNPs and genes associated to reproductive seasonality in sheep by using the Illumina OvineSNP50 Beadchip. A total of 239 adult Rasa Aragonesa breed ewes from one flock were controlled from January to August. Three reproductive seasonality traits were considered: the total days of anoestrus (TDA), based on weekly individual plasma progesterone levels and defined as the sum of days in anoestrus, considering anoestrus those periods with three or more consecutive P4 concentrations lower than 0.5 ng/ml; the progesterone cycling months (P4CM), defined for each ewe as the rate of cycling months between January and August based on progesterone determinations and the oestrus cycling months (OCM), defined for each ewe as the rate of months cycling between January and August based on oestrus records. Genotyping of 123 ewes was performed with the OvineSNP50 Infinium Beadchip. After the quality control (QC) performed on the raw genotypes, a total of 47,206 SNPs distributed over the 27 ovine chromosomes and 110 ewes were included in subsequent analyses. Principal component analysis revealed a substructure within the total dataset and identified 4 principal clusters in the experimental flock. None of the SNPs overcame the genome-wide significance level (P = 1.06 × 10-6). However, the SNPs OAR4_66002395 (9.41E-6), and OAR8_25877010 (1.86E-5) reached the genome-wide suggestive significance level (set to 2.32 × 10-5) for TDA and P4CM traits, respectively, while OAR23_14608581 was significant for both TDA (2.02E-5) and P4CM (1.05E-5) traits. Five SNPs evidenced association at chromosome-wise level: SNPs OAR4_66002395, OAR23_14608581 and s20800 (DTA), and OAR8_25877010, OAR23_14608581 and s48474 (P4CM). Several genes related to circadian and circannual rhythms were found close to these SNPs: NPSR1 (SNP OAR4_66002395), HS3ST5 (SNP OAR8_25877010), RPTOR (SNP s48474), and NPTX1 (SNP s48474) and could be considered as candidate gene related to TDA and P4CM traits.

Ovine HSP90AA1 expression rate is affected by several SNPs at the promoter under both basal and heat stress conditions.

The aim of this work was to investigate the association between polymorphisms located at the HSP90AA1 ovine gene promoter and gene expression rate under different environmental conditions, using a mixed model approach. Blood samples from 120 unrelated rams of the Manchega sheep breed were collected at three time points differing in environmental conditions. Rams were selected on the basis of their genotype for the transversion G/C located 660 base pairs upstream the gene transcription initiation site. Animals were also genotyped for another set of 6 SNPs located at the gene promoter. Two SNPs, G/C-660 and A/G-444, were associated with gene overexpression resulting from heat stress. The composed genotype CC-660-AG-444 was the genotype having the highest expression rates with fold changes ranging from 2.2 to 3.0. The genotype AG-522 showed the highest expression levels under control conditions with a fold change of 1.4. Under these conditions, the composed genotype CC-601-TT-524-AG-522-TT-468 is expected to be correlated with higher basal expression of the gene according to genotype frequencies and linkage disequilibrium values. Some putative transcription factors were predicted for binding sites where the SNPs considered are located. Since the expression rate of the gene under alternative environmental conditions seems to depend on the composed genotype of several SNPs located at its promoter, a cooperative regulation of the transcription of the HSP90AA1 gene could be hypothesized. Nevertheless epigenetic regulation mechanisms cannot be discarded.

Time trends, environmental factors and genetic basis of semen traits collected in Holstein bulls under commercial conditions.

The fact that results of artificial insemination (AI) are declining in highly selected dairy cattle populations has added a renewed interest to the evaluation of male fertility. Data from 42,348 ejaculates collected from 1990 to 2007 on 502 Holstein bulls were analysed in a Bayesian framework to provide estimates of the evolution of semen traits routinely collected in AI centres throughout the last decades of intense selection for production traits and estimate genetic parameters. The traits under consideration were volume (VOL), concentration (CONC), number of spermatozoa per ejaculate (NESPZ), mass motility score (MM), individual motility (IM), and post-thawing motility (PTM). The environmental factors studied were year-season and week of collection, which account for changes in environmental and technical conditions along time, age at collection, ejaculate order, time from previous collection (TPC) and time between collection and freezing (TCF) (only for PTM). Bull's inbreeding coefficient (Fi), bull's permanent environmental and additive genetic effects were also considered. The use of reduced models was evaluated using the Bayes factor. For all the systematic effects tested, strong or very strong evidence in favour of including the effect in the model was obtained, except for Fi for motility traits and TCF for PTM. No systematic time trends for environment or bull effects were observed, except for PTM, which showed an increasing environmental trend, associated with improvements in freezing-thawing protocols. Heritability estimates were moderate (0.16-0.22), except for IM, which presented a low value (0.07). Genetic correlations among motilities and between motilities and CONC were large and positive [0.38-0.87], VOL showed a negative correlation with CONC (-0.13) but with ample HPD 95%. The magnitude of heritabilities would allow an efficient selection if required and grants the use of these traits as indicators of the sperm viability component of bulls breeding soundness.

Polymorphisms in the HSP90AA1 5' flanking region are associated with scrapie incubation period in sheep.

Susceptibility to scrapie is mainly controlled by point mutations at the PRNP locus. However, additional quantitative trait loci (QTL) have been identified across the genome including a region in OAR18. The gene which encodes the inducible form of the cytoplasmic Hsp90 chaperone (HSP90AA1) maps within this region and seems to be associated with the resistance/susceptibility to scrapie in sheep. Here, we have analyzed several polymorphisms which were previously described in the ovine HSP90AA1 5' flanking region and in intron 10 in two naturally scrapie infected Romanov sheep populations. First, we have studied 58 ARQ/VRQ animals pertaining to the sire family where the QTL influencing scrapie incubation period in OAR18 was detected. We have found a significant association between polymorphisms localized at -660 and -528 in the HSP90AA1 5' flanking region and the scrapie incubation period. These two polymorphisms have also been studied in a second sample constituted by 62 VRQ/VRQ sheep showing an extreme incubation period. Results are concordant with the first dataset. Finally, we have studied the HSP90AA1 expression in scrapie and control animals (N = 41) with different HSP90AA1 genotypes by real time PCR on blood samples. The HSP90AA1 expression rate was equivalent in CC(-600)AA(-528) and CG(-600)AG(-528) scrapie resistant animals (ARR/ARR) and was higher in their CC(-600)AA(-528) than in their CG(-600)AG(-528) scrapie susceptible counterparts (VRQ/VRQ). Our results support the hypothesis that the ovine HSP90AA1 gene acts as a modulator of scrapie susceptibility, contributing to the observed differences in the incubation period of scrapie infected animals with the same PRNP genotype.

A SNP in the HSP90AA1 gene 5' flanking region is associated with the adaptation to differential thermal conditions in the ovine species.

Molecular chaperones have long been understood to be preferentially transcribed in response to multiple perturbations of the cellular homeostasis. In this study, several polymorphisms in the gene encoding the inducible form of the cytoplasmic Hsp90 (HSP90AA1) were addressed in 24 sheep breeds reared in different climatic regions of Europe, Africa, and Asia. Significant differences in the genotype frequencies for a C/G single nucleotide polymorphism (SNP) located at position -660 in the HSP90AA1 5'flanking region were found between the different breeds. Regression analyses reflected significant correlations (from 0.41 to 0.62) between the alternative genotypes of this polymorphism and several climatic and geographic variables characteristic of the regions where these breeds are reared. Real-time analysis revealed that animals bearing the CC(-660) genotype presented higher expression levels than those presenting the CG(-660) or GG(-660) in summer, but not in spring. Mutation at -660 site seems to affect HSP90AA1 transcription rates which could have important effects on the adaptation to different environmental conditions in sheep. Thus, the variability found in the genotype frequencies for the SNP at -660 in the ovine HSP90AA1 locus could be the result of the different environmental pressures occurring in the regions where these breed are maintained.

Structural and functional analysis of the HSP90AA1 gene: distribution of polymorphisms among sheep with different responses to scrapie.

Scrapie is a transmissible spongiform encephalopathy in sheep and goats. Susceptibility to this neurodegenerative disease is mainly controlled by point mutations at the PRNP locus. Other genes, apart from PRNP, have been reported to modulate resistance/susceptibility to scrapie. On the basis of several studies in Alzheimer and different transmissible spongiform encephalopathy models, HSP90AA1 was chosen as a putative positional and functional candidate gene that might be involved in the polygenic variance mentioned above. In the present work, the ovine HSP90AA1 gene including the promoter and other regulatory regions has been isolated and characterized. Several sequence polymorphisms have also been identified. FISH-mapping localized the HSP90AA1 gene on ovine chromosome OAR19q24dist, which was confirmed by linkage analysis. This chromosome region has been shown to include a quantitative trait loci (QTL) for scrapie incubation period in sheep. Expression analyses were carried out in spleen and cerebellum samples. No differences in the expression of the HSP90AA1 gene were found in any of these tissues (p > 0.05) between control and infected animal samples. Nevertheless, association analyses revealed that several polymorphisms in the 5' and 3' regions of the HSP90AA1 gene were differentially distributed among animals with different responses to scrapie infection. Thus, results presented here support the hypothesis that HSP90AA1 could be a positional and functional candidate gene modulating the response to scrapie in sheep.

Structural and functional analysis of the ovine laminin receptor gene (RPSA): Possible involvement of the LRP/LR protein in scrapie response.

Scrapie is a prion disease affecting sheep and goats. Susceptibility to this neurodegenerative disease shows polygenic variance. The involvement of the laminin receptor (LRP/LR) in the metabolism and propagation of prions has previously been demonstrated. In the present work, the ovine laminin receptor gene (RPSA) was isolated, characterized, and mapped to ovine chromosome OAR19q13. Real-time RT-PCR revealed a significant decrease in RPSA mRNA in cerebellum after scrapie infection. Conversely, no differences were detected in other brain regions such as diencephalon and medulla oblongata. Association analysis showed that a polymorphism reflecting the presence of a RPSA pseudogene was overrepresented in a group of sheep resistant to scrapie infection. No amino acid change in the LRP/LR protein was found in the 126 sheep analyzed. However, interesting amino acid positions (241, 272, and 290), which could participate in the species barrier to scrapie and maybe to other transmissible spongiform encephalopathies, were identified by comparing LRP/LR sequences from various mammals with variable levels of resistance to scrapie.

IL-1 family members as candidate genes modulating scrapie susceptibility in sheep: localization, partial characterization, and expression.

Scrapie (SC) is a transmissible spongiform encephalopathy (TSE) in sheep and goats. Susceptibility to this neurodegenerative disease is controlled mainly by point mutations at the PRNP locus. Other genes, apart from PRNP, have been reported to modulate resistance/susceptibility to SC. On the basis of several studies on Alzheimer's disease and different TSE models, and of requirement for correct homeostasis of cytokines in brain, IL1B and IL1RN were chosen as putative positional and functional candidate genes that might be involved in the polygenic variance mentioned above. In the present work, ovine IL1B and IL1RN genes were partially isolated and characterized, including promoter and other regulatory regions. In addition, several sequence polymorphisms were identified. Furthermore, their cytogenetic positions on sheep chromosomes were determined by FISH and confirmed by linkage analysis, localizing both genes in OAR3p22, a region previously described as carrying a QTL for SC incubation period in sheep. Finally, expression analyses were carried out in eight naturally SC-infected and five uninfected sheep with the same genotype for PRNP (ARQ/ARQ). This comparison was performed using real-time RT-PCR in samples of spleen and cerebellum. Results showed differences in the expression of both cytokines in cerebellum (p < 0.05) but not in spleen (p > 0.05).