Interindividual Variability In Expression Research Articles

Cytochrome P450 Transcriptional Regulation by Testis-Specific Y-Encoded-Like Protein: Identification of Novel Upstream Transcription Factors.

Cytochrome P450s (CYPs) display significant inter-individual variation in expression, much of which remains unexplained by known CYP single-nucleotide polymorphisms (SNPs). Testis-specific Y-encoded-like proteins (TSPYLs) are transcriptional regulators for several drug-metabolizing CYPs including CYP3A4 However, transcription factors (TFs) that might influence CYP expression through an effect on TSPYL expression are unknown. Therefore, we studied regulators of TSPYL expression in hepatic cell lines and their possible SNP-dependent variation. Specifically, we identified candidate TFs that might influence TSPYL expression using the ENCODE ChIPseq database. Subsequently, the expression of TSPYL1/2/4 as well as that of selected CYP targets for TSPYL regulation were assayed in hepatic cell lines before and after knockdown of TFs that might influence CYP expression through TSPYL-dependent mechanisms. Those results were confirmed by studies of TF binding to TSPYL1/2/4 gene promoter regions. In hepatic cell lines, knockdown of the REST and ZBTB7A TFs resulted in decreased TSPYL1 and TSPYL4 expression and increased CYP3A4 expression, changes reversed by TSPYL1/4 overexpression. Potential binding sites for REST and ZBTB7A on the promoters of TSPYL1 and TSPYL4 were confirmed by chromatin immunoprecipitation. Finally, common SNP variants in upstream binding sites on the TSPYL1/4 promoters were identified and luciferase reporter constructs confirmed SNP-dependent modulation of TSPYL1/4 gene transcription. In summary, we identified REST and ZBTB7A as regulators of the expression of TSPYL genes which themselves can contribute to regulation of CYP expression and-potentially-of drug metabolism. SNP-dependent modulation of TSPYL transcription may contribute to individual variation in both CYP expression and-downstream-drug response phenotypes. SIGNIFICANCE STATEMENT: Testis-specific Y-encoded-like proteins (TSPYLs) are transcriptional regulators of cytochrome P450 (CYP) gene expression. Here, we report that variation in TSPYL expression as a result of the effects of genetically regulated TSPYL transcription factors is an additional factor that could result in downstream variation in CYP expression and potentially, as a result, variation in drug biotransformation.

Analysis of polymorphisms of canine Cytochrome P 450-CYP2D15.

Cytochrome P450 (CYP) proteins constitute a large ancient family of oxidative enzymes essential for the efficient elimination of a wide variety of clinically used drugs. Polymorphic variants of human CYP2D6 are associated with the conversion rate and efficacy of several drugs such as antidepressants. Polymorphisms of the canine orthologue CYP2D15 are of interest because these antidepressants are also used in dogs with behavioral problems and the outcome of the treatment is variable. However, the annotated CYP2D15 gene is incomplete and inaccurately assembled in CanFam3.1, hampering DNA sequence analysis of the gene in individual dogs. We elucidated the complete exon–intron structure of CYP2D15 to enable comprehensive genotyping of the gene using genomic DNA. We surveyed variations of the gene in four diverse dog breeds and identified novel polymorphisms in exon 2 in border collies. Further investigation to establish the impact of these canine CYP2D15 polymorphisms on interindividual variability in expression and function of this metabolizing enzyme is now feasible. Further knowledge of CYP pharmacogenetics will help individualize therapy and thereby increase therapeutic efficacy, especially in the use of antidepressants in veterinary behavioral medicine.

Gene expression variation and parental allele inheritance in a Xiphophorus interspecies hybridization model.

