Exploiting meta-analysis of genome-wide interaction with serum 25-hydroxyvitamin D to identify novel genetic loci associated with pulmonary function

BACKGROUND

Higher 25-hydroxyvitamin D (25(OH)D) concentrations in serum has a positive association with pulmonary function. Investigating genome-wide interactions with 25(OH)D may reveal new biological insights into pulmonary function.

OBJECTIVES

We aimed to identify novel genetic variants associated with pulmonary function by accounting for 25(OH)D interactions.

METHODS

We included 211,264 participants from the observational United Kingdom Biobank study with pulmonary function tests (PFTs), genome-wide genotypes, and 25(OH)D concentrations from 4 ancestral backgrounds-European, African, East Asian, and South Asian. Among PFTs, we focused on forced expiratory volume in the first second (FEV1) and forced vital capacity (FVC) because both were previously associated with 25(OH)D. We performed genome-wide association study (GWAS) analyses that accounted for variant×25(OH)D interaction using the joint 2 degree-of-freedom (2df) method, stratified by participants' smoking history and ancestry, and meta-analyzed results. We evaluated interaction effects to determine how variant-PFT associations were modified by 25(OH)D concentrations and conducted pathway enrichment analysis to examine the biological relevance of our findings.

RESULTS

Our GWAS meta-analyses, accounting for interaction with 25(OH)D, revealed 30 genetic variants significantly associated with FEV1 or FVC (P2df <5.00×10-8) that were not previously reported for PFT-related traits. These novel variant signals were enriched in lung function-relevant pathways, including the p38 MAPK pathway. Among variants with genome-wide-significant 2df results, smoking-stratified meta-analyses identified 5 variants with 25(OH)D interactions that influenced FEV1 in both smoking groups (never smokers P1df interaction<2.65×10-4; ever smokers P1df interaction<1.71×10-5); rs3130553, rs2894186, rs79277477, and rs3130929 associations were only evident in never smokers, and the rs4678408 association was only found in ever smokers.

CONCLUSION

Genetic variant associations with lung function can be modified by 25(OH)D, and smoking history can further modify variant×25(OH)D interactions. These results expand the known genetic architecture of pulmonary function and add evidence that gene-environment interactions, including with 25(OH)D and smoking, influence lung function.

Evaluating 17 methods incorporating biological function with GWAS summary statistics to accelerate discovery demonstrates a tradeoff between high sensitivity and high positive predictive value

Where sufficiently large genome-wide association study (GWAS) samples are not currently available or feasible, methods that leverage increasing knowledge of the biological function of variants may illuminate discoveries without increasing sample size. We comprehensively evaluated 17 functional weighting methods for identifying novel associations. We assessed the performance of these methods using published results from multiple GWAS waves across each of five complex traits. Although no method achieved both high sensitivity and positive predictive value (PPV) for any trait, a subset of methods utilizing pleiotropy and expression quantitative trait loci nominated variants with high PPV (>75%) for multiple traits. Application of functionally weighting methods to enhance GWAS power for locus discovery is unlikely to circumvent the need for larger sample sizes in truly underpowered GWAS, but these results suggest that applying functional weighting to GWAS can accurately nominate additional novel loci from available samples for follow-up studies.

Change in plasma α-tocopherol associations with attenuated pulmonary function decline and with CYP4F2 missense variation

BACKGROUND

Vitamin E (vitE) is hypothesized to attenuate age-related decline in pulmonary function.

OBJECTIVES

We investigated the association between change in plasma vitE (∆vitE) and pulmonary function decline [forced expiratory volume in the first second (FEV1)] and examined genetic and nongenetic factors associated with ∆vitE.

METHODS

We studied 1144 men randomly assigned to vitE in SELECT (Selenium and Vitamin E Cancer Prevention Trial). ∆vitE was the difference between baseline and year 3 vitE concentrations measured with GC-MS. FEV1 was measured longitudinally by spirometry. We genotyped 555 men (vitE-only arm) using the Illumina Expanded Multi-Ethnic Genotyping Array (MEGAex). We used mixed-effects linear regression modeling to examine the ∆vitE-FEV1 association.

RESULTS

Higher ∆vitE was associated with lower baseline α-tocopherol (α-TOH), higher baseline γ-tocopherol, higher baseline free cholesterol, European ancestry (as opposed to African) (all P < 0.05), and the minor allele of a missense variant in cytochrome P450 family 4 subfamily F member 2 (CYP4F2) (rs2108622-T; 2.4 µmol/L higher ∆vitE, SE: 0.8 µmol/L; P = 0.0032). Higher ∆vitE was associated with attenuated FEV1 decline, with stronger effects in adherent participants (≥80% of supplements consumed): a statistically significant ∆vitE × time interaction (P = 0.014) indicated that a 1-unit increase in ∆vitE was associated with a 2.2-mL/y attenuation in FEV1 decline (SE: 0.9 mL/y). The effect size for 1 SD higher ∆vitE (+4 µmol/mmol free-cholesterol-adjusted α-TOH) was roughly one-quarter of the effect of 1 y of aging, but in the opposite direction. The ∆vitE-FEV1 association was similar in never smokers (2.4-mL/y attenuated FEV1 decline, SE: 1.0 mL/y; P = 0.017, n = 364), and current smokers (2.8-mL/y, SE: 1.6 mL/y; P = 0.079, n = 214), but there was little to no effect in former smokers (-0.64-mL/y, SE: 0.9 mL/y; P = 0.45, n = 564).

