Summary-data based Mendelian randomization identifies gene expression regulatory polymorphisms associated with bovine paratuberculosis by modulation of the nuclear factor Kappa ß (NF-κß)-mediated inflammatory response.

Genome-wide association studies (GWAS) have identified host genetic variants associated with paratuberculosis (PTB) susceptibility. Most of the GWAS-identified SNPs are in non-coding regions. Connecting these non-coding variants and downstream affected genes is a challenge and, up to date, only a few functional mutations or expression quantitative loci (cis-eQTLs) associated with PTB susceptibility have been identified. In the current study, the associations between imputed whole-genome sequence genotypes and whole RNA-Sequencing data from peripheral blood (PB) and ileocecal valve (ICV) samples of Spanish Holstein cows (N = 16) were analyzed with TensorQTL. This approach allowed the identification of 88 and 37 cis-eQTLs regulating the expression levels of 90 and 37 genes in PB and ICV samples, respectively (False discorey rate, FDR ≤ 0.05). Next, we applied summary-based data Mendelian randomization (SMR) to integrate the cis-eQTL dataset with GWAS data obtained from a cohort of 813 culled cattle that were classified according to the presence or absence of PTB-associated histopathological lesions in gut tissues. After multiple testing corrections (FDR ≤ 0.05), we identified two novel cis-eQTLs affecting the expression of the early growth response factor 4 (EGR4) and the bovine neuroblastoma breakpoint family member 6-like protein isoform 2 (MGC134040) that showed pleiotropic associations with the presence of multifocal and diffuse lesions in gut tissues; P = 0.002 and P = 0.017, respectively. While EGR4 acts as a brake on T-cell proliferation and cytokine production through interaction with the nuclear factor Kappa ß (NF-κß), MGC134040 is a target gene of NF-κß. Our findings provide a better understanding of the genetic factors influencing PTB outcomes, confirm that the multifocal lesions are localized/confined lesions that have different underlying host genetics than the diffuse lesions, and highlight regulatory SNPs and regulated-gene targets to design future functional studies.

Double immunofluorescence staining of whole-mount small intestinal mucosa samples as a tool for characterization of three-dimensional paratuberculosis granulomas.

Bovine paratuberculosis (PTB) is a chronic granulomatous enteritis, caused by Mycobacterium avium subsp. paratuberculosis (Map). The progression of PTB from subclinical to the clinical stage of the disease is determined locally at the level of the granuloma, a host defence hallmark against mycobacterial infection. Therefore, in-depth characterization of distinct cell populations controlling granuloma formation is critical to understanding PTB progression. Confocal laser scanning microscopy (CLSM) has been extensively used to visualize two or more proteins of interest concomitantly within a variety of cellular structures. As such, it is an invaluable tool for the correct identification and characterization of different cell populations. In this study, a novel approach, CLSM of whole-mount small intestinal mucosa samples, is used to characterize three-dimensional (3-D) paratuberculosis granulomas and epithelioid macrophages. Detailed optimized procedures to perform CLSM in whole mount small intestinal mucosa samples and also in formalin fixed paraffin embedded (FFPE) intestinal tissue sections of Holstein Friesian cows presenting different types of PTB-associated histological lesions are described.

Bovine Intelectin 2 Expression as a Biomarker of Paratuberculosis Disease Progression.

Paratuberculosis (PTB), a chronic granulomatous enteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP), is responsible for important economic losses in the dairy industry. Our previous RNA-sequencing (RNA-Seq) analysis showed that bovine intelectin 2 (ITLN2) precursor gene was overexpressed in ileocecal valve (ICV) samples of animals with focal (log2 fold-change = 10.6) and diffuse (log2 fold-change = 6.8) PTB-associated lesions compared to animals without lesions. This study analyzes the potential use of ITLN2, a protein that has been described as fundamental in the innate immune response to infections, as a biomarker of MAP infection. The presence of ITLN2 was investigated by quantitative immunohistochemical analysis of ICV samples of 20 Holstein Friesian cows showing focal (n = 5), multifocal (n = 5), diffuse (n = 5) and no histological lesions (n = 5). Significant differences were observed in the mean number of ITLN2 immunostained goblet and Paneth cells between the three histopathological types and the control. The number of immunolabelled cells was higher in the focal histopathological type (116.9 ± 113.9) followed by the multifocal (108.7 ± 140.5), diffuse (76.5 ± 97.8) and control types (41.0 ± 81.3). These results validate ITLN2 as a post-mortem biomarker of disease progression.

