Variance components for bovine tuberculosis infection and multi-breed genome-wide association analysis using imputed whole genome sequence data

Integrative and comparative genomic analyses of mammalian macrophage responses to intracellular mycobacterial pathogens

Mycobacterium tuberculosis, the causative agent of human tuberculosis (hTB), is a close evolutionary relative of Mycobacterium bovis, which causes bovine tuberculosis (bTB), one of the most damaging infectious diseases to livestock agriculture. Previous studies have shown that the pathogenesis of bTB disease is comparable to hTB disease, and that the bovine and human alveolar macrophage (bAM and hAM, respectively) transcriptomes are extensively reprogrammed in response to infection with these intracellular mycobacterial pathogens. In this study, a multi-omics integrative approach was applied with functional genomics and GWAS data sets across the two primary hosts (Bos taurus and Homo sapiens) and both pathogens (M. bovis and M. tuberculosis). Four different experimental infection groups were used: 1) bAM infected with M. bovis, 2) bAM infected with M. tuberculosis, 3) hAM infected with M. tuberculosis, and 4) human monocyte-derived macrophages (hMDM) infected with M. tuberculosis. RNA-seq data from these experiments 24 h post-infection (24 hpi) was analysed using three computational pipelines: 1) differentially expressed genes, 2) differential gene expression interaction networks, and 3) combined pathway analysis. The results were integrated with high-resolution bovine and human GWAS data sets to detect novel quantitative trait loci (QTLs) for resistance to mycobacterial infection and resilience to disease. This revealed common and unique response macrophage pathways for both pathogens and identified 32 genes (12 bovine and 20 human) significantly enriched for SNPs associated with disease resistance, the majority of which encode key components of the NF-κB signalling pathway and that also drive formation of the granuloma.

Genome-Wide Association Study Reveals Quantitative Trait Loci and Candidate Genes Associated with High Interferon-gamma Production in Holstein Cattle Naturally Infected with Mycobacterium Bovis

Mycobacterium bovis (Mb) is the causative agent of bovine tuberculosis (bTb). Genetic selection aiming to identify less susceptible animals has been proposed as a complementary measure in ongoing programs toward controlling Mb infection. However, individual animal phenotypes for bTb based on interferon-gamma (IFNɣ) and its use in bovine selective breeding programs have not been explored. In the current study, IFNɣ production was measured using a specific IFNɣ ELISA kit in bovine purified protein derivative (bPPD)-stimulated blood samples collected from Holstein cattle. DNA isolated from the peripheral blood samples collected from the animals included in the study was genotyped with the EuroG Medium Density bead Chip, and the genotypes were imputed to whole-genome sequences. A genome-wide association analysis (GWAS) revealed that the IFNɣ in response to bPPD was associated with a specific genetic profile (heritability = 0.23) and allowed the identification of 163 SNPs, 72 quantitative trait loci (QTLs), 197 candidate genes, and 8 microRNAs (miRNAs) associated with this phenotype. No negative correlations between this phenotype and other phenotypes and traits included in the Spanish breeding program were observed. Taken together, our results define a heritable and distinct immunogenetic profile associated with strong production of IFNɣ in response to Mb.

Will we ever eradicate animal tuberculosis?

Two characteristics of the Mycobacterium tuberculosis complex (MTC) are particularly relevant for tuberculosis (TB) epidemiology and control, namely the ability of this group of pathogens to survive in the environment and thereby facilitate indirect transmission via water or feed, and the capacity to infect multiple host species including human beings, cattle, wildlife, and domestic animals other than cattle. As a consequence, rather than keeping the focus on certain animal species regarded as maintenance hosts, we postulate that it is time to think of complex and dynamic multi-host MTC maintenance communities where several wild and domestic species and the environment contribute to pathogen maintenance. Regarding the global situation of animal TB, many industrialized countries have reached the Officially Tuberculosis Free status. However, infection of cattle with M. bovis still occurs in most countries around the world. In low- and middle-income countries, human and animal TB infection is endemic and bovine TB control programs are often not implemented because standard TB control through testing and culling, movement control and slaughterhouse inspection is too expensive or ethically unacceptable. In facing increasingly complex epidemiological scenarios, modern integrated disease control should rely on three main pillars: (1) a close involvement of farmers including collaborative decision making, (2) expanding the surveillance and control targets to all three host categories, the environment, and their interactions, and (3) setting up new control schemes or upgrading established ones switching from single tool test and cull approaches to integrated ones including farm biosafety and vaccination.

Selective breeding can contribute to bovine tuberculosis control and eradication

