Identification and characterization of miRNAs in spleens of sheep subjected to repetitive vaccination.

Accumulative evidence has shown that short non-coding RNAs such as miRNAs can regulate the innate and adaptive immune responses. Aluminium hydroxide is a commonly used adjuvant in human and veterinary vaccines. Despite its extended use, its mechanism of action is not fully understood and very few in vivo studies have been done to enhance understanding at the molecular level. In this work, we took advantage of a previous long-term experiment in which lambs were exposed to three different treatments by parallel subcutaneous inoculations with aluminium-containing commercial vaccines, an equivalent dose of aluminium or mock injections. Spleen samples were used for miRNA-seq. A total of 46 and 16 miRNAs were found differentially expressed when animals inoculated with commercial vaccines or the adjuvant alone were compared with control animals, respectively. Some miRNAs previously related to macrophage polarization were found dysregulated exclusively by the commercial vaccine treatment but not in the aluminium inoculated animals. The dysregulated miRNAs in vaccine group let-7b-5p, miR-29a-3p, miR-27a and miR-101-3p are candidates for further research, since they may play key roles in the immune response induced by aluminium adjuvants added to vaccines. Finally, protein-protein interaction network analysis points towards leucocyte transendothelial migration as a specific mechanism in animals receiving adjuvant only.

The sheep miRNAome: Characterization and distribution of miRNAs in 21 tissues.

MicroRNAs (miRNAs) are evolutionarily conserved small non-coding RNAs that regulate gene expression. In livestock species, miRNAs also show great potential as biomarkers for animal health and product quality. In sheep, few miRNAs have been described in comparison with other livestock species or model organisms. We uniformly analyzed 172 public ovine small RNA sequencing datasets from 21 different tissues in order to predict conserved and novel miRNA precursors and profile their expression patterns. In addition to the 106 annotated sheep miRNAs, 1047 previously unannotated miRNA precursor sequences were detected and 41 % of them were assigned an orthologue from other close species. Considering expression levels, a set of miRNAs with high sequence conservation were detected in all tissues, while 733 mature miRNAs were robustly expressed in at least one tissue. 270 miRNAs showed high tissue specificity index values. Brain, male reproductive tissues and PBMCs showed the most distinct expression patterns. Strikingly, over one hundred precursors from the ruminant specific family of mir-2284/mir-2285 miRNAs were found, which were enriched in immune related tissues. This work supports the known high conservation of many miRNAs, but also highlights the potential of clade-specific innovations in ruminant evolution.

Transcriptomic meta-analysis reveals unannotated long non-coding RNAs related to the immune response in sheep.

Long non-coding RNAs (lncRNAs) are involved in several biological processes, including the immune system response to pathogens and vaccines. The annotation and functional characterization of lncRNAs is more advanced in humans than in livestock species. Here, we take advantage of the increasing number of high-throughput functional experiments deposited in public databases in order to uniformly analyse, profile unannotated lncRNAs and integrate 422 ovine RNA-seq samples from the ovine immune system. We identified 12302 unannotated lncRNA genes with support from independent CAGE-seq and histone modification ChIP-seq assays. Unannotated lncRNAs showed low expression levels and sequence conservation across other mammal species. There were differences in expression levels depending on the genomic location-based lncRNA classification. Differential expression analyses between unstimulated and samples stimulated with pathogen infection or vaccination resulted in hundreds of lncRNAs with changed expression. Gene co-expression analyses revealed immune gene-enriched clusters associated with immune system activation and related to interferon signalling, antiviral response or endoplasmic reticulum stress. Besides, differential co-expression networks were constructed in order to find condition-specific relationships between coding genes and lncRNAs. Overall, using a diverse set of immune system samples and bioinformatic approaches we identify several ovine lncRNAs associated with the response to an external stimulus. These findings help in the improvement of the ovine lncRNA catalogue and provide sheep-specific evidence for the implication in the general immune response for several lncRNAs.

Identification of sheep lncRNAs related to the immune response to vaccines and aluminium adjuvants.