Understanding the genetic mechanisms underlying segregation of phenotypic variation through successive generations is important for understanding physiological changes and disease risk. Tracing the etiology of variation in gene expression enables identification of genetic interactions, and may uncover molecular mechanisms leading to the phenotypic expression of a trait, especially when utilizing model organisms that have well-defined genetic lineages. There are a plethora of studies that describe relationships between gene expression and genotype, however, the idea that global variations in gene expression are also controlled by genotype remains novel. Despite the identification of loci that control gene expression variation, the global understanding of how genome constitution affects trait variability is unknown. To study this question, we utilized Xiphophorus fish of different, but tractable genetic backgrounds (inbred, F1 interspecies hybrids, and backcross hybrid progeny), and measured each individual’s gene expression concurrent with the degrees of inter-individual expression variation. We found, (a) F1 interspecies hybrids exhibited less variability than inbred animals, indicting gene expression variation is not affected by the fraction of heterozygous loci within an individual genome, and (b), that mixing genotypes in backcross populations led to higher levels of gene expression variability, supporting the idea that expression variability is caused by heterogeneity of genotypes of cis or trans loci. In conclusion, heterogeneity of genotype, introduced by inheritance of different alleles, accounts for the largest effects on global phenotypical variability.

Environmental Canalization of Life Span and Gene Expression in Caenorhabditis elegans.

Animals, particularly poikilotherms, exhibit distinct physiologies at different environmental temperatures. Here, we hypothesized that temperature-based differences in physiology could affect the amount of variation in complex quantitative traits. Specifically, we examined, in Caenorhabditis elegans, how different temperatures (15°C, 20°C, and 25°C) affected the amount of interindividual variation in life span and also expression of three reporter genes-transcriptional reporters for vit-2, gpd-2, and hsp-16.2 (a life-span biomarker). We found the expected inverse relationship between temperature and average life span. Surprisingly, we found that at the highest temperature, there were fewer differences between individuals in life span and less interindividual variation in expression of all three reporters. We suggest that growth at 25°C might canalize (reduce interindividual differences in) life span and expression of some genes by eliciting a small constitutive heat shock response. Growth at 25°C requires wild-type hsf-1, which encodes the main heat shock response transcriptional activator. We speculate that increased chaperone activity at 25°C may reduce interindividual variation in gene expression by increasing protein folding efficiency. We hypothesize that reduced variation in gene expression may ultimately cause reduced variation in life span.

Caenorhabditis elegans Genes Affecting Interindividual Variation in Life-span Biomarker Gene Expression.

Genetically identical organisms grown in homogenous environments differ in quantitative phenotypes. Differences in one such trait, expression of a single biomarker gene, can identify isogenic cells or organisms that later manifest different fates. For example, in isogenic populations of young adult Caenorhabditis elegans, differences in Green Fluorescent Protein (GFP) expressed from the hsp-16.2 promoter predict differences in life span. Thus, it is of interest to determine how interindividual differences in biomarker gene expression arise. Prior reports showed that the thermosensory neurons and insulin signaling systems controlled the magnitude of the heat shock response, including absolute expression of hsp-16.2. Here, we tested whether these regulatory signals might also influence variation in hsp-16.2 reporter expression. Genetic experiments showed that the action of AFD thermosensory neurons increases interindividual variation in biomarker expression. Further genetic experimentation showed the insulin signaling system acts to decrease interindividual variation in life-span biomarker expression; in other words, insulin signaling canalizes expression of the hsp-16.2-driven life-span biomarker. Our results show that specific signaling systems regulate not only expression level, but also the amount of interindividual expression variation for a life-span biomarker gene. They raise the possibility that manipulation of these systems might offer means to reduce heterogeneity in the aging process.

Long non-coding RNAs display higher natural expression variation than protein-coding genes in healthy humans.

