Pharmacogenomics of HIV Therapy: Summary of a Workshop Sponsored by the National Institute of Allergy and Infectious Diseases

Approximately 1 million people in the United States and over 30 million worldwide are living with human immunodeficiency virus type 1 (HIV-1). While mortality from untreated infection approaches 100%, survival improves markedly with use of contemporary antiretroviral therapies (ART). In the United States, 25 drugs are approved for treating HIV-1, and increasing numbers are available in resource-limited countries. Safe and effective ART is a cornerstone in the global struggle against the acquired immunodeficiency syndrome. Variable responses to ART are due at least in part to human genetic variants that affect drug metabolism, drug disposition, and off-site drug targets. Defining effects of human genetic variants on HIV treatment toxicity, efficacy, and pharmacokinetics has far-reaching implications. In 2010, the National Institute of Allergy and Infectious Diseases sponsored a workshop entitled, Pharmacogenomics A Path Towards Personalized HIV Care. This article summarizes workshop objectives, presentations, discussions, and recommendations derived from this meeting.

Tandem repeat sequence variation as causative cis-eQTLs for protein-coding gene expression variation: the case of CSTB.

Association studies have revealed expression quantitative trait loci (eQTLs) for a large number of genes. However, the causative variants that regulate gene expression levels are generally unknown. We hypothesized that copy-number variation of sequence repeats contribute to the expression variation of some genes. Our laboratory has previously identified that the rare expansion of a repeat c.-174CGGGGCGGGGCG in the promoter region of the CSTB gene causes a silencing of the gene, resulting in progressive myoclonus epilepsy. Here, we genotyped the repeat length and quantified CSTB expression by quantitative real-time polymerase chain reaction in 173 lymphoblastoid cell lines (LCLs) and fibroblast samples from the GenCord collection. The majority of alleles contain either two or three copies of this repeat. Independent analysis revealed that the c.-174CGGGGCGGGGCG repeat length is strongly associated with CSTB expression (P = 3.14 × 10(-11)) in LCLs only. Examination of both genotyped and imputed single-nucleotide polymorphisms (SNPs) within 2 Mb of CSTB revealed that the dodecamer repeat represents the strongest cis-eQTL for CSTB in LCLs. We conclude that the common two or three copy variation is likely the causative cis-eQTL for CSTB expression variation. More broadly, we propose that polymorphic tandem repeats may represent the causative variation of a fraction of cis-eQTLs in the genome.

Molecular Proxies for Climate Maladaptation in a Long-Lived Tree (Pinus pinaster Aiton, Pinaceae).

Understanding adaptive genetic responses to climate change is a main challenge for preserving biological diversity. Successful predictive models for climate-driven range shifts of species depend on the integration of information on adaptation, including that derived from genomic studies. Long-lived forest trees can experience substantial environmental change across generations, which results in a much more prominent adaptation lag than in annual species. Here, we show that candidate-gene SNPs (single nucleotide polymorphisms) can be used as predictors of maladaptation to climate in maritime pine (Pinus pinaster Aiton), an outcrossing long-lived keystone tree. A set of 18 SNPs potentially associated with climate, 5 of them involving amino acid-changing variants, were retained after performing logistic regression, latent factor mixed models, and Bayesian analyses of SNP-climate correlations. These relationships identified temperature as an important adaptive driver in maritime pine and highlighted that selective forces are operating differentially in geographically discrete gene pools. The frequency of the locally advantageous alleles at these selected loci was strongly correlated with survival in a common garden under extreme (hot and dry) climate conditions, which suggests that candidate-gene SNPs can be used to forecast the likely destiny of natural forest ecosystems under climate change scenarios. Differential levels of forest decline are anticipated for distinct maritime pine gene pools. Geographically defined molecular proxies for climate adaptation will thus critically enhance the predictive power of range-shift models and help establish mitigation measures for long-lived keystone forest trees in the face of impending climate change.

Evolutionary trajectories of primate genes involved in HIV pathogenesis.

