Evidence for cattle major histocompatibility complex (BoLA) class IIDQA1gene heterozygote advantage against clinical mastitis caused byStreptococciandEscherichiaspecies

Genetic factors underlying host resistance to Rhipicephalus microplus tick infestation in Braford cattle: a systems biology perspective

Approximately 80% of the world's cattle are raised in regions with a high risk of tick-borne diseases, resulting in significant economic losses due to parasitism by Rhipicephalus (Boophilus) microplus. However, the lack of a systemic biology approach hampers a comprehensive understanding of tick-host interactions that mediate tick resistance phenotypes. Here, we conducted a genome-wide association study (GWAS) of 2933 Braford cattle and found 340 single-nucleotide polymorphisms (SNPs) associated with tick counts. Gene expression analyses were performed on skin samples obtained from previously tick-exposed heifers with extremely high or low estimated breeding values for R. microplus counts. Evaluations were performed both before and after artificial infestation with ticks. Differentially expressed genes were found within 1-Mb windows centered at significant SNPs from GWAS. A total of 330 genes were related to the breakdown of homeostasis that was induced by larval attachment to bovine skin. Enrichment analysis pointed to a key role of proteolysis and signal transduction via JAK/STAT, NFKB and WNT/beta catenin signaling pathways. Integrative analysis on matrixEQTL revealed two cis-eQTLs and four significant SNPs in the genes peptidyl arginine deiminase type IV (PADI4) and LOC11449251. The integration of genomic data from QTL maps and transcriptome analyses has identified a set of twelve key genes that show significant associations with tick loads. These genes could be key candidates to improve the accuracy of genomic predictions for tick resistance in Braford cattle.

BoLA-DRB3*15:01 allele is associated with susceptibility to early enzootic bovine leukosis onset in Holstein-Friesian and Japanese Black cattle

