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.

Guided construction of single cell reference for human and mouse lung.

Accurate cell type identification is a key and rate-limiting step in single-cell data analysis. Single-cell references with comprehensive cell types, reproducible and functionally validated cell identities, and common nomenclatures are much needed by the research community for automated cell type annotation, data integration, and data sharing. Here, we develop a computational pipeline utilizing the LungMAP CellCards as a dictionary to consolidate single-cell transcriptomic datasets of 104 human lungs and 17 mouse lung samples to construct LungMAP single-cell reference (CellRef) for both normal human and mouse lungs. CellRefs define 48 human and 40 mouse lung cell types catalogued from diverse anatomic locations and developmental time points. We demonstrate the accuracy and stability of LungMAP CellRefs and their utility for automated cell type annotation of both normal and diseased lungs using multiple independent methods and testing data. We develop user-friendly web interfaces for easy access and maximal utilization of the LungMAP CellRefs.

A review of the genomics of neonatal abstinence syndrome.

Neonatal abstinence syndrome (NAS) is a constellation of signs of withdrawal occurring after birth following in utero exposure to licit or illicit opioids. Despite significant research and public health efforts, NAS remains challenging to diagnose, predict, and manage due to highly variable expression. Biomarker discovery in the field of NAS is crucial for stratifying risk, allocating resources, monitoring longitudinal outcomes, and identifying novel therapeutics. There is considerable interest in identifying important genetic and epigenetic markers of NAS severity and outcome that can guide medical decision making, research efforts, and public policy. A number of recent studies have suggested that genetic and epigenetic changes are associated with NAS severity, including evidence of neurodevelopmental instability. This review will provide an overview of the role of genetics and epigenetics in short and longer-term NAS outcomes. We will also describe novel research efforts using polygenic risk scores for NAS risk stratification and salivary gene expression to understand neurobehavioral modulation. Finally, emerging research focused on neuroinflammation from prenatal opioid exposure may elucidate novel mechanisms that could lead to development of future novel therapeutics.

Polygenic risk scores and the need for pharmacotherapy in neonatal abstinence syndrome.

BACKGROUND: The aim of this study was to identify genetic variants associated with NAS through a genome-wide association study (GWAS) and estimate a Polygenic Risk Score (PRS) model for NAS. METHODS: A prospective case-control study included 476 in utero opioid-exposed term neonates. A GWAS of 1000 genomes-imputed genotypes was performed to identify variants associated with need for pharmacotherapy for NAS. PRS models for estimating genetic predisposition were generated via a nested cross-validation approach using 382 neonates of European ancestry. PRS predictive ability, discrimination, and calibration were assessed. RESULTS: Cross-ancestry GWAS identified one intergenic locus on chromosome 7 downstream of SNX13 exhibiting genome-wide association with need for pharmacotherapy. PRS models derived from the GWAS for a subset of the European ancestry neonates reliably discriminated between need for pharmacotherapy using cis variant effect sizes within validation sets of European and African American ancestry neonates. PRS were less effective when applying variant effect sizes across datasets and in calibration analyses. CONCLUSIONS: GWAS has the potential to identify genetic loci associated with need for pharmacotherapy for NAS and enable development of clinically predictive PRS models. Larger GWAS with additional ancestries are needed to confirm the observed SNX13 association and the accuracy of PRS in NAS risk prediction models. IMPACT: Genetic associations appear to be important in neonatal abstinence syndrome. This is the first genome-wide association in neonates with neonatal abstinence syndrome. Polygenic risk scores can be developed examining single-nucleotide polymorphisms across the entire genome. Polygenic risk scores were higher in neonates receiving pharmacotherapy for treatment of their neonatal abstinence syndrome. Future studies with larger cohorts are needed to better delineate these genetic associations.

LungMAP Portal Ecosystem: Systems-Level Exploration of the Lung.

An improved understanding of the human lung necessitates advanced systems models informed by an ever-increasing repertoire of molecular omics, cellular, imaging, and pathological datasets. To centralize and standardize information across broad lung research efforts we expanded the LungMAP.net website into a new gateway portal. This portal connects a broad spectrum of research networks, bulk and single-cell multi-omics data and a diverse collection of image data that span mammalian lung development, and disease. The data are standardized across species and technologies using harmonized data and metadata models that leverage recent advances including those from the Human Cell Atlas, diverse ontologies, and the LungMAP CellCards initiative. To cultivate future discoveries, we have aggregated a diverse collection of single-cell atlases for multiple species (human, rhesus, mouse), to enable consistent queries across technologies, cohorts, age, disease, and drug treatment. These atlases are provided as independent and integrated queryable datasets, with an emphasis on dynamic visualization, figure generation, re-analysis, cell-type curation, and automated reference-based classification of user-provided single-cell genomics datasets (Azimuth). As this resource grows, we intend to increase the breadth of available interactive interfaces, supported data types, data portals and datasets from LungMAP and external research efforts.