Identification of loci associated with susceptibility to bovine paratuberculosis and with the dysregulation of the MECOM, eEF1A2, and U1 spliceosomal RNA expression.

Although genome-wide association studies have identified single nucleotide polymorphisms (SNPs) associated with the susceptibility to Mycobacterium avium subsp. paratuberculosis (MAP) infection, only a few functional mutations for bovine paratuberculosis (PTB) have been characterized. Expression quantitative trait loci (eQTLs) are genetic variants typically located in gene regulatory regions that alter gene expression in an allele-specific manner. eQTLs can be considered as functional links between genomic variants, gene expression, and ultimately phenotype. In the current study, peripheral blood (PB) and ileocecal valve (ICV) gene expression was quantified by RNA-Seq from fourteen Holstein cattle with no lesions and with PTB-associated histopathological lesions in gut tissues. Genotypes were generated from the Illumina LD EuroG10K BeadChip. The associations between gene expression levels (normalized read counts) and genetic variants were analyzed by a linear regression analysis using R Matrix eQTL 2.2. This approach allowed the identification of 192 and 48 cis-eQTLs associated with the expression of 145 and 43 genes in the PB and ICV samples, respectively. To investigate potential relationships between these cis-eQTLs and MAP infection, a case-control study was performed using the genotypes for all the identified cis-eQTLs and phenotypical data (histopathology, ELISA for MAP-antibodies detection, tissue PCR, and bacteriological culture) of 986 culled cows. Our results suggested that the heterozygous genotype in the cis-eQTL-rs43744169 (T/C) was associated with the up-regulation of the MDS1 and EVI1 complex (MECOM) expression, with positive ELISA, PCR, and bacteriological culture results, and with increased risk of progression to clinical PTB. As supporting evidence, the presence of the minor allele was associated with higher MECOM levels in plasma samples from infected cows and with increased MAP survival in an ex-vivo macrophage killing assay. Moreover, the presence of the two minor alleles in the cis-eQTL-rs110345285 (C/C) was associated with the dysregulation of the eukaryotic elongation factor 1-α2 (eEF1A2) expression and with increased ELISA (OD) values. Finally, the presence of the minor allele in the cis-eQTL rs109859270 (C/T) was associated with the up-regulation of the U1 spliceosomal RNA expression and with an increased risk of progression to clinical PTB. The introduction of these novel functional variants into marker-assisted breeding programs is expected to have a relevant effect on PTB control.

Detection of latent forms of Mycobacterium avium subsp. paratuberculosis infection using host biomarker-based ELISAs greatly improves paratuberculosis diagnostic sensitivity.

Bovine paratuberculosis (PTB) is a chronic granulomatous enteritis, caused by Mycobacterium avium subsp. paratuberculosis (MAP), responsible for important economic losses in the dairy industry. Current diagnostic methods have low sensitivities for detection of latent forms of MAP infection, defined by focal granulomatous lesions and scarce humoral response or MAP presence. In contrast, patent infections correspond to multifocal and diffuse types of enteritis where there is increased antibody production, and substantial mycobacterial load. Our previous RNA-Seq analysis allowed the selection of five candidate biomarkers overexpressed in peripheral blood of MAP infected Holstein cows with focal (ABCA13 and MMP8) and diffuse (FAM84A, SPARC and DES) lesions vs. control animals with no detectable PTB-associated lesions in intestine and regional lymph nodes. The aim of the current study was to assess the PTB diagnostic potential of commercial ELISAs designed for the specific detection of these biomarkers. The ability of these ELISAs to identify animals with latent and/or patent forms of MAP infection was investigated using serum from naturally infected cattle (n = 88) and non-infected control animals (n = 67). ROC analysis revealed that the ABCA13-based ELISA showed the highest diagnostic accuracy for the detection of infected animals with focal lesions (AUC 0.837, sensitivity 79.25% and specificity 88.06%) and with any type of histological lesion (AUC 0.793, sensitivity 69.41% and specificity 86.57%) improving on the diagnostic performance of the popular IDEXX ELISA and other conventional diagnostic methods. SPARC and MMP8 showed the highest diagnostic accuracy for the detection of animals with multifocal (AUC 0.852) and diffuse lesions (AUC 0.831), respectively. In conclusion, our results suggest that quantification of ABCA13, SPARC and MMP8 by ELISA has the potential for implementation as a diagnostic tool to reliably identify MAP infection, greatly improving early detection of MAP latent infections when antibody responses and fecal shedding are undetectable using conventional diagnostic methods.