Bovine tuberculosis (bTB) persists in many countries having a significant impact on public health and livestock industry finances. The incidence and prevalence of new cases in parts of the UK and elsewhere over the past decades warrant intensified efforts towards achieving Officially Tuberculosis Free (OTF) status in the respective regions. Genetic selection aiming to identify and remove inherently susceptible animals from breeding has been proposed as an additional measure in ongoing programmes towards controlling the disease. The presence of genetic variation among individual animals in their capacity to respond to Mycobacterium bovis exposure has been documented and heritability estimates of 0.06-0.18 have been reported. Despite their moderate magnitude, these estimates suggest that host resistance to bTB is amenable to improvement with selective breeding. Although relatively slow, genetic progress can be constant, cumulative and permanent, thereby complementing ongoing disease control measures. Importantly, mostly no antagonistic genetic correlations have been found between bTB resistance and other animal traits suggesting that carefully incorporating the former in breeding decisions should not adversely affect bovine productivity. Simulation studies have demonstrated the potential impact of genetic selection on reducing the probability of a breakdown to occur or the duration and severity of a breakdown that has already been declared. Furthermore, research on the bovine genome has identified multiple genomic markers and genes associated with bTB resistance. Nevertheless, the combined outcomes of these studies suggest that host resistance to bTB is a complex, polygenic trait, with no single gene alone explaining the inherent differences between resistant and susceptible animals. Such results support the development of accurate genomic breeding values that duly capture the collective effect of multiple genes to underpin selective breeding programmes. In addition to improving host resistance to bTB, scientists and practitioners have considered the possibility of reducing host infectivity. Ongoing studies have suggested the presence of genetic variation for infectivity and confirmed that bTB eradication would be accelerated if selective breeding considered both host resistance and infectivity traits. In conclusion, research activity on bTB genetics has generated knowledge and insights to support selective breeding as an additional measure towards controlling and eradicating the disease.

Genomic diversity and relationship analyses of endangered German Black Pied cattle (DSN) to 68 other taurine breeds based on whole-genome sequencing

German Black Pied cattle (Deutsches Schwarzbuntes Niederungsrind, DSN) are an endangered dual-purpose cattle breed originating from the North Sea region. The population comprises about 2,500 cattle and is considered one of the ancestral populations of the modern Holstein breed. The current study aimed at defining the breeds closest related to DSN cattle, characterizing their genomic diversity and inbreeding. In addition, the detection of selection signatures between DSN and Holstein was a goal. Relationship analyses using fixation index (F_ST), phylogenetic, and admixture analyses were performed between DSN and 68 other breeds from the 1000 Bull Genomes Project. Nucleotide diversity, observed heterozygosity, and expected heterozygosity were calculated as metrics for genomic diversity. Inbreeding was measured as excess of homozygosity (F_Hom) and genomic inbreeding (F_RoH) through runs of homozygosity (RoHs). Region-wide F_ST and cross-population-extended haplotype homozygosity (XP-EHH) between DSN and Holstein were used to detect selection signatures between the two breeds, and RoH islands were used to detect selection signatures within DSN and Holstein. DSN showed a close genetic relationship with breeds from the Netherlands, Belgium, Northern Germany, and Scandinavia, such as Dutch Friesian Red, Dutch Improved Red, Belgian Red White Campine, Red White Dual Purpose, Modern Angler, Modern Danish Red, and Holstein. The nucleotide diversity in DSN (0.151%) was higher than in Holstein (0.147%) and other breeds, e.g., Norwegian Red (0.149%), Red White Dual Purpose (0.149%), Swedish Red (0.149%), Hereford (0.145%), Angus (0.143%), and Jersey (0.136%). The F_Hom and F_RoH values in DSN were among the lowest. Regions with high F_ST between DSN and Holstein, significant XP-EHH regions, and RoH islands detected in both breeds harbor candidate genes that were previously reported for milk, meat, fertility, production, and health traits, including one QTL detected in DSN for endoparasite infection resistance. The selection signatures between DSN and Holstein provide evidence of regions responsible for the dual-purpose properties of DSN and the milk type of Holstein. Despite the small population size, DSN has a high level of diversity and low inbreeding. F_ST supports its relatedness to breeds from the same geographic origin and provides information on potential gene pools that could be used to maintain diversity in DSN.

Genome-wide analysis of runs of homozygosity in Italian Mediterranean buffalo