Multi-trait genome-wide association study of opioid addiction: OPRM1 and beyond.

Opioid addiction (OA) is moderately heritable, yet only rs1799971, the A118G variant in OPRM1, has been identified as a genome-wide significant association with OA and independently replicated. We applied genomic structural equation modeling to conduct a GWAS of the new Genetics of Opioid Addiction Consortium (GENOA) data together with published studies (Psychiatric Genomics Consortium, Million Veteran Program, and Partners Health), comprising 23,367 cases and effective sample size of 88,114 individuals of European ancestry. Genetic correlations among the various OA phenotypes were uniformly high (rg > 0.9). We observed the strongest evidence to date for OPRM1: lead SNP rs9478500 (p = 2.56 × 10-9). Gene-based analyses identified novel genome-wide significant associations with PPP6C and FURIN. Variants within these loci appear to be pleiotropic for addiction and related traits.

GAWMerge expands GWAS sample size and diversity by combining array-based genotyping and whole-genome sequencing.

Genome-wide association studies (GWAS) have made impactful discoveries for complex diseases, often by amassing very large sample sizes. Yet, GWAS of many diseases remain underpowered, especially for non-European ancestries. One cost-effective approach to increase sample size is to combine existing cohorts, which may have limited sample size or be case-only, with public controls, but this approach is limited by the need for a large overlap in variants across genotyping arrays and the scarcity of non-European controls. We developed and validated a protocol, Genotyping Array-WGS Merge (GAWMerge), for combining genotypes from arrays and whole-genome sequencing, ensuring complete variant overlap, and allowing for diverse samples like Trans-Omics for Precision Medicine to be used. Our protocol involves phasing, imputation, and filtering. We illustrated its ability to control technology driven artifacts and type-I error, as well as recover known disease-associated signals across technologies, independent datasets, and ancestries in smoking-related cohorts. GAWMerge enables genetic studies to leverage existing cohorts to validly increase sample size and enhance discovery for understudied traits and ancestries.

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.

Genetically predicted serum vitamin D and COVID-19: a Mendelian randomisation study.

OBJECTIVES: To investigate causality of the association of serum vitamin D with the risk and severity of COVID-19 infection. DESIGN: Two-sample Mendelian randomisation study. SETTING: Summary data from genome-wide analyses in the population-based UK Biobank and SUNLIGHT Consortium, applied to meta-analysed results of genome-wide analyses in the COVID-19 Host Genetics Initiative. PARTICIPANTS: 17 965 COVID-19 cases including 11 085 laboratory or physician-confirmed cases, 7885 hospitalised cases and 4336 severe respiratory cases, and 1 370 547 controls, primarily of European ancestry. EXPOSURES: Genetically predicted variation in serum vitamin D status, instrumented by genome-wide significant single nucleotide polymorphisms (SNPs) associated with serum vitamin D or risk of vitamin D deficiency/insufficiency. MAIN OUTCOME MEASURES: Susceptibility to and severity of COVID-19 infection, including severe respiratory infection and hospitalisation. RESULTS: Mendelian randomisation analysis, sufficiently powered to detect effects comparable to those seen in observational studies, provided little to no evidence for an effect of genetically predicted serum vitamin D on susceptibility to or severity of COVID-19 infection. Using SNPs in loci related to vitamin D metabolism as genetic instruments for serum vitamin D concentrations, the OR per SD higher serum vitamin D was 1.04 (95% CI 0.92 to 1.18) for any COVID-19 infection versus population controls, 1.05 (0.84 to 1.31) for hospitalised COVID-19 versus population controls, 0.96 (0.64 to 1.43) for severe respiratory COVID-19 versus population controls, 1.15 (0.99 to 1.35) for COVID-19 positive versus COVID-19 negative and 1.44 (0.75 to 2.78) for hospitalised COVID-19 versus non-hospitalised COVID-19. Results were similar in analyses using SNPs with genome-wide significant associations with serum vitamin D (ie, including SNPs in loci with no known relationship to vitamin D metabolism) and in analyses using SNPs with genome-wide significant associations with risk of vitamin D deficiency or insufficiency. CONCLUSIONS: These findings suggest that genetically predicted differences in long-term vitamin D nutritional status do not causally affect susceptibility to and severity of COVID-19 infection, and that associations observed in previous studies may have been driven by confounding. These results do not exclude the possibility of low-magnitude causal effects or causal effects of acute responses to therapeutic doses of vitamin D.

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.

Ontology-guided segmentation and object identification for developmental mouse lung immunofluorescent images.