The current availability of five complete genomes of different primate species allows the analysis of genetic divergence over the last 40 million years of evolution. We hypothesized that the interspecies differences observed in susceptibility to HIV-1 would be influenced by the long-range selective pressures on host genes associated with HIV-1 pathogenesis. We established a list of human genes (n = 140) proposed to be involved in HIV-1 biology and pathogenesis and a control set of 100 random genes. We retrieved the orthologous genes from the genome of humans and of four nonhuman primates (Pan troglodytes, Pongo pygmaeus abeli, Macaca mulatta, and Callithrix jacchus) and analyzed the nucleotide substitution patterns of this data set using codon-based maximum likelihood procedures. In addition, we evaluated whether the candidate genes have been targets of recent positive selection in humans by analyzing HapMap Phase 2 single-nucleotide polymorphisms genotyped in a region centered on each candidate gene. A total of 1,064 sequences were used for the analyses. Similar median K(A)/K(S) values were estimated for the set of genes involved in HIV-1 pathogenesis and for control genes, 0.19 and 0.15, respectively. However, genes of the innate immunity had median values of 0.37 (P value = 0.0001, compared with control genes), and genes of intrinsic cellular defense had K(A)/K(S) values around or greater than 1.0 (P value = 0.0002). Detailed assessment allowed the identification of residues under positive selection in 13 proteins: AKT1, APOBEC3G, APOBEC3H, CD4, DEFB1, GML, IL4, IL8RA, L-SIGN/CLEC4M, PTPRC/CD45, Tetherin/BST2, TLR7, and TRIM5alpha. A number of those residues are relevant for HIV-1 biology. The set of 140 genes involved in HIV-1 pathogenesis did not show a significant enrichment in signals of recent positive selection in humans (intraspecies selection). However, we identified within or near these genes 24 polymorphisms showing strong signatures of recent positive selection. Interestingly, the DEFB1 gene presented signatures of both interspecies positive selection in primates and intraspecies recent positive selection in humans. The systematic assessment of long-acting selective pressures on primate genomes is a useful tool to extend our understanding of genetic variation influencing contemporary susceptibility to HIV-1.

Common variants in the ATP2B1 gene are associated with susceptibility to hypertension: the Japanese Millennium Genome Project.

Hypertension is one of the most common complex genetic disorders. We have described previously 38 single nucleotide polymorphisms (SNPs) with suggestive association with hypertension in Japanese individuals. In this study we extend our previous findings by analyzing a large sample of Japanese individuals (n=14 105) for the most associated SNPs. We also conducted replication analyses in Japanese of susceptibility loci for hypertension identified recently from genome-wide association studies of European ancestries. Association analysis revealed significant association of the ATP2B1 rs2070759 polymorphism with hypertension (P=5.3×10(-5); allelic odds ratio: 1.17 [95% CI: 1.09 to 1.26]). Additional SNPs in ATP2B1 were subsequently genotyped, and the most significant association was with rs11105378 (odds ratio: 1.31 [95% CI: 1.21 to 1.42]; P=4.1×10(-11)). Association of rs11105378 with hypertension was cross-validated by replication analysis with the Global Blood Pressure Genetics consortium data set (odds ratio: 1.13 [95% CI: 1.05 to 1.21]; P=5.9×10(-4)). Mean adjusted systolic blood pressure was highly significantly associated with the same SNP in a meta-analysis with individuals of European descent (P=1.4×10(-18)). ATP2B1 mRNA expression levels in umbilical artery smooth muscle cells were found to be significantly different among rs11105378 genotypes. Seven SNPs discovered in published genome-wide association studies were also genotyped in the Japanese population. In the combined analysis with replicated 3 genes, FGF5 rs1458038, CYP17A1, rs1004467, and CSK rs1378942, odds ratio of the highest risk group was 2.27 (95% CI: 1.65 to 3.12; P=4.6×10(-7)) compared with the lower risk group. In summary, this study confirmed common genetic variation in ATP2B1, as well as FGF5, CYP17A1, and CSK, to be associated with blood pressure levels and risk of hypertension.

The recent breakthroughs in the understanding of host genomics in hepatitis C.