BACKGROUND: Long non-coding RNAs (lncRNAs) are involved in several immune processes, including the immune response to vaccination, but most of them remain uncharacterised in livestock species. The mechanism of action of aluminium adjuvants as vaccine components is neither not fully understood. RESULTS: We built a transcriptome from sheep PBMCs RNA-seq data in order to identify unannotated lncRNAs and analysed their expression patterns along protein coding genes. We found 2284 novel lncRNAs and assessed their conservation in terms of sequence and synteny. Differential expression analysis performed between animals inoculated with commercial vaccines or aluminium adjuvant alone and the co-expression analysis revealed lncRNAs related to the immune response to vaccines and adjuvants. A group of co-expressed genes enriched in cytokine signalling and production highlighted the differences between different treatments. A number of differentially expressed lncRNAs were correlated with a divergently located protein-coding gene, such as the OSM cytokine. Other lncRNAs were predicted to act as sponges of miRNAs involved in immune response regulation. CONCLUSIONS: This work enlarges the lncRNA catalogue in sheep and puts an accent on their involvement in the immune response to repetitive vaccination, providing a basis for further characterisation of the non-coding sheep transcriptome within different immune cells.

Gene Expression Profiling Reveals New Pathways and Genes Associated with Visna/Maedi Viral Disease.

Visna/Maedi virus (VMV) is a lentivirus that infects the cells of the monocyte/macrophage lineage in sheep, goats and wild ruminants. Infection with VMV causes a multisystemic inflammatory disorder, which includes pneumonia, encephalitis, mastitis or arthritis. The immune response to VMV infection is complex, and the infection and pathogenesis of this virus are not totally characterized yet. In this work, a gene expression microarray was used to identify the differentially expressed genes in VMV infection and disease development by comparing sheep with different serologic status and with presence of VM-characteristic clinical lesions. The expression profile analysis has revealed many interesting genes that may be associated with the viral infection process. Among them, the OXT gene appeared significantly up-regulated, so the oxytocin-secreting system could play an essential role in VM disease. Moreover, some of the most significantly enriched functions in up-regulated genes appeared the complement pathway, which (in combination with the Toll-like receptor signaling network) could compose a mechanism in the VMV pathogenesis. Identifying the host genetic factors associated with VMV infection can be applied to develop strategies for preventing infection and develop effective vaccines that lead to therapeutic treatments.

The dystrophia myotonica WD repeat-containing protein DMWD and WDR20 differentially regulate USP12 deubiquitinase.

Despite its potential clinical relevance, the product of the DMWD (dystrophia myotonica, WD repeat containing) gene is a largely uncharacterized protein. The DMWD amino acid sequence is similar to that of WDR20, a known regulator of the USP12 and USP46 deubiquitinases (DUBs). Here, we apply a combination of in silico and experimental methods to investigate several aspects of DMWD biology. Molecular evolution and phylogenetic analyses reveal that WDR20 and DMWD, similar to USP12 and USP46, arose by duplication of a common ancestor during the whole genome duplication event in the vertebrate ancestor lineage. The analysis of public human gene expression datasets suggests that DMWD expression is positively correlated with USP12 expression in normal tissues and negatively correlated with WDR20 expression in tumors. Strikingly, a survey of the annotated interactome for DMWD and WDR20 reveals a largely nonoverlapping set of interactors for these proteins. Experimentally, we first confirmed that DMWD binds both USP12 and USP46 through direct coimmunoprecipitation of epitope-tagged proteins. We found that DMWD and WDR20 share the same binding interface in USP12, suggesting that their interaction with the DUB may be mutually exclusive. Finally, we show that both DMWD and WDR20 promote USP12 enzymatic activity, but they differentially modulate the subcellular localization of the DUB. Altogether, our findings suggest a model whereby mutually exclusive binding of DMWD and WDR20 to USP12 may lead to formation of deubiquitinase complexes with distinct subcellular localization, potentially targeting different substrate repertoires.

Novel circRNA discovery in sheep shows evidence of high backsplice junction conservation.

Circular RNAs (circRNAs) are covalently closed circular non-coding RNAs. Due to their structure, circRNAs are more stable and have longer half-lives than linear RNAs making them good candidates for disease biomarkers. Despite the scientific relevance of these molecules, the study of circRNAs in non-model organisms is still in its infancy. Here, we analyse total RNA-seq data to identify circRNAs in sheep from peripheral blood mononuclear cells (PBMCs) and parietal lobe cortex. Out of 2510 and 3403 circRNAs detected in parietal lobe cortex and in PBMCs, a total of 1379 novel circRNAs were discovered. Remarkably, around 63% of all detected circRNAs were found to be completely homologous to a circRNA annotated in human. Functional enrichment analysis was conducted for both tissues based on GO terms and KEGG pathways. The enriched terms suggest an important role of circRNAs from encephalon in synaptic functions and the involvement of circRNAs from PBMCs in basic immune system functions. In addition to this, we investigated the role of circRNAs in repetitive vaccination experiments via differential expression analysis and did not detect any significant relationship. At last, our results support both the miRNA sponge and the miRNA shuttle functions of CDR1-AS in sheep brain. To our knowledge, this is the first study on circRNA annotation in sheep PBMCs or parietal lobe cortex samples.

