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.

Characterization of the ovine ribosomal protein SA gene and its pseudogenes.

BACKGROUND: The ribosomal protein SA (RPSA), previously named 37-kDa laminin receptor precursor/67-kDa laminin receptor (LRP/LR) is a multifunctional protein that plays a role in a number of pathological processes, such as cancer and prion diseases. In all investigated species, RPSA is a member of a multicopy gene family consisting of one full length functional gene and several pseudogenes. Therefore, for studies on RPSA related pathways/pathologies, it is important to characterize the whole family and to address the possible function of the other RPSA family members. The present work aims at deciphering the RPSA family in sheep. RESULTS: In addition to the full length functional ovine RPSA gene, 11 other members of this multicopy gene family, all processed pseudogenes, were identified. Comparison between the RPSA transcript and these pseudogenes shows a large variety in sequence identities ranging from 99% to 74%. Only one of the 11 pseudogenes, i.e. RPSAP7, shares the same open reading frame (ORF) of 295 amino acids with the RPSA gene, differing in only one amino acid. All members of the RPSA family were annotated by comparative mapping and fluorescence in situ hybridization (FISH) localization. Transcription was investigated in the cerebrum, cerebellum, spleen, muscle, lymph node, duodenum and blood, and transcripts were detected for 6 of the 11 pseudogenes in some of these tissues. CONCLUSIONS: In the present work we have characterized the ovine RPSA family. Our results have revealed the existence of 11 ovine RPSA pseudogenes and provide new data on their structure and sequence. Such information will facilitate molecular studies of the functional RPSA gene taking into account the existence of these pseudogenes in the design of experiments. It remains to be investigated if the transcribed members are functional as regulatory non-coding RNA or as functional proteins.

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.

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).