Runs of homozygosity (ROH) are a powerful tool to explore patterns of genomic inbreeding in animal populations and detect signatures of selection. The present study used ROH analysis to evaluate the genome-wide patterns of homozygosity, inbreeding levels, and distribution of ROH islands using the SNP data sets from 899 Mediterranean buffaloes. A total of 42,433 ROH segments were identified, with an average of 47.20 segments per individual. The ROH comprising mostly shorter segments (1-4 Mb) accounted for approximately 72.29% of all ROH. In contrast, the larger ROH (>8 Mb) class accounted for only 7.97% of all ROH segments. Estimated inbreeding coefficients from ROH (F_ROH) ranged from 0.0201 to 0.0371. Pearson correlations between F_ROH and genomic relationship matrix increased with the increase of ROH length. We identified ROH hotspots in 12 genomic regions, located on chromosomes 1, 2, 3, 5, 17, and 19, harboring a total of 122 genes. Protein-protein interaction (PPI) analysis revealed the clustering of these genes into 7 PPI networks. Many genes located in these regions were associated with different production traits. In addition, 5 ROH islands overlapped with cattle quantitative trait loci that were mainly associated with milk traits. These findings revealed the genome-wide autozygosity patterns and inbreeding levels in Mediterranean buffalo. Our study identified many candidate genes related to production traits that could be used to assist in selective breeding for genetic improvement of buffalo.

Host Serum Proteins as Potential Biomarkers of Bovine Tuberculosis Resistance Phenotype

Eradication of bovine tuberculosis (bTB) continues to be a worldwide challenge. The lack of reliable vaccines dampens the control and eradication programs of Mycobacterium bovis infection and spread. Selection and breeding of cattle resistant to M. bovis infection would greatly enhance the effectiveness of bTB eradication programs. Here, we have evaluated the potential of serum proteins as biomarkers of cattle resistance to bTB in Holstein-Friesian cows, 6-8-year-old, born and raised in similar conditions in herds with bTB prevalence >30%. Serum proteins obtained from uninfected cows (bTB-resistant; R) were compared to those from infected cows (bTB-susceptible; S), defined by a negative or positive bTB diagnosis, respectively. bTB diagnosis included: (i) single intradermal (caudal fold) tuberculin test, (ii) whole blood IFN-gamma test, (iii) gross visible lesions in lymph nodes and lungs by inspection at the abattoir, and (iv) a bacteriological culture for M. bovis. Using 2D-GE and LC-ESI-MS/MS, we found higher expression levels of primary amine oxidase (AO), complement component 5 (C5), and serotransferrin (TF) in R cattle than S cattle. In-house developed and standardized ELISAs for these novel biomarkers showed the best sensitivities of 72, 77, 77%, and specificities of 94, 94, 83%, for AO, C5, and TF, respectively. AUC-ROC (95% CI) values of 0.8935 (0.7906-0.9964), 0.9290 (0.8484-1.010), and 0.8580 (0.7291-0.9869) were obtained at cut-off points of 192.0, 176.5 ng/ml, and 2.1 mg/ml for AO, C5, and TF, respectively. These proteins are involved in inflammatory/immunomodulatory responses to infections and may provide a novel avenue of research to determine the mechanisms of protection against bTB. Overall, our results indicate that these proteins could be novel biomarkers to help identify cattle resistant to bTB, which in turn could be used to strengthen the effectiveness of existing eradication programs against bTB.

Taurine and Indicine Haplotype Representation in Advanced Generation Individuals From Three American Breeds

Development of the American Breeds of beef cattle began in the 1920s as breeders and U. S. Experiment Station researchers began to create Bos taurus taurus × Bos taurus indicus hybrids using Brahman as the B. t. indicus source. By 1954, U.S. Breed Associations had been formed for Brangus (5/8 Angus × 3/8 Brahman), Beefmaster (½ Brahman × ¼ Shorthorn × ¼ Hereford), and Santa Gertrudis (5/8 Shorthorn × 3/8 Brahman). While these breeds were developed using mating designs expected to create base generation animals with the required genome contributions from progenitor breeds, each association has now registered advanced generation animals in which selection or drift may have caused the realized genome compositions to differ from initial expected proportions. The availability of high-density SNP genotypes for 9,161 Brangus, 3,762 Beefmaster, and 1,942 Santa Gertrudis animals allowed us to compare the realized genomic architectures of breed members to the base generation expectations. We used RFMix to estimate local ancestry and identify genomic regions in which the proportion of Brahman ancestry differed significantly from a priori expectations. For all three breeds, lower than expected levels of Brahman composition were found genome-wide, particularly in early-generation animals where we demonstrate that selection on beef production traits was likely responsible for the taurine enrichment. Using a proxy for generation number, we also contrasted the genomes of early- and advanced-generation animals and found that the indicine composition of the genome has increased with generation number likely due to selection on adaptive traits. Many of the most-highly differentiated genomic regions were breed specific, suggesting that differences in breeding objectives and selection intensities exist between the breeds. Global ancestry estimation is commonly performed in admixed animals to control for stratification in association studies. However, local ancestry estimation provides the opportunity to investigate the evolution of specific chromosomal segments and estimate haplotype effects on trait variation in admixed individuals. Investigating the genomic architecture of the American Breeds not only allows the estimation of indicine and taurine genome proportions genome-wide, but also the locations within the genome where either taurine or indicine alleles confer a selective advantage.