Local Lung Immune Response to Mycobacterium bovis Challenge after BCG and M. bovis Heat-Inactivated Vaccination in European Badger (Meles meles).

Tuberculosis (TB) vaccination could be used as a key part of integrated strategies for the disease's control if an effective and safe vaccine under field conditions is obtained. Recent studies in Spain have evaluated the protective efficacy of two oral vaccines against experimental challenge with live intra-bronchial Mycobacterium bovis in captive badgers: the live-attenuated M. bovis BCG vaccine (Danish strain) and a heat-inactivated M. bovis (HIMB) vaccine. With the objective of increasing the knowledge of the cellular development progress of infection and generating further tools to discriminate between mild and severe TB lesions between and within animals, the immunopathology of tuberculous lesions was studied to characterize the local immune response (cell type profile) within lung granulomas from control (non-vaccinated), BCG vaccinated and HIMB-vaccinated experimentally infected badgers with M. bovis. Four immunohistochemical protocols, for the specific detection of macrophages, T lymphocytes, B lymphocytes and plasma cells within TB granulomas in formalin fixed sections of the right middle lung lobe (lobe targeted for the M. bovis delivery), were performed. Immunolabelled sections were scanned and five randomly selected areas were analyzed with digital image analysis software. The results were expressed as the proportion of the positively immunolabelled area within the total area of the selected site. Data was analyzed using the statistical analysis software (SAS). In the three treatment groups, macrophages were the most abundant inflammatory cells within the granulomas, followed by B lymphocytes and plasma cells. T lymphocyes were absent in those granulomas. This would suggest a predominance of a non-specific innate response mediated by phagocytic cells over an adaptative humoral immune response. The proportion of macrophages and plasma cells was higher in BCG and HIMB-vaccinated badgers, respectively, suggesting the establishment of an adaptative humoral response in HIMB-vaccinated badgers. The lower bacterial load at the lung level, as well as the volume of lesions in lungs using magnetic resonance imaging in badgers with the HIMB vaccine in relation with local immune response presented, must be highlighted, since it would be an advantage in favor of its use under field conditions in terms of reducing TB transmission and environmental contamination.

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.

Myxoma virus jumps species to the Iberian hare.

The study of myxoma virus (MYXV) infections in the European rabbit (Oryctolagus cuniculus) has produced one of the most accepted host-pathogen evolutionary models. To date, myxomatosis has been limited to the European rabbit with sporadic reports in hares. However, reports of widespread mortalities in the Iberian hare (Lepus granatensis) with myxomatosis-like clinical signs indicate a potential species jump has occurred. The presence of MYXV DNA was confirmed by PCR in 244 samples received from regional veterinary services, animal health laboratories, hunters or rangers over a 5-month period. PCR analysis of 4 MYXV positive hare samples revealed a 2.8 kb insertion located within the M009 gene with respect to MYXV. The presence of this insertion was subsequently confirmed in 20 samples from 18 Spanish provinces. Sanger sequencing and subsequent analysis show that the insert contained 4 ORFs which are phylogenetically related to MYXV genes M060, M061, M064 and M065. The complete MYXV genome from hare tissue was sequenced using Ion torrent next-generation technology and a summary of the data presented here. With the exception of the inserted region, the virus genome had no large scale modifications and 110 mutations with respect to the MYXV reference strain Lausanne were observed. The next phase in the evolution of MYXV has taken place as a host species jump from the European rabbit to the Iberian hare an occurrence which could have important effects on this naïve population.

Vaccination of rabbits with immunodominant antigens from Sarcoptes scabiei induced high levels of humoral responses and pro-inflammatory cytokines but confers limited protection.

