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.

Adaptive evolution during an ongoing range expansion: the invasive bank vole (Myodes glareolus) in Ireland

Range expansions are extremely common, but have only recently begun to attract attention in terms of their genetic consequences. As populations expand, demes at the wave front experience strong genetic drift, which is expected to reduce genetic diversity and potentially cause ‘allele surfing’, where alleles may become fixed over a wide geographical area even if their effects are deleterious. Previous simulation models show that range expansions can generate very strong selective gradients on dispersal, reproduction, competition and immunity. To investigate the effects of range expansion on genetic diversity and adaptation, we studied the population genomics of the bank vole (Myodes glareolus) in Ireland. The bank vole was likely introduced in the late 1920s and is expanding its range at a rate of ~2.5 km/year. Using genotyping-by-sequencing, we genotyped 281 bank voles at 5979 SNP loci. Fourteen sample sites were arranged in three transects running from the introduction site to the wave front of the expansion. We found significant declines in genetic diversity along all three transects. However, there was no evidence that sites at the wave front had accumulated more deleterious mutations. We looked for outlier loci with strong correlations between allele frequency and distance from the introduction site, where the direction of correlation was the same in all three transects. Amongst these outliers, we found significant enrichment for genic SNPs, suggesting the action of selection. Candidates for selection included several genes with immunological functions and several genes that could influence behaviour.

Associations between prefrontal γ-aminobutyric acid concentration and the tryptophan hydroxylase isoform 2 gene, a panic disorder risk allele in women

Associations between the central serotonergic and γ-aminobutyric acid (GABA) systems play key roles in the prefrontal cortical regulation of emotion and cognition and in the pathophysiology and pharmacotherapy of highly prevalent psychiatric disorders. The goal of this study was to test the effects of common variants of the tryptophan hydroxylase isoform 2 (TPH2) gene on GABA concentration in the prefrontal cortex (PFC) using magnetic resonance spectroscopy. In this study involving 64 individuals, we examined the associations between prefrontal cortical GABA concentration and 12 single nucleotide polymorphisms (SNPs) spanning the TPH2 gene, including rs4570625 (−703 G/T SNP), a potentially functional TPH2 polymorphism that has been associated with decreased TPH2 mRNA expression and panic disorder. Our results revealed a significant association between increased GABA concentration in the PFC and the T-allele frequencies of two TPH2 SNPs, namely rs4570625 (−703 G/T) and rs2129575 (p≤0.0004) and the C-allele frequency of one TPH2 SNP, namely rs1386491 (p = 0.0003) in female subjects. We concluded that rs4570625 (−703 G/T), rs2129575 and rs1386491 play a significant role in GABAergic neurotransmission and may contribute to the sex-specific dysfunction of the GABAergic system in the PFC.

A chromosomal region on ECA13 is associated with maxillary prognathism in horses.

Hereditary variations in head morphology and head malformations are known in many species. The most common variation encountered in horses is maxillary prognathism. Prognathism and brachygnathism are syndromes of the upper and lower jaw, respectively. The resulting malocclusion can negatively affect teeth wear, and is considered a non-desirable trait in breeding programs. We performed a case-control analysis for maxillary prognathism in horses using 96 cases and 763 controls. All horses had been previously genotyped with a commercially available 50 k SNP array. We analyzed the data with a mixed-model considering the genomic relationships in order to account for population stratification. Two SNPs within a region on the distal end of chromosome ECA 13 reached the Bonferroni corrected genome-wide significance level. There is no known prognathism candidate gene located within this region. Therefore, our findings in the horse offer the possibility of identifying a novel gene involved in the complex genetics of prognathism that might also be relevant for humans and other livestock species.

Analysis of Variation in Clubfoot Candidate Genes

