Acceleration of Age-Associated Methylation Patterns in HIV-1-Infected Adults

Patients with treated HIV-1-infection experience earlier occurrence of aging-associated diseases, raising speculation that HIV-1-infection, or antiretroviral treatment, may accelerate aging. We recently described an age-related co-methylation module comprised of hundreds of CpGs; however, it is unknown whether aging and HIV-1-infection exert negative health effects through similar, or disparate, mechanisms. We investigated whether HIV-1-infection would induce age-associated methylation changes. We evaluated DNA methylation levels at >450,000 CpG sites in peripheral blood mononuclear cells (PBMC) of young (20-35) and older (36-56) adults in two separate groups of participants. Each age group for each data set consisted of 12 HIV-1-infected and 12 age-matched HIV-1-uninfected samples for a total of 96 samples. The effects of age and HIV-1 infection on methylation at each CpG revealed a strong correlation of 0.49, p<1 x10^-200 and 0.47, p<1x10^-200. Weighted gene correlation network analysis (WGCNA) identified 17 co-methylation modules; module 3 (ME3) was significantly correlated with age (cor=0.70) and HIV-1 status (cor=0.31). Older HIV-1⁺ individuals had a greater number of hypermethylated CpGs across ME3 (p=0.015). In a multivariate model, ME3 was significantly associated with age and HIV status (Data set 1: β_age= 0.007088, p=2.08 x 10^-9; β_HIV= 0.099574, p=0.0011; Data set 2: β_age= 0.008762, p=1.27x 10^-5; β_HIV= 0.128649, p= 0.0001). Using this model, we estimate that HIV-1 infection accelerates age-related methylation by approximately 13.7 years in data set 1 and 14.7 years in data set 2. The genes related to CpGs in ME3 are enriched for polycomb group target genes known to be involved in cell renewal and aging. The overlap between ME3 and an aging methylation module found in solid tissues is also highly significant (Fisher-exact p=5.6 x 10^-6, odds ratio=1.91). These data demonstrate that HIV-1 infection is associated with methylation patterns that are similar to age-associated patterns and suggest that general aging and HIV-1 related aging work through some common cellular and molecular mechanisms. These results are an important first step for finding potential therapeutic targets and novel clinical approaches to mitigate the detrimental effects of both HIV-1-infection and aging.     

Six Novel Loci Associated with Circulating VEGF Levels Identified by a Meta-analysis of Genome-Wide Association Studies

Vascular endothelial growth factor (VEGF) is an angiogenic and neurotrophic factor, secreted by endothelial cells, known to impact various physiological and disease processes from cancer to cardiovascular disease and to be pharmacologically modifiable. We sought to identify novel loci associated with circulating VEGF levels through a genome-wide association meta-analysis combining data from European-ancestry individuals and using a dense variant map from 1000 genomes imputation panel. Six discovery cohorts including 13,312 samples were analyzed, followed by in-silico and de-novo replication studies including an additional 2,800 individuals. A total of 10 genome-wide significant variants were identified at 7 loci. Four were novel loci (5q14.3, 10q21.3, 16q24.2 and 18q22.3) and the leading variants at these loci were rs114694170 (MEF2C, P = 6.79x10^-13), rs74506613 (JMJD1C, P = 1.17x10^-19), rs4782371 (ZFPM1, P = 1.59x10^-9) and rs2639990 (ZADH2, P = 1.72x10^-8), respectively. We also identified two new independent variants (rs34528081, VEGFA, P = 1.52x10^-18; rs7043199, VLDLR-AS1, P = 5.12x10^-14) at the 3 previously identified loci and strengthened the evidence for the four previously identified SNPs (rs6921438, LOC100132354, P = 7.39x10^-1467; rs1740073, C6orf223, P = 2.34x10^-17; rs6993770, ZFPM2, P = 2.44x10^-60; rs2375981, KCNV2, P = 1.48x10^-100). These variants collectively explained up to 52% of the VEGF phenotypic variance. We explored biological links between genes in the associated loci using Ingenuity Pathway Analysis that emphasized their roles in embryonic development and function. Gene set enrichment analysis identified the ERK5 pathway as enriched in genes containing VEGF associated variants. eQTL analysis showed, in three of the identified regions, variants acting as both cis and trans eQTLs for multiple genes. Most of these genes, as well as some of those in the associated loci, were involved in platelet biogenesis and functionality, suggesting the importance of this process in regulation of VEGF levels. This work also provided new insights into the involvement of genes implicated in various angiogenesis related pathologies in determining circulating VEGF levels. The understanding of the molecular mechanisms by which the identified genes affect circulating VEGF levels could be important in the development of novel VEGF-related therapies for such diseases.     