CONCLUSIONS

Greater response to vitE supplementation was associated with attenuated FEV1 decline. The response to supplementation differed by rs2108622 such that individuals with the C allele, compared with the T allele, may need a higher dietary intake to reach the same plasma vitE concentration.

Integration of evidence across human and model organism studies: A meeting report

The National Institute on Drug Abuse and Joint Institute for Biological Sciences at the Oak Ridge National Laboratory hosted a meeting attended by a diverse group of scientists with expertise in substance use disorders (SUDs), computational biology, and FAIR (Findability, Accessibility, Interoperability, and Reusability) data sharing. The meeting's objective was to discuss and evaluate better strategies to integrate genetic, epigenetic, and 'omics data across human and model organisms to achieve deeper mechanistic insight into SUDs. Specific topics were to (a) evaluate the current state of substance use genetics and genomics research and fundamental gaps, (b) identify opportunities and challenges of integration and sharing across species and data types, (c) identify current tools and resources for integration of genetic, epigenetic, and phenotypic data, (d) discuss steps and impediment related to data integration, and (e) outline future steps to support more effective collaboration-particularly between animal model research communities and human genetics and clinical research teams. This review summarizes key facets of this catalytic discussion with a focus on new opportunities and gaps in resources and knowledge on SUDs.

Expanding the genetic architecture of nicotine dependence and its shared genetics with multiple traits

Cigarette smoking is the leading cause of preventable morbidity and mortality. Genetic variation contributes to initiation, regular smoking, nicotine dependence, and cessation. We present a Fagerström Test for Nicotine Dependence (FTND)-based genome-wide association study in 58,000 European or African ancestry smokers. We observe five genome-wide significant loci, including previously unreported loci MAGI2/GNAI1 (rs2714700) and TENM2 (rs1862416), and extend loci reported for other smoking traits to nicotine dependence. Using the heaviness of smoking index from UK Biobank (N = 33,791), rs2714700 is consistently associated; rs1862416 is not associated, likely reflecting nicotine dependence features not captured by the heaviness of smoking index. Both variants influence nearby gene expression (rs2714700/MAGI2-AS3 in hippocampus; rs1862416/TENM2 in lung), and expression of genes spanning nicotine dependence-associated variants is enriched in cerebellum. Nicotine dependence (SNP-based heritability = 8.6%) is genetically correlated with 18 other smoking traits (rg = 0.40-1.09) and co-morbidities. Our results highlight nicotine dependence-specific loci, emphasizing the FTND as a composite phenotype that expands genetic knowledge of smoking.

Genome-Wide Meta-Analyses of FTND and TTFC Phenotypes

INTRODUCTION

FTND (Fagerstrӧm test for nicotine dependence) and TTFC (time to smoke first cigarette in the morning) are common measures of nicotine dependence (ND). However, genome-wide meta-analysis for these phenotypes has not been reported.

METHODS

Genome-wide meta-analyses for FTND (N = 19,431) and TTFC (N = 18,567) phenotypes were conducted for adult smokers of European ancestry from 14 independent cohorts.

RESULTS

We found that SORBS2 on 4q35 (p = 4.05 × 10-8), BG182718 on 11q22 (p = 1.02 × 10-8), and AA333164 on 14q21 (p = 4.11 × 10-9) were associated with TTFC phenotype. We attempted replication of leading candidates with independent samples (FTND, N = 7010 and TTFC, N = 10 061), however, due to limited power of the replication samples, the replication of these new loci did not reach significance. In gene-based analyses, COPB2 was found associated with FTND phenotype, and TFCP2L1, RELN, and INO80C were associated with TTFC phenotype. In pathway and network analyses, we found that the interconnected interactions among the endocytosis, regulation of actin cytoskeleton, axon guidance, MAPK signaling, and chemokine signaling pathways were involved in ND.

CONCLUSIONS

Our analyses identified several promising candidates for both FTND and TTFC phenotypes, and further verification of these candidates was necessary. Candidates supported by both FTND and TTFC (CHRNA4, THSD7B, RBFOX1, and ZNF804A) were associated with addiction to alcohol, cocaine, and heroin, and were associated with autism and schizophrenia. We also identified novel pathways involved in cigarette smoking. The pathway interactions highlighted the importance of receptor recycling and internalization in ND.