Whole transcriptome approach to evaluate the effect of aluminium hydroxide in ovine encephalon.

Aluminium hydroxide adjuvants are crucial for livestock and human vaccines. Few studies have analysed their effect on the central nervous system in vivo. In this work, lambs received three different treatments of parallel subcutaneous inoculations during 16 months with aluminium-containing commercial vaccines, an equivalent dose of aluminium hydroxide or mock injections. Brain samples were sequenced by RNA-seq and miRNA-seq for the expression analysis of mRNAs, long non-coding RNAs and microRNAs and three expression comparisons were made. Although few differentially expressed genes were identified, some dysregulated genes by aluminium hydroxide alone were linked to neurological functions, the lncRNA TUNA among them, or were enriched in mitochondrial energy metabolism related functions. In the same way, the miRNA expression was mainly disrupted by the adjuvant alone treatment. Some differentially expressed miRNAs had been previously linked to neurological diseases, oxidative stress and apoptosis. In brief, in this study aluminium hydroxide alone altered the transcriptome of the encephalon to a higher degree than commercial vaccines that present a milder effect. The expression changes in the animals inoculated with aluminium hydroxide suggest mitochondrial disfunction. Further research is needed to elucidate to which extent these changes could have pathological consequences.

RNA-Seq analysis of ileocecal valve and peripheral blood from Holstein cattle infected with Mycobacterium avium subsp. paratuberculosis revealed dysregulation of the CXCL8/IL8 signaling pathway.

Paratuberculosis is chronic granulomatous enteritis of ruminants caused by Mycobacterium avium subsp. paratuberculosis (MAP). Whole RNA-sequencing (RNA-Seq) is a promising source of novel biomarkers for early MAP infection and disease progression in cattle. Since the blood transcriptome is widely used as a source of biomarkers, we analyzed whether it recapitulates, at least in part, the transcriptome of the ileocecal valve (ICV), the primary site of MAP colonization. Total RNA was prepared from peripheral blood (PB) and ICV samples, and RNA-Seq was used to compare gene expression between animals with focal or diffuse histopathological lesions in gut tissues versus control animals with no detectable signs of infection. Our results demonstrated both shared, and PB and ICV-specific gene expression in response to a natural MAP infection. As expected, the number of differentially expressed (DE) genes was larger in the ICV than in the PB samples. Among the DE genes in the PB and ICV samples, there were some common genes irrespective of the type of lesion including the C-X-C motif chemokine ligand 8 (CXCL8/IL8), apolipoprotein L (APOLD1), and the interferon inducible protein 27 (IFI27). The biological processes (BP) enriched in the PB gene expression profiles from the cows with diffuse lesions included the killing of cells of other organism, defense response, immune response and the regulation of neutrophil chemotaxis. Two of these BP, the defense and immune response, were also enriched in the ICV from the cows with diffuse lesions. Metabolic analysis of the DE genes revealed that the N-glycan biosynthesis, bile secretion, one-carbon pool by folate and purine metabolism were significantly enriched in the ICV from the cows with focal lesions. In the ICV from cows with diffuse lesions; the valine, leucine and isoleucine degradation route, purine metabolism, vitamin digestion and absorption and the cholesterol routes were enriched. Some of the identified DE genes, BP and metabolic pathways will be studied further to develop novel diagnostic tools, vaccines and immunotherapeutics.

Expression analysis of lung miRNAs responding to ovine VM virus infection by RNA-seq.

