Open Access

Identification of single nucleotide polymorphisms in Toll-like receptor candidate genes associated with tuberculosis infection in water buffalo (Bubalus bubalis)

  • Flora Alfano1Email author,
  • Simone Peletto2,
  • Maria Gabriella Lucibelli1,
  • Giorgia Borriello1,
  • Giovanna Urciuolo1,
  • Maria Grazia Maniaci2,
  • Rosanna Desiato2,
  • Michela Tarantino3,
  • Amalia Barone4,
  • Paolo Pasquali3,
  • Pier Luigi Acutis2 and
  • Giorgio Galiero1
BMC Genetics201415:139

https://doi.org/10.1186/s12863-014-0139-y

Received: 12 March 2014

Accepted: 27 November 2014

Published: 14 December 2014

Abstract

Background

Toll-like receptors play a key role in innate immunity by recognizing pathogens and activating appropriate responses. Pathogens express several signal molecules (pathogen-associated molecular patterns, PAMPs) essential for survival and pathogenicity. Recognition of PAMPs triggers an array of anti-microbial immune responses through the induction of various inflammatory cytokines. The objective of this work was to perform a case-control study to characterize the distribution of polymorphisms in three candidate genes (toll-like receptor 2, toll-like receptor 4, toll-like receptor 9) and to test their role as potential risk factors for tuberculosis infection in water buffalo (Bubalus bubalis).

Results

The case-control study included 184 subjects, 59 of which resulted positive to both intradermal TB test and Mycobacterium bovis isolation (cases) and 125 resulted negative to at least three consecutive intradermal TB tests. The statistical analysis indicated that two polymorphisms exhibited significant differences in allelic frequencies between cases and controls. Indeed, the TT genotype at TLR9 2340 C > T locus resulted significantly associated with susceptibility to bovine tuberculosis (P = 0.030, OR = 3.31, 95% CI = 1.05-10.40). One polymorphism resulted significantly associated with resistance to the disease, and included the CC genotype, at the TLR4 672 A > C locus (P = 0.01, OR = 0.26, 95% CI = 0.08-0.80). Haplotype reconstruction of the TLR2 gene revealed one haplotype (CTTACCAGCGGCCAGTCCC) associated with disease resistance (P = 0.04, OR = 0.51, 95% CI = 0.27–0.96), including the allelic variant associated with disease resistance.

Conclusions

The work describes novel mutations in bubaline TLR2, TLR4 and TLR9 genes and presents their association with M. bovis infection. These results will enhance our ability to determine the risk of developing the disease by improving the knowledge of the immune mechanisms involved in host response to mycobacterial infection, and will allow the creation of multiple layers of disease resistance in herds by selective breeding.

Keywords

Bubalus bubalis TLRs Genetic resistance Case–control study

Background

Tuberculosis is a major zoonosis that causes serious economic losses in livestock industry. Among domestic ruminants, an economically important role, in many parts of the world and particularly in Southern Italy, is played by water buffalo (Bubalus bubalis). Milk from this animal species is used for the production of the worldwide famous “mozzarella di bufala” cheese, whose product specifications allow the use of raw milk [1]. It is therefore mandatory the use of milk collected from brucellosis- and tuberculosis-free herds. Tuberculosis in Southern Italy has not yet been eradicated (prevalence of 0.7% in bovine and water buffalo herds in 2012), and there is an increase of new cases of the disease (incidence of 0.65%) (data from the veterinary epidemiology monitoring field office, Istituto Zooprofilattico Sperimentale del Mezzogiorno www.izsmportici.it). The Italian National plan for tuberculosis control is based on a test-and-slaughter approach and does not allow vaccination. The genetic variation within the host may play a crucial role in the immunity to infections and resistance or susceptibility to disease, and selective breeding for disease-resistant genotypes might represent an emerging approach supporting disease control [2],[3]. The role of genes in protection against bacterial infections also shown in water buffalo [4] prompted the search for polymorphisms conferring resistance to mycobacterial infection in this species.

Microbial infection initiates complex interactions between the pathogen and the host. Pathogens express several molecules, known as pathogen-associated molecular patterns (PAMPs), which are essential for survival and pathogenicity. Recognition of PAMPs triggers an array of anti-microbial immune responses through the induction of various inflammatory cytokines [5].

Toll-like receptors (TLRs) are reported to be involved in immune responses [6]-[8] and their polymorphisms have been associated with mutated susceptibility to mycobacteria, the causative agent of bovine tuberculosis, both in humans and different animal species [9]-[15]. TLRs are transmembrane proteins that play a key role in innate immunity by recognizing pathogens and subsequently activating appropriate responses. They are characterized by an extracellular N-terminal domain constituted by 16 to 28 leucine rich repeats (LRR) involved in ligand recognition [16], and an intracellular C-terminal domain known as the toll/IL-1 receptor (TIR) domain, required for the interaction and recruitment of various adaptor molecules to activate the downstream signalling pathways [17]. In water buffalo, each cell contains 25 chromosome pairs (Bubalus bubalis 2n = 50) [18], and the TLRs 2 and 9 loci are localized on chromosomes 17 and 21, respectively, while the TLR 4 locus is located on chromosome 3 [19].