Identification of loci associated with pathological outcomes in Holstein cattle infected with Mycobacterium avium subsp. paratuberculosis using whole-genome sequence data

Bovine paratuberculosis (PTB), caused by Mycobacterium avium subsp. paratuberculosis (MAP), is a chronic granulomatous enteritis that affects cattle worldwide. According to their severity and extension, PTB-associated histological lesions have been classified into the following groups; focal, multifocal, and diffuse. It is unknown whether these lesions represent sequential stages or divergent outcomes. In the current study, the associations between host genetic and pathology were explored by genotyping 813 Spanish Holstein cows with no visible lesions (N = 373) and with focal (N = 371), multifocal (N = 33), and diffuse (N = 33) lesions in gut tissues and regional lymph nodes. DNA from peripheral blood samples of these animals was genotyped with the bovine EuroG MD Bead Chip, and the corresponding genotypes were imputed to whole-genome sequencing (WGS) data using the 1000 Bull genomes reference population. A genome-wide association study (GWAS) was performed using the WGS data and the presence or absence of each type of histological lesion in a case–control approach. A total of 192 and 92 single nucleotide polymorphisms (SNPs) defining 13 and 9 distinct quantitative trait loci (QTLs) were highly-associated (P ≤ 5 × 10⁻⁷) with the multifocal (heritability = 0.075) and the diffuse (heritability = 0.189) lesions, respectively. No overlap was seen in the SNPs controlling these distinct pathological outcomes. The identified QTLs overlapped with some QTLs previously associated with PTB susceptibility, bovine tuberculosis susceptibility, clinical mastitis, somatic cell score, bovine respiratory disease susceptibility, tick resistance, IgG level, and length of productive life. Pathway analysis with candidate genes overlapping the identified QTLs revealed a significant enrichment of the keratinization pathway and cholesterol metabolism in the animals with multifocal and diffuse lesions, respectively. To test whether the enrichment of SNP variants in candidate genes involved in the cholesterol metabolism was associated with the diffuse lesions; the levels of total cholesterol were measured in plasma samples of cattle with focal, multifocal, or diffuse lesions or with no visible lesions. Our results showed reduced levels of plasma cholesterol in cattle with diffuse lesions. Taken together, our findings suggested that the variation in MAP-associated pathological outcomes might be, in part, genetically determined and indicative of distinct host responses.

Identification of loci associated with susceptibility to Mycobacterium avium subsp. paratuberculosis infection in Holstein cattle using combinations of diagnostic tests and imputed whole-genome sequence data

Bovine paratuberculosis (PTB) is a chronic inflammatory disease caused by Mycobacterium avium susbp. paratuberculosis (MAP). Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) significantly associated with susceptibility to bovine PTB. The main objective of this study was to identify quantitative trait loci (QTLs) associated with MAP infection in Spanish Holstein cows (N = 983) using combinations of diagnostic tests and imputed whole-genome sequence (WGS) data. The infection status of these animals was defined by three diagnostic methods including ELISA for MAP-antibodies detection, and tissue culture and PCR for MAP detection. The 983 cows included in this study were genotyped with the Bovine MD SNP50 Bead Chip, and the corresponding genotypes were imputed to WGS using the 1,000 Bull genomes reference population. In total, 33.77 million SNP variants per animal were identified across the genome. Linear mixed models were used to calculate the heritability (h²) estimates for each diagnostic test and test combinations. Next, we performed a case-control GWAS using the imputed WGS datasets and the phenotypes and combinations of phenotypes with h² estimates > 0.080. After performing the GWAS, the test combinations that showed SNPs with a significant association (P_FDR ≤ 0.05), were the ELISA-tissue PCR-tissue culture, ELISA-tissue culture, and ELISA-tissue PCR. A total of twelve quantitative trait loci (QTLs) highly associated with MAP infection status were identified on the Bos taurus autosomes (BTA) 4, BTA5, BTA11, BTA12, BTA14, BTA23, BTA24, and BTA28, and some of these QTLs were linked to immune-modulating genes. The identified QTLs on BTA23 spanning from 18.81 to 22.95 Mb of the Bos taurus genome overlapped with several QTLs previously found to be associated with PTB susceptibility, bovine tuberculosis susceptibility, and clinical mastitis. The results from this study provide more clues regarding the molecular mechanisms underlying susceptibility to PTB infection in cattle and might be used to develop national genetic evaluations for PTB in Spain.