BACKGROUND: Vaccination is an attractive ecological alternative to the use of acaricides for parasite control. However, effective anti-parasite vaccines against sarcoptic mange have not yet been developed. The purpose of this study was first to identify Sarcoptes scabiei immunodominant antigens and second to evaluate them as vaccine candidates in a rabbit/S. scabiei var. cuniculi model. METHODS: The S. scabiei Ssλ15 immunodominant antigen was selected by immunoscreening of a S. scabiei var. hominis cDNA. The full-length cDNA was sequenced and cloned into the pGEX vector and the recombinant protein expressed in BL21 (DE3) cells and purified. A vaccination trial was performed consisting of a test group (n = 8) immunised with recAgs (a mix of two recombinant antigens, Ssλ15 and the previously described Ssλ20∆B3) and a control group (n = 8) immunised with PBS. All analyses were performed with R Statistical Environment with α set at 0.050. RESULTS: The full-length open reading frame of the 1,821 nt cloned cDNA encodes a 64 kDa polypeptide, the sequence of which had 96 % identity with a hypothetical protein of S. scabiei. Ssλ15 was localised by immunostaining of skin sections in the tegument surrounding the mouthparts and the coxa in the legs of mites. Rabbit immunisation with recAgs induced high levels of specific IgG (P < 0.010) and increased levels of total IgEs. However, no significant clinical protection against S. scabiei challenge was detected. Unexpectedly, the group immunised with the recAgs mix had significantly higher lesion scores (P = 0.050) although lower mean mite densities than those observed in the control group. These results might indicate that the lesions in the recAgs group were due not only to the mites density but also to an exacerbated immunological response after challenge, which is in agreement with the specific high levels of pro-inflammatory cytokines (IL-1 and TNFα) detected after challenge in this group. CONCLUSIONS: The selected antigens delivered as recombinant proteins had no clinical protective efficacy against S. scabiei infestation although immunisation reduced mite density. However, these results pave the way for future studies on alternative production systems, adjuvants, delivery methods and combinations of antigens in order to manage stimulation of clinical protective immune responses.

Influence of age and gender in the sensory analysis of balance control.

Postural control is achieved through the integration at the central nervous system level of information obtained by the visual, somatosensory and vestibular systems. Computerized dynamic posturography and the Sway Star system are both used to carry out sensory analysis. The purpose of this study was to determine the influence of sex and age on sensory analysis, measured with these two systems, and to compare their results. A prospective trial was conducted with 70 healthy individuals (average age: 44.9 years) uniformly distributed in seven age groups, who underwent postural study with both systems. We used SPSS 16.0 for statistical study: comparison of means test for influence of gender and age and Pearson's correlation test (p < 0.05). Gender variable had no influence. The influence of age in vestibular input was found to be significant with both posturography systems, while visual input was only found to be significant with the Sway Star. The results with the two systems were not comparable. Sensory contribution does not remain stable throughout life. Visual information decreases with age, reaching a minimum at 40-49 years, and may correspond to the deterioration of eyesight with age. Propioceptive information showed no statistically significant changes, and several forms of treatment might correct the deterioration of this system. Vestibular information reaches a maximum in the 40-49 years age group in an attempt to compensate for visual deterioration, and decreases again in subsequent decades. This may be due to aging of the vestibular system and the difficulty in its correction.

Temporal stability in the genetic structure of Sarcoptes scabiei under the host-taxon law: empirical evidences from wildlife-derived Sarcoptes mite in Asturias, Spain.

BACKGROUND: Implicitly, parasite molecular studies assume temporal genetic stability. In this study we tested, for the first time to our knowledge, the extent of changes in genetic diversity and structure of Sarcoptes mite populations from Pyrenean chamois (Rupicapra pyrenaica) in Asturias (Spain), using one multiplex of 9 microsatellite markers and Sarcoptes samples from sympatric Pyrenean chamois, red deer (Cervus elaphus), roe deer (Capreolus capreolus) and red fox (Vulpes vulpes). RESULTS: The analysis of an 11-years interval period found little change in the genetic diversity (allelic diversity, and observed and expected heterozygosity). The temporal stability in the genetic diversity was confirmed by population structure analysis, which was not significantly variable over time. Population structure analysis revealed temporal stability in the genetic diversity of Sarcoptes mite under the host-taxon law (herbivore derived- and carnivore derived-Sarcoptes mite) among the sympatric wild animals from Asturias. CONCLUSIONS: The confirmation of parasite temporal genetic stability is of vital interest to allow generalizations to be made, which have further implications regarding the genetic structure, epidemiology and monitoring protocols of the ubiquitous Sarcoptes mite. This could eventually be applied to other parasite species.

Identification and heterologous expression of a Sarcoptes scabiei cDNA encoding a structural antigen with immunodiagnostic potential.

The mite Sarcoptes scabiei causes sarcoptic mange (or scabies), a disease of considerable human and veterinary significance. An S. scabiei cDNA clone of about 2 kb was isolated from a S. scabiei var. hominis expression library by immunological screening using blood serum from a naturally infected chamois (Rupicapra rupicapra). The nucleotide sequence of the identified cDNA contains an open reading frame of 1930 bp that encodes a 642 amino acid polypeptide. This polypeptide shows tandem repeats of a glycine-serine rich 20 residue sequence followed by a unique C-terminal glutamate rich 54 residue sequence. The cDNA or the deduced polypeptide did not show significant similarities to any of the sequences in the databases. A carboxyl-terminal fragment of this polypeptide (residues 380 to 642) was efficiently expressed in Escherichia coli as a fusion with Glutathione S-transferase and then was used to produce a specific antiserum. The antigen encoded by the cDNA was located at the integument of the mite's epidermis and the cavities surrounding its vital organs. Western blot analysis of mite extracts using the specific antiserum against the recombinant protein identified antigens larger that 60 kDa indicating that the isolated cDNA did not contain the full ORF. Moreover, we designed a diagnostic assay based on the carboxyl-terminal fragment of the antigen for the identification of infected animals.

