Microarray analysis of tick-infested skin in resistant and susceptible cattle confirms the role of inflammatory pathways in immune activation and larval rejection

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

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

Application of quantitative proteomics to discover biomarkers for tick resistance in cattle

Introduction

Breeding for tick resistance is a sustainable alternative to control cattle ticks due to widespread resistance to acaricidal drugs and the lack of a protective vaccine. The most accurate method used to characterise the phenotype for tick resistance in field studies is the standard tick count, but this is labour-intensive and can be hazardous to the operator. Efficient genetic selection requires reliable phenotyping or biomarker(s) for accurately identifying tick-resistant cattle. Although breed-specific genes associated with tick resistance have been identified, the mechanisms behind tick resistance have not yet been fully characterised.

Methods

This study applied quantitative proteomics to examine the differential abundance of serum and skin proteins using samples from naïve tick-resistant and -susceptible Brangus cattle at two-time points following tick exposure. The proteins were digested into peptides, followed by identification and quantification using sequential window acquisition of all theoretical fragment ion mass spectrometry.

Results

Resistant naïve cattle had a suite of proteins associated with immune response, blood coagulation and wound healing that were significantly (adjusted P < 10- 5) more abundant compared with susceptible naïve cattle. These proteins included complement factors (C3, C4, C4a), alpha-1-acid glycoprotein (AGP), beta-2-glycoprotein-1, keratins (KRT1 & KRT3) and fibrinogens (alpha & beta). The mass spectrometry findings were validated by identifying differences in the relative abundance of selected serum proteins with ELISA. The proteins showing a significantly different abundance in resistant cattle following early and prolonged tick exposures (compared to resistant naïve) were associated with immune response, blood coagulation, homeostasis, and wound healing. In contrast, susceptible cattle developed some of these responses only after prolonged tick exposure.

Discussion

Resistant cattle were able to transmigrate immune-response related proteins towards the tick bite sites, which may prevent tick feeding. Significantly differentially abundant proteins identified in this research in resistant naïve cattle may provide a rapid and efficient protective response to tick infestation. Physical barrier (skin integrity and wound healing) mechanisms and systemic immune responses were key contributors to resistance. Immune response-related proteins such as C4, C4a, AGP and CGN1 (naïve samples), CD14, GC and AGP (post-infestation) should be further investigated as potential biomarkers for tick resistance.

Serum proteomes of Santa Gertrudis cattle before and after infestation with Rhipicephalus australis ticks

Previous studies have applied genomics and transcriptomics to identify immune and genetic markers as key indicator traits for cattle tick susceptibility/resistance; however, results differed between breeds, and there is lack of information on the use of host proteomics. Serum samples from Santa Gertrudis cattle (naïve and phenotyped over 105 days as tick-resistant [TR] or tick-susceptible [TS]) were used to conduct differential abundance analyses of protein profiles. Serum proteins were digested into peptides followed by identification and quantification using sequential window acquisition of all instances of theoretical fragment ion mass spectrometry. Before tick infestation, abundance of 28 proteins differed significantly (adjusted P < 10^-5 ) between TR and TS. These differences were also observed following tick infestation (TR vs TS) with a further eight differentially abundant proteins in TR cattle, suggesting possible roles in adaptive responses. The intragroup comparisons (TS-0 vs TS and TR-0 vs TR) showed that tick infestation elicited quite similar responses in both groups of cattle, but with relatively stronger responses in TR cattle. Many of the significantly differentially abundant proteins in TR Santa Gertrudis cattle (before and after tick infestation) were associated with immune responses including complement factors, chemotaxis for immune cells and acute-phase responses.

Temporal analysis of the bovine lymph node transcriptome during cattle tick (Rhipicephalus microplus) infestation