BACKGROUND: Immunofluorescent confocal microscopy uses labeled antibodies as probes against specific macromolecules to discriminate between multiple cell types. For images of the developmental mouse lung, these cells are themselves organized into densely packed higher-level anatomical structures. These types of images can be challenging to segment automatically for several reasons, including the relevance of biomedical context, dependence on the specific set of probes used, prohibitive cost of generating labeled training data, as well as the complexity and dense packing of anatomical structures in the image. The use of an application ontology helps surmount these challenges by combining image data with its metadata to provide a meaningful biological context, modeled after how a human expert would make use of contextual information to identify histological structures, that constrains and simplifies the process of segmentation and object identification. RESULTS: We propose an innovative approach for the semi-supervised analysis of complex and densely packed anatomical structures from immunofluorescent images that utilizes an application ontology to provide a simplified context for image segmentation and object identification. We describe how the logical organization of biological facts in the form of an ontology can provide useful constraints that facilitate automatic processing of complex images. We demonstrate the results of ontology-guided segmentation and object identification in mouse developmental lung images from the Bioinformatics REsource ATlas for the Healthy lung database of the Molecular Atlas of Lung Development (LungMAP1) program CONCLUSION: We describe a novel ontology-guided approach to segmentation and classification of complex immunofluorescence images of the developing mouse lung. The ontology is used to automatically generate constraints for each image based on its biomedical context, which facilitates image segmentation and classification.

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.

Spatial distribution of marker gene activity in the mouse lung during alveolarization.

This data is a curated collection of visual images of gene expression patterns from the pre- and post-natal mouse lung, accompanied by associated mRNA probe sequences and RNA-Seq expression profiles. Mammalian lungs undergo significant growth and cellular differentiation before and after the transition to breathing air. Documenting normal lung development is an important step in understanding abnormal lung development, as well as the challenges faced during a preterm birth. Images in this dataset indicate the spatial distribution of mRNA transcripts for over 500 different genes that are active during lung development, as initially determined via RNA-Seq. Images were systematically acquired using high-throughput in situ hybridization with non-radioactive digoxigenin-labeled mRNA probes across mouse lungs from developmental time points E16.5, E18.5, P7, and P28. The dataset was produced as part of The Molecular Atlas of Lung Development Program (LungMAP) and is hosted at https://lungmap.net. This manuscript describes the nature of the data and the protocols for generating the dataset.

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.

Genomewide association study of HLA alloimmunization in previously pregnant blood donors.

BACKGROUND: Alloimmunization through blood transfusion, transplantation, or circulating fetal cells during pregnancy is a significant concern. Some exposed individuals make alloantibodies while others do not, implying variation in genetic risk factors. STUDY DESIGN AND METHODS: We conducted a genomewide association study (GWAS) of 9,427,497 single-nucleotide polymorphisms (SNPs) to identify genetic variants for HLA alloimmunization in previously pregnant blood donors with (n = 752) and without (n = 753) HLA Class I or II alloantibodies. RESULTS: A SNP in the neurexophilin 2 (NXPH2) gene surpassed genome-wide significance (p = 2.06 × 10-8 ), with multiple adjacent markers p < 10-6 , for women with anti-Class I alloantibodies only. Little is currently known about the function of NXPH2, although gene family members have been shown to impact immunity. SNPs in the E2F7 gene, a transcription factor related to cell cycle control and cellular proliferation, also approached genomewide significance (p = 2.5 × 10-7 ). CONCLUSION: Further work to extend the GWAS approach and to characterize variants in NXPH2 and E2F7 in the context of alloantibody formation is warranted.

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.

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.

Associations of common variants in the BST2 region with HIV-1 acquisition in African American and European American people who inject drugs.

OBJECTIVE: The bone marrow stromal cell antigen 2 (BST2) gene encodes a host restriction factor that acts as an innate immune sensor of HIV-1 exposure and suppresses release of HIV-1 particles. We aimed to identify associations of variants in the BST2 gene region with HIV-1 acquisition and disease progression. DESIGN/METHODS: Using HIV+ cases and HIV- controls from the Urban Health Study (n=3136 African Americans and European Americans who inject drugs), we tested 470 variants in BST2 and its flanking regions for association with HIV-1 acquisition and log-transformed viral load. RESULTS: We found that the single nucleotide polymorphism (SNP) rs113189798 surpassed the P value threshold corrected for multiple testing. The rs113189798-G allele (frequency=16% in African Americans, 4% in European Americans) was associated with increased HIV-1 acquisition risk (meta-analysis P=1.43 × 10): odds ratio (95% confidence interval) of 1.22 (1.01-1.49) in African Americans and 2.17 (1.43-3.33) in European Americans. We also found that the previously reported rs12609479-A allele (frequency=35% in African Americans, 81% in European Americans) was nominally associated with decreased risk of acquiring HIV-1 in our study (meta-analysis P=0.036). Rs12609479-A is predicted to increase BST2 expression and thereby decrease risk of acquiring HIV-1. Rs113189798 and rs12609479 were only weakly correlated [square of the correlation coefficient (r)=0.2-0.4] and represented distinct association signals. None of our tested variants were significantly associated with log-transformed viral load among the HIV-infected cases. CONCLUSION: Our findings support BST2 as a genetic susceptibility factor for HIV-1 acquisition: identifying a novel SNP association for rs13189798 and linking the previously reported regulatory SNP rs12609479 to HIV-1 acquisition.