Kawasaki disease (KD) is a diffuse and acute small-vessel vasculitis observed in children and has genetic and autoimmune components. We genotyped 112 case-parent trios of European decent (confirmed by AIMS) using the ImmunoChip array and performed association analyses with susceptibility to KD and IVIG non-response. KD susceptibility was assessed using the transmission disequilibrium test whereas IVIG non-response was evaluated using multivariable logistic regression analysis. We replicated SNPs in three gene regions (FCGR, CD40/CDH22, and HLA-DQB2/HLA-DOB) that have been previously associated with KD and provide support to other findings of several novel SNPs in genes with potential pathway in KD pathogenesis. SNP rs838143 in the 3′ UTR of FUT1 gene (2.7×10-5) and rs9847915 in the intergenic region of LOC730109 ∣ BRD7P2 (6.81×10-7) were the top hits for KD susceptibility in additive and dominant models, respectively. The top hits for IVIG responsiveness were rs1200332 in the intergenic region of BAZ1A ∣ C14orf19 (1.4×10-4) and rs4889606 in the intron of the STX1B gene (6.95×10-5) in additive and dominant models, respectively. Our study suggests that genes and biological pathways involved in autoimmune diseases play an important role in the pathogenesis of KD and IVIG response mechanism.
Mucocutaneous lymph node syndrome; Kawasaki disease; intravenous immunoglobulins therapy; Immunogenetics; Immune-related genes
Kawasaki disease (KD), response to intravenous immunoglobulin (IVIG) therapy, and associated coronary artery disease progression have been associated with genetic polymorphisms in Fc gamma receptor (FcγR) genes. However, it is not known whether the existing gene copy number (GCN) variability relates to KD treatment response, susceptibility, or associated sequelae.
The copy number of individuals with KD (n = 510) and their family members (n = 808) for three variable FcγRs was assessed using pyrosequencing. We performed the transmission disequilibrium test to examine the association of GCN for FcγRs (FcγR2C, FcγR3A, and FcγR3B) with susceptibility and used logistic regression models to determine its association with IVIG treatment outcomes.
FcγR2C and FcγR3B GCN were significantly associated with KD susceptibility. IVIG response was associated with GCN variations of FcγR3B in Whites and FcγR2C in Hispanics, and gene risk score based on single nucleotide polymorphism and GCN in FcγRs were significantly different between IVIG responders and nonresponders among Whites. We found no significant associations between coronary artery disease and any of the FcγR copy numbers.
GCN of FcγR2C and FcγR3B influences IVIG treatment response and predisposes individuals to KD, providing potential insights into understanding the mechanism of the FcγR gene family in the IVIG pathway.
copy number variation; FcγR; genetic risk; intravenous immunoglobulin; Kawasaki disease
Single nucleotide polymorphisms (SNPs) in the Ryanodine receptor 3 (RYR3) gene are associated with common carotid intima media thickness (CCA cIMT) in HIV-infected men. We evaluated SNPs in the RYR3 gene among HIV-infected women participating in Women’s Interagency HIV Study (WIHS).
CCA cIMT was measured using B-mode ultrasound and the 838 SNPs in the RYR3 gene region were genotyped using the Illumina HumanOmni2.5-quad beadchip. The CCA cIMT genetic association was assessed using linear regression analyses among 1213 women and also separately among White (n=139), Black (n=720) and Hispanic (n=354) women after adjusting for confounders. A summary measure of pooled association was estimated using a meta-analytic approach by combining the effect estimates from the three races. Haploblocks were inferred using Gabriel’s method and haplotype association analyses were conducted among the three races separately.
SNP rs62012610 was associated with CCA cIMT among the Hispanics (p=4.41× 10−5), rs11856930 among Whites (p=5.62× 10−4), and rs2572204 among Blacks (p=2.45× 10−3). Meta-analysis revealed several associations of SNPs in the same direction and of similar magnitude, particularly among Blacks and Hispanics. Additionally, several haplotypes within three haploblocks containing SNPs previously related with CCA cIMT were also associated in Whites and Hispanics.
Consistent with previous research among HIV-infected men, SNPs within the RYR3 region were associated with subclinical atherosclerosis among HIV-infected women. Allelic heterogeneity observed across the three races suggests that the contribution of the RYR3 gene to CCA cIMT is complex, and warrants future studies to better understand regional SNP function.
RYR3; single nucleotide polymorphisms; HIV infection; CCA; cIMT; subclinical atherosclerosis
The major histocompatibility complex (MHC) region on chromosome 6p21.3 is suspected to host susceptibility loci for HIV-related Kaposi’s sarcoma (HIV-KS). A nested case-control study in the Multicenter AIDS Cohort Study was designed to conduct fine genetic association mapping across central MHC. Individuals co-infected with HIV-1 and HHV-8 who later developed KS were defined as cases (n=354) and were matched 1:1 with co-infected KS-free controls.
We report data for new independent MHC class II and III susceptibility loci. In particular, class II HLA-DMB emerged as a strong candidate, with the intronic variant rs6902982 A>G associated with a 4-fold increase of risk (OR= 4.09; 95% CI: 1.90–8.80; p= 0.0003). A striking multiplicative effect on the estimated risk was associated with further carriage of two non-synonymous variants, rs1800453 A>G (Asp697Gly) and rs4148880 A>G (Ile393Val), in the linked TAP1 gene (OR=10.5; 95% CI: 2.54–43.6; p=0.0012). The class III susceptibility variant is moderately associated with HIV-KS and lies within a 120 Kb-long haplotype (OR=1.52; 95% CI: 1.01–2.28; p=0.047) formed by rs7029 A>G (GPANK1 3’UTR), rs1065356 G>A (LY6G6C), rs3749953 A>G (MSH5-SAPCD1 readthrough) and rs707926 G>A (VARS). Our data suggest that antigen processing by MHC class II molecules is a target pathway in the pathogenesis of HIV-KS.
Like other members of the γ-herpesvirus family, human herpes virus 8 (HHV-8), the etiologic agent of classic and HIV-related Kaposi’s sarcoma (HIV-KS) acquired and evolved several human genes with key immune modulatory and cellular growth control functions. The encoded viral homologs substitute for their human counterparts but escape cellular regulation, leading to uncontrolled cell proliferation. We postulated that DNA variants in the human homologs of viral genes that potentially alter the expression or the binding of the encoded factors controlling the antiviral response may facilitate viral interference. To test whether cellular homologs are candidate susceptibility genes, we evaluated the association of DNA variants in 92 immune-related genes including 7 cellular homologs with the risk for HIV-KS in a matched case and control study nested in the Multicenter AIDS Cohort Study. Low- and high-risk gene-by-gene interactions were estimated by multifactor dimensionality reduction and used as predictors in conditional logistic models. Among the most significant gene interactions at risk (OR=2.84–3.92; Bonferroni-adjusted p= 9.9×10−3−2.6×10−4), three comprised human homologs of two latently expressed viral genes, cyclin D1 (CCND1) and interleukin-6 (IL-6), in conjunction with angiogenic genes (VEGF, EDN-1 and EDNRB). At lower significance thresholds (adjusted p < 0.05), human homologs related to apoptosis (CFLAR) and chemotaxis (CCL2) emerged as candidates. This “proof of concept” study identified human homologs involved in the regulation of type I interferon-induced signaling, cell cycle and apoptosis potentially as important determinants of HIV-KS
Kaposi’s sarcoma; Immunodeficiency; Herpes Virus 8; Multifactor Dimensionality Reduction; Polymorphism; Genetic association
Our work aimed to examine the potential influence of variants in interleukin/interleukin receptors genes on high-risk (HR-HPV) HPV clearance. Clearance of genital HR-HPV infection was evaluated for 134 HIV-1 seropositive African-American female adolescents from the Reaching for Excellence in Adolescent Care and Health (REACH) cohort. Genotyping targeted 225 single nucleotide polymorphisms (SNPs) within the exons, 5′ untranslated region (UTR) and 3′ UTR sequences of 27 immune-related candidate genes encoding interleukin family of cytokines. Cox proportional hazard models were used to determine the association of type- specific HPV clearance adjusting for time-varying CD4+ T-cell count and low-risk (LR-HPV) HPV co-infections. HR-HPV clearance rates were significantly (p< 0.001) associated with five SNPs (rs228942, rs419598, rs315950, rs7737000, rs9292618) mapped to coding and regulatory regions in three genes (IL2RB, IL1RN, and IL7R). These data suggest that the analyzed genetic variants in interleukin family of cytokines modulate HR-HPV clearance in HIV-1 seropositive African-Americans that warrants replication.
HPV clearance; genetic association; interleukins; HIV-1 seropositive; African American adolescents
Nepal has one of the highest cervical cancer rates in South Asia. Only a few studies in populations from urban areas have investigated type specific distribution of human papillomavirus (HPV) in Nepali women. Data on high-risk HPV (HR-HPV) types are not currently available for rural populations in Nepal. We aimed to assess the distribution of HR- HPV among rural Nepali women while assessing self-collected and clinician-collected cervico-vaginal specimens as sample collection methods for HPV screening.
Study participants were recruited during a health camp conducted by Nepal Fertility Care Center in Achham District of rural far western Nepal. Women of reproductive age completed a socio-demographic and clinical questionnaire, and provided two specimens; one cervical-vaginal specimen using a self-collection method and another cervical specimen collected by health camp auxiliary nurse midwives during a pelvic examination. All samples were tested for 14 different HR-HPV mRNA and also specific for HPV16/18/45 mRNA.
Of 261 women with both clinician- and self-collected cervical samples, 25 tested positive for HR-HPV, resulting in an overall HR-HPV prevalence of 9.6% (95% confidence Interval [CI]: 6.3–13.8). The overall Kappa value assessing agreement between clinician- and self-collected tests was 0.62 (95% CI: 0.43–0.81), indicating a “good” level of agreement. Abnormal cytology was reported for 8 women. One woman identified with squamous cell carcinoma (SCC), and 7 women with high grade squamous intraepithelial lesions (HSIL). Seven of the 8 women tested positive for HR-HPV (87.5%) in clinician-collected samples and 6 in self-collected samples (75.0%).
This is the first study to assess HR-HPV among rural Nepali women. Self-collected sampling methods should be the subject of additional research in Nepal for screening HR-HPV, associated with pre-cancer lesions and cancer, in women in rural areas with limited access to health services.
Following the publication of the ENCODE project results, there has been increasing interest in investigating different areas of the chromosome and evaluating the relative contribution of each area to expressed phenotypes. This study aims to evaluate the contribution of variants, classified by minor allele frequency and gene annotation, to the observed interindividual differences. In this study, we fitted Bayesian linear regression models to data from Genetic Analysis Workshop 18 (n = 395) to estimate the variance of standardized and log-transformed systolic blood pressure that can be explained by subsets of genetic markers. Rare and very rare variants explained an overall higher proportion of the variance, as did markers located within a gene rather than flanking regions. The proportion of variance explained by rare and very rare variants decreased when we controlled for the number of markers, suggesting that the number of contributing rare alleles plays an important role in the genetic architecture of chronic disease traits. Our findings lend support to the "common disease, rare variant" hypothesis for systolic blood pressure and highlight allele frequency and functional annotation of a polymorphism as potentially crucial considerations in whole genome study designs.
Persistent high-risk human papillomavirus (HR-HPV) is a necessary and causal factor of cervical cancer. Most women naturally clear HPV infections; however, the biological mechanisms related to HPV pathogenesis have not been clearly elucidated. Host genetic factors that specifically regulate immune response could play an important role. All HIV-positive women in the HIV Epidemiology Research Study (HERS) with a HR-HPV infection and at least one follow-up biannual visit were included in the study. Cervicovaginal lavage samples were tested for HPV using type-specific HPV hybridization assays. Type-specific HPV clearance was defined as two consecutive HPV-negative tests after a positive test. DNA from participants was genotyped for 196,524 variants within 186 known immune related loci using the custom ImmunoChip microarray. To assess the influence of each single-nucleotide polymorphism (SNP) with HR-HPV clearance, the Cox proportional hazards model with the Wei-Lin-Weissfeld approach was used, adjusting for CD4+ count, low risk HPV (LR-HPV) co-infection, and relevant confounders. Three analytical models were performed: race-specific (African Americans (n = 258), European Americans (n = 87), Hispanics (n = 55), race-adjusted combined analysis, and meta-analysis of pooled independent race-specific analyses. Women were followed for a median time of 1,617 days. Overall, three SNPs (rs1112085, rs11102637, and rs12030900) in the MAGI-3 gene and one SNP (rs8031627) in the SMAD3 gene were associated with HR-HPV clearance (p<10−6). A variant (rs1633038) in HLA-G were also significantly associated in African American. Results from this study support associations of immune-related genes, having potential biological mechanism, with differential cervical HR-HPV infection outcomes.
The Centers for Disease Control (CDC) testing recommendations suggest universal opt-out testing in all health care settings, including cancer clinics. The incidence of non-AIDS-defining cancers (NADCs) is on the rise among HIV patients. However, to date, no data exist on the prevalence of HIV infection among NADC patients in the United States. Knowledge of HIV infection may affect clinical management, prognosis, and overall patient survival and decrease new infections in the population. The purpose of this study was to determine the point seroprevalence of HIV infection in cancer patients being seen in medical oncology clinics. A total of 634 individuals (mean age=53.2 years) participated and were tested for HIV. None of the participants tested positive for HIV in any of the three clinics. Using a futility analysis, the upper end of the 95% confidence interval for prevalence of undiagnosed HIV in cancer patients was less than 0.3%. Most participants were female (59.2%) and non-Hispanic (96.6%). The majority of study participants were white (76.5%) or African-American (17.7%). Breast cancer (19.7%), colon cancer (10.3%), and melanoma (9.7%) were the most commonly reported non-AIDS-defining cancers. While our study suggested that there was no occurrence of undiagnosed HIV among NADC patients, it is important to note that our population was largely white, females with insurance and with a different distribution of cancer than the most prevalent NADC among HIV patients. Furthermore, one-third of the patients did not consent to participate and further studies are needed to assess reasons for their unwillingness along with other populations, specifically minorities and individuals with low socioeconomic status (SES).
Nearly one-third of adults in the U.S. have hypertension, which is associated with increased cardiovascular disease (CVD) morbidity and mortality. The goal of antihypertensive pharmacogenetic research is to enhance understanding of drug response based on the interaction of individual genetic architecture and antihypertensive therapy to improve blood pressure control and ultimately prevent CVD outcomes. In the context of the Genetics of Hypertension Associated Treatment (GenHAT) study and using a case-only design, we examined whether single nucleotide polymorphisms in RYR3 interact with four classes of antihypertensive drugs, particularly the calcium channel blocker amlodipine versus other classes, to modify the risk of coronary heart disease (CHD; fatal CHD and non-fatal myocardial infarction combined) and heart failure in high-risk hypertensive individuals. RYR3 mediates the mobilization of stored Ca+2 in cardiac and skeletal muscle to initiate muscle contraction. There was suggestive evidence of pharmacogenetic effects on heart failure, the strongest of which was for rs877087, with the smallest p-value =.0005 for the codominant model when comparing amlodipine versus all other treatments. There were no pharmacogenetic effects observed for CHD. The findings reported here for the case-only analysis of the antihypertensive pharmacogenetic effect of RYR3 among 3,058 CHD cases and 1,940 heart failure cases show that a hypertensive patient’s genetic profile may help predict which medication(s) might better lower cardiovascular disease risk.
RYR3 gene; calcium channel blocker; hypertension; coronary heart disease; heart failure; genetic interaction
Genetic changes occurring in different stages of pre-cancer lesions reflect causal events initiating and promoting the progression to cancer. Co-existing pre-cancerous lesions including low- and high-grade squamous intraepithelial lesion (LGSIL and HGSIL), and adjacent “normal” cervical epithelium from six formalin-fixed paraffin-embedded samples were selected. Tissues from these 18 samples were isolated using laser-capture microdissection, RNA was extracted and sequenced. RNA-sequencing generated 2.4 billion raw reads in 18 samples, of which ~50.1% mapped to known and annotated genes in the human genome. There were 40 genes up-regulated and 3 down-regulated (normal to LGSIL) in at least one-third of the sample pairs (same direction and FDR p < 0.05) including S100A7 and KLK6. Previous studies have shown that S110A7 and KLK7 are up-regulated in several other cancers, whereas CCL18, CFTR, and SLC6A14, also differentially expressed in two samples, are up-regulated specifically in cervical cancer. These differentially expressed genes in normal to LGSIL progression were enriched in pathways related to epithelial cell differentiation, keratinocyte differentiation, peptidase, and extracellular activities. In progression from LGSIL to HGSIL, two genes were up-regulated and five down-regulated in at least two samples. Further investigations using co-existing samples, which account for all internal confounders, will provide insights to better understand progression of cervical pre-cancer.
RNA-sequencing; gene expression; cervical dysplasia; co-existing lesions; human genome
Bacterial vaginosis (BV) is a common vaginal disorder in women of reproductive age, especially among women with HIV-1 infection. Several bacterial products including lipopolysaccharides (LPS), lipoteichoic acids (LTA), and peptidoglycans (PGN) are stimulatory ligands for Toll-like receptors (TLRs), and recent evidence indicates the important role of variation in TLR genes for permitting overgrowth of gram negative and BV-type flora. We assessed whether genetic polymorphisms in five TLR genes (TLR1, TLR2, TLR4, TLR6, and TLR9) could be determinants of differential host immune responses to BV in 159 HIV-1-positive African American adolescents enrolled in the Reaching for Excellence in Adolescent Care and Health (REACH) study. BV was assessed biannually and diagnosed either by a Nugent Score of at least 7 of 10, or using the Amsel Criteria. Cox-proportional hazards regression models, adjusted for concurrent Chlamydia and Gonorrhea infections, douching, and absolute CD4 cell count, were used to identify host genetic factors associated with BV. Two SNPs were associated with BV as diagnosed by the Nugent Score and the combined criteria: a minor allele G of rs4986790 (frequency=0.07), which encodes a His to Tyr substitution in TLR4 (HR=1.47, 95% CI 1.15–1.87) and rs187084 (frequency=0.24) on TLR9. The minor allele of rs1898830 (frequency=0.13) was associated with an increased hazard of BV defined by the Amsel criteria (HR=1.86, 95%CI 1.17–2.95). Further studies are warranted to confirm the associations of TLR gene variants and also to understand the underlying pathways and immunogenetic correlates in the context of HIV-1 infection.
HIV-1; Bacterial Vaginosis; Toll Like Receptors
In the United States (US), kidney dysfunction is prevalent in almost 30% of HIV-infected patients and is an independent predictor of mortality. Proteinuria and elevated serum cystatin C (eCysC) are used as markers of kidney disease in the general population; however, the prevalence of these markers in HIV-infected adolescents is largely unknown.
This study includes 304 HIV-infected adolescents from the Reaching for Excellence in Adolescent Care and Health (REACH) cohort, an observational study of adolescents recruited from 13 US cities. Clinical and demographic characteristics of participants were evaluated as correlates of proteinuria, a urine protein to creatinine ratio (UP/Cr) of ≥200 mg/g. Select univariate predictors were assessed to determine the association with urinary protein excretion and serum cystatin C in multivariable linear regression models and proteinuria and elevated serum cystatin C (eCysC ≥ 75th percentile) in multivariable logistic regression models.
Overall, 19.1% of the participants had proteinuria while 23.7% had an eCysC. Low CD4+ T-lymphocyte counts (<200 cells/mm3) were significantly associated with a greater UP/Cr in linear models and with proteinuria in logistic regression models. CD4+ T-lymphocyte counts <500 cells/mm3 were significantly associated with a greater serum cystatin C concentration in linear models and with eCysC in logistic regression models.
Proteinuria among HIV-infected adolescents in REACH was approximately two-fold greater than healthy US adolescents. Both proteinuria and eCysC are associated with CD4+ T-lymphocyte counts. Further studies investigating early markers of kidney disease and the association with immune status and inflammation in adolescents are needed.
HIV; adolescent; kidney; CD4+ T-cells; proteinuria; serum cystatin C
Both the prophylactic human papillomavirus (HPV) vaccines, Gardasil® and Cervarix®, are licensed for the prevention of cervical cancer in females, and Gardasil is also licensed for the prevention of genital warts and anal cancer in both males and females. This review focuses on the uptake of these vaccines in adolescent males and females in the USA and the barriers associated with vaccine initiation and completion. In the USA in 2009, approximately 44.3% of adolescent females aged 13–17 years had received at least one dose of the HPV vaccine, but only 26.7% had received all three doses. In general, the Northeast and Midwest regions of the USA have the highest rates of HPV vaccine initiation in adolescent females, while the Southeast has the lowest rates of vaccine initiation. Uptake of the first dose of the HPV vaccine in adolescent females did not vary by race/ethnicity; however, completion of all three doses is lower among African Americans (23.1%) and Latinos (23.4%) compared with Caucasians (29.3%). At present, vaccination rates among adolescent females are lower than expected, and thus vaccine models suggest that it is more cost-effective to vaccinate both adolescent males and females. Current guidelines for HPV vaccination in adolescent males is recommended only for “permissive use,” which leaves this population out of routine vaccination for HPV. The uptake of the vaccine is challenged by the high cost, feasibility, and logistics of three-dose deliveries. The biggest impact on acceptability of the vaccine is by adolescents, physicians, parents, and the community. Future efforts need to focus on HPV vaccine education among adolescents and decreasing the barriers associated with poor vaccine uptake and completion in adolescents before their sexual debut, but Papanicolau screening should remain routine among adults and those already infected until a therapeutic vaccine can be developed.
human papillomavirus; vaccine uptake; adolescent health
Populations of African ancestry continue to account for a disproportionate burden of human immunodeficiency virus type 1 (HIV-1) epidemic in the US. We investigated the effects of human leukocyte antigen (HLA) class I markers in association with virologic and immunologic control of HIV-1 infection among 338 HIV-1 subtype B-infected African Americans in two cohorts: REACH (Reaching for Excellence in Adolescent Care and Health) and HERS (HIV Epidemiology Research Study). One-year treatment-free interval measurements of HIV-1 RNA viral loads and CD4+ T-cells were examined both separately and combined to represent three categories of HIV-1 disease control (76 “controllers,” 169 “intermediates,” and 93 “non-controllers”). Certain previously or newly implicated HLA class I alleles (A*32, A*36, A*74, B*14, B*1510, B*3501, B*45, B*53, B*57, Cw*04, Cw*08, Cw*12, and Cw*18) were associated with one or more of the endpoints in univariate analyses. After multivariable adjustments for other genetic and non-genetic risk factors of HIV-1 progression, the subset of alleles more strongly or consistently associated with HIV-1 disease control included A*32, A*74, B*14, B*45, B*53, B*57, and Cw*08. Carriage of infrequent HLA-B but not HLA-A alleles was associated with more favorable disease outcomes. Certain HLA class I associations with control of HIV-1 infection span the boundaries of race and viral subtype; while others appear confined within one or the other of those boundaries.
HLA class I; Allele frequency; HIV-1 control; African American
Genetic variants in the inhibiting FcγRIIB mediate anti-inflammatory responses and influence IVIG refractoriness (IVIG-R). However, these variants are rare in Asian and Hispanic populations so other genes in the pathway could be potentially involved. IVIG is ineffective in mice lacking SIGN-R1, a related molecule to human DC-SIGN. Further, DC-SIGN is a known receptor for sialylated Fc, the component responsible for the anti-inflammatory action of IVIG. Thus, we hypothesized that DC-SIGN would also be involved in the pathway of IVIG response in Kawasaki Disease (KD) patients.
A case-control approach was performed to examine the differential distribution of five single nucleotide polymorphisms (SNPs) in DC-SIGN promoter with IVIG-R among White (158 vs. 62), Asian (64 vs. 12) and Hispanic (55 vs. 20) KD patients. Distinct differences in allele frequency distributions of several variants in the DC-SIGN promoter were observed in the three ethnic groups. Further, Asians with the major allele “A” in rs2287886 were more likely (OR = 1.76, p = 0.04) to be IVIG non-responder, but this allele is a minor allele in other two ethnic groups, where the association was not apparent.
DC-SIGN can potentially complement the role of FcγRIIB in the anti-inflammatory cascade involved in the IVIG response mechanism.
Kawasaki disease; IVIG treatment response; FcγR; Coronary artery disease; DC-SIGN
To replicate the association of variants in RYR3 gene with common carotid intima-media thickness (cIMT), a surrogate marker of atherosclerosis, we genotyped single nucleotide polymorphisms (SNPs) rs2229116 and rs7177922 in a sub-population of 244 HIV-positive and HIV-negative men. SNP rs2229116 was associated with common cIMT in HIV infected white men after adjusting for age and use of stavudine (d4T). The association was more evident at younger ages and decreased among older individuals.
Several lines of evidence have supported a host genetic contribution to vaccine response, but genome-wide assessments for specific determinants have been sparse. Here we describe a genome-wide association study (GWAS) of protective antigen-specific antibody (AbPA) responses among 726 European-Americans who received Anthrax Vaccine Adsorbed (AVA) as part of a clinical trial. After quality control, 736,996 SNPs were tested for association with the AbPA response to 3 or 4 AVA vaccinations given over a 6-month period. No SNP achieved the threshold of genome-wide significance (p=5x10−8), but suggestive associations (p<1x10−5) were observed for SNPs in or near the class II region of the major histocompatibility complex (MHC), in the promoter region of SPSB1, and adjacent to MEX3C. Multivariable regression modeling suggested that much of the association signal within the MHC corresponded to previously identified HLA DR-DQ haplotypes involving component HLA-DRB1 alleles of *15:01, *01:01, or *01:02. We estimated the proportion of additive genetic variance explained by common SNP variation for the AbPA response after the 6 month vaccination. This analysis indicated a significant, albeit imprecisely estimated, contribution of variation tagged by common polymorphisms (p=0.032). Future studies will be required to replicate these findings in European Americans and to further elucidate the host genetic factors underlying variable immune response to AVA.
Anthrax vaccines; Bacillus anthracis; bacterial vaccines; vaccination; Genome-wide association study
A functional polymorphism in the inhibitory IgG-Fc receptor FcγRIIB influences intravenous immunoglobulin (IVIG) response in Kawasaki Disease (KD) a vasculitis preferentially affecting the coronary arteries in children. We tested the hypothesis that the polymorphisms in the activating receptors (Fcγ RIIA, Fcγ RIIIA and Fcγ RIIIB) also influence susceptibility, IVIG treatment response, and coronary artery disease (CAD) in KD patients.
Methods and Results
We genotyped polymorphisms in the activating FcγRIIA, FcγRIIIA and FcγRIIIB genes using pyrosequencing in 443 KD patients, including 266 trios and 150 single parent-child pairs, in northwest US and genetically determined race with 155 ancestry information markers. We used the FBAT program to test for transmission disequilibrium and further generated pseudo-sibling controls for comparisons to the cases. The FcγRIIA-131H variant showed an association with KD (p = 0.001) with ORadditive = 1.51 [1.16–1.96], p = 0.002) for the primary combined population, which persisted in both Caucasian (p = .04) and Asian (p = .01) subgroups and is consistent with the recent genome-wide association study. We also identified over-transmission of FcγRIIIB-NA1 among IVIG non-responders (p = 0.0002), and specifically to Caucasian IVIG non-responders (p = 0.007). Odds ratios for overall and Caucasian non-responders were respectively 3.67 [1.75–7.66], p = 0.0006 and 3.60 [1.34–9.70], p = 0.01. Excess NA1 transmission also occurred to KD with CAD (ORadditive = 2.13 [1.11–4.0], p = 0.02).
A common variation in FcγRIIA is associated with increased KD susceptibility. The FcγRIIIB-NA1, which confers higher affinity for IgG compared to NA2, is a determining factor for treatment response. These activating FcγRs play an important role in KD pathogenesis and mechanism of IVIG anti-inflammatory.
coronary disease; pediatrics; Kawasaki disease; IVIG treatment response; FcγR
Anthrax and its etiologic agent remain a biological threat. Anthrax vaccine is highly effective, but vaccine-induced IgG antibody responses vary widely following required doses of vaccinations. Such variation can be related to genetic factors, especially genomic copy number variants (CNVs) that are known to be enriched among genes with immunologic function. We have tested this hypothesis in two study populations from a clinical trial of anthrax vaccination.
We performed CNV-based genome-wide association analyses separately on 794 European Americans and 200 African-Americans. Antibodies to protective antigen were measured at week 8 (early response) and week 30 (peak response) using an enzyme-linked immunosorbent assay. We used DNA microarray data (Affymetrix 6.0) and two CNV detection algorithms, hidden markov model (PennCNV) and circular binary segmentation (GeneSpring) to determine CNVs in all individuals. Multivariable regression analyses were used to identify CNV-specific associations after adjusting for relevant non-genetic covariates.
Within the 22 autosomal chromosomes, 2,943 non-overlapping CNV regions were detected by both algorithms. Genomic insertions containing HLA-DRB5, DRB1 and DQA1/DRA genes in the major histocompatibility complex (MHC) region (chromosome 6p21.3) were moderately associated with elevated early antibody response (β = 0.14, p = 1.78×10−3) among European Americans, and the strongest association was observed between peak antibody response and a segmental insertion on chromosome 1, containing NBPF4, NBPF5, STXMP3, CLCC1, and GPSM2 genes (β = 1.66, p = 6.06×10−5). For African-Americans, segmental deletions spanning PRR20, PCDH17 and PCH68 genes on chromosome 13 were associated with elevated early antibody production (β = 0.18, p = 4.47×10−5). Population-specific findings aside, one genomic insertion on chromosome 17 (containing NSF, ARL17 and LRRC37A genes) was associated with elevated peak antibody response in both populations.
Multiple CNV regions, including the one consisting of MHC genes that is consistent with earlier research, can be important to humoral immune responses to anthrax vaccine adsorbed.
Several studies have assessed risk factors associated with herpes simplex virus-2 (HSV-2) prevalence in adults; however, few have focused on HSV-2 incidence, particularly in adolescents. The objective of this study was to determine HSV-2 prevalence and incidence and associated risk factors in a HIV-1-positive and at risk HIV-1-negative adolescent population.
Sera were tested for HSV-2 antibodies in 518 adolescents in the Reaching for Excellence in Adolescent Care and Health (REACH) cohort at baseline and again at the final follow-up visit. Prevalence at baseline and incidence (per person years) of HSV-2 infection were calculated. Furthermore, among HIV-1-positive individuals, a subgroup analysis was performed to assess risk factors for HSV-2 infection. Conditional logistic regression was used to estimate odds ratios (OR) and p-values (p) for associations between CD4+ T-cell (CD4+) count, HIV-1 viral load (VL), and HSV-2 acquisition, adjusting for antiretroviral therapy use, other sexually transmitted infections, gender, race, and number of sexual partners.
At baseline, 179 (35%) subjects were HSV-2 positive, with an additional 47 (16%) new cases being identified during a median follow-up time of 1.95 years and an incidence rate of 7.35 cases per 100 person years (py). Several risk factors were associated with HSV-2 prevalence (being female, non-Hispanic, uncertainty of sexual preference, and HIV-1 positive) and incidence (using drugs, alcohol, and number of new sexual partners). Among HIV-1 positives, an increase in CD4+ count by 50 cell/mm3 (OR, 1.17; 95% CI 1.04–1.31, p=0.008) was associated with HSV-2 acquisition.
The high prevalence and incidence of HSV-2 infection among adolescents, compared to the general population at this age group suggests a critical need for screening and preventive programs among this targeted group.
HIV-1; HSV-2; CD4+ count; adolescents
Several CC-motif chemokine ligands (CCLs) can block HIV-1 binding sites on CC-motif chemokine receptor 5 (CCR5) and inhibit viral entry. We studied single nucleotide polymorphisms (SNPs) in genes encoding three CCR5 ligands [CCL3 (MIP-1α), CCL4 (MIP-1β), and CCL5 (RANTES)] along with an adjacent gene encoding a CCR2 ligand [CCL2 (MCP-1)] to identify candidate markers for HIV-1 infection and pathogenesis. Analyses of 567 HIV-1 serodiscordant Zambian couples revealed that rs5029410C (in CCL3 intron 2) was associated with lower viral load (VL) in seroconverters, adjusted for gender and age (regression β=−0.57 log10, P=4×10−6). In addition, rs34171309A in CCL3 exon 3 was associated with increased risk of HIV-1 acquisition in exposed seronegatives (hazard ratio=1.52, P=0.006 when adjusted for donor VL and genital ulcer/inflammation). The CCL3 exon 3 SNP, encoding a conservative Glu-to-Asp substitution, and five neighboring SNPs in tight linkage disequilibrium all showed similar associations with HIV-1 acquisition. How these multiple CCL3 SNPs may alter the occurrence or course of HIV-1 infection remains to be determined.
HIV-1 transmission; CCL2; CCL3; CCL4; CCL5; SNP
Genome-wide association (GWA) studies have become a standard approach for discovering and validating genomic polymorphisms putatively associated with phenotypes of interest. Accounting for population structure in GWA studies is critical to attain unbiased parameter measurements and control Type I error. One common approach to accounting for population structure is to include several principal components derived from the entire autosomal dataset, which reflects population structure signal. However, knowing which components to include is subjective and generally not conclusive. We examined how phylogenetic signal from mitochondrial DNA (mtDNA) and chromosome Y (chr:Y) markers is concordant with principal component data based on autosomal markers to determine whether mtDNA and chr:Y phylogenetic data can help guide principal component selection. Using HAPMAP and other original data from individuals of multiple ancestries, we examined the relationships of mtDNA and chr:Y phylogenetic signal with the autosomal PCA using best subset logistic regression. We show that while the two approaches agree at times, this is independent of the component order and not completely represented in the Eigen values. Additionally, we use simulations to demonstrate that our approach leads to a slightly reduced Type I error rate compared to the standard approach. This approach provides preliminary evidence to support the theoretical concept that mtDNA and chr:Y data can be informative in locating the PCs that are most associated with evolutionary history of populations that are being studied, although the utility of such information will depend on the specific situation.
phylogeny; PCA; Y chromosome; mitochondria; population sub-structure
Genetic studies may help explain abnormalities of fat distribution in HIV-infected patients treated with antiretroviral therapy (ARV).
Subcutaneous adipose tissue (SAT) volume measured by magnetic resonance imaging (MRI) in leg, lower trunk, upper trunk, and arm was examined in 192 HIV-infected Caucasian men, ARV-treated from the Fat Redistribution and Metabolic Change in HIV infection (FRAM) study. Single nucleotide polymorphisms (SNPs) were assayed using the Illumina HumanCNV370-quad beadchip. Multivariate and univariate genome wide association analyses of the four SAT depots were implemented in PLINK software adjusted for age and ARV duration. Functional annotation analysis (FAA) using Ingenuity Systems Pathway Analysis tool (IPA) was carried out for markers with P<10-3 near known genes identified by multivariate analysis.
Loci (rs10504906, rs13267998, rs921231) in or near the anion exchanger solute carrier family 26, member 7 isoform a (SLC26A7) were strongly associated with upper trunk and arm SAT (9.8*10-7≤P<7.8*10-6). Loci (rs193139, rs7523050, rs1761621) in and near a gene rich region including G-protein-signaling modulator 2 (GPSM2) and syntaxin binding protein 3 (STXBP3) were significantly associated with lower body SAT depots (9.9*10-7≤P<9.5*10-6). GPSM2 is associated with cell division and cancer while STXBP3 is associated with glucose metabolism in adipoctyes. IPA identified atherosclerosis, mitochondrial function and T-Cell mediated apoptosis as processes related to SAT volume in HIV-infected individuals (P<5*10-3).
Our results are limited by the small sample size and replication is needed, however this genomic scan uncovered new genes associated with metabolism and inflammatory pathways that may affect SAT volume in ARV-treated HIV-infected patients.
HIV; HAART; GWAS; Subcutaneous Fat; SAT