IMPLICATIONS

Understanding the genetic architecture of cigarette smoking and ND is critical to develop effective prevention and treatment. Our study identified novel candidates and biological pathways involved in FTND and TTFC phenotypes, and this will facilitate further investigation of these candidates and pathways.

Omega-3 Fatty Acids and Genome-Wide Interaction Analyses Reveal DPP10-Pulmonary Function Association

RATIONALE

Omega-3 polyunsaturated fatty acids (n-3 PUFAs) have anti-inflammatory properties that could benefit adults with comprised pulmonary health.

OBJECTIVE

To investigate n-3 PUFA associations with spirometric measures of pulmonary function tests (PFTs) and determine underlying genetic susceptibility.

METHODS

Associations of n-3 PUFA biomarkers (α-linolenic acid, eicosapentaenoic acid, docosapentaenoic acid [DPA], and docosahexaenoic acid [DHA]) were evaluated with PFTs (FEV1, FVC, and FEV1/FVC) in meta-analyses across seven cohorts from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium (N = 16,134 of European or African ancestry). PFT-associated n-3 PUFAs were carried forward to genome-wide interaction analyses in the four largest cohorts (N = 11,962) and replicated in one cohort (N = 1,687). Cohort-specific results were combined using joint 2 degree-of-freedom (2df) meta-analyses of SNP associations and their interactions with n-3 PUFAs.

RESULTS

DPA and DHA were positively associated with FEV1 and FVC (P < 0.025), with evidence for effect modification by smoking and by sex. Genome-wide analyses identified a novel association of rs11693320-an intronic DPP10 SNP-with FVC when incorporating an interaction with DHA, and the finding was replicated (P2df = 9.4 × 10-9 across discovery and replication cohorts). The rs11693320-A allele (frequency, ∼80%) was associated with lower FVC (PSNP = 2.1 × 10-9; βSNP = -161.0 ml), and the association was attenuated by higher DHA levels (PSNP×DHA interaction = 2.1 × 10-7; βSNP×DHA interaction = 36.2 ml).

CONCLUSIONS

We corroborated beneficial effects of n-3 PUFAs on pulmonary function. By modeling genome-wide n-3 PUFA interactions, we identified a novel DPP10 SNP association with FVC that was not detectable in much larger studies ignoring this interaction.

Genome-wide association study across European and African American ancestries identifies a SNP in DNMT3B contributing to nicotine dependence

Cigarette smoking is a leading cause of preventable mortality worldwide. Nicotine dependence, which reduces the likelihood of quitting smoking, is a heritable trait with firmly established associations with sequence variants in nicotine acetylcholine receptor genes and at other loci. To search for additional loci, we conducted a genome-wide association study (GWAS) meta-analysis of nicotine dependence, totaling 38,602 smokers (28,677 Europeans/European Americans and 9925 African Americans) across 15 studies. In this largest-ever GWAS meta-analysis for nicotine dependence and the largest-ever cross-ancestry GWAS meta-analysis for any smoking phenotype, we reconfirmed the well-known CHRNA5-CHRNA3-CHRNB4 genes and further yielded a novel association in the DNA methyltransferase gene DNMT3B. The intronic DNMT3B rs910083-C allele (frequency=44-77%) was associated with increased risk of nicotine dependence at P=3.7 × 10-8 (odds ratio (OR)=1.06 and 95% confidence interval (CI)=1.04-1.07 for severe vs mild dependence). The association was independently confirmed in the UK Biobank (N=48,931) using heavy vs never smoking as a proxy phenotype (P=3.6 × 10-4, OR=1.05, and 95% CI=1.02-1.08). Rs910083-C is also associated with increased risk of squamous cell lung carcinoma in the International Lung Cancer Consortium (N=60,586, meta-analysis P=0.0095, OR=1.05, and 95% CI=1.01-1.09). Moreover, rs910083-C was implicated as a cis-methylation quantitative trait locus (QTL) variant associated with higher DNMT3B methylation in fetal brain (N=166, P=2.3 × 10-26) and a cis-expression QTL variant associated with higher DNMT3B expression in adult cerebellum from the Genotype-Tissue Expression project (N=103, P=3.0 × 10-6) and the independent Brain eQTL Almanac (N=134, P=0.028). This novel DNMT3B cis-acting QTL variant highlights the importance of genetically influenced regulation in brain on the risks of nicotine dependence, heavy smoking and consequent lung cancer.

Large-scale association analysis identifies new lung cancer susceptibility loci and heterogeneity in genetic susceptibility across histological subtypes