Ancient autozygous segments subject to positive selection suggest adaptive immune responses in West African cattle

Small-sized and trypanotolerant West African taurine (Bos taurus) cattle are a unique case of human-mediated process of adaptation to a challenging environment. Extensive gene flow with Sahelian zebu (B. indicus), bigger and with some resistance to tick attack, occurred for centuries and allowed the apparition of stable crossbred populations (sanga) having intermediate characteristics. Up to 237 individuals belonging to 10 different taurine, zebu and sanga cattle populations sampled in Benin, Burkina Faso and Niger were typed using the BovineHD BeadChip of Illumina to identify signatures of selection, assessed using three different Extended-Haplotype-Homozygosity-based statistics, overlapping with ancient, originated 1024 or 2048 generations ago, Homozygosity-By-Descent segments in the cattle genome. Candidate genomic regions were defined ensuring their importance within cattle type and using zebu as reference. Functional annotation analysis identified four statistically significant Annotation Clusters in taurine cattle (from ACt1 to ACt4), one (ACs1) in sanga, and another (ACz1) in zebu cattle, fitting well with expectations. ACt1 included genes primarily associated with innate immunity; ACt2 involved bitter taste receptor genes of importance to adaptation to changing environments; ACt3 included 68 genes coding ATP-binding proteins, some of them located on trypanotolerance-related QTL regions, that can partially underlie immune response and the additive mechanism of trypanotolerance; ACt4 was associated with growth and small size (NPPC gene); ACs1 included genes involved in immune response; and ACz1 is related with ectoparasite resistance. Our results provide a new set of genomic areas and candidate genes giving new insights on the genomic impact of adaptation in West African cattle.

Molecular epidemiology and drug susceptibility profiles of Mycobacterium tuberculosis complex isolates from Northern Ghana

OBJECTIVE

We conducted a cross-sectional study in the five administrative regions of Northern Ghana to determine the diversity of Mycobacterium tuberculosis complex (MTBC) sub/lineages and their susceptibility to isoniazid (INH) and rifampicin (RIF).

METHODS

Sputum specimens were collected and cultured from 566 pulmonary tuberculosis patients reporting to 17 health facilities from 2015 to 2019. Mycobacterial isolates obtained from solid cultures were confirmed as members of the MTBC by PCR amplification of IS6110 and rpoß and assigned lineages and sub-lineages using spoligotyping.

RESULTS

Of 294 mycobacterial isolates recovered, MTBC species identified were: M. tuberculosis sensu stricto (Mtbss) 241 (82.0%), M. africanum 41 (13.9%) and M. bovis four (1.4%) with eight (2.7%) unidentified. The human-adapted lineages (L) identified (N=279) were L1 (8/279, 2.9%), L2 (15/279, 5.4%), L3 (7/279, 2.5%), L4 (208/279, 74.5%), L5 (13/279, 4.7%) and L6 (28/279, 10.0%) with three unidentified lineages. Among the 208 L4, the dominant sub-lineages in the region were the Cameroon 120/208 (57.7%) and Ghana 50/208 (24.0%). We found 4.4% (13/294) and 0.7% (2/294) of the patients infected with MTBC isolates resistant to INH only and RIF only, respectively, with 2.4% (7/294) being infected with MDR strains. Whereas L6 was associated with the elderly, we identified that the Ghana sub-lineage of L4 was associated with both INH and MDR (p<0.05), making them important TB pathogens in Northern Ghana and a growing public health concern.

Mycobacterial Infection of Precision-Cut Lung Slices Reveals Type 1 Interferon Pathway Is Locally Induced by Mycobacterium bovis but Not M. tuberculosis in a Cattle Breed