Toll-like receptors 2 and 4 mediate their effects by recognizing bacteria, including mycobacteria. Toll-like receptor 2 is sensitive to peptidoglycan and to lipoarabinomannan. Toll-like receptor 4 is involved in bacterial lipopolysaccharide recognition [20], but its sensitivity to mycobacterial antigens has also been reported [21]. Toll-like receptor 9 is essential for responses to bacterial DNA and in particular to unmethylated CpG dinucleotides (CpG DNA) [22].

The aim of this study was (i) to identify single nucleotide polymorphisms (SNPs) in the coding sequence (CDS) of three candidate genes: toll-like receptor 2 (TLR 2), toll-like receptor 4 (TLR 4), toll-like receptor 9 (TLR 9), (ii) to characterize the distribution of the detected polymorphisms and (iii) to perform a case-control study to test their role as potential risk factors for tuberculosis infection in water buffalo.

Results

Identification of SNPs

The water buffalo TLR2 and TLR9 genes consist of two exons, while the TLR4 consists of three exons. The analysis of the entire CDSs enabled us to identify 29 new SNPs. Specifically, 18 were identified in TLR2 (five non-synonymous and 13 synonymous), nine in TLR4 (five non-synonymous and four synonymous), and two in TLR9 (both synonymous). All the found SNPs were bi-allelic with the exception of a tri-allelic SNP (572 A > CG) found in TLR4. In addition we identified a dinucleotidic SNP (482/483 GC > CT) in TLR2. Table 1 shows all the SNPs, their positions in the CDS, the encoded protein codon and the GenBank reference number.
Table 1

Detected SNPs in bubaline TLRs2 , 4 and 9 genes

Gene

SNP a

AA change b

Protein domain c

dbSNP ID

TLR2

42 C > T

Silent 14

- d

rs1388116474:C > T

 

53 C > T

M18T

-

rs1388116475:C > T

 

108 C > T

Silent 36

LRR 8

rs1388116476:C > T

 

153 G > A

Silent 51

LRR_RI

rs1388116477:G > A

 

156 C > T

Silent 52

LRR_RI

rs1388116478:C > T

 

374 T > C

A125V

LRR 8

rs1388116479:T > C

 

381 A > G

Silent 127

LRR 8

rs1388116480:A > G

 

482/483 GC > CT

S161T

LRR 8

rs1388116481:GC > CT

 

519 G > C

Silent 173

LRR_RI

rs1388116482:G > C

 

1034 A > G

S345N

-

rs1388116483:A > G

 

1375 T > C

Silent 459

LRR4

rs1388116484:T > C

 

1407 C > T

Silent 469

LRR4

rs1388116485:C > T

 

1650 G > A

Silent 550

LRR_CT

rs1388116486:G > A

 

1678 A > G

A560T

LRR_CT

rs1388116487:A > G

 

1707 C > T e

Silent 569

LRR_CT

rs1388116488:C > T

 

1731 C > T

Silent 577

LRR_CT

rs1388116489:C > T

 

1740 C > T

Silent 580

LRR_CT

rs1388116490:C > T

 

2064 T > C

Silent 688

TIR

rs1388116491:T > C

TLR4

572 A > C f

Y191S

LRR8

rs1388116492:A > CG

 

572 A > G f

silent

LRR8

rs1388116492:A > CG

 

574 C > T

Q192W

LRR8

rs1388116493:C > T

 

575 A > G

Q192W

LRR8

rs1388116494:A > G

 

576 T > G

Silent192

LRR8

rs1388116495:T > G

 

577 G > A

E193K

LRR8

rs1388116496:G > A

 

579 A > G

Silent193

LRR8

rs1388116497:A > G

 

647 G > A

Silent216

LRR_RI

rs1388116498:G > A

 

662 G > A

Silent221

-

rs1388116499:G > A

 

672 A > C

Silent224

-

rs1388116500:A > C

TLR9

2340 C > T

Silent780

-

rs1388116501:C > T

 

2475 A > G

Silent825

-

rs1388116502:A > G

a SNPs positions were calculated taking the ATG start codon as position 1 based on the sequences GenBank: HM756161 (TLR2 gene), GenBank: JN786600 (TLR4 gene), GenBank: HQ242778 (TLR9 gene).

b Amino acid positions are given according to the ATG start codon.

c Protein domain has been predicted from the cds sequence by using the Conserved Domain Database available at http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi [23],[24].

d Protein position with unknown function.

e This SNP is conserved also in the bovine species [9].

f This polymorphic site exhibits a SNP with three different alleles (572 A > CG).

Case-control study

Calculation of the allelic and genotypic frequencies, and application of the χ2 test, revealed that both populations (cases and controls) conformed to the Hardy-Weinberg equilibrium with regard to all the polymorphic loci (P > 0.05) of the TLRs 2, 4 and 9 genes.

As a first step of the genetic association study, Pearson’s chi-squared test was carried out for each SNP both independently and as genotype. As many SNPs exhibited a strong linkage disequilibrium a multiple testing correction was applied by permutation analysis. Within the 29 analysed polymorphic sites, three SNPs exhibited statistically significant differences in frequency distribution of one or more associated genotypes between cases and controls (Table 2). In particular, one SNP was located in the TLR2 gene (1650 G > A ), one SNP in TLR4 (672 A > C) and one in TLR9 (2340 C > T). Then, in order to test inheritance model, a deviation test from additivity was performed with plink version 1.07, suggesting in our case an additive model of inheritance. SNPs significantly associated to the disease (P < 0.05) were used for further analysis: a logistic regression model was performed to calculate the Odds Ratio (OR) with the wild type category used as a baseline (Table 3).
Table 2