Isolated clubfoot, a common birth defect occurring in more than 135,000 livebirths worldwide each year, is associated with significant health care and financial burdens. Clubfoot is defined by forefoot adduction, hindfoot varus, midfoot cavus and hindfoot equinus. Isolated clubfoot, which is the focus of these studies, is distinct from syndromic clubfoot because there are no other associated malformations. Population, family, twin and segregation analysis studies provide evidence that genetic and environmental factors play an etiologic role in isolated clubfoot. The studies described in this thesis were performed to define the role of genetic variation in isolated clubfoot. Interrogation of a deletion region associated with syndromic clubfoot, suggested that CASP8 and CASP10, two apoptotic genes, play a role in isolated clubfoot. To explore the role of apoptotic genes in clubfoot, SNPs spanning genes involved in the apoptotic pathway in the six chromosomal deletion regions, and limb patterning genes, HOXD and HOXA, were interrogated. SNPs in mitochondrial mediated apoptotic genes and several SNPs in HOXA and HOXD genes were modestly associated with clubfoot with the most significant SNP, rs3801776, located in the basal promoter of HOXA9. Several significant associations were found with SNPs in NFAT2 and TNIP2. Significant gene interactions were detected between SNPs in HOX and apoptotic genes. These findings suggest a model for clubfoot in which variation in one gene is not sufficient to cause the malformation but requires variation several genes to perturb protein expression sufficiently to alter muscle and foot development. These results significantly impact our knowledge base by delineating underlying mechanisms causing clubfoot.

Xenobiotic Metabolism Genes and Clubfoot

Idiopathic or isolated clubfoot is a common orthopedic birth defect that affects approximately 135,000 children worldwide. It is characterized by equinus, varus and adductus deformities of the ankle and foot. Correction of clubfoot involves months of serial manipulations, castings and bracing, with surgical correction needed in forty percent of cases. Multifactorial etiology has been suggested in numerous studies with both environmental and genetic factors playing an etiologic role. Maternal smoking during pregnancy is the only common environmental factor that has consistently been shown to increase the risk for clubfoot. Moreover, a positive family history of clubfoot and maternal smoking increases the risk of clubfoot twenty fold. These findings suggest that genetic variation in smoking metabolism genes may increase susceptibility to clubfoot. Based on this reasoning, we interrogated eight candidate genes, chosen based on their involvement in phase 1 and 2 cigarette smoke metabolism. Twenty-two SNPs and two null alleles in eight genes (CYP1A1, CYP1A2, CYP1B1, CYP2A6, EPHX1, NAT2, GSTM1 and GSTT1) were genotyped in a dataset composed of nonHispanic white and Hispanic multiplex and simplex families. Only one SNP in CYP1A1, rs1048943, had significantly altered transmission in the aggregate and multiplex NHW datasets (p=0.003 and p=0.009). Perturbation of CYP1A1 by rs1048943 polymorphism causes an increase in the amount of harmful, adduct forming metabolic intermediates. A significant gene interaction between EPHX1 and NAT2 was also found (p=0.007). This interaction may affect the metabolism of harmful metabolic intermediates. Additionally, marginal interactions were found for other xenobiotic genes and these interactions may play a contributory role in clubfoot. Importantly, for CYP1A2, significant maternal (p=0.03; RR=1.24; 95% CI: 1.04-1.44) and fetal (p=0.01; RR=1.33; 95% CI: 1.13-1.54) genotypic effects were identified suggesting that both maternal and fetal genotypes impact normal limb development. No association was found for maternal smoking status and tobacco metabolism genes. Together, these results suggest that xenobiotic metabolism genes may play a contributory role in the etiology of clubfoot regardless of maternal smoking status and may impact foot development through perturbation of tobacco metabolic pathways.

Reduced neutrophil count in people of African descent is due to a regulatory variant in the Duffy antigen receptor for chemokines gene.

Persistently low white blood cell count (WBC) and neutrophil count is a well-described phenomenon in persons of African ancestry, whose etiology remains unknown. We recently used admixture mapping to identify an approximately 1-megabase region on chromosome 1, where ancestry status (African or European) almost entirely accounted for the difference in WBC between African Americans and European Americans. To identify the specific genetic change responsible for this association, we analyzed genotype and phenotype data from 6,005 African Americans from the Jackson Heart Study (JHS), the Health, Aging and Body Composition (Health ABC) Study, and the Atherosclerosis Risk in Communities (ARIC) Study. We demonstrate that the causal variant must be at least 91% different in frequency between West Africans and European Americans. An excellent candidate is the Duffy Null polymorphism (SNP rs2814778 at chromosome 1q23.2), which is the only polymorphism in the region known to be so differentiated in frequency and is already known to protect against Plasmodium vivax malaria. We confirm that rs2814778 is predictive of WBC and neutrophil count in African Americans above beyond the previously described admixture association (P = 3.8 x 10(-5)), establishing a novel phenotype for this genetic variant.

Multiple independent genetic factors at NOS1AP modulate the QT interval in a multi-ethnic population.