Livestock production is a fundamental source of revenue and nutrition, wherein cattle-farming constitutes one of the major agricultural industries. Vectors and vector-borne diseases constitute one of the major factors that decrease the livelihood of all farming communities, more so in resource-poor communities and developing countries. Understanding the immunological responses during tick infestation in cattle is instrumental in the development of novel and improved tick control strategies, such as vaccines. In this study, gene expression patterns were compared within the lymph nodes of three cattle breeds at different life stages of the cattle tick, Rhipicephalus microplus. For Bonsmara (5/8Bos taurus indicus × 3/8B. t. taurus) cattle specifically, some 183 genes were found to be differentially expressed within the lymph nodes during larval and adult tick feeding, relative to uninfested cattle. Overall, the data provides evidence for a transcriptional regulatory network that is activated during immature tick infestation, but is down-regulated towards basal transcriptional levels when adult ticks are feeding. Specific processes in the lymph nodes of Bonsmara cattle were found to be differentially regulated on a transcriptional level. These include: (1) Leukocyte recruitment to the lymph node via chemokines and chemotaxis, (2) Trans-endothelial and intranodal movement on the reticular network, (3) Active regulation of cellular transcription and translation in the lymph node (including leukocyte associated cellular regulatory networks) and (4) Chemokine receptors regulating the movement of cells out of the lymph node. This work provides a first transcriptome analysis of bovine lymph node responses in tick-infested cattle. Findings show a dynamic immune response to tick infestation for the Bonsmara cattle breed, and that suppression of the maturation of the cattle hosts' immunity is especially evident during the larval feeding stages.

Network analysis uncovers putative genes affecting resistance to tick infestation in Braford cattle skin

Background

Genetic resistance in cattle is considered a suitable way to control tick burden and its consequent losses for livestock production. Exploring tick-resistant (R) and tick-susceptible (S) hosts, we investigated the genetic mechanisms underlying the variation of Braford resistance to tick infestation. Skin biopsies from four-times-artificially infested R (n = 20) and S (n = 19) hosts, obtained before the first and 24 h after the fourth tick infestation were submitted to RNA-Sequencing. Differential gene expression, functional enrichment, and network analysis were performed to identify genetic pathways and transcription factors (TFs) affecting host resistance.

Results

Intergroup comparisons of hosts before (Rpre vs. Spre) and after (Rpost vs. Spost) tick infestation found 51 differentially expressed genes (DEGs), of which almost all presented high variation (TopDEGs), and 38 were redundant genes. Gene expression was consistently different between R and S hosts, suggesting the existence of specific anti-tick mechanisms. In the intragroup comparisons, Rpost vs. Rpre and Spost vs. Spre, we found more than two thousand DEGs in response to tick infestation in both resistance groups. Redundant and non-redundant TopDEGs with potential anti-tick functions suggested a role in the development of different levels of resistance within the same breed. Leukocyte chemotaxis was over-represented in both hosts, whereas skin degradation and remodeling were only found in TopDEGs from R hosts. Also, these genes indicated the participation of cytokines, such as IL6 and IL22, and the activation of Wingless (WNT)-signaling pathway. A central gene of this pathway, WNT7A, was consistently modulated when hosts were compared. Moreover, the findings based on a genome-wide association study (GWAS) corroborate the prediction of the WNT-signaling pathway as a candidate mechanism of resistance. The regulation of immune response was the most relevant pathway predicted for S hosts. Members of Ap1 and NF-kB families were the most relevant TFs predicted for R and S, respectively.

Conclusion

This work provides indications of genetic mechanisms presented by Braford cattle with different levels of resistance in response to tick infestation, contributing to the search of candidate genes for tick resistance in bovine.

Towards a new phenotype for tick resistance in beef and dairy cattle: a review

About 80% of the world’s cattle are affected by ticks and tick-borne diseases, both of which cause significant production losses. Cattle host resistance to ticks is the most important factor affecting the economics of tick control, but it is largely neglected in tick-control programs due to technical difficulties and costs associated with identifying individual-animal variation in resistance. The present paper reviews the scientific literature to identify factors affecting resistance of cattle to ticks and the biological mechanisms of host tick resistance, to develop alternative phenotype(s) for tick resistance. If new cost-effective phenotype(s) can be developed and validated, then tick resistance of cattle could be genetically improved using genomic selection, and incorporated into breeding objectives to simultaneously improve cattle productive attributes and tick resistance. The phenotype(s) could also be used to improve tick control by using cattle management. On the basis of the present review, it is recommended that three possible phenotypes (haemolytic analysis; measures of skin hypersensitivity reactions; simplified artificial tick infestations) be further developed to determine their practical feasibility for consistently, cost-effectively and reliably measuring cattle tick resistance in thousands of individual animals in commercial and smallholder farmer herds in tropical and subtropical areas globally. During evaluation of these potential new phenotypes, additional measurements should be included to determine the possibility of developing a volatile-based resistance phenotype, to simultaneously improve cattle resistance to both ticks and biting flies. Because the current measurements of volatile chemistry do not satisfy the requirements of a simple, cost-effective phenotype for use in commercial cattle herds, consideration should also be given to inclusion of potentially simpler measures to enable indirect genetic selection for volatile-based resistance to ticks.