BACKGROUND: Hepatitis C Virus (HCV) infection is spontaneously resolved in about 30% of acutely infected individuals. In those who progress to chronic hepatitis C, HCV therapy permanently eradicates infection in about 40% of cases. It has long been suspected that host genetic factors are key determinants for the control of HCV infection. DESIGN: We will review in this study four genome-wide association studies (GWAS) and two large candidate gene studies that assessed the role of host genetic variation for the natural and treatment-induced control of HCV infection. RESULTS: The studies consistently identified genetic variation in interleukin 28B (IL28B) as the strongest predictor for the control of HCV infection. Importantly, single nucleotide polymorphisms (SNPs) in IL28B strongly predicted both spontaneous and treatment-induced HCV recovery. IL28B is located on chromosome 19 and encodes interferon-λ, a type III interferon with antiviral activity, which is mediated through the JAK-STAT pathway by inducing interferon-stimulated genes. The SNPs identified in the GWAS are in high linkage disequilibrium with coding or functional non-coding SNPs that might modulate function and/or expression of IL28B. The role of the different IL28B alleles on gene expression and cytokine function has not yet been established. CONCLUSIONS: These findings provide strong genetic evidence for the influence of interferon-λ for both the natural and treatment-induced control of HCV infection, and support the further investigation of interferon-λ for the treatment of chronic hepatitis C. Furthermore, genetic testing before HCV therapy could provide important information towards an individualized HCV treatment.

Depressive disorder moderates the effect of the FTO gene on body mass index.

There is evidence that obesity-related disorders are increased among people with depression. Variation in the FTO (fat mass and obesity associated) gene has been shown to contribute to common forms of human obesity. This study aimed to investigate the genetic influence of polymorphisms in FTO in relation to body mass index (BMI) in two independent samples of major depressive disorder (MDD) cases and controls. We analysed 88 polymorphisms in the FTO gene in a clinically ascertained sample of 2442 MDD cases and 809 controls (Radiant Study). In all, 8 of the top 10 single-nucleotide polymorphisms (SNPs) showing the strongest associations with BMI were followed-up in a population-based cohort (PsyCoLaus Study) consisting of 1292 depression cases and 1690 controls. Linear regression analyses of the FTO variants and BMI yielded 10 SNPs significantly associated with increased BMI in the depressive group but not the control group in the Radiant sample. The same pattern was found in the PsyCoLaus sample. We found a significant interaction between genotype and affected status in relation to BMI for seven SNPs in Radiant (P<0.0057), with PsyCoLaus giving supportive evidence for five SNPs (P-values between 0.03 and 0.06), which increased in significance when the data were combined in a meta-analysis. This is the first study investigating FTO and BMI within the context of MDD, and the results indicate that having a history of depression moderates the effect of FTO on BMI. This finding suggests that FTO is involved in the mechanism underlying the association between mood disorders and obesity.

Host genetic determinants of spontaneous hepatitis C clearance

Acute infection with the hepatitis C virus (HCV) induces a wide range of innate and adaptive immune responses. A total of 20-50% of acutely HCV-infected individuals permanently control the virus, referred to as 'spontaneous hepatitis C clearance', while the infection progresses to chronic hepatitis C in the majority of cases. Numerous studies have examined host genetic determinants of hepatitis C infection outcome and revealed the influence of genetic polymorphisms of human leukocyte antigens, killer immunoglobulin-like receptors, chemokines, interleukins and interferon-stimulated genes on spontaneous hepatitis C clearance. However, most genetic associations were not confirmed in independent cohorts, revealed opposing results in diverse populations or were limited by varying definitions of hepatitis C outcomes or small sample size. Coordinated efforts are needed in the search for key genetic determinants of spontaneous hepatitis C clearance that include well-conducted candidate genetic and genome-wide association studies, direct sequencing and follow-up functional studies.

Familiality and SNP heritability of age at onset and episodicity in major depressive disorder.

BACKGROUND: Strategies to dissect phenotypic and genetic heterogeneity of major depressive disorder (MDD) have mainly relied on subphenotypes, such as age at onset (AAO) and recurrence/episodicity. Yet, evidence on whether these subphenotypes are familial or heritable is scarce. The aims of this study are to investigate the familiality of AAO and episode frequency in MDD and to assess the proportion of their variance explained by common single nucleotide polymorphisms (SNP heritability). METHOD: For investigating familiality, we used 691 families with 2-5 full siblings with recurrent MDD from the DeNt study. We fitted (square root) AAO and episode count in a linear and a negative binomial mixed model, respectively, with family as random effect and adjusting for sex, age and center. The strength of familiality was assessed with intraclass correlation coefficients (ICC). For estimating SNP heritabilities, we used 3468 unrelated MDD cases from the RADIANT and GSK Munich studies. After similarly adjusting for covariates, derived residuals were used with the GREML method in GCTA (genome-wide complex trait analysis) software. RESULTS: Significant familial clustering was found for both AAO (ICC = 0.28) and episodicity (ICC = 0.07). We calculated from respective ICC estimates the maximal additive heritability of AAO (0.56) and episodicity (0.15). SNP heritability of AAO was 0.17 (p = 0.04); analysis was underpowered for calculating SNP heritability of episodicity. CONCLUSIONS: AAO and episodicity aggregate in families to a moderate and small degree, respectively. AAO is under stronger additive genetic control than episodicity. Larger samples are needed to calculate the SNP heritability of episodicity. The described statistical framework could be useful in future analyses.