BackgroundLong non-coding RNAs (lncRNAs) are increasingly implicated as gene regulators and may ultimately be more numerous than protein-coding genes in the human genome. Despite large numbers of reported lncRNAs, reference annotations are likely incomplete due to their lower and tighter tissue-specific expression compared to mRNAs. An unexplored factor potentially confounding lncRNA identification is inter-individual expression variability. Here, we characterize lncRNA natural expression variability in human primary granulocytes.ResultsWe annotate granulocyte lncRNAs and mRNAs in RNA-seq data from 10 healthy individuals, identifying multiple lncRNAs absent from reference annotations, and use this to investigate three known features (higher tissue-specificity, lower expression, and reduced splicing efficiency) of lncRNAs relative to mRNAs. Expression variability was examined in seven individuals sampled three times at 1- or more than 1-month intervals. We show that lncRNAs display significantly more inter-individual expression variability compared to mRNAs. We confirm this finding in two independent human datasets by analyzing multiple tissues from the GTEx project and lymphoblastoid cell lines from the GEUVADIS project. Using the latter dataset we also show that including more human donors into the transcriptome annotation pipeline allows identification of an increasing number of lncRNAs, but minimally affects mRNA gene number.ConclusionsA comprehensive annotation of lncRNAs is known to require an approach that is sensitive to low and tight tissue-specific expression. Here we show that increased inter-individual expression variability is an additional general lncRNA feature to consider when creating a comprehensive annotation of human lncRNAs or proposing their use as prognostic or disease markers.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-016-0873-8) contains supplementary material, which is available to authorized users.

Changes in correlation between promoter methylation and gene expression in cancer.

BackgroundMethylation of high-density CpG regions known as CpG Islands (CGIs) has been widely described as a mechanism associated with gene expression regulation. Aberrant promoter methylation is considered a hallmark of cancer involved in silencing of tumor suppressor genes and activation of oncogenes. However, recent studies have also challenged the simple model of gene expression control by promoter methylation in cancer, and the precise mechanism of and role played by changes in DNA methylation in carcinogenesis remains elusive.ResultsUsing a large dataset of 672 matched cancerous and healthy methylomes, gene expression, and copy number profiles accross 3 types of tissues from The Cancer Genome Atlas (TCGA), we perform a detailed meta-analysis to clarify the interplay between promoter methylation and gene expression in normal and cancer samples. On the one hand, we recover the existence of a CpG island methylator phenotype (CIMP) with prognostic value in a subset of breast, colon and lung cancer samples, where a common subset of promoter CGIs hypomethylated in normal samples become hypermethylated. However, this hypermethylation is not accompanied by a decrease in expression of the corresponding genes, which are already lowly expressed in the normal genes. On the other hand, we identify tissue-specific sets of genes, different between normal and cancer samples, whose inter-individual variation in expression is significantly correlated with the variation in methylation of the 3’ flanking regions of the promoter CGIs. These subsets of genes are not the same in the different tissues, nor between normal and cancerous samples, but transcription factors are over-represented in all subsets.ConclusionOur results suggest that epigenetic reprogramming in cancer does not contribute to cancer development via direct inhibition of gene expression through promoter hypermethylation. It may instead modify how the expression of a few specific genes, particularly transcription factors, are associated with DNA methylation variations in a tissue-dependent manner.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1994-2) contains supplementary material, which is available to authorized users.

Aberrant RNA splicing in cancer; expression changes and driver mutations of splicing factor genes.

Alternative splicing is a widespread process contributing to structural transcript variation and proteome diversity. In cancer, the splicing process is commonly disrupted, resulting in both functional and non-functional end-products. Cancer-specific splicing events are known to contribute to disease progression; however, the dysregulated splicing patterns found on a genome-wide scale have until recently been less well-studied. In this review, we provide an overview of aberrant RNA splicing and its regulation in cancer. We then focus on the executors of the splicing process. Based on a comprehensive catalog of splicing factor encoding genes and analyses of available gene expression and somatic mutation data, we identify cancer-associated patterns of dysregulation. Splicing factor genes are shown to be significantly differentially expressed between cancer and corresponding normal samples, and to have reduced inter-individual expression variation in cancer. Furthermore, we identify enrichment of predicted cancer-critical genes among the splicing factors. In addition to previously described oncogenic splicing factor genes, we propose 24 novel cancer-critical splicing factors predicted from somatic mutations.

Interindividual variation in gene expression responses and metabolite formation in acetaminophen-exposed primary human hepatocytes.

