Profiling of differentially expressed genes in sheep T lymphocytes response to an artificial primary Haemonchus contortus infection
© Yang et al.; licensee BioMed Central. 2015
Received: 28 January 2015
Accepted: 6 April 2015
Published: 18 April 2015
Haemonchus contortus is a common bloodsucking nematode causing widespread economic loss in agriculture. Upon H. contortus infection, a series of host responses is elicited, especially those related to T lymphocyte immunity. Existing studies mainly focus on the general immune responses of sheep T lymphocyte to H. contortus, lacking investigations at the molecular level. The objective of this study was to obtain a systematic transcriptional profiling of the T lymphocytes in H. contortus primary-infected sheep.
Nematode-free sheep were orally infected once with H. contortus L3s. T lymphocyte samples were collected from the peripheral blood of 0, 3, 30 and 60 days post infection (dpi) infected sheep. Microarrays were used to compare gene transcription levels between samples. Quantitative RT-PCR was employed to validate the microarray data. Gene Ontology and KEGG pathway analysis were utilized for the annotation of differentially expressed genes.
Our microarray data was consistent with qPCR results. From microarrays, 853, 242 and 42 differentially expressed genes were obtained in the 3d vs. 0d, 30d vs. 0d and 60d vs. 0d comparison groups, respectively. Gene Ontology and KEGG pathway analysis indicated that these genes were involved in metabolism, signaling, cell growth and immune system processes. Functional analysis of significant differentially expressed genes, such as SLC9A3R2, ABCB9, COMMD4, SUGT1, FCER1G, GSK3A, PAK4 and FCER2, revealed a crucial association with cellular homeostasis maintenance and immune response. Our data suggested that maintaining both effective immunological response and natural cellular activity are important for T lymphocytes in fighting against H. contortus infection.
Our results provide a substantial list of candidate genes in sheep T lymphocytes response to H. contortus infection, and contribute novel insights into a general immune response upon infection.
Haemonchus contortus is an important ruminant gastrointestinal nematode, causing diarrhea, anemia, emaciation, and death in serious infection . Haemonchosis results in a significant global economic loss in the farming industry every year. Whilst chemical control is the main strategy , evolution of drug resistance and residues severely reduces the success of anthelmintic treatment programs [3,4]. New control strategies, such as immunological protective method, are urgently in need and are being developed rapidly. At the same time, approaches to better understand the underlying mechanisms of host immune response to H. contortus infection require rapid implementation.
In general, the developmental process of H. contortus in the host involves a few steps . Infective third-stage larvae (L3s) become exsheathed and develop into L4s after oral infection by the host. After that L4s arrive at the abomasums ready to parasitize at 3–5 days post infection (dpi). Before reaching adulthood, they need to escape from the host’s immune system, avoiding host’s rejection. Subsequently the adults begin to lay eggs at approximately 18 dpi and the parasitic load peak by 25–30 dpi. Finally, H. contortus reaches a stable parasitic level upon achieving immune evasion in its host.
To expel the nematodes during such infection, the host relies mainly on T lymphocytes, especially T helper 2 (TH2) cells [6,7]. The TH2-type immune response induces the production of the various cytokines, such as IL-4, IL-5, IL-10, IL-13, IL-25 and IL-31 ; contributes to B cell differentiation leading to the expression of antibodies, such as IgE, IgG1, IgG4 and IgA [9,10]; and gathers eosinophils to target and eradicate the nematodes . In addition, the TH2-type immune response can reduce T helper 1 (TH1)-type immune response-mediated pathological inflammation through cross suppression, a reaction that further challenges nematode survival [12-14]. During such immune responses, T lymphocytes also mobilize factors involved in cell homeostasis to maintain a stable environment for the cells, ensuring the immune system in optimum working conditions.
A number of studies have described the general aspects of sheep T lymphocyte immune responses to H. contortus [15,16], while Nicholas et al. examined gastrointestinal tract and lymph node tissues and gene expression associated with H. contortus resistance in sheep . Thus, we speculated that obvious changes in gene expression profiles of T lymphocytes upon H. contortus infection must exist. We hypothesized that these differentially expressed genes are likely to participate in T lymphocyte immune responses and/or cell homeostasis maintenance. This study therefore aims to obtain the profiles of differentially expressed genes in the sheep T lymphocytes in response to H. contortus infection to increase our current understanding of sheep immune response to H. contortus.
Animals and parasites
Five six-month-old female Hu sheep were raised in nematode-free conditions and handled humanely in strict accordance with the Zhejiang University Guidance for care and use of experimental animals (Zhejiang University, Code ZJU201308-1-10-072). Four sheep were orally infected with 17,000 H. contortus L3s, and one was kept as an uninfected control. Infective larvae of H. contortus ZJ strain were cultured from eggs using standard methods. The eggs were collected from naturally singly-infected sheep.
Peripheral blood samples of each sheep were collected in sodium citrate coated tubes at 0, 3, 30 and 60 dpi. T lymphocytes were then separated using lymphocyte separation medium (Huadong Medicine), purified through nylon wool fiber [18-20] and resuspended in TRIzol reagent (Invitrogen). Samples were preserved in liquid nitrogen until RNA preparation.
Total RNA of T lymphocyte samples were obtained using TRIzol reagent (Invitrogen) following the manufacturer’s instructions and further purified using RNeasy Mini kit (Qiagen). RNA quantity and quality were determined using an Infinigen SSP-3300 ultramicro-spectrophotometer and an Agilent 2100 Bioanalyzer respectively. T lymphocytes RNA samples from three sheep were chosen for microarray experiments based on their good quantity and quality measurement. cRNA of these twelve RNA samples was then synthesized, labeled with Cy3 and used for hybridization.
Sheep draft genome sequences and supplementary sequences from other species such as Bos taurus  were acquired from GenBank to prepare a microarray . An Agilent sheep microarray containing 19,075 UniGene probes, 25,865 expressed sequence tags (ESTs, blast against UniGene, e value > 0.001) and 1,417 control probes were used in this experiment (Biostar Genechip INC., Shanghai). The GEO accession number of the microarray platform is GPL16283 , while that of the microarray series is GSE42302 . The various Cy3 labeled cRNA samples were hybridized to the respective microarrays . Fluorescence signals were scanned by Agilent scanistor with scan resolution 5 μm, PMT 100%.
Data statistics and analysis
Feature Extraction software was used to read the raw data. The raw data files were processed using the Agi4x44Preprocess package based on linear models for microarray data (limma) package developed within the Bioconductor project in the R statistical programming environment. Background correction and normalization, probes filtering and replicated probes merger were included in data pre-processing to generate microarray results for gene expression profiling. In order to eliminate errors caused by one-sample analysis, the microarray results of three infected sheep of the same stage were combined to proceed to multi-sample analysis using limma package (P < 0.05), which were designated in groups as the 0d group, 3d group, 30d group and 60d group. Genes with absolute log2 fold change greater than 1 (P < 0.05) in pair-wise comparison between groups were scored as differentially expressed genes in response to H. contortus infection.
Functional categories of differentially expressed genes were classified with the bioinformatics analysis resource Database for Annotation, Visualization and Integrated Discovery (DAVID, http://david.abcc.ncifcrf.gov/). Gene Ontology (GO), KEGG pathway and PANTHER pathway were used in DAVID for gene function clustering .
Quantitative RT-PCR (qPCR) validation
Primers of ovine beta-actin and five immune-relevant genes selected for quantitative RT-PCR validation
Primer sequence (5′ – 3′)
Product size (bp)
Fc fragment of IgE, high affinity I, receptor for; gamma polypeptide
immunoglobulin lambda-like polypeptide 1
chemokine (C-X-C motif) ligand 2
MHC class I antigen
Analysis of global gene expression
General situation of differentially expressed genes in six comparisons
3d vs. 0d
30d vs. 0d
60d vs. 0d
30d vs. 3d
60d vs. 3d
60d vs. 30d
KEGG Pathway annotation
PANTHER Pathway annotation
Common thirty-five genes differentially expressed in the 3d vs. 0d, 30d vs. 0d and 60d vs. 0d comparisons
Log fold change
3d vs. 0d
3d vs. 0d
3d vs. 0d
similar to A kinase (PRKA) anchor protein (yotiao) 9
solute carrier family 9 (sodium/hydrogen exchanger), member 3 regulator 2
sodium channel, voltage-gated, type VIII, alpha
similar to Chromosome 1 open reading frame 2
retinoid X receptor, beta
nudix (nucleoside diphosphate linked moiety X)-type motif 14
UDP-Gal: betaGlcNAc beta 1,4- galactosyltransferase, polypeptide 2
ATP-binding cassette, sub-family B (MDR/TAP), member 9
kallikrein related-peptidase 10
chromosome 5 open reading frame 24 ortholog
similar to prion-like protein doppel
mast cell proteinase-1
SMAD family member 4
zinc finger protein 330
COMM domain containing 4
GO enrichment analysis
KEGG pathway analysis
The correlative genes involved in the immune-relevant pathways from the results of KEGG pathway analysis in the six comparisons
3d vs. 0d
30d vs. 0d
60d vs. 0d
30d vs. 3d
60d vs. 3d
Hematopoietic cell lineage
Complement and coagulation cascades
C5AR1↑, C1R↑, CD59↑
Toll-like receptor signaling pathway
IRF7↓, MAP3K7IP2↓, MAPK9↓
IRAK1↑, MAP3K7IP2↑, MAPK9↑, IRF7↑
PIK3R2↑, MAPK9↑, IRF7↑
NOD-like receptor signaling pathway
SUGT1↓, MAP3K7IP2↓, MAPK9↓
MAP3K7IP2↑, MAPK9↑, SUGT1↑
RIG-I-like receptor signaling pathway
IRF7↓, LOC782671↓, MAPK9↓
LOC782671↑, MAPK9↑, IRF7↑
LOC782671↑, MAPK9↑, IRF7↑
Cytosolic DNA-sensing pathway
Natural killer cell mediated cytotoxicity
PPP3R1↑, SHC1↑, FCER1G↑
PIK3R2↑, NFATC1↑, PPP3R1↑, FCER1G↑
T cell receptor signaling pathway
PPP3R1↓, NCK2↑, PAK4↑
PAK4↓, NCK2↓, PIK3R2↑, NFATC1↑, PPP3R1↑
B cell receptor signaling pathway
PIK3R2↑, NFATC1↑, PPP3R1↑
Fc epsilon RI signaling pathway
PIK3R2↑, MAPK9↑, FCER1G↑
Leukocyte transendothelial migration
CYBA↓, CYBB↑, MSN↓, CTNNB1↓, MYL9↑
CXCR4↓, MYL9↓, RHOH↓, CYBB↑
MYL9↓, RHOH↓, PIK3R2↑, CYBB↓, CYBA↑
Intestinal immune network for IgA production
TGFB1↓, PIGR↑, TGFB2↑
CXCR4↓, TGFB2↓, TGFB1↑
Cemokine signaling pathway
FOXO3↓, GNB2↓, GSK3A↓
CXCR4↓, GNB2↑, FOXO3↑, GSK3A↑, SHC1↑
GNG10↓, PIK3R2↑, GSK3A↑, FOXO3↑, GNB2↑
Fc gamma R-mediated phagocytosis
Top up- or down-regulated genes
The top 10 up- or down-regulated genes in each of the six comparisons are listed in Additional file 3. Approximately half of these genes were uncharacterized, while the other half consisted of genes of known or noted hypothetical performances that include roles in immune response, protein metabolic process, sugar metabolic process and signaling transduction. Of particular interest are SUGT1, FCER1G, GSK3A, PAK4, SHC1 and FCER2, which are immune-related genes.
Reverse transcription and quantitative RT-PCR (qPCR) validation
In this experiment, we infected sheep with a single bolus dose of H. contortus larvae to simulate a primary infection. While a trickle infection challenge would be a better strategy to mirror field-acquired immunity, mimicking a mature individual developing immune response at a natural speed, our approach and data, on the other hand reflect a much more direct response to infection in sheep, including initial innate immune responses, followed by adaptive immune response to surviving worms.
As the sheep genome has yet to be sequenced, we used microarray composed of ovine sequences and useful sequences from other species, to identify genes that were significantly differentially expressed in the various stages in sheep T lymphocytes upon H. contortus infection. Our results show that the number of modulated genes decreased with post infection day, inferring that H. contortus had eventually achieved homeostasis with the host via immune evasion. This inference is supported by the fecal egg count (FEC) during H. contortus infection and the final abomasal necropsy, where FEC peaked at 2,400 eggs per gram (data not shown) while only a few adult nematodes were found at the point of abomasal necropsy, even though each sheep was infected with as many as 17,000 L3s. This drop was a significant reduction.
The bulk of modulated genes in the 3d vs. 0d, 60d vs. 0d, and 60d vs. 30d comparisons were down-regulated genes. However majority of those in the 30d vs. 0d, 30d vs. 3d and 60d vs. 3d comparisons were up-regulated (Table 2). When sheep are being challenged by a large number of L3s, the larvae migrate, grow up to blood-feeding L4s and begin to invade the host system within a very short time. During this early stage of infection, larvae infiltrated the abomasal crypts and released a substantial arsenal of secreted products . These products with known roles in tissues degradation and possible immune-modulatory roles are likely to circulate through systematic circulation, destroying immune cells including T lymphocytes, in an effort to assist larvae development and evasion . In addition, the dose of 17,000 worms is considered to be an enormous amount for nematode-naive sheep and represents a dose that can lead to host fatality. As such, it is possible that the initial transcript patterns reflect the physical and physiological trauma that the animal was experiencing. Thus, down-regulated genes accounted for most of the modulated genes in the 3d vs. 0d comparison. At 25–30 dpi, the H. contortus larvae that failed to be expelled from the host system would now present as adults. This status depends delicately on the interactions between the parasites and their host . Female adults were fully gravid by now after copulation with the males. While the adults continued to grow, the sheep began to mount a stronger T lymphocyte immune response against these parasites and the millions of egg produced . Thus, differentially expressed genes related to immune system were mostly up-regulated genes in the 30d vs. 0d and 30d vs. 3d comparisons. However, more down, than up-regulated genes were observed at 60 dpi. This observation could be due to a strong wave of immune response elicited in the middle of the infection that has brought about a decreased worm burden during the later stage of infection. Only about one hundred adults were (data not shown) found in the abomasums when the sheep were autopsied at 60 dpi. These remaining parasites managed to survive from the host immune responses by immune-escape and eventually achieving homeostasis within its host.
To investigate the modulated genes in greater detail, we applied GO annotation for gene clustering. About half of the genes are clustered into one or more than one GO terms. The most frequent GO terms are binding, cellular process, biological regulation and metabolic process. To function optimally, T lymphocytes perform certain housekeeping tasks such as protecting cell structure from parasite attack and maintaining intracellular biological regulation . In association with these protective measures against parasite invasion, many genes involved in cellular homeostasis maintenance were differentially expressed throughout the infection (Additional file 1). Several genes help in the communication between T lymphocytes and other immune cells and/or antigens while others assist in cell metabolism. This clustering data revealed a strong relationship between sheep immune response and cellular homeostasis.
Interestingly, there are only fourteen known genes among the thirty-five genes shared in the 3d vs. 0d, 30d vs. 0d and 60d vs. 0d comparisons, including SLC9A3R2, SCN8A, B4GALT2, ZNF330, ABCB9, MCP1, RXRB, C7H5ORF24, COMMD4, NUDT14, KLK10, SLC28A1, SMAD4, and SEPT7. SLC9A3R2 encodes the sodium-hydrogen exchange regulator confactor NHE-RF2, which generally interacts with a sodium/hydrogen exchanger NHE3 that is primarily responsible for transepithelial sodium balance to maintain endothelial homeostasis . Decreased transcription of SLC9A3R2 at 3 dpi may reflect an initial physical trauma experienced from an enormous infection dose, which gave rise to an imbalance intracellular sodium level. A subsequent increased transcriptional activity at 30 dpi and 60 dpi helped to recover the disturbed salt balance and provided signals to stimulate chemokine generation and NF-κB activation for the regulation of a series of immunological responses . SCN8A encodes for a protein that modulates sodium balance as well. It is the alpha subunit of a voltage-gated sodium channel, forming ion conduction pore to produce a functional channel. The beta and alpha subunits work together to regulate channel gating [32,33]. Like SLC9A3R2, SCN8A down and up-regulation at 3 dpi and the two later time points respectively also indicated sodium balance regulation during infection. ZNF330, a highly conserved protein, was first characterized immunologically as a human autoantigen . It has been demonstrated that this zinc finger protein is capable of mediating mitochondrial apoptotic pathway and directing specific proteins to the nucleolus to regulate nucleolar transcription . ZNF330 displayed a similar expression profile and function as B4GALT2. B4GALT2 is a member of the beta-1,4-galactosyltransferase gene family. As membrane-bound glycoproteins, members in the family direct proteins to the Golgi apparatus and remain as a transmembrane anchor . It is undeniable that H. contortus infection induces cellular damage in sheep. Reduced transcription of B4GALT2 and ZNF330 at 3 dpi was consistent with the inability of protein transportation between or within T lymphocytes, later where cells tried to activate apoptosis to achieve homeostasis by increasing B4GALT2 and ZNF330 levels at 30 and 60 dpi. . ABCB9 is a member of the MDR/TAP subfamily of ATP-binding cassette (ABC) transporters. Though the function of the transporter has yet to be determined, ABCB1, another member of the MDR/TAP subfamily, shows correlation with migration of immune cells and adaptive immunity [38,39]. Fujimoto speculated that more ABCB9 might help to regulate antigen presentation and multidrug resistance . Increased expression of MCP1 encoded mast cell protease 1 would contribute to the clearance of gastrointestinal nematode . Mast cell proteases are protective mediators of inflammation . Up-regulated levels of ABCB9 and MCP1 during infection help the host to mobilize immune cells and present antigen against H. contortus infection. A higher transcription level of RXRB encoding retinoid X receptor beta, a nuclear receptor, modulates cell growth and differentiation through increasing transcriptional function and DNA binding . The gene maps to the MHC class II regions . C7H5ORF24, a bovine ortholog of human chromosome 5 open reading frame, acts as a cancer marker to address DNA damage response [45,46]. The up-regulation of RXRB and C7H5ORF24 at the three time points indicated modulation of cell development and repair during T lymphocytes homeostasis maintenance. COMMD4, down-regulated in the three comparisons, is a protein widely conserved throughout evolution and involved in NF-κB termination [47-49]. NF-κB plays a key role in cellular responses to stimuli, as a rapid-acting primary transcription factor. Κ light chains of NF-κB are also key components of immunoglobulins, making NF-κB an important humoral immune modulator. The factor was activated to up-regulate genes associated with T-cell development and humoral immune response when the level of COMMD4 is low. Misregulation of NF-κB could possibly lead to cancer, viral infection and inadequate immune development. In brief, the protein complex finally manifests as a cytokine production and cell survival factor . Decreased expression of NUDT14 encoding for UDP-glucose pyrophosphatases (UGPPase) would influence the synthesis of glycogen from glucose (UDP-glucose metabolism) . Methylation of KLK10 would be a novel prognostic marker of cancers . A lower level of KLK10 limited the activity of tumorigenesis. CNT1, expressed by SLC28A1, is down-regulated that would affect nucleotide biosynthesis in nucleoside salvage pathways and pyrimidine-nucleoside transportation . SMAD4, down-regulated during infection, binds receptor-regulated SMADs to serve as a transcription activator that regulates TGF-beta receptor-mediated signaling . SEPT7 encoded protein is highly similar to yeast CDC10, which is localized to the cytoplasmic membrane and is an absolute requirement for the completion of cytokinesis. However, decreased SEPT7 expression has no effect on mitosis in T lymphocytes, which suggest that this septin may not be the key molecule even though it is able to destabilize microtubule . Lower transcription levels of NUDT14, KLK10, SLC28A1, SMAD4 and SEPT7 during a series of stages of infection may reflect the inabilities of these genes to participate in glucose metabolism, tumorigenesis, nucleoside biosynthesis and transport, signaling and microtubule stabilization respectively. However other supplementary genes, such as C7H5ORF24 mentioned above, were able to perform similar functions to maintain T lymphocyte homeostasis. Down-regulation of SLC9A3R2, SCN8A, B4GALT2, ZNF330, NUDT14, KLK10, SLC28A1, SMAD4 and SEPT7 in the 3d vs. 0d comparison suggested a strong homeostasis regulation during the early stage of infection. This is likely to be contributed by the initial trauma induced by an enormous infection dose used in our experiments. Twenty-one genes remained unknown. Our data implies that maintaining both an effective immunological response and a constant cellular environment within the host is of equal importance during T lymphocytes battle against H. contortus infection.
Differentially expressed genes were also subjected to KEGG pathway analysis to obtain a more systematic functional annotation. In addition to signaling, cell communication and cell growth and death, KEGG annotations locate more genes to immune system. These genes contribute to almost all immunological processes. Of particular interest are SUGT1, PAK4, FCER1G, GSK3A, SHC1 and FCER2. These genes are also amongst the top ten differentially expressed genes. PAK4 serves as an enzyme in the family of serine/threonine p21-activating kinases. A higher expression of PAK4 at 3 dpi and 30dpi may suggest T cell activation and inflammation mediated by chemokine and cytokine signaling pathway [56,57]. PAK proteins play a role in cytoskeleton reorganization and nuclear signaling and have been implicated in a large-scale of biological activities. PAK4 is also a critical effector that link Rho GTPases to actin cytoskeleton reorganization. Furthermore, this kinase may act as a factor to support cytoskeleton organization of T lymphocytes. The up-regulation of SUGT1 enhanced NOD-like receptor activation . The NOD-like receptors are involved in inflammation regulation through the activation of IL-1 and the NF-κB signaling pathway, meanwhile the proteins also regulate apoptotic responses. Depletion of SUGT1 prevents multiple cellular responses associated with Nod1 activation, but did not affect Nod1 protein stability. Heat shock proteins are binding partners of the SUGT1 protein, involved in innate immune responses . Stronger inflammatory responses are suggested in T lymphocytes at 30 dpi by greater SUGT1 and PAK4 expression compared to that at 0 dpi. GSK3A is a multi-task serine kinase with a role in phosphoinositide 3-kinase (especially PIK3CG, which regulates cytotoxicity in NK cells) signaling pathway and chemokine signaling pathway [60,61]. The up-regulated glycogen synthase kinase modulates various transcription factors, as well as embryo development. However the kinase seems to be non-essential to vitality, as GSK3A mutant mice showed no significant abnormalities. There are three SHC1 protein isoforms. These isoforms differ in activity and domain structures, showing vital functions in various aspects. A higher level of SHC1 delivers a strong signal for chemokine and natural killer cell mediated cytotoxicity signaling pathway [62-64]. When SHC1 acts as a dynamic scaffold protein in cell surface receptor, they may regulate cell invasion and cytoskeleton reorganization. The proteins also operate in the regulation of apoptosis and tumorigenicity in mammalian cells. GSK3A and SHC1 levels were increased at 30 dpi and 60 dpi. They co-regulate chemokine signaling pathways while mediate other immune-relevant signaling pathways separately. The four genes mentioned above exhibit roles in both immunological and cellular processes in T lymphocytes. FCER1G is the gamma chain of high-affinity receptor for the Fc region of IgE. Up-regulated FCER1G directs mast cell and coupling allergens to suppress inflammatory and immediate hypersensitivity responses in immunity to parasites and in allergic diseases [65-67]. Another receptor for IgE, FCER2 (CD23), also known as the low-affinity receptor, is involved in resistance to parasites and is also found to be up-regulated in this study. CD23 molecules regulate IgE levels, and then bind to the IgE immune complexes on B cells to activate macrophages and eosinophils . Upon H. contortus infection, secretions from nematodes enter the blood stream of the sheep and are captured by IgE. The antigen is then transferred from CD23+ B cells to CD4+ T cells via CD11c + antigen presenting cells, leading to an intensive antibody response . FCER1G cooperates with FCER2 to enhance host humoral immune responses, which are reflected by their up-regulations at 30 dpi and 60 dpi.
Analysis of the top ten differentially expressed genes revealed a big increase in the expression of IL-13 in the 60d vs. 30d comparison. Most cytokine kinetics studies have reported early up-regulation in IL-13 mRNA levels by 3dpi, peak at 7dpi followed by a gradual decline until 30dpi. IL-13, one of the key cytokines in Th2-type immune response, induces IgE secretion that in turn, guiding immune cells to kill nematodes . Yet a striking discovery in our experiment shows a recovery in IL-13 mRNA levels by 60 dpi. We speculate that this interesting observation may imply a final “mop-up” phenomenon to eliminate all nematodes at the final stage of sheep immunity against H. contortus infection.
The data reported in our manuscript was collected from animals. The ethics was approved by the experimental animal ethics committee of Zhejiang University and the reference number is ZJU201308-1-10-072.
In summary, the present study sought to describe changes in gene expression across various time points during the infection period of sheep T lymphocytes by H. contortus using microarrays and modern bioinformatics analysis. A list of differentially expressed genes having roles in homeostasis maintenance and immune response was obtained, together with many other genes of unknown function. Changes occurring in the sheep T lymphocytes at the molecular level correspond well to the general description of nematode immunity at the cellular level. The elaborate list of candidate genes obtained here provides new research directions for the study of long-term resistance to H. contortus infection in sheep. Future work includes the identification of novel H. contortus infection resistance genes in sheep.
We thank Zhejiang Academy of Agricultural Sciences for their assistance in collecting H. contortus eggs. We thank Biostar Genechip INC., Shanghai for their assistance in microarray experiments, original data processing and corresponding bioinformatics analysis. The project was supported by grant from the National Natural Science Foundation of China (No. 30972174, 31372426), and the National Key Basic Research Program (973 Program) of China (No. 2015CB150300).
- Zhou QJ, Zhang HL, Jiang XL, Du AF. The gene structure and promoter region of the vaccine target aminopeptidase H11 from the blood-sucking nematode parasite of ruminants, Haemonchus contortus. Funct Integr Genomics. 2010;10(4):589–601.View ArticlePubMedGoogle Scholar
- Freeman AS, Nghiem C, Li J, Ashton FT, Guerrero J, Shoop WL, et al. Amphidial structure of ivermectin-resistant and susceptible laboratory and field strains of Haemonchus contortus. Vet Parasitol. 2003;110(3–4):217–26.View ArticlePubMedGoogle Scholar
- Sutherland IA, Leathwick DM. Anthelmintic resistance in nematode parasites of cattle: a global issue? Trends Parasitol. 2011;27(4):176–81.View ArticlePubMedGoogle Scholar
- James CE, Hudson AL, Davey MW. Drug resistance mechanisms in helminths: is it survival of the fittest? Trends Parasitol. 2009;25(7):328–35.View ArticlePubMedGoogle Scholar
- Nikolaou S, Gasser RB. Prospects for exploring molecular developmental processes in Haemonchus contortus. Int J Parasitol. 2006;36(8):859–68.View ArticlePubMedGoogle Scholar
- Lacroux C, Nguyen TH, Andreoletti O, Prevot F, Grisez C, Bergeaud JP, et al. Haemonchus contortus (Nematoda: Trichostrongylidae) infection in lambs elicits an unequivocal Th2 immune response. Vet Res. 2006;37(4):607–22.View ArticlePubMedGoogle Scholar
- Shakya KP, Miller JE, Horohov DW. A Th2 type of immune response is associated with increased resistance to Haemonchus contortus in naturally infected Gulf Coast Native lambs. Vet Parasitol. 2009;163(1–2):57–66.View ArticlePubMedGoogle Scholar
- Meeusen EN, Balic A, Bowles V. Cells, cytokines and other molecules associated with rejection of gastrointestinal nematode parasites. Vet Immunol Immunopathol. 2005;108(1–2):121–5.View ArticlePubMedGoogle Scholar
- Bambou JC, De la Chevrotiere C, Varo H, Arquet R, Kooyman FN, Mandonnet N. Serum antibody responses in Creole kids experimentally infected with Haemonchus contortus. Vet Parasitol. 2008;158(4):311–8.View ArticlePubMedGoogle Scholar
- De la Chevrotiere C, Bambou JC, Arquet R, Jacquiet P, Mandonnet N. Genetic analysis of the potential role of IgA and IgE responses against Haemonchus contortus in parasite resistance of Creole goats. Vet Parasitol. 2012;186(3–4):337–43.View ArticlePubMedGoogle Scholar
- Terefe G, Lacroux C, Prevot F, Grisez C, Bergeaud JP, Bleuart C, et al. Eosinophils in Haemonchus contortus-infected resistant and susceptible breeds of sheep: abomasal tissue recruitment and in vitro functional state. Vet Parasitol. 2009;165(1–2):161–4.View ArticlePubMedGoogle Scholar
- Fox JG, Beck P, Dangler CA, Whary MT, Wang TC, Shi HN, et al. Concurrent enteric helminth infection modulates inflammation and gastric immune responses and reduces helicobacter-induced gastric atrophy. Nat Med. 2000;6(5):536–42.View ArticlePubMedGoogle Scholar
- Zaccone P, Fehervari Z, Jones FM, Sidobre S, Kronenberg M, Dunne DW, et al. Schistosoma mansoni antigens modulate the activity of the innate immune response and prevent onset of type 1 diabetes. Eur J Immunol. 2003;33(5):1439–49.View ArticlePubMedGoogle Scholar
- Chen F, Liu Z, Wu W, Rozo C, Bowdridge S, Millman A, et al. An essential role for TH2-type responses in limiting acute tissue damage during experimental helminth infection. Nat Med. 2012;18(2):260–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Miller JE, Horohov DW. Immunological aspects of nematode parasite control in sheep. J Anim Sci. 2006;84(Suppl):E124–32.PubMedGoogle Scholar
- Shakya KP, Miller JE, Lomax LG, Burnett DD. Evaluation of immune response to artificial infections of Haemonchus contortus in Gulf Coast Native compared with Suffolk lambs. Vet Parasitol. 2011;181(2–4):239–47.View ArticlePubMedGoogle Scholar
- Andronicos N, Hunt P, Windon R. Expression of genes in gastrointestinal and lymphatic tissues during parasite infection in sheep genetically resistant or susceptible to Trichostrongylus colubriformis and Haemonchus contortus. Int J Parasitol. 2010;40(4):417–29.View ArticlePubMedGoogle Scholar
- Takeuchi A, Shiota M, Tatsugami K, Yokomizo A, Kuroiwa K, Dejima T, et al. YB-1 suppression induces STAT3 proteolysis and sensitizes renal cancer to interferon-alpha. Cancer Immunol Immunother. 2013;62(3):517–27.View ArticlePubMedGoogle Scholar
- Kaburaki T, Zhang Q, Jin X, Uchiyama M, Fujino Y, Nakahara H, et al. Effects of Japanese herbal medicine Sairei-to on murine experimental autoimmune uveitis. Graefes Arch Clin Exp Ophthalmol. 2013;251(12):2733–9.View ArticlePubMedGoogle Scholar
- Trizio D, Cudkowicz G. Separation of T and B lymphocytes by nylon wool columns: evaluation of efficacy by functional assays in vivo. J Immunol. 1974;113(4):1093–7.PubMedGoogle Scholar
- Bar-Or C, Czosnek H, Koltai H. Cross-species microarray hybridizations: a developing tool for studying species diversity. Trends Genet. 2007;23(4):200–7.View ArticlePubMedGoogle Scholar
- Du AF: Gene Expression Omnlbus: Agilent-028803 Custom Whole Sheep Genome Oligo Microarray. 2012. http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GPL16283. Accessed 16 Nov 2012.
- Yang Y, Du AF: Gene Expression Omnlbus: Profiling of differentially expressed genes in sheep T lymphocytes response to Haemonchus contortus infection. 2012. http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE42302. Accessed 16 Nov 2012.
- Webster LM, Paterson S, Mougeot F, Martinez-Padilla J, Piertney SB. Transcriptomic response of red grouse to gastro-intestinal nematode parasites and testosterone: implications for population dynamics. Mol Ecol. 2011;20(5):920–31.View ArticlePubMedGoogle Scholar
- Xiang LX, He D, Dong WR, Zhang YW, Shao JZ. Deep sequencing-based transcriptome profiling analysis of bacteria-challenged Lateolabrax japonicus reveals insight into the immune-relevant genes in marine fish. BMC Genomics. 2010;11:472.View ArticlePubMed CentralPubMedGoogle Scholar
- Castilho PC, Buckley BA, Somero G, Block BA. Heterologous hybridization to a complementary DNA microarray reveals the effect of thermal acclimation in the endothermic bluefin tuna (Thunnus orientalis). Mol Ecol. 2009;18(10):2092–102.View ArticlePubMedGoogle Scholar
- MacKinnon KM, Burton JL, Zajac AM, Notter DR. Microarray analysis reveals difference in gene expression profiles of hair and wool sheep infected with Haemonchus contortus. Vet Immunol Immunopathol. 2009;130(3–4):210–20.View ArticlePubMedGoogle Scholar
- Schwarz EM, Korhonen PK, Campbell BE, Young ND, Jex AR, Jabbar A, et al. The genome and developmental transcriptome of the strongylid nematode Haemonchus contortus. Genome Biol. 2013;14(8):R89.View ArticlePubMed CentralPubMedGoogle Scholar
- Sprent J, Surh CD. Normal T cell homeostasis: the conversion of naive cells into memory-phenotype cells. Nat Immunol. 2011;12(6):478–84.View ArticlePubMed CentralPubMedGoogle Scholar
- Bhattacharya R, Wang E, Dutta SK, Vohra PK EG, Prakash YS, Mukhopadhyay D. NHERF-2 maintains endothelial homeostasis. Blood. 2012;119(20):4798–806.View ArticlePubMed CentralPubMedGoogle Scholar
- Subramanian H, Gupta K, Ali H. Roles for NHERF1 and NHERF2 on the regulation of C3a receptor signaling in human mast cells. PLoS One. 2012;7(12), e51355.View ArticlePubMed CentralPubMedGoogle Scholar
- Isom LL. Sodium channel beta subunits: anything but auxiliary. Neuroscientist. 2001;7(1):42–54.View ArticlePubMedGoogle Scholar
- Calhoun JD, Isom LL. The role of non-pore-forming beta subunits in physiology and pathophysiology of voltage-gated sodium channels. Handb Exp Pharmacol. 2014;221:51–89.View ArticlePubMedGoogle Scholar
- Bolivar J, Garcia-Cozar FJ, Astola A, Iglesias C, Pendon C, Valdivia MM. Genomic structure and chromosome location of the human gene encoding the zinc finger autoantigen ZNF330. Cytogenet Cell Genet. 2001;93(3–4):234–8.PubMedGoogle Scholar
- De Melo IS, Iglesias C, Benitez-Rondan A, Medina F, Martinez-Barbera JP, Bolivar J. NOA36/ZNF330 is a conserved cystein-rich protein with proapoptotic activity in human cells. Biochim Biophys Acta. 2009;1793(12):1876–85.View ArticlePubMedGoogle Scholar
- Hennet T. The galactosyltransferase family. Cell Mol Life Sci. 2002;59(7):1081–95.View ArticlePubMedGoogle Scholar
- Sato T, Furukawa K. Expression and Transcriptional Regulation of β4-Galactosyltransferase Genes in Cancer. In: Endo T, Seeberger PH, Hart GW, Wong C-H, Taniguchi N, editors. Glycoscience: Biology and Medicine. Japan: Springer; 2014. p. 1–5.Google Scholar
- Bangert I, Tumulka F, Abele R. The lysosomal polypeptide transporter TAPL: more than a housekeeping factor? Biol Chem. 2011;392(1–2):61–6.PubMedGoogle Scholar
- Pohl A, Devaux PF, Herrmann A. Function of prokaryotic and eukaryotic ABC proteins in lipid transport. Biochim Biophys Acta. 2005;1733(1):29–52.View ArticlePubMedGoogle Scholar
- Fujimoto Y, Kamakura A, Motohashi Y, Ohashi-Kobayashi A, Maeda M. Transporter associated with antigen processing-like (ABCB9) stably expressed in Chinese hamster ovary-K1 cells is sorted to the microdomains of lysosomal membranes. Biol Pharm Bull. 2011;34(1):36–40.View ArticlePubMedGoogle Scholar
- Lawrence CE, Paterson YY, Wright SH, Knight PA, Miller HR. Mouse mast cell protease-1 is required for the enteropathy induced by gastrointestinal helminth infection in the mouse. Gastroenterology. 2004;127(1):155–65.View ArticlePubMedGoogle Scholar
- Caughey GH. Mast cell proteases as protective and inflammatory mediators. Adv Exp Med Biol. 2011;716:212–34.View ArticlePubMed CentralPubMedGoogle Scholar
- Gianni M, Ponzanelli I, Mologni L, Reichert U, Rambaldi A, Terao M, et al. Retinoid-dependent growth inhibition, differentiation and apoptosis in acute promyelocytic leukemia cells. Expression and activation of caspases. Cell Death Differ. 2000;7(5):447–60.View ArticlePubMedGoogle Scholar
- Rajsbaum R, Fici D, Fraser PA, Flores-Villanueva PO, Awdeh ZL. Polymorphism of the human retinoid X receptor beta and linkage disequilibrium with HLA-DPB1. Tissue Antigens. 2001;58(1):24–9.View ArticlePubMedGoogle Scholar
- Jia WH, Zhang B, Matsuo K, Shin A, Xiang YB, Jee SH, et al. Genome-wide association analyses in East Asians identify new susceptibility loci for colorectal cancer. Nat Genet. 2013;45(2):191–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Ye F, Zhang Y, Liu Y, Yamada K, Tso JL, Menjivar JC, et al. Protective properties of radio-chemoresistant glioblastoma stem cell clones are associated with metabolic adaptation to reduced glucose dependence. PLoS One. 2013;8(11):e80397.View ArticlePubMed CentralPubMedGoogle Scholar
- Maine GN, Burstein E. COMMD proteins: COMMing to the scene. Cell Mol Life Sci. 2007;64(15):1997–2005.View ArticlePubMed CentralPubMedGoogle Scholar
- Jin P, Gao Y, Chen L, Ma F. Cloning and characterization of a COMMD4 gene from amphioxus (Branchiostoma belcheri): an insight into the function and evolution of COMMD4. Immunol Lett. 2012;148(2):110–6.View ArticlePubMedGoogle Scholar
- Bartuzi P, Hofker MH, Van de Sluis B. Tuning NF-kappaB activity: a touch of COMMD proteins. Biochim Biophys Acta. 2013;1832(12):2315–21.View ArticlePubMedGoogle Scholar
- Burstein E, Hoberg JE, Wilkinson AS, Rumble JM, Csomos RA, Komarck CM, et al. COMMD proteins, a novel family of structural and functional homologs of MURR1. J Biol Chem. 2005;280(23):22222–32.View ArticlePubMedGoogle Scholar
- Heyen CA, Tagliabracci VS, Zhai L, Roach PJ. Characterization of mouse UDP-glucose pyrophosphatase, a Nudix hydrolase encoded by the Nudt14 gene. Biochem Biophys Res Commun. 2009;390(4):1414–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Talieri M, Alexopoulou DK, Scorilas A, Kypraios D, Arnogiannaki N, Devetzi M, et al. Expression analysis and clinical evaluation of kallikrein-related peptidase 10 (KLK10) in colorectal cancer. Tumour Biol. 2011;32(4):737–44.View ArticlePubMedGoogle Scholar
- Ho MF, Rose’Meyer RB. Vascular adenosine receptors; potential clinical applications. Curr Vasc Pharmacol. 2013;11(3):327–37.View ArticlePubMedGoogle Scholar
- Warmflash A, Zhang Q, Sorre B, Vonica A, Siggia ED, Brivanlou AH. Dynamics of TGF-beta signaling reveal adaptive and pulsatile behaviors reflected in the nuclear localization of transcription factor Smad4. Proc Natl Acad Sci U S A. 2012;109(28):E1947–56.View ArticlePubMed CentralPubMedGoogle Scholar
- Menon MB, Sawada A, Chaturvedi A, Mishra P, Schuster-Gossler K, Galla M, et al. Genetic deletion of SEPT7 reveals a cell type-specific role of septins in microtubule destabilization for the completion of cytokinesis. PLoS Genet. 2014;10(8):e1004558.View ArticlePubMed CentralPubMedGoogle Scholar
- Abo A, Qu J, Cammarano MS, Dan C, Fritsch A, Baud V, et al. PAK4, a novel effector for Cdc42Hs, is implicated in the reorganization of the actin cytoskeleton and in the formation of filopodia. EMBO J. 1998;17(22):6527–40.View ArticlePubMed CentralPubMedGoogle Scholar
- Dan C, Kelly A, Bernard O, Minden A. Cytoskeletal changes regulated by the PAK4 serine/threonine kinase are mediated by LIM kinase 1 and cofilin. J Biol Chem. 2001;276(34):32115–21.View ArticlePubMedGoogle Scholar
- Da Silva CJ, Miranda Y, Leonard N, Ulevitch R. SGT1 is essential for Nod1 activation. Proc Natl Acad Sci U S A. 2007;104(16):6764–9.View ArticleGoogle Scholar
- Mayor A, Martinon F, De Smedt T, Petrilli V, Tschopp J. A crucial function of SGT1 and HSP90 in inflammasome activity links mammalian and plant innate immune responses. Nat Immunol. 2007;8(5):497–503.View ArticlePubMedGoogle Scholar
- Ali A, Hoeflich KP, Woodgett JR. Glycogen synthase kinase-3: properties, functions, and regulation. Chem Rev. 2001;101(8):2527–40.View ArticlePubMedGoogle Scholar
- Aparicio IM, Garcia-Herreros M, Fair T, Lonergan P. Identification and regulation of glycogen synthase kinase-3 during bovine embryo development. Reproduction. 2010;140(1):83–92.View ArticlePubMedGoogle Scholar
- Pelicci G, Lanfrancone L, Grignani F, McGlade J, Cavallo F, Forni G, et al. A novel transforming protein (SHC) with an SH2 domain is implicated in mitogenic signal transduction. Cell. 1992;70(1):93–104.View ArticlePubMedGoogle Scholar
- Shih HJ, Chen HH, Chen YA, Wu MH, Liou GG, Chang WW, et al. Targeting MCT-1 oncogene inhibits Shc pathway and xenograft tumorigenicity. Oncotarget. 2012;3(11):1401–15.PubMed CentralPubMedGoogle Scholar
- Zheng Y, Zhang C, Croucher DR, Soliman MA, St-Denis N, Pasculescu A, et al. Temporal regulation of EGF signalling networks by the scaffold protein Shc1. Nature. 2013;499(7457):166–71.View ArticlePubMedGoogle Scholar
- Von Bubnoff D, Novak N, Kraft S, Bieber T. The central role of FcepsilonRI in allergy. Clin Exp Dermatol. 2003;28(2):184–7.View ArticlePubMedGoogle Scholar
- Siraganian RP. Mast cell signal transduction from the high-affinity IgE receptor. Curr Opin Immunol. 2003;15(6):639–46.View ArticlePubMedGoogle Scholar
- Kadarmideen HN, Watson-Haigh NS, Andronicos NM. Systems biology of ovine intestinal parasite resistance: disease gene modules and biomarkers. Mol Biosyst. 2011;7(1):235–46.View ArticlePubMedGoogle Scholar
- Kijimoto-Ochiai S. CD23 (the low-affinity IgE receptor) as a C-type lectin: a multidomain and multifunctional molecule. Cell Mol Life Sci. 2002;59(4):648–64.View ArticlePubMedGoogle Scholar
- Henningsson F, Ding Z, Dahlin JS, Linkevicius M, Carlsson F, Gronvik KO, et al. IgE-mediated enhancement of CD4+ T cell responses in mice requires antigen presentation by CD11c + cells and not by B cells. PLoS One. 2011;6(7):e21760.View ArticlePubMed CentralPubMedGoogle Scholar
- Liang HE, Reinhardt RL, Bando JK, Sullivan BM, Ho IC, Locksley RM. Divergent expression patterns of IL-4 and IL-13 define unique functions in allergic immunity. Nat Immunol. 2012;13(1):58–66.View ArticleGoogle Scholar
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.