Results of permutation analysis

Chromosome

SNP

Empirical P -value

17

42-C > T

0.9000

17

53-C > T

0.1967

17

108-C > T

0.1714

17

153-G > A

0.2195

17

156-C > T

0.9452

17

374-T > C

0.1229

17

381-A > G

0.9788

17

482/483-GC > CT

0.8222

17

519-G > C

0.5175

17

1034-A > G

0.2191

17

1375-T > C

0.1302

17

1407-C > T

0.8738

17

1650-G > A

0.0138*

17

1678-A > G

0.1505

17

1707-C > T

0.2203

17

1731-C > T

0.9286

17

1740-C > T

0.8985

17

2064-T > C

0.1296

3

572 C > A G

0.4909

3

574 C > T

0.4288

3

575 A > G

0.4380

3

576 T > G

0.2490

3

577 G > A

0.6373

3

579 A > G

0.2551

3

647 G > A

0.1735

3

662 G > A

0.1006

3

672 A > C

0.0183*

21

2340 C > T

0.0521*

21

2475 A > G

0.1254

*statistically significant.

Table 3

Polymorphic sites including genotypes with statistically significant differences in frequency distribution between cases and controls

Gene

SNPs

Genotype

P -value

OR

CI 95%

TLR2

1650 G > A

G/A

0.166

0.612

0.26 - 1.01

  

A/A

0.021

0.000

--

TLR4

672 A > C

A/C

0.689

0.865

0.42 - 1.76

  

C/C

0.012

0.258

0.08 - 0.80

TLR9

2340 C > T

C/T

0.048

2.917

0.96 - 8.89

  

T/T

0.030

3.306

1.05 - 10.40

Bold lines highlight genotypes with statistically significant P-values and ORs.

In the TLR2 gene, the 1650 G > A locus exhibited a χ2 test of independence with a significant P-value (P = 0.018) while the OR for the AA genotype, calculated with the Woolf’s method, exhibited a not statistically significant value (P = 0.0512, OR = 0.10, 95% CI = 0.01-1.74). However, in this case, a significant test for trend (P = 0.01) in calculating OR values in logistic regression indicates a plausible strong correlation with resistance to the disease.

For the TLR4 gene, the SNP 672 A > C showed a χ2 test with a significant P-value (P = 0.031); the OR demonstrated a significant P-value (P = 0.012) for the CC genotype at the 672 A > C locus and a value lower than 1 (OR = 0.26, 95% CI = 0.08 to 0.80), thus indicating a significant association with resistance to tuberculosis.

In the TLR9 gene, the analysis revealed one polymorphic site (2340 C > T) with a significant P-value (0.053) at the χ2 test including one genotype (T/T) associated with susceptibility to the disease (P = 0.030 and OR = 3.31, 95% CI = 1.05-10.40).

Haplotype reconstruction based on TLR2 polymorphisms performed by PHASE software, generated 25 possible haplotypes, seven of which displayed a frequency greater than 0.01% in the whole sample. The software PHASE can reveal any significant imbalance in haplotype distribution between groups by performing a permutation test; the application of this analysis to the detected haplotypes revealed the presence of one haplotype (CTTACCAGCGGCCAGTCCC) associated with disease resistance (P = 0.041 and OR = 0.51, 95% CI = 0.27-0.96). This haplotype contained the allelic variant that resulted possibly associated with resistance to the disease.

In the TLR4 gene, 34 haplotypes were observed in the two tested groups, six of which displayed a frequency greater than 0.01% in the whole population. None of the predicted haplotypes showed a significant association with resistance/susceptibility to tuberculosis.

Finally, four haplotypes were observed in the TLR9 gene, none of which displayed a significant association with resistance/susceptibility to tuberculosis.

Discussion

The polymorphisms of single nucleotides constitute excellent genetic markers for various population studies as they can, for instance, reveal traces of natural selection. Moreover, it is now known that the contribution of genes to the incidence of or predisposition to diseases can be determined by comparing inter-individual genetic differences.

Each individual, to a greater or lesser degree than others, is genetically prone to develop certain diseases. In most cases, the susceptibility or resistance to a given pathology does not of course mean that the individual will necessarily be affected by the disease or remain disease-free; it does, however, indicate a higher or lower risk than that of the general population. Indeed, the onset of symptoms is due to an interaction between genetic factors and environmental factors. Nevertheless, cases have been described in which genetic factors are sufficient to determine the course of infection. In scrapie, for instance, which is a transmissible disease of sheep and goats caused by a prion, there are reports of genotypes that are defined as “resistant”, which are practically refractory to natural infection [25],[26].

Our analysis of the water buffalo TLR2, TLR4 and TLR9 genes enabled us to identify 29 SNPs. In the three genes, the 76% of the SNPs fell within the LRR protein domains, regions responsible for ligand recognition. This finding is in agreement with conservation of the transmembrane and functional TIR domains already observed in human, bovine and murine TLRs genes [9]-[13].