Genetic Determinants of Pelvic Organ Prolapse among African American and Hispanic Women in the Women’s Health Initiative

Current evidence suggests a multifactorial etiology to pelvic organ prolapse (POP), including genetic predisposition. We conducted a genome-wide association study of POP in African American (AA) and Hispanic (HP) women from the Women’s Health Initiative Hormone Therapy study. Cases were defined as any POP (grades 1–3) or moderate/severe POP (grades 2–3), while controls had grade 0 POP. We performed race-specific multiple logistic regression analyses between SNPs imputed to 1000 genomes in relation to POP (grade 0 vs 1–3; grade 0 vs 2–3) adjusting for age at diagnosis, body mass index, parity, and genetic ancestry. There were 1274 controls and 1427 cases of any POP and 317 cases of moderate/severe POP. Although none of the analyses reached genome-wide significance (p<5x10^-8), we noted variants in several loci that met p<10⁻⁶. In race-specific analysis of grade 0 vs 2–3, intronic SNPs in the CPE gene (rs28573326, OR:2.14; 95% CI 1.62–2.83; p = 1.0x10^-7) were associated with POP in AAs, and SNPs in the gene AL132709.5 (rs1950626, OR:2.96; 95% CI 1.96–4.48, p = 2.6x10^-7) were associated with POP in HPs. Inverse variance fixed-effect meta-analysis of the race-specific results showed suggestive signals for SNPs in the DPP6 gene (rs11243354, OR:1.36; p = 4.2x10^-7) in the grade 0 vs 1–3 analyses and for SNPs around PGBD5 (rs740494, OR:2.17; p = 8.6x10^-7) and SHC3 (rs2209875, OR:0.60; p = 9.3x10^-7) in the grade 0 vs 2–3 analyses. While we did not identify genome-wide significant findings, we document several SNPs reaching suggestive statistical significance. Further interrogation of POP in larger minority samples is warranted.     

Replication of EPHA1 and CD33 associations with late-onset Alzheimer's disease: a multi-centre case-control study

Background

A recently published genome-wide association study (GWAS) of late-onset Alzheimer's disease (LOAD) revealed genome-wide significant association of variants in or near MS4A4A, CD2AP, EPHA1 and CD33. Meta-analyses of this and a previously published GWAS revealed significant association at ABCA7 and MS4A, independent evidence for association of CD2AP, CD33 and EPHA1 and an opposing yet significant association of a variant near ARID5B. In this study, we genotyped five variants (in or near CD2AP, EPHA1, ARID5B, and CD33) in a large (2,634 LOAD, 4,201 controls), independent dataset comprising six case-control series from the USA and Europe. We performed meta-analyses of the association of these variants with LOAD and tested for association using logistic regression adjusted by age-at-diagnosis, gender, and APOE ε4 dosage.

Results

We found no significant evidence of series heterogeneity. Associations with LOAD were successfully replicated for EPHA1 (rs11767557; OR = 0.87, p = 5 × 10^-4) and CD33 (rs3865444; OR = 0.92, p = 0.049), with odds ratios comparable to those previously reported. Although the two ARID5B variants (rs2588969 and rs494288) showed significant association with LOAD in meta-analysis of our dataset (p = 0.046 and 0.008, respectively), the associations did not survive adjustment for covariates (p = 0.30 and 0.11, respectively). We had insufficient evidence in our data to support the association of the CD2AP variant (rs9349407, p = 0.56).

Conclusions

Our data overwhelmingly support the association of EPHA1 and CD33 variants with LOAD risk: addition of our data to the results previously reported (total n > 42,000) increased the strength of evidence for these variants, providing impressive p-values of 2.1 × 10^-15 (EPHA1) and 1.8 × 10^-13 (CD33).     

Parent-of-Origin Effects of the APOB Gene on Adiposity in Young Adults

