Association between BoLA-DRB3 and somatic cell count in Holstein cattle from Argentina

More than Three Decades of Bm86: What We Know and Where to Go

Tick and tick-borne disease control have been a serious research focus for many decades. In a global climate of increasing acaricide resistance, host immunity against tick infestation has become a much-needed complementary strategy to common chemical control. From the earliest acquired resistance studies in small animal models to proof of concept in large production animals, it was the isolation, characterization, and final recombinant protein production of the midgut antigen Bm86 from the Australian cattle tick strain of Rhipicephalus (Boophilus) microplus (later reinstated as R. (B.) australis) that established tick subunit vaccines as a viable alternative in tick and tick-borne disease control. In the past 37 years, this antigen has spawned numerous tick subunit vaccines (either Bm86-based or novel), and though we are still describing its molecular structure and function, this antigen remains the gold standard for all tick vaccines. In this paper, advances in tick vaccine development over the past three decades are discussed alongside the development of biotechnology, where existing gaps and future directives in the field are highlighted.

Candidate genes for mastitis resistance in dairy cattle: a data integration approach

BACKGROUND

Inflammation of the mammary tissue (mastitis) is one of the most detrimental health conditions in dairy ruminants and is considered the most economically important infectious disease of the dairy sector. Improving mastitis resistance is becoming an important goal in dairy ruminant breeding programmes. However, mastitis resistance is a complex trait and identification of mastitis-associated alleles in livestock is difficult. Currently, the only applicable approach to identify candidate loci for complex traits in large farm animals is to combine different information that supports the functionality of the identified genomic regions with respect to a complex trait.

METHODS

To identify the most promising candidate loci for mastitis resistance we integrated heterogeneous data from multiple sources and compiled the information into a comprehensive database of mastitis-associated candidate loci. Mastitis-associated candidate genes reported in association, expression, and mouse model studies were collected by searching the relevant literature and databases. The collected data were integrated into a single database, screened for overlaps, and used for gene set enrichment analysis.

RESULTS

The database contains candidate genes from association and expression studies and relevant transgenic mouse models. The 2448 collected candidate loci are evenly distributed across bovine chromosomes. Data integration and analysis revealed overlaps between different studies and/or with mastitis-associated QTL, revealing promising candidate genes for mastitis resistance.

CONCLUSION

Mastitis resistance is a complex trait influenced by numerous alleles. Based on the number of independent studies, we were able to prioritise candidate genes and propose a list of the 22 most promising. To our knowledge this is the most comprehensive database of mastitis associated candidate genes and could be helpful in selecting genes for functional validation studies.

High throughput analysis of MHC-I and MHC-DR diversity of Brazilian cattle populations

The major histocompatibility complex (MHC) contains many genes that play key roles in initiating and regulating immune responses. This includes the polymorphic MHCI and MHCII genes that present epitopes to CD8+ and CD4+ T-cells, respectively. Consequently, the characterisation of the repertoire of MHC genes is an important component of improving our understanding of the genetic variation that determines the outcomes of immune responses. In cattle, MHC (BoLA) research has predominantly focused on Holstein-Friesian animals (as the most economically important breed globally), although the development of high-throughput approaches has allowed the BoLA-DRB3 repertoire to be studied in a greater variety of breeds. In a previous study we reported on the development of a MiSeq-based method to enable high-throughput and high-resolution analysis of bovine MHCI repertoires. Herein, we report on the expansion of this methodology to incorporate analysis of the BoLA-DRB3 and its application to analyse MHC diversity in a large cohort of cattle from Brazil (>500 animals), including representatives from the three major Bos indicus breeds present in Brazil - Guzerat, Gir and Nelore. This large-scale description of paired MHCI-DRB3 repertoires in Bos indicus cattle has identified a small number of novel DRB3 alleles, a large number of novel MHCI alleles and haplotypes, and provided novel insights into MHCI-MHCII association - further expanding our knowledge of bovine MHC diversity.

Copy number variations in BOLA-DQA2, BOLA-DQB, and BOLA-DQA5 show the genomic architecture and haplotype frequency of major histocompatibility complex class II genes in Holstein cows