Wonders of tick saliva

Saliva of ticks is arguably the most complex saliva of any animal. This is particularly the case for ixodid species that feed for many days firmly attached to the same skin site of their obliging host. Sequencing and spectrometry technologies combined with bioinformatics are enumerating ingredients in the saliva cocktail. The dynamic and expanding saliva recipe is helping decipher the wonderous activities of tick saliva, revealing how ticks stealthily hide from their hosts while satisfying their gluttony and sharing their individual resources. This review takes a tick perspective on the composition and functions of tick saliva, covering water balance, gasket and holdfast, control of host responses, dynamics, individuality, mate guarding, saliva-assisted transmission, and redundancy. It highlights areas sometimes overlooked - feeding aggregation and sharing of sialomes, and the contribution of salivary gland storage granules - and questions whether the huge diversity of tick saliva molecules is 'redundant' or more a reflection on the enormous adaptability wonderous saliva confers on ticks.

Genome-wide association studies for tick resistance in Bos taurus × Bos indicus crossbred cattle: A deeper look into this intricate mechanism

Rhipicephalus (Boophilus) microplus is the main cattle ectoparasite in tropical areas. Gir × Holstein crossbred cows are well adapted to different production systems in Brazil. In this context, we performed genome-wide association study (GWAS) and post-GWAS analyses for R. microplus resistance in an experimental Gir × Holstein F₂ population. Single nucleotide polymorphisms (SNP) identified in GWAS were used to build gene networks and to investigate the breed of origin for its alleles. Tick artificial infestations were performed during the dry and rainy seasons. Illumina BovineSNP50 BeadChip (Illumina Inc., San Diego, CA) and single-step BLUP procedure was used for GWAS. Post-GWAS analyses were performed by gene ontology terms enrichment and gene transcription factors networks, generated from enriched transcription factors, identified from the promoter sequences of selected gene sets. The genetic origin of marker alleles in the F₂ population was assigned using the breed of origin of alleles approach. Heritability estimates for tick counts were 0.40 ± 0.11 in the rainy season and 0.54 ± 0.11 in the dry season. The top ten 0.5-Mbp windows with the highest percentage of genetic variance explained by SNP markers were found in chromosomes 10 and 23 for both the dry and rainy seasons. Gene network analyses allowed the identification of genes involved with biological processes relevant to immune system functions (TREM1, TREM2, and CD83). Gene-transcription factors network allowed the identification of genes involved with immune functions (MYO5A, TREML1, and PRSS16). In resistant animals, the average proportion of animals showing significant SNPs with paternal and maternal alleles originated from Gir breed was 44.8% whereas the proportion of animals with both paternal and maternal alleles originated from Holstein breed was 11.3%. Susceptible animals showing both paternal and maternal alleles originated from Holstein breed represented 44.6% on average, whereas both paternal and maternal alleles originated from Gir breed animals represented 9.3%. This study allowed us to identify candidate genes for tick resistance in Gir × Holstein crossbreds in both rainy and dry seasons. According to the origin of alleles analysis, we found that most animals classified as resistant showed 2 alleles from Gir breed, while the susceptible ones showed alleles from Holstein. Based on these results, the identified genes may be thoroughly investigated in additional experiments aiming to validate their effects on tick resistance phenotype in cattle.

Bovine Immune Factors Underlying Tick Resistance: Integration and Future Directions

The mechanisms underlying tick resistance within and between cattle breeds have been studied for decades. Several previous papers on bovine immune parameters contributing to tick resistance discussed findings across DNA, RNA, protein, cellular, and tissue levels. However, the differences between bovine host species, tick species and the experimental layouts were not always taken into account. This review aims to (a) give a comprehensive summary of studies investigating immune marker differences between cattle breeds with varying degrees of tick resistance, and (b) to integrate key findings and suggest hypotheses on likely immune-regulated pathways driving resistance. Experimental issues, which may have skewed conclusions, are highlighted. In future, improved experimental strategies will enable more focused studies to identify and integrate immune markers and/or pathways. Most conclusive thus far is the involvement of histamine, granulocytes and their associated pathways in the tick-resistance mechanism. Interestingly, different immune markers might be involved in the mechanisms within a single host breed in contrast to between breeds. Also, differences are evident at each tick life stage, limiting the level to which datasets can be compared. Future studies to further elucidate immune molecule dynamics across the entire tick life cycle and in-depth investigation of promising markers and pathways on both molecular and cellular level are in dire need to obtain a scientifically sound hypothesis on the drivers of tick resistance.