Tuberculosis exacts a terrible toll on human and animal health. While Mycobacterium tuberculosis (Mtb) is restricted to humans, Mycobacterium bovis (Mb) is present in a large range of mammalian hosts. In cattle, bovine TB (bTB) is a noticeable disease responsible for important economic losses in developed countries and underestimated zoonosis in the developing world. Early interactions that take place between mycobacteria and the lung tissue early after aerosol infection govern the outcome of the disease. In cattle, these early steps remain poorly characterized. The precision-cut lung slice (PCLS) model preserves the structure and cell diversity of the lung. We developed this model in cattle in order to study the early lung response to mycobacterial infection. In situ imaging of PCLS infected with fluorescent Mb revealed bacilli in the alveolar compartment, in adjacent or inside alveolar macrophages, and in close contact with pneumocytes. We analyzed the global transcriptional lung inflammation signature following infection of PCLS with Mb and Mtb in two French beef breeds: Blonde d'Aquitaine and Charolaise. Whereas, lungs from the Blonde d'Aquitaine produced high levels of mediators of neutrophil and monocyte recruitment in response to infection, such signatures were not observed in the Charolaise in our study. In the Blonde d'Aquitaine lung, whereas the inflammatory response was highly induced by two Mb strains, AF2122 isolated from cattle in the UK and Mb3601 circulating in France, the response against two Mtb strains, H37Rv, the reference laboratory strain, and BTB1558, isolated from zebu in Ethiopia, was very low. Strikingly, the type I interferon pathway was only induced by Mb but not Mtb strains, indicating that this pathway may be involved in mycobacterial virulence and host tropism. Hence, the PCLS model in cattle is a valuable tool to deepen our understanding of early interactions between lung host cells and mycobacteria. It revealed striking differences between cattle breeds and mycobacterial strains. This model could help in deciphering biomarkers of resistance vs. susceptibility to bTB in cattle as such information is still critically needed for bovine genetic selection programs and would greatly help the global effort to eradicate bTB.

Integrative genomics of the mammalian alveolar macrophage response to intracellular mycobacteria

Background

Bovine TB (bTB), caused by infection with Mycobacterium bovis, is a major endemic disease affecting global cattle production. The key innate immune cell that first encounters the pathogen is the alveolar macrophage, previously shown to be substantially reprogrammed during intracellular infection by the pathogen. Here we use differential expression, and correlation- and interaction-based network approaches to analyse the host response to infection with M. bovis at the transcriptome level to identify core infection response pathways and gene modules. These outputs were then integrated with genome-wide association study (GWAS) data sets to enhance detection of genomic variants for susceptibility/resistance to M. bovis infection.

Results

The host gene expression data consisted of RNA-seq data from bovine alveolar macrophages (bAM) infected with M. bovis at 24 and 48 h post-infection (hpi) compared to non-infected control bAM. These RNA-seq data were analysed using three distinct computational pipelines to produce six separate gene sets: 1) DE genes filtered using stringent fold-change and P-value thresholds (DEG-24: 378 genes, DEG-48: 390 genes); 2) genes obtained from expression correlation networks (CON-24: 460 genes, CON-48: 416 genes); and 3) genes obtained from differential expression networks (DEN-24: 339 genes, DEN-48: 495 genes). These six gene sets were integrated with three bTB breed GWAS data sets by employing a new genomics data integration tool—gwinteR. Using GWAS summary statistics, this methodology enabled detection of 36, 102 and 921 prioritised SNPs for Charolais, Limousin and Holstein-Friesian, respectively.

Conclusions

The results from the three parallel analyses showed that the three computational approaches could identify genes significantly enriched for SNPs associated with susceptibility/resistance to M. bovis infection. Results indicate distinct and significant overlap in SNP discovery, demonstrating that network-based integration of biologically relevant transcriptomics data can leverage substantial additional information from GWAS data sets. These analyses also demonstrated significant differences among breeds, with the Holstein-Friesian breed GWAS proving most useful for prioritising SNPS through data integration. Because the functional genomics data were generated using bAM from this population, this suggests that the genomic architecture of bTB resilience traits may be more breed-specific than previously assumed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07643-w.

Cross-sectional analyses of a national database to determine if superior genetic merit translates to superior dairy cow performance

Various studies have validated that genetic divergence in dairy cattle translates to phenotypic differences; nonetheless, many studies that consider the breeding goal, or associated traits, have generally been small scale, often undertaken in controlled environments, and they lack consideration for the entire suite of traits included in the breeding goal. Therefore, the objective of the present study was to fill this void, and in doing so, provide producers with confidence that the estimated breeding values (EBV) included in the breeding goal do (or otherwise) translate to desired changes in performance among commercial cattle; an additional outcome of such an approach is the identification of potential areas for improvements. Performance data on 536,923 Irish dairy cows (and their progeny) from 13,399 commercial spring-calving herds were used. Association analyses between the cow's EBV of each trait included in the Irish total merit index for dairy cows (which was derived before her own performance data accumulated) and her subsequent performance were undertaken using linear mixed models; milk production, fertility, calving, maintenance (i.e., liveweight), beef, health, and management traits were all considered in the analyses. Results confirm that excelling in EBV for individual traits, as well as on the total merit index, generally delivers superior phenotypic performance; examples of the improved performance for genetically elite animals include a greater yield and concentration of both milk fat and milk protein, despite a lower milk volume, superior reproductive performance, better survival, improved udder and hoof health, lighter cows, and fewer calving complications; all these gains were achieved with minimal to no effect on the beef merit of the dairy cow's progeny. The associated phenotypic change in each performance trait per unit change in its respective EBV was largely in line with the direction and magnitude of expectation, the exception being for calving interval. Per unit change in calving interval EBV, the direction of phenotypic response was as anticipated but the magnitude of the response was only half of what was expected. Despite the deviation from expectation between the calving interval EBV and its associated phenotype, a superior total merit index or a superior fertility EBV was indeed associated with an improvement in all detailed fertility performance phenotypes investigated. Results substantiate that breeding is a sustainable strategy of improving phenotypic performance in commercial dairy cattle and, by extension, profit.