Bovine major histocompatibility complex (MHC) class II region contains many genes. The bovine leukocyte antigen (BoLA)-DRB3 was reportedly associated with susceptibility of various phenotypes of infections including bovine leukemia virus-induced lymphoma. However, the association of the remaining genes with various phenotypes has not been clarified due to the complicated genomic structure of the MHC class II region. Thus, in this study, we elucidated the MHC class II genomic structure, including the novel alleles and copy number variations (CNVs). We determined the copy numbers of BOLA-DQA2 (DQA2), BOLA-DQB (DQB2), BOLA-DQA5 (DQA5), BLA-DQB (DQB1), LOC100848815 (DQA1), and BOLA-DRB3 (DRB3) in 127 unrelated Holstein cows by TaqMan copy number assay. The genomes were sequenced using target next-generation sequencing (NGS) based on multiplex polymerase chain reaction. Combining the results of the copy numbers and alleles, we identified the BoLA alleles directly without haplotype estimation. Pairwise linkage disequilibrium (LD) analysis between alleles and genes were performed. The CNVs of DQA2, DQB2, and DQA5 in Holstein cows were detected. The frequency of the whole gene deletion in DQA2, DQB2, and DQA5 was 35.4%, 93.7%, and 93.7%, respectively. After target NGS, we identified 37 alleles in the six genes. Fifteen novel alleles (40.5%) were not registered in the IPD-MHC Database. LD analysis showed strong LD among the DQB2*deletion, DQA5*deletion, and DRB3*27:03 alleles. Our findings will provide important insights into the identification of the BoLA genes associated with various infection-related phenotypes.

Genetic mechanisms regulating the host response during mastitis

Mastitis is a very costly and common disease in the dairy industry. The study of the transcriptome from healthy and mastitic milk somatic cell samples using RNA-Sequencing technology can provide measurements of transcript levels associated with the immune response to the infection. The objective of this study was to characterize the Holstein milk somatic cell transcriptome from 6 cows to determine host response to intramammary infections. RNA-Sequencing was performed on 2 samples from each cow from 2 separate quarters, one classified as healthy (n = 6) and one as mastitic (n = 6). In total, 449 genes were differentially expressed between the healthy and mastitic quarters (false discovery rate <0.05, fold change >±2). Among the differentially expressed genes, the most expressed genes based on reads per kilobase per million mapped reads (RPKM) in the healthy group were associated with milk components (CSN2 and CSN3), and in the mastitic group they were associated with immunity (B2M and CD74). In silico functional analysis was performed using the list of 449 differentially expressed genes, which identified 36 significantly enriched metabolic pathways (false discovery rate <0.01), some of which were associated with the immune system, such as cytokine-cytokine interaction and cell adhesion molecules. Seven functional candidate genes were selected, based on the criteria of being highly differentially expressed between healthy and mastitic groups and significantly enriched in metabolic pathways that are relevant to the inflammatory process (GLYCAM1, B2M, CD74, BoLA-DRA, FCER1G, SDS, and NFKBIA). Last, we identified the differentially expressed genes that are located in quantitative trait locus regions previously known to be associated with mastitis, specifically clinical mastitis, somatic cell count, and somatic cell score. It was concluded that multiple genes within quantitative trait locus regions could potentially affect host response to mastitis-causing agents, making some cows more susceptible to intramammary infections. The identification of potential candidate genes with functional, statistical, biological, and positional relevance associated with host defense to infection will contribute to a better understanding of the underlying genetic architecture associated with mastitis. This in turn will improve the sustainability of agricultural practices by facilitating the selection of cows with improved host defense leading to increased resistance to mastitis.

Molecular characterization of MHC II DRB3 gene of swamp- and riverine-type water buffaloes

Objective

Major histocompatibility complex (MHC) is a set of molecular proteins on the surface of antigen presenting cells encoded by a large gene family which are important parts of the immune system. This study was conducted to convey information on the genetic characteristics of the MHC II DRB3 gene in riverine and swamp buffaloes.

Materials and Methods

Characterization of MHC II DRB3 gene was carried out using polymerase chain reaction (PCR)-based assay. Thirty-milliliter milk samples were collected from 10 swamp-type and 10 riverine-type buffaloes. RNA from milk samples were extracted using Trizol and then followed by reverse transcription-PCR (RT-PCR).

Results

The phylogenetic analysis with 1,000 bootstrap replications clearly showed complex parsimony in MHC II DRB3 gene between 10 riverine- and 10 swamp-type but also confirmed that the samples are similar to Bubalus bubalis. Aligned sequences of the 20 water buffaloes were compared with three other ruminants (Bos taurus, Ovis aries, and Capra hircus) and non-ruminant (Sus scrofa) that serve as an outgroup. MHC sequences from GenBank show that there was an average of 705 identical pairs, with 22 transitional pairs and 30 transversional pairs with a ratio of 0.7.

Conclusion

Based on the molecular data, the current study conforms to other works of literature that this gene is highly polymorphic which can be due to its function in the immune responsiveness and disease resistance. Further study on the immunological response of MHC II DRB3 to infection may elucidate its underlying function and role in the protection against specific disease of animals.