The case-control study carried out revealed the presence of associations between some of the identified polymorphisms and resistance/susceptibility to the disease caused by M. bovis. Specifically, in TLR4, we observed a correlation between resistance to disease and the genotype C/C at the 672 A > C locus (P = 0.012, OR = 0.26, 95% CI = 0.08-0.80) when a logistic regression model was performed. A likewise correlation with weak statistical evidence was seen in TLR2 at the 1650 A > G locus when a logistic regression model was performed between wild type genotype and mutated genotype (P = 0.0512, OR = 0.10, 95% CI = 0.01-1.74). Moreover, a correlation was seen between susceptibility to tuberculosis and the T/T genotype at the 2340 C > T locus in the TLR9 gene (P = 0.030, OR = 3.31, 95% CI = 1.05-10.40).

Haplotype reconstruction based on TLR2 polymorphisms enabled us to ascertain whether the interaction of the single alleles exerted an influence on the animal’s response to the disease; indeed, this analysis revealed a haplotype (CTTACCAGCGGCCAGTCCC) that was associated with disease resistance (P = 0.041 and OR = 0.51, 95% CI = 0.27-0.96). This haplotype contained the allelic variant which resulted possibly associated with resistance to the disease. However, as this variant is not exclusive of the associated haplotype, we may hypothesise that other loci exert a modulatory effect.

Haplotype reconstruction and the permutation test conducted on the polymorphisms identified in TLR4 did not retrieve any disease-associated haplotype.

Previous studies have demonstrated the influence of genetic factors and some gene polymorphisms on susceptibility to several bacterial diseases, including TB. TLRs have been found to be associated with both the innate immune and adaptive immune responses, and play a crucial role in the immune recognition of M. tuberculosis. Mutations and polymorphisms in TLRs and TLR signalling molecules have revealed the importance of TLRs both in human and animal defence against disease [14]-[22],[25]-[28]. One nsSNP in TLR2 (Arg753Gln) has been shown to increase human predisposition to staphylococcal infection, tuberculosis, rheumatic fever and urinary tract infection [10]-[15]. A nsSNP in the bovine TLR2 gene (Val220Met) has been shown to decrease the response to M. avium subsp. paratuberculosis (Map) [11], while another SNP in the same gene, C1903T (Silent569), was found to be associated with resistance to Map infection in cattle [9]. Several polymorphisms have been found in the bovine TLR4 gene, 12 SNPs affecting ligand binding domain [29] and four SNPs associated with susceptibility to Map infection [11]. Studies conducted on the bovine TLR9 gene demonstrated a high degree of genetic variability in this gene [14]-[22],[25]-[30], even though no significant association with bovine tuberculosis could be inferred from mutations in this gene [14]. A SNP in the human TLR9 gene (rs352139) was found to be associated with susceptibility to TB (P = 0.01, OR = 2.37) [31]. In our study, among the polymorphisms identified in water buffalo, only one (C1707T in the TLR2 gene) showed correspondence to a polymorphism reported in cattle (T1903C) [9]. Specifically, Koets and coworkers [9] found this SNP associated with resistance to Map, while in water buffalo the same polymorphism did not show any significant association with susceptibility/resistance to TB.

To date, it is estimated that about 50 human diseases are caused entirely or partly by synonymous mutations, and it has also been demonstrated that they play a role in disease resistance [32],[33]. In the human TLR2, for example, synonymous mutations have been identified which are associated with resistance to liver cell carcinoma [17]. A recent meta-analysis identified a significant association between an intron mutation of TLR9 and tuberculosis in humans [34].

Moreover, it has been observed that silent mutations can interfere with various stages of protein “construction”, from DNA transcription to the translation of mRNA into protein.

Silent changes in codons located in the exonic splicing enhancers (ESE) can prevent introns from being properly eliminated. A notable example of the damage that can be caused by mutations in a splicing enhancer has been documented for the CFTR (Cystic Fibrosis Transmembrane Conductance Regulator) gene [35],[36].

Finally, synonymous mutations can alter the stability of the transcript by reducing its affinity for RNA binding proteins. Indeed, in humans, a synonymous mutation of the gene coding for corneodesmosine (CDSN) has been associated with the onset of psoriasis in various ethnic groups [37].

An alternative hypothesis could be that some of the detected synonymous mutations are not directly related to disease susceptibility/resistance, but rather that one or more are linked to other functional mutations on the same chromosome which have not yet been identified. Functional studies will be crucial in order to understand how these polymorphisms act.

Conclusions

The work describes novel mutations in the bubaline TLR2, TLR4 and TLR9 genes and describes their association with M. bovis infection. Studies of the genetic factors involved in complex diseases have not yet provided clear explanations for the onset of such diseases, though they may help to identify their underlying physio-pathological pathways. A thorough knowledge of these mechanisms will improve our ability to determine the risk of developing the disease and to understand it in its entirety. This will pave the way to create multiple layers of disease resistance in herds by selective breeding and also to identify and synthesize innovative drugs.

Methods

Study design