Cattle Tick Rhipicephalus microplus-Host Interface: A Review of Resistant and Susceptible Host Responses

Ticks are able to transmit tick-borne infectious agents to vertebrate hosts which cause major constraints to public and livestock health. The costs associated with mortality, relapse, treatments, and decreased production yields are economically significant. Ticks adapted to a hematophagous existence after the vertebrate hemostatic system evolved into a multi-layered defense system against foreign invasion (pathogens and ectoparasites), blood loss, and immune responses. Subsequently, ticks evolved by developing an ability to suppress the vertebrate host immune system with a devastating impact particularly for exotic and crossbred cattle. Host genetics defines the immune responsiveness against ticks and tick-borne pathogens. To gain an insight into the naturally acquired resistant and susceptible cattle breed against ticks, studies have been conducted comparing the incidence of tick infestation on bovine hosts from divergent genetic backgrounds. It is well-documented that purebred and crossbred Bos taurus indicus cattle are more resistant to ticks and tick-borne pathogens compared to purebred European Bos taurus taurus cattle. Genetic studies identifying Quantitative Trait Loci markers using microsatellites and SNPs have been inconsistent with very low percentages relating phenotypic variation with tick infestation. Several skin gene expression and immunological studies have been undertaken using different breeds, different samples (peripheral blood, skin with tick feeding), infestation protocols and geographic environments. Susceptible breeds were commonly found to be associated with the increased expression of toll like receptors, MHC Class II, calcium binding proteins, and complement factors with an increased presence of neutrophils in the skin following tick feeding. Resistant breeds had higher levels of T cells present in the skin prior to tick infestation and thus seem to respond to ticks more efficiently. The skin of resistant breeds also contained higher numbers of eosinophils, mast cells and basophils with up-regulated proteases, cathepsins, keratins, collagens and extracellular matrix proteins in response to feeding ticks. Here we review immunological and molecular determinants that explore the cattle tick Rhipicephalus microplus-host resistance phenomenon as well as contemplating new insights and future directions to study tick resistance and susceptibility, in order to facilitate interventions for tick control.

Peripheral cellular and humoral responses to infestation with the cattle tick Rhipicephalus microplus in Santa Gertrudis cattle

Resistance to cattle tick infestation in single-host ticks is primarily manifested against the larval stage and results in the immature tick failing to attach successfully and obtain a meal. This study was conducted to identify immune responses that characterize the tick-resistant phenotype in cattle. Thirty-five tick-naïve Santa Gertrudis heifers were used in this study, thirty of which were artificially infested for thirteen weeks with tick larvae while five animals remained at a tick-free quarantine property to serve as a control group. Following thirteen weeks of tick infestation, the animals in this trial exhibited highly divergent tick-resistant phenotypes. Blood samples collected throughout the trial were used to measure peripheral immune parameters: haematology, the percentage of cellular subsets comprising the peripheral blood mononuclear cell (PBMC) population, tick-specific IgG₁ and IgG₂ antibody titres, IgG1 avidity for tick antigens and the ability of PBMC to recognize and proliferate in response to stimulation with tick antigens in vitro. The tick-susceptible cattle developed significantly higher tick-specific IgG₁ antibody titres compared to the tick-resistant animals. These results suggest that the heightened antibody response either does not play a role in resistance or might contribute to increased susceptibility to infestation.

Analyses of reaction norms reveal new chromosome regions associated with tick resistance in cattle