Sequence diversity of major histo-compatibility complex class II DQA1 in Indian Tharparkar cattle: novel alleles and in-silico analysis

Exon 2 of MHC class II gene codes for the first domain of the molecule that forms the peptide-binding groove and its polymorphism partly explains functional MHC diversity. A 850 bp DQA1 gene fragment spanning from intron I to exon III was typed by sequencing of 40 Tharparkar cattle of various agro-climatic zones of northern India along with 10 Tharparkar crossbreds. On analysis of nucleotide sequences, a total of 30 polymorphic sites (1 insertion and 29 SNPs) were identified in 14 MHC alleles leading to amino acid changes in 5 places in 249 bp (exon 2). Five new BoLa DQA1 alleles were identified and reported. The within group mean distance was highest in Tharparkar herd of Bikaner (0.045) and lowest (0.020) in that of Surathgarh (breeding tract) whereas, between groups mean distance was highest in Bikaner Tharparkar-Suratgarh Tharparkar pair. There was excess of nonsynonymous over synonymous nucleotide substitutions in the present study. The effects of these substitutions were predicted using I-Mutant and Panther online resources. The mean ratio of dN/dS was found to be >1.0 at 12 codons with two mutation hotspots at 13th codon (P = 0.002) and 64th codon (P = 0.01). The phylo-geographic analysis revealed that alleles 5, 7 and 13 formed a different cluster with alleles 7 and 13 grouped by the most frequent allele (BoLa-DQA*1401).

Prevalence of bovine dermatophilosis and disease-associated alleles in zebu Goudali cattle and their Italian Simmental crosses ranching in the western highland plateau savannah of Cameroon

Abundance of native pastures makes Cameroon's western highland savannah (WHS) a hotspot for low-input beef-type cattle. Dumbo Ranch is central to cattle seed stock multiplication in WHS and holds that Dermatophilus congolensis infection undermines production. The bovine BoLA-DRB3 has been variously demonstrated as the principal gene of the major histocompatibility locus associated with immunity and resistance to dermatophilosis in cattle. We studied the profile of dermatophilosis prevalence in zebu Goudali (G) and its Simmental composite, SimGoud (SG), at Dumbo Ranch and determined the distribution of a dermatophilosis-associated susceptibility allele of the BoLA-DRB3 gene by allele-specific polymerase chain reaction (PCR). We recorded a 42 % prevalence of dermatophilosis in the studied cohort (337 animals). Dermatophilosis was more common in older cattle than in cattle ≤36 months (p ≤ 0.05). G was more affected compared to SG, because of the prevalence of the disease in the oldest animals and the age distribution of the experimental subjects. No susceptible homozygote was observed. About 85 and 15 % of the cohort carried the homozygous resistant and heterozygous condition, respectively. This genotype distribution was not affected by cattle type. The study confirms the presence of dermatophilosis among G and SG cattle in WHS. However, there was no correlation between the presence of the disease-associated susceptible allele considered and clinical manifestation. Screening for this dermatophilosis resistance-associated allele of BoLA-DRB3 gene appeared not useful for selection of G and SG in WHS.

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.

Locus BoLA-DRB3 is just an ordinary site of the polygene when explaining genetic variance of somatic cell count and milk yield

The study aimed at clarifying the problem of the hitherto contradictory results regarding usefulness of BoLA-DRB3 locus as a marker in selection against mastitis and for milk yield. Treating the BoLA-DRB3 locus effect as random was proposed in place of considering it fixed. Somatic cell counts and milk yields recorded monthly on a test day (22 424) of 619 Polish Holstein cows genotyped for BoLA-DRB3 were analysed with an animal model including a random effect for genotype at this locus. The BoLA-DRB3 alleles were defined as restriction patterns obtained with three endonucleases. Two alternative BoLA-DRB3 additive genotype (co)variance structures were constructed for 161 genotypes recorded. One was based on the allelic similarity of the genotypes resulting in element values of 0 (no common allele), 0·5 (one allele in common), and 1 (diagonal). The other considered restriction site similarity (up to 3 in 1 allele) giving element values of 0 (no common restriction sites) and then increasingly in steps of 1/6 up to 6/6 (diagonal), where the numerator represents the number of common sites between genotypes. The DRB3 variance component for the natural logarithm of somatic cell count did not exceed 0·006 of the polygenic additive component or 0·003 for milk yield. Hence, unless we fail to detect the causative site or to properly define traits being the projection of a site, the effect of the genotype at the BoLA-DRB3 locus does not explain variation in somatic cell count and milk yield at a degree expected of a genetic marker.