Although several lung cancer susceptibility loci have been identified, much of the heritability for lung cancer remains unexplained. Here 14,803 cases and 12,262 controls of European descent were genotyped on the OncoArray and combined with existing data for an aggregated genome-wide association study (GWAS) analysis of lung cancer in 29,266 cases and 56,450 controls. We identified 18 susceptibility loci achieving genome-wide significance, including 10 new loci. The new loci highlight the striking heterogeneity in genetic susceptibility across the histological subtypes of lung cancer, with four loci associated with lung cancer overall and six loci associated with lung adenocarcinoma. Gene expression quantitative trait locus (eQTL) analysis in 1,425 normal lung tissue samples highlights RNASET2, SECISBP2L and NRG1 as candidate genes. Other loci include genes such as a cholinergic nicotinic receptor, CHRNA2, and the telomere-related genes OFBC1 and RTEL1. Further exploration of the target genes will continue to provide new insights into the etiology of lung cancer.

KAT2B polymorphism identified for drug abuse in African Americans with regulatory links to drug abuse pathways in human prefrontal cortex

Drug abuse is a common and heritable set of disorders, but the underlying genetic factors are largely unknown. We conducted genome-wide association studies of drug abuse using 7 million imputed single nucleotide polymorphisms (SNPs) and insertions/deletions in African Americans (AAs; n = 3742) and European Americans (EAs; n = 6845). Cases were drawn from the Urban Health Study of street-recruited people, who injected drugs and reported abusing opioids, cocaine, marijuana, stimulants and/or other drugs 10 or more times in the past 30 days, and were compared with population controls. Independent replication testing was conducted in 755 AAs and 1131 EAs from the Genetic Association Information Network. An intronic SNP (rs9829896) in the K(lysine) acetyltransferase 2B (KAT2B) gene was significantly associated with drug abuse in AAs (P = 4.63 × 10^-8 ) and independently replicated in AAs (P = 0.0019). The rs9829896-C allele (frequency = 12%) had odds ratios of 0.68 and 0.53 across the AA cohorts: meta-analysis P = 3.93 × 10^-10 . Rs9829896-C was not associated with drug abuse across the EA cohorts: frequency = 36% and meta-analysis P = 0.12. Using dorsolateral prefrontal cortex data from the BrainCloud cohort, we found that rs9829896-C was associated with reduced KAT2B expression in AAs (n = 113, P = 0.050) but not EAs (n = 110, P = 0.39). KAT2B encodes a transcriptional regulator in the cyclic adenosine monophosphate and dopamine signaling pathways, and rs9829896-C was associated with expression of genes in these pathways: reduced CREBBP expression (P = 0.011) and increased OPRM1 expression (P = 0.016), both in AAs only. Our study identified the KAT2B SNP rs9829896 as having novel and biologically plausible associations with drug abuse and gene expression in AAs but not EAs, suggesting ancestry-specific effects.

Replication of ZNF804A gene variant associations with risk of heroin addiction

Heroin addiction is heritable, but few specific genetic variants have been reproducibly associated with this disease. The zinc finger protein 804A (ZNF804A) gene is a biologically plausible susceptibility gene for heroin addiction, given its function as a transcription factor in human brain. Novel associations of two common ZNF804A single nucleotide polymorphisms (SNPs), rs7597593 and rs1344706, with heroin addiction have been reported in Han Chinese. Both SNPs have also been implicated for regulating ZNF804A expression in human brain, including the addiction-relevant dorsolateral prefrontal cortex. In this independent replication study, we tested the rs7597593 and rs1344706 SNP genotypes and their corresponding haplotypes for association with heroin addiction using cases drawn from the Urban Health Study and population controls: total N = 10 757 [7095 European Americans (EAs) and 3662 African Americans (AAs)]. We independently replicated both ZNF804A SNP associations in EAs: the rs7597593-T (P = 0.016) and rs1344706-A (P = 0.029) alleles both being associated with increased risk of heroin addiction, consistent with the prior report. Neither SNP was associated in AAs alone, but meta-analysis across both ancestry groups resulted in significant associations for rs1344706-A [P = 0.016, odds ratio (95% confidence interval) = 1.13 (1.02-1.25)] and its haplotype with rs7597593-T [P = 0.0067, odds ratio (95% confidence interval) = 1.16 (1.04-1.29)]. By showing consistent associations across independent studies and diverse ancestry groups, our study provides evidence that these two ZNF804A SNPs and their risk haplotype are among the few replicable genetic associations with heroin addiction.

Genome-wide meta-analysis reveals common splice site acceptor variant in CHRNA4 associated with nicotine dependence

