Skip to main content
  • Research article
  • Open access
  • Published:

Characterisation and expression profile of the bovine cathelicidin gene repertoire in mammary tissue

Abstract

Background

Cathelicidins comprise a major group of host-defence peptides. Conserved across a wide range of species, they have several functions related to host defence. Only one cathelicidin has been found in humans but several cathelicidin genes occur in the bovine genome. We propose that these molecules may have a protective role against mastitis. The aim of this study was to characterise the cathelicidin gene-cluster in the bovine genome and to identify sites of expression in the bovine mammary gland.

Results

Bioinformatic analysis of the bovine genome (BosTau7) revealed seven protein-coding cathelicidin genes, CATHL 1-7, including two identical copies of CATHL4, as well as three additional putative cathelicidin genes, all clustered on the long arm of chromosome 22. Six of the seven protein-coding genes were expressed in leukocytes extracted from milk of high somatic cell count (SCC) cows. CATHL5 was expressed across several sites in the mammary gland, but did not increase in response to Staphylococcus aureus infection.

Conclusions

Here, we characterise the bovine cathelicidin gene cluster and reconcile inconsistencies in the datasets of previous studies. Constitutive cathelicidin expression in the mammary gland suggests a possible role for these host defence peptides its protection.

Background

Rapid and effective responses to pathogenic challenge are essential for the survival of all living organisms. Production of host-defence peptides (HDPs), important effector molecules of the innate immune response, is key to effective anti-microbial activity at many sites, in many eukaryotic species. Cathelicidins are a family of host-defence peptides found in a diverse range of species, including hagfish [1], amphibians [2], fish [3], birds [4], snakes [5] and mammals ([6] and reviewed in [7]). Initially named antimicrobial peptides (AMPs) and categorised by their ability to act as endogenous antibiotics by disrupting microbial cell membranes [8], it is now becoming clear that their biological activity is multifunctional and includes chemotactic and immunoregulatory activities [9, 10]. For example, the human cathelicidin peptide LL-37 (37 amino acids in length, beginning with two leucines), is chemotactic for neutrophils, monocytes and T cells [11] and has been shown to influence dendritic cell and monocyte function as well as TLR signalling [9].

Cathelicidins are so called because of a highly conserved N-terminal-coding region of the precursor protein known as the cathelin domain, followed by a highly variable domain which codes for peptides with antimicrobial activity (Figure 1A) [12]. The cathelin domain contains two disulfide bonds between cysteine residues C85-C96 and C107-C124 (Figure 1) and was given its name based on sequence similarity to a protein called cathelin, a member of the cystatin superfamily of cysteine proteinase inhibitors (cathe-l-in is an acronym for cathepsin L inhibitor) [13].

Figure 1
figure 1

Cathelicidin structure. A) Cathelicidin genes are approximately 2 kilobases in size with a conserved four exon – three intron arrangement. Exons one, two and three code for the highly conserved N-terminal-coding region of the precursor protein known as the cathelin domain. Exon four codes for a highly variable domain with antimicrobial properties in the mature peptide. B) Alignment of bovine cathelicidin precursor proteins with those from human (Homo sapiens, hCAP-18), mouse (Mus musculus, CAMP), horse (Equus caballus, ecath-2), dog (Canis lupus familiaris, CAMP) and chicken (Gallus gallus, cathl2). (See Additional file 1 for details of the mature peptide encoded by each cathelicidin gene.) Conserved residues are indicated by the asterisk. Exons 1-3 are highly conserved between species with considerable variation in sequence within the exon 4 sequence within and between species. These peptides contain two disulfide bonds between cysteine residues C85-C96 and C107-C124.

Although ubiquitous in mammals, the number of cathelicidin genes in the genome of any single species varies. For instance, while humans and mice have a single cathelicidin gene which codes for an α-helical mature peptide, multiple cathelicidin genes have been found in cattle and sheep [7, 14–16]. In these latter species, the expanded cathelicidin repertoire also includes linear peptides whose mechanism of action may be different to the ancestral α-helical structures [7].

The bovine lineage was one of the first in which cathelicidins were discovered during studies of the antimicrobial activity of bovine neutrophil lysates [17]. In particular, indolicidin, which was first discovered in bovine neutrophils [18] has been shown to exert antimicrobial activity against a number of well-known pathogens, including Staphylococcus aureus[19]. Additional cathelicidin genes have been discovered in the bovine genome using molecular cloning strategies [20–23]. To date, at least seven distinct protein-coding cathelicidin genes have been identified in the bovine genome and are found in a single cluster on chromosome 22 [15]. Southern blot analysis of bovine genomic libraries identified a second copy of CATHL4 and of CATHL1, as well as two related genes that were found to contain premature stop codons [23].

When tested for in vitro antimicrobial activity, cathelicidins rapidly killed a wide range of microorganisms (reviewed in [7, 16]). BMAP-27 and BMAP-28 peptides, encoded by CATHL6 and CATHL5 respectively, display significant levels of antimicrobial activity against E. coli when assayed in milk from cows with mastitis [24] (see Additional file 1 for names of mature antimicrobial peptides encoded by the CATHL genes). The importance of cathelicidins in vivo can be inferred from their localisation at sites which are exposed to microbial invasion as well as in leukocytes. LL-37 has been detected in human neutrophils, mast cells, monocytes and macrophages [25–28]. Epithelial cells in skin, gut, lung, epididymis and mammary gland are other important sources of LL-37 [29–33], which is therefore thought to have a key role in host defence at sites that are in contact with the external environment [34]. CATHL5 transcript has been detected in the mammary gland of healthy cows, but a significant increase in transcript levels in animals with a naturally-occurring intramammary infection was not observed [24].