Extremes of electrocardiographic QT interval are associated with increased risk for sudden cardiac death (SCD); thus, identification and characterization of genetic variants that modulate QT interval may elucidate the underlying etiology of SCD. Previous studies have revealed an association between a common genetic variant in NOS1AP and QT interval in populations of European ancestry, but this finding has not been extended to other ethnic populations. We sought to characterize the effects of NOS1AP genetic variants on QT interval in the multi-ethnic population-based Dallas Heart Study (DHS, n = 3,072). The SNP most strongly associated with QT interval in previous samples of European ancestry, rs16847548, was the most strongly associated in White (P = 0.005) and Black (P = 3.6 x 10(-5)) participants, with the same direction of effect in Hispanics (P = 0.17), and further showed a significant SNP x sex-interaction (P = 0.03). A second SNP, rs16856785, uncorrelated with rs16847548, was also associated with QT interval in Blacks (P = 0.01), with qualitatively similar results in Whites and Hispanics. In a previously genotyped cohort of 14,107 White individuals drawn from the combined Atherosclerotic Risk in Communities (ARIC) and Cardiovascular Health Study (CHS) cohorts, we validated both the second locus at rs16856785 (P = 7.63 x 10(-8)), as well as the sex-interaction with rs16847548 (P = 8.68 x 10(-6)). These data extend the association of genetic variants in NOS1AP with QT interval to a Black population, with similar trends, though not statistically significant at P<0.05, in Hispanics. In addition, we identify a strong sex-interaction and the presence of a second independent site within NOS1AP associated with the QT interval. These results highlight the consistent and complex role of NOS1AP genetic variants in modulating QT interval.

A genome-wide association study of resistance to HIV infection in highly exposed uninfected individuals with hemophilia A.

Human genetic variation contributes to differences in susceptibility to HIV-1 infection. To search for novel host resistance factors, we performed a genome-wide association study (GWAS) in hemophilia patients highly exposed to potentially contaminated factor VIII infusions. Individuals with hemophilia A and a documented history of factor VIII infusions before the introduction of viral inactivation procedures (1979-1984) were recruited from 36 hemophilia treatment centers (HTCs), and their genome-wide genetic variants were compared with those from matched HIV-infected individuals. Homozygous carriers of known CCR5 resistance mutations were excluded. Single nucleotide polymorphisms (SNPs) and inferred copy number variants (CNVs) were tested using logistic regression. In addition, we performed a pathway enrichment analysis, a heritability analysis, and a search for epistatic interactions with CCR5 Δ32 heterozygosity. A total of 560 HIV-uninfected cases were recruited: 36 (6.4%) were homozygous for CCR5 Δ32 or m303. After quality control and SNP imputation, we tested 1 081 435 SNPs and 3686 CNVs for association with HIV-1 serostatus in 431 cases and 765 HIV-infected controls. No SNP or CNV reached genome-wide significance. The additional analyses did not reveal any strong genetic effect. Highly exposed, yet uninfected hemophiliacs form an ideal study group to investigate host resistance factors. Using a genome-wide approach, we did not detect any significant associations between SNPs and HIV-1 susceptibility, indicating that common genetic variants of major effect are unlikely to explain the observed resistance phenotype in this population.

Gene Discovery in Nonsyndromic Cleft Lip with or without Cleft Palate

Nonsyndromic cleft lip with or without cleft palate (NSCLP), a common, complex orofacial birth defect that affects approximately 4,000 newborns each year in the United States, is caused by both genetic and environmental factors. Orofacial clefts affect the mouth and nose, causing severe deformity of the face, which require medical, dental and speech therapies. Despite having substantial genetic liability, less than 25% of the genetic contribute to NSCLP has been identified. The studies described in this thesis were performed to identify genes that contribute to NSCLP and to demonstrate the role of these genes in normal craniofacial development. Using genome scan and candidate gene approaches, novel associations with NSCLP were identified. These include MYH9 (7 SNPs, 0.009≤p<0.05), Wnt3A (4 SNPs, 0.001≤p≤0.005), Wnt11 (2 SNPs, 0.001≤p≤0.01) and CRISPLD2 (4 SNPs, 0.001≤p<0.05). The most interesting findings were for CRISPLD2. This gene is expressed in the fused mouse palate at E17.5. In zebrafish, crispld2 localized to the craniofacial region by one day post fertilization. Morpholino knockdown of crispld2 resulted in a lower survival rates and altered neural crest cell (NCC) clustering. Because NCCs form the tissues that populate the craniofacies, this NCC abnormality resulted in cartilage abnormalities of the jaw including fewer ceratobranchial cartilages forming the lower jaw (three pairs compared to five) and broader craniofacies compared to wild-type zebrafish. These findings suggest that the CRISPLD2 gene plays an important role in normal craniofacial development and perturbation of this gene in humans contributes to orofacial clefting. Overall, these results are important because they contribute to our understanding of normal craniofacial development and orofacial clefting etiology, information that can be used to develop better methods to diagnose, counsel and potentially treat NSCLP patients.