BACKGROUND: MicroRNAs (miRNAs) are short endogenous, single-stranded, noncoding small RNA molecules of approximately 22 nucleotides in length. They regulate gene expression posttranscriptionally by silencing mRNA expression, thus orchestrating many physiological processes. The Small Ruminant Lentiviruses (SRLV) group includes the Visna Maedi Virus (VMV) and Caprine Arthritis Encephalitis (CAEV) viruses, which cause a disease in sheep and goats characterized by pneumonia, mastitis, arthritis and encephalitis. Their main target cells are from the monocyte/macrophage lineage. To date, there are no studies on the role of miRNAs in this viral disease. RESULTS: Using RNA-seq technology and bioinformatics analysis, the expression levels of miRNAs during different clinical stages of infection were studied. A total of 212 miRNAs were identified, of which 46 were conserved sequences in other species but found for the first time in sheep, and 12 were completely novel. Differential expression analysis comparing the uninfected and seropositive groups showed changes in several miRNAs; however, no significant differences were detected between seropositive asymptomatic and diseased sheep. The robust increase in the expression level of oar-miR-21 is consistent with its increased expression in other viral diseases. Furthermore, the target prediction of the dysregulated miRNAs revealed that they control genes involved in proliferation-related signalling pathways, such as the PI3K-Akt, AMPK and ErbB pathways. CONCLUSIONS: To the best of our knowledge, this is the first study reporting miRNA profiling in sheep in response to SRLV infection. The known functions of oar-miR-21 as a regulator of inflammation and proliferation appear to be a possible cause of the lesions caused in the sheep's lungs. This miRNA could be an indicator for the severity of the lung lesions, or a putative target for therapeutic intervention.

Molecular Signature of Aluminum Hydroxide Adjuvant in Ovine PBMCs by Integrated mRNA and microRNA Transcriptome Sequencing.

There have been few in vivo studies on the effect of aluminum hydroxide adjuvant and its influence on the immune response to vaccination. In this study, lambs received a parallel subcutaneous treatment with either commercial vaccines containing aluminum hydroxide or an equivalent dose of this compound only with the aim of identifying the activated molecular signature. Blood samples were taken from each animal at the beginning and at the end of the experiment and PBMCs isolated. Total RNA and miRNA libraries were prepared and sequenced. After alignment to the Oar3.1 reference genome and differential expression with 3 programs, gene enrichment modeling was performed. For miRNAs, miRBase and RNAcentral databases were used for detection and characterization. Three expression comparisons were made: vaccinated animals at the beginning and at the end of the treatment, adjuvanted animals at the same times, and animals of both treatments at the end of the experiment. After exposure to both treatments, a total of 2,473; 2,980 and 429 differentially expressed genes were identified in vaccinated animals, adjuvanted animals and animals at the end of both treatments, respectively. In both adjuvant and vaccine treated animals the NF-κB signaling pathway was enriched. On the other hand, it can be observed a downregulation of cytokines and cytokine receptors in the adjuvanted group compared to the vaccinated group at the final time, suggesting a milder induction of the immune response when the adjuvant is alone. As for the miRNA analysis, 95 miRNAs were detected: 64 previously annotated in Ovis aries, 11 annotated in Bos taurus and 20 newly described. Interestingly, 6 miRNAs were differentially expressed in adjuvant treated animals, and 3 and 1 in the other two comparisons. Lastly, an integrated miRNA-mRNA expression profile was developed, in which a miRNA-mediated regulation of genes related to DNA damage stimulus was observed. In brief, it seems that aluminum-containing adjuvants are not simple delivery vehicles for antigens, but also induce endogenous danger signals that can stimulate the immune system. Whether this contributes to long-lasting immune activation or to the overstimulation of the immune system remains to be elucidated.

Genome-wide reexamination of endogenous retroviruses in Rattus norvegicus.

Endogenous retroviruses (ERVs) are remnants of retroviral infections that are present in a large number of vertebrate genomes. Based on the proposal that the rat could act as a reservoir of retroviruses, rat ERVs were analysed in silico using a whole-genome approach. To enrich the detected ERV groups, we applied an upgraded approach based on the hidden Markov model. We found 2637 elements that were classified into the following groups: 9 groups of Class I; 15 of Class II, 7 of them previously described; 1 of Class III; and 3 groups whose classification was unclear but were distantly related to Class I. Sixteen ERV groups seemed to be specific to rat. The high number of rat-specific groups might be related to the contact of rats with retroviruses and their role as a reservoir. In addition, the env gene of the more extended groups seemed to be undetectable.

Lack of association between polymorphic copies of endogenous Jaagsiekte sheep retrovirus (enJSRVs) and Ovine Pulmonary Adenocarcinoma.