Bioinformatics-driven identification and examination of candidate genes for non-alcoholic fatty liver disease.

ObjectiveCandidate genes for non-alcoholic fatty liver disease (NAFLD) identified by a bioinformatics approach were examined for variant associations to quantitative traits of NAFLD-related phenotypes.Research Design and MethodsBy integrating public database text mining, trans-organism protein-protein interaction transferal, and information on liver protein expression a protein-protein interaction network was constructed and from this a smaller isolated interactome was identified. Five genes from this interactome were selected for genetic analysis. Twenty-one tag single-nucleotide polymorphisms (SNPs) which captured all common variation in these genes were genotyped in 10,196 Danes, and analyzed for association with NAFLD-related quantitative traits, type 2 diabetes (T2D), central obesity, and WHO-defined metabolic syndrome (MetS).Results273 genes were included in the protein-protein interaction analysis and EHHADH, ECHS1, HADHA, HADHB, and ACADL were selected for further examination. A total of 10 nominal statistical significant associations (P<0.05) to quantitative metabolic traits were identified. Also, the case-control study showed associations between variation in the five genes and T2D, central obesity, and MetS, respectively. Bonferroni adjustments for multiple testing negated all associations.ConclusionsUsing a bioinformatics approach we identified five candidate genes for NAFLD. However, we failed to provide evidence of associations with major effects between SNPs in these five genes and NAFLD-related quantitative traits, T2D, central obesity, and MetS.

Biological, clinical and population relevance of 95 loci for blood lipids.

Plasma concentrations of total cholesterol, low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides are among the most important risk factors for coronary artery disease (CAD) and are targets for therapeutic intervention. We screened the genome for common variants associated with plasma lipids in >100,000 individuals of European ancestry. Here we report 95 significantly associated loci (P < 5 x 10(-8)), with 59 showing genome-wide significant association with lipid traits for the first time. The newly reported associations include single nucleotide polymorphisms (SNPs) near known lipid regulators (for example, CYP7A1, NPC1L1 and SCARB1) as well as in scores of loci not previously implicated in lipoprotein metabolism. The 95 loci contribute not only to normal variation in lipid traits but also to extreme lipid phenotypes and have an impact on lipid traits in three non-European populations (East Asians, South Asians and African Americans). Our results identify several novel loci associated with plasma lipids that are also associated with CAD. Finally, we validated three of the novel genes-GALNT2, PPP1R3B and TTC39B-with experiments in mouse models. Taken together, our findings provide the foundation to develop a broader biological understanding of lipoprotein metabolism and to identify new therapeutic opportunities for the prevention of CAD.

Common variants in mendelian kidney disease genes and their association with renal function.

Many common genetic variants identified by genome-wide association studies for complex traits map to genes previously linked to rare inherited Mendelian disorders. A systematic analysis of common single-nucleotide polymorphisms (SNPs) in genes responsible for Mendelian diseases with kidney phenotypes has not been performed. We thus developed a comprehensive database of genes for Mendelian kidney conditions and evaluated the association between common genetic variants within these genes and kidney function in the general population. Using the Online Mendelian Inheritance in Man database, we identified 731 unique disease entries related to specific renal search terms and confirmed a kidney phenotype in 218 of these entries, corresponding to mutations in 258 genes. We interrogated common SNPs (minor allele frequency >5%) within these genes for association with the estimated GFR in 74,354 European-ancestry participants from the CKDGen Consortium. However, the top four candidate SNPs (rs6433115 at LRP2, rs1050700 at TSC1, rs249942 at PALB2, and rs9827843 at ROBO2) did not achieve significance in a stage 2 meta-analysis performed in 56,246 additional independent individuals, indicating that these common SNPs are not associated with estimated GFR. The effect of less common or rare variants in these genes on kidney function in the general population and disease-specific cohorts requires further research.