Despite single nucleotide polymorphism (SNP) availability and frequent cost reduction has allowed genome-wide association studies even in complex traits as tick resistance, the use of this information source in SNP by environment interaction context is unknown for many economically important traits in cattle. We aimed at identifying putative genomic regions explaining differences in tick resistance in Hereford and Braford cattle under SNP by environment point of view as well as to identify candidate genes derived from outliers/significant markers. The environment was defined as contemporary group means of tick counts, since they seemed to be the most appropriate entities to describe the environmental gradient in beef cattle. A total of 4363 animals having tick counts (n=10 673) originated from 197 sires and 3966 dams were used. Genotypes were acquired on 3591 of these cattle. From top 1% SNPs (410) having the greatest effects in each environment, 75 were consistently relevant in all environments, which indicated SNP by environment interaction. The outliers/significant SNPs were mapped on chromosomes 1, 2, 5, 6, 7, 9, 11, 13, 14, 15, 16, 18, 21, 23, 24, 26 and 28, and potential candidate genes were detected across environments. The presence of SNP by environment interaction for tick resistance indicates that genetic expression of resistance depends upon tick burden. Markers with major portion of genetic variance explained across environments appeared to be close to genes with different direct or indirect functions related to immune system, inflammatory process and mechanisms of tissue destruction/repair, such as energy metabolism and cell differentiation.

The Essential Role of Tick Salivary Glands and Saliva in Tick Feeding and Pathogen Transmission

As long-term pool feeders, ticks have developed myriad strategies to remain discreetly but solidly attached to their hosts for the duration of their blood meal. The critical biological material that dampens host defenses and facilitates the flow of blood-thus assuring adequate feeding-is tick saliva. Saliva exhibits cytolytic, vasodilator, anticoagulant, anti-inflammatory, and immunosuppressive activity. This essential fluid is secreted by the salivary glands, which also mediate several other biological functions, including secretion of cement and hygroscopic components, as well as the watery component of blood as regards hard ticks. When salivary glands are invaded by tick-borne pathogens, pathogens may be transmitted via saliva, which is injected alternately with blood uptake during the tick bite. Both salivary glands and saliva thus play a key role in transmission of pathogenic microorganisms to vertebrate hosts. During their long co-evolution with ticks and vertebrate hosts, microorganisms have indeed developed various strategies to exploit tick salivary molecules to ensure both acquisition by ticks and transmission, local infection and systemic dissemination within the vertebrate host.

Immune and biochemical responses in skin differ between bovine hosts genetically susceptible and resistant to the cattle tick Rhipicephalus microplus

Background

Ticks attach to and penetrate their hosts’ skin and inactivate multiple components of host responses in order to acquire a blood meal. Infestation loads with the cattle tick, Rhipicephalus microplus, are heritable: some breeds carry high loads of reproductively successful ticks, whereas in others, few ticks feed and reproduce efficiently.

Methods

In order to elucidate the mechanisms that result in the different outcomes of infestations with cattle ticks, we examined global gene expression and inflammation induced by tick bites in skins from one resistant and one susceptible breed of cattle that underwent primary infestations with larvae and nymphs of R. microplus. We also examined the expression profiles of genes encoding secreted tick proteins that mediate parasitism in larvae and nymphs feeding on these breeds.

Results

Functional analyses of differentially expressed genes in the skin suggest that allergic contact-like dermatitis develops with ensuing production of IL-6, CXCL-8 and CCL-2 and is sustained by HMGB1, ISG15 and PKR, leading to expression of pro-inflammatory chemokines and cytokines that recruit granulocytes and T lymphocytes. Importantly, this response is delayed in susceptible hosts. Histopathological analyses of infested skins showed inflammatory reactions surrounding tick cement cones that enable attachment in both breeds, but in genetically tick-resistant bovines they destabilized the cone. The transcription data provided insights into tick-mediated activation of basophils, which have previously been shown to be a key to host resistance in model systems. Skin from tick-susceptible bovines expressed more transcripts encoding enzymes that detoxify tissues. Interestingly, these enzymes also produce volatile odoriferous compounds and, accordingly, skin rubbings from tick-susceptible bovines attracted significantly more tick larvae than rubbings from resistant hosts. Moreover, transcripts encoding secreted modulatory molecules by the tick were significantly more abundant in larval and in nymphal salivary glands from ticks feeding on susceptible bovines.

Conclusions

Compared with tick-susceptible hosts, genes encoding enzymes producing volatile compounds exhibit significantly lower expression in resistant hosts, which may render them less attractive to larvae; resistant hosts expose ticks to an earlier inflammatory response, which in ticks is associated with significantly lower expression of genes encoding salivary proteins that suppress host immunity, inflammation and coagulation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1945-z) contains supplementary material, which is available to authorized users.