The sequencing of the bovine genome has provided an opportunity for detailed study of the cathelicidins in a ruminant species. Initial examination of the complete cathelicidin cluster by Elsik et al. [35] indicated there were three additional cathelicidin genes, which they called cathelicidin 8, cathelicidin 9 and cathelicidin 10. Later Dawson et al.[36] searched for immune gene family expansions in artiodactyls and also identified three potential novel bovine cathelicidin genes, which they named CATHL2L, CATHL3L1 and CATHL3L2. However although both groups identified 10 cathelicidins in the bovine genome, the two sets of results were not completely congruent, with differences in the positioning of the novel sequences relative to the known cathelicidins in the cluster. Here, we attempt to reconcile these differences by conducting a comprehensive annotation of the cathelicidin locus in Bos taurus using homology-based search methods of the BLAST family of programs [37] and the more sensitive Hidden Markov Models (HMM) approach [38]. We examined cathelicidin expression in somatic cells isolated from milk samples from animals with high somatic cell counts (SCC), a metric commonly used to detect intramammary infection (mastitis). We also used an in vivo model of mastitis to examine cathelicidin expression across several regions of the mammary gland.

Results

Characterisation of the bovine cathelicidin gene cluster

The BosTau7 assembly of the bovine genome provided us with the opportunity to reconstruct the full cathelicidin region for this species. A HMM profile was constructed based on the alignment of protein sequences corresponding to the seven known bovine cathelicidins. This was then used to search the sequenced bovine genome which had been translated in all six reading frames. Evidence for the presence of all seven protein-coding cathelicidin genes was found on chromosome 22 (BTA22) of the bovine genome BosTau7 assembly (Figure 2). Two copies of CATHL4 were detected within this locus, as were second copies of exons 3 and 4 of CATHL1. This is in agreement with the results obtained by Scocchi et al.[21] using Southern blotting and indicates that our findings are not the result of genome assembly errors. All cathelicidin genes are clustered in an approximately 100 kb length of DNA, located on the long arm of chromosome 22 (22q24). The genes are localized on one strand and all are transcribed in the same direction.

Figure 2
figure 2

Genomic organisation of the bovine cathelicidin gene cluster on bovine chromosome 22 (BTA22). Comparison of recent annotations of the bovine cathelicidin cluster [33, 34]. Our results are based on the BosTau7 (Baylor) genomic assembly, while the Elsik et al. and Dawson et al. results are based on the alternative BosTau6 (University of Maryland) assembly. Both assemblies are based on sequence data from the same Hereford animal. Genes predicted to be pseudogenised contain premature stop codons.

This search led to the identification of 4 cathelin-containing motifs additional to the seven protein expressing genes already annotated: one intact gene, two genes containing premature stop codons and one partial gene duplication. Chromosomal location and strand orientation of the identified cathelicidins was determined using the BLAST-like Alignment Tool (BLAT) at the University of California, Santa Cruz genome browser (Figure 2).