Loci identified in genome-wide association studies (GWAS) of cardio-metabolic traits account for a small proportion of the traits'' heritability. To date, most association studies have not considered parent-of-origin effects (POEs). Here we report investigation of POEs on adiposity and glycemic traits in young adults. The Jerusalem Perinatal Family Follow-Up Study (JPS), comprising 1250 young adults and their mothers was used for discovery. Focusing on 18 genes identified by previous GWAS as associated with cardio-metabolic traits, we used linear regression to examine the associations of maternally- and paternally-derived offspring minor alleles with body mass index (BMI), waist circumference (WC), fasting glucose and insulin. We replicated and meta-analyzed JPS findings in individuals of European ancestry aged ≤50 belonging to pedigrees from the Framingham Heart Study, Family Heart Study and Erasmus Rucphen Family study (total N≅4800). We considered p<2.7x10^-4 statistically significant to account for multiple testing. We identified a common coding variant in the 4^th exon of APOB (rs1367117) with a significant maternally-derived effect on BMI (β = 0.8; 95%CI:0.4,1.1; p = 3.1x10^-5) and WC (β = 2.7; 95%CI:1.7,3.7; p = 2.1x10^-7). The corresponding paternally-derived effects were non-significant (p>0.6). Suggestive maternally-derived associations of rs1367117 were observed with fasting glucose (β = 0.9; 95%CI:0.3,1.5; p = 4.0x10^-3) and insulin (ln-transformed, β = 0.06; 95%CI:0.03,0.1; p = 7.4x10^-4). Bioinformatic annotation for rs1367117 revealed a variety of regulatory functions in this region in liver and adipose tissues and a 50% methylation pattern in liver only, consistent with allelic-specific methylation, which may indicate tissue-specific POE. Our findings demonstrate a maternal-specific association between a common APOB variant and adiposity, an association that was not previously detected in GWAS. These results provide evidence for the role of regulatory mechanisms, POEs specifically, in adiposity. In addition this study highlights the benefit of utilizing family studies for deciphering the genetic architecture of complex traits.     

A Genome Wide Association Study Links Glutamate Receptor Pathway to Sporadic Creutzfeldt-Jakob Disease Risk

We performed a genome-wide association (GWA) study in 434 sporadic Creutzfeldt-Jakob disease (sCJD) patients and 1939 controls from the United Kingdom, Germany and The Netherlands. The findings were replicated in an independent sample of 1109 sCJD and 2264 controls provided by a multinational consortium. From the initial GWA analysis we selected 23 SNPs for further genotyping in 1109 sCJD cases from seven different countries. Five SNPs were significantly associated with sCJD after correction for multiple testing. Subsequently these five SNPs were genotyped in 2264 controls. The pooled analysis, including 1543 sCJD cases and 4203 controls, yielded two genome wide significant results: rs6107516 (p-value=7.62x10^-9) a variant tagging the prion protein gene (PRNP); and rs6951643 (p-value=1.66x10^-8) tagging the Glutamate Receptor Metabotropic 8 gene (GRM8). Next we analysed the data stratifying by country of origin combining samples from the pooled analysis with genotypes from the 1000 Genomes Project and imputed genotypes from the Rotterdam Study (Total n=12967). The meta-analysis of the results showed that rs6107516 (p-value=3.00x10^-8) and rs6951643 (p-value=3.91x10^-5) remained as the two most significantly associated SNPs. Rs6951643 is located in an intronic region of GRM8, a gene that was additionally tagged by a cluster of 12 SNPs within our top100 ranked results. GRM8 encodes for mGluR8, a protein which belongs to the metabotropic glutamate receptor family, recently shown to be involved in the transduction of cellular signals triggered by the prion protein. Pathway enrichment analyses performed with both Ingenuity Pathway Analysis and ALIGATOR postulates glutamate receptor signalling as one of the main pathways associated with sCJD. In summary, we have detected GRM8 as a novel, non-PRNP, genome-wide significant marker associated with heightened disease risk, providing additional evidence supporting a role of glutamate receptors in sCJD pathogenesis.     

Callus induction and plantlet regeneration in damiana (Turnera diffusa,Willd.)

Experiments were performed to determine optimal conditions for callus induction and plantlet regeneration from damiana (Turnera diffusa, Willd.) leaf explants. MS media [7] and B5 media [4] were used; both of them were supplemented with different combinations of 2,4-dichlorophenoxyacetic, acid (2,4-D) and 6-benzyl-adenine (BA). In all cases the samples were incubated under continous fluorescent light, 100 mol.m^-2.sec^-1, at 25±2°C for 60 days. Highest callus percentage was obtained in MS media supplemented with 2,4-D 10^–5 M plus BA 10^-8M. Maximum percentage and multiple shoots and leaves were produced in MS media with BA 10^-5M and BA 10^-5M plus 2,4-D 10^-7M. These explants were transferred for rooting to MS medium without plant growth regulators. Later, rooted plantlets were successfully established in soil.     