Genetic Diversity of Cameroon Cattle and a Putative Genomic Map for Resistance to Bovine Tuberculosis

Bovine Tuberculosis (bTB) caused by Mycobacterium bovis is a livestock disease of global economic and public health importance. There are currently no effective vaccines available for livestock and so control relies on animal level surveillance and pasteurization of dairy products. A new alternative control approach is to exploit the genetic variability of the host; recent studies have demonstrated that breeding European taurine cattle, such as Holsteins for increased resistance to bTB is feasible. The utility of such an approach is still unknown for African cattle populations. This study aims to assess genetic variation in bTB resistance and the underlying genomic architecture of cattle in Cameroon. We conducted a cross-sectional study of 2,346 slaughter cattle in Cameroon. Retropharyngeal lymph node samples were collected and cultured on Lowenstein Jensen media and the BACTEC MGIT 960 system, and M. bovis was identified using the Hain® Genotype kits. A total of 153 cattle were positive for M. bovis and were archived along with a random selection of negative samples. In this study, we genotyped archived samples from 212 cattle. Their genomic diversity was characterized using PCA, hierarchical clustering and admixture analysis. We assessed genetic variation in bTB resistance using heritability analysis and compared quantitative trait loci. Previous research on this study population have shown that Fulani cattle are more susceptible to bTB than mixed breeds. However, here we show that these apparent phenotypic differences in breeds are not reflected by clear genomic differences. At the genetic level, both the Fulani and mixed cattle show similar patterns of admixture with evidence of both taurine and indicine ancestry. There was little European taurine introgression within the studied population. Hierarchical clustering showed clusters of cattle that differed in their susceptibility to bTB. Our findings allude to bTB resistance being polygenic in nature. This study highlights the potential for genetic control of bTB in Africa and the need for further research into the genetics of bTB resistance within African cattle populations.

The impact of human single nucleotide polymorphisms on Bacillus Calmette-Guérin responses

The influence of genetic variability on human immune responses has major implications for the understanding of disease mechanisms and host-pathogen interactions. Bacillus Calmette-Guérin (BCG) vaccine, which is given globally to protect against tuberculosis, has high variability in its protective efficacy against mycobacteria and its beneficial off-target (heterologous) effects. Single nucleotide polymorphisms (SNPs) are major cause of genetic variation and have been strongly associated with susceptibility to tuberculosis and outcomes following BCG immunotherapy for cancer. This review discusses the contribution of SNPs to the variability in mycobacterial-specific and off-target BCG responses, and the implications for this on development of novel TB vaccines and strategies to harness the beneficial off-target effects of BCG.

Analysis of Genetic Variation in the Bovine SLC11A1 Gene, Its Influence on the Expression of NRAMP1 and Potential Association With Resistance to Bovine Tuberculosis

Bovine tuberculosis (bTB), caused by Mycobacterium bovis, is a chronic zoonotic disease where host genetics is thought to contribute to susceptibility or resistance. One of the genes implicated is the SLC11A1 gene, that encodes for the natural resistance-associated macrophage protein 1 (NRAMP1). The aim of this study was to identify SLC11A1 polymorphisms and to investigate any resulting functional differences in NRAMP1 expression that might be correlated with resistance/susceptibility to M. bovis infection. Sequencing of the SLC11A1 gene in cDNA isolated from Brown Swiss, Holstein Friesian, and Sahiwal cattle identified five single nucleotide polymorphisms (SNPs) in the coding region, but only one of these (SNP4, c.1066C>G, rs109453173) was present in all three cattle breeds and therefore warranted further investigation. Additionally, variations of 10, 11, and 12 GT repeats were identified in a microsatellite (MS1) in the SLC11A1 3′UTR. Measurement of NRAMP1 expression in bovine macrophages by ELISA showed no differences between cells generated from the different breeds. Furthermore, variations in the length of the MS1 microsatellite did not impact on NRAMP1 protein expression as analyzed by luciferase reporter assay. However, further analysis of the ELISA data identified that the presence of the alternative G allele at SNP4 was associated with increased expression of NRAMP1 in bovine macrophages. Since NRAMP1 has been shown to influence the survival of intracellular pathogens such as M. bovis through the sequestering of iron, it is possible that cattle expressing the alternative G allele might have an increased resistance to bTB through increased NRAMP1 expression in their macrophages.

