Comparative gene expression profiling in two congenic mouse strains following Bordetella pertussis infection
© Banus et al; licensee BioMed Central Ltd. 2007
Received: 15 June 2007
Accepted: 12 October 2007
Published: 12 October 2007
Susceptibility to Bordetella pertussis infection varies widely. These differences can partly be explained by genetic host factors. HcB-28 mice are more resistant to B. pertussis infection than C3H mice, which could partially be ascribed to the B. pertussis susceptibility locus-1 (Bps1) on chromosome 12. The presence of C57BL/10 genome on this locus instead of C3H genome resulted in a decreased number of bacteria in the lung. To further elucidate the role of host genetic factors, in particular in the Bps1 locus, in B. pertussis infection, and to identify candidate genes within in this region, we compared expression profiles in the lungs of the C3H and HcB-28 mouse strains following B. pertussis inoculation. Twelve and a half percent of the genomes of these mice are from a different genetic background.
Upon B. pertussis inoculation 2,353 genes were differentially expressed in the lungs of both mouse strains. Two hundred and six genes were differentially expressed between the two mouse strains, but, remarkably, none of these were up- or down-regulated upon B. pertussis infection. Of these 206 genes, 17 were located in the Bps1 region. Eight of these genes, which showed a strong difference in gene expression between the two mouse strains, map to the immunoglobulin heavy chain complex (Igh).
Gene expression changes upon B. pertussis infection are highly identical between the two mouse strains despite the differences in the course of B. pertussis infection. Because the genes that were differentially regulated between the mouse strains only showed differences in expression before infection, it appears likely that such intrinsic differences in gene regulation are involved in determining differences in susceptibility to B. pertussis infection. Alternatively, such genetic differences in susceptibility may be explained by genes that are not differentially regulated between these two mouse strains. Genes in the Igh complex, among which Igh-1a/b, are likely candidates to explain differences in susceptibility to B. pertussis. Thus, by microarray analysis we significantly reduced the number of candidate susceptibility genes within the Bps1 locus. Further work should establish the role of the Igh complex in B. pertussis infection.
Bordetella pertussis is a gram-negative bacterium that can cause the respiratory disease known as pertussis or whooping cough in humans. Susceptibility to this disease and its course vary widely between individuals . We have previously shown that genetically divergent mouse strains differ in their response to B. pertussis infection, underlining that infection is influenced by host genetic factors[2, 3]. In addition, a role of several host genetic loci in the course of B. pertussis infection has been indicated, such as the toll-like receptor 4 (Tlr4) gene [2–5], the interferon gamma receptor gene , and three novel loci, B. pertussis susceptibility locus 1, 2, and 3 (Bps1, 2, and 3)  that showed linkage with the severity of infection.
Although Bps-1 has not yet been validated, we believe that the significance of this locus warrants further study. To further elucidate the role of host genetic factors, in particular the Bps1 locus, in B. pertussis infection, and to identify candidate genes within this region, we studied expression profiles in the lungs of mice following B. pertussis inoculation. The traditional approach for mapping genes in susceptibility loci is a combination of positional cloning and linkage analysis [11, 12]. Although this strategy has proven to be effective [7, 13], the approach is quite costly and animal-consuming. Previously we identified changes in gene expression in the lungs of C3H/DISnA mice after B. pertussis infection, and we especially focused on differentially expressed genes in the lungs of infected and non-infected mice located in Bps1,-2, and -3 . We found that the expression of 1,841 genes was significantly changed upon B. pertussis infection. These genes are involved in various immune-related processes, such as the acute-phase response, antigen presentation, cytokine production, inflammation, and apoptosis. Nine of the differentially expressed genes are located in Bps1, 13 are located in Bps-2, and 62 are located in Bps-3.
In the present study we compared the gene expression profiles in the lungs of two mouse strains, i.e. HcB-28/Dem and C3H, which showed a different course of B. pertussis infection, in order to further identify candidate susceptibility genes without the need for positional cloning. We hypothesized that the phenotypic differences displayed by these mice in the response to B. pertussis can partly be explained by a different gene expression profile between the mouse strains, and that this approach could lead to the identification of candidate genes affecting the course of infection. Using this approach we reduced the number of candidate susceptibility genes within the Bps1 locus.
Forty-eight female HcB-28/Dem and 48 female C3H/DISnA (the background strain of the HcB-28/Dem) mice were intranasally inoculated with 2*107 colony forming units (CFU) of the streptomycin-resistant Tohama strain of B. pertussis (B213) in 40 μl Verwey medium (The Netherlands Vaccine Institute, Bilthoven, the Netherlands), or as a control with Verwey Medium only (total of 96 mice). One, three, and five days after inoculation, 8 infected and 8 control mice were euthanized. To remove blood from the lungs, mice were perfused with phosphate-buffered saline (PBS, Tritium Microbiology, Veldhoven, the Netherlands). Subsequently, the lungs and trachea were collected [15, 16]. For RNA extraction, the right lung was collected in RNA stabilization reagent (RNAlater, Qiagen, Venlo, the Netherlands). For histological examination, the left lung was fixed intratracheally using 4% formalin.
The number of viable B. pertussis bacteria was determined in the trachea to confirm a proper infection  (Note that tracheal counts are not representative for bacterial clearance from the lungs.) To this end, approximately one centimeter of the trachea was collected in 500 μl Verwey medium. Bacterial suspensions were diluted in Verwey medium and the number of CFU was determined by plating on Bordet Gengou agar supplemented with 15% sheep blood and 30 μg/ml streptomycin (Tritium Microbiology). Plates were incubated for 4 days at 35°C and the resulting colonies were counted using a ProtoCOL Colony counter (Synbiosis, Cambridge, United Kingdom).
All B. pertussis-infected mice had between 102 and 105 CFU in the trachea during the first five days after inoculation, confirming an actual infection (data not shown).
For the analysis of intrinsic immunological and cellular differences between the two mouse strains, independent from infection, blood and spleen samples were collected from five untreated female mice of both strains.
Breeding pairs of C3H/DISnA and HcB-28/Dem were supplied by the Department of Molecular and Cellular Biology, Roswell Park Cancer Institute, Buffalo, New York. Female mice were bred to the appropriate numbers at our animal testing facility in Bilthoven, the Netherlands. Mice were acclimatized for at least one week before the start of the experiments. Mice received standard laboratory food (SRM-A, Hope Farms, Woerden, the Netherlands) and tap water ad libitum. All animal experiments were approved by the Institute's Animal Ethics Committee and were performed according to NIH guidelines  and Dutch legislation.
Clinical and pathological examinations
Mice were weighed before inoculation, and subsequently every day after inoculation to determine the relative change in weight. Lung weights were determined post mortem as quantitative parameter for lung inflammation. Formalin-fixed lungs were embedded in paraplast (Monoject Inc., St Louis, MO). Sections (5 μm) were stained with hematoxylin-eosin. Lung lesions were examined for infiltration of inflammatory cells in the peribronchiolar space (peribronchiolitis), infiltration of inflammatory cells in the alveoli (alveolitis), infiltration of inflammatory cells in the perivascular space (perivasculitis), hypertrophy of mucus-producing glands, free protein (exudate), and eosinophilia. Lung lesions were scored semi-quantitatively as absent, minimal, slight, moderate, marked, or strong, as previously described .
Enumeration of cellular subsets
Spleens were harvested and weighed. Cell suspensions were prepared and the number of nucleated cells per spleen was determined using a Coulter Counter Z2 (Beckman Coulter, Mijdrecht, the Netherlands). The percentage of B- and T-cells, as well as T-cell subsets were determined using a fluorescence-activated cell sorter (FACS-Calibur, BD Biosciences, Alphen aan den Rijn, the Netherlands). T-cells were detected with phycoerythrin-labeled (PE) anti-CD3ε antibodies (Molecular probes, Invitrogen, Breda, the Netherlands). CD4+ T-cells were determined using allophycocyanin-labeled (APC) antibodies (Molecular probes) and CD8+T-cells were determined using fluorescein isothiocyanate-labeled (FITC) antibodies (Molecular probes). B-cells were determined using PE-labeled anti-CD19 antibodies (Molecular probes).
Splenocyte proliferation test
The cell suspensions were used to measure the spontaneous (medium only) or stimulus-induced 3H-thymidine uptake. The cells were stimulated by adding 10 μg/ml Concavalin A (Con A, T-cell stimulus), 70 μg/ml lipopolysacharide (LPS, B-cell stimulus) or 15 μg/ml Lectin (B- and T-cell stimulus) to the culture medium. Cells were incubated for 72 hr at 37°C. 3H-thymidine was present during the final 24 hr of culture. See  for details.
The morphologic hemogram was determined using an ADVIA 120 Hematology system (Multispecies analyzer, Bayer, Mijdrecht, the Netherlands).
Blood was collected in Vacuette Minicollect tubes (Greiner bio-one, Alphen a/d Rijn, the Netherlands), and serum was removed by centrifugation. Immunoglobulins were determined using the multiplex Beadlyte Mouse Immunoglobulin Isotyping Kit (Millipore, Billerica, MA) for the Luminex platform (Luminex, Oosterhoud, the Netherlands), as prescribed by the manufacturer.
Microarray analysis experiments were performed as described previously . Briefly, total RNA was extracted from lungs and amplified using the Amino Allyl MessageAmp II aRNA kit (Ambion Inc., Austin, TX). RNA amplification was performed to obtain more nucleic acid for labeling. This results in a stronger fluorescence signal and a better signal/noise ratio. Because of this, less experimental samples (and therefore animals) are needed to obtain sufficient statistical power. RNA samples from individual mice were labeled with Cy3. A common reference containing a RNA pool of all samples isolated was labeled with Cy5.
Microarray slides were spotted at the Microarray Department of the University of Amsterdam. The slides contain 21,997 65-mer oligonucleotides from the Sigma-Compugen Mouse oligonucleotide library, 192 additional 65-mer oligonucleotides of genes of which 185 are located in the Bps1 region spotted in duplicate (Isogen Life science, IJsselstein, the Netherlands), as well as appropriate controls and blank spots. The Bps-1 oligonucleotides were annotated according to NCBI mouse genome build 34.1; all other oligonucleotides were annotated according to the Sigma-Compugen Mouse oligonucleotide library.
At least 5 slides per group were used. Raw microarray signal data were normalized and analyzed in the R statistical software environment [21–23]. Significance of differences in gene expression between the experimental groups was calculated in R with ANOVA. The false discovery rate (FDR) was calculated according to Benjamini and Hochberg . Gene expression data were visualized by hierarchical clustering (using Euclidian distance and Ward linkage) using GeneMaths (Applied Maths, St-Martens-Latem, Belgium). Gene categories were defined by Gene Ontology (GO) , Classification and enrichment according to GO categories or location were determined using DAVID/EASE [26, 27], using EASE-scores as P- value for enrichment.
Response of C3H and HcB mice to infection
We have previously examined clearance of B. pertussis from the lungs of infected RCS mice, and established that HcB-28 mice contained lower numbers of bacteria in their lungs seven days post inoculation compared to C3H mice (Figure 1b) . To examine whether the genetic differences between "resistant" HcB-28 mice and "susceptible" C3H mice also result in a different transcriptional response to B. pertussis infection, we compared gene expression profiles of the lungs of infected and uninfected HcB-28 and C3H mice [2, 14]. mRNA transcription profiles were determined in at least five animals per group.
The 206 genes that were differentially expressed between the two strains were all unaffected by infection, but displayed a difference in basal expression level. Of the 206 strain-dependent genes, 107 genes were expressed at a higher level in C3H/DISnA mice (group A) and 99 genes were expressed at a higher level in the HcB-28/Dem mice (group B). Of the 2,353 infection-regulated genes, 963 were down-regulated upon infection (group C), and 1,390 were up-regulated upon infection (group D). All genes that were affected by B. pertussis infection followed a similar trend in time. One day post-inoculation differential gene expression was already observed, but the effect was stronger on days three and five post-inoculation. Thus there was no identifiable subset of genes that displayed a different kinetic trend in gene expression.
Gene-ontology based classification of genes with different expression
Total on array
Mus musculus 12 a
Bps-1 (Mus musculus 12)
Listed in any of the above
Candidate susceptibility genes
RIKEN cDNA 2610204M08 gene
cDNA sequence BC022687
similar to monoclonal antibody heavy chain
similar to immunoglobulin heavy chain
similar to Ig H-chain
similar to IgE antibody heavy chain (VDJ)
similar to Igh-VJ558 protein
similar to heavy chain V region VH558 A1/A4 precursor
similar to Ig heavy chain variable region precursor
RIKEN cDNA 1700001K19 gene
protein phosphatase 1, regulatory (inhibitor) subunit 13B
adenylosuccinate synthetase like 1
thymoma viral proto-oncogene 1
expressed sequence AI450948
immunoglobulin heavy chain 1a (serum IgG2a)
Of the 2,353 genes that were differentially expressed after B. pertussis infection, 1,702 genes were identical to the genes we have described previously to be regulated upon B. pertussis infection . Thus by examining additional mice, we identified approximately 650 additional genes regulated by B. pertussis, most of which were weakly up- or down-regulated. Most up-regulated genes upon B. pertussis infection (group D) are involved in immune- and inflammation-related processes or in generic processes, while most down-regulated genes (group C) are involved in non-immune processes such as muscle contraction.
Hematologic and immunologic analyses in non-infected C3H and HcB-28 mice
To investigate possible differences in the basal immunological parameters between the two strains, we analyzed hematologic parameters, mitogen-induced splenocyte proliferation, and lymphocyte subset distributions. Because the two mouse strains displayed clear differences in basal gene expression, we wanted to examine if these differences led to alterations in blood cell composition or basal immune status. We observed a slightly higher number of platelets (P = 0.0003) in the C3H/DISnA mice (1.1*1012/l) compared to the HcB-28/Dem mice (9.6*1011/l), and a slightly higher number of white blood cells (P = 0.006) in the C3H/DISnA mice (7.3*109/l) compared to the HcB-28/Dem mice (6.5*109/l), but no differences in the numbers of red blood cells, reticulocytes, lymphocytes, neutrophils, eosinophils, and basophils (data not shown).
Splenocyte proliferation was determined by mitogen-induced 3H-thymidine uptake. There was no difference in proliferation between the two strains after ex vivo stimulation with Con A (T-cell stimulus), LPS (B-cell stimulus) or Lectin (B- and T-cell stimulus) (data not shown). We observed no differences in the number of CD19+ (B), CD3ε+ (T), CD4+(Th) and CD8+ (CTL) cells (data not shown).
Clinical and pathological findings in B. pertussis-infected C3H and HcB-28 mice
As a quantitative endpoint for inflammation, we determined lung weights relative to body weight. All mock-inoculated mice had relative lung weights of approximately one percent five days after inoculation, while B. pertussis- infected mice had relative lung weights of up to three percent. This difference in relative lung weights due to B. pertussis infection is significant from day one till day five post-inoculation, but was similar in both mouse strains (data not shown).
Studying genetic differences in susceptibility to B. pertussis infection may point to novel insights in the pathogenesis of this infection. We have previously identified B. pertussis susceptibility locus-1 (Bps1) in HcB-28 mice . The Bps1 locus is located on chromosome 12, spanning a region of 185 genes and has a dominant positive effect on the clearance of B. pertussis from the lung. In this study we examined gene expression profiles in HcB-28 and C3H mice, which differ in their susceptibility to B. pertussis. Twelve and a half percent of the genomes of these mice are from a different genetic background [9, 10]. The traditional approach for identification of relevant genes in susceptibility loci is a combination of positional cloning and linkage analysis [11, 12]. This method has proven to be effective [7, 13], but has disadvantages. We have therefore chosen for an alternative strategy. Using this approach we attempted to identify candidate susceptibility genes that control the difference between these two mouse strains. We hypothesized that the difference in susceptibility to B. pertussis infection could (partly) be explained by a different gene expression profile between the mouse strains.
We have previously shown that B. pertussis infection in C3H mice induces a wide transcriptional response, which appears to be partly specific for B. pertussis and partly non-specific . This study revealed that 1,841 genes are differentially expressed in the lungs of mice after B. pertussis inoculation, and most up-regulated genes are involved in immune- and inflammation-related processes or in generic processes, while most down-regulated genes are involved in non-immune processes. In the present study, we found that HcB-28 and C3H mice showed a similar gene expression profile upon infection and identified approximately 650 additional genes regulated by B. pertussis, most of which were weakly up- or down-regulated. However, the gene expression profiles and enrichment for GO categories were identical as described in the previous study. A substantial number of genes and pathways suggest a central role of PMN recruitment and activation in the pathogenesis of B. pertussis infection. The transcriptional profiles further indicate in particular the significance of TLR activation and apoptosis . The reason for the detection of the 650 additional genes, besides the strain differences, is that by adding an extra mouse strain the number of samples doubled increasing the power of detection. Hundred and thirty-nine genes which we have described previously to be regulated upon B. pertussis infection were not detected in the present study. These genes were borderline significantly regulated (median FDR of 0.03) and only slightly induced (median 1.4-fold). These genes are therefore probably less important in the host response to B. pertussis infection. The finding that this list of 139 genes does not show significant enrichment for any GO-term, including immunological terms, corroborates this.
Although the two mouse strains differ in 12.5% of their genome (12.5% of the genomes of these mice are from a different genetic background), we observed no marked differences in their phenotypical characteristics other than the previously observed difference in bacterial numbers in the lungs after infection (Figure 1b). C3H mice did have slightly higher numbers of circulating platelets and white blood cells compared to the HcB mice, but the cellular proportions, as well as the proliferation of splenocytes was identical for both strains. There was also no significant difference observed in body weight, lung weight and histopathological findings in response to B. pertussis infection between the two mouse strains. The major difference is that HcB-28 mice did not have detectable IgG2a serum levels.
We observed 206 genes that were differentially expressed between the two mouse strains, but these genes were identically expressed in mock- or B. pertussis-inoculated mice. The majority of these genes (65%) are unannotated. These unannotated genes tend to have no GO functional annotation, because they are not "regular" protein-coding genes and many of them are not (sufficiently) mapped to a chromosomal locus to warrant including them under chromosome 12 or Bps-1. The mentioned 23 genes, which were mapped to chromosome 12, are significant at P = 5.29e-014 (Fisher exact probability, Bonferroni correction for multiple testing). The same test applied to all other mouse chromosomes yielded P values > 0.05. Twenty-three out of the 206 genes were located on chromosome 12, which can be explained by the fact that approximately one fifth of the genetic variation between the two mouse strains is due to variation on chromosome 12 . Because the genes that were differentially regulated between the two mouse strains only showed differences in expression before infection, it appears likely that such intrinsic differences in gene regulation are involved in determining differences in susceptibility to B. pertussis infection. Alternatively, such genetic differences may be explained by genes that are not differentially regulated between these two strains of mice, or by processes at present not fully characterized and possibly involving differential expression of genes by mechanisms such as microRNA's. Remarkably, 17 of these genes were located in the Bps1 region, 8 of which mapped to the Igh complex. Among these 8 genes were the Igh-1 gene and genes that encode for Ig heavy chain variable regions. We observed a significantly higher expression (up to 2.8-fold) of the gene variant of Igh-1 isotype b in HcB mice compared to C3H mice. The Igh locus is genetically polymorphic and very complex . The Igh-1 gene exists in 2 major genetic variants (Igh-1a and Igh-1b) with 83.8% similarity . The Igh-1a allele codes for the heavy chain of IgG2a while the Igh-1b allele codes for the heavy chain of IgG2c . Mouse strains such as C57BL/6 and C57BL/10 (the donor strain of the HcB mice) only contain the gene variant Igh-1b and are therefore incapable of producing IgG2a, while mouse strains such as BALB/c only contain the gene variant of Igh-1a and are therefore incapable of producing IgG2c [29, 30]. The oligo for Igh-1 spotted on the microarray, was designed based on accession number XM_484178 annotated as Igh-1a. Because this annotation is based on the reference C57BL/6 mice, this oligo can be considered as Igh-1b annotated. In the sera of HcB mice we detected no IgG2a while C3H mice had significant titers of IgG2a, thereby confirming the expression results. It is tempting to speculate if and how genes within the Igh complex may affect differences early in the course of B. pertussis infection. Possible mechanisms may include differences in transcriptional gene regulation affecting immune responsiveness, different function of the IgG2a or c isotypes, or different usage of V chains. This latter possibility might imply the existence of "natural antibodies" reacting with B. pertussis epitopes. It has previously been shown that genes within the Igh-1 locus are predominantly associated with the course of a herpes simplex virus type-1 (HSV-1) infection in mice by an unknown mechanism [31–33]. Pro-inflammatory cytokines such as IL-1β, IL-4, IL-6 and IL-7 participate in this infection . Interestingly, natural killer cell activity appears to be regulated by the Igh-1 locus but could not simply explain the differences in HSV-1 susceptibility . Igh-linked genes have further been implicated in T suppressor cell activity [35, 36].
In this study we reduced the number of candidate susceptibility genes within the Bps1 locus by microarray analysis. Gene expression changes upon B. pertussis infection appear highly identical between C3H and HcB-28 mouse strains despite the different course of B. pertussis infection in these strains. Because the genes that were differentially regulated between the mouse strains only showed differences in expression before infection, it appears likely that such intrinsic differences in gene regulation are involved in determining differences in susceptibility to B. pertussis infection. Alternatively, such genetic differences in susceptibility may be explained by genes that are not differentially regulated between these two mouse strains or by processes other than differential gene expression. Genes in the Igh complex, among which Igh-1, may be likely candidates to explain differences in susceptibility to B. pertussis. Further work should establish the role of the Igh complex in B. pertussis infection and determine its mode of action.
Raw data as well as the detailed description of the experiment was uploaded to the freely accessible online database ArrayExpress . Expression levels of all 2,559 genes are presented in the additional files.
We would like to thank Yvonne Wallbrink, Liset de la Fonteyne and Sisca de Vlugt- van den Koedijk for the hematology, FACS analysis and histotechnical contributions to this paper. Furthermore we would like to thank all biotechnicians of our animal facility, for facilitating and performing the animal experiment.
- Carbonetti NH: Immunomodulation in the pathogenesis of Bordetella pertussis infection and disease. Curr Opin Pharmacol. 2007Google Scholar
- Banus S, van Kranen HJ, Mooi FR, Hoebee B, Nagelkerke NJ, Demant P, Kimman TG: Genetic Control of Bordetella pertussis Infection: Identification of Susceptibility Loci Using Recombinant Congenic Strains of Mice. Infect Immun. 2005, 73: 741-747. 10.1128/IAI.73.2.741-747.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Banus S, Vandebriel RJ, de Ruiter H, Dormans JA, Nagelkerke NJ, Mooi FR, Hoebee B, van Kranen HJ, Kimman TG: Host Genetics of Bordetella pertussis Infection in Mice: Significance of Toll-Like Receptor 4 in Genetic Susceptibility and Pathobiology. Infect Immun. 2006, 74: 2596-2605. 10.1128/IAI.74.5.2596-2605.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Mann PB, Kennett MJ, Harvill ET: Toll-Like Receptor 4 Is Critical to Innate Host Defense in a Murine Model of Bordetellosis. J Infect Dis. 2004, 189: 833-836. 10.1086/381898.View ArticlePubMedGoogle Scholar
- Mann PB, Wolfe D, Latz E, Golenbock D, Preston A, Harvill ET: Comparative toll-like receptor 4-mediated innate host defense to Bordetella infection. Infect Immun. 2005, 73: 8144-8152. 10.1128/IAI.73.12.8144-8152.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Mahon BP, Sheahan BJ, Griffin F, Murphy G, Mills KH: Atypical disease after Bordetella pertussis respiratory infection of mice with targeted disruptions of interferon-gamma receptor or immunoglobulin mu chain genes. J Exp Med. 1997, 186: 1843-1851. 10.1084/jem.186.11.1843.PubMed CentralView ArticlePubMedGoogle Scholar
- Demant P: Cancer susceptibility in the mouse: genetics, biology and implications for human cancer. Nat Rev Genet. 2003, 4: 721-734. 10.1038/nrg1157.View ArticlePubMedGoogle Scholar
- Demant P, Hart AA: Recombinant congenic strains--a new tool for analyzing genetic traits determined by more than one gene. Immunogenetics. 1986, 24: 416-422. 10.1007/BF00377961.View ArticlePubMedGoogle Scholar
- Groot PC, Moen CJ, Dietrich W, Stoye JP, Lander ES, Demant P: The recombinant congenic strains for analysis of multigenic traits: genetic composition. FASEB J. 1992, 6: 2826-2835.PubMedGoogle Scholar
- Stassen AP, Groot PC, Eppig JT, Demant P: Genetic composition of the recombinant congenic strains. Mamm Genome. 1996, 7: 55-58. 10.1007/s003359900013.View ArticlePubMedGoogle Scholar
- Rannala B: Finding genes influencing susceptibility to complex diseases in the post-genome era. Am J Pharmacogenomics. 2001, 1: 203-221. 10.2165/00129785-200101030-00005.View ArticlePubMedGoogle Scholar
- Ruivenkamp CA, van Wezel T, Zanon C, Stassen AP, Vlcek C, Csikos T, Klous AM, Tripodis N, Perrakis A, Boerrigter L, Groot PC, Lindeman J, Mooi WJ, Meijjer GA, Scholten G, Dauwerse H, Paces V, Van Zandwijk N, Van Ommen GJ, Demant P: Ptprj is a candidate for the mouse colon-cancer susceptibility locus Scc1 and is frequently deleted in human cancers. Nat Genet. 2002, 31: 295-300. 10.1038/ng903.View ArticlePubMedGoogle Scholar
- Ruivenkamp C, Hermsen M, Postma C, Klous A, Baak J, Meijer G, Demant P: LOH of PTPRJ occurs early in colorectal cancer and is associated with chromosomal loss of 18q12-21. Oncogene. 2003, 22: 3472-3474. 10.1038/sj.onc.1206246.View ArticlePubMedGoogle Scholar
- Banus S, Pennings J, Vandebriel R, Wester P, Breit T, Mooi F, Hoebee B, Kimman T: Lung response to Bordetella pertussis infection in mice identified by gene-expression profiling. Immunogenetics. 2007, 59: 555-564. 10.1007/s00251-007-0227-5.PubMed CentralView ArticlePubMedGoogle Scholar
- Willems RJ, Kamerbeek J, Geuijen CA, Top J, Gielen H, Gaastra W, Mooi FR: The efficacy of a whole cell pertussis vaccine and fimbriae against Bordetella pertussis and Bordetella parapertussis infections in a respiratory mouse model. Vaccine. 1998, 16: 410-416. 10.1016/S0264-410X(97)80919-X.View ArticlePubMedGoogle Scholar
- King AJ, Berbers G, van Oirschot HF, Hoogerhout P, Knipping K, Mooi FR: Role of the polymorphic region 1 of the Bordetella pertussis protein pertactin in immunity. Microbiology. 2001, 147: 2885-2895.View ArticlePubMedGoogle Scholar
- Verwey WF, Thiele EH, Sage DN, Suchardt LT: A simplified liquid culture medium for the growth of Haemophilus pertussis. J Bacteriol. 1949, 58: 127-134.PubMed CentralGoogle Scholar
- NIH: Revised guide for the care and use of laboratory animals. NIH GUIDE. 1996, 25:Google Scholar
- Barends M, van Oosten M, De Rond CG, Dormans JA, Osterhaus AD, Neijens HJ, Kimman TG: Timing of infection and prior immunization with respiratory syncytial virus (RSV) in RSV-enhanced allergic inflammation. J Infect Dis. 2004, 189: 1866-1872. 10.1086/386341.View ArticlePubMedGoogle Scholar
- Vos JG, de Klerk A, Krajnc EI, Kruizinga W, van Ommen B, Rozing J: Toxicity of bis(tri-n-butyltin)oxide in the rat. II. Suppression of thymus-dependent immune responses and of parameters of nonspecific resistance after short-term exposure. Toxicol Appl Pharmacol. 1984, 75: 387-408. 10.1016/0041-008X(84)90177-7.View ArticlePubMedGoogle Scholar
- Team RDC: R: A language and enviroment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2005, ISBN 3-900051-07-0: [http://www.r-project.org]Google Scholar
- Smyth GK, Speed T: Normalization of cDNA microarray data. Methods. 2003, 31: 265-273. 10.1016/S1046-2023(03)00155-5.View ArticlePubMedGoogle Scholar
- Pennings JLA, Heisterkamp SH: Normal probability plots for microarray experiments. Proceedings of the 12th International Conference on Intelligent Systems for Molecular Biology and the 3rd European Conference on Computational Biology 143. 2004Google Scholar
- Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: a Practical and Powerful Approach to Multiple Testing. J R Statist Soc B. 1995, 57: 288-300.Google Scholar
- GO: Gene Ontology. 2007, [http://www.geneontology.org]Google Scholar
- Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: 3-10.1186/gb-2003-4-5-p3.View ArticleGoogle Scholar
- Hosack DA, Dennis G, Sherman BT, Lane HC, Lempicki RA: Identifying biological themes within lists of genes with EASE. Genome Biol. 2003, 4: R70-10.1186/gb-2003-4-10-r70.PubMed CentralView ArticlePubMedGoogle Scholar
- Solin ML, Kaartinen M: Allelic polymorphism of mouse Igh-J locus, which encodes immunoglobulin heavy chain joining (JH) segments. Immunogenetics. 1992, 36: 306-313. 10.1007/BF00215659.View ArticlePubMedGoogle Scholar
- Martin RM, Silva A, Lew AM: The Igh-1 sequence of the non-obese diabetic (NOD) mouse assigns it to the IgG2c isotype. Immunogenetics. 1997, 46: 167-168. 10.1007/s002510050258.View ArticlePubMedGoogle Scholar
- Martin RM, Brady JL, Lew AM: The need for IgG2c specific antiserum when isotyping antibodies from C57BL/6 and NOD mice. J Immunol Methods. 1998, 212: 187-192. 10.1016/S0022-1759(98)00015-5.View ArticlePubMedGoogle Scholar
- Foster CS, Opremcak EM, Rice B, Wells P, Chung H, Thompson P, Fong LP, Raizman M: Clinical, pathologic, and immunopathologic characteristics of experimental murine herpes simplex virus stromal keratitis and uveitis is controlled by gene products from the Igh-1 locus on chromosome 12. Trans Am Ophthalmol Soc. 1987, 85: 293-311.PubMed CentralPubMedGoogle Scholar
- Opremcak EM, Wells PA, Thompson P, Daigle JA, Rice BA, Millin JA, Foster CS: Immunogenetic influence of Igh-1 phenotype on experimental herpes simplex virus type-1 corneal infection. Invest Ophthalmol Vis Sci. 1988, 29: 749-754.PubMedGoogle Scholar
- Tamesis RR, Foster CS: Natural killer cellular cytotoxicity against herpes simplex virus-infected cells in Igh-1-disparate mice. Invest Ophthalmol Vis Sci. 1990, 31: 2224-2229.PubMedGoogle Scholar
- Arrunategui-Correa V, Baltatzis S, Foster CS: The role of cytokines in experimental herpes simplex keratitis. Acta Virol. 1999, 43: 325-329.PubMedGoogle Scholar
- Lake JP, Kapp JA, Pierce CW: Characterization of L-glutamic acid60-L-alanine30-L-tyrosine10-specific suppressor T cells in responder mice restricted by Igh-C-linked genes. J Immunol. 1988, 140: 3296-3302.PubMedGoogle Scholar
- O'Hara RM, Sherr DH, Dorf ME: In vitro generation of suppressor T cells. Induction of CD3+, IgH-restricted suppressor cells. J Immunol. 1988, 141: 2935-2942.PubMedGoogle Scholar
- ArrayExpress. 2007, [http://www.ebi.ac.uk/arrayexpress/]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.