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.

MMP13 mutations are the cause of recessive metaphyseal dysplasia, Spahr type.

Metaphyseal dysplasia, Spahr type (MDST; OMIM 250400) was described in 1961 based on the observation of four children in one family who had rickets-like metaphyseal changes but normal blood chemistry and moderate short stature. Its molecular basis and nosologic status remained unknown. We followed up on those individuals and diagnosed the disorder in an additional member of the family. We used exome sequencing to ascertain the underlying mutation and explored its consequences on three-dimensional models of the affected protein. The MDST phenotype is associated with moderate short stature and knee pain in adults, while extra-skeletal complications are not observed. The sequencing showed that MDST segregated with a c.619T>G single nucleotide transversion in MMP13. The predicted non-conservative amino acid substitution, p.Trp207Gly, disrupts a crucial hydrogen bond in the calcium-binding region of the catalytic domain of the matrix metalloproteinase, MMP13. The MDST phenotype is associated with recessive MMP13 mutations, confirming the importance of this metalloproteinase in the metaphyseal growth plate. Dominant MMP13 mutations have been associated with metaphyseal anadysplasia (OMIM 602111), while a single child homozygous for a MMP13 mutation had been previously diagnosed as "recessive metaphyseal anadysplasia," that we conclude is the same nosologic entity as MDST. Molecular confirmation of MDST allows distinction of it from dominant conditions (e.g., metaphyseal dysplasia, Schmid type; OMIM # 156500) and from more severe multi-system conditions (such as cartilage-hair hypoplasia; OMIM # 250250) and to give precise recurrence risks and prognosis. © 2014 Wiley Periodicals, Inc.

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.

Molecular mechanisms of penicillin-induced killing and penicillin tolerance in Streptococcus gordonii

Abstract The main thesis topic relates to the 'molecular mechanisms of penicillin-induced bacterial death. Indeed, bacteria have developed two principal mechanisms to escape the killing effect of ß-lactam antibiotics: resistance and tolerance. Resistant bacteria are characterized by their ability to grow in the presence of drug concentrations higher than the one inhibiting the growth of susceptible members of the same species. Hence, resistant bacteria have an increased minimal inhibitory concentration (MIC) of the drug. Nevertheless, when exposed to antibiotic concentrations exceeding their new MIC, resistant bacteria remain sensitive to the antibiotic killing effect. In contrast, tolerant bacteria have an unchanged MIC. However, they have a considerably increased ability to survive drug-induced killing, even at concentrations exceeding their MIC by several orders of magnitude. In other words, in the presence of the antibiotic, tolerant bacteria become persister cells which stop growing but are not killed. In the present thesis, it is shown that the survival phenotype of a tolerant Streptococcus gordonii strain depends on two components belonging to sugar metabolism pathways. First, the transcription factor CcpA which mediates a global regulatory mechanism allowing bacteria to utilize the most efficient sugar source for their growth. We show that the inactivation of the ccpA gene leads to a partial loss of penicillin tolerance both in vitro and in a rat model of experimental endocarditis. Second, the Enzyme I of the phosphotransferase system which is involved in the uptake and phosphorylation of sugars. Here, we -show that a single nucleotide mutation in ptsI, the gene encoding the Enzyme I, is sufficient to confer a fully tolerant phenotype in S. gordonii both in vivo and in vivo. The mutation results in a radical proline to arginine substitution in the C-terminal domain of the protein, probably leading to a decrease in its homodimerization and subsequent activity. Taken together our results prove that tolerance is a global survival mechanism linked to sugar metabolism. We hypothesize that, in the presence of the antibiotic, the already altered metabolic processes of the tolerant strain are completely inactivated. Hence, bacteria may enter in a dormant state and become insensitive to the bactericidal effect of ß-lactams, which depends on actively dividing cells. This thesis manuscript also contains two other side-projects. The first one establishes that the ability to form a biofilm is not a requisite for the successful establishment of endocarditis due to S. gordonii. The second one characterizes the S. gordonii a-phosphoglucomutase gene, and shows that its inactivation results in a loss of in vitro fitness and in vivo virulence. Résumé Le sujet principal de cette thèse concerne les mécanismes moléculaires de la mort bactérienne induite par la pénicilline. En effet, les bactéries ont développé deux mécanismes principaux pour échapper à l'effet bactéricide des ß-lactamines : la résistance et la tolérance. Les bactéries résistantes sont caractérisées par leur capacité de croître en présence de concentration d'antibiotiques plus élevées que celles inhibant la croissance des organismes sensibles de la même espèce. Les bactéries résistantes ont donc une augmentation de leur concentration minimale inhibitrice (CMI) à l'antibiotique. Néanmoins, quand elles sont exposées à des concentrations dépassant leur nouvelle CMI, elles restent sensibles à l'effet bactéricide. Au contraire, les bactéries tolérantes ont une CMI inchangée. Toutefois, elles ont une très importante capacité à survivre à l'effet bactéricide des ß-lactamines, ceci même à des concentrations excédant leur CMI de plusieurs ordres de grandeur. En d'autres termes, en présence de l'antibiotique, les bactéries tolérantes deviennent des cellules persistantes qui arrêtent leur croissance mais ne sont pas tuées. Dans la présente thèse, il est montré que le phénotype de survie d'un Streptococcus gordonii tolérant dépend de deux composants appartenant aux voies du métabolisme des sucres. Premièrement, le facteur de transcription CcpA qui contrôle un système global de régulation permettant à la bactérie d'utiliser les sources de sucre les plus efficaces pour sa croissance. Il est montré que l'inactivation du gène ccpA résulte en la perte partielle de la tolérance à la pénicilline aussi bien in vitro que dans un modèle d'endocardite expérimentale chez le rat. Deuxièmement, l'Enzyme I du système de phosphotransfert impliqué dans l'import et la phosphorylation des sucres. Nous montrons qu'une mutation ponctuelle d'un nucléotide dans ptsl, le gène codant pour l'Enzyme I, suffit à complètement conférer un phénotype tolérant chez S. gordonii aussi bien in vitro qu'in vivo. La mutation induit la substitution radicale d'une proline en une arginine dans le domaine C-terminal de la protéine, résultant probablement en une diminution de sa capacité d'homodimérisation et donc d'activité. Dans leur ensemble, nos résultats prouvent que la tolérance est un mécanisme global de survie lié au métabolisme des sucres. Nous présentons l'hypothèse que, en présence de l'antibiotique, les processus métaboliques déjà altérés de la souche tolérante deviennent complètement inactifs. En conséquence, les bactéries entreraient dans un état dormant nonréplicatif, devenant ainsi insensibles à l'effet bactéricide des ß-lactamines qui nécessite des cellules en cours de division active. Le manuscrit de cette thèse contient également deux projets secondaires. Le premier montre que la capacité de former un biofilm n'est pas un prérequis pour le succès de l'initiation de l'endocardite à S. gordonii. Le second caractérise le gène de l'a-phosphoglucomutase de S. gordonii et montre que son inactivation résulte en une perte de fitness in vitro et de virulence in vivo.

Monitoring of malaria parasite resistance to chloroquine and sulphadoxine-pyrimethamine in the Solomon Islands by DNA microarray technology.

BACKGROUND: Little information is available on resistance to anti-malarial drugs in the Solomon Islands (SI). The analysis of single nucleotide polymorphisms (SNPs) in drug resistance associated parasite genes is a potential alternative to classical time- and resource-consuming in vivo studies to monitor drug resistance. Mutations in pfmdr1 and pfcrt were shown to indicate chloroquine (CQ) resistance, mutations in pfdhfr and pfdhps indicate sulphadoxine-pyrimethamine (SP) resistance, and mutations in pfATPase6 indicate resistance to artemisinin derivatives. METHODS: The relationship between the rate of treatment failure among 25 symptomatic Plasmodium falciparum-infected patients presenting at the clinic and the pattern of resistance-associated SNPs in P. falciparum infecting 76 asymptomatic individuals from the surrounding population was investigated. The study was conducted in the SI in 2004. Patients presenting at a local clinic with microscopically confirmed P. falciparum malaria were recruited and treated with CQ+SP. Rates of treatment failure were estimated during a 28-day follow-up period. In parallel, a DNA microarray technology was used to analyse mutations associated with CQ, SP, and artemisinin derivative resistance among samples from the asymptomatic community. Mutation and haplotype frequencies were determined, as well as the multiplicity of infection. RESULTS: The in vivo study showed an efficacy of 88% for CQ+SP to treat P. falciparum infections. DNA microarray analyses indicated a low diversity in the parasite population with one major haplotype present in 98.7% of the cases. It was composed of fixed mutations at position 86 in pfmdr1, positions 72, 75, 76, 220, 326 and 356 in pfcrt, and positions 59 and 108 in pfdhfr. No mutation was observed in pfdhps or in pfATPase6. The mean multiplicity of infection was 1.39. CONCLUSION: This work provides the first insight into drug resistance markers of P. falciparum in the SI. The obtained results indicated the presence of a very homogenous P. falciparum population circulating in the community. Although CQ+SP could still clear most infections, seven fixed mutations associated with CQ resistance and two fixed mutations related to SP resistance were observed. Whether the absence of mutations in pfATPase6 indicates the efficacy of artemisinin derivatives remains to be proven.

SNP2GO: functional analysis of genome-wide association studies.

Genome-wide association studies (GWAS) are designed to identify the portion of single-nucleotide polymorphisms (SNPs) in genome sequences associated with a complex trait. Strategies based on the gene list enrichment concept are currently applied for the functional analysis of GWAS, according to which a significant overrepresentation of candidate genes associated with a biological pathway is used as a proxy to infer overrepresentation of candidate SNPs in the pathway. Here we show that such inference is not always valid and introduce the program SNP2GO, which implements a new method to properly test for the overrepresentation of candidate SNPs in biological pathways.

Human leukocyte antigen class II variants and adult-onset asthma: does occupational allergen exposure play a role?

Recently, a locus centred on rs9273349 in the HLA-DQ region emerged from genome-wide association studies of adult-onset asthma. We aimed to further investigate the role of human leukocyte antigen (HLA) class II in adult-onset asthma and a possible interaction with occupational exposures. We imputed classical HLA-II alleles from 7579 single-nucleotide polymorphisms in 6025 subjects (1202 with adult-onset asthma) from European cohorts: ECRHS, SAPALDIA, EGEA and B58C, and from surveys of bakers and agricultural workers. Based on an asthma-specific job-exposure matrix, 2629 subjects had ever been exposed to high molecular weight (HMW) allergens. We explored associations between 23 common HLA-II alleles and adult-onset asthma, and tested for gene-environment interaction with occupational exposure to HMW allergens. Interaction was also tested for rs9273349. Marginal associations of classical HLA-II alleles and adult-onset asthma were not statistically significant. Interaction was detected between the DPB1*03:01 allele and exposure to HMW allergens (p = 0.009), in particular to latex (p = 0.01). In the unexposed group, the DPB1*03:01 allele was associated with adult-onset asthma (OR 0.67, 95%CI 0.53-0.86). HMW allergen exposures did not modify the association of rs9273349 with adult-onset asthma. Common classical HLA-II alleles were not marginally associated with adult-onset asthma. The association of latex exposure and adult-onset asthma may be modified by DPB1*03:01.

A common biological basis of obesity and nicotine addiction.

Smoking influences body weight such that smokers weigh less than non-smokers and smoking cessation often leads to weight increase. The relationship between body weight and smoking is partly explained by the effect of nicotine on appetite and metabolism. However, the brain reward system is involved in the control of the intake of both food and tobacco. We evaluated the effect of single-nucleotide polymorphisms (SNPs) affecting body mass index (BMI) on smoking behavior, and tested the 32 SNPs identified in a meta-analysis for association with two smoking phenotypes, smoking initiation (SI) and the number of cigarettes smoked per day (CPD) in an Icelandic sample (N=34,216 smokers). Combined according to their effect on BMI, the SNPs correlate with both SI (r=0.019, P=0.00054) and CPD (r=0.032, P=8.0 × 10(-7)). These findings replicate in a second large data set (N=127,274, thereof 76,242 smokers) for both SI (P=1.2 × 10(-5)) and CPD (P=9.3 × 10(-5)). Notably, the variant most strongly associated with BMI (rs1558902-A in FTO) did not associate with smoking behavior. The association with smoking behavior is not due to the effect of the SNPs on BMI. Our results strongly point to a common biological basis of the regulation of our appetite for tobacco and food, and thus the vulnerability to nicotine addiction and obesity.

Genetic Association of Recovery from Eating Disorders: The Role of GABA Receptor SNPs.

Follow-up studies of eating disorders (EDs) suggest outcomes ranging from recovery to chronic illness or death, but predictors of outcome have not been consistently identified. We tested 5151 single-nucleotide polymorphisms (SNPs) in approximately 350 candidate genes for association with recovery from ED in 1878 women. Initial analyses focused on a strictly defined discovery cohort of women who were over age 25 years, carried a lifetime diagnosis of an ED, and for whom data were available regarding the presence (n=361 ongoing symptoms in the past year, ie, 'ill') or absence (n=115 no symptoms in the past year, ie, 'recovered') of ED symptoms. An intronic SNP (rs17536211) in GABRG1 showed the strongest statistical evidence of association (p=4.63 × 10(-6), false discovery rate (FDR)=0.021, odds ratio (OR)=0.46). We replicated these findings in a more liberally defined cohort of women age 25 years or younger (n=464 ill, n=107 recovered; p=0.0336, OR=0.68; combined sample p=4.57 × 10(-6), FDR=0.0049, OR=0.55). Enrichment analyses revealed that GABA (γ-aminobutyric acid) SNPs were over-represented among SNPs associated at p<0.05 in both the discovery (Z=3.64, p=0.0003) and combined cohorts (Z=2.07, p=0.0388). In follow-up phenomic association analyses with a third independent cohort (n=154 ED cases, n=677 controls), rs17536211 was associated with trait anxiety (p=0.049), suggesting a possible mechanism through which this variant may influence ED outcome. These findings could provide new insights into the development of more effective interventions for the most treatment-resistant patients.

Common variants in UMOD associate with urinary uromodulin levels: a meta-analysis.

Uromodulin is expressed exclusively in the thick ascending limb and is the most abundant protein excreted in normal urine. Variants in UMOD, which encodes uromodulin, are associated with renal function, and urinary uromodulin levels may be a biomarker for kidney disease. However, the genetic factors regulating uromodulin excretion are unknown. We conducted a meta-analysis of urinary uromodulin levels to identify associated common genetic variants in the general population. We included 10,884 individuals of European descent from three genetic isolates and three urban cohorts. Each study measured uromodulin indexed to creatinine and conducted linear regression analysis of approximately 2.5 million single nucleotide polymorphisms using an additive model. We also tested whether variants in genes expressed in the thick ascending limb associate with uromodulin levels. rs12917707, located near UMOD and previously associated with renal function and CKD, had the strongest association with urinary uromodulin levels (P<0.001). In all cohorts, carriers of a G allele of this variant had higher uromodulin levels than noncarriers did (geometric means 10.24, 14.05, and 17.67 μg/g creatinine for zero, one, or two copies of the G allele). rs12446492 in the adjacent gene PDILT (protein disulfide isomerase-like, testis expressed) also reached genome-wide significance (P<0.001). Regarding genes expressed in the thick ascending limb, variants in KCNJ1, SORL1, and CAB39 associated with urinary uromodulin levels. These data indicate that common variants in the UMOD promoter region may influence urinary uromodulin levels. They also provide insights into uromodulin biology and the association of UMOD variants with renal function.

Genetic polymorphisms of the main transcription factors for adiponectin gene promoter in regulation of adiponectin levels: association analysis in three European cohorts.

Adiponectin serum concentrations are an important biomarker in cardiovascular epidemiology with heritability etimates of 30-70%. However, known genetic variants in the adiponectin gene locus (ADIPOQ) account for only 2%-8% of its variance. As transcription factors are thought to play an under-acknowledged role in carrying functional variants, we hypothesized that genetic polymorphisms in genes coding for the main transcription factors for the ADIPOQ promoter influence adiponectin levels. Single nucleotide polymorphisms (SNPs) at these genes were selected based on the haplotype block structure and previously published evidence to be associated with adiponectin levels. We performed association analyses of the 24 selected SNPs at forkhead box O1 (FOXO1), sterol-regulatory-element-binding transcription factor 1 (SREBF1), sirtuin 1 (SIRT1), peroxisome-proliferator-activated receptor gamma (PPARG) and transcription factor activating enhancer binding protein 2 beta (TFAP2B) gene loci with adiponectin levels in three different European cohorts: SAPHIR (n = 1742), KORA F3 (n = 1636) and CoLaus (n = 5355). In each study population, the association of SNPs with adiponectin levels on log-scale was tested using linear regression adjusted for age, sex and body mass index, applying both an additive and a recessive genetic model. A pooled effect size was obtained by meta-analysis assuming a fixed effects model. We applied a significance threshold of 0.0033 accounting for the multiple testing situation. A significant association was only found for variants within SREBF1 applying an additive genetic model (smallest p-value for rs1889018 on log(adiponectin) = 0.002, β on original scale = -0.217 µg/ml), explaining ∼0.4% of variation of adiponectin levels. Recessive genetic models or haplotype analyses of the FOXO1, SREBF1, SIRT1, TFAPB2B genes or sex-stratified analyses did not reveal additional information on the regulation of adiponectin levels. The role of genetic variations at the SREBF1 gene in regulating adiponectin needs further investigation by functional studies.

Molecular substitution rate increases in myrmecophilous lycaenid butterflies (Lepidoptera)

Is species diversification driven by neutral- or niche-based processes? Butterflies of the Lycaenidae family have developed mutualistic interactions with ants. This biotic requirement increased the spatial fragmentation of populations of lower effective population size (Ne) compared with autonomous species. The nearly neutral theory predicts that species with smaller Ne should fix more mutations because of the increased strength of drift. Taking into account the phylogenetic relatedness among species, this study shows that species with a stronger dependence on ants displayed more intra-specific Single Nucleotide Polymorphisms compared with species with low or no myrmecophily. This phenomenon can cause more pronounced genetic differentiation between populations and could ultimately promote speciation in a similar manner as on physical islands. The large species diversity observed in this family could be the consequence of this neutral process enhancing the diversification of lineages.

Viral genotype-specific role of PNPLA3, PPARG, MTTP, and IL28B in hepatitis C virus-associated steatosis.

BACKGROUND & AIMS: Steatosis is a prominent feature of hepatitis C, especially in patients infected with genotype 3. The analysis of genetic polymorphisms influencing steatosis in chronic hepatitis C has been limited by the studies' small sample size, and important single nucleotide polymorphisms (SNPs), such as those in the patatin-like phospholipase family 3 protein (PNPLA3), were never evaluated. METHODS: We analyzed the role of SNPs, from 19 systematically selected candidate genes, on steatosis in 626 Caucasian hepatitis C virus (HCV) infected patients. SNPs were extracted from a genome-wide association-generated dataset. Associations of alleles with the presence and/or different severity of steatosis were evaluated by univariate and multivariate logistic regression, accounting for all relevant covariates. RESULTS: The risk of steatosis was increased by carriage of I148M in PNPLA3, but only in patients with HCV genotypes non-3 (odds ratio [OR]=1.9, 95% confidence interval [CI]=1.6-2.3, p<0.001) and similar, albeit weaker associations were found for SNPs in peroxisome proliferator-activated receptor-γ (PPARG) and interleukin-28B (IL28B). Carriage of a SNP in the microsomal triglyceride transfer protein (MTTP) increased the risk of steatosis, but only in patients with HCV genotype 3 (rs1800803, OR=3.4, 95% CI=2.4-4.9, p=0.001). CONCLUSIONS: The rs738409 SNP in PNPLA3 is associated with an increased risk of steatosis in patients infected with HCV genotypes non-3. Host genes affect steatosis depending on the infecting HCV genotype, suggesting their interaction with viral factors.

Pharmacogenetics-based population pharmacokinetic analysis of etravirine in HIV-1 infected individuals.

OBJECTIVES: Etravirine (ETV) is metabolized by cytochrome P450 (CYP) 3A, 2C9, and 2C19. Metabolites are glucuronidated by uridine diphosphate glucuronosyltransferases (UGT). To identify the potential impact of genetic and non-genetic factors involved in ETV metabolism, we carried out a two-step pharmacogenetics-based population pharmacokinetic study in HIV-1 infected individuals. MATERIALS AND METHODS: The study population included 144 individuals contributing 289 ETV plasma concentrations and four individuals contributing 23 ETV plasma concentrations collected in a rich sampling design. Genetic variants [n=125 single-nucleotide polymorphisms (SNPs)] in 34 genes with a predicted role in ETV metabolism were selected. A first step population pharmacokinetic model included non-genetic and known genetic factors (seven SNPs in CYP2C, one SNP in CYP3A5) as covariates. Post-hoc individual ETV clearance (CL) was used in a second (discovery) step, in which the effect of the remaining 98 SNPs in CYP3A, P450 cytochrome oxidoreductase (POR), nuclear receptor genes, and UGTs was investigated. RESULTS: A one-compartment model with zero-order absorption best characterized ETV pharmacokinetics. The average ETV CL was 41 (l/h) (CV 51.1%), the volume of distribution was 1325 l, and the mean absorption time was 1.2 h. The administration of darunavir/ritonavir or tenofovir was the only non-genetic covariate influencing ETV CL significantly, resulting in a 40% [95% confidence interval (CI): 13-69%] and a 42% (95% CI: 17-68%) increase in ETV CL, respectively. Carriers of rs4244285 (CYP2C19*2) had 23% (8-38%) lower ETV CL. Co-administered antiretroviral agents and genetic factors explained 16% of the variance in ETV concentrations. None of the SNPs in the discovery step influenced ETV CL. CONCLUSION: ETV concentrations are highly variable, and co-administered antiretroviral agents and genetic factors explained only a modest part of the interindividual variability in ETV elimination. Opposing effects of interacting drugs effectively abrogate genetic influences on ETV CL, and vice-versa.