Acetaminophen (APAP) is a readily available over-the-counter drug and is one of the most commonly used analgesics/antipyretics worldwide. Large interindividual variation in susceptibility toward APAP-induced liver failure has been reported. However, the exact underlying factors causing this variability in susceptibility are still largely unknown. The aim of this study was to better understand this variability in response to APAP by evaluating interindividual differences in gene expression changes and APAP metabolite formation in primary human hepatocytes (PHH) from several donors (n = 5) exposed in vitro to a non-toxic to toxic APAP dose range. To evaluate interindividual variation, gene expression data/levels of metabolites were plotted against APAP dose/donor. The correlation in APAP dose response between donors was calculated by comparing data points from one donor to the data points of all other donors using a Pearson-based correlation analysis. From that, a correlation score/donor for each gene/metabolite was defined, representing the similarity of the omics response to APAP in PHH of a particular donor to all other donors. The top 1 % highest variable genes were selected for further evaluation using gene set overrepresentation analysis. The biological processes in which the genes with high interindividual variation in expression were involved include liver regeneration, inflammatory responses, mitochondrial stress responses, hepatocarcinogenesis, cell cycle, and drug efficacy. Additionally, the interindividual variation in the expression of these genes could be associated with the variability in expression levels of hydroxyl/methoxy-APAP and C8H13O5N-APAP-glucuronide. The before-mentioned metabolites or their derivatives have also been reported in blood of humans exposed to therapeutic APAP doses. Possibly these findings can contribute to elucidating the causative factors of interindividual susceptibility toward APAP. Electronic supplementary materialThe online version of this article (doi:10.1007/s00204-015-1545-2) contains supplementary material, which is available to authorized users.

Allele Frequency Distributions of the Drug Metabolizer Genes <i>CYP2C9*2</i>, <i>CYP2C9*3</i>, and <i>CYP2C19*17</i> in the Buginese Population of Indonesia

