Genome-wide association study of antibody response to Newcastle disease virus in chicken
- Chenglong Luo†1, 3,
- Hao Qu†1, 3,
- Jie Ma1, 3,
- Jie Wang1, 3,
- Chunyu Li1, 3,
- Chunfen Yang1, 3,
- Xiaoxiang Hu2,
- Ning Li2 and
- Dingming Shu1, 3Email author
© Luo et al.; licensee BioMed Central Ltd. 2013
Received: 1 September 2012
Accepted: 6 May 2013
Published: 10 May 2013
Since the first outbreak in Indonesia in 1926, Newcastle disease has become one of the most common and contagious bird diseases throughout the world. To date, enhancing host antibody response by vaccination remains the most efficient strategy to control outbreaks of Newcastle disease. Antibody response plays an important role in host resistance to Newcastle disease, and selection for antibody response can effectively improve disease resistance in chickens. However, the molecular basis of the variation in antibody response to Newcastle disease virus (NDV) is not clear. The aim of this study was to detect genes modulating antibody response to NDV by a genome-wide association study (GWAS) in chickens.
To identify genes or chromosomal regions associated with antibody response to NDV after immunization, a GWAS was performed using 39,833 SNP markers in a chicken F2 resource population derived from a cross between two broiler lines that differed in their resistance. Two SNP effects reached 5% Bonferroni genome-wide significance (P<1.26×10-6). These two SNPs, rs15354805 and rs15355555, were both on chicken (Gallus gallus) chromosome 1 and spanned approximately 600 Kb, from 100.4 Mb to 101.0 Mb. Rs15354805 is in intron 7 of the chicken Roundabout, axon guidance receptor, homolog 2 (ROBO2) gene, and rs15355555 is located about 243 Kb upstream of ROBO2. Rs15354805 explained 5% of the phenotypic variation in antibody response to NDV, post immunization, in chickens. Rs15355555 had a similar effect as rs15354805 because of its linkage disequilibrium with rs15354805 (r2=0.98).
The region at about 100 Mb from the proximal end of chicken chromosome 1, including the ROBO1 and ROBO2 genes, has a strong effect on the antibody response to the NDV in chickens. This study paves the way for further research on the host immune response to NDV.
KeywordsChicken Newcastle disease Antibody response Genome-wide association study
Newcastle disease is a highly contagious viral disease of birds affecting over 250 domestic and wild avian species. Since Newcastle disease was first discovered in Indonesia in 1926, it has spread throughout the world . Newcastle disease is one of the most serious problems affecting poultry industries in many countries because its contagiousness and devastating effects [2–5]. Newcastle disease virus (NDV) has been classified into velogenic, mesogenic and lentogenic strains on the basis of their pathogenesis and virulence. Chickens with Newcastle disease have severe neurological and respiratory signs and show decreased egg quality and production [1, 6, 7]. There is no effective treatment for Newcastle disease; however, the use of prophylactic vaccines and maintenance of strict biosecurity measures can reduce the likelihood of outbreaks. Thus, the ability of chickens to mount an antibody response to NDV plays a key role in controlling Newcastle disease outbreaks, why understanding the molecular basis of immune response to NDV is important for the control of avian Newcastle disease. The antibody response to the same virus differs between chicken breeds , and selection for an antibody response may improve disease resistance in chickens . However, there have been few studies on controlling Newcastle disease outbreaks from the perspective of host resistance. Ten quantitative trait loci (QTLs) have been reported to be linked to the antibody response to NDV on chicken (Gallus gallus) chromosomes (GGA) 2, 3, 4, 5, 9, 13, 16, 18, 22 and Z [10, 11], but only a few causative genes have been identified because of low map resolution. The chicken 60K single nucleotide polymorphism (SNP) chip provides an opportunity to increase precision of QTL mapping by increasing genome coverage and map resolution . Using this chicken 60K SNP chip, some QTLs and candidate genes affecting chicken growth have been identified at 172–175 Mb on GGA1  and 71.6–80.2 Mb on GGA4 .
The aim of this study was to identify QTLs for the antibody response to NDV in chickens. A genome-wide association study (GWAS) was performed using the chicken 60K SNP chip, based on the data from an F2 resource population that was derived from the cross between two broiler lines with differing resistances.
Characterization of the antibody response to NDV in the F2population
The antibody response to NDV was positive in all 511 F2 birds at 41 days after the second immunization, as assessed by an enzyme linked immunosorbent assay (ELISA) test, and were highly variable. The mean ± standard error of the adjusted antibody responses to NDV was 3.38 ± 0.05, its coefficient of variance was 34.0%, and its distribution fitted a normal distribution (P=0.5527, Additional file 1: Figure S1). The heritability of the secondary antibody response to NDV was estimated as 0.24 ± 0.08.
Genome-wide association analysis
Genome-wide significant* SNPs for antibody response to the Newcastle disease virus
243 Kb upstream of ROBO2
Evaluation of the candidate region comprising rs15354805 and rs15355555
Bioinformatic analysis of candidate regions
As shown in Additional file 3: Table S2, there are 194 genes or non-coding RNAs in the candidate regions of 1-Mb windows (SNP position ± 0.5 Mb) surrounding each selected SNP (P<1.00×10-4). These genes participate in at least 25 KEGG pathways, and comprise at least 108 protein domains in the Interpro database of chickens. Most of the genes in the significantly (P<0.1) enriched pathways and protein domains are related to the immune system (Additional file 5: Table S3). Among these, CD2 and CD80, which contain nearly significant (P<1.00×10-4) SNPs (rs13904539, rs13904683 and GGaluGA032608) for the antibody response to NDV, belong to the cell adhesion molecules pathway (gga04514), which involves the major histocompatibility complex (MHC) and T cell receptor genes. LOC418424, junction adhesion molecule 2, CD200, CD80, ROBO1 and ROBO2 on GGA1, and Aggrecan on GGA10 all possess a conserved protein domain called the immunoglobulin-like fold (IPR013783, Additional file 5: Table S3), suggesting that these proteins could be members of the complex gene network that modulates the immune response. The Cbl proto-oncogene, E3 ubiquitin protein ligase B, interleukin 10 receptor, beta, interferon gamma receptor 1 and interferon gamma receptor 2 genes, whose proteins participate in the Jak-STAT signaling pathway (gga04630), are also implicated in the antibody response to NDV.
The host immune response to viruses is a complex process. The antibody response to NDV can be considered to be a quantitative trait under polygenic control, but with some QTLs [10, 11]. In the current study, GWAS identified a region located approximately 100 Mb from the proximal end of GGA1, in which rs15354805 was most significantly (P<1.26×10-6) associated with, and accounted for 5% of the phenotypic variation of, the antibody response to NDV (Table 1). Moderately high heritability (0.41 ± 0.17) of the chicken antibody response to NDV was observed three decades ago . Heritabilities of the primary antibody response to NDV two weeks post-vaccination were 0.27 ± 0.06 and 0.29 ± 0.05 in two Tanzania chicken ecotypes (Kuchi and Tanzania Medium), respectively . Lwelamira also reported recently that the heritability of the chicken primary antibody response to NDV at two weeks post-immunization was 0.22 ± 0.08 in the offspring of Kuchi chickens , a figure comparable to the 0.24 ± 0.08 estimated here for the secondary antibody response. The QTL linked to rs15354805 therefore accounted for >20% (0.05/0.24) of the additive genetic variance of this trait. It is reasonable to consider that this QTL would be important for improving the ability of the immune response to NDV in chickens by use of marker-assisted selection. In addition, the effect of combining rs15354805 on GGA1 and rs14008095 on GGA10 is almost same as the sum of the effects of all SNPs with P values less than 1.00×10-4 for the antibody response to NDV (7% vs. 8% of the phenotypic variation) in this study. On that basis, there are probably only two divergent QTLs for this trait on GGA1 and GGA10, while other SNP effects might result from LD.
It should be noted that the region found in the present study to be significantly (P<1.26×10-6) associated with the antibody response to NDV differs from those found in earlier reports [10, 11]. Yonash et al. reported that QTLs for the antibody response to NDV were located on GGA2 and GGA18 , whereas Biscarini et al. found 13 QTLs associated with the antibody response to NDV on GGA3, GGA4, GGA5, GGA9, GGA13, GGA16, GGA22 and GGAZ . There are many possible reasons for the differing results among these studies. The population specificity of QTLs for the antibody response to NDV is probably the most important reason . The current study was based on a broiler resource population, while the other two studies were based on nine different laying lines , and a chicken resource population built from the intercross of broilers with high and low antibody response to Escherichia coli vaccine . In addition, Yonash et al. focused on the primary antibody response to NDV , while the current study examined the secondary antibody response to NDV. In the primary antibody response, the predominant class of antibody produced is immunoglobulin M, and in the secondary antibody response it is immunoglobulin Y; therefore, the ranking of birds on antibody response to the primary immunization may differ from that after the secondary boosting resulting in different QTLs being detected. In fact, the QTLs for the antibody response to NDV in this study may reflect the ability of the memory cell pool to respond to NDV. Interestingly, we observed that certain immunity-related genes, such as CD2, CD80 and CD200, on the detected QTLs (Additional file 3: Table S2). However, some detected QTLs in the current study were similar to those in two other studies, if the significance threshold set as P<0.05 (Figure 1). For example, rs14520447 (GGA5: 21181095 bp, P=0.0031), rs14520528 (GGA5: 21329838 bp, P=0.0194) and GGAluGA278144 (GGA5: 21355363 bp, P=0.0194) were to the strongest of the QTLs detected by Biscarini et al. , and rs10727400 (GGA18: 4572183 bp, P=0.0336), rs15821136 (GGA18: 4587942 bp, P=0.0171) and GGAluGA120237 (GGA18: 4594180 bp, P=0.0200) were similar to the QTLs detected on GGA18 in the report of Yonash et al. .
The MHC plays an important role in animal immunity [19–22]. However, the present study found that the most reliable QTL linked to antibody response to NDV could be in the region that included rs15354805 and rs15355555 on GGA1, whereas it did not detect any significant SNPs on GGA16, which harbors the chicken MHC. This may be because very few SNP markers on GGA16 could be used for QTL mapping in the current study (Additional file 2: Table S1). Further studies are necessary in order to detect QTLs for the antibody response to NDV on GGA16.
The 2 Mb region of 99–101 Mb region containing three coding genes on GGA1 had the most significant effects for the secondary antibody response to NDV (Figure 2). Within this region, LOC769250 is a predicted gene with no clear function, and ROBO1 and ROBO2 belong to the ROBO family, which are members of the immunoglobulin superfamily . The encoded products of both genes are integral membrane protein receptors for the SLIT-family, and have similar functions in axon guidance and migration of neuronal precursor cells [24–27]. Many molecules, such as Semaphorin 7A, Neurokinin A and Plexin-B2, participate in nervous system and immune system processes [28–32]. Perhaps ROBO1 and ROBO2 also affect the immune system. In fact, SLIT/ROBO signaling can influence the expression of Rac1 and CDC42 . Rac1 plays a crucial role in T cell development  while CDC42 is essential for memory T cell growth and differentiation [35, 36]. Therefore, SLIT/ROBO signaling could contribute to adaptive immunity through Rac1 or/and CDC42. Consistently, the focus of QTL mapping in this study was the difference of secondary antibody response to NDV that would be mainly related to activation, growth and differentiation of the memory T cells after the second immunization. In addition, SLIT/ROBO signaling plays important roles in leukocyte chemotaxis [37–39]. Therefore, ROBO1 and ROBO2 could well be candidate genes for modulating the antibody response to NDV. However, determining the functional role of these QTLs in the antibody response to NDV needs requires further studies. Studies of chicken network pathways involving SLIT/ROBO signaling need to be replicated to validate the present findings.
In addition to ROBO1 and ROBO2, some of the other genes that are close to the QTL region (P<1.00×10-4, Additional file 3: Table S2) in the current study might modulate antibody responses to NDV, based on their biological function and networks. These genes could comprise a complex network (Additional file 6: Figure S3). On GGA10, interleukin 16, whose encoded protein is linked to ROBO2 by three pathways (Additional file 6: Figure S3), participates in leukocyte chemotaxis and positive regulation of responses to external stimuli , while B-cell CLL/lymphoma 2 - related protein A1 is essential for lymphocyte activation and cell survival . In fact, seven genes in the QTL region on GGA1 were enriched in the cell adhesion molecules pathway, which affects the immune system (Additional file 5: Table S3). Among these genes, CD2 and CD80 were initially considered to be candidate genes. CD2 not only plays an important role in the generation of nanotubes between natural killer cells and their targets , but also provides signals to control the development and differentiation of T cells [43–45]. CD80 is expressed on B cells, but also induces T cell proliferation and cytokine production [46–48]. In addition, CD200, which is about 180 Kb from rs14856616 (P=6.08×10-6), is expressed on the surface of immune cells and affects the adaptive immune system by stimulating proliferation of T cells . Interestingly, evidence shows that dendritic cells and T cells show a stronger response in the absence of CD200 inhibitory signaling . In addition, there is a relationship between CD200 and influenza virus infection . SAM domain, SH3 domain and nuclear localization signals 1 is another immune-related protein that participates in B cell activation and differentiation . CD247, IFNGR2 and IFNAR1 belong to the natural killer cell mediated cytotoxicity pathway (Additional file 5: Table S3). Obviously, a large number of genes related to innate and adaptive immunity are in the 89.7-111.0 Mb region of GGA1. This QTL region also overlaps the QTL for the antibody response to sheep red blood cells . Thus this genomic region is of great importance for the chicken immune response and, possibly, disease resistance, and has been impacted by intense artificial selection (e.g., ZHp=−6.2 in ROBO2) during chicken domestication .
This study has exposed a novel QTL region associated with the antibody response to NDV. The QTL region was located approximately 100 Mb from the proximal end of GGA1. This region may play an important role in the chicken immune response. ROBO1 and ROBO2 were found to be important prospective candidate genes for modulating the antibody response to NDV in chickens; however, the roles of ROBO1 and ROBO2 require further study, such as exploring the effects of ROBO1 and/or ROBO2 silencing and overexpression in vitro and in vivo on immune gene networks.
The Animal Care Committee of the Institute of Animal Science, Guangdong Academy of Agricultural Sciences (Guangzhou, P. R. China) approved this study (approval number GAAS-IAS-2009-73). Animals involved in the current study were humanely euthanized, as necessary, to reduce any suffering.
The F2 resource population was built from the full-sib intercross of two divergent lines (23 P and 51 F1). The first line was a fast-growing Chinese yellow broiler raised at Guangdong Wiz Agricultural Science & Technology Co. Ltd., which had undergone more than ten generations of selection for high growth rate and meat quality tailored to Chinese tastes. The second line was the Huiyang beard chicken, a Chinese local breed with a low growth rate and high meat quality tailored to Chinese tastes. The resistance of the second line is better than that of the first line (Data not shown). The F2 population comprised 511 individuals from six hatches. Before 40 days of age, each hatch was maintained in a group cage. From 40 days to 91 days, each chicken was reared in an individual cage with an individual trough. The house was equipped with water curtain systems and supplied with 24-hour light per day. A starter diet (200 g total protein and 2,900 kcal ME/kg) from hatch to 35 days of age and a grower diet (180 g total protein and 2,950 kcal ME/kg) from 43 days to 91 days of age were provided ad libitum.
The chickens were immunized with a commercial NDV live vaccine of the LaSota strain (Intervet International B.V., Boxmeer, Netherlands), using the standard dose given in the instructions of the vaccine, by eye drop at 25 days of age and at 50 days of age, respectively. At 41 days after the second immunization (91 days of age), serum samples were collected from the chickens. The antibody response to NDV was determined using an indirect ELISA and expressed as the S/P value of corresponding dilutions, according to the instruction of the commercial ELISA kit for testing the antibody response to NDV (BioCheck, Inc., Foster City, CA, USA). The ELISA kit measures NDV-specific immunoglobulin, Y but not immunoglobulin A nor immunoglobulin M.
Genomic DNA extraction from venous blood was performed using a phenol/chloroform method. For all 511 F2 individuals, the quality and concentration of genomic DNA fulfilled the requirements of the Illumina Infinium SNP genotyping platform. Genotyping using the Illumina 60K Chicken SNP Beadchip  was carried out at the Illumina-certified service provider, DNA LandMarks Inc., Saint-Jean-sur-Richelieu, Quebec, Canada. Quality control was assessed in GenomeStudio v2008.1 . Six samples were excluded because they had more than 5% missing SNP genotypes. The final SNP set used in this study consisted of 39,833 SNP markers, meeting the following selection criteria: call frequency >95%, heterozygosity cluster intensity and separation value >0.4, and minor allele frequency >0.1. Information about the SNPs on each chromosome is summarized in Additional file 2: Table S1. Gene locations were based on the Ensembl Genome Browser  and the National Center for Biotechnology Information (NCBI) .
where A is the additive relationship matrix between individuals in the pedigree of all birds; is animal additive genetic variance; I is an identity matrix; and is the residual variance. Heritability of antibody response to NDV was estimated as , and , thus, the ratio of the QTL effect compared with the additive genetic variance was calculated as , where , and R2 are the phenotypic variance, the QTL effect variance, and coefficient of correction for the regression equation with one or more QTL as independent variable and the antibody response to NDV as a dependent variable, which indicated the ratio of the QTL effect compared with the phenotypic variance of antibody response to NDV, respectively.
where Yijkl is the phenotypic value for the antibody response to NDV, μ is the overall mean, SEXi is the effect of the ith sex, Hj is the effect of the jth hatch, Sk is the effect of the kth sire, Dkl is the effect of the lth dam within kth sire, and eijkl is the residual effect. The adjusted phenotypic value was calculated as μ + eijkl.
Second, the adjusted antibody response to NDV, which followed a normal distribution (Additional file 1: Figure S1), was used to carry out genome-wide association analysis in PLINK using linear regression analyses without any covariates other than SNPs . The threshold P value of the 5% Bonferroni genome-wide significance was 1.26×10-6 (0.05/39833), and the threshold P value for the “suggestive linkage” significance, which allowed one false positive effect in a genome-wide test, was 2.51×10-6 (1/39833). Using the “gap” package in R v2.12.0 (http://www.r-project.org), a Manhattan plot was drawn to describe the P values of all SNPs. This study calculated the fraction of the phenotypic variance explained by the significantly associated SNPs according to a previously described method . Haploview software  was used to calculate LD as an r-square value and build haplotype blocks under the four gamete rule (the 4th gamete must be observed at frequency >0.01) in the candidate region associated with antibody response to NDV in the F2 population.
Names of genes in the candidate regions with 1-Mb windows (SNP position ± 0.5 Mb) surrounding each significant SNP (P<1.00×10-4) were obtained from the Ensembl Genome Browser  and the NCBI . A pathway and protein domain analysis was performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID) [62, 63]. A network of gene correlation was constructed by integrating the pathways from the human database of the Kyoto Encyclopedia of Genes and Genomes (KEGG)  (http://www.genome.jp/kegg/), the Gene Map Annotator and Pathway Profiler (GenMAPP)  (genmapp.org/) and the BioCarta (http://www.biocarta.com/genes/index.asp).
The authors would like to thank Ms Zheya Sheng (College of Biological Science, China Agricultural University, Beijing, China) for her technical assistance in SNP genotyping data handling. Dr W. Bruce Currie, Emeritus Professor, Cornell University offered suggestions on presentation. This study was supported by three grants. The first was the Natural Science Foundation of Guangdong Province (Grant No. S2011010001471). The second was from the Guangdong Key Scientific and Technological Projects (Grant No. 2011A020201007). The third was from the National Natural Science Foundation of China (Grant No. U0831003).
- Alexander DJ: Newcastle disease. Br Poult Sci. 2001, 42 (1): 5-22. 10.1080/713655022.View ArticlePubMedGoogle Scholar
- Ababneh MM, Dalab AE, Alsaad SR: Al-Zghoul MB. 2012, Al-Natour MQ: Molecular characterization of a recent Newcastle disease virus outbreak in Jordan. Res Vet SciGoogle Scholar
- Aldous EW, Manvell RJ, Cox WJ, Ceeraz V, Harwood DG, Shell W, Alexander DJ, Brown IH: Outbreak of Newcastle disease in pheasants (Phasianus colchicus) in south-east England in July 2005. Vet Rec. 2007, 160 (14): 482-484. 10.1136/vr.160.14.482.View ArticlePubMedGoogle Scholar
- Irvine RM, Aldous EW, Manvell RJ, Cox WJ, Ceeraz V, Fuller CM, Wood AM, Milne JC, Wilson M, Hepple RG: Outbreak of Newcastle disease due to pigeon paramyxovirus type 1 in grey partridges (Perdix perdix) in Scotland in October 2006. Vet Rec. 2009, 165 (18): 531-535. 10.1136/vr.165.18.531.View ArticlePubMedGoogle Scholar
- Dortmans JC, Peeters BP, Koch G: Newcastle disease virus outbreaks: vaccine mismatch or inadequate application?. Vet Microbiol. 2012, 160 (1–2): 17-22.View ArticlePubMedGoogle Scholar
- Miller PJ, Decanini EL, Afonso CL: Newcastle disease: evolution of genotypes and the related diagnostic challenges. Infect Genet Evol. 2010, 10 (1): 26-35. 10.1016/j.meegid.2009.09.012.View ArticlePubMedGoogle Scholar
- Villegas P: Viral diseases of the respiratory system. Poult Sci. 1998, 77 (8): 1143-1145.View ArticlePubMedGoogle Scholar
- Pitcovski J, Cahaner A, Heller ED, Zouri T, Gutter B, Gotfried Y, Leitner G: Immune response and resistance to infectious bursal disease virus of chicken lines selected for high or low antibody response to Escherichia coli. Poult Sci. 2001, 80 (7): 879-884.View ArticlePubMedGoogle Scholar
- Gross WG, Siegel PB, Hall RW, Domermuth CH, DuBoise RT: Production and persistence of antibodies in chickens to sheep erythrocytes. 2. Resistance to infectious diseases. Poult Sci. 1980, 59 (2): 205-210. 10.3382/ps.0590205.View ArticlePubMedGoogle Scholar
- Biscarini F, Bovenhuis H, van Arendonk JA, Parmentier HK, Jungerius AP, van der Poel JJ: Across-line SNP association study of innate and adaptive immune response in laying hens. Anim Genet. 2010, 41 (1): 26-38. 10.1111/j.1365-2052.2009.01960.x.View ArticlePubMedGoogle Scholar
- Yonash N, Cheng HH, Hillel J, Heller DE, Cahaner A: DNA microsatellites linked to quantitative trait loci affecting antibody response and survival rate in meat-type chickens. Poult Sci. 2001, 80 (1): 22-28.View ArticlePubMedGoogle Scholar
- Groenen MA, Megens HJ, Zare Y, Warren WC, Hillier LW, Crooijmans RP, Vereijken A, Okimoto R, Muir W, Cheng HH: The development and characterization of a 60K SNP chip for chicken. BMC Genomics. 2011, 12 (1): 274-10.1186/1471-2164-12-274.PubMed CentralView ArticlePubMedGoogle Scholar
- Xie L, Luo C, Zhang C, Zhang R, Tang J, Nie Q, Ma L, Hu X, Li N, Da Y: Genome-Wide Association Study Identified a Narrow Chromosome 1 Region Associated with Chicken Growth Traits. PLoS One. 2012, 7 (2): e30910-10.1371/journal.pone.0030910.PubMed CentralView ArticlePubMedGoogle Scholar
- Gu X, Feng C, Ma L, Song C, Wang Y, Da Y, Li H, Chen K, Ye S, Ge C: Genome-wide association study of body weight in chicken F2 resource population. PLoS One. 2011, 6 (7): e21872-10.1371/journal.pone.0021872.PubMed CentralView ArticlePubMedGoogle Scholar
- Soller M, Heller D, Peleg B, Ron-Kuper N, Hornstein K: Genetic and phenotypic correlations between immune response to Escherichia coli and to Newcastle disease virus vaccines. Poult Sci. 1981, 60 (1): 49-53. 10.3382/ps.0600049.View ArticlePubMedGoogle Scholar
- Lwelamira J, Kifaro GC, Gwakisa PS: Genetic parameters for body weights, egg traits and antibody response against Newcastle Disease Virus (NDV) vaccine among two Tanzania chicken ecotypes. Trop Anim Health Prod. 2009, 41 (1): 51-59. 10.1007/s11250-008-9153-2.View ArticlePubMedGoogle Scholar
- Lwelamira J: Phenotypic and genetic parameters for body weights and antibody response against Newcastle disease virus (NDV) vaccine for Kuchi chicken ecotype of Tanzania under extensive management. Trop Anim Health Prod. 2012, 44 (7): 1529-1534. 10.1007/s11250-012-0099-z.View ArticlePubMedGoogle Scholar
- Ioannidis JP, Thomas G, Daly MJ: Validating, augmenting and refining genome-wide association signals. Nat Rev Genet. 2009, 10 (5): 318-329.View ArticlePubMedGoogle Scholar
- Chapman SJ, Hill AV: Human genetic susceptibility to infectious disease. Nat Rev Genet. 2012, 13 (3): 175-188.PubMedGoogle Scholar
- Fruci D, Benevolo M, Cifaldi L, Lorenzi S, Lo ME, Tremante E, Giacomini P: Major histocompatibility complex class i and tumour immuno-evasion: how to fool T cells and natural killer cells at one time. Curr Oncol. 2012, 19 (1): 39-41.PubMed CentralView ArticlePubMedGoogle Scholar
- Westerdahl H, Asghar M, Hasselquist D, Bensch S: Quantitative disease resistance: to better understand parasite-mediated selection on major histocompatibility complex. Proc Biol Sci. 2012, 279 (1728): 577-584. 10.1098/rspb.2011.0917.PubMed CentralView ArticlePubMedGoogle Scholar
- Sumners LH, Cox CM, Kim S, Salevsky JE, Siegel PB, Dalloul RA: Immunological responses to Clostridium perfringens alpha-toxin in two genetically divergent lines of chickens as influenced by major histocompatibility complex genotype. Poult Sci. 2012, 91 (3): 592-603. 10.3382/ps.2011-01775.View ArticlePubMedGoogle Scholar
- National Center for Biotechnology Information: 2011, Entrez gene, Available: http://www.ncbi.nlm.nih.gov
- Long H, Sabatier C, Ma L, Plump A, Yuan W, Ornitz DM, Tamada A, Murakami F, Goodman CS, Tessier-Lavigne M: Conserved roles for Slit and Robo proteins in midline commissural axon guidance. Neuron. 2004, 42 (2): 213-223. 10.1016/S0896-6273(04)00179-5.View ArticlePubMedGoogle Scholar
- Jia L, Cheng L, Raper J: Slit/Robo signaling is necessary to confine early neural crest cells to the ventral migratory pathway in the trunk. Dev Biol. 2005, 282 (2): 411-421. 10.1016/j.ydbio.2005.03.021.View ArticlePubMedGoogle Scholar
- Dugan JP, Stratton A, Riley HP, Farmer WT, Mastick GS: Midbrain dopaminergic axons are guided longitudinally through the diencephalon by Slit/Robo signals. Mol Cell Neurosci. 2011, 46 (1): 347-356. 10.1016/j.mcn.2010.11.003.PubMed CentralView ArticlePubMedGoogle Scholar
- Mastick GS, Farmer WT, Altick AL, Nural HF, Dugan JP, Kidd T, Charron F: Longitudinal axons are guided by Slit/Robo signals from the floor plate. Cell Adh Migr. 2010, 4 (3): 337-341. 10.4161/cam.4.3.11219.PubMed CentralView ArticlePubMedGoogle Scholar
- Pasterkamp RJ, Peschon JJ, Spriggs MK, Kolodkin AL: Semaphorin 7A promotes axon outgrowth through integrins and MAPKs. Nature. 2003, 424 (6947): 398-405. 10.1038/nature01790.View ArticlePubMedGoogle Scholar
- Helmuth L: NEUROSCIENCE: Immune Molecules Prune Neural Links. Science. 2000, 290 (5499): 2051a-10.1126/science.290.5499.2051a.View ArticleGoogle Scholar
- Ransohoff RM, Brown MA: Innate immunity in the central nervous system. J Clin Invest. 2012, 122 (4): 1164-1171. 10.1172/JCI58644.PubMed CentralView ArticlePubMedGoogle Scholar
- Kitamura H, Kobayashi M, Wakita D, Nishimura T: Neuropeptide signaling activates dendritic cell-mediated type 1 immune responses through neurokinin-2 receptor. J Immunol. 2012, 188 (9): 4200-4208. 10.4049/jimmunol.1102521.View ArticlePubMedGoogle Scholar
- Roney KE, O'Connor BP, Wen H, Holl EK, Guthrie EH, Davis BK, Jones SW, Jha S, Sharek L, Garcia-Mata R: Plexin-B2 negatively regulates macrophage motility, Rac, and Cdc42 activation. PLoS One. 2011, 6 (9): e24795-10.1371/journal.pone.0024795.PubMed CentralView ArticlePubMedGoogle Scholar
- Ypsilanti AR, Zagar Y, Chedotal A: Moving away from the midline: new developments for Slit and Robo. Development. 2010, 137 (12): 1939-1952. 10.1242/dev.044511.View ArticlePubMedGoogle Scholar
- Guo F, Cancelas JA, Hildeman D, Williams DA, Zheng Y: Rac GTPase isoforms Rac1 and Rac2 play a redundant and crucial role in T-cell development. Blood. 2008, 112 (5): 1767-1775. 10.1182/blood-2008-01-132068.PubMed CentralView ArticlePubMedGoogle Scholar
- Guo F, Zhang S, Tripathi P, Mattner J, Phelan J, Sproles A, Mo J, Wills-Karp M, Grimes HL, Hildeman D: Distinct roles of Cdc42 in thymopoiesis and effector and memory T cell differentiation. PLoS One. 2011, 6 (3): e18002-10.1371/journal.pone.0018002.PubMed CentralView ArticlePubMedGoogle Scholar
- Smits K, Iannucci V, Stove V, Van Hauwe P, Naessens E, Meuwissen PJ, Arien KK, Bentahir M, Plum J, Verhasselt B: Rho GTPase Cdc42 is essential for human T-cell development. Haematologica. 2010, 95 (3): 367-375. 10.3324/haematol.2009.006890.PubMed CentralView ArticlePubMedGoogle Scholar
- Fernandis AZ, Ganju RK: Slit: a roadblock for chemotaxis. Sci STKE. 2001, 2001 (91): e1-Google Scholar
- Wu JY, Feng L, Park HT, Havlioglu N, Wen L, Tang H, Bacon KB, Jiang Z, Zhang X, Rao Y: The neuronal repellent Slit inhibits leukocyte chemotaxis induced by chemotactic factors. Nature. 2001, 410 (6831): 948-952. 10.1038/35073616.PubMed CentralView ArticlePubMedGoogle Scholar
- Tole S, Mukovozov IM, Huang YW, Magalhaes MA, Yan M, Crow MR, Liu GY, Sun CX, Durocher Y, Glogauer M: The axonal repellent, Slit2, inhibits directional migration of circulating neutrophils. J Leukoc Biol. 2009, 86 (6): 1403-1415. 10.1189/jlb.0609391.View ArticlePubMedGoogle Scholar
- Nimmanapalli R, Sharmila C, Reddy PG: Immunomodulation of caprine lentiviral infection by interleukin-16. Comp Immunol Microbiol Infect Dis. 2010, 33 (6): 529-536. 10.1016/j.cimid.2009.09.003.View ArticlePubMedGoogle Scholar
- Rasooly R, Schuster GU, Gregg JP, Xiao JH, Chandraratna RA, Stephensen CB: Retinoid x receptor agonists increase bcl2a1 expression and decrease apoptosis of naive T lymphocytes. J Immunol. 2005, 175 (12): 7916-7929.View ArticlePubMedGoogle Scholar
- Comerci CJ, Mace EM, Banerjee PP, Orange JS: CD2 Promotes Human Natural Killer Cell Membrane Nanotube Formation. PLoS One. 2012, 7 (10): e47664-10.1371/journal.pone.0047664.PubMed CentralView ArticlePubMedGoogle Scholar
- Kaizuka Y, Douglass AD, Vardhana S, Dustin ML, Vale RD: The coreceptor CD2 uses plasma membrane microdomains to transduce signals in T cells. J Cell Biol. 2009, 185 (3): 521-534. 10.1083/jcb.200809136.PubMed CentralView ArticlePubMedGoogle Scholar
- Kalland ME, Oberprieler NG, Vang T, Tasken K, Torgersen KM: T cell-signaling network analysis reveals distinct differences between CD28 and CD2 costimulation responses in various subsets and in the MAPK pathway between resting and activated regulatory T cells. J Immunol. 2011, 187 (10): 5233-5245. 10.4049/jimmunol.1101804.View ArticlePubMedGoogle Scholar
- Meinl E, Lengenfelder D, Blank N, Pirzer R, Barata L, Hivroz C: Differential requirement of ZAP-70 for CD2-mediated activation pathways of mature human T cells. J Immunol. 2000, 165 (7): 3578-3583.View ArticlePubMedGoogle Scholar
- Good-Jacobson KL, Song E, Anderson S, Sharpe AH, Shlomchik MJ: CD80 expression on B cells regulates murine T follicular helper development, germinal center B cell survival, and plasma cell generation. J Immunol. 2012, 188 (9): 4217-4225. 10.4049/jimmunol.1102885.PubMed CentralView ArticlePubMedGoogle Scholar
- Lim TS, Goh JK, Mortellaro A, Lim CT, Hammerling GJ, Ricciardi-Castagnoli P: CD80 and CD86 differentially regulate mechanical interactions of T-cells with antigen-presenting dendritic cells and B-cells. PLoS One. 2012, 7 (9): e45185-10.1371/journal.pone.0045185.PubMed CentralView ArticlePubMedGoogle Scholar
- Perez N, Karumuthil-Melethil S, Li R, Prabhakar BS, Holterman MJ, Vasu C: Preferential costimulation by CD80 results in IL-10-dependent TGF-beta1(+) -adaptive regulatory T cell generation. J Immunol. 2008, 180 (10): 6566-6576.PubMed CentralView ArticlePubMedGoogle Scholar
- Wright GJ, Jones M, Puklavec MJ, Brown MH, Barclay AN: The unusual distribution of the neuronal/lymphoid cell surface CD200 (OX2) glycoprotein is conserved in humans. Immunology. 2001, 102 (2): 173-179. 10.1046/j.1365-2567.2001.01163.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Holt PG, Strickland DH: The CD200-CD200R axis in local control of lung inflammation. Nat Immunol. 2008, 9 (9): 1011-1013. 10.1038/ni0908-1011.View ArticlePubMedGoogle Scholar
- Rygiel TP, Rijkers ES, de Ruiter T, Stolte EH, van der Valk M, Rimmelzwaan GF, Boon L, van Loon AM, Coenjaerts FE, Hoek RM: Lack of CD200 enhances pathological T cell responses during influenza infection. J Immunol. 2009, 183 (3): 1990-1996. 10.4049/jimmunol.0900252.View ArticlePubMedGoogle Scholar
- Zhu YX, Benn S, Li ZH, Wei E, Masih-Khan E, Trieu Y, Bali M, McGlade CJ, Claudio JO, Stewart AK: The SH3-SAM adaptor HACS1 is up-regulated in B cell activation signaling cascades. J Exp Med. 2004, 200 (6): 737-747. 10.1084/jem.20031816.PubMed CentralView ArticlePubMedGoogle Scholar
- Dorshorst BJ, Siegel PB, Ashwell CM: Genomic regions associated with antibody response to sheep red blood cells in the chicken. Anim Genet. 2011, 42 (3): 300-308. 10.1111/j.1365-2052.2010.02146.x.View ArticlePubMedGoogle Scholar
- Rubin CJ, Zody MC, Eriksson J, Meadows JR, Sherwood E, Webster MT, Jiang L, Ingman M, Sharpe T, Ka S: Whole-genome resequencing reveals loci under selection during chicken domestication. Nature. 2010, 464 (7288): 587-591. 10.1038/nature08832.View ArticlePubMedGoogle Scholar
- Illumina: Infinium® genotyping data analysis: 2010, http://www.illumina.com/Documents/products/technotes/technote_infinium_genotyping_data_analysis.pdf,
- Ensembl Genome Browser: Ensembl release 63. 2011, Available: http://useast.ensembl.org/index.htmlGoogle Scholar
- Gilmour AR, Thompson R, Cullis BR: Average Information REML: An Efficient Algorithm for Variance Parameter Estimation in Linear Mixed Models. Biometrics. 1995, 51 (4): 1440-1450. 10.2307/2533274.View ArticleGoogle Scholar
- Madsen P, SO Rensen P, Su G, Damgaard LH, Thomsen H, Labouriau R: DMU - a package for analyzing multivariate mixed models. 2006, Belo Horizonte, Minas Gerais, Brazil: Proc. 8th World Congr. Genet. Appl. Livest. Prod, 11-27.Google Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81 (3): 559-575. 10.1086/519795.PubMed CentralView ArticlePubMedGoogle Scholar
- Gudbjartsson DF, Walters GB, Thorleifsson G, Stefansson H, Halldorsson BV, Zusmanovich P, Sulem P, Thorlacius S, Gylfason A, Steinberg S: Many sequence variants affecting diversity of adult human height. Nat Genet. 2008, 40 (5): 609-615. 10.1038/ng.122.View ArticlePubMedGoogle Scholar
- Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.View ArticlePubMedGoogle Scholar
- Dennis GJ, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4 (5): P3-10.1186/gb-2003-4-5-p3.View ArticlePubMedGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4 (1): 44-57.View ArticleGoogle Scholar
- Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28 (1): 27-30. 10.1093/nar/28.1.27.PubMed CentralView ArticlePubMedGoogle Scholar
- Dahlquist KD, Salomonis N, Vranizan K, Lawlor SC, Conklin BR: GenMAPP, a new tool for viewing and analyzing microarray data on biological pathways. Nat Genet. 2002, 31 (1): 19-20. 10.1038/ng0502-19.View ArticlePubMedGoogle Scholar
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.