We conducted a 1000 Genomes-imputed genome-wide association study (GWAS) meta-analysis for nicotine dependence, defined by the Fagerström Test for Nicotine Dependence in 17 074 ever smokers from five European-ancestry samples. We followed up novel variants in 7469 ever smokers from five independent European-ancestry samples. We identified genome-wide significant association in the alpha-4 nicotinic receptor subunit (CHRNA4) gene on chromosome 20q13: lowest P=8.0 × 10(-9) across all the samples for rs2273500-C (frequency=0.15; odds ratio=1.12 and 95% confidence interval=1.08-1.17 for severe vs mild dependence). rs2273500-C, a splice site acceptor variant resulting in an alternate CHRNA4 transcript predicted to be targeted for nonsense-mediated decay, was associated with decreased CHRNA4 expression in physiologically normal human brains (lowest P=7.3 × 10(-4)). Importantly, rs2273500-C was associated with increased lung cancer risk (N=28 998, odds ratio=1.06 and 95% confidence interval=1.00-1.12), likely through its effect on smoking, as rs2273500-C was no longer associated with lung cancer after adjustment for smoking. Using criteria for smoking behavior that encompass more than the single 'cigarettes per day' item, we identified a common CHRNA4 variant with important regulatory properties that contributes to nicotine dependence and smoking-related consequences.

Cis-Expression Quantitative Trait Loci Mapping Reveals Replicable Associations with Heroin Addiction in OPRM1

BACKGROUND

No opioid receptor, mu 1 (OPRM1) gene polymorphisms, including the functional single nucleotide polymorphism (SNP) rs1799971, have been conclusively associated with heroin/other opioid addiction, despite their biological plausibility. We used evidence of polymorphisms altering OPRM1 expression in normal human brain tissue to nominate and then test associations with heroin addiction.

METHODS

We tested 103 OPRM1 SNPs for association with OPRM1 messenger RNA expression in prefrontal cortex from 224 European Americans and African Americans of the BrainCloud cohort. We then tested the 16 putative cis-expression quantitative trait loci (cis-eQTL) SNPs for association with heroin addiction in the Urban Health Study and two replication cohorts, totaling 16,729 European Americans, African Americans, and Australians of European ancestry.

RESULTS

Four putative cis-eQTL SNPs were significantly associated with heroin addiction in the Urban Health Study (smallest p = 8.9 × 10(-5)): rs9478495, rs3778150, rs9384169, and rs562859. Rs3778150, located in OPRM1 intron 1, was significantly replicated (p = 6.3 × 10(-5)). Meta-analysis across all case-control cohorts resulted in p = 4.3 × 10(-8): the rs3778150-C allele (frequency = 16%-19%) being associated with increased heroin addiction risk. Importantly, the functional SNP allele rs1799971-A was associated with heroin addiction only in the presence of rs3778150-C (p = 1.48 × 10(-6) for rs1799971-A/rs3778150-C and p = .79 for rs1799971-A/rs3778150-T haplotypes). Lastly, replication was observed for six other intron 1 SNPs that had prior suggestive associations with heroin addiction (smallest p = 2.7 × 10(-8) for rs3823010).

CONCLUSIONS

Our findings show that common OPRM1 intron 1 SNPs have replicable associations with heroin addiction. The haplotype structure of rs3778150 and nearby SNPs may underlie the inconsistent associations between rs1799971 and heroin addiction.

A multiancestry study identifies novel genetic associations with CHRNA5 methylation in human brain and risk of nicotine dependence

Nicotine dependence is influenced by chromosome 15q25.1 single nucleotide polymorphisms (SNPs), including the missense SNP rs16969968 that alters function of the α5 nicotinic acetylcholine receptor (CHRNA5) and noncoding SNPs that regulate CHRNA5 mRNA expression. We tested for cis-methylation quantitative trait loci (cis-meQTLs) using SNP genotypes and DNA methylation levels measured across the IREB2-HYKK-PSMA4-CHRNA5-CHRNA3-CHRNB4 genes on chromosome 15q25.1 in the BrainCloud and Brain QTL cohorts [total N = 175 European-Americans and 65 African-Americans (AAs)]. We identified eight SNPs that were significantly associated with CHRNA5 methylation in prefrontal cortex: P ranging from 6.0 × 10(-10) to 5.6 × 10(-5). These SNP-methylation associations were also significant in frontal cortex, temporal cortex and pons: P ranging from 4.8 × 10(-12) to 3.4 × 10(-3). Of the eight cis-meQTL SNPs, only the intronic CHRNB4 SNP rs11636753 was associated with CHRNA5 methylation independently of the known SNP effects in prefrontal cortex, and it was the most significantly associated SNP with nicotine dependence across five independent cohorts (total N = 7858 European ancestry and 3238 AA participants): P = 6.7 × 10(-4), odds ratio (OR) [95% confidence interval (CI)] = 1.11 (1.05-1.18). The rs11636753 major allele (G) was associated with lower CHRNA5 DNA methylation, lower CHRNA5 mRNA expression and increased nicotine dependence risk. Haplotype analyses showed that rs11636753-G and the functional rs16969968-A alleles together increased risk of nicotine dependence more than each variant alone: P = 3.1 × 10(-12), OR (95% CI) = 1.32 (1.22-1.43). Our findings identify a novel regulatory SNP association with nicotine dependence and connect, for the first time, previously observed differences in CHRNA5 mRNA expression and nicotine dependence risk to underlying DNA methylation differences.