Accuracy of whole-genome sequence imputation using hybrid peeling in large pedigreed livestock populations

BACKGROUND

The coupling of appropriate sequencing strategies and imputation methods is critical for assembling large whole-genome sequence datasets from livestock populations for research and breeding. In this paper, we describe and validate the coupling of a sequencing strategy with the imputation method hybrid peeling in real animal breeding settings.

METHODS

We used data from four pig populations of different size (18,349 to 107,815 individuals) that were widely genotyped at densities between 15,000 and 75,000 markers genome-wide. Around 2% of the individuals in each population were sequenced (most of them at 1× or 2× and 37-92 individuals per population, totalling 284, at 15-30×). We imputed whole-genome sequence data with hybrid peeling. We evaluated the imputation accuracy by removing the sequence data of the 284 individuals with high coverage, using a leave-one-out design. We simulated data that mimicked the sequencing strategy used in the real populations to quantify the factors that affected the individual-wise and variant-wise imputation accuracies using regression trees.

RESULTS

Imputation accuracy was high for the majority of individuals in all four populations (median individual-wise dosage correlation: 0.97). Imputation accuracy was lower for individuals in the earliest generations of each population than for the rest, due to the lack of marker array data for themselves and their ancestors. The main factors that determined the individual-wise imputation accuracy were the genotyping status, the availability of marker array data for immediate ancestors, and the degree of connectedness to the rest of the population, but sequencing coverage of the relatives had no effect. The main factors that determined variant-wise imputation accuracy were the minor allele frequency and the number of individuals with sequencing coverage at each variant site. Results were validated with the empirical observations.

CONCLUSIONS

We demonstrate that the coupling of an appropriate sequencing strategy and hybrid peeling is a powerful strategy for generating whole-genome sequence data with high accuracy in large pedigreed populations where only a small fraction of individuals (2%) had been sequenced, mostly at low coverage. This is a critical step for the successful implementation of whole-genome sequence data for genomic prediction and fine-mapping of causal variants.

Candidate genes associated with the heritable humoral response to Mycobacterium avium ssp. paratuberculosis in dairy cows have factors in common with gastrointestinal diseases in humans

Infection of cattle with bovine paratuberculosis (i.e., Johne's disease) is caused by Mycobacterium avium ssp. paratuberculosis (MAP) and results in a chronic incurable gastroenteritis. This disease, which has economic ramifications for the cattle industry, is increasing in detected prevalence globally; subclinically infected animals can silently shed the bacterium into the environment for years, exposing contemporaries and hampering disease-control programs. The objective of the present study was to first quantify the genetic parameters for humoral response to MAP in dairy cattle. This was followed by a genome-based association analysis and subsequent downstream bioinformatic analyses from imputed whole genome sequence SNP data. After edits, ELISA test records were available on 136,767 cows; analyses were also undertaken on a subset of 33,818 of these animals from herds with at least 5 MAP ELISA-positive cows, with at least 1 of those positive cows being homebred. Variance components were estimated using univariate animal and sire linear mixed models. The heritability calculated from the animal model for humoral response to MAP using alternative phenotype definitions varied from 0.02 (standard error = 0.003) to 0.05 (standard error = 0.008). The genome-based associations were undertaken within a mixed model framework using weighted deregressed estimated breeding values as a dependent variable on 1,883 phenotyped animals that were ≥87.5% Holstein-Friesian. Putative susceptibility quantitative trait loci (QTL) were identified on Bos taurus autosome 1, 3, 5, 6, 8, 9, 10, 11, 13, 14, 18, 21, 23, 25, 26, 27, and 29; mapping the most significant SNP to genes within and overlapping these QTL revealed that the most significant associations were with the 10 functional candidate genes KALRN, ZBTB20, LPP, SLA2, FI3A1, LRCH3, DNAJC6, ZDHHC14, SNX1, and HAS2. Pathway analysis failed to reveal significantly enriched biological pathways, when both bovine-specific pathway data and human ortholog data were taken into account. The existence of genetic variation for MAP susceptibility in a large data set of dairy cows signifies the potential of breeding programs for reducing MAP susceptibility. Furthermore, the identification of susceptible QTL facilitates greater biological understanding of bovine paratuberculosis and potential therapeutic targets for future investigation. The novel molecular similarities identified between bovine paratuberculosis and human inflammatory bowel disease suggest potential for human therapeutic interventions to be translated to veterinary medicine and vice versa.