In the present study, genetic analyses of diversity and differentiation were performed on four Basque-Navarrese semiferal native horse breeds. In total, 417 animals were genotyped for 12 microsatellite markers. Mean heterozygosity was higher than in other horse breeds, surely as a consequence of management. Although the population size of some of these breeds has declined appreciably in the past century, no genetic bottleneck was detected in any of the breeds, possibly because it was not narrow enough to be detectable. In the phylogenetic tree, the Jaca Navarra breed was very similar to the Pottoka, but appeared to stand in an intermediate position between this and the meat breeds. Assuming that Pottoka is the breed less affected by admixture, the others gradually distanced themselves from it through varying influences from outside breeds, among other factors. In a comparative study with other breeds, the French breeds Ardanais, Comtois, and Breton were the closest to the four native breeds. Three different approaches for evaluating the distribution of genetic diversity were applied. The high intrabreed variability of Euskal Herriko Mendiko Zaldia (EHMZ) was pointed out in these analyses. In our opinion, cultural, economic, and scientific factors should also be considered in the management of these horse breeds.

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.

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.

Endogenous retroviruses (ERVs) are the proviral phase of exogenous retroviruses that become integrated into a host germ line. They can play an important role in the host genome. Bioinformatic tools have been used to detect ERVs in several vertebrates, primarily primates and rodents. Less information is available regarding ERVs in other mammalian groups, and the source of this information is basically experimental. We analyzed the genome of the cow (Bos taurus) using three different methods. A BLAST-based method detected 928 possible ERVs, LTR_STRUC detected 4,487 elements flanked by long terminal repeats (LTRs), and Retrotector detected 9,698 ERVs. The ERVs were not homogeneously distributed across chromosomes; the number of ERVs was positively correlated with chromosomal size and negatively correlated with chromosomal GC content. The bovine ERVs (BoERVs) were classified into 24 putative families, with 20 of them not previously described. One of these new families, BoERV1, was the most abundant family and appeared to be specific to ruminants. An analysis of representatives of ERV families from rodents, primates, and ruminants showed a phylogenetic relationship following their hosts' relationships. This study demonstrates the importance of using multiple methods when trying to identify new ERVs and shows that the number of bovine ERV families is not as limited as previously thought.

No abstract

BackgroundHuman 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.ObjectivesWe 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.MethodsFrom 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.ResultsGenetic 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.ConclusionsMultiple 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.

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.

Ovine pulmonary adenocarcinoma (OPA) and Maedi-Visna (Maedi) are two chronic respiratory diseases of retroviral origin which occur worldwide. It is known that different host genetic factors influence the outcome of viral infections. To determine if variation in the Mhc-DRB1 gene was associated with progression to these ovine diseases, sheep lungs with and without OPA and Maedi lesions were collected. A sequence-based method was applied and 40 different alleles were detected in the sample analysed. In the allele-by-allele association analysis, allele DRB1*0325 had a significant association with susceptibility to Maedi (P = 0.045). For OPA, DRB1*0143 and DRB1*0323 were significantly associated with susceptibility (P = 0.024 and P = 0.029), and allele DRB1*0702 was significantly associated with resistance (P = 0.012). Based on these results, the Mhc-DRB1 alleles were classified by effect in three categories-susceptible (S), resistant (R) and neutral (N)-and animals were reassigned the genotypes as S/S, S/R, S/N, R/R, R/N and N/N. In a second analysis, penalised logistic regression models including a flock effect were run. In Maedi, significant association was detected for the N/S heterozygote (P = 0.0007), but not for the S/S homozygote, probably as a result of the low number of S/S animals. In OPA, association was detected for both the S/S and R/R homozygotes (P = 0.005 and P = 0.047). This allele grouping method may be applied in association studies with highly variable genes. This is the first study demonstrating significant associations between sheep Mhc-DRB1 alleles and susceptibility to OPA and Maedi. Therefore, both diseases are suitable candidates for more comprehensive genetic studies.