Novel genetic locus implicated for HIV-1 acquisition with putative regulatory links to HIV replication and infectivity: a genome-wide association study

Fifty percent of variability in HIV-1 susceptibility is attributable to host genetics. Thus identifying genetic associations is essential to understanding pathogenesis of HIV-1 and important for targeting drug development. To date, however, CCR5 remains the only gene conclusively associated with HIV acquisition. To identify novel host genetic determinants of HIV-1 acquisition, we conducted a genome-wide association study among a high-risk sample of 3,136 injection drug users (IDUs) from the Urban Health Study (UHS). In addition to being IDUs, HIV- controls were frequency-matched to cases on environmental exposures to enhance detection of genetic effects. We tested independent replication in the Women’s Interagency HIV Study (N=2,533). We also examined publicly available gene expression data to link SNPs associated with HIV acquisition to known mechanisms affecting HIV replication/infectivity. Analysis of the UHS nominated eight genetic regions for replication testing. SNP rs4878712 in FRMPD1 met multiple testing correction for independent replication (P=1.38x10^-4), although the UHS-WIHS meta-analysis p-value did not reach genome-wide significance (P=4.47x10^-7 vs. P<5.0x10^-8) Gene expression analyses provided promising biological support for the protective G allele at rs4878712 lowering risk of HIV: (1) the G allele was associated with reduced expression of FBXO10 (r=-0.49, P=6.9x10^-5); (2) FBXO10 is a component of the Skp1-Cul1-F-box protein E3 ubiquitin ligase complex that targets Bcl-2 protein for degradation; (3) lower FBXO10 expression was associated with higher BCL2 expression (r=-0.49, P=8x10^-5); (4) higher basal levels of Bcl-2 are known to reduce HIV replication and infectivity in human and animal in vitro studies. These results suggest new potential biological pathways by which host genetics affect susceptibility to HIV upon exposure for follow-up in subsequent studies.

Imputation across genotyping arrays for genome-wide association studies: assessment of bias and a correction strategy

A great promise of publicly sharing genome-wide association data is the potential to create composite sets of controls. However, studies often use different genotyping arrays, and imputation to a common set of SNPs has shown substantial bias: a problem which has no broadly applicable solution. Based on the idea that using differing genotyped SNP sets as inputs creates differential imputation errors and thus bias in the composite set of controls, we examined the degree to which each of the following occurs: (1) imputation based on the union of genotyped SNPs (i.e., SNPs available on one or more arrays) results in bias, as evidenced by spurious associations (type 1 error) between imputed genotypes and arbitrarily assigned case/control status; (2) imputation based on the intersection of genotyped SNPs (i.e., SNPs available on all arrays) does not evidence such bias; and (3) imputation quality varies by the size of the intersection of genotyped SNP sets. Imputations were conducted in European Americans and African Americans with reference to HapMap phase II and III data. Imputation based on the union of genotyped SNPs across the Illumina 1M and 550v3 arrays showed spurious associations for 0.2 % of SNPs: ~2,000 false positives per million SNPs imputed. Biases remained problematic for very similar arrays (550v1 vs. 550v3) and were substantial for dissimilar arrays (Illumina 1M vs. Affymetrix 6.0). In all instances, imputing based on the intersection of genotyped SNPs (as few as 30 % of the total SNPs genotyped) eliminated such bias while still achieving good imputation quality.

Assessment of genotype imputation performance using 1000 Genomes in African American studies

Genotype imputation, used in genome-wide association studies to expand coverage of single nucleotide polymorphisms (SNPs), has performed poorly in African Americans compared to less admixed populations. Overall, imputation has typically relied on HapMap reference haplotype panels from Africans (YRI), European Americans (CEU), and Asians (CHB/JPT). The 1000 Genomes project offers a wider range of reference populations, such as African Americans (ASW), but their imputation performance has had limited evaluation. Using 595 African Americans genotyped on Illumina's HumanHap550v3 BeadChip, we compared imputation results from four software programs (IMPUTE2, BEAGLE, MaCH, and MaCH-Admix) and three reference panels consisting of different combinations of 1000 Genomes populations (February 2012 release): (1) 3 specifically selected populations (YRI, CEU, and ASW); (2) 8 populations of diverse African (AFR) or European (AFR) descent; and (3) all 14 available populations (ALL). Based on chromosome 22, we calculated three performance metrics: (1) concordance (percentage of masked genotyped SNPs with imputed and true genotype agreement); (2) imputation quality score (IQS; concordance adjusted for chance agreement, which is particularly informative for low minor allele frequency [MAF] SNPs); and (3) average r2hat (estimated correlation between the imputed and true genotypes, for all imputed SNPs). Across the reference panels, IMPUTE2 and MaCH had the highest concordance (91%-93%), but IMPUTE2 had the highest IQS (81%-83%) and average r2hat (0.68 using YRI+ASW+CEU, 0.62 using AFR+EUR, and 0.55 using ALL). Imputation quality for most programs was reduced by the addition of more distantly related reference populations, due entirely to the introduction of low frequency SNPs (MAF≤2%) that are monomorphic in the more closely related panels. While imputation was optimized by using IMPUTE2 with reference to the ALL panel (average r2hat = 0.86 for SNPs with MAF>2%), use of the ALL panel for African American studies requires careful interpretation of the population specificity and imputation quality of low frequency SNPs.