Enzootic bovine leukosis (EBL) is typically observed in cattle older than 3 years, but some cases of onset in cattle younger than 3 years have been reported in Japan. BoLA-DRB3 polymorphisms are associated with susceptibility to EBL onset. However, little is known about the relationship between the polymorphisms and EBL onset in young cattle. In the present study, we performed BoLA-DRB3 genotyping in 59 EBL cattle younger than 3 years (25 Holstein-Friesian and 34 Japanese Black) and compared the results with those of 69 EBL cattle older than 3 years (38 Holstein-Friesian and 31 Japanese Black). The BoLA-DRB3*15:01 allele was detected at a frequency of 37.3 % (48.0 % and 29.4 % in Holstein-Friesian and Japanese Black, respectively) and was identified as an early EBL onset susceptibility allele. Nine EBL cattle younger than 3 years (5 Holstein-Friesian and 4 Japanese Black), but only 1 EBL cattle older than 3 years (1 Holstein-Friesian), had a BoLA-DRB3*15:01/*15:01 homozygous genotype. The frequency of the BoLA-DRB3*15:01 allele occurring with a different allele (BoLA-DRB3*015:01/other) in cattle younger than 3 years was 44.1 % (56.0 % Holstein-Friesian and 35.3 % Japanese Black) and significantly higher than that in cattle older than 3 years (28.9 % Holstein-Friesian and 9.7 % Japanese Black) (P = 0.0013). These results suggest that BoLA-DRB3*15:01/*15:01 and BoLA-DRB3*15:01/other genotypes are early EBL onset susceptibility genotypes. The present findings may contribute to cattle breeding selection.

Influence of BoLA-DRB3 Polymorphism and Bovine Leukemia Virus (BLV) Infection on Dairy Cattle Productivity

Enzootic bovine leukosis caused by the bovine leukemia virus (BLV) results in substantial damage to the livestock industry; however, we lack an effective cure or vaccine. BoLA-DRB3 polymorphism in BLV-infected cattle is associated with the proviral load (PVL), infectivity in the blood, development of lymphoma, and in utero infection of calves. Additionally, it is related to the PVL, infectivity, and anti-BLV antibody levels in milk. However, the effects of the BoLA-DRB3 allele and BLV infection on dairy cattle productivity remain poorly understood. Therefore, we investigated the effect of BLV infection and BoLA-DRB3 allele polymorphism on dairy cattle productivity in 147 Holstein dams raised on Japanese dairy farms. Our findings suggested that BLV infection significantly increased milk yield. Furthermore, the BoLA-DRB3 allele alone, and the combined effect of BLV infection and the BoLA-DRB3 allele had no effect. These results indicate that on-farm breeding and selection of resistant cattle, or the preferential elimination of susceptible cattle, does not affect dairy cattle productivity. Additionally, BLV infection is more likely to affect dairy cattle productivity than BoLA-DRB3 polymorphism.

BoLA-DRB3 Polymorphism Controls Proviral Load and Infectivity of Bovine Leukemia Virus (BLV) in Milk

Bovine leukemia virus (BLV), which causes enzootic bovine leukosis, is transmitted to calves through the milk of BLV-infected dams. Bovine leukocyte antigen (BoLA)-DRB3 is a polymorphic gene associated with BLV infectivity and proviral load (PVL). However, the effect of BoLA-DRB3 polymorphism on the infectivity and PVL of milk from BLV-infected dams remains unknown. This study examined milk from 259 BLV-infected dams, including susceptible dams carrying at least one BoLA-DRB3*012:01 or *015:01 allele with high PVL, resistant dams carrying at least one BoLA-DRB3*002:01, *009:02, or *014:01:01 allele with low PVL, and neutral dams carrying other alleles. The detection rate of BLV provirus and PVL were significantly higher in milk from susceptible dams than in that from resistant dams. This result was confirmed in a three-year follow-up study in which milk from susceptible dams showed a higher BLV provirus detection rate over a longer period than that from resistant dams. The visualization of infectivity of milk cells using a luminescence syncytium induction assay showed that the infectious risk of milk from BLV-infected dams was markedly high for susceptible dams compared to resistant ones. This is the first report confirming that BoLA-DRB3 polymorphism affects the PVL and infectivity of milk from BLV-infected dams.

Comparing Class II MHC DRB3 Diversity in Colombian Simmental and Simbrah Cattle Across Worldwide Bovine Populations

The major histocompatibility complex (MHC) exerts great influence on responses to infectious diseases and vaccination due to its fundamental role in the adaptive immune system. Knowledge about MHC polymorphism distribution among breeds can provide insights into cattle evolution and diversification as well as population-based immune response variability, thus guiding further studies. Colombian Simmental and Simbrah cattle’s BoLA-DRB3 genetic diversity was compared to that of taurine and zebuine breeds worldwide to estimate functional diversity. High allele richness was observed for Simmental and Simbrah cattle; nevertheless, high homozygosity was associated with individual low sequence variability in both the β1 domain and the peptide binding region (PBR), thereby implying reduced MHC-presented peptide repertoire size. There were strong signals of positive selection acting on BoLA-DRB3 in all populations, some of which were poorly structured and displayed common alleles accounting for their high genetic similarity. PBR sequence correlation analysis suggested that, except for a few populations exhibiting some divergence at PBR, global diversity regarding potential MHC-presented peptide repertoire could be similar for the cattle populations analyzed here, which points to the retention of functional diversity in spite of the selective pressures imposed by breeding.

Diversity of the BoLA-DRB3 gene in cattle breeds from tropical and subtropical regions of Argentina

Bovine leukocyte antigens (BoLA) have been widely studied because of their primary function in the recognition of pathogens by the immune system. To date, however, the characterization of the BoLA-DRB3 gene in Latin American Zebu and mixed zebuine breeds is scarce. By a sequence-based typing method, here we sequenced exon 2 of BoLA class II DRB3 gene in 264 animals from the five most commonly used breeds in northern Argentina (Creole, Brahman, Braford, Brangus, and Nellore).The Bos taurus, Bos indicus, and mixed breeds analyzed here contained 61 previously reported alleles. Genetic diversity was high at both allelic and nucleotide sequence levels, particularly in the mixed breeds Braford and Brangus. In contrast to previous reports on DRB3 diversity, no evidence of balancing selection was found in our data. Differentiation among breeds was highly significant, as shown by FST (F_ST = 0.052, P < 0.001) and cluster analyses. In accordance with historical origin of the breeds, UPGMA trees and metric multidimensional scaling (MDS) analyses showed that Creole is distantly related to the other zebuine breeds. Among them, Brahman, Braford, and Brangus exhibited the closest affiliations. Despite the overall differentiation of the breeds, analysis of the peptide binding regions at the aminoacid level revealed that the key aminoacids involved in peptide recognition are greatly conserved suggesting little influence of domestication and breeding in functional MHC variability. In sum, this is the first report of BoLA-DRB3 diversity in pure and mixed Bos indicus cattle breeds from Argentina. Knowledge of BoLA-DRB3 variability in breeds adapted to tropical and subtropical environments contributes not only to the characterization of MHC diversity but also to the design of peptide-based vaccines.

Characterization of BoLA class II DQA and DQB by PCR-RFLP, cloning, and sequencing reveals sequence diversity in crossbred cattle

The BoLA class II DQA and DQB genes in crossbred cattle were studied using PCR-RFLP, cloning, and sequencing techniques. Seventy-two crossbred cattle (Vrindavani) were used in the current study. HaeIII and XbaI restriction enzymes digested DQA exon 2-3, revealing seven (HaeIII-A-G) and three (XbaI A-C) motifs, respectively. The BoLA-DQB gene was analyzed using PCR-RFLP with PstI and TaqI restriction enzymes, yielding five restriction motifs for each restriction enzyme (PstI-A-E and TaqI-A-E). In crossbred cattle, addition, deletion, and substitutions were observed in distinct sequences, resulting in variations in overall gene length. Changes in nucleotides at positions 64-80, 110-200, and 207-264 were largely responsible for polymorphism in DQA exon 2. The phylogenetic analysis predicted a high degree of nucleotide and amino acid changes in DQA exon 2-3 and DQB exon 2. DQA genes had a nucleotide dissimilarity of 0.3-25.4 percent, while DQB genes had a nucleotide dissimilarity of 1.5-14.3 percent. We cloned and sequenced 20 genotypes based on PCR-RFLP of the DQA and DQB genes. The current study observed variation in the DQA and DQB genes and will serve as a foundation for future research on the BoLA DQA and DQB genes.

Integrative Systems Biology Analysis Elucidates Mastitis Disease Underlying Functional Modules in Dairy Cattle

Background: Mastitis is the most prevalent disease in dairy cattle and one of the most significant bovine pathologies affecting milk production, animal health, and reproduction. In addition, mastitis is the most common, expensive, and contagious infection in the dairy industry. Methods: A meta-analysis of microarray and RNA-seq data was conducted to identify candidate genes and functional modules associated with mastitis disease. The results were then applied to systems biology analysis via weighted gene coexpression network analysis (WGCNA), Gene Ontology, enrichment analysis for the Kyoto Encyclopedia of Genes and Genomes (KEGG), and modeling using machine-learning algorithms. Results: Microarray and RNA-seq datasets were generated for 2,089 and 2,794 meta-genes, respectively. Between microarray and RNA-seq datasets, a total of 360 meta-genes were found that were significantly enriched as "peroxisome," "NOD-like receptor signaling pathway," "IL-17 signaling pathway," and "TNF signaling pathway" KEGG pathways. The turquoise module (n = 214 genes) and the brown module (n = 57 genes) were identified as critical functional modules associated with mastitis through WGCNA. PRDX5, RAB5C, ACTN4, SLC25A16, MAPK6, CD53, NCKAP1L, ARHGEF2, COL9A1, and PTPRC genes were detected as hub genes in identified functional modules. Finally, using attribute weighting and machine-learning methods, hub genes that are sufficiently informative in Escherichia coli mastitis were used to optimize predictive models. The constructed model proposed the optimal approach for the meta-genes and validated several high-ranked genes as biomarkers for E. coli mastitis using the decision tree (DT) method. Conclusion: The candidate genes and pathways proposed in this study may shed new light on the underlying molecular mechanisms of mastitis disease and suggest new approaches for diagnosing and treating E. coli mastitis in dairy cattle.

BoLA-DRB3 gene haplotypes show divergence in native Sudanese cattle from taurine and indicine breeds

Autochthonous Sudanese cattle breeds, namely Baggara for beef and Butana and Kenana for dairy, are characterized by their adaptive characteristics and high performance in hot and dry agro-ecosystems. They are thus used largely by nomadic and semi-nomadic pastoralists. We analyzed the diversity and genetic structure of the BoLA-DRB3 gene, a genetic locus linked to the immune response, for the indigenous cattle of Sudan and in the context of the global cattle repository. Blood samples (n = 225) were taken from three indigenous breeds (Baggara; n = 113, Butana; n = 60 and Kenana; n = 52) distributed across six regions of Sudan. Nucleotide sequences were genotyped using the sequence-based typing method. We describe 53 alleles, including seven novel alleles. Principal component analysis (PCA) of the protein pockets implicated in the antigen-binding function of the MHC complex revealed that pockets 4 and 9 (respectively) differentiate Kenana-Baggara and Kenana-Butana breeds from other breeds. Venn analysis of Sudanese, Southeast Asian, European and American cattle breeds with 115 alleles showed 14 were unique to Sudanese breeds. Gene frequency distributions of Baggara cattle showed an even distribution suggesting balancing selection, while the selection index (ω) revealed the presence of diversifying selection in several amino acid sites along the BoLA-DRB3 exon 2 of these native breeds. The results of several PCA were in agreement with clustering patterns observed on the neighbor joining (NJ) trees. These results provide insight into their high survival rate for different tropical diseases and their reproductive capacity in Sudan's harsh environment.

MHC Haplotyping of SARS-CoV-2 Patients: HLA Subtypes Are Not Associated with the Presence and Severity of COVID-19 in the Israeli Population

HLA haplotypes were found to be associated with increased risk for viral infections or disease severity in various diseases, including SARS. Several genetic variants are associated with COVID-19 severity. Studies have proposed associations, based on a very small sample and a large number of tested HLA alleles, but no clear association between HLA and COVID-19 incidence or severity has been reported. We conducted a large-scale HLA analysis of Israeli individuals who tested positive for SARS-CoV-2 infection by PCR. Overall, 72,912 individuals with known HLA haplotypes were included in the study, of whom 6413 (8.8%) were found to have SARS-CoV-2 by PCR. A total of 20,937 subjects were of Ashkenazi origin (at least 2/4 grandparents). One hundred eighty-one patients (2.8% of the infected) were hospitalized due to the disease. None of the 66 most common HLA loci (within the five HLA subgroups: A, B, C, DQB1, DRB1) was found to be associated with SARS-CoV-2 infection or hospitalization in the general Israeli population. Similarly, no association was detected in the Ashkenazi Jewish subset. Moreover, no association was found between heterozygosity in any of the HLA loci and either infection or hospitalization. We conclude that HLA haplotypes are not a major risk/protecting factor among the Israeli population for SARS-CoV-2 infection or severity. Our results suggest that if any HLA association exists with the disease it is very weak, and of limited effect on the pandemic.

Nucleotide sequence variation of the major histocompatibility complex class II DQA1 gene in different cattle breeds from Nigeria and New Zealand

The major histocompatibility complex (MHC) plays a role in immune response. Among other activities, the bovine MHC genes (BoLA) trigger immune responses, including the activation of antibody-producing B-cells. In this study, White Fulani (n = 24), Red Bororo (n = 5) and Holstein-White × Fulani-cross (n = 11) cattle from Nigeria, and New Zealand Holstein-Friesian × Jersey-cross (n = 40) cattle were used to investigate variability in exon 2 of BoLA-DQA1. Ten alleles were identified using a PCR-Single Strand Conformation Polymorphism (SSCP) approach and their nucleotide sequences confirmed by DNA sequencing. A total of 12.60 % of all nucleotide positions analysed were revealed to be variable and two novel BoLA-DQA1 alleles are reported here for the first time.

Bovine major histocompatibility complex (BoLA) heterozygote advantage against the outcome of bovine leukemia virus infection

Bovine leukemia virus (BLV) causes enzootic bovine leucosis. Host genetic heterozygosity at the major histocompatibility complex can enhance the ability to combat infectious diseases. However, heterozygote advantage is loci specific and depends on disease type. Bovine leukocyte antigen (BoLA)-DRB3 polymorphisms are related with BLV-infection outcome; however, whether BoLA-DRB3 heterozygotes have an advantage against BLV-induced lymphoma and proviral load (PVL) remains unclear. By analyzing 1567 BLV-infected individuals, we found that BoLA-DRB3 heterozygous status was significantly associated with lymphoma resistance irrespective of cattle breeds (p < 0.0001). Similarly, decreased PVL was observed in BoLA-DRB3 heterozygotes (p = 0.0407 for Holstein cows; p = 0.0889 for Japanese Black cattle). Our report provides first evidence of BoLA-DRB3 heterozygote advantage against BLV infection outcome.

Characterization of bovine MHC DRB3 diversity in global cattle breeds, with a focus on cattle in Myanmar

Background

Myanmar cattle populations predominantly consist of native cattle breeds (Pyer Sein and Shwe), characterized by their geographical location and coat color, and the Holstein-Friesian crossbreed, which is highly adapted to the harsh tropical climates of this region. Here, we analyzed the diversity and genetic structure of the BoLA-DRB3 gene, a genetic locus that has been linked to the immune response, in Myanmar cattle populations.

Methods

Blood samples (n = 294) were taken from two native breeds (Pyer Sein, n = 163 and Shwe Ni, n = 69) and a cattle crossbreed (Holstein-Friesian, n = 62) distributed across six regions of Myanmar (Bago, n = 38; Sagaing, n = 77; Mandalay, n = 46; Magway, n = 46; Kayin, n = 43; Yangon, n = 44). In addition, a database that included 2428 BoLA-DRB3 genotypes from European (Angus, Hereford, Holstein, Shorthorn, Overo Negro, Overo Colorado, and Jersey), Zebuine (Nellore, Brahman and Gir), Asian Native from Japan and Philippine and Latin-American Creole breeds was also included. Furthermore, the information from the IPD–MHC database was also used in the present analysis. DNA was genotyped using the sequence-based typing method. DNA electropherograms were analyzed using the Assign 400ATF software.

Results

We detected 71 distinct alleles, including three new variants for the BoLA-DRB3 gene. Venn analysis showed that 11 of these alleles were only detected in Myanmar native breeds and 26 were only shared with Asian native and/or Zebu groups. The number of alleles ranged from 33 in Holstein-Friesians to 58 in Pyer Seins, and the observed versus unbiased expected heterozygosity were higher than 0.84 in all the three the populations analyzed. The F_ST analysis showed a low level of genetic differentiation between the two Myanmar native breeds (F_ST = 0.003), and between these native breeds and the Holstein-Friesians (F_ST <  0.021). The average F_ST value for all the Myanmar Holstein-Friesian crossbred and Myanmar native populations was 0.0136 and 0.0121, respectively. Principal component analysis (PCA) and tree analysis showed that Myanmar native populations grouped in a narrow cluster that diverged clearly from the Holstein-Friesian populations. Furthermore, the BoLA-DRB3 allele frequencies suggested that while some Myanmar native populations from Bago, Mandalay and Yangon regions were more closely related to Zebu breeds (Gir and Brahman), populations from Kayin, Magway and Sagaing regions were more related to the Philippines native breeds. On the contrary, PCA showed that the Holstein-Friesian populations demonstrated a high degree of dispersion, which is likely the result of the different degrees of native admixture in these populations.

Conclusion

This study is the first to report the genetic diversity of the BoLA-DRB3 gene in two native breeds and one exotic cattle crossbreed from Myanmar. The results obtained contribute to our understanding of the genetic diversity and distribution of BoLA-DRB3 gene alleles in Myanmar, and increases our knowledge of the worldwide variability of cattle BoLA-DRB3 genes, an important locus for immune response and protection against pathogens.

Diversity analysis at MHC class II DQA locus in buffalo (Bubalus bubalis) indicates extensive duplication and trans-species evolution

Variation at MHC Class II-DQA locus in riverine and swamp buffaloes (Bubu) has been explored in this study. Through sequencing of buffalo DQA, 48 nucleotide variants identified from 17 individuals, reporting 42 novel alleles, including one pseudogene. Individual animal displayed two to seven variants, suggesting the presence of more than two Bubu-DQA loci, as an evidence of extensive duplication. dN values were found to be higher than dS values at peptide binding sites, separately for riverine and swamp buffaloes, indicating locus being under positive selection. Evolutionary analysis revealed numerous trans-species polymorphism with alleles from water buffalo assigned to at least three different loci (Bubu-DQA1, DQA2, DQA3). Alleles of both the sub-species intermixed within the cluster, showing convergent evolution of MHC alleles in bovines. The results thus suggest that both riverine and swamp buffaloes share con-current arrangement of DQA region, comparable to cattle in terms of copy number and population polymorphism.

Whole genome characterization of autochthonous Bos taurus brachyceros and introduced Bos indicus indicus cattle breeds in Cameroon regarding their adaptive phenotypic traits and pathogen resistance

Background

African indigenous taurine cattle display unique adaptive traits shaped by husbandry management, regional climate and exposure to endemic pathogens. They are less productive with respect to milk and meat production which has been associated with amongst others, small size, traditional beliefs, husbandry practices, limited feed resources, disease burden and lack of sustained breeding for trait improvement. This resulted in the severe dwindling of their population size rendering them vulnerable to extinction.

The Namchi taurine cattle breed is referred to as [Namchi (Doayo)] and shows resistance traits against trypanosome infection and exposure to tick infestation. Nonetheless, the historically later introduced Zebu cattle are the main cattle breeds in Africa today, even though they suffer more from locally prevailing pathogens.

By using a whole genome sequencing approach, we sequenced with high depth for the first time the genomes of five cattle breeds from Cameroon in order to provide a valuable genetic resource for future African cattle breeding: the Namchi, an endangered trypano-tolerant taurine breed, the Kapsiki, an indigenous trypano-susceptible taurine breed, and three Zebu (Bos indicus indicus) breeds: Ngaoundere Gudali, White Fulani and Red Fulani.

Results

Approximately 167 Gigabases of raw sequencing data were generated for each breed and mapped to the cattle reference genomes ARS-UCD1.2 and UMD3.1.The coverage was 103 to 140-fold when aligning the reads to ARS-UCD1.2 with an average mapping rate of ~ 99%, and 22 to 30-fold when aligning the reads to UMD3.1 with an average mapping rate of ~ 64%. The single nucleotide polymorphisms (SNPs) obtained from analysis using the genome ARS-UCD1.2 were compared with reference genomes of European Bos taurus Holstein, the Asian Bos indicus Brahman, and the African trypanotolerant N’Dama breeds.

A total of ~ 100 million (M) SNPs were identified and 7.7 M of those were breed-specific. An approximately 11.1 M constituted of small insertions and deletions. By using only breed-specific non-synonymous variants we identified genes as genetic signatures and associated Gene Ontology (GO) terms that could explain certain cattle-breed specific phenotypes such as increased tolerance against trypanosome parasites in the Namchi breed and heat tolerance in the Kapsiki breed. Phylogenetic analysis grouped, except for Namchi, the Bos taurus breeds Kapsiki, N’Dama and Holstein together while the B. indicus breeds White and Red Fulani, Gudali and Brahman clustered separately. The deviating result for Namchi indicates a hybrid status of the selected animal with a recent introgression of Zebu genes into its genome.

Conclusions

The findings provide the first comprehensive set of genome-wide variant data of the most important Cameroonian cattle breeds. The genomic data shall constitute a foundation for breed amelioration whilst exploiting the heritable traits and support conservation efforts for the endangered local cattle breeds.

HLA Heterozygote Advantage against HIV-1 Is Driven by Quantitative and Qualitative Differences in HLA Allele-Specific Peptide Presentation

Pathogen-mediated balancing selection is regarded as a key driver of host immunogenetic diversity. A hallmark for balancing selection in humans is the heterozygote advantage at genes of the human leukocyte antigen (HLA), resulting in improved HIV-1 control. However, the actual mechanism of the observed heterozygote advantage is still elusive. HLA heterozygotes may present a broader array of antigenic viral peptides to immune cells, possibly resulting in a more efficient cytotoxic T-cell response. Alternatively, heterozygosity may simply increase the chance to carry the most protective HLA alleles, as individual HLA alleles are known to differ substantially in their association with HIV-1 control. Here, we used data from 6,311 HIV-1-infected individuals to explore the relative contribution of quantitative and qualitative aspects of peptide presentation in HLA heterozygote advantage against HIV. Screening the entire HIV-1 proteome, we observed that heterozygous individuals exhibited a broader array of HIV-1 peptides presented by their HLA class I alleles. In addition, viral load was negatively correlated with the breadth of the HIV-1 peptide repertoire bound by an individual's HLA variants, particularly at HLA-B. This suggests that heterozygote advantage at HLA-B is at least in part mediated by quantitative peptide presentation. We also observed higher HIV-1 sequence diversity among HLA-B heterozygous individuals, suggesting stronger evolutionary pressure from HLA heterozygosity. However, HLA heterozygotes were also more likely to carry certain HLA alleles, including the highly protective HLA-B*57:01 variant, indicating that HLA heterozygote advantage ultimately results from a combination of quantitative and qualitative effects in antigen presentation.

Major Histocompatibility Complex Class II (DRB3) Genetic Diversity in Spanish Morucha and Colombian Normande Cattle Compared to Taurine and Zebu Populations

Bovine leukocyte antigens (BoLA) have been used as disease markers and immunological traits in cattle due to their primary role in pathogen recognition by the immune system. A higher MHC allele diversity in a population will allow presenting a broader peptide repertoire. However, loss of overall diversity due to domestication process can decrease a population's peptide repertoire. Within the context of zebu and taurine cattle populations, BoLA-DRB3 genetic diversity in Spanish Morucha and Colombian Normande cattle was analyzed and an approach to estimate functional diversity was performed. Sequence-based typing was used for identifying 29, 23, 27, and 28 alleles in Spanish Morucha, Nariño-, Boyacá-, and Cundinamarca-Normande cattle, respectively. These breeds had remarkably low heterozygosity levels and the Hardy-Weinberg principle revealed significant heterozygote deficiency. F_ST and D_A genetic distance showed that Colombian Normande populations had greater variability than other phenotypically homogeneous breeds, such as Holstein. It was also found that Spanish Morucha cattle were strongly differentiated from other cattle breeds. Spanish Morucha had greater divergence in the peptide-binding region regarding other cattle breeds. However, peptide-binding region covariation indicated that the potential peptide repertoire seemed equivalent among cattle breeds. Despite the genetic divergence observed, the extent of the potential peptide repertoire in the cattle populations studied appears to be similar and thus their pathogen recognition potential should be equivalent, suggesting that functional diversity might persist in the face of bottlenecks imposed by domestication and breeding.

Probiotics and Prebiotics for the Amelioration of Type 1 Diabetes: Present and Future Perspectives

Type 1-diabetes (T1D) is an autoimmune disease characterized by immune-mediated destruction of pancreatic beta (β)-cells. Genetic and environmental interactions play an important role in immune system malfunction by priming an aggressive adaptive immune response against β-cells. The microbes inhabiting the human intestine closely interact with the enteric mucosal immune system. Gut microbiota colonization and immune system maturation occur in parallel during early years of life; hence, perturbations in the gut microbiota can impair the functions of immune cells and vice-versa. Abnormal gut microbiota perturbations (dysbiosis) are often detected in T1D subjects, particularly those diagnosed as multiple-autoantibody-positive as a result of an aggressive and adverse immunoresponse. The pathogenesis of T1D involves activation of self-reactive T-cells, resulting in the destruction of β-cells by CD8⁺ T-lymphocytes. It is also becoming clear that gut microbes interact closely with T-cells. The amelioration of gut dysbiosis using specific probiotics and prebiotics has been found to be associated with decline in the autoimmune response (with diminished inflammation) and gut integrity (through increased expression of tight-junction proteins in the intestinal epithelium). This review discusses the potential interactions between gut microbiota and immune mechanisms that are involved in the progression of T1D and contemplates the potential effects and prospects of gut microbiota modulators, including probiotic and prebiotic interventions, in the amelioration of T1D pathology, in both human and animal models.

Detection of quantitative trait loci and putative causal variants affecting somatic cell score in dairy sheep by using a 50K SNP chip and whole-genome sequencing

This study presents a scan of the ovine genome to identify quantitative trait loci (QTL) influencing the somatic cell score (SCS), a classical indicator of subclinical mastitis in sheep, and a subsequent high-resolution analysis of one of the identified QTL regions based on the analysis of whole-genome sequence data sets. A half-sib commercial population of Churra sheep genotyped with a 50K SNP chip was analyzed using linkage analysis (LA) and combined linkage and linkage disequilibrium analysis (LDLA). By LA, 2 5% chromosome-wide significant QTL on OAR5 and OAR25 and one 5% genome-wide significant QTL on ovine chromosome 20 (OAR20) were detected, whereas 22 significant associations were identified by LDLA. Two of the associations detected by LDLA replicated LA-detected effects (OAR20, OAR25). We compared the detected associations with previously reported QTL in sheep and cattle, and functional candidate genes were identified within the estimated confidence intervals. We then performed a high-resolution analysis of the OAR20 QTL region, the most significant QTL region identified by LA that replicated a QTL previously described in Churra sheep for SCS using microsatellite markers. For that, 2 segregating trios of 2 segregating families for the OAR20 QTL (each including the Qq sire and 2 daughters, QQ and qq) were selected for whole-genome sequencing. The bioinformatic analysis of the 6 sequenced samples performed across the genomic interval considered (14.2-41.7 Mb) identified a total of 227,030 variants commonly identified by 2 independent software packages. For the 3 different concordance tests considered, due to discrepancies regarding the QTL peak in the segregating families, the list of mutations concordant with the QTL segregating pattern was processed to identify the variants identified in immune-related genes that show a moderate/high impact on the encoded protein function. Among a list of 85 missense variants concordant with the QTL segregation pattern that were within candidate immune-related genes, 13 variants distributed across 7 genes [PKHD1, NOTCH4, AGER, ENSOARG00000009395 (HLA-C, Homo sapiens), ENSOARG00000015002 (HLA-B, H. sapiens), MOG, and ENSOARG00000018075 (BoLA, Bos taurus, orthologous to human HLA-A] were predicted to cause deleterious effects on protein function. Future studies should assess the possible associations of the candidate variants identified herein in commercial populations with indicator traits of udder inflammation (SCS, clinical mastitis).

Transcriptome analysis and identification of significantly differentially expressed genes in Holstein calves subjected to severe thermal stress

RNA-Seq analysis was used to characterize transcriptome response of Holstein calves to thermal stress. A total of eight animals aged between 2 and 3 months were randomly selected and subjected to thermal stress corresponding to a temperature humidity index of 95 in an environmentally controlled house for 12 h consecutively for 3 days. A set of 15,787 unigenes were found to be expressed and after a threshold of threefold change, and a Q value <0.05; 502, 394, and 376 genes were found to be differentially expressed on days 1, 2, and 3 out of which 343, 261 and 256 genes were upregulated and 159, 133, and 120 genes were downregulated. Only 356 genes out of these were expressed on all 3 days, and only they were considered as significantly differentially expressed. KEGG pathway analysis revealed that ten pathways were significantly enriched; the top two among them were protein processing in endoplasmic reticulum and MAPK signaling pathways. These results suggest that thermal stress triggered a complex response in Holstein calves and the animals adjusted their physiological and metabolic processes to survive. Many of the genes identified in this study have not been previously reported to be involved in thermal stress response. The results of this study extend our understanding of the animal's response to thermal stress and some of the identified genes may prove useful in the efforts to breed Holstein cattle with superior thermotolerance, which might help in minimizing production loss due to thermal stress.

Effect of swine leukocyte antigen-DQA gene variation on diarrhea in Large White, Landrace, and Duroc piglets

Piglet diarrhea is one of the most common factors that affects the benefits of the swine industry. Although recent studies have shown that exon 2 of SLA-DQA is associated with piglet resistance to diarrhea, contributions of genetic variation in the additional exon coding regions of this gene remain unclear. Here, we investigated variation in exons 1, 3 and 4 of the SLA-DQA gene and evaluated their effects on diarrheal infection in 425 suckling piglets. No variation was identified in exon 1. In exon 3, there were eight alleles detected, generated by 14 single nucleotide polymorphisms (SNPs) and three nucleotide deletions, eight SNPs being newly identified. Four allele sequences and three SNPs were identified in exon 4, only one SNP being newly identified. Statistical analysis showed that the genotypes of exon 3 are significantly associated with piglet diarrhea; indeed, genotypes DQA*wb01/wb02 and wb04/wb05 are clearly associated with resistance to piglet diarrhea, as they have the lowest probabilities of infection (P < 0.05). However, no significant association was found between the genotypes of exon 4 and diarrhea (P > 0.05). These results provide important new information concerning the level of genetic diversity at the SLA-DQA locus and suggest that further genetic association studies of piglet diarrhea resistance should include analyses of both exons 2 and 3 of this locus.

Characterization of bovine MHC class II DRB3 diversity in South American Holstein cattle populations

Holstein cattle dominate the global milk production industry because of their outstanding milk production, however, this breed is susceptible to tropical endemic pathogens and suffers from heat stress and thus fewer Holstein populations are raised in tropical areas. The bovine major histocompatibility complex (BoLA)-DRB3 class II gene is used as a marker for disease and immunological traits, and its polymorphism has been studied extensively in Holstein cattle from temperate and cold regions. We studied the genetic diversity of the BoLA-DRB3 gene in South American Holstein populations to determine whether tropical populations have diverged from those bred in temperate and cold regions by selection and/or crossbreeding with local native breeds. We specifically studied Exon 2 of this gene from 855 South American Holstein individuals by a polymerase chain reaction (PCR) sequence-based typing method. We found a high degree of gene diversity at the allelic (Na > 20 and He > 0.87) and molecular (π > 0.080) levels, but a low degree of population structure (FST = 0.009215). A principal components analysis and tree showed that the Bolivian subtropical population had the largest genetic divergence compared with Holsteins bred in temperate or cold regions, and that this population was closely related to Bolivian Creole cattle. Our results suggest that Holstein genetic divergence can be explained by selection and/or gene introgression from local germplasms. This is the first examination of BoLA-DRB3 in Holsteins adapted to tropical environments, and contributes to an ongoing effort to catalog bovine MHC allele frequencies by breed and location.

Short Communication: New alleles at the BoLA-DQA1 locus in Holstein–Fresian cattle

Kułajaj, D., Pokorska, J., Ormian, M. and Dusza, M. 2015. Short Communication: New alleles at the BoLA-DQA1 locus in Holstein Fresian cattle. Can. J. Anim. Sci. 95: 161–164. The aim of the study was to analyze genetic polymorphism at locus BoLA-DQA1 in Polish Holstein-Friesian cattle, using next generation sequencing technology (Miseq sequencer – Illumina). In the studied cattle population we identified 14 alleles, of which three were determined for the first time (BoLA-DQA1*0106, BoLA-DQA1*3201, BoLA-DQA1*3301). The frequencies of these alleles were respectively: 0.0423, 0.0282, 0.0282. At locus BoLA-DQA1 the highest frequency was found for the BoLA-DQA1*10011 and BoLA-DQA1*0101 alleles and their frequencies were respectively: 0.3592 and 0.2606. The frequencies of other alleles identified in the studied population ranged from 0.0070 to 0.0704. The results obtained in this study indicate that the BoLA-DQA1 locus is highly polymorphic and new alleles are still identifiable.

Discovering a vaccine against neosporosis using computers: is it feasible?

A vaccine is urgently needed to prevent cattle neosporosis. This infectious disease is caused by the parasite Neospora caninum, a complex biological system with multifaceted life cycles. An in silico vaccine discovery approach attempts to transform digital abstractions of this system into adequate knowledge to predict candidates. Researchers need current information to implement such an approach, such as understanding evasion mechanisms of the immune system, type of immune response to elicit, availability of data and prediction programs, and statistical models to analyze predictions. Taken together, an in silico approach involves assembly of an intricate jigsaw of interdisciplinary and interdependent knowledge. In this review, we focus on the approach influencing vaccine development against Neospora caninum, which can be generalized to other pathogenic apicomplexans.

The great diversity of major histocompatibility complex class II genes in Philippine native cattle

Bovine leukocyte antigens (BoLA) are extensively used as markers for bovine disease and immunological traits. However, none of the BoLA genes in Southeast Asian breeds have been characterized by polymerase chain reaction (PCR)-sequence-based typing (SBT). Therefore, we sequenced exon 2 of the BoLA class II DRB3 gene from 1120 individual cows belonging to the Holstein, Sahiwal, Simbrah, Jersey, Brahman, and Philippine native breeds using PCR-SBT. Several cross-breeds were also examined. BoLA-DRB3 PCR-SBT identified 78 previously reported alleles and five novel alleles. The number of BoLA-DRB3 alleles identified in each breed from the Philippines was higher (71 in Philippine native cattle, 58 in Brahman, 46 in Holstein × Sahiwal, and 57 in Philippine native × Brahman) than that identified in breeds from other countries (e.g., 23 alleles in Japanese Black and 35 in Bolivian Yacumeño cattle). A phylogenetic tree based on the D_A distance calculated from the BoLA-DRB3 allele frequency showed that Philippine native cattle from different Philippine islands are closely related, and all of them are closely similar to Philippine Brahman cattle but not to native Japanese and Latin American breeds. Furthermore, the BoLA-DRB3 allele frequency in Philippine native cattle from Luzon Island, located in the Northern Philippines was different from that in cattle from Iloilo, Bohol, and Leyte Islands, which are located in the Southern Philippines. Therefore, we conclude that Philippine native cattle can be divided into two populations, North and South areas. Moreover, a neutrality test revealed that Philippine native cattle from Leyte showed significantly greater genetic diversity, which may be maintained by balancing selection. This study shows that Asian breeds have high levels of BoLA-DRB3 polymorphism. This finding, especially the identification of five novel BoLA-DRB3 alleles, will be helpful for future SBT studies of BoLA-DRB3 alleles in East Asian cattle.

Haplotype determination of the upstream regulatory region and the second exon of the BoLA-DRB3 gene in Holstein cattle

Polymorphisms of the BoLA-DRB3 gene are located primarily in the second exon [antigen binding site (ABS)] and, to a lesser extent, in the upstream regulatory region (URR). It can be hypothesised that exon 2 and the URR are under different types of natural selection. The aim of this work was to determine the URR-exon 2 haplotypes; 34 Holstein samples were genotyped by direct sequencing. A total of 7 URR alleles and 23 exon 2 alleles were detected, and 3 of the URR alleles were novel. Our results may suggest that no relationship exists between the URR and exon 2 of the BoLA-DRB3 gene (linkage disequilibrium P value > 0.05), most likely due to recombination over time. Our results also suggest that both regions of class II genes may be included in the development of new genotyping methods based on next-generation DNA sequencing technologies.

A study of the association between chronic superficial keratitis and polymorphisms in the upstream regulatory regions of DLA-DRB1, DLA-DQB1 and DLA-DQA1

Canine chronic superficial keratitis (CSK) is an inflammatory corneal disease that primarily occurs in German shepherd dogs (GSDs). Several studies support the hypothesis that CSK is an immune-mediated disease. To investigate the genetic factors associated with CSK development, the upstream regulatory regions (URRs) of the DLA-DRB, -DQA and -DQB genes were genotyped in 60 dogs, including 32 CSK animals. LD analysis identified two blocks (r(2)≤45), with two DLA-DRB1 and five DLA-DQB1 haplotypes. Analysis of DLA-URR alleles/haplotypes showed a significant association between DQB1*-154 [C/T] (p=0.016) and CSK, suggesting that the T variant may increase the risk for developing CSK disease (OR=3, 95% CI=1.25-7.68). When haplotype associations were performed, the URR-DQB*CATT haplotype was significantly associated with CSK (p=0.016), increasing the risk of develop this disease over two-fold (OR=3, 95%, CI=1.25-7.68). These results showed that dogs homozygous at DRB1*69 [C/T] had a risk for developing CSK disease that was over four times the risk for heterozygotes. This genetic association supports the previous clinical, histological and pharmacological studies that suggest that CSK is an immune-mediated disease, and this association could potentially be used to identify susceptible animals.

BLV-CoCoMo-qPCR: a useful tool for evaluating bovine leukemia virus infection status

BACKGROUND

Bovine leukemia virus (BLV) is associated with enzootic bovine leukosis, which is the most common neoplastic disease of cattle. BLV infects cattle worldwide, imposing a severe economic impact on the dairy cattle industry. Recently, we developed a new quantitative real-time polymerase chain reaction (PCR) method using Coordination of Common Motifs (CoCoMo) primers to measure the proviral load of known and novel BLV variants in BLV-infected animals. Indeed, the assay was highly effective in detecting BLV in cattle from a range of international locations. This assay enabled us to demonstrate that proviral load correlates not only with BLV infection capacity as assessed by syncytium formation, but also with BLV disease progression. In this study, we compared the sensitivity of our BLV-CoCoMo-qPCR method for detecting BLV proviruses with the sensitivities of two real-time PCR systems, and also determined the differences of proviral load with serotests.

RESULTS

BLV-CoCoMo-qPCR was found to be highly sensitive when compared with the real-time PCR-based TaqMan MGB assay developed by Lew et al. and the commercial TaKaRa cycleave PCR system. The BLV copy number determined by BLV-CoCoMo-qPCR was only partially correlated with the positive rate for anti-BLV antibody as determined by the enzyme-linked immunosorbent assay, passive hemagglutination reaction, or agar gel immunodiffusion. This result indicates that, although serotests are widely used for the diagnosis of BLV infection, it is difficult to detect BLV infection with confidence by using serological tests alone. Two cattle were experimentally infected with BLV. The kinetics of the provirus did not precisely correlate with the change in anti-BLV antibody production. Moreover, both reactions were different in cattle that carried different bovine leukocyte antigen (BoLA)-DRB3 genotypes.

CONCLUSIONS

Our results suggest that the quantitative measurement of proviral load by BLV-CoCoMo-qPCR is useful tool for evaluating the progression of BLV-induced disease. BLV-CoCoMo-qPCR allows us to monitor the spread of BLV infection in different viewpoint compared with classical serotest.

Major histocompatibility complex class II compatibility, but not class I, predicts mate choice in a bird with highly developed olfaction

Mate choice for major histocompatibility complex (MHC) compatibility has been found in several taxa, although rarely in birds. MHC is a crucial component in adaptive immunity and by choosing an MHC-dissimilar partner, heterozygosity and potentially broad pathogen resistance is maximized in the offspring. The MHC genotype influences odour cues and preferences in mammals and fish and hence olfactory-based mate choice can occur. We tested whether blue petrels, Halobaena caerulea, choose partners based on MHC compatibility. This bird is long-lived, monogamous and can discriminate between individual odours using olfaction, which makes it exceptionally well suited for this analysis. We screened MHC class I and II B alleles in blue petrels using 454-pyrosequencing and quantified the phylogenetic, functional and allele-sharing similarity between individuals. Partners were functionally more dissimilar at the MHC class II B loci than expected from random mating (p = 0.033), whereas there was no such difference at the MHC class I loci. Phylogenetic and non-sequence-based MHC allele-sharing measures detected no MHC dissimilarity between partners for either MHC class I or II B. Our study provides evidence of mate choice for MHC compatibility in a bird with a high dependency on odour cues, suggesting that MHC odour-mediated mate choice occurs in birds.

Identification and diversity of bovine major histocompatibility complex class II haplotypes in Japanese Black and Holstein cattle in Japan

Bovine leukocyte antigen (BoLA), the major histocompatibility complex of cattle, is one of the most polymorphic gene clusters. We genotyped a population of 109 Japanese Black and 39 Holstein cattle to analyze their BoLA class II haplotypes, BoLA-DRB3 locus, 5 BoLA-DQA loci, and 5 BoLA-DQB loci. We identified 26 previously reported DRB3 alleles, 22 previously reported and 3 new DQA alleles, and 24 previously reported and 6 new DQB alleles. A dendrogram was constructed based on the predicted amino acid sequences of the α1 or β1 domains encoded by BoLA-DQA or -DQB alleles, which revealed that DQA alleles were clustered into 5 loci, whereas DQB alleles could not be clearly assigned to specific DQB loci. The BoLA-DRB3-DQA-DQB haplotypes were sorted by sequential analytical processes, and 42 distinct haplotypes, including 11 previously published haplotypes and 31 novel haplotypes, were defined. Strong linkage disequilibrium was present in the BoLA genes. We also compared DRB3-DQA1 haplotype frequencies between 507 Japanese Black and 143 Holstein cattle. Thirty-nine DRB3-DQA1 haplotypes were identified, including 29 haplotypes from Japanese Black and 22 haplotypes from Holstein cattle. The majority of the haplotypes could be identified in both breeds, although several haplotypes were identified in only a single breed. This is the first report presenting a detailed study of the BoLA class II haplotype in Japanese Black and Holstein cattle in Japan.

Expression of genes associated with the antigen presentation and processing pathway are consistently regulated in early Mycobacterium avium subsp. paratuberculosis infection

The primary objective of this study was to evaluate early gene expression profiles associated with paratuberculosis in cattle exposed to known infectious doses of Mycobacterium avium sub-species paratuberculosis (MAP). A Johne's disease experimental infection field trial was conducted on a mixed population of Holstein and Holstein Red cattle. Blood samples from four MAP exposed and four unexposed cattle, selected based on IFNγ expressions were taken at 9, 13 and 21 weeks and RNA processed to Affymetrix GeneChip™ Bovine Genome arrays. Ontological analysis revealed consistent differences in gene expression between MAP exposed and control animals. A stark variation was observed in expression of a number of genes along antigen presentation pathways, suggesting that MAP exposure potentially results in the host immune response switching to a CD8(+) biased antigen presentation profile. This requires further in-depth analysis since it exposes a hitherto unconfirmed association between MAP exposure and in vivo MHC gene modulation.

Functional adaptations of the transcriptome to mastitis-causing pathogens: the mammary gland and beyond

Application of microarrays to the study of intramammary infections in recent years has provided a wealth of fundamental information on the transcriptomics adaptation of tissue/cells to the disease. Due to its heavy toll on productivity and health of the animal, in vivo and in vitro transcriptomics works involving different mastitis-causing pathogens have been conducted on the mammary gland, primarily on livestock species such as cow and sheep, with few studies in non-ruminants. However, the response to an infectious challenge originating in the mammary gland elicits systemic responses in the animal and encompasses tissues such as liver and immune cells in the circulation, with also potential effects on other tissues such as adipose. The susceptibility of the animal to develop mastitis likely is affected by factors beyond the mammary gland, e.g. negative energy balance as it occurs around parturition. Objectives of this review are to discuss the use of systems biology concepts for the holistic study of animal responses to intramammary infection; providing an update of recent work using transcriptomics to study mammary and peripheral tissue (i.e. liver) as well as neutrophils and macrophage responses to mastitis-causing pathogens; discuss the effect of negative energy balance on mastitis predisposition; and analyze the bovine and murine mammary innate-immune responses during lactation and involution using a novel functional analysis approach to uncover potential predisposing factors to mastitis throughout an animal's productive life.

Eca20 microsatellite polymorphisms in equine viral arteritis-infected horses from Argentina

We investigated the association of equine arteritis virus (EAV) infection and three short tandem repeat (STR) polymorphisms located within or in close proximity to equine lymphocyte antigen (ELA) region. We used a case-control design as a first approach before proceeding to select candidate genes. One hundred and sixty-five Silla Argentino horses were taken in 2002 from positive serological detections of EAV in Argentina, to determine whether STR genotypes were correlated to genetic susceptibility to EVA. Allele frequency distribution did not show significant differences between both groups (P = 0.0781). However, in particular alleles, Fisher exact test and odds ratio calculations showed significant values >1 for TKY08 and LEX52, and <1 for UM011, TKY08, LEX52 and VHL20. Interestingly, TKY08 STR is located in ELA class I region.

A new method for typing bovine major histocompatibility complex class II DRB3 alleles by combining two established PCR sequence-based techniques

Recently, two polymerase chain reaction sequence-based typing (PCR-SBT) methods were reported for the genotyping of the bovine leukocyte antigen (BoLA)-DRB3. One technique is a single PCR-SBT (sPCR-SBT) method that generates heterozygous sequences that are subsequently analyzed by the haplofinder program, while the other technique is a nested PCR-SBT (nPCR-SBT) method that allows the analysis of heterozygous sequences using the assign 400ATF software. In this study, these techniques were compared and then integrated to produce an improved genotyping method. The primer set used for sPCR-SBT was more accurate than those used for nPCR-SBT. Combining sPCR-SBT with the assign 400ATF software previously reported for nPCR-SBT enables rapid and accurate genotyping of a large number of DNA samples.

Escherichia coli infection induces distinct local and systemic transcriptome responses in the mammary gland

Background

Coliform bacteria are the most common etiologic agents in severe mastitis of cows. Escherichia coli infections are mostly restricted to a single udder quarter whereas neighboring quarters stay clinically inapparent, implicating the presence of a systemic defense reaction. To address its underlying mechanism, we performed a transcriptome study of mammary tissue from udder quarters inoculated with E. coli (6 h and 24 h post infection), from neighboring quarters of the same animals, and from untreated control animals.

Results

After 6 h 13 probe sets of differentially expressed genes (DEG) were detected in infected quarters versus control animals. Eighteen hours later 2154 and 476 DEG were found in infected and in neighboring quarters vs. control animals. Cluster analysis revealed DEG found only in infected quarters (local response) and DEG detected in both infected and neighboring quarters (systemic response). The first group includes genes mainly involved in immune response and inflammation, while the systemic reaction comprises antigen processing and presentation, cytokines, protein degradation and apoptosis. Enhanced expression of antimicrobial genes (S100A8, S100A9, S100A12, CXCL2, GNLY), acute phase genes (LBP, SAA3, CP, BF, C6, C4BPA, IF), and indicators of oxidative stress (GPX3, MT1A, MT2A, SOD2) point to an active defense reaction in infected and neighboring healthy quarters. Its early onset is indicated by increased transcription of NFIL3 at 6 h. NFIL3 is a predicted regulator of many genes of the systemic response at 24 h. The significance of our transcriptome study was evidenced by some recent findings with candidate gene based approaches.

Conclusions

The discovery and holistic analysis of an extensive systemic reaction in the mammary gland significantly expands the knowledge of host-pathogen interactions in mastitis which may be relevant for the development of novel therapies and for genetic selection towards mastitis resistance.