A Phenome-Wide Association Study of genes associated with COVID-19 severity reveals shared genetics with complex diseases in the Million Veteran Program

The study aims to determine the shared genetic architecture between COVID-19 severity with existing medical conditions using electronic health record (EHR) data. We conducted a Phenome-Wide Association Study (PheWAS) of genetic variants associated with critical illness (n = 35) or hospitalization (n = 42) due to severe COVID-19 using genome-wide association summary data from the Host Genetics Initiative. PheWAS analysis was performed using genotype-phenotype data from the Veterans Affairs Million Veteran Program (MVP). Phenotypes were defined by International Classification of Diseases (ICD) codes mapped to clinically relevant groups using published PheWAS methods. Among 658,582 Veterans, variants associated with severe COVID-19 were tested for association across 1,559 phenotypes. Variants at the ABO locus (rs495828, rs505922) associated with the largest number of phenotypes (n_rs495828 = 53 and n_rs505922 = 59); strongest association with venous embolism, odds ratio (OR_rs495828 1.33 (p = 1.32 x 10⁻¹⁹⁹), and thrombosis OR_rs505922 1.33, p = 2.2 x10^-265. Among 67 respiratory conditions tested, 11 had significant associations including MUC5B locus (rs35705950) with increased risk of idiopathic fibrosing alveolitis OR 2.83, p = 4.12 × 10⁻¹⁹¹; CRHR1 (rs61667602) associated with reduced risk of pulmonary fibrosis, OR 0.84, p = 2.26× 10⁻¹². The TYK2 locus (rs11085727) associated with reduced risk for autoimmune conditions, e.g., psoriasis OR 0.88, p = 6.48 x10^-23, lupus OR 0.84, p = 3.97 x 10⁻⁰⁶. PheWAS stratified by ancestry demonstrated differences in genotype-phenotype associations. LMNA (rs581342) associated with neutropenia OR 1.29 p = 4.1 x 10⁻¹³ among Veterans of African and Hispanic ancestry but not European. Overall, we observed a shared genetic architecture between COVID-19 severity and conditions related to underlying risk factors for severe and poor COVID-19 outcomes. Differing associations between genotype-phenotype across ancestries may inform heterogenous outcomes observed with COVID-19. Divergent associations between risk for severe COVID-19 with autoimmune inflammatory conditions both respiratory and non-respiratory highlights the shared pathways and fine balance of immune host response and autoimmunity and caution required when considering treatment targets.     

Genetic Variants in the Bone Morphogenic Protein Gene Family Modify the Association between Residential Exposure to Traffic and Peripheral Arterial Disease

There is a growing literature indicating that genetic variants modify many of the associations between environmental exposures and clinical outcomes, potentially by increasing susceptibility to these exposures. However, genome-scale investigations of these interactions have been rarely performed particularly in the case of air pollution exposures. We performed race-stratified genome-wide gene-environment interaction association studies on European-American (EA, N = 1623) and African-American (AA, N = 554) cohorts to investigate the joint influence of common single nucleotide polymorphisms (SNPs) and residential exposure to traffic (“traffic exposure”)—a recognized vascular disease risk factor—on peripheral arterial disease (PAD). Traffic exposure was estimated via the distance from the primary residence to the nearest major roadway, defined as the nearest limited access highways or major arterial. The rs755249-traffic exposure interaction was associated with PAD at a genome-wide significant level (P = 2.29x10^-8) in European-Americans. Rs755249 is located in the 3’ untranslated region of BMP8A, a member of the bone morphogenic protein (BMP) gene family. Further investigation revealed several variants in BMP genes associated with PAD via an interaction with traffic exposure in both the EA and AA cohorts; this included interactions with non-synonymous variants in BMP2, which is regulated by air pollution exposure. The BMP family of genes is linked to vascular growth and calcification and is a novel gene family for the study of PAD pathophysiology. Further investigation of BMP8A using the Genotype Tissue Expression Database revealed multiple variants with nominally significant (P < 0.05) interaction P-values in our EA cohort were significant BMP8A eQTLs in tissue types highlight relevant for PAD such as rs755249 (tibial nerve, eQTL P = 3.6x10^-6) and rs1180341 (tibial artery, eQTL P = 5.3x10^-6). Together these results reveal a novel gene, and possibly gene family, associated with PAD via an interaction with traffic air pollution exposure. These results also highlight the potential for interactions studies, particularly at the genome scale, to reveal novel biology linking environmental exposures to clinical outcomes.     

APLP2 Regulates Refractive Error and Myopia Development in Mice and Humans

Myopia is the most common vision disorder and the leading cause of visual impairment worldwide. However, gene variants identified to date explain less than 10% of the variance in refractive error, leaving the majority of heritability unexplained (“missing heritability”). Previously, we reported that expression of APLP2 was strongly associated with myopia in a primate model. Here, we found that low-frequency variants near the 5’-end of APLP2 were associated with refractive error in a prospective UK birth cohort (n = 3,819 children; top SNP rs188663068, p = 5.0 × 10⁻⁴) and a CREAM consortium panel (n = 45,756 adults; top SNP rs7127037, p = 6.6 × 10⁻³). These variants showed evidence of differential effect on childhood longitudinal refractive error trajectories depending on time spent reading (gene x time spent reading x age interaction, p = 4.0 × 10⁻³). Furthermore, Aplp2 knockout mice developed high degrees of hyperopia (+11.5 ± 2.2 D, p < 1.0 × 10⁻⁴) compared to both heterozygous (-0.8 ± 2.0 D, p < 1.0 × 10⁻⁴) and wild-type (+0.3 ± 2.2 D, p < 1.0 × 10⁻⁴) littermates and exhibited a dose-dependent reduction in susceptibility to environmentally induced myopia (F(2, 33) = 191.0, p < 1.0 × 10⁻⁴). This phenotype was associated with reduced contrast sensitivity (F(12, 120) = 3.6, p = 1.5 × 10⁻⁴) and changes in the electrophysiological properties of retinal amacrine cells, which expressed Aplp2. This work identifies APLP2 as one of the “missing” myopia genes, demonstrating the importance of a low-frequency gene variant in the development of human myopia. It also demonstrates an important role for APLP2 in refractive development in mice and humans, suggesting a high level of evolutionary conservation of the signaling pathways underlying refractive eye development.     

The Virological and Immunological Characteristics of the HIV-1-Infected Population in Brazil: From Initial Diagnosis to Impact of Antiretroviral Use

BackgroundImmunological and virological status of HIV-infected individuals entering the Brazilian public system over time was analyzed. We evaluated the impact of ART on virological, immunological and antiretroviral resistance over time.MethodsCD4+ T cell counts, viral loads and genotypes from patients over 13 years old from 2001–2011 were analyzed according to demographic data. We compared groups using parametric t-tests and linear regression analysis in the R statistical software language.ResultsMean baseline CD4+ T cell counts varied from 348 (2003) to 389 (2009) and was higher among women (p = 1.1 x 10⁻⁸), lower in older patients (p< 1 x 10⁻⁸) and lower in less developed regions (p = 1.864 x 10⁻⁵). Percentage of treated patients with undetectable viral loads increased linearly from 46% (2001) to 77% (2011), was lower among women (p = 2.851 x 10⁻⁶), younger ages (p = 1 x 10⁻³), and in less developed regions (p = 1.782 x 10⁻⁴). NRTI acquired resistance was 86% in 2001–3 and decreased over time. NNRTI resistance increased from 2001-3(50%) to 2006–9 (60%), PI resistance decreased from 2001–3 (60%) to 2009 (40%), and 3-class resistance was stable over time around 25%. Subtype prevalence comprised B (75.3%), B/F recombinants (12.2%), C (5.7%), F (5.3%) and B/C recombinants (1.5%), with regional variations. Three-class resistance was 26.5% among Bs, 22.4% among Fs and 17.2% among Cs.ConclusionsHIV diagnosis occurs late, especially among elderly Brazilians. Younger individuals need special attention due to poor virological response to treatment. Antiretroviral Resistance profile is subtype related.     

A genome-wide association study identifies 5 loci associated with frozen shoulder and implicates diabetes as a causal risk factor

Frozen shoulder is a painful condition that often requires surgery and affects up to 5% of individuals aged 40–60 years. Little is known about the causes of the condition, but diabetes is a strong risk factor. To begin to understand the biological mechanisms involved, we aimed to identify genetic variants associated with frozen shoulder and to use Mendelian randomization to test the causal role of diabetes. We performed a genome-wide association study (GWAS) of frozen shoulder in the UK Biobank using data from 10,104 cases identified from inpatient, surgical and primary care codes. We used data from FinnGen for replication and meta-analysis. We used one-sample and two-sample Mendelian randomization approaches to test for a causal association of diabetes with frozen shoulder. We identified five genome-wide significant loci. The most significant locus (lead SNP rs28971325; OR = 1.20, [95% CI: 1.16–1.24], p = 5x10^-29) contained WNT7B. This variant was also associated with Dupuytren’s disease (OR = 2.31 [2.24, 2.39], p<1x10^-300) as were a further two of the frozen shoulder associated variants. The Mendelian randomization results provided evidence that type 1 diabetes is a causal risk factor for frozen shoulder (OR = 1.03 [1.02–1.05], p = 3x10^-6). There was no evidence that obesity was causally associated with frozen shoulder, suggesting that diabetes influences risk of the condition through glycemic rather than mechanical effects. We have identified genetic loci associated with frozen shoulder. There is a large overlap with Dupuytren’s disease associated loci. Diabetes is a likely causal risk factor. Our results provide evidence of biological mechanisms involved in this common painful condition.     

Molecular Marker Based Coefficient of Parentage Analysis for Establishing Distinctness in Indian Rice (Oryza sativa L) Varieties

Molecular characterization of 32 Indian rice varieties of different agro-climatic zones resulted in mean heterozygosity value of 0.622, 0.819 and 0.890 over polymorphic loci and marker index value of 1.00, 6.75 and 4.16, respectively for RAPD, ISSR and STMS primers. The three marker systems resulted in 201 polymorphic bands (94.36%) out of a total of 213 bands. The probability of a chance identical match between two varieties was very low (2.08x10¹⁰) in combined molecular marker analysis as compared to individual marker system (RAPD, 7.5 x 10^-4; ISSR, 1.5 x 10^-3 and STMS, 3.9 x 10^-6). The combined average genetic similarity for molecular markers and coefficient of parentage (COP) across all 496 pairwise combinations revealed a non-significant relationship (r = 0.215).     

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.     

A GWAS Study on Liver Function Test Using eMERGE Network Participants

Introduction

Liver enzyme levels and total serum bilirubin are under genetic control and in recent years genome-wide population-based association studies have identified different susceptibility loci for these traits. We conducted a genome-wide association study in European ancestry participants from the Electronic Medical Records and Genomics (eMERGE) Network dataset of patient medical records with available genotyping data in order to identify genetic contributors to variability in serum bilirubin levels and other liver function tests and to compare the effects between adult and pediatric populations.

Methods

The process of whole genome imputation of eMERGE samples with standard quality control measures have been described previously. After removing missing data and outliers based on principal components (PC) analyses, 3294 samples from European ancestry were used for the GWAS study. The association between each single nucleotide polymorphism (SNP) and total serum bilirubin and other liver function tests was tested using linear regression, adjusting for age, gender, site, platform and ancestry principal components (PC).

Results

Consistent with previous results, a strong association signal has been detected for UGT1A gene cluster (best SNP rs887829, beta = 0.15, p = 1.30x10^-118) for total serum bilirubin level. Indeed, in this region more than 176 SNPs (or indels) had p<10⁻⁸ spanning 150Kb on the long arm of chromosome 2q37.1. In addition, we found a similar level of magnitude in a pediatric group (p = 8.26x10^-47, beta = 0.17). Further imputation using sequencing data as a reference panel revealed association of other markers including known TA7 repeat indels (rs8175347) (p = 9.78x10^-117) and rs111741722 (p = 5.41x10^-119) which were in proxy (r2 = 0.99) with rs887829. Among rare variants, two Asian subjects homozygous for coding SNP rs4148323 (G71R) were identified. Additional known effects for total serum bilirubin were also confirmed including organic anion transporters SLCO1B1-SLCO1B3, TDRP and ZMYND8 at FDR<0.05 with no gene-gene interaction effects. Phenome-wide association studies (PheWAS) suggest a protective effect of TA7 repeat against cerebrovascular disease in an adult cohort (OR = 0.75, p = 0.0008). Among other liver function tests, we also confirmed the previous effect of the ABO blood group locus for variation in serum alkaline phosphatase (rs579459, p = 9.44x10^-15).

Conclusions

Taken together, our data present interesting findings with strong confirmation of previous effects by simply using the eMERGE electronic health record phenotyping. In addition, our findings indicate that similar to the adult population, the UGT1A1 is the main locus responsible for normal variation of serum bilirubin in pediatric populations.