GATA transcription factor required for immunity to bacterial and fungal pathogens

In the past decade, Caenorhabditis elegans has been used to dissect several genetic pathways involved in immunity; however, little is known about transcription factors that regulate the expression of immune effectors. C. elegans does not appear to have a functional homolog of the key immune transcription factor NF-κB. Here we show that that the intestinal GATA transcription factor ELT-2 is required for both immunity to Salmonella enterica and expression of a C-type lectin gene, clec-67, which is expressed in the intestinal cells and is a good marker of S. enterica infection. We also found that ELT-2 is required for immunity to Pseudomonas aeruginosa, Enterococcus faecalis, and Cryptococcus neoformans. Lack of immune inhibition by DAF-2, which negatively regulates the FOXO transcription factor DAF-16, rescues the hypersusceptibility to pathogens phenotype of elt-2(RNAi) animals. Our results indicate that ELT-2 is part of a multi-pathogen defense pathway that regulates innate immunity independently of the DAF-2/DAF-16 signaling pathway.

Further investigation of simian immunodeficiency virus Vif function in human cells

Primate lentivirus Vif proteins function by suppressing the antiviral activity of the cell-encoded apolipoprotein B mRNA-editing enzyme-catalytic polypeptide-like (APOBEC) proteins APOBEC3G and APOBEC3F. It has been hypothesized that species-specific susceptibilities of APOBEC proteins to Vif proteins may help govern the transmission of primate lentiviruses to new host species. Consistent with this view and with previous results, we report that the Vif proteins of several diverse simian immunodeficiency viruses (SIVs) that are not known to infect humans are not effective inhibitors of human APOBEC3G or APOBEC3F when assessed in transient-transfection experiments. Unexpectedly, this lack of SIV Vif function did not prevent the replication of two vif-deficient SIVs (SIVtan and SIVmnd1; isolated from tantalus monkeys and mandrills, respectively) in a human T-cell line, HUT78, that expresses both APOBEC 3G and APOBEC3F, a finding which demonstrates that some SIVs are partially resistant to the antiretroviral effects of these enzymes irrespective of Vif function. Additional virus replication studies also revealed that the Vif protein of SIVtan is, in fact, active in human T cells, as it substantially enhanced the replication of its cognate virus and human immunodeficiency virus type 1. In sum, we now consider it improbable that species-specific restrictions to SIV Vif function can explain the lack of human infection with certain SIVs. Instead, our data reveal that the species-specific modulation of Vif function is more complex than previously envisioned and that additional (as-yet-unidentified) viral or host factors may be involved in regulating this dynamic interaction between host and pathogen.

The antiretroviral enzyme APOBEC3G is degraded by the proteasome in response to HIV-1 Vif

The human protein apolipoprotein B mRNA-editing enzyme-catalytic polypeptide-like-3G (APOBEC3G), also known as CEM-15, mediates a newly described form of innate resistance to retroviral infection by catalyzing the deamination of deoxycytidine to deoxyuridine in viral cDNA replication intermediates. Because DNA deamination takes place after virus entry into target cells, APOBEC3G function is dependent on its association with the viral nucleoprotein complexes that synthesize cDNA and must therefore be incorporated into virions as they assemble in infected cells. Here we show that the HIV-1 virion infectivity factor (Vif) protein protects the virus from APOBEC3G-mediated inactivation by preventing its incorporation into progeny virions, thus allowing the ensuing infection to proceed without DNA deamination. In addition to helping exclude APOBEC3G from nascent virions, Vif also removes APOBEC3G from virus-producing cells by inducing its ubiquitination and subsequent degradation by the proteasome. Our findings indicate that pharmacologic strategies aimed at stabilizing APOBEC3G in HIV-1 infected cells should be explored as potential HIV/AIDS therapeutics.

Comprehensive investigation of the molecular defect in vif-deficient human immunodeficiency virus type 1 virions