The present study is part of the genetic mapping of Indonesia focusing on drug metabolizing enzymes, which started with the Buginese population of Makassar, South Sulawesi. The two CYP450 gene subfamilies, i.e. CYP2C9 and CYP2C19 are of interest as they exhibit wide inter-individual variation in expression, which influence the drug metabolism capacity. The CYP2C9 alleles of interest in this study were CYP2C9*2 and *3, and of CYP2C19 was CYP2C19*17. The study aimed to determine the frequencies of the CYP2C9 genotype, which contains *1, *2 and *3 alleles, and the CYP2C19 genotype, which comprises the *1 and *17 alleles in the Buginese. Ninety six Buginese subjects, comprising 48 males and 48 females were studied. CYP2C9 and CYP2C19 alleles were detected by a PCR-RFLP assay method. Results showed that there was no CYP2C9*2 allele present, while the frequencies of CYP2C9*3 and CYP2C19*17 overall were 1.56 % and 4.68 %, respectively. The frequency of the CYP2C9*3 allele in females was 2.08%, and not statistically different from that in males (1.05%). The frequency of the CYP2C19*17 allele in females (8.33%), was significantly different (P<0.05) from that in males (1.05%). No subject carried the CYP2C9*2/*2, CYP2C9*3/*3, CYP2C19*17/*17, or CYP2C9*3/CYP2C19*17 genotype. The study is the first to describe the drug metabolizing enzyme polymorphisms, CYP2C9 and CYP2C19, in the Indonesian Buginese population.

Ultrasensitive proteome analysis using paramagnetic bead technology.

In order to obtain a systems-level understanding of a complex biological system, detailed proteome information is essential. Despite great progress in proteomics technologies, thorough interrogation of the proteome from quantity-limited biological samples is hampered by inefficiencies during processing. To address these challenges, here we introduce a novel protocol using paramagnetic beads, termed Single-Pot Solid-Phase-enhanced Sample Preparation (SP3). SP3 provides a rapid and unbiased means of proteomic sample preparation in a single tube that facilitates ultrasensitive analysis by outperforming existing protocols in terms of efficiency, scalability, speed, throughput, and flexibility. To illustrate these benefits, characterization of 1,000 HeLa cells and single Drosophila embryos is used to establish that SP3 provides an enhanced platform for profiling proteomes derived from sub-microgram amounts of material. These data present a first view of developmental stage-specific proteome dynamics in Drosophila at a single-embryo resolution, permitting characterization of inter-individual expression variation. Together, the findings of this work position SP3 as a superior protocol that facilitates exciting new directions in multiple areas of proteomics ranging from developmental biology to clinical applications.

Developmental trajectory of intestinal MDR1/ABCB1 mRNA expression in children

Developmental trajectory of intestinal <scp>MDR1/ABCB1</scp> mRNA expression in children

How to get the most from microarray data: advice from reverse genomics

BackgroundWhole-genome profiling of gene expression is a powerful tool for identifying cancer-associated genes. Genes differentially expressed between normal and tumorous tissues are usually considered to be cancer associated. We recently demonstrated that the analysis of interindividual variation in gene expression can be useful for identifying cancer associated genes. The goal of this study was to identify the best microarray data–derived predictor of known cancer associated genes.ResultsWe found that the traditional approach of identifying cancer genes—identifying differentially expressed genes—is not very efficient. The analysis of interindividual variation of gene expression in tumor samples identifies cancer-associated genes more effectively. The results were consistent across 4 major types of cancer: breast, colorectal, lung, and prostate. We used recently reported cancer-associated genes (2011–2012) for validation and found that novel cancer-associated genes can be best identified by elevated variance of the gene expression in tumor samples.ConclusionsThe observation that the high interindividual variation of gene expression in tumor tissues is the best predictor of cancer-associated genes is likely a result of tumor heterogeneity on gene level. Computer simulation demonstrates that in the case of heterogeneity, an assessment of variance in tumors provides a better identification of cancer genes than does the comparison of the expression in normal and tumor tissues. Our results thus challenge the current paradigm that comparing the mean expression between normal and tumorous tissues is the best approach to identifying cancer-associated genes; we found that the high interindividual variation in expression is a better approach, and that using variation would improve our chances of identifying cancer-associated genes.

Genome-Wide Variation of Cytosine Modifications Between European and African Populations and the Implications for Complex Traits

Elucidating cytosine modification differences between human populations can enhance our understanding of ethnic specificity in complex traits. In this study, cytosine modification levels in 133 HapMap lymphoblastoid cell lines derived from individuals of European or African ancestry were profiled using the Illumina HumanMethylation450 BeadChip. Approximately 13% of the analyzed CpG sites showed differential modification between the two populations at a false discovery rate of 1%. The CpG sites with greater modification levels in European descent were enriched in the proximal regulatory regions, while those greater in African descent were biased toward gene bodies. More than half of the detected population-specific cytosine modifications could be explained primarily by local genetic variation. In addition, a substantial proportion of local modification quantitative trait loci exhibited population-specific effects, suggesting that genetic epistasis and/or genotype × environment interactions could be common. Distinct correlations were observed between gene expression levels and cytosine modifications in proximal regions and gene bodies, suggesting epigenetic regulation of interindividual expression variation. Furthermore, quantitative trait loci associated with population-specific modifications can be colocalized with expression quantitative trait loci and single nucleotide polymorphisms previously identified for complex traits with known racial disparities. Our findings revealed abundant population-specific cytosine modifications and the underlying genetic basis, as well as the relatively independent contribution of genetic and epigenetic variations to population differences in gene expression.

Genome-wide survey of interindividual differences of RNA stability in human lymphoblastoid cell lines

The extent to which RNA stability differs between individuals and its contribution to the interindividual expression variation remain unknown. We conducted a genome-wide analysis of RNA stability in seven human HapMap lymphoblastoid cell lines (LCLs) and analyzed the effect of DNA sequence variation on RNA half-life differences. Twenty-six percent of the expressed genes exhibited RNA half-life differences between LCLs at a false discovery rate (FDR) < 0.05, which accounted for ~ 37% of the gene expression differences between individuals. Nonsense polymorphisms were associated with reduced RNA half-lives. In genes presenting interindividual RNA half-life differences, higher coding GC3 contents (G and C percentages at the third-codon positions) were correlated with increased RNA half-life. Consistently, G and C alleles of single nucleotide polymorphisms (SNPs) in protein coding sequences were associated with enhanced RNA stability. These results suggest widespread interindividual differences in RNA stability related to DNA sequence and composition variation.

Evidence for widespread changes in promoter methylation profile in human placenta in response to increasing gestational age and environmental/stochastic factors

BackgroundThe human placenta facilitates the exchange of nutrients, gas and waste between the fetal and maternal circulations. It also protects the fetus from the maternal immune response. Due to its role at the feto-maternal interface, the placenta is subject to many environmental exposures that can potentially alter its epigenetic profile. Previous studies have reported gene expression differences in placenta over gestation, as well as inter-individual variation in expression of some genes. However, the factors contributing to this variation in gene expression remain poorly understood.ResultsIn this study, we performed a genome-wide DNA methylation analysis of gene promoters in placenta tissue from three pregnancy trimesters. We identified large-scale differences in DNA methylation levels between first, second and third trimesters, with an overall progressive increase in average methylation from first to third trimester. The most differentially methylated genes included many immune regulators, reflecting the change in placental immuno-modulation as pregnancy progresses. We also detected increased inter-individual variation in the third trimester relative to first and second, supporting an accumulation of environmentally induced (or stochastic) changes in DNA methylation pattern. These highly variable genes were enriched for those involved in amino acid and other metabolic pathways, potentially reflecting the adaptation of the human placenta to different environments.ConclusionsThe identification of cellular pathways subject to drift in response to environmental influences provide a basis for future studies examining the role of specific environmental factors on DNA methylation pattern and placenta-associated adverse pregnancy outcomes.

Canalization effect in the coagulation cascade and the interindividual variability of oral anticoagulant response. a simulation Study

BackgroundIncreasing the predictability and reducing the rate of side effects of oral anticoagulant treatment (OAT) requires further clarification of the cause of about 50% of the interindividual variability of OAT response that is currently unaccounted for. We explore numerically the hypothesis that the effect of the interindividual expression variability of coagulation proteins, which does not usually result in a variability of the coagulation times in untreated subjects, is unmasked by OAT.ResultsWe developed a stochastic variant of the Hockin-Mann model of the tissue factor coagulation pathway, using literature data for the variability of coagulation protein levels in the blood of normal subjects. We simulated in vitro coagulation and estimated the Prothrombin Time and the INR across a model population. In a model of untreated subjects a "canalization effect" can be observed in that a coefficient of variation of up to 33% of each protein level results in a simulated INR of 1 with a clinically irrelevant dispersion of 0.12. When the mean and the standard deviation of vitamin-K dependent protein levels were reduced by 80%, corresponding to the usual Warfarin treatment intensity, the simulated INR was 2.98 ± 0.48, a clinically relevant dispersion, corresponding to a reduction of the canalization effect.Then we combined the Hockin-Mann stochastic model with our previously published model of population response to Warfarin, that takes into account the genetical and the phenotypical variability of Warfarin pharmacokinetics and pharmacodynamics. We used the combined model to evaluate the coagulation protein variability effect on the variability of the Warfarin dose required to reach an INR target of 2.5. The dose variance when removing the coagulation protein variability was 30% lower. The dose was mostly related to the pretreatment levels of factors VII, X, and the tissue factor pathway inhibitor (TFPI).ConclusionsIt may be worth exploring in experimental studies whether the pretreatment levels of coagulation proteins, in particular VII, X and TFPI, are predictors of the individual warfarin dose, even though, maybe due to a canalization-type effect, their effect on the INR variance in untreated subjects appears low.

The endothelial microenvironment in the venous valvular sinus: thromboresistance trends and inter-individual variation

The valve sinuses of the deep venous system are frequent sites of venous thrombus initiation. We previously reported that, in comparison with the non-valvular lumenal endothelium, the valve sinus endothelium had decreased expression of von Willebrand factor (vWF) and increased expression of endothelial protein C receptor (EPCR) and thrombomodulin (TM), suggesting alteration in the procoagulant/anticoagulant balance. We hypothesized that increased stasis in the deeper recesses of the venous valves would be associated with a gradient of increased thromboresistance. Expression of EPCR, TM, and vWF was analyzed via quantitative confocal immunofluorescence in residual saphenous veins collected following coronary artery bypass procedures. In agreement with our hypothesis, endothelial expression of vWF in the valve sinus decreased from the uppermost to the deepest region of the valve sinus. In contrast to our hypothesis, EPCR expression decreased from the uppermost to the deepest region of the valve sinus (p<0.001) and TM expression remained unchanged throughout the valve sinus. Comparison of the non-valvular lumenal endothelium with the valve sinus endothelium demonstrated significantly decreased vWF expression (p<0.001) in the valvular sinus consistent with our previous report; however, we did not observe statistically significant differences in EPCR or TM expression in this comparison. In addition, remarkable inter-individual variation in expression of these three proteins was also observed. These findings suggest that the genesis of these observations is more complex than predicted by our initial hypothesis, likely due, at least in part, to the complex rheology of the valvular sinus microenvironment.

Long-term balancing selection maintains trans-specific polymorphisms in the human TRIM5 gene

The human TRIM5 genes encodes a retroviral restriction factor (TRIM5α). Evolutionary analyses of this gene in mammals have revealed a complex and multifaceted scenario, suggesting that TRIM5 has been the target of exceptionally strong selective pressures, possibly exerted by recurrent waves of retroviral infections. TRIM5 displays inter-individual expression variability in humans and high levels of TRIM5 mRNA have been associated with a reduced risk of HIV-1 infection. We resequenced TRIM5 in chimpanzees and identified two polymorphisms in intron 1 that are shared with humans. Analysis of the gene region encompassing the two trans-specific variants in human populations identified exceptional nucleotide diversity levels and an excess of polymorphism compared to fixed divergence. Most tests rejected the null hypothesis of neutral evolution for this region and haplotype analysis revealed the presence of two deeply separated clades. Calculation of the time to the most recent common ancestor (TMRCA) for TRIM5 haplotypes yielded estimates ranging between 4 and 7 million years. Overall, these data indicate that long-term balancing selection, an extremely rare process outside MHC genes, has maintained trans-specific polymorphisms in the first intron of TRIM5. Bioinformatic analyses indicated that variants in intron 1 may affect transcription factor-binding sites and, therefore, TRIM5 transcriptional activity. Data herein confirm an extremely complex evolutionary history of TRIM5 genes in primates and open the possibility that regulatory variants in the gene modulate the susceptibility to HIV-1.

Functional analysis of genetic variations in the 5′-flanking region of the human MDR1 gene

P-glycoprotein (P-gp), the product of the MDR1 gene, shows large interindividual variations in expression, which leads to differences in the pharmacokinetics of the substrate drugs. The functions of single nucleotide polymorphisms located in the nuclear receptor-responsive element of the 5′-flanking region in the human MDR1 gene were analyzed in order to clarify the mechanism underlying the interindividual variation in P-gp expression. Electrophoretic mobility shift assays revealed that the − 7833C > T substitution in the nuclear receptor-responsive region of MDR1 decreases the binding affinities of four nuclear receptors to their responsive elements: vitamin D receptor (VDR), thyroid hormone receptor (TR), constitutive androstane receptor (CAR), and pregnane X receptor (PXR). A reporter gene assay revealed that the C-to-T substitution at − 7833 also reduces the transcriptional activation of MDR1 by VDR, TRβ, CAR, and PXR. However, another SNP (− 1211T > C substitution), which results in the formation of a xenobiotic responsive element-like sequence and a hypoxia responsive element-like sequence, failed to affect the aryl hydrocarbon receptor-dependent and hypoxia-induced transcriptional activation of MDR1. Although the frequency of the − 7833C > T substitution in MDR1 is relatively low, the SNP is crucial because it may alter the pharmacokinetics of P-gp substrates in a small subset of the population.