Upregulation of defensins in burn sheep small intestine.

OBJECTIVE: The aim of this study was to visualize and localize the sheep antimicrobials, beta-defensins 1, 2, and 3, (SBD-1, SBD-2, SBD-3), sheep neutrophil defensin alpha (SNP-1), and the cathelicidin LL-37 in sheep small intestine after burn injury, our hypothesis being that these compounds would be upregulated in an effort to overcome a compromised endothelial lining. Response to burn injury includes the release of proinflammatory cytokines and systemic immune suppression that, if untreated, can progress to multiple organ failure and death, so protective mechanisms have to be initiated and implemented. METHODS: Tissue sections were probed with antibodies to the antimicrobials and then visualized with fluorescently labeled secondary antibodies and subjected to fluorescence deconvolution microscopy and image reconstruction. RESULTS: In both the sham and burn samples, all the aforementioned antimicrobials were seen in each of the layers of small intestine, the highest concentration being localized to the epithelium. SBD-2, SBD-3, and SNP-1 were upregulated in both enterocytes and Paneth cells, while SNP-1 and LL-37 showed increases in both the inner circular and outer longitudinal muscle layers of the muscularis externa following burn injury. Each of the defensins, except SBD-1, was also seen in between the muscle layers of the externa and while burn caused slight increases of SBD-2, SBD-3, and SNP-1 in this location, LL-37 content was significantly decreased. CONCLUSION: That while each of these human antimicrobials is present in multiple layers of sheep small intestine, SBD-2, SBD-3, SNP-1, and LL-37 are upregulated in the specific layers of the small intestine.

HLA-DPB1 and DPB2 are genetic loci for systemic sclerosis: a genome-wide association study in Koreans with replication in North Americans.

OBJECTIVE: To identify systemic sclerosis (SSc) susceptibility loci via a genome-wide association study. METHODS: A genome-wide association study was performed in 137 patients with SSc and 564 controls from Korea using the Affymetrix Human SNP Array 5.0. After fine-mapping studies, the results were replicated in 1,107 SSc patients and 2,747 controls from a US Caucasian population. RESULTS: The single-nucleotide polymorphisms (SNPs) (rs3128930, rs7763822, rs7764491, rs3117230, and rs3128965) of HLA-DPB1 and DPB2 on chromosome 6 formed a distinctive peak with log P values for association with SSc susceptibility (P=8.16x10(-13)). Subtyping analysis of HLA-DPB1 showed that DPB1*1301 (P=7.61x10(-8)) and DPB1*0901 (P=2.55x10(-5)) were the subtypes most susceptible to SSc in Korean subjects. In US Caucasians, 2 pairs of SNPs, rs7763822/rs7764491 and rs3117230/rs3128965, showed strong association with SSc patients who had either circulating anti-DNA topoisomerase I (P=7.58x10(-17)/4.84x10(-16)) or anticentromere autoantibodies (P=1.12x10(-3)/3.2x10(-5)), respectively. CONCLUSION: The results of our genome-wide association study in Korean subjects indicate that the region of HLA-DPB1 and DPB2 contains the loci most susceptible to SSc in a Korean population. The confirmatory studies in US Caucasians indicate that specific SNPs of HLA-DPB1 and/or DPB2 are strongly associated with US Caucasian patients with SSc who are positive for anti-DNA topoisomerase I or anticentromere autoantibodies.

Altered transmission of HOX and apoptotic SNPs identify a potential common pathway for clubfoot.