Replication of human immunodeficiency virus type 1 (HIV-1) in primary blood lymphocytes, certain T-cell lines (nonpermissive cells), and most likely in vivo is highly dependent on the virally encoded Vif protein. Evidence suggests that Vif acts late in the viral life cycle during assembly, budding, and/or maturation to counteract the antiviral activity of the CEM15 protein and possibly other antiviral factors. Because HIV-1 virions produced in the absence of Vif are severely restricted at a postentry, preintegration step of infection, it is presumed that such virions differ from wild-type virions in some way. In the present study, we established a protocol for producing large quantities of vif-deficient HIV-1 (HIV-1/Delta vif) from an acute infection of nonpermissive T cells and performed a thorough examination of the defect in these virions. Aside from the expected lack of Vif, we observed no apparent abnormalities in the packaging, modification, processing, or function of proteins in Delta vif virions. In addition, we found no consistent defect in the ability of Delta vif virions to perform intravirion reverse transcription under a variety of assay conditions, suggesting that the reverse transcription complexes in these particles can behave normally under cell-free conditions. Consistent with this finding, neither the placement of the primer tRNA3Lys nor its ability to promote reverse transcription in an in vitro assay was affected by a lack of Vif. Based on the inability of this comprehensive analysis to uncover molecular defects in Delta vif virions, we speculate that such defects are likely to be subtle and/or rare.

Isolation of a human gene that inhibits HIV-1 infection and is suppressed by the viral Vif protein

Viruses have developed diverse non-immune strategies to counteract host-mediated mechanisms that confer resistance to infection. The Vif (virion infectivity factor) proteins are encoded by primate immunodeficiency viruses, most notably human immunodeficiency virus-1 (HIV-1). These proteins are potent regulators of virus infection and replication and are consequently essential for pathogenic infections in vivo. HIV-1 Vif seems to be required during the late stages of virus production for the suppression of an innate antiviral phenotype that resides in human T lymphocytes. Thus, in the absence of Vif, expression of this phenotype renders progeny virions non-infectious. Here, we describe a unique cellular gene, CEM15, whose transient or stable expression in cells that do not normally express CEM15 recreates this phenotype, but whose antiviral action is overcome by the presence of Vif. Because the Vif:CEM15 regulatory circuit is critical for HIV-1 replication, perturbing the circuit may be a promising target for future HIV/AIDS therapies.

Unusual polymorphisms in human immunodeficiency virus type 1 associated with nonprogressive infection

Factors accounting for long-term nonprogression may include infection with an attenuated strain of human immunodeficiency virus type 1 (HIV-1), genetic polymorphisms in the host, and virus-specific immune responses. In this study, we examined eight individuals with nonprogressing or slowly progressing HIV-1 infection, none of whom were homozygous for host-specific polymorphisms (CCR5-Delta32, CCR2-64I, and SDF-1-3'A) which have been associated with slower disease progression. HIV-1 was recovered from seven of the eight, and recovered virus was used for sequencing the full-length HIV-1 genome; full-length HIV-1 genome sequences from the eighth were determined following amplification of viral sequences directly from peripheral blood mononuclear cells (PBMC). Longitudinal studies of one individual with HIV-1 that consistently exhibited a slow/low growth phenotype revealed a single amino acid deletion in a conserved region of the gp41 transmembrane protein that was not seen in any of 131 envelope sequences in the Los Alamos HIV-1 sequence database. Genetic analysis also revealed that five of the eight individuals harbored HIV-1 with unusual 1- or 2-amino-acid deletions in the Gag sequence compared to subgroup B Gag consensus sequences. These deletions in Gag have either never been observed previously or are extremely rare in the database. Three individuals had deletions in Nef, and one had a 4-amino-acid insertion in Vpu. The unusual polymorphisms in Gag, Env, and Nef described here were also found in stored PBMC samples taken 3 to 11 years prior to, or in one case 4 years subsequent to, the time of sampling for the original sequencing. In all, seven of the eight individuals exhibited one or more unusual polymorphisms; a total of 13 unusual polymorphisms were documented in these seven individuals. These polymorphisms may have been present from the time of initial infection or may have appeared in response to immune surveillance or other selective pressures. Our results indicate that unusual, difficult-to-revert polymorphisms in HIV-1 can be found associated with slow progression or nonprogression in a majority of such cases.

Evidence for a newly discovered cellular anti-HIV-1 phenotype

Animal cells have developed many ways to suppress viral replication, and viruses have evolved diverse strategies to resist these. Here we provide evidence that the virion infectivity factor protein of human immunodeficiency virus type 1 (HIV-1) functions to counteract a newly discovered activity in human cells that otherwise inhibits virus replication. This anti-viral phenotype is shown by human T cells, the principal in vivo targets for HIV-1, and, based on our present understanding of virion infectivity factor action, is presumed to act by interfering with a late step(s) in the virus life cycle. These observations indicate that the inhibition of virion infectivity factor function in vivo may prevent HIV-1 replication by 'unmasking' an innate anti-viral phenotype.