Ovine Pulmonary Adenocarcinoma (OPA) is a retrovirus-induced lung tumor of sheep, goat and mouflon, and its etiologic agent, Jaagsiekte sheep retrovirus (JSRV) is the only virus known to cause a naturally occurred lung adenocarcinoma. The oncogenic JSRV has several endogenous counterparts termed enJSRVs, some of which have been shown to interfere with JSRV replication at early and late stages of the retroviral cycle inhibiting JSRV exit from the cell, and thus, protecting sheep against the infection. In this work, Latxa sheep breed animals were classified depending on the presence/absence of OPA-characteristic clinical lesions in the lung. Using a PCR genotyping method and a logistic regression-based association study, five polymorphic enJSRV copies were analyzed in 49 OPA positive sheep and 124 control individuals. Our results showed that the frequency of the provirus enJSRV-16 is much higher in Latxa sheep breed than in other breeds, suggesting a recent proliferation of this provirus in the studied breed. However, no polymorphic enJSRV was found to be statistically associated with the susceptibility/resistance to OPA development.

Characterization and genotyping of the DRB1 gene of the major histocompatibility complex (MHC) in the Marmota monax, animal model of hepatitis B.

The major histocompatibility complex (MHC)-containing genes are among the most polymorphic in vertebrates. MHC genes code for proteins that are critical in the immune system response. In this study, the polymorphism of the second exon of the MHC class II DRB gene was characterized in the Eastern woodchuck (Marmota monax). Woodchucks chronically infected with the woodchuck hepatitis virus (WHV) represent the best available animal model for the study of chronic hepatitis B infection in humans. In the genotyped animals we found fifteen alleles, which were expressed in two independent loci and that were named DRB1A and DRB1B in this work. The 15 alleles investigated showed an elevated divergence. A significant excess of non-synonymous substitutions was detected, which could indicate that a historical positive selection is acting in the woodchuck DRB1 genes. This hypothesis was confirmed in our study by the high variability in or near the antigen binding sites (ABS) and by the results obtained in sequence variability analyses. This analysis identified the presence of a microsatellite sequence that is located at the start of the second intron, which could further allow the development of a fast and cheap semiautomatic sequencing method.

Genomic environment and digital expression of bovine endogenous retroviruses.

In this work, we attempt to identify the location and expression of bovine ERVs and the nature of nearby genes. As many as 1610 bovine genes contain a bovine ERV (BoERV) inserted into their introns. Most of these BoERVs present an antisense orientation, which could be a consequence of the detrimental effects of the sense insertion. Based on the overrepresentation of Gene Ontology terms, we determined that some genes located in the vicinity of BoERVs are related to viral response and chromatin assembly. In addition, we identified some genes that belonged to IFNs that are inserted in or between BoERVs, pointing out a possible role of BoERVs in some gene duplication. Based on EST database mining, significant expression was detected for BoERV genes in the thyroid and in embryos only. Transcription factor binding sites for Rxra, a steroid and thyroid hormone receptor, were detected in three of the most expressed BoERVs in the thyroid glands.

Evolutionary history of bovine endogenous retroviruses in the Bovidae family.

BACKGROUND: Endogenous retroviruses (ERVs) are genomic elements of retroviral origin that are present in the genomes of almost all vertebrates. In cattle, more than 13,000 elements related to ERVs have been detected, and based on the pol gene, 24 families or groups of bovine ERVs have been described. However, information about ERVs in other bovids and the presence of families of related bovine ERVs in different species of the Bovidae family is scarce. RESULTS: The 24 families of bovine ERVs previously detected in cattle (Bos taurus) were also detected in zebus (Bos indicus) and yaks (Bos grunniens). In addition, six new families, named BoERV25 to BoERV30, were detected in the three Bos species. Five more ruminant species were screened for related ERVs: 26 families were detected in these species, but four families (BoERV24, BoERV26, BoERV28 and BoERV29) were specific to cattle, zebus, yaks and buffalo. An analysis of the homology of the ERVs of cattle, zebus and yaks revealed that the level of LTR divergence was similar between ERVs from cattle and zebus but was less similar between with ERVs from cattle and yaks. In addition, purifying selection was detected in the genes and retroviral regions of clusters of ERVs of cattle, zebus and yaks. CONCLUSIONS: In this work, the 24 ERV families previously identified in cattle were also found in two other species in the Bos genus. In addition, six new bovine ERV families were detected. Based on LTR divergence, the most recently inserted families are from Class II. The divergence of the LTR, used as an indirect estimate of the ERV insertion time, seemed to be influenced by the differences in genome evolution since the divergence of the species. In addition, purifying selection could be acting on clusters of ERVs from different species.

Retroviral infections in sheep and goats: small ruminant lentiviruses and host interaction.

