Figures
Abstract
Background
Viral attachment and cell entry host factors are important for viral replication, pathogenesis, and the generation and sustenance of immune responses after infection and/or vaccination, and are plausible genetic regulators of vaccine-induced immunity.
Methods
Using a tag-SNP approach in candidate gene study, we assessed the role of selected cell surface receptor genes, attachment factor-related genes, along with other immune genes in the genetic control of immune response variations after live rubella vaccination in two independent study cohorts.
Citation: Haralambieva IH, Lambert ND, Ovsyannikova IG, Kennedy RB, Larrabee BR, Pankratz VS, et al. (2014) Associations between Single Nucleotide Polymorphisms in Cellular Viral Receptors and Attachment Factor-Related Genes and Humoral Immunity to Rubella Vaccination. PLoS ONE 9(6): e99997. https://doi.org/10.1371/journal.pone.0099997
Editor: Tjeerd Kimman, Wageningen University and Research Centre, Netherlands
Received: February 7, 2014; Accepted: May 21, 2014; Published: June 19, 2014
Copyright: © 2014 Haralambieva et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The research reported in this publication was supported by the National Institute of Allergy and Infectious Diseases of the National Institutes of Health under Award Number R37 AI048793-11 (which recently received a MERIT award). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: GP is the chair of a Safety Evaluation Committee for novel non-rubella investigational vaccine trials being conducted by Merck Research Laboratories. GP offers consultative advice on vaccine development to Merck & Co. Inc., CSL Biotherapies, Avianax, Sanofi Pasteur, Dynavax, Novartis Vaccines and Therapeutics, PAXVAX Inc, and Emergent Biosolutions. GP and IO hold two patents related to vaccinia peptide research: U.S. Patent “Vaccinia Virus Polypeptides,” U.S. Patent Number: 13/222,862 (issued 31 August 2011), and U.S. Patent “Naturally Processed Measles Virus Peptides from Class II HLA Molecules,” U.S. Patent Number US 2010/0119536 A1 (issued 13 May 2010). These activities have been reviewed by the Mayo Clinic Conflict of Interest Review Board and are conducted in compliance with Mayo Clinic Conflict of Interest policies. This research has been reviewed by the Mayo Clinic Conflict of Interest Review Board and was conducted in compliance with Mayo Clinic Conflict of Interest policies. This does not alter adherence to PLOS ONE policies on sharing data and materials.
Introduction
Host genetic determinants play an important role in the generation and regulation of humoral and cellular immune responses after vaccination and/or infection [1], [2]. Studying measles-mumps-rubella (MMR) vaccine immunogenicity in twins, we have previously demonstrated a high level of heritability, approximately 46%, of humoral immune response variance following rubella vaccination [3]. Multiple host genes operate at a gene and gene-network level to shape and regulate the quality, duration and magnitude of rubella vaccine-induced humoral and cellular immune responses [4], [5], [6], [7], [8], [9], [10].
Viral attachment, cell entry, fusion with the cell membrane, and viral genome translocation into target cells are essential early stages initiating the viral infectious cycle. These are key steps in viral replication and dissemination, as well as virus-host interactions involving the generation and maintenance of the immune response. Recent discoveries point to the role of newly discovered cellular receptors and attachment factors for several important human viruses, including rubella virus (Myelin Oligodendrocyte Glycoprotein, MOG), measles virus (measles virus epithelial cell receptor PVRL4, poliovirus receptor-related 4 gene, Nectin-4; and the transmembrane C-type lectin DC-SIGN/CD209), Rift Valley Fever virus/RVFV (DC-SIGN/CD209), poliovirus (Poliovirus receptor, PVR/CD155, Nectin-5) and herpesvirus (poliovirus receptor-related protein 2, PVRL2/CD112, Nectin-2) in infection, disease pathogenesis and immunity [11], [12], [13], [14], [15], [16], [17], [18], [19], [20]. Host cell surface receptors and attachment factors are plausible genetic regulators of measles vaccine-induced immunity [21], [22], [23], [24], but their influence on immunity after rubella vaccination has never been examined.
In the current candidate gene association study, we follow up on previously found genetic associations and also make use of recent major discoveries in the virology field, and explore the plausible role of selected cell surface receptor-, and attachment factor-related genes, such as MOG and poliovirus receptor-related gene family members, in the genetic control of immune response variations after live rubella vaccination. Our results from two independent study cohorts (discovery and replication) strongly suggest that genetic variants from these genes play a role in modulating humoral immune responses following rubella vaccination.
Methods
The methods described herein are similar or identical to those published for our previous studies [6], [7], [8], [9], [25], [26], [27].
Study Participants
The study cohort was a large population-based sample of 2,250 healthy children, older adolescents, and healthy adults (age 11 to 40 years), residing in Rochester, MN, and San Diego, CA, with clinical and demographic characteristics previously reported [7], [9], [21], [28], [29], [30].
The Rochester cohort comprised a sample of 1,145 individuals from three independent age-stratified random cohorts of healthy schoolchildren and young adults from all socioeconomic strata from Olmsted County, MN, enrolled between 2001 and 2009, as published elsewhere [7], [9], [21], [28], [29], [30]. Eleven hundred and one parents agreed to let their children join the current rubella vaccine study.
Between July 2005 and September 2006, we enrolled an additional 1,076 healthy older adolescents and healthy adults (age 18 to 40 years, San Diego cohort) from armed forces personnel at the Naval Health Research Center (NHRC) in San Diego, CA, as previously described [30]. All subjects included in the current rubella vaccine study had a documented receipt of measles-mumps-rubella (MMR) vaccine. The Institutional Review Boards of the Mayo Clinic and the NHRC approved the study, and written informed consent was obtained from each subject, from the parents of all children who participated in the study, as well as written assent from age-appropriate participants.
Rubella Virus-specific Neutralization Assay (sICNA)
We used a modified version of the immunocolorimetric-based neutralization method described by Chen et al. [31], optimized to a high-throughput micro-format, to measure rubella virus-specific neutralizing antibodies [25], [26]. Subjects’ sera were heat-inactivated for 1 hour at 56°C. Sera were serially diluted in two-fold, in triplicate for each dilution, beginning from 1∶12.5 through 1∶100 (to yield a final volume of 30 µL per dilution), using phosphate-buffered saline (PBS, pH 7.4) supplemented with 1% fetal bovine serum (FBS). Rubella virus stock (vaccine virus HPV77) was diluted to a working concentration of 1.2×103 plaque-forming units (PFU)/mL, and was added (30 µL) to an equal volume of diluted serum (or diluent as in the case of virus-only control), yielding a final serum dilution series of 1∶25 through 1∶200. The plate was incubated for 1.5 hour at 37°C, 5% CO2. Fifty microliters of each mixture were used to inoculate confluent Vero cell monolayers (cultured in flat-bottom 96-well plates) and the cells were incubated for 1 hour at 37°C, 5% CO2. After the incubation period, DMEM supplemented with 5% FBS and 50 µg/mL Gentamicin (Gibco; Invitrogen, CA, USA) was added to each well and the plate was further incubated for 72 hours at 37°C, 5% CO2. The plates were developed using an indirect immunocolorimetric method and rubella E1-specific monoclonal antibody (Centers for Disease Control and Prevention/CDC, GA, USA), as previously described [31]. Optical density (OD) values were measured using a measurement/reference wavelength pair of 450 nm/630 nm. The neutralization titer was calculated as the highest dilution at which the input virus signal was reduced by at least 50% within the dilution series (neutralization titer 50/NT50). The Loess method of statistical interpolation was used to estimate interpolated NT50 values from observed values [25], [32]. The intra-class correlation coefficient (ICC), based on log-transformed estimates from samples with repeated NT50 measurements was, 0.89, demonstrating a high degree of reproducibility in the assay [25], [26].
SNP Selection and Candidate SNP Genotyping
The SNP selection and genotyping methods described herein are similar or identical to those published for our previous genetic association studies [7], [8], [9], [28], [29]. The genotyping effort included follow-up on 197 previously identified genetic associations (p<0.05) between SNPs in immune response genes (cytokine and cytokine receptor genes, Toll-like receptor and signaling genes, vitamin D and vitamin A receptor family genes, antiviral effector and other innate genes) and immune measures after rubella vaccination. A total of 571 additional SNPs on the GoldenGate 768-plex panel were filled with candidate tagSNPs selected from eight newly discovered viral receptor genes (MOG, PVR, PVRL1, PVRL2, PVRL3, PVRL4, BTN2A1, BTN3A1) using the tagSNP selection approach based on linkage disequilibrium (LD), as previously described [7], [8], [9], [28], [29], [33].
The 768 SNPs were genotyped using a custom designed 768-plex Illumina GoldenGate™ assay (Illumina Inc., San Diego, CA) following the manufacturer’s instructions. The BeadStudio 2 software was used to call genotypes. The genotyping quality was high, with a mean SNP call rate of 99.7% and 99.8% (for the Rochester cohort and the San Diego cohort, respectively), and a mean subject call rate of 99.7% and 99.8%. The estimated concordance of duplicates was 96%–100%.
For the Rochester cohort, after removing samples that failed genotyping and/or QC checks, there were 1,039 samples and 606 SNPs available for analysis, including 887 samples and 555 SNPs for the Caucasian-only analysis. For the San Diego cohort, we had data (available for analysis) on 989 subjects and 611 SNPs, including 542 samples and 565 SNPs for the Caucasian-only analysis.
Statistical Methods/Analysis
Summaries of the subject characteristics were obtained within the two study cohorts, counts and percentages for categorical features and medians and interquartile ranges (25th and 75th percentiles) for continuous features. Principal components analysis of a large collection of independent SNPs was applied to estimate and define racial groups, as previously described [34], [35]. As the Rochester cohort was predominantly Caucasian, the data for both cohorts was subset to this race. Hardy–Weinberg equilibrium (HWE) was tested within the Caucasian subsets of each cohort. Principal component analysis was repeated via EIGENSTRAT [36] on a collection of SNPs except those that were in low linkage disequilibrium (r2<0.1) and that had HWE p-values>1×10−3. Non-autosomal SNPs, or SNPs that had a minor allele frequency less than 0.01, were not included in these population stratification analyses, which produced additional race and cohort-specific population stratification eigenvectors.
Linear regression models, which modeled the quantitative neutralizing antibody outcome as a linear combination of covariates and SNP genotype variables, were used to assess associations between neutralizing antibody levels and SNP genotypes. An additive genetic model, which coded SNP genotypes into a single variable reflecting the number of minor alleles carried by each individual, was used for the primary tests of significance. As sensitivity analyses, dominant and recessive genetic effects were also examined for their associations with the outcome, but for most of the SNPs these alternative models did not prove to be superior compared to the additive genetic model used in our primary assessment of significance. Neutralizing antibody (NT50) levels were analyzed as the phenotype of interest after logarithmic transformation in order to conform to linear models assumptions. SNP associations were assessed while adjusting for gender, assay batch effect, quartiles of age at enrollment, immunization age, time since last immunization to enrollment, and race-specific population stratification eigenvectors.
Associations were tested independently within each cohort, and the results were compared between the two cohorts. In order to determine which SNPs were potentially associated with neutralizing antibody titers across the two study cohorts, a fixed effects meta-analysis was performed for each SNP [37]. These meta-analyses used estimates of effect from the two study cohorts, and combined them using a fixed effects paradigm, to obtain a pooled estimate of effect, along with its pooled standard error and test of significance. Additionally, the meta-analyses provided a test of heterogeneity between cohorts. SNPs that did not display significant heterogeneity between the two cohorts (p>0.1) and that had meta-analysis p-values less than 0.05 were considered to be potentially associated with neutralizing antibody levels, and meriting further follow-up in future studies.
Analyses were carried out using the SAS software system (SAS Institute, Inc., Cary, NC) and the PLINK whole genome data analysis software package.
Results
Demographic Characteristics and Immune Variables of the Study Subjects
The demographic characteristics of the study cohorts have been reported previously [7], [9], [25], [26], [27], [28], [29], [30], [38], and are summarized in Table 1. Briefly, out of the Rochester (discovery) cohort, 1,029 subjects met all inclusion, exclusion and QC criteria, and had genotyping and rubella neutralizing antibody data available for this study. Of these 1,029 subjects, 887 (85.4%) were Caucasians, and were included in the final analysis. Out of the San Diego (replication) cohort, 985 subjects met all inclusion, exclusion and QC criteria, and had genotyping and rubella neutralizing antibody data available. Of these 985 subjects, 542 (55.03%) were Caucasians, and were included in the final analysis. The median age at enrollment was 15 years (IQR, 13, 17) for the Rochester cohort, and 23 years (IQR, 22, 27) for the San Diego cohort. The median neutralizing antibody level was 55.5 NT50 for the Rochester cohort and 62.8 NT50 for the San Diego cohort (Table 1) Detailed characterization of the immune response variables in the study cohorts and their interrelationships has been published previously [25], [26], [27].
Associations between SNPs and Antibody Titers in Two Distinct Cohorts
We found three replicated genetic associations between SNPs tagging the poliovirus receptor gene (PVR, nectin-like protein 5, CD155) in high linkage disequilibrium/LD (rs1550551, rs2165538 and rs73936843; r2 = 1, Fig. 1), and rubella-specific neutralizing antibodies (meta-analysis p-value = 0.008, Table 2). Although relatively rare, the heterozygous genotype of these PVR SNPs (no homozygous minor allele genotype was observed) was consistently associated with a significant decrease (> 2-fold) in neutralizing antibody titers after rubella vaccination in both cohorts. Similarly, the heterozygous genotype of two SNPs (rs78245864 and rs41290122, in LD; r2 = 0.66, Fig. 1) in the PVRL2 gene (poliovirus receptor-related 2, herpesvirus entry mediator B, nectin 2, CD112) was associated with a significant increase (∼ 30%) in rubella-specific neutralizing antibody titers (meta-analysis p-values of 0.038 and 0.04 for rs78245864 and rs41290122, respectively, Table 2). Furthermore, five additional PVRL2 genetic variants (some in LD, Fig. 1) demonstrated some evidence for association with rubella-specific NT50 titers post vaccination (meta-analysis p-value<0.044, Table 2), although some of the findings (association with high/low response) were not consistent between the two cohorts. Our meta-analysis also revealed potential associations between SNPs in the promoter regions of CD209/DC-SIGN (rs2287886) and IL6 (rs1880241) genes and rubella-specific neutralizing antibodies (meta-analysis p-values of 0.008 and 0.027, for rs2287886 and rs1880241, respectively, Table 2). We found a replicated genetic association for a SNP (rs1153600) in intron 3 of the vitamin A receptor gene RARB that demonstrated a significant decrease in neutralizing antibody titers with the representation of the minor allele (meta-analysis p-value = 0.037, Table 2). The homozygous minor allele genotype of a SNP (rs16895223) in the putative rubella virus receptor gene MOG (myelin oligodendrocyte glycoprotein) was associated with a decrease in rubella-specific antibodies (meta-analysis p-value = 0.046, Table 2). In addition, BTN2A1 SNP rs1977198, IRF9 SNP rs17256713 and EIF2AK2 SNPrs4648212 also demonstrated evidence for association with rubella-specific NT50 titers (p<0.049, Table 2).
A. Haplotype block structure of the PVR and PVRL2 genetic variants, analyzed using Haploview software, version 4.2 (all SNPs presented were directly genotyped in the study). The r2 color scheme is: white (r2 = 0), shades of grey (0<r2<1), black (r2 = 1). The numbers report the r2 value multiplied by 100. B. Locus zoom plot of PVR/PVRL2 genetic region on 19q13, showing multiple SNPs associated with rubella-specific neutralizing antibody levels. P-value is represented on the left Y-axis. The recombination rate is represented on the right Y-axis and as a line at the bottom of the graph. SNPs are represented as circles, color-coded for LD with respect to rs73936843 (the top three PVR SNPs rs73936843, rs1550551 and rs2165538 are in perfect LD and overlap, as represented on the graph).
The analysis results for all SNPs assessed for association with neutralizing antibody levels in our current study are presented in Table S1.
Discussion
Much knowledge has been accumulated on the role of genetic factors in the inter-individual immune response variation following immunization with measles-mumps-rubella (MMR) vaccine. A genome-wide transcriptional study (using mRNA-Seq) in high and low antibody responders to rubella vaccine identified HLA genes, and several innate immunity and inflammation-related genes (MEFV, Mediterranean fever gene, EMR3, EGF-like module containing, mucin-like, hormone receptor-like 3 gene, etc.) that discriminated between high and low humoral immune response following rubella immunization [4]. We and others have highlighted the importance of HLA polymorphisms, SNPs in cytokine and cytokine receptor genes, antiviral genes, toll-like receptors and pathway signaling genes, vitamin A and D receptor and cellular viral receptor genes (CD46 and SLAM for measles virus), and other important innate immune genes (TRIM genes, DC-SIGN, etc.) for inter-individual immune response variation after measles and rubella vaccination [2], [5], [6], [7], [8], [9], [10], [21], [22], [23], [24], [28], [29], [39], [40], [41], [42]. Because the immune response encompasses a system of interrelated components and biological processes, it is likely that multiple genetic factors and gene networks influence host response and vaccine-induced immunity (involving antibody response), including measles and rubella vaccines [5], [40].
The purpose of the current study was to refine our previous findings through replication, and to examine, in two independent cohorts, the role of newly discovered host viral receptors, attachment factors and immune regulators in the adaptive immune response induction and/or maintenance following vaccination with a live rubella vaccine.
Our meta-analysis results provide evidence for the involvement of multiple genetic variants (some in LD) within the PVR-PVRL2 gene region (Fig. 1) on chromosome 19 (19q13.31–32) in the genetic control of neutralizing antibody response to rubella vaccine. The poliovirus receptor (PVR, CD155) and poliovirus receptor-related 2 (PVRL2, CD112) genes encode transmembrane glycoproteins/nectins of the immunoglobulin (Ig) superfamily, which are components of adherens junctions and serve as cellular entry receptors for poliovirus (PVR), pseudorabies virus, and certain mutant strains of herpes simplex virus (PRRL2) [17], [18], [19]. As virus-specific cell surface receptors, these proteins have important functional downstream effects on processes related to viral entry, such as viral replication, cell-to-cell spread, viral tropism, pathogenesis and antiviral immunity. For example, the Ala67Thr mutation in the poliovirus receptor gene (PVR) was previously associated with a higher risk for developing vaccine-induced and wild-type virus-induced poliomyelitis [43].
Furthermore, both PVR (CD155) and PVRL2 (CD112) are directly involved in immunity and regulate T cell (CTL) and natural killer/NK cell activation/cytotoxicity by functioning as specific ligands for the immune costimulatory receptors DNAX Accessory Molecule-1/DNAM-1 (CD226) and CD96/Tactile and by promoting IFNγ production [15], [44], [45], [46]. Importantly, CD155/PVR may play a role in Th2 differentiation and humoral immunity by polarizing naïve CD4+T cells to a Th2 phenotype [16]. In particular, TLR agonists are known to: upregulate CD155/PVR expression (dependent on MYD88, TRIF and NF-κB); suppress Th1 cell differentiation; and modulate antigen/(OVA)-specific IgG titers in a mouse model [16]. Consistent with these findings, Maier et al. demonstrated less efficient humoral immune response (IgG and IgA) to orally administered antigens in mice deficient of CD155 [47].
In summary, based on the literature and our findings, we speculate that one or more of the multiple PVR and PVRL2 genetic variants, associated with variation in rubella vaccine-specific neutralizing antibody levels, or other tagged causal SNP/SNPs in the 19q13 region, influence the gene expression and/or protein function of CD155 (PVR) and CD112 (PVRL2). In turn, this may have downstream effects on rubella virus-specific humoral immunity by either modulating the unspecific (or specific) virus attachment/entry, replication and immune response priming or, more likely, by modulating Th1/Th2 cell differentiation and promoting more (or less) efficient antibody response [16].
The only cellular receptor for rubella virus identified to date is the myelin oligodendrocyte glycoprotein (MOG), an adhesion molecule that is important in nerve myelination and is implicated as a target antigen in the pathogenesis of several autoimmune demyelinating diseases (including multiple sclerosis/MS), and in the central nervous system (CNS) damage in congenital rubella syndrome [11]. However, rubella virus cell tropism and infectivity suggest that other, yet undiscovered, cellular receptors and adhesion factors capable of mediating infection exist for this pathogen. MOG is expressed predominantly in the central nervous system and is hardly detectable in immune system organs/tissues and cells such as spleen, thymus and PBMCs [11]. Our analysis identified two intronic polymorphisms (rs16895223, p = 0.046 and rs1977198, p = 0.027, Table 2) in the MOG gene and the highly homologous butyrophilin BTN2A1 gene, both on chromosome 6p22, that were associated with variations in rubella-specific neutralizing antibody response. Host genetic variation in these or other tagged functional MOG/BTN2A1 SNPs may potentially impair or enhance viral entry into susceptible cells, alter virus replication/propagation and antigen abundance and modulate subsequent immune response. Alternatively, these genes may interfere in immune regulation, as demonstrated for the butyrophilins. The latter are co-stimulatory molecules that are implicated in T cell inhibition (negative regulation of T cell proliferation, activation, cytokine production), immune signaling and interaction with DC-SIGN on dendritic cells (BTN2A1 is identified as a ligand for DC-SIGN), and genetic variation in these genes have been associated with inflammatory diseases [48], [49].
Of note, our results also demonstrate an association between a promoter SNP in DC-SIGN (rs2287886, p = 0.008, Table 2) and humoral immune response variations following rubella vaccination. DC-SIGN (Cell-Specific Intercellular adhesion molecule-3-Grabbing Non-integrin, CD209) is a C-type lectin that mediates dendritic cell function and activation of CD4+ T cells and binds multiple microorganisms, including viruses, by recognizing mannose type glycoproteins [20]. This molecule is exploited by many viruses for attachment and/or cellular entry: measles virus (MV); rift valley fever virus (RVFV); influenza A viruses; human immunodeficiency virus (HIV-1); hepatitis C virus (HCV); herpes simplex virus (HSV1); human cytomegalovirus (HCMV); Ebola virus; SARS coronavirus and dengue virus [20]. Polymorphisms in DC-SIGN (and, in particular, the promoter SNP rs2287886) have been associated with the clinical course and outcome of dengue virus infection, cytomegalovirus infection, tick-borne encephalitis virus (TBEV) infection and pulmonary aspergillosis [50], [51], [52]. Furthermore, the promoter SNP rs2287886 is functional and was demonstrated to influence DC-SIGN gene expression and HCMC infection in dendritic cells [52]. Interestingly, we found the same DC-SIGN promoter SNP to be associated with variations in measles vaccine-induced TNFα secretion in African-Americans and several other DC-SIGN SNPs that were associated with variations in measles-specific neutralizing antibody levels in Caucasians and African-Americans [21]. Therefore, it is likely that DC-SIGN, PVR, PVRL2 or other attachment factors and/or cellular viral receptors and pathway-related genes are involved in the genetic control of rubella vaccine-induced immune response heterogeneity after vaccination.
Another important finding is the observed genetic association between a SNP (rs1153600) in intron 3 of the retinoic acid receptor, beta (RARB) gene and rubella-specific neutralizing antibody response (the homozygous minor allele and the heterozygous genotypes were associated with lower antibody titers in both cohorts). RARB is a receptor for retinoic acid, the biologically active metabolite of vitamin A, with involvement in cell growth, differentiation, cell signaling and gene regulation. Interestingly, vitamin A has direct effects on the immune system and regulates antigen presentation, lymphocyte homing, proliferation of B and T cells, T-helper differentiation, cytokine production, T cell activation and cytotoxicity, as well as enhancing the antibody response to vaccination, including the measles vaccine, oral polio vaccine, and tetanus and diphtheria toxoids [53], [54]. Furthermore, our previous studies demonstrated associations between: the same RARB SNP (rs1153600, p = 0.034) and two other RARB SNPs and rubella (whole) virus-specific IgG levels (as measured by a chemiluminescent immunoassay) in 738 children/young adults following two doses of MMR vaccine [8]; and RARB SNPs/haplotypes and variations in measles-specific neutralizing antibody levels and cellular immune outcomes after MMR vaccination [41]. The findings from the current study and previous reports add evidence to the importance of vitamin A receptors and pathway in immune response regulation after immunization.
Other potential findings include SNP associations between polymorphisms in the IL6 gene (promoter region SNP rs1880241), IRF9 gene (3′intergenic SNP rs17256713) and EIF2AK2 gene (intronic SNP rs4648212), which we previously found to be associated with whole rubella virus-specific IgG levels using single-SNP and/or multigenic assessment analysis [5], [7]. Other SNP associations identified in our previous studies were not confirmed in the current study, possibly due to differences in study design, sample size and cohorts’ characteristics, analytical approach, differences in immune response measurement (rubella whole virus-specific IgG levels vs. rubella virus-specific neutralizing antibody levels), and/or false-positive findings [4], [5], [6], [7], [8], [9], [10].
While our efforts have identified a collection of genetic variants that are of high interest as potentially playing a role in variability of neutralizing antibody response after rubella vaccination, it must be noted that some of the genetic associations may not be true positives. The limitations of our study reflect the fact that the characteristics of the study participants were different between the two study cohorts. This made it difficult to compare and/or combine the results from the two cohorts using a method that would typically be applied to control for false-positive findings. For this reason, we applied meta-analytic approaches to evaluate the degree of evidence of association present when evaluating combined results across the two cohorts, with the understanding that false-positive results would be less likely to show consistent associations between cohorts. Given the signals that we observed, and the biological plausibility, it is likely that at least several of the identified associations/genes do indeed play a role in influencing neutralizing antibody levels observed in vaccinated individuals. In addition, our analysis focused on additive genetic effects, and while the support for other genetic models over the additive model was not compellingly strong in our sensitivity analyses, it is possible that for some of the genetic variants an alternative genetic model may prove superior. Other limitations include the influence of unpredictable factors confounding analysis results (e.g., potential wild-type rubella virus exposure in deployed military personnel from the San Diego cohort).
The strengths of our study include the use of a clinically relevant immune outcome to assess vaccine response (i.e., rubella-specific functional/neutralizing antibody levels measured using a state-of-the-art, standardized, high-throughput immune assay). A significant strength is the genetic association analysis approach, adjusting for known confounding variables, and the use of two independent cohorts (subsetted to Caucasian subjects only) to minimize the false-positive findings. The identified genes and genetic variants will be further fine mapped in order to identify candidate functional variants in the genomic regions tagged by the original SNP (if not functional). The p-values, and magnitudes of effect, obtained from the fine-mapping analyses will provide important insights into which of the variants is most likely to be causing the genetic association that was originally observed. Finally, the identification of likely causal genetic variants will be complemented with functional studies to reveal the SNP immediate functional effect and the downstream consequences/biological mechanisms for humoral immune response variation.
In summary, our results identified multiple genetic variants, primarily in genes related to viral attachment and entry, and/or immune regulation that are associated with inter-individual differences in rubella-specific neutralizing antibody response after vaccination. Such findings may indicate evidence for novel receptors used by the rubella virus for cell entry, and will assist in development of improved vaccines and vaccine formulations for achieving optimal immune response after vaccination.
Supporting Information
Table S1.
SNPs assessed for association with neutralizing antibody levels after rubella vaccination in two different cohorts.
https://doi.org/10.1371/journal.pone.0099997.s001
(DOCX)
Author Contributions
Conceived and designed the experiments: IH NL IO RK GAP. Performed the experiments: IH NL. Analyzed the data: BL VSP IH. Contributed reagents/materials/analysis tools: IO RK GAP BL VSP. Wrote the paper: IH NL IO RK GAP VSP.
References
- 1. Poland GA, Ovsyannikova IG, Kennedy RB, Haralambieva IH, Jacobson RM (2011) Vaccinomics and a New Paradigm for the Development of Preventive Vaccines Against Viral Infections. Omics 15: 625–636.
- 2. Haralambieva IH, Ovsyannikova IG, Pankratz VS, Kennedy RB, Jacobson RM, et al. (2013) The genetic basis for interindividual immune response variation to measles vaccine: new understanding and new vaccine approaches. Expert Review of Vaccines 12: 57–70.
- 3. Tan PL, Jacobson RM, Poland GA, Jacobsen SJ, Pankratz SV (2001) Twin studies of immunogenicity - determining the genetic contribution to vaccine failure. Vaccine 19: 2434–2439.
- 4. Haralambieva IH, Oberg AL, Ovsyannikova IG, Kennedy RB, Grill DE, et al. (2013) Genome-wide characterization of transcriptional patterns in high and low antibody responders to rubella vaccination. PLos ONE 8: e62149.
- 5. Pankratz VS, Vierkant RA, O'Byrne MM, Ovsyannikova IG, Poland GA (2010) Associations between SNPs in candidate immune-relevant genes and rubella antibody levels: a multigenic assessment. BMC Immunol 11: 48.
- 6. Dhiman N, Haralambieva IH, Kennedy RB, Vierkant RA, O'Byrne MM, et al. (2010) SNP/haplotype associations in cytokine and cytokine receptor genes and immunity to rubella vaccine. Immunogenetics 62: 197–210.
- 7. Haralambieva IH, Dhiman N, Ovsyannikova IG, Vierkant RA, Pankratz VS, et al. (2010) 2′-5′-Oligoadenylate synthetase single-nucleotide polymorphisms and haplotypes are associated with variations in immune responses to rubella vaccine. Hum Immunol 71: 383–391.
- 8. Ovsyannikova IG, Haralambieva IH, Dhiman N, O'Byrne MM, Pankratz VS, et al. (2010) Polymorphisms in the vitamin A receptor and innate immunity genes influence the antibody response to rubella vaccination. J Infect Dis 201: 207–213.
- 9. Ovsyannikova IG, Dhiman N, Haralambieva IH, Vierkant RA, O'Byrne MM, et al. (2010) Rubella vaccine-induced cellular immunity: evidence of associations with polymorphisms in the Toll-like, vitamin A and D receptors, and innate immune response genes. Human Genetics 127: 207–221.
- 10. Ovsyannikova IG, Jacobson RM, Vierkant RA, O'Byrne MM, Poland GA (2009) Replication of rubella vaccine population genetic studies: validation of HLA genotype and humoral response associations. Vaccine 27: 6926–6931.
- 11. Cong H, Jiang Y, Tien P (2011) Identification of the myelin oligodendrocyte glycoprotein as a cellular receptor for rubella virus. Journal of Virology 85: 11038–11047.
- 12. Noyce RS, Bondre DG, Ha MN, Lin LT, Sisson G, et al. (2011) Tumor Cell Marker PVRL4 (Nectin 4) Is an Epithelial Cell Receptor for Measles Virus. PLoS Pathog 7: e1002240.
- 13. Muhlebach MD, Mateo M, Sinn PL, Prufer S, Uhlig KM, et al. (2011) Adherens junction protein nectin-4 is the epithelial receptor for measles virus. Nature 480: 530–533.
- 14. Lozach PY, Kuhbacher A, Meier R, Mancini R, Bitto D, et al. (2011) DC-SIGN as a receptor for phleboviruses. Cell Host & Microbe 10: 75–88.
- 15. Bottino C, Castriconi R, Pende D, Rivera P, Nanni M, et al. (2003) Identification of PVR (CD155) and Nectin-2 (CD112) as cell surface ligands for the human DNAM-1 (CD226) activating molecule. The Journal of Experimental Medicine 198: 557–567.
- 16. Kamran N, Takai Y, Miyoshi J, Biswas SK, Wong JS, et al. (2013) Toll-like receptor ligands induce expression of the costimulatory molecule CD155 on antigen-presenting cells. PLos ONE 8: e54406.
- 17. He Y, Bowman VD, Mueller S, Bator CM, Bella J, et al. (2000) Interaction of the poliovirus receptor with poliovirus. Proceedings of the National Academy of Sciences of the United States of America 97: 79–84.
- 18. Warner MS, Geraghty RJ, Martinez WM, Montgomery RI, Whitbeck JC, et al. (1998) A cell surface protein with herpesvirus entry activity (HveB) confers susceptibility to infection by mutants of herpes simplex virus type 1, herpes simplex virus type 2, and pseudorabies virus. Virology 246: 179–189.
- 19. Struyf F, Martinez WM, Spear PG (2002) Mutations in the N-terminal domains of nectin-1 and nectin-2 reveal differences in requirements for entry of various alphaherpesviruses and for nectin-nectin interactions. Journal of Virology 76: 12940–12950.
- 20. Backovic M, Rey FA (2012) Virus entry: old viruses, new receptors. Current Opinion in Virology 2: 4–13.
- 21. Ovsyannikova IG, Haralambieva IH, Vierkant RA, O'Byrne MM, Jacobson RM, et al. (2011) The Association of CD46, SLAM, and CD209 Cellular Receptor Gene SNPs with Variations in Measles Vaccine-Induced Immune Responses–A Replication Study and Examination of Novel Polymorphisms. Human Heredity 72: 206–223.
- 22. Dhiman N, Cunningham JM, Jacobson RM, Vierkant RA, Wu Y, et al. (2007) Variations in measles vaccine-specific humoral immunity by polymorphisms in SLAM and CD46 measles virus receptors. Journal of Allergy and Clinical Immunology 120: 666–672.
- 23. Clifford HD, Hayden CM, Khoo SK, Zhang G, Le Souef PN, et al. (2012) CD46 measles virus receptor polymorphisms influence receptor protein expression and primary measles vaccine responses in naive Australian children. ClinVaccine Immunol 19: 704–710.
- 24. Clifford HD, Richmond P, Khoo SK, Zhang G, Yerkovich ST, et al. (2011) SLAM and DC-SIGN measles receptor polymorphisms and their impact on antibody and cytokine responses to measles vaccine. Vaccine 29: 5407–5413.
- 25. Lambert ND, Pankratz VS, Larrabee BR, Ogee-Nwankwo A, Chen MH, et al. (2014) High-throughput Assay Optimization and Statistical Interpolation of Rubella-Specific Neutralizing Antibody Titers. Clinical and Vaccine Immunology 21: 340–346.
- 26. Lambert ND, Haralambieva IH, Ovsyannikova IG, Larrabee BR, Pankratz VS, et al. (2013) Characterization of humoral and cellular immunity to rubella vaccine in four distinct cohorts. Immunologic Research 58: 1–8.
- 27. Haralambieva IH, Salk HM, Lambert ND, Ovsyannikova IG, Kennedy RB, et al. (2014) Associations between race, sex and immune response variations to rubella vaccination in two independent cohorts. Vaccine 32: 1946–1953.
- 28. Haralambieva IH, Ovsyannikova IG, Umlauf BJ, Vierkant RA, Pankratz SV, et al. (2011) Genetic polymorphisms in host antiviral genes: associations with humoral and cellular immunity to measles vaccine. Vaccine 29: 8988–8997.
- 29. Haralambieva IH, Ovsyannikova IG, Kennedy RB, Vierkant RA, Pankratz SV, et al. (2011) Associations between single nucleotide polymorphisms and haplotypes in cytokine and cytokine receptor genes and immunity to measles vaccination. Vaccine 29: 7883–7895.
- 30. Haralambieva IH, Ovsyannikova IG, Dhiman N, Kennedy RB, O'Byrne M, et al. (2011) Common SNPs/Haplotypes in IL18R1 and IL18 Genes are Associated with Variations in Humoral Immunity to Smallpox Vaccination in Caucasians and African-Americans. Journal of Infectious Diseases 204: 433–441.
- 31. Chen MH, Zhu Z, Zhang Y, Favors S, Xu WB, et al. (2007) An indirect immunocolorimetric assay to detect rubella virus infected cells. Journal of Virological Methods 146: 414–418.
- 32.
Cleveland WS, Gross E, Shyu WM (1992) Local regression models. In: Chambers JM, Hastie TJ, editors. Statistical Models: Wadsworth & Brooks/Cole.
- 33. Carlson CS, Eberle MA, Rieder MJ, Yi Q, Kruglyak L, et al. (2004) Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium. Am J Hum Genet 74: 106–120.
- 34. Kennedy RB, Ovsyannikova IG, Pankratz VS, Haralambieva IH, Vierkant RA, et al. (2012) Genome-wide genetic associations with IFNgamma response to smallpox vaccine. Human Genetics 131: 1433–1451.
- 35. Haralambieva IH, Ovsyannikova IG, Kennedy RB, Larrabee BR, Shane Pankratz V, et al. (2013) Race and sex-based differences in cytokine immune responses to smallpox vaccine in healthy individuals. Human Immunology 74: 1263–1266.
- 36. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, et al. (2006) Principal components analysis corrects for stratification in genome-wide association studies. Nature Genetics 38: 904–909.
- 37. Greenland S, O'Rourke K (2008) Meta-Analysis. In: Rothman KJ, Greeland S, Lash T, editors. Modern Epidemiology. 3rd ed: Lippincott Williams and Wilkins. 652.
- 38. Haralambieva IH, Ovsyannikova IG, O'Byrne M, Pankratz VS, Jacobson RM, et al. (2011) A large observational study to concurrently assess persistence of measles specific B-cell and T-cell immunity in individuals following two doses of MMR vaccine. Vaccine 29: 4485–4491.
- 39. Ovsyannikova IG, Haralambieva IH, Vierkant RA, O'Byrne MM, Poland GA (2013) Associations between polymorphisms in the antiviral TRIM genes and measles vaccine immunity. Human Immunology 74: 768–774.
- 40. Kennedy RB, Ovsyannikova IG, Haralambieva IH, O'Byrne MM, Jacobson RM, et al. (2012) Multigenic control of measles vaccine immunity mediated by polymorphisms in measles receptor, innate pathway, and cytokine genes. Vaccine 30: 2159–2167.
- 41. Ovsyannikova IG, Haralambieva IH, Vierkant RA, O'Byrne MM, Jacobson RM, et al. (2012) Effects of Vitamin A and D Receptor Gene Polymophisms/Haplotypes on Immune Responses to Measles Vaccine. Pharmacogenetics and Genomics 22: 20–31.
- 42. Ovsyannikova IG, Haralambieva IH, Vierkant RA, Pankratz VS, Poland GA (2011) The role of polymorphisms in toll-like receptors and their associated intracellular signaling genes in measles vaccine immunity. Human Genetics 130: 547–561.
- 43. Kindberg E, Ax C, Fiore L, Svensson L (2009) Ala67Thr mutation in the poliovirus receptor CD155 is a potential risk factor for vaccine and wild-type paralytic poliomyelitis. Journal of Medical Virology 81: 933–936.
- 44. Shibuya A, Campbell D, Hannum C, Yssel H, Franz-Bacon K, et al. (1996) DNAM-1, a novel adhesion molecule involved in the cytolytic function of T lymphocytes. Immunity 4: 573–581.
- 45. Fuchs A, Cella M, Giurisato E, Shaw AS, Colonna M (2004) Cutting edge: CD96 (tactile) promotes NK cell-target cell adhesion by interacting with the poliovirus receptor (CD155). Journal of Immunology 172: 3994–3998.
- 46. Tahara-Hanaoka S, Shibuya K, Onoda Y, Zhang H, Yamazaki S, et al. (2004) Functional characterization of DNAM-1 (CD226) interaction with its ligands PVR (CD155) and nectin-2 (PRR-2/CD112). International Immunology 16: 533–538.
- 47. Maier MK, Seth S, Czeloth N, Qiu Q, Ravens I, et al. (2007) The adhesion receptor CD155 determines the magnitude of humoral immune responses against orally ingested antigens. European Journal of Immunology 37: 2214–2225.
- 48. Abeler-Dorner L, Swamy M, Williams G, Hayday AC, Bas A (2012) Butyrophilins: an emerging family of immune regulators. Trends in Immunology 33: 34–41.
- 49. Malcherek G, Mayr L, Roda-Navarro P, Rhodes D, Miller N, et al. (2007) The B7 homolog butyrophilin BTN2A1 is a novel ligand for DC-SIGN. Journal of Immunology 179: 3804–3811.
- 50. Alagarasu K, Damle IM, Bachal RV, Mulay AP, Shah PS, et al. (2013) Association of promoter region polymorphisms of CD209 gene with clinical outcomes of dengue virus infection in Western India. Infection, Genetics and Evolution: Journal of Molecular Epidemiology and Evolutionary Genetics in Infectious Diseases 17: 239–242.
- 51. Barkhash AV, Perelygin AA, Babenko VN, Brinton MA, Voevoda MI (2012) Single nucleotide polymorphism in the promoter region of the CD209 gene is associated with human predisposition to severe forms of tick-borne encephalitis. Antiviral Research 93: 64–68.
- 52. Mezger M, Steffens M, Semmler C, Arlt EM, Zimmer M, et al. (2008) Investigation of promoter variations in dendritic cell-specific ICAM3-grabbing non-integrin (DC-SIGN) (CD209) and their relevance for human cytomegalovirus reactivation and disease after allogeneic stem-cell transplantation. Clinical microbiology and infection: the official publication of the European Society of Clinical Microbiology and Infectious Diseases 14: 228–234.
- 53. Villamor E, Fawzi WW (2005) Effects of vitamin a supplementation on immune responses and correlation with clinical outcomes. Clin Microbiol Rev 18: 446–464.
- 54. Mora JR, Iwata M, Von Andrian UH (2008) Vitamin effects on the immune system: vitamins A and D take centre stage. Nat Rev Immunol 8: 685–698.