Genetic association between SNPs and TB was studied by a case-control approach on 184 water buffaloes from 22 herds in the Campania region. Sample collection was performed within the Italian National Plan for the control of tuberculosis, which imposes culling of all the subjects positive to the delayed hypersensitivity test [38] and microbiological processing of the lymphoid tissue for isolation of M. bovis [38]. Our analyses, performed during the period 2009-2011, revealed the presence of 59 animals positive to both intradermal TB test and M. bovis microbiological isolation. These 59 animals were therefore classified as cases for the case-control study. The choice of the microbiological test to confirm cases identification was considered in order to exclude the possibility of false-positive subjects diagnosed by the delayed hypersensitivity test. The control subjects were collected among the animals which tested negative to repeated intradermal TB tests (at least three tests) during the entire period of the study. The control animals were at least five years-old, and were born and raised in the same herds of cases to ensure an equal level of exposure to infection. The number of controls was estimated taking into account an odds ratio value of 2.0 as a threshold for significance, with a ratio for cases and controls of 1:2 (power = 80%, confidence = 95%) [12]. Based on these parameters, with an imposed number of cases (59), the required number of controls was 108 for a one-tailed analysis. The present case-control study therefore included 59 cases and 125 controls.

TLRs sequencing

DNA was extracted from blood (control animals) and lymph node (cases) samples using the QIAamp DNA mini kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. PCR primers for amplification of the CDS of the genes TLR2, TLR4 and TLR9 were designed using the Web interface Primer3; they are shown in Table 4. The TLR2 coding sequence was amplified using two primer pairs according to the following touchdown PCR protocol: initial step of 7 cycles of 1 min at 95°C, 1 min at 65°C (decreasing by 1°C after each cycle), 1 min 30 s at 72°C, followed by 35 cycles of 1 min at 95°C, 1 min at 58°C, 1 min 30 s with a final extension of 10 min at 72°C. The TLR4 CDS was amplified using four primers sets with the following thermal profile: 35 cycles of 30 s at 95°C, 30 s at Tann and 45 s at 72°C, with a final extension of 10 min at 72°C. Tann is the annealing temperature specifically indicated for each primers pair in Table 4. The TLR9 CDS was amplified using six primer sets with the following thermal profile: 5 cycles of 30 s at 95°C, 20 s at 60°C, and 15 s at 72°C followed by 30 additional cycles of 20 s at 95°C, 20 s at Tann, and 30 s at 72°C, with a final extension of 10 min at 72°C. All thermal profiles included an initial step at 95°C for 15 min for Taq DNA polymerase activation. PCR reactions were performed with the HotStar HiFidelity Polymerase Kit (Qiagen), and the reaction mixture included 50 ng of genomic DNA and 0.5 μM of each primer (forward and reverse) in a final volume of 25 μl. Amplicons were purified (QiaQuick purification kit, Qiagen), bi-directionally sequenced using the Big Dye Terminator cycle sequencing kit v.1.1 or v.3.1 (Life Technologies) and purified using the DyeEx spin kit (Qiagen). Samples (5μl) were denatured with 10 μl of Hi-Di formamide (Life Technologies) at 95°C for 5 min, and separated by capillary electrophoresis on either ABI PRISM 310 or 3130 sequencers (Life Technologies). Sequencing data were manually inspected by the Sequencing Analysis software v5.4 (Life Technologies).
Table 4

PCR primers, annealing temperatures and amplicon length of amplified bubaline TLRs 2, 4 and 9 genes

Oligo

Sequence (5′ 3′)

PCR T ann(°C)

Amplicon (bp)

TLR2 A

For: TTTGTAGGTCAAATCACTGGACA

58 a

1331 bp

Rev: TCCTGGCCACTGACAAGTTT

TLR2 B

For: GCCCTTCCTTCAAACCTTG

58 a

1234 bp

Rev: CACCACCAGACCAAGACTGA

TLR4 A

For: GTGTGGAGACCTAGATGACTGG

60

706 bp

Rev: GTACGCTATCCGGAATTGTTCA

TLR4 B

For: CTTTCCTGGAGGGACTGTGC

60

435 bp

Rev: CCACGAAGTTTGAACCTAAGGTAA

TLR4 C

For: CTACCAAGCCTTCAGTATCTAG

60

743 bp

Rev: GGCATGTCCTCCATATCTAAAG

TLR4 D

For: AAGGACCAGAGGCAGCTCTT

58

801 bp

Rev: TAACTGAACACGCCCTGCAT

TLR9 A

For: CCAGCCTCTCCTTAATCTCC

54

718 bp

Rev: CGGAACCAATCTTTCTCTAGTT

TLR9 B

For: CCTGACACCTTCAGTCACCT

55

651 bp

Rev: GCGGGTAAACATCTCTTGCT

TLR9 C

For: CGTCAGCTCAAAGGACTTCA

56

546 bp

Rev: AGGGTGTGCAGATGGTTCTC

TLR9 D

For: GGGAGACCTCTATCTCTGCTTT

56

428 bp

Rev: CGCTCACGTCTAGGATTTTC

TLR9 E

For: CTTCAGAAGCTGGACGTGAG

55

685 bp

Rev: TCTTGCGGCTGCTGTAGAC

TLR9 F

For: TGCTCTATGATGCCTTCGTG

55

424 bp

Rev: AGGTTGGCCCAGAAACTACC

a For these primer pairs a touchdown PCR thermal profile was used.

SNP selection and genotyping