Genomic DNA corresponding to these putative cathelicidins was used for prediction of intron/exon boundaries with GenScan software (Figure 3) before the additional predicted cathelicidin sequences were compared with those proposed by previous authors [35, 36] (Figure 2). ClustalW, in MEGA 5.2 [39], was used to align the additional cathelicidin sequences found in our analysis with the CATHL8, CATHL9 and CATHL10 genes described by Elsik et al.[35], as downloaded from the Bovine Genome Database (http://www.bovinegenome.org[40]). CATHL8 corresponds with the partial copy of the single intact gene we found. However the third and fourth exon of our predicted intact gene share a high level of sequence similarity with other bovine cathelicidins (Figure 1B), yet differ from that predicted by Elsik et al. Another copy of the first two exons of this putative gene lies adjacent. CATHL9 is identical to one of the pseudogenes, with the exception of a tryptophan residue in place of the stop codon we have predicted, while our analysis also predicts a different fourth exon. CATHL10 appears to be identical to the duplication of CATHL4 we have annotated.

Figure 3
figure 3

Predicted structure of CATHL8 . Genomic DNA corresponding to the putative cathelicidin, CATHL8, was retrieved using BLAST-like Alignment Tool (BLAT) at the University of California, Santa Cruz genome browser (http://genome.ucsc.edu) and used for prediction of intron/ exon boundaries using GenScan (http://genes.mit.edu/GENSCAN.html). The coding sequence is in uppercase letters, noncoding sequences are in lowercase letters. The deduced amino acid sequence of the open reading frame is indicated in single letter code and the stop codon is indicated by an x. The TATA-box signal is in bold and underlined and the polyadenylation signal is in bold and double underlined. Line numbers represent genomic co-ordinates on the BosTau7 assembly.

The same procedure was used to align the additional cathelicidin sequences found in our analysis with the Ensembl sequence IDs referenced by Dawson et al.[36] for CATHL2L, CATHL3L1 and CATHL3L2. Based on these alignments, CATHL2L is 100% identical to the partial duplication of CATHL8, while CATHL3L1 is 100% identical to the pseudogenised version of CATHL9 we have predicted. The remaining pseudogene we found shares full sequence similarity with CATHL3L2. Based on these results, we named the predicted additional cathelicidin genes CATHL8, CATHL3L2 and CATHL9, following the rule of priority (Figure 2). See Additional files 2 and 3 for the predicted structures of CATHL9 and CATHL3L2 respectively.

CATHL8 displayed all the characteristics of a functional cathelicidin including 2 kilobase size and a conserved four exon – three intron arrangement (Figure 3) with a TATA-box just upstream from the transcription start site and a polyadenylation signal located 54 bp from the stop codon (Figure 3). In addition, several potential recognition sites for transcription factors involved in the transcription of immune-related genes were found in the 5′ flanking region of this predicted gene (Additional file 4). However, alternative 5′ splice junctions were predicted, changing the 3′ boundary of exon 2 and exon 3 when compared to cathelicidins 1-7 (Figure 4A and B). Therefore, the predicted sequences and their intron-exon boundaries are not consistent with cathelicidins 1-7, and as a result, CATHL8 does not display the canonical cysteine spacing that is invariably conserved among other cathelicidin family members [7].

Figure 4
figure 4

Alternative 5′ splice junctions predicted for exon 2 and exon 3 of CATHL8 and CATHL3L2 . Genomic DNA corresponding to exon 2 and exon 3 of CATHL 1-7, CATHL8, CATHL3L2 and CATHL9 was retrieved from the University of California, Santa Cruz genome browser (http://genome.ucsc.edu). These sequences were aligned using ClustalW, in MEGA5.2 [39]. A) Exon 2 alignment with donor sites (splice junctions) highlighted in grey. B) Exon 3 alignment with donor sites highlighted in grey. Conserved nucleotides are indicated by an asterisk.

Phylogenetic analysis of all protein-coding and predicted bovine cathelicidins revealed that all the predicted putative cathelicidin genes had highest nucleotide identity to CATHL3 (Figure 5). Bootstrap support is relatively low at all nodes because the cathelin-like domain, which is encoded by three out of the four exons, is so highly conserved, meaning the evolutionary signal is mostly coming from the relatively short but variable fourth exon.

Figure 5
figure 5

Phylogram indicating evolutionary relationship between the bovine cathelicidin genes. Maximum likelihood phylogram based on an alignment of the nucleotide sequences of the exons from the 10 bovine cathelicidin genes which have all four exons. Sequences were obtained from the BosTau7 genome assembly using the UCSC genome browser. A Tamura 3-parameter model with a gamma distribution of substitution rates among sites was selected as being the most appropriate using jModelTest2. Numbers at the nodes are the results of 10,000 bootstrap replicates. Scale bar is number of substitutions per site.

Expression of cathelicidin genes 1-8

Cathelicidins 1-7 were constitutively expressed in a range of tissues from multiple body systems (Table 1). However, we were unable to detect expression of the putative genes within cDNA of these tissues.

Table 1 qRT-PCR expression of cathelicidin genes 1-7 in a range of healthy bovine tissues

Cathelicidin gene expression in tissues from healthy mammary gland

We used quantitative real-time PCR (qRT-PCR) to examine the tissue expression patterns of cathelicidins 1-9 in bovine alveolar, ductal, gland cistern and teat canal tissues from the mammary gland of 3 healthy animals and 3 animals which had been experimentally-infected with Staphylococcus aureus. Only CATHL5 mRNA was detected. This was present across all regions in each animal, but did not increase in the infected animals (Figure 6). Once again, none of the samples was positive for CATHL8 or CATHL9 using any of the primer sets we designed.

Figure 6
figure 6

Expression of CATHL5 across several regions in the bovine mammary gland. CATHL5 mRNA is constitutively expressed in alveolar, ductal, gland cistern and teat canal regions of the bovine mammary gland. Both healthy (and Staphylococcus aureus- infected) animals are included; CATHL5 was not significantly upregulated in the infected group. Samples were collected 48 hours after infection. Results are expressed as a ratio of reference gene expression (H3F3A) in each sample, with mean values indicated by horizontal lines.

Cathelicidin gene expression in milk cells from animals with a high somatic cell count

Expression of cathelicidin genes 1-9 was analysed in somatic cells extracted from milk samples taken from 5 cows with a high somatic cell count (Table 2). While six of the seven genes already known to be protein-coding were expressed, expression patterns were highly variable between animals, with no one cathelicidin gene being expressed in all samples. Neither CATHL8 nor CATHL9 expression was detected.

Table 2 Expression of cathelicidin genes 1-7 in 5 Holstein Friesian cows with high somatic cell count (SCC)

Discussion

In this study, we used the most recent assembly of the bovine genome to characterise the cathelicidin gene cluster with a view to reconciling older functional studies with more recent genomic analyses whose results were not fully congruent. It should be noted that these analyses were each carried out on one particular assembly of the genome, of which there have been several, with each being slightly different with respect to the cathelicidin cluster (Figure 2). This search identified seven genes already known to code for functional proteins (cathelicidins 1-7). We also found three additional putative genes whose function is unclear, along with one complete and two incomplete gene duplications clustered on chromosome 22 (Figure 2).

Phylogenetic analysis (Figure 5) indicates CATHL8, CATHL9 and CATHL3L2 are most similar to CATHL3, which gives rise to a linear protein whose proline and arginine-rich motifs are responsible for its antimicrobial function [41]. Gennaro et al. found that cathelicidin peptides of this type need at least 16 residues to be functional. CATHL9 fulfils this criterion, with an antimicrobial domain of 36 residues which contains several proline-arginine repeats [42] (Additional file 2). Despite the sequencing results of this gene indicating a stop codon in exon 1, Scocchi et al. were able to detect a transcript in bovine bone marrow RNA via Northern blotting using an oligonucleotide probe based on the antimicrobial domain (exon 4) of this gene [42]. On further examination, this transcript was found to consist of exons 1, 2 and 3 of CATHL3, which is adjacent to CATHL9, with exons 2, 3 and 4 of CATHL9. A peptide based on the antimicrobial domain of CATHL9 was synthesised and was found to have antimicrobial properties against E. coli at low salt concentrations [42]. Therefore CATHL9 (and perhaps other cathelicidin genes) may extend the bovine cathelicidin repertoire through unusual transcription mechanisms.

Based on the gene expression profiling carried out here, the 7 known protein-coding cathelicidin genes are constitutively expressed in a broad range of tissue types (Table 1), including the mammary gland. Cathelicidin proteins are a major component of the neutrophil secondary granule and increased levels of these proteins are observed in a range of inflammatory conditions [7]. Extensive neutrophil recruitment from the circulation to the lumen of the mammary gland is a hallmark of the early immune response to mammary infection [43]. As an initial step to evaluate the role of the bovine cathelicidins in mammary gland defence, the expression pattern of cathelicidin genes in somatic cells isolated from milk samples was examined. While mRNAs of six of the seven protein-coding cathelicidins were detected, there was no single pattern of expression among the five animals studied (Table 2).

Pryor et al.[44] used a pan-cathelicidin antibody targeting the cathelin-like domain along with antibodies targeting cathelicidin proteins 1, 2 and 4 to carry out Western blotting on milk samples from 35 cows with naturally-occurring mastitis involving a range of pathogens. In addition to finding considerable inter-animal variation in the pattern of individual cathelicidin proteins secretion, several samples tested positive with the consensus antibody, yet tested negative with the three specific antibodies. The authors posited that this may be due to the secretion of additional cathelicidin proteins besides those tested. Our results using leukocytes extracted from milk indicate this is probably the case.

The same group have also induced an intramammary Staphylococcus uberis infection and measured levels of cathelicidin protein in milk using the antibody specific for targeting the cathelicidin consensus sequence, comparing these results to those from animals with naturally occurring infection [45]. While there was a strong correlation between cathelicidin protein concentration and the progression of infection as measured by somatic cell count in the experimentally-infected animals, the results obtained from animals with naturally-occurring mastitis were less conclusive, as a quarter of animals testing positive for bacterial infection had no detectable cathelicidin protein in their milk.

We hypothesised that our in vivo mastitis model may similarly build upon existing data on cathelicidin expression in the mammary gland gathered from naturally infected animals [24]. This model has previously demonstrated that the expression of several innate immune molecules both increases in Staphylococcus aureus infected animals and varies across different regions in the gland [46]. However this was not the case for the cathelicidins. CATHL5 was the only gene whose transcript was detected. While it is constitutively expressed it was not upregulated 48 hours after infection (Figure 6). More extensive analyses of all the cathelicidins across the full spectrum of infection of the mammary gland will be required to determine their effect in this important clinical condition.

It is possible that additional reagents for each HDP may reveal more subtle differences in expression for each cathelicidin. Our demonstration of significant expression of CATHL5 in mammary tissue as well as milk somatic cells suggests that other cells in addition to immune cells are capable of secreting these potent immunomodulatory molecules.

Conclusion

In this study, we used the recently released BosTau7 assembly of the bovine genome to reconstruct the full cathelicidin gene repertoire for this species. In addition to cathelicidins 1-7 which are known to code for distinct peptides, three putative cathelicidin genes cluster on chromosome (BTA22); two of these are pseudogenes and one a potentially functional cathelicidin gene. We also identified two CATHL4 genes adjacent to each other on BTA22. Our in vivo data indicate that CATHL5 is constitutively expressed in the healthy mammary gland, while transcripts of several other protein-coding genes can be detected in leukocytes extracted from the milk of animals with a high somatic cell count.

Methods

Characterisation of the bovine cathelicidin gene cluster

Publicly available protein sequences corresponding to the 7 known bovine cathelicidins were retrieved from Uniprot (http://www.uniprot.org/); (Additional file 1A). To carry out Hidden Markov Model [38] searches of the bovine genome for novel cathelicidin genes, the entire bovine genome was translated in all six reading frames using a purpose-written PERL script. To generate accurate HMM models representing the cathelicidin family, the protein sequences of the seven known bovine cathelicidins were aligned using ClustalW, in MEGA5.2 [39]. These sequences were used in the construction of the HMM by the hmmbuild program in HMMER 2.1.1 (http://www.biology.wustl.edu/gcg/hmmerbuild.html). The generated HMM profile was then searched against the translated genome to identify putative cathelicidin-like regions. Chromosomal location and strand orientation of the identified cathelicidins were determined using the BLAST-like Alignment Tool (BLAT) at the University of California, Santa Cruz (UCSC) genome browser (http://genome.ucsc.edu). Genomic DNA corresponding to putative cathelicidins was retrieved using BLAT and used for prediction of intron/ exon boundaries using GenScan (http://genes.mit.edu/GENSCAN.html; [47]).

Comparison of results with published genomic analyses of the bovine cathelicidin cluster

The three predicted putative cathelicidin sequences located bioinformatically were compared with those published in previous bioinformatic analyses of this gene cluster. Predicted genes described by Elsik et al.[35] were downloaded from the Bovine Genome Project database (http://www.bovinegenome.org; [40]). DNA sequences for the putative CATHL2L, CATHL3L1 and CATHL3L2 genes published by Dawson et al.[36] were downloaded from the UCSC genome browser using the co-ordinates for the referenced ENSEMBL sequence tags. These sequences were aligned with our results using ClustalW, in MEGA 5.2 [39], and sequences were defined as being equivalent when there was 100% agreement between two sequences.

Phylogenetic analysis

The exons from CATHL1-7, as well as those for the putative complete cathelicidins, CATHL8, CATHL9 and CATHL3L2, were downloaded from the BosTau7 genome assembly on the UCSC genome browser. MEGA5.2 [39] was used to produce a maximum likelihood phylogram. jModelTest2 [48] was used to identify the most appropriate model of sequence evolution for the dataset. A Tamura 3-parameter model with a gamma distribution of substitution rates among sites was chosen. 10,000 bootstrap replicates were run to estimate the reliability of each grouping in the final tree.

Bovine tissue panel collection

An extensive range of tissues was collected at a local abattoir from recently euthanized healthy cattle and was immediately flash frozen in liquid nitrogen. Each tissue type was collected from an individual animal. These included tissues from the lung, rumen, small intestine, large intestine, testis, uterus, mammary gland, spleen, liver and lymph node. Peripheral blood mononuclear cells were extracted from whole blood of uninfected cattle using a Percoll™ gradient (GE Healthcare UK, Buckinghamshire, UK) and previously described methods [49]. All experimental procedures were carried out under license from the Irish Department of Health and Children in accordance with the European Community Directive 86-609-EC and were approved by the Animal Research Ethics Committee, University College Dublin.

Bovine mammary gland tissue panel

Our group has previously used an in vivo model to study the effects of S. aureus infection on the innate immune gene expression in the bovine mammary gland [46, see [50] for detailed description of the model]. Tissues from the unchallenged control animals and from the quarter which was infected 48 hours prior to sample collection in treated animals were used in this study. Bacteriological examination of milk samples taken before the experiment indicated these animals were not infected with any major mastitis-causing pathogen prior to the experiment.

Extraction of somatic cells from milk samples

Milk samples were collected from a nearby dairy farm. 100 mL volumes of milk were hand-milked from a single quarter in each of 5 Holstein-Friesian cows identified as having a high somatic cell count at the most recent routine count. All samples were collected midway through milking. Samples were processed following the methodology used by Piepers et al.[51]. 25 mL volumes of milk were diluted 50:50 in cold phosphate-buffered saline (PBS) (Life Technologies, Dublin, Ireland) before centrifuging at 300 × g for 15 minutes at 4°C. Supernatants and fat layers were discarded and cell pellets resuspended in 10 mL of cold PBS. Samples were centrifuged at 4°C for 10 minutes at 200 × g before washing the cells once more. Cells were counted using a Z1 Coulter Counter (Beckman Coulter, Co. Clare, Ireland).

RNA isolation

RNA was extracted from bovine tissue and somatic cells isolated from milk samples using the RNeasy® Mini Kit (Qiagen), and DNAse digested according to the manufacturer’s instructions.

Quantitative real-time PCR (qRT-PCR)

One μg of total RNA from each sample was reverse transcribed into cDNA using OmniScript™ III first strand synthesis kit with oligo (dT) primers according to the manufacturer’s instructions (Life Technologies). The cDNA was quantified using the ND-1000 NanoDrop® spectrophotometer (Thermo Fisher Scientific, Dublin, Ireland) and then diluted to a 20 ng/μl working concentration. Gene specific primers for quantitative real-time PCR were designed using Primer3 express software. These were designed to traverse introns and then commercially synthesized (Life Technologies) (Additional file 1B). qRT-PCR was performed using the Sybr Green-based fluorescent method and the MX3000P® quantitative PCR system (Agilent Technologies, Cork, Ireland), as follows: 95°C for 10 min, followed by 40 cycles of 95°C for 30 sec, 60°C for 1 min and 72°C for 30 sec, and finally amplicon dissociation at 95°C for 1 min, 55°C for 30 sec and 95°C for 30 sec. Each reaction had a total volume of 25μl: 12.5μl of Brilliant II 2X qPCR low ROX master mix (Agilent Technologies), 10.5μl of the primer-H20 and 2μl of cDNA. Optimal primer concentrations were determined by titrating 100, 250 and 800 nM final concentrations and dissociation curves were examined for the presence of a single product. Amplicons were also assessed on 1.5% native agarose gels to ensure they were the predicted length. Only reactions which yielded a single product with a single band of the correct size were used to determine relative expression levels.

Data analysis

For differential expression analysis, qRT-PCR data (Cq values) were converted to gene expression fold changes using the 2-ΔΔCq (Cq represents the quantification cycle) method [52], and recorded relative to control samples. GAPDH was used as a reference gene for the multi-organ panel, while H3F3A was previously found to be the most stably expressed from a panel of reference genes tested in the mammary gland samples [46] as determined using Genorm [53].

For baseline expression (expression in control samples) analysis, expression levels of the gene of interest (GOI) was determined as a ratio of the level of expression of the house keeping gene (HK) using the formula (2-Cq(GOI)/2-Cq(HK)).

Statistical analysis of qRT-PCR results was carried out using the non-parametric Mann–Whitney U test as implemented in version 5.01 of GraphPad Prism (GraphPad Software, San Diego, CA). P values of <0.05 were considered statistically significant.

References

  1. Uzzell T, Stolzenberg ED, Shinnar AE, Zasloff M: Hagfish intestinal antimicrobial peptides are ancient cathelicidins. Peptides. 2003, 24: 1655-1667. 10.1016/j.peptides.2003.08.024.

    Article  CAS  PubMed  Google Scholar 

  2. Hao X, Yang H, Wei L, Yang S, Zhu W, Ma D, Yu H, Lai R: Amphibian cathelicidin fills the evolutionary gap of cathelicidin in vertebrate. Amino Acids. 2012, 43: 677-685. 10.1007/s00726-011-1116-7.

    Article  CAS  PubMed  Google Scholar 

  3. Chang C-I, Pleguezuelos O, Zhang Y-A, Zou J, Secombes CJ: Identification of a novel cathelicidin gene in the rainbow trout, Oncorhynchus mykiss. Infect Immun. 2005, 73: 5053-5064. 10.1128/IAI.73.8.5053-5064.2005.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  4. Lynn DJ, Higgs R, Gaines S, Tierney J, James T, Lloyd AT, Fares MA, Mulcahy G, O’Farrelly C: Bioinformatic discovery and initial characterisation of nine novel antimicrobial peptide genes in the chicken. Immunogenetics. 2004, 56: 170-177. 10.1007/s00251-004-0675-0.

    Article  CAS  PubMed  Google Scholar 

  5. Zhao H, Gan T-X, Liu X-D, Jin Y, Lee W-H, Shen J-H, Zhang Y: Identification and characterization of novel reptile cathelicidins from elapid snakes. Peptides. 2008, 29: 1685-1691. 10.1016/j.peptides.2008.06.008.

    Article  CAS  PubMed  Google Scholar 

  6. Leonard BC, Chu H, Johns JL, Gallo RL, Moore PF, Marks SL, Bevins CL: Expression and activity of a novel cathelicidin from domestic cats. PLoS One. 2011, 6: e18756-10.1371/journal.pone.0018756.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  7. Zanetti M: The role of cathelicidins in the innate host defenses of mammals. Curr Issues Mol Biol. 2005, 7: 179-196.

    CAS  PubMed  Google Scholar 

  8. Brogden KA: Antimicrobial peptides: pore formers or metabolic inhibitors in bacteria?. Nat Rev Micro. 2005, 3: 238-250. 10.1038/nrmicro1098.

    Article  CAS  Google Scholar 

  9. Lai Y, Gallo RL: AMPed up immunity: how antimicrobial peptides have multiple roles in immune defense. Trends Immunol. 2009, 30: 131-141. 10.1016/j.it.2008.12.003.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  10. Wuerth K, Hancock REW: New insights into cathelicidin modulation of adaptive immunity. Eur J Immunol. 2011, 41: 2817-2819. 10.1002/eji.201142055.

    Article  CAS  PubMed  Google Scholar 

  11. De Y, Chen Q, Schmidt AP, Anderson GM, Wang JM, Wooters J, Oppenheim JJ, Chertov O: LL-37, the neutrophil granule- and epithelial cell-derived cathelicidin, utilizes formyl peptide receptor-like 1 (FPRL1) as a receptor to chemoattract human peripheral blood neutrophils, monocytes, and T cells. J Exp Med. 2000, 192: 1069-1074. 10.1084/jem.192.7.1069.

    Article  Google Scholar 

  12. Zanetti M, Gennaro R, Romeo D: Cathelicidins: a novel protein family with a common proregion and a variable C-terminal antimicrobial domain. FEBS Lett. 1995, 374: 1-5. 10.1016/0014-5793(95)01050-O.

    Article  CAS  PubMed  Google Scholar 

  13. Ritonja A, Kopitar M, Jerala R, Turk V: Primary structure of a new cysteine proteinase inhibitor from pig leucocytes. FEBS Lett. 1989, 255: 211-214. 10.1016/0014-5793(89)81093-2.

    Article  CAS  PubMed  Google Scholar 

  14. Zaiou M, Gallo RL: Cathelicidins, essential gene-encoded mammalian antibiotics. J Mol Med. 2002, 80: 549-561. 10.1007/s00109-002-0350-6.

    Article  CAS  PubMed  Google Scholar 

  15. Tomasinsig L, Zanetti M: The cathelicidins- structure, function and evolution. Cur Protein Pept Sci. 2005, 6: 23-34. 10.2174/1389203053027520.

    Article  CAS  Google Scholar 

  16. Zanetti M: Cathelicidins, multifunctional peptides of the innate immunity. J Leukoc Biol. 2004, 75: 39-48.

    Article  PubMed  Google Scholar 

  17. Romeo D, Skerlavaj B, Bolognesi M, Gennaro R: Structure and bactericidal activity of an antibiotic dodecapeptide purified from bovine neutrophils. J Biol Chem. 1988, 263: 9573-9575.

    CAS  PubMed  Google Scholar 

  18. Selsted ME, Novotny MJ, Morris WL, Tang YQ, Smith W, Cullor JS: Indolicidin, a novel bactericidal tridecapeptide amide from neutrophils. J Biol Chem. 1992, 267: 4292-4295.

    CAS  PubMed  Google Scholar 

  19. van Abel RJ, Tang YQ, Rao VS, Dobbs CH, Tran D, Barany G, Selsted ME: Synthesis and characterization of indolicidin, a tryptophan-rich antimicrobial peptide from bovine neutrophils. Int J Pept Protein Res. 1995, 45: 401-409.

    Article  CAS  PubMed  Google Scholar 

  20. Zanetti M, Del Sal G, Storici P, Schneider C, Romeo D: The cDNA of the neutrophil antibiotic Bac5 predicts a pro-sequence homologous to a cysteine proteinase inhibitor that is common to other neutrophil antibiotics. J Biol Chem. 1993, 268: 522-526.

    CAS  PubMed  Google Scholar 

  21. Gennaro R, Skerlavaj B, Romeo D: Purification, composition, and activity of two bactenecins, antibacterial peptides of bovine neutrophils. Infect Immun. 1989, 57: 3142-3146.

    CAS  PubMed Central  PubMed  Google Scholar 

  22. Storici P, Del Sal G, Schneider C, Zanetti M: cDNA sequence analysis of an antibiotic dodecapeptide from neutrophils. FEBS Lett. 1992, 314: 187-190. 10.1016/0014-5793(92)80971-I.

    Article  CAS  PubMed  Google Scholar 

  23. Scocchi M, Wang S, Zanetti M: Structural organization of the bovine cathelicidin gene family and identification of a novel member. FEBS Lett. 1997, 417: 311-315. 10.1016/S0014-5793(97)01310-0.

    Article  CAS  PubMed  Google Scholar 

  24. Tomasinsig L, De Conti G, Skerlavaj B, Piccinini R, Mazzilli M, D’Este F, Tossi A, Zanetti M: Broad-spectrum activity against bacterial mastitis pathogens and activation of mammary epithelial cells support a protective role of neutrophil cathelicidins in bovine mastitis. Infect Immun. 2010, 78: 1781-1788. 10.1128/IAI.01090-09.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  25. Di Nardo A, Vitiello A, Gallo RL: Cutting edge: mast cell antimicrobial activity is mediated by expression of cathelicidin antimicrobial peptide. J Immunol. 2003, 170: 2274-2278.

    Article  CAS  PubMed  Google Scholar 

  26. Agerberth B, Gunne H, Odeberg J, Kogner P, Boman HG, Gudmundsson GH: FALL-39, a putative human peptide antibiotic, is cysteine-free and expressed in bone marrow and testis. Proc Natl Acad Sci. 1995, 92: 195-199. 10.1073/pnas.92.1.195.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  27. Cowland JB, Johnsen AH, Borregaard N: hCAP-18, a cathelin/pro-bactenecin-like protein of human neutrophil specific granules. FEBS Lett. 1995, 368: 173-176. 10.1016/0014-5793(95)00634-L.

    Article  CAS  PubMed  Google Scholar 

  28. Koon HW, Shih DQ, Chen J, Bakirtzi K, Hing TC, Law I, Ho S, Ichikawa R, Zhao D, Xu H, Gallo R, Dempsey P, Cheng G, Targan SR, Pothoulakis C: Cathelicidin Signaling via the Toll-like receptor protects against colitis in mice. Gastroenterology. 2011, 141: 1852.e3-1863.e3.

    Article  Google Scholar 

  29. Murakami M, Lopez-Garcia B, Braff M, Dorschner RA, Gallo RL: Postsecretory processing generates multiple cathelicidins for enhanced topical antimicrobial defense. J Immunol. 2004, 172: 3070-3077.

    Article  CAS  PubMed  Google Scholar 

  30. Malm J, Sørensen O, Persson T, Frohm-Nilsson M, Johansson B, Bjartell A, Lilja H, Ståhle-Bäckdahl M, Borregaard N, Egesten A: The human cationic antimicrobial protein (hCAP-18) is expressed in the epithelium of human epididymis, is present in seminal plasma at high concentrations, and is attached to spermatozoa. Infect Immun. 2000, 68: 4297-4302. 10.1128/IAI.68.7.4297-4302.2000.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  31. Bals R, Wang X, Zasloff M, Wilson J: The peptide antibiotic LL-37/hCAP-18 is expressed in epithelia of the human lung where it has broad antimicrobial activity at the airway surface. Proc Natl Acad Sci USA. 1998, 95: 9541-9546. 10.1073/pnas.95.16.9541.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  32. Hase K, Eckmann L, Leopard JD, Varki N, Kagnoff MF: Cell differentiation is a key determinant of cathelicidin LL-37/human cationic antimicrobial protein 18 expression by human colon epithelium. Infect Immun. 2002, 70: 953-963. 10.1128/IAI.70.2.953-963.2002.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  33. Murakami M, Dorschner R, Stern L, Lin K, Gallo R: Expression and secretion of cathelicidin antimicrobial peptides in murine mammary glands and human milk. Pediatr Res. 2005, 57: 10-15. 10.1203/01.PDR.0000148068.32201.50.

    Article  CAS  PubMed  Google Scholar 

  34. Zasloff M: Defending the epithelium. Nature Med. 2006, 12: 607-608. 10.1038/nm0606-607.

    Article  CAS  PubMed  Google Scholar 

  35. Sequencing TBG, Consortium A, Elsik CG, Tellam RL, Worley KC: The genome sequence of taurine cattle: a window to ruminant biology and evolution. Science. 2009, 324: 522-528.

    Article  Google Scholar 

  36. Dawson HD, Loveland JE, Pascal G, Gilbert JGR, Uenishi H, Mann KM, Sang Y, Zhang J, Carvalho-Silva D, Hunt T, Hardy M, Hu Z, Zhao SH, Anselmo A, Shinkai H, Chen C, Badaoui B, Berman D, Amid C, Kay M, Lloyd D, Snow C, Morozumi T, Cheng RPY, Bystrom M, Kapetanovic R, Schwartz JC, Kataria R, Astley M, Fritz E, Steward C, Thomas M, Wilming L, Toki D, Archibald AL, Bed’Hom B, Beraldi D, Huang TH, Ait-Ali T, Blecha F, Botti S, Freeman TC, Giuffra E, Hume DA, Lunney JK, Murtaugh MP, Reecy JM, Harrow JL, Rogel-Gaillard C, Tuggle CK: Structural and functional annotation of the porcine immunome. BMC Genom. 2013, 14: 332-10.1186/1471-2164-14-332.

    Article  CAS  Google Scholar 

  37. Altschul S, Madden T, Schaffer A, Zhang J, Zhang Z, Miller W, Lipman D: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  38. Eddy SR: Profile hidden Markov models. Bioinformatics. 1998, 14: 755-763. 10.1093/bioinformatics/14.9.755.

    Article  CAS  PubMed  Google Scholar 

  39. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  40. Childers CP, Reese JT, Sundaram JP, Vile DC, Dickens CM, Childs KL, Salih H, Bennett AK, Hagen DE, Adelson DL, Elsik CG: Bovine Genome Database: integrated tools for genome annotation and discovery. Nucleic Acids Res. 2011, 39: D830-D834. 10.1093/nar/gkq1235.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  41. Gennaro R, Scocchi M, Skerlavaj B, Tossi A, Romeo D: Synthetic approach to the identification of the antibacterial domain of bactenecins, Pro/Arg-rich peptides from bovine neutrophils. Peptides: Chemistry, Structure and Biology: Proceedings of the Thirteenth American Peptide Symposium: 20-25 June 1993; Edmonton. Edited by: Hodges RS, Smith JA. 1994, Leiden: ESCOM, 461-463.

    Google Scholar 

  42. Scocchi M, Wang S, Gennaro R, Zanetti M: Cloning and analysis of a transcript derived from two contiguous genes of the cathelicidin family. Biochim Biophys Acta. 1998, 1398: 393-396. 10.1016/S0167-4781(98)00076-1.

    Article  CAS  PubMed  Google Scholar 

  43. Thomas L, Haider W, Hill A, Cook R: Pathologic findings of experimentally induced Streptococcus uberis infection in the mammary gland of cows. Am J Vet Res. 1994, 55: 1723-

    CAS  PubMed  Google Scholar 

  44. Pryor SM, Smolenski GA, Wieliczko RJ, Broadhurst MK, Stelwagen K, Wheeler TT, Haigh BJ: Cathelicidin levels in milk from cows infected with a range of mastitis causing pathogens. Proc New Zeal Soc An. 2010, 70: 243-245.

    Google Scholar 

  45. Smolenski G, Wieliczko R, Pryor S, Broadhurst M, Wheeler T, Haigh B: The abundance of milk cathelicidin proteins during bovine mastitis. Vet Immunol Immunopathol. 2011, 143: 125-130. 10.1016/j.vetimm.2011.06.034.

    Article  CAS  PubMed  Google Scholar 

  46. Whelehan CJ, Meade KG, Eckersall DP, Young FJ, O’Farrelly C: Experimental Staphylococcus aureus infection of the mammary gland induces region-specific changes in innate immune gene expression. Vet Immunol Immunopathol. 2011, 140: 181-189. 10.1016/j.vetimm.2010.11.013.

    Article  CAS  PubMed  Google Scholar 

  47. Burge C, Karlin S: Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997, 268: 78-94. 10.1006/jmbi.1997.0951.

    Article  CAS  PubMed  Google Scholar 

  48. Darriba D, Taboada GL, Doallo R, Posada D: jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012, 9: 772-

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  49. Ulmer A, Scholz W, Ernst M, Brandt E, Flad HD: Isolation and subfractionation of human peripheral blood mononuclear cells (PBMC) by density gradient centrifugation on Percoll. Immunobiology. 1984, 166: 238-250. 10.1016/S0171-2985(84)80042-X.

    Article  CAS  PubMed  Google Scholar 

  50. Eckersall PD, Young FJ, Nolan AM, Knight CH, McComb C, Waterston MM: Acute phase proteins in bovine milk in an experimental model of Staphylococcus aureus subclinical mastitis. J Dairy Sci. 2006, 89: 1488-1501. 10.3168/jds.S0022-0302(06)72216-0.

    Article  CAS  PubMed  Google Scholar 

  51. Piepers S, De Vliegher S, Demeyere K, Lambrecht BN, de Kruif A, Opsomer G: Technical note: Flow cytometric identification of bovine milk neutrophils and simultaneous quantification of their viability. J Dairy Sci. 2009, 92: 626-631. 10.3168/jds.2008-1393.

    Article  CAS  PubMed  Google Scholar 

  52. Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008, 3: 1101-1108. 10.1038/nprot.2008.73.

    Article  CAS  PubMed  Google Scholar 

  53. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: 7-

    Article  Google Scholar 

Download references

Acknowledgements

Thanks to Michael O’Flaherty for facilitating the collection of milk samples at his farm, to Margaret Murray for assistance with milk sample processing. This work was supported by an Irish Department of Agriculture, Fisheries and Food Research Stimulus Fund Grant (RSF 06-340).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Cliona O’Farrelly.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

ABR, COF and CJW wrote the manuscript. CJW, ATL and ABR performed the bioinformatic analysis. CJW, ABR and FN prepared the RNA for qRT-PCR. CJW and ABR performed the qRT-PCR experiments and data analysis. PDE provided mammary tissues from the in vivo mastitis model. FN, AC and KGM assisted with sample collection and manuscript editing. COF proposed and supervised the study. All authors read and approved the final manuscript.

Cormac J Whelehan, Anne Barry-Reidy contributed equally to this work.

Electronic supplementary material

12864_2013_5761_MOESM1_ESM.docx

Additional file 1:A: Cathelicidin Uniprot Accession numbers. B: Gene-specific oligonucleotide primers used for qRT-PCR. (DOCX 18 KB)

12864_2013_5761_MOESM2_ESM.pdf

Additional file 2:Predicted structure of CATHL9. This figure shows the predicted coding sequence of CATHL9, the deduced amino acid sequence and genomic co-ordinates on the Baylor Btau_4.6.1 (BosTau7) assembly. Genomic DNA corresponding to the putative cathelicidin was retrieved using BLAST-like Alignment Tool (BLAT) at the University of California, Santa Cruz genome browser (http://genome.ucsc.edu) and used for prediction of intron/exon boundaries using GenScan (http://genes.mit.edu/GENSCAN.html). The predicted coding sequence is in uppercase letters, noncoding sequences are in lowercase letters. The deduced amino acid sequence of the open reading frame is indicated in single letter code and the stop codon is indicated by an X. Line numbers represent genomic co-ordinates on the BosTau7 assembly. // denotes a break in the sequence. (PDF 147 KB)

12864_2013_5761_MOESM3_ESM.pdf

Additional file 3:Predicted structure of CATHL3L2. This figure shows the predicted coding sequence of CATHL3L2, the deduced amino acid sequence and genomic co-ordinates on the Baylor Btau_4.6.1 (BosTau7) assembly. Genomic DNA corresponding to the putative cathelicidin was retrieved using BLAST-like Alignment Tool (BLAT) at the University of California, Santa Cruz genome browser (http://genome.ucsc.edu) and used for prediction of intron/exon boundaries using GenScan (http://genes.mit.edu/GENSCAN.html). The predicted coding sequence is in uppercase letters, noncoding sequence are in lowercase letters. The deduced amino acid sequence of the open reading frame is indicated in single letter code and the stop codon is indicated by an X. Line numbers represent genomic co-ordinates on the BosTau7 assembly. // denotes a break in the sequence. (PDF 151 KB)

12864_2013_5761_MOESM4_ESM.pdf

Additional file 4:DNA sequence of the CATHL8 promoter region. This figure shows several potential recognition sites for transcription factors involved in the transcription of immune-related genes in the 5′ flanking of CATHL8. Genomic DNA corresponding to the CATHL8 promoter region (866 base pairs upstream of the translation start codon) was retrieved using BLAST-like Alignment Tool (BLAT) at the University of California, Santa Cruz genome browser (http://genome.ucsc.edu). Predicted transcription binding sites were identified using alibaba2 (http://www.gene-regulation.com/pub/programs.html#alibaba2) and TFSEARCH (http://www.cbrc.jp/research/db/TFSEARCH.html) programs. Several potential recognition sites for transcription factors involved in the transcription of immune-related genes are underlined. Translation start codon ATG is in bold and TATA-box is in italics. (PDF 202 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Whelehan, C.J., Barry-Reidy, A., Meade, K.G. et al. Characterisation and expression profile of the bovine cathelicidin gene repertoire in mammary tissue. BMC Genomics 15, 128 (2014). https://doi.org/10.1186/1471-2164-15-128

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-15-128

Keywords