Clubfoot is a common birth defect that affects 135,000 newborns each year worldwide. It is characterized by equinus deformity of one or both feet and hypoplastic calf muscles. Despite numerous study approaches, the cause(s) remains poorly understood although a multifactorial etiology is generally accepted. We considered the HOXA and HOXD gene clusters and insulin-like growth factor binding protein 3 (IGFBP3) as candidate genes because of their important roles in limb and muscle morphogenesis. Twenty SNPs from the HOXA and HOXD gene clusters and 12 SNPs in IGFBP3 were genotyped in a sample composed of non-Hispanic white and Hispanic multiplex and simplex families (discovery samples) and a second sample of non-Hispanic white simplex trios (validation sample). Four SNPs (rs6668, rs2428431, rs3801776, and rs3779456) in the HOXA cluster demonstrated altered transmission in the discovery sample, but only rs3801776, located in the HOXA basal promoter region, showed altered transmission in both the discovery and validation samples (P = 0.004 and 0.028). Interestingly, HOXA9 is expressed in muscle during development. An SNP in IGFBP3, rs13223993, also showed altered transmission (P = 0.003) in the discovery sample. Gene-gene interactions were identified between variants in HOXA, HOXD, and IGFBP3 and with previously associated SNPs in mitochondrial-mediated apoptotic genes. The most significant interactions were found between CASP3 SNPS and variants in HOXA, HOXD, and IGFBP3. These results suggest a biologic model for clubfoot in which perturbation of HOX and apoptotic genes together affect muscle and limb development, which may cause the downstream failure of limb rotation into a plantar grade position.

The Oxytricha trifallax macronuclear genome: a complex eukaryotic genome with 16,000 tiny chromosomes.

The macronuclear genome of the ciliate Oxytricha trifallax displays an extreme and unique eukaryotic genome architecture with extensive genomic variation. During sexual genome development, the expressed, somatic macronuclear genome is whittled down to the genic portion of a small fraction (∼5%) of its precursor "silent" germline micronuclear genome by a process of "unscrambling" and fragmentation. The tiny macronuclear "nanochromosomes" typically encode single, protein-coding genes (a small portion, 10%, encode 2-8 genes), have minimal noncoding regions, and are differentially amplified to an average of ∼2,000 copies. We report the high-quality genome assembly of ∼16,000 complete nanochromosomes (∼50 Mb haploid genome size) that vary from 469 bp to 66 kb long (mean ∼3.2 kb) and encode ∼18,500 genes. Alternative DNA fragmentation processes ∼10% of the nanochromosomes into multiple isoforms that usually encode complete genes. Nucleotide diversity in the macronucleus is very high (SNP heterozygosity is ∼4.0%), suggesting that Oxytricha trifallax may have one of the largest known effective population sizes of eukaryotes. Comparison to other ciliates with nonscrambled genomes and long macronuclear chromosomes (on the order of 100 kb) suggests several candidate proteins that could be involved in genome rearrangement, including domesticated MULE and IS1595-like DDE transposases. The assembly of the highly fragmented Oxytricha macronuclear genome is the first completed genome with such an unusual architecture. This genome sequence provides tantalizing glimpses into novel molecular biology and evolution. For example, Oxytricha maintains tens of millions of telomeres per cell and has also evolved an intriguing expansion of telomere end-binding proteins. In conjunction with the micronuclear genome in progress, the O. trifallax macronuclear genome will provide an invaluable resource for investigating programmed genome rearrangements, complementing studies of rearrangements arising during evolution and disease.

Dilution of candidates: the case of iron-related genes in restless legs syndrome

Restless legs syndrome (RLS) is a common multifactorial disease. Some genetic risk factors have been identified. RLS susceptibility also has been related to iron. We therefore asked whether known iron-related genes are candidates for association with RLS and, vice versa, whether known RLS-associated loci influence iron parameters in serum. RLS/control samples (n = 954/1814 in the discovery step, 735/736 in replication 1, and 736/735 in replication 2) were tested for association with SNPs located within 4 Mb intervals surrounding each gene from a list of 111 iron-related genes using a discovery threshold of P = 5 × 10(-4). Two population cohorts (KORA F3 and F4 with together n = 3447) were tested for association of six known RLS loci with iron, ferritin, transferrin, transferrin-saturation, and soluble transferrin receptor. Results were negative. None of the candidate SNPs at the iron-related gene loci was confirmed significantly. An intronic SNP, rs2576036, of KATNAL2 at 18q21.1 was significant in the first (P = 0.00085) but not in the second replication step (joint nominal P-value = 0.044). Especially, rs1800652 (C282Y) in the HFE gene did not associate with RLS. Moreover, SNPs at the known RLS loci did not significantly affect serum iron parameters in the KORA cohorts. In conclusion, the correlation between RLS and iron parameters in serum may be weaker than assumed. Moreover, in a general power analysis, we show that genetic effects are diluted if they are transmitted via an intermediate trait to an end-phenotype. Sample size formulas are provided for small effect sizes.