Sequences were analyzed by multiple alignment using BioEdit v.7.1.3 [39] and SeqManII 5.00 software (DNASTAR Inc.). The 184 water buffalo sequences of the TLR2, TLR4 and TLR9 genes were compared with the publicly available sequences GenBank:HM756161 (Toll-like receptor 2 gene, complete cds), GenBank:JN786600 (Toll-like receptor 4 gene, complete cds), GenBank:HQ242778 (Toll-like receptor 9 gene, complete cds), respectively. Each forward and reverse sequence from a single DNA sample was compared with the other in order to generate a consensus sequence and to identify polymorphisms among the samples.

Genetic association analysis

Genotype frequencies were tested for deviation from Hardy-Weinberg equilibrium (HWE) by using Fisher Exact test and Chi-square test (χ2). Pearson’s χ2 test was carried out for each SNP both independently and by using a multiple testing correction based on permutation [40],[41].

SNPs that were significantly associated with the disease were entered in an univariate logistic regression model. All statistical analyses were performed using Stata statistical software, release 10.0 (Stata Corp., College Station, TX, USA) and PLINK software (PLINK 1.07, Shaun Purcell http://pngu.mgh.harvard.edu/purcell/plink/).

Haplotype reconstruction was performed using PHASE software, version 2.1 [42].

Declarations

Acknowledgements

We thank Cristiana Maurella, DVM MSE (BEAR – Biostatistica Epidemiologia e Analisi del Rischio, Istituto Zooprofilattico Sperimentale del Piemonte, Liguria e Valle d’Aosta, Turin, Italy) for helpful advices on the statistical analysis.

Funding for the project was provided by the Italian Ministry of Health (grant IZS ME 02/09 RC).

The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Authors’ Affiliations

(1)
Istituto Zooprofilattico Sperimentale del Mezzogiorno
(2)
Istituto Zooprofilattico Sperimentale del Piemonte Liguria e Valle d’Aosta
(3)
Dipartimento di Sanità Pubblica Veterinaria e Sicurezza Alimentare, Istituto Superiore di Sanità
(4)
Dipartimento di Agraria, Università degli Studi di Napoli “Federico II”

References

  1. Coppola S, Parente E, Dumontet S, La Peccerella A: The microflora of natural whey cultures utilized as starters in the manufacture of Mozzarella cheese from water-buffalo milk. Lait. 1988, 68: 295-310. 10.1051/lait:1988319.View ArticleGoogle Scholar
  2. Adams LG, Templeton JW: Genetic resistance to bacterial diseases of animals. Rev Sci Tech. 1998, 17: 200-219.PubMedGoogle Scholar
  3. Bishop SC, MacKenzie KM: Genetic management strategies for controlling infectious diseases in livestock populations. Genet Sel Evol. 2003, 35 (Suppl 1): S3-S17. 10.1186/1297-9686-35-S1-S3.PubMedPubMed CentralView ArticleGoogle Scholar
  4. Capparelli R, Alfano F, Amoroso MG, Borriello G, Fenizia D, Bianco A, Roperto S, Roperto F, Iannelli D: Protective effect of the Nramp1 BB genotype against Brucella abortus in the water buffalo (Bubalus bubalis). Infect Immun. 2007, 75: 988-996. 10.1128/IAI.00948-06.PubMedPubMed CentralView ArticleGoogle Scholar
  5. Trinchieri G, Sher A: Cooperation of Toll-like receptor signals in innate immune defence. Nat Rev Immunol. 2007, 7: 179-190. 10.1038/nri2038.PubMedView ArticleGoogle Scholar
  6. Akira S, Takeda K, Kaisho T: Toll-like receptors: critical proteins linking innate and acquired immunity. Nat Immunol. 2001, 2: 675-680. 10.1038/90609.PubMedView ArticleGoogle Scholar
  7. Takeda K, Akira S: TLR signaling pathways. Semin Immunol. 2004, 16: 3-9. 10.1016/j.smim.2003.10.003.PubMedView ArticleGoogle Scholar
  8. Takeuchi O, Akira S: TOLL-like receptors; their physiological role and signal transduction system. Int Immunopharmacol. 2001, 1: 625-635. 10.1016/S1567-5769(01)00010-8.PubMedView ArticleGoogle Scholar
  9. Koets A, Santema W, Mertens H, Oostenrijk D, Keestra M, Overdijk M, Labouriau R, Franken P, Frijters A, Nielen M, Rutten V: Susceptibility to paratuberculosis infection in cattle is associated with single nucleotide polymorphisms in Toll-like receptor 2 which modulate immune responses against Mycobacterium avium subspecies paratuberculosis . Prev Vet Med. 2010, 93: 305-315. 10.1016/j.prevetmed.2009.11.008.PubMedView ArticleGoogle Scholar
  10. Lorenz E, Mira JP, Cornish KL, Arbour NC, Schwartz DA: A novel polymorphism in the toll-like receptor 2 gene and its potential association with staphylococcal infection. Infect Immun. 2000, 68: 6398-6401. 10.1128/IAI.68.11.6398-6401.2000.PubMedPubMed CentralView ArticleGoogle Scholar
  11. Mucha R, Bhide MR, Chakurkar EB, Novak M, Mikula I: Toll-like receptors TLR1, TLR2 and TLR4 gene mutations and natural resistance to Mycobacterium avium subsp. paratuberculosis infection in cattle. Vet Immunol Immunopathol. 2009, 128: 381-388. 10.1016/j.vetimm.2008.12.007.PubMedView ArticleGoogle Scholar
  12. Pinedo PJ, Buergelt CD, Donovan GA, Melendez P, Morel L, Wu R, Langaee TY, Rae DO: Candidate gene polymorphisms (BoIFNG, TLR4, SLC11A1) as risk factors for paratuberculosis infection in cattle. Prev Vet Med. 2009, 91: 189-196. 10.1016/j.prevetmed.2009.05.020.PubMedView ArticleGoogle Scholar
  13. Smirnova I, Poltorak A, Chan EK, McBride C, Beutler B: Phylogenetic variation and polymorphism at the toll-like receptor 4 locus (TLR4). Genome Biol. 2000, 1: 1-10.1186/gb-2000-1-1-research002.View ArticleGoogle Scholar
  14. Sun L, Song Y, Riaz H, Yang H, Hua G, Guo A, Yang L: Polymorphisms in toll-like receptor 1 and 9 genes and their association with tuberculosis susceptibility in Chinese Holstein cattle. Vet Immunol Immunopathol. 2012, 147: 195-201. 10.1016/j.vetimm.2012.04.016.PubMedView ArticleGoogle Scholar
  15. Tabel Y, Berdeli A, Mir S: Association of TLR2 gene Arg753Gln polymorphism with urinary tract infection in children. Int J Immunogenet. 2007, 34: 399-405. 10.1111/j.1744-313X.2007.00709.x.PubMedView ArticleGoogle Scholar
  16. Wei T, Gong J, Jamitzky F, Heckl WM, Stark RW, Rössle SC: Homology modeling of human Toll-like receptors TLR7, 8, and 9 ligand-binding domains. Protein Sci. 2009, 18: 1684-1691. 10.1002/pro.186.PubMedPubMed CentralView ArticleGoogle Scholar
  17. Medzhitov R: Toll-like receptors and innate immunity. Nat Rev Immunol. 2001, 1: 135-145. 10.1038/35100529.PubMedView ArticleGoogle Scholar
  18. Amaral ME, Grant JR, Riggs PK, Stafuzza NB, Filho EA, Goldammer T, Weikard R, Brunner RM, Kochan KJ, Greco AJ, Jeong J, Cai Z, Lin G, Prasad A, Kumar S, Saradhi GP, Mathew B, Kumar MA, Miziara MN, Mariani P, Caetano AR, Galvão SR, Tantia MS, Vijh RK, Mishra B, Kumar ST, Pelai VA, Santana AM, Fornitano LC, Jones BC, et al: A first generation whole genome RH map of the river buffalo with comparison to domestic cattle. BMC Genomics. 2008, 9: 631-10.1186/1471-2164-9-631.PubMedPubMed CentralView ArticleGoogle Scholar
  19. Mitra M, Taraphder S, Sonawane GS, Verma A: Nucleotide sequencing and SNP detection of Toll-like receptor-4 gene in Murah Buffalo (Bubalus bubalis). Int ScholRes Net. 2012, 7: 659513-Art. IDGoogle Scholar
  20. Poltorak A, He X, Smirnova I, Liu MY, Huffel CV, Du X, Birdwell D, Alejos E, Silva M, Galanos C, Freudenberg M, Ricciardi-Castagnoli P, Layton B, Beutler B: Defective LPS signaling in C3H/HeJ and C57BL/10ScCr mice: mutations in TLR4 gene. Science. 1998, 282: 2085-2088. 10.1126/science.282.5396.2085.PubMedView ArticleGoogle Scholar
  21. Quesniaux V, Fremond C, Jacobs M, Parida S, Nicolle D, Yeremeev V, Bihl F, Erard F, Botha T, Drennan M, Soler MN, Le Bert M, Schnyder B, Ryffel B: Toll-like receptor pathways in the immune responses to mycobacteria. Microbes Infect. 2004, 6: 946-959. 10.1016/j.micinf.2004.04.016.PubMedView ArticleGoogle Scholar
  22. Griebel PJ, Brownlie R, Manuja A, Nichani A, Mookherjee N, Popowych Y, Mutwiri G, Hecker R, Babiuk LA: Bovine toll-like receptor 9: A comparative analysis of molecular structure, function and expression. Vet Immunol Immunopathol. 2005, 108: 11-16. 10.1016/j.vetimm.2005.07.012.PubMedView ArticleGoogle Scholar
  23. Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, Gwadz M, Hurwitz DI, Jackson JD, Ke Z, Lanczycki CJ, Lu F, Marchler GH, Mullokandov M, Omelchenko MV, Robertson CL, Song JS, Thanki N, Yamashita RA, Zhang D, Zhang N, Zheng C, Bryant SH: CDD: specific functional annotation with the Conserved Domain Database. Nucleic Acids Res. 2009, 37: 205-210. 10.1093/nar/gkn845.View ArticleGoogle Scholar
  24. Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, Gwadz M, Hurwitz DI, Jackson JD, Ke Z, Lanczycki CJ, Lu F, Marchler GH, Mullokandov M, Omelchenko MV, Robertson CL, Song JS, Thanki N, Yamashita RA, Zhang D, Zhang N, Zheng C, Bryant SH, CDD: A Conserved domain database for the functional annotation of proteins. Nucleic Acids Res. 2011, 39: 225-229. 10.1093/nar/gkq1189.View ArticleGoogle Scholar
  25. Acutis PL, Martucci F, D’Angelo A, Peletto S, Colussi S, Maurella C, Porcario C, Iulini B, Mazza M, Dell’Atti L, Zuccon F, Corona C, Martinelli N, Casalone C, Caramelli M, Lombardi G: Resistance to classical scrapie in experimentally challenged goats carrying mutation K222 of the prion protein gene. Vet Res. 2012, 43: 8-10.1186/1297-9716-43-8.PubMedPubMed CentralView ArticleGoogle Scholar
  26. Goldmann W, Hunter N, Smith G, Foster J, Hope J: PrP genotype and agent effects in scrapie: change in allelic interaction with different isolates of agent in sheep, a natural host of scrapie. J Gen Virol. 1994, 75: 989-995. 10.1099/0022-1317-75-5-989.PubMedView ArticleGoogle Scholar
  27. Novák K: Functional polymorphisms in Toll-like receptor genes for innate immunity in farm animals. Vet Immunol Immunopathol. 2014, 157: 1-11. 10.1016/j.vetimm.2013.10.016.PubMedView ArticleGoogle Scholar
  28. Werling D, Piercy J, Coffey TJ: Expression of TOLL-like receptors (TLR) by bovine antigen-presenting cells-potential role in pathogen discrimination?. Vet Immunol Immunopathol. 2006, 112: 2-11. 10.1016/j.vetimm.2006.03.007.PubMedView ArticleGoogle Scholar
  29. White SN, Taylor KH, Abbey CA, Gill CA, Womack JE: Haplotype variation in bovine Toll-like receptor 4 and computational prediction of a positively selected ligand-binding domain. Proc Natl Acad Sci U S A. 2003, 100: 10364-10369. 10.1073/pnas.1333957100.PubMedPubMed CentralView ArticleGoogle Scholar
  30. Cargill EJ, Womack JE: Detection of polymorphisms in bovine toll-like receptors 3, 7, 8, and 9. Genomics. 2007, 89: 745-755. 10.1016/j.ygeno.2007.02.008.PubMedView ArticleGoogle Scholar
  31. Torres-García D, Cruz-Lagunas A, García-Sancho Figueroa MC, Fernández-Plata R, Baez-Saldaña R, Mendoza-Milla C, Barquera R, Carrera-Eusebio A, Ramírez-Bravo S, Campos L, Angeles J, Vargas-Alarcón G, Granados J, Gopal R, Khader SA, Yunis EJ, Zuñiga J: Variants in toll-like receptor 9 gene influence susceptibility to tuberculosis in a Mexican population. J Transl Med. 2013, 11: 220-10.1186/1479-5876-11-220.PubMedPubMed CentralView ArticleGoogle Scholar
  32. Chamary JV, Hurst LD: The price of silent mutations. Sci Am. 2009, 300: 46-53. 10.1038/scientificamerican0609-46.PubMedView ArticleGoogle Scholar
  33. Sauna ZE, Kimchi-Sarfaty C, Ambudkar SV, Gottesman MM: The sounds of silence: synonymous mutations affect function. Pharmacogenomics. 2007, 8: 527-532. 10.2217/14622416.8.6.527.PubMedView ArticleGoogle Scholar
  34. Kobayashi K, Yuliwulandari R, Yanai H, Naka I, Lien LT, Hang NTL, Hijikata M, Keicho N, Tokunaga K: Association of TLR polymorphisms with development of tuberculosis in Indonesian females. Tissue Antigens. 2012, 79: 190-197. 10.1111/j.1399-0039.2011.01821.x.PubMedView ArticleGoogle Scholar
  35. Pagani F, Buratti E, Stuani C, Baralle FE: Missense, nonsense, and neutral mutations define juxtaposed regulatory elements of splicing in cystic fibrosis transmembrane regulator exon 9. J Biol Chem. 2003, 278: 26580-26588. 10.1074/jbc.M212813200.PubMedView ArticleGoogle Scholar
  36. Pagani F, Raponi M, Baralle FE: Synonymous mutations in CFTR exon 12 affect splicing and are not neutral in evolution. Proc Natl Acad Sci U S A. 2005, 102: 6368-6372. 10.1073/pnas.0502288102.PubMedPubMed CentralView ArticleGoogle Scholar
  37. Capon F, Trembath RC, Barker JN: An update on the genetics of psoriasis. Dermatol Clin. 2004, 22: 339-347. 10.1016/S0733-8635(03)00125-6.PubMedView ArticleGoogle Scholar
  38. OIE: Bovine tuberculosis. In Manual of Diagnostic Tests and Vaccines for Terrestrial Animals. 2013. Chapter 2.4.7.Google Scholar
  39. Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleis Acids Symposium Series. 1999, 41: 95-98.Google Scholar
  40. Gao X, Starmer J, Martin ER: A multiple testing correction method for genetic association studies using correlated single nucleotide polymorphisms. Genet Epidemiol. 2008, 32: 361-369. 10.1002/gepi.20310.PubMedView ArticleGoogle Scholar
  41. Nyholt DR: A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other. Am J Hum Genet. 2004, 74: 765-769. 10.1086/383251.PubMedPubMed CentralView ArticleGoogle Scholar
  42. Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001, 68: 978-989. 10.1086/319501.PubMedPubMed CentralView ArticleGoogle Scholar

Copyright

© Alfano et al.; licensee BioMed Central. 2014

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.