Small ruminant lentiviruses (SRLV) are members of the Retrovirus family comprising the closely related Visna/Maedi Virus (VMV) and the Caprine Arthritis-Encephalitis Virus (CAEV), which infect sheep and goats. Both infect cells of the monocyte/macrophage lineage and cause lifelong infections. Infection by VMV and CAEV can lead to Visna/Maedi (VM) and Caprine Arthritis-Encephalitis (CAE) respectively, slow progressive inflammatory diseases primarily affecting the lungs, nervous system, joints and mammary glands. VM and CAE are distributed worldwide and develop over a period of months or years, always leading to the death of the host, with the consequent economic and welfare implications. Currently, the control of VM and CAE relies on the control of transmission and culling of infected animals. However, there is evidence that host genetics play an important role in determining Susceptibility/Resistance to SRLV infection and disease progression, but little work has been performed in small ruminants. More research is necessary to understand the host-SRLV interaction.

Expression analysis of 13 ovine immune response candidate genes in Visna/Maedi disease progression.

Visna/Maedi virus (VMV) is a lentivirus that infects cells of the monocyte/macrophage lineage in sheep. Infection with VMV may lead to Visna/Maedi (VM) disease, which causes a multisystemic inflammatory disorder causing pneumonia, encephalitis, mastitis and arthritis. The role of ovine immune response genes in the development of VM disease is not fully understood. In this work, sheep of the Rasa Aragonesa breed were divided into two groups depending on the presence/absence of VM-characteristic clinical lesions in the aforementioned organs and the relative levels of candidate gene expression, including cytokines and innate immunity loci were measured by qPCR in the lung and udder. Sheep with lung lesions showed differential expression in five target genes: CCR5, TLR7, and TLR8 were up regulated and IL2 and TNFα down regulated. TNFα up regulation was detected in the udder.

Human papillomavirus (HPV) genotype 18 variants in patients with clinical manifestations of HPV related infections in Bilbao, Spain.

BACKGROUND: Human papillomavirus (HPV) variants differ in their biological and chemical properties, and therefore, may present differences in pathogenicity. Most authors classified variants based on the phylogenetic analysis of L1 region. Nevertheless, recombination in HPV samples is becoming a usual finding and thus, characterizing genetic variability in other regions should be essential. OBJECTIVES: We aimed to characterize the genetic variability of HPV 18 in 5 genomic regions: E6, E7, E4, L1 and the Upstream Regulatory Region (URR), working with both single infection and multiple HPV infection samples. Furthermore, we aimed to assess the prevalence of HPV 18 variants in our region and look for possible existence of recombination as well as analyze the relationship between these variants and the type of lesion. METHODS: From 2007 to 2010, Clinical Microbiology and Infection Control Department analyzed 44 samples which were positive for HPV 18. Genetic variability was determined in PCR products and variants were assigned to European, Asian-amerindian or African lineage. Recombination and association of variants with different types of lesion was studied. RESULTS: Genetic analysis of the regions revealed a total of 56 nucleotide variations. European, African and Asian-amerindian variants were found in 25/44 (56.8%), 10/44 (22.7%) and 5/44 (11.4%) samples, respectively. We detected the presence of recombinant variants in 2/44 (4.5%) cases. Samples taken from high-grade squamous intraepithelial lesions (H-SIL) only presented variants with specific-african substitutions. CONCLUSIONS: Multiple HPV infection, non-european HPV variants prevalence and existence of recombination are considered risk factors for HPV persistence and progression of intraepithelial abnormalities, and therefore, should be taken into consideration in order to help to design and optimize diagnostics protocols as well as improve epidemiologic studies.Our study is one of the few studies in Spain which analyses the genetic variability of HPV18 and we showed the importance of characterizing more than one genomic region in order to detect recombination and classify HPV variants properly.

Amino acid signatures in the Ovar-DRB1 peptide-binding pockets are associated with Ovine Pulmonary Adenocarcinoma susceptibility/resistance.

Potential relationships between amino acid motifs of various alleles of the ovine major histocompatibility complex DR (Ovar-DR) molecule and occurrence of clinical OPA caused by JSRV were investigated in a case-control study. Latxa sheep (n=132) screened for presence/absence of pulmonary OPA lesions were typed for their Ovar-DRB1 2nd exon alleles by PCR and sequence-based typing (PCR-SBT). The polymorphic amino acid residues derived from the obtained 34 DRB1 protein variants were subjected to a logistic regression-based association study. The amino acids at several positions showed significant associations with the presence/absence of pulmonary OPA lesions; some of the residues were located within the peptide binding cleft of the DRB molecule, including pockets P1, P4, P7 and P9.