Dyslipidemia in HIV-infected individuals: from pharmacogenetics to pharmacogenomics.

HIV-infected individuals may have accelerated atherogenesis and an increased risk for premature coronary artery disease. Dyslipidemia represents a key pro-atherogenic mechanism. In HIV-infected patients, dyslipidemia is typically attributed to the adverse effects of antiretroviral therapy. Nine recent genome-wide association studies have afforded a comprehensive, unbiased inventory of common SNPs at 36 genetic loci that are reproducibly associated with dyslipidemia in the general population. Genome-wide association study-validated SNPs have now been demonstrated to contribute to dyslipidemia in the setting of HIV infection and antiretroviral therapy. In a Swiss HIV-infected study population, a similar proportion of serum lipid variability was explained by antiretroviral therapy and by genetic background. In the individual patient, both antiretroviral therapy and the cumulative effect of SNPs contribute to the risk of high low-density lipoprotein cholesterol, low high-density lipoprotein cholesterol and hypertriglyceridemia. Genetic variants presumably contribute to additional major metabolic complications in HIV-infected individuals, including diabetes mellitus and coronary artery disease. In an effort to explain an increasing proportion of the heritability of complex metabolic traits, ongoing large-scale gene resequencing studies are focusing on the effects of rare SNPs and structural genetic variants.

Genetic vulnerability factor for schizophrenia: association with glutamate cysteine ligase modifier gene and abnormal enzyme activity

Background: We previously reported in schizophrenia patients a decreased level of glutathione ([GSH]), the principal non-protein antioxidant and redox regulator, both in cerebrospinal-fluid and prefrontal cortex. To identify possible genetic causation, we studied genes involved in GSH metabolism. Methods: Genotyping: mass spectrometry analysis of polymerase chain reaction (PCR) amplified DNA fragments purified from peripheral blood. Gene expression: real-time PCR of total RNA isolated from fibroblast cultures derived from skin of patients (DSM-IV) and healthy controls (DIGS). Results: Case-control association study of single nucleotide polymorphisms (SNP) from the GSH key synthesizing enzyme glutamate-cysteine-ligase (GCL) modifier subunit (GCLM) was performed in two populations: Swiss (patients/controls: 40/31) and Danish (349/348). We found a strong association of SNP rs2301022 in GCLM gene (Danish: c2=3.2; P=0.001 after correction for multiple testing). Evidence for GCLM as a risk factor was confirmed in linkage study of NIMH families. Moreover, we observed a decrease in GCLM mRNA levels in patient fibroblasts, consistently with the association study. Interestingly, Dalton and collaborators reported in GCLM knock-out mice an increased feedback inhibition of GCL activity, resulting in 60% decrease of brain [GSH], a situation analogous to patients. These mice also exhibited an increased sensitivity to oxidative stress. Similarly, under oxidative stress conditions, GCL enzymatic activity was also decreased in patient fibroblasts. Conclusions: These results at the genetic and functional levels, combined with observations that GSH deficient models reveal morphological, electrophysiological, and behavioral anomalies analogous to those observed in patients, suggest that GCLM allelic variant is a vulnerability factor for schizophrenia.

Novel Approach Identifies SNPs in SLC2A10 and KCNK9 with Evidence for Parent-of-Origin Effect on Body Mass Index.

The phenotypic effect of some single nucleotide polymorphisms (SNPs) depends on their parental origin. We present a novel approach to detect parent-of-origin effects (POEs) in genome-wide genotype data of unrelated individuals. The method exploits increased phenotypic variance in the heterozygous genotype group relative to the homozygous groups. We applied the method to >56,000 unrelated individuals to search for POEs influencing body mass index (BMI). Six lead SNPs were carried forward for replication in five family-based studies (of ∼4,000 trios). Two SNPs replicated: the paternal rs2471083-C allele (located near the imprinted KCNK9 gene) and the paternal rs3091869-T allele (located near the SLC2A10 gene) increased BMI equally (beta = 0.11 (SD), P<0.0027) compared to the respective maternal alleles. Real-time PCR experiments of lymphoblastoid cell lines from the CEPH families showed that expression of both genes was dependent on parental origin of the SNPs alleles (P<0.01). Our scheme opens new opportunities to exploit GWAS data of unrelated individuals to identify POEs and demonstrates that they play an important role in adult obesity.