Generation of a recombinant avian coronavirus infectious bronchitis virus using transient dominant selection.

A reverse genetics system for the avian coronavirus infectious bronchitis virus (IBV) has been described in which a full-length cDNA, corresponding to the IBV (Beaudette-CK) genome, was inserted into the vaccinia virus genome following in vitro assembly of three contiguous cDNAs [Casais, R., Thiel, V., Siddell, S.G., Cavanagh, D., Britton, P., 2001. Reverse genetics system for the avian coronavirus infectious bronchitis virus. J. Virol. 75, 12359-12369]. The method has subsequently been used to generate a recombinant IBV expressing a chimaeric S gene [Casais, R., Dove, B., Cavanagh, D., Britton, P., 2003. Recombinant avian infectious bronchitis virus expressing a heterologous spike gene demonstrates that the spike protein is a determinant of cell tropism. J. Virol. 77, 9084-9089]. Use of vaccinia virus as a vector for the full-length cDNA of the IBV genome has the advantage that modifications can be made to the IBV cDNA using homologous recombination, a method frequently used to insert and delete sequences from the vaccinia virus genome. We describe the use of homologous recombination as a method for modifying the Beaudette full-length cDNA, within the vaccinia virus genome, without the requirement for in vitro assembly of the IBV cDNA. To demonstrate the feasibility of the method we exchanged the ectodomain of the Beaudette spike gene for the corresponding region from IBV M41 and generated two recombinant infectious bronchitis viruses (rIBVs) expressing the chimaeric S protein, validating the method as an alternative way for generating rIBVs.

Recombinant infectious bronchitis coronavirus Beaudette with the spike protein gene of the pathogenic M41 strain remains attenuated but induces protective immunity.

We have replaced the ectodomain of the spike (S) protein of the Beaudette strain (Beau-R; apathogenic for Gallus domesticus chickens) of avian infectious bronchitis coronavirus (IBV) with that from the pathogenic M41 strain to produce recombinant IBV BeauR-M41(S). We have previously shown that this changed the tropism of the virus in vitro (R. Casais, B. Dove, D. Cavanagh, and P. Britton, J. Virol. 77:9084-9089, 2003). Herein we have assessed the pathogenicity and immunogenicity of BeauR-M41(S). There were no consistent differences in pathogenicity between the recombinant BeauR-M41(S) and its apathogenic parent Beau-R (based on snicking, nasal discharge, wheezing, watery eyes, rales, and ciliostasis in trachea), and both replicated poorly in trachea and nose compared to M41; the S protein from the pathogenic M41 had not altered the apathogenic nature of Beau-R. Both Beau-R and BeauR-M41(S) induced protection against challenge with M41 as assessed by absence of recovery of challenge virus and nasal exudate. With regard to snicking and ciliostasis, BeauR-M41(S) induced greater protection (seven out of nine chicks [77%]; assessed by ciliostasis) than Beau-R (one out of nine; 11%) but less than M41 (100%). The greater protection induced by BeauR-M41(S) against M41 may be related to the ectodomain of the spike protein of Beau-R differing from that of M41 by 4.1%; a small number of epitopes on the S protein may play a disproportionate role in the induction of immunity. The results are promising for the prospects of S-gene exchange for IBV vaccine development.

Recombinant avian infectious bronchitis virus expressing a heterologous spike gene demonstrates that the spike protein is a determinant of cell tropism.

A recombinant infectious bronchitis virus (IBV), BeauR-M41(S), was generated using our reverse genetics system (R. Casais, V. Thiel, S. G. Siddell, D. Cavanagh, and P. Britton, J. Virol. 75:12359-12369, 2001), in which the ectodomain region of the spike gene from IBV M41-CK replaced the corresponding region of the IBV Beaudette genome. BeauR-M41(S) acquired the same cell tropism phenotype as IBV M41-CK in four different cell types, demonstrating that the IBV spike glycoprotein is a determinant of cell tropism.