A signature of balancing selection in the region upstream to the human UGT2B4 gene and implications for breast cancer risk

UDP-glucuronosyltransferase 2 family, polypeptide B4 (UGT2B4) is an important metabolizing enzyme involved in the clearance of many xenobiotics and endogenous substrates, especially steroid hormones and bile acids. The HapMap data show that numerous SNPs upstream of UGT2B4 are in near-perfect linkage disequilibrium with each other and occur at intermediate frequency, indicating that this region might contain a target of natural selection. To investigate this possibility, we chose three regions (4.8 kb in total) for resequencing and observed a striking excess of intermediate-frequency alleles that define two major haplotypes separated by many mutation events and with little differentiation across populations, thus suggesting that the variation pattern upstream UGT2B4 is highly unusual and may be the result of balancing selection. We propose that this pattern is due to the maintenance of a regulatory polymorphism involved in the fine tuning of UGT2B4 expression so that heterozygous genotypes result in optimal enzyme levels. Considering the important role of steroid hormones in breast cancer susceptibility, we hypothesized that variation in this region could predispose to breast cancer. To test this hypothesis, we genotyped tag SNP rs13129471 in 1,261 patients and 825 normal women of African ancestry from three populations. The frequency comparison indicated that rs13129471 was significantly associated with breast cancer after adjusting for ethnicity [P = 0.003; heterozygous odds ratio (OR) 1.02, 95% confidence interval (CI) 0.81-1.28; homozygous OR 1.50, 95% CI 1.15-1.95]. Our results provide new insights into UGT2B4 sequence variation and indicate that a signal of natural selection may lead to the identification of disease susceptibility variants.

Differential allelic expression of c.1568C > A at UGT2B15 is due to variation in a novel cis-regulatory element in the 3'UTR

Differential allelic expression (DAE) is a powerful tool to identify cis-regulatory elements for gene expression. The UDP-glucuronosyltransferase 2 family, polypeptide B15 (UGT2B15), is an important enzyme involved in the metabolism of multiple endobiotics and xenobiotics. In the present study, we measured the relative expression of two alleles at SNP c.1568C>A (rs4148269) in this gene, which causes an amino acid substitution (T523K). An excess of the C over the A allele was consistently observed in both liver (P=0.0021) and breast (P=0.012) samples, suggesting that SNP(s) in strong linkage disequilibrium (LD) with c.1568C>A can regulate UGT2B15 expression in both tissues. By resequencing, one such SNP, c.1761T>C (rs3100) in 3' untranslated region (UTR), was identified. Reporter gene assays showed that the 1761T allele results in a significantly higher gene expression level than the 1761C allele in HepG2, MCF-7, LNCaP, and Caco-2 cell lines (all P<0.001), thus indicating that this variation can regulate UGT2B15 gene expression in liver, breast, colon, and prostate tissues. Considering its location, we postulated that this SNP is within an unknown microRNA binding site and can influence microRNA targeting. Considering the importance of UGT2B15 in metabolism, we proposed that this SNP might contribute to multiple cancer risk and variability in drug response.

SNP discovery, expression and cis-regulatory variation in the UGT2B genes

UGT2B enzymes metabolize multiple endogenous and exogenous molecules, including steroid hormones and clinical drugs. However, little is known about the inter-individual variation in gene expression and its determinants. We re-sequenced candidate regulatory regions and the partial coding regions (41.1 kb) of UGT2B genes and identified 332 genetic variants. We measured gene expression in normal breast and liver samples and observed different patterns. The expression levels varied greatly across individuals in both tissues and were significantly correlated with each other in liver. Genotyping of tagging single-nucleotide polymorphisms (SNPs) in the same samples and association tests between genotype and transcript levels identified 62 variants that were associated with at least one UGT2B mRNA levels in either tissue. Most of these cis-regulatory SNPs were not shared between tissues, suggesting that this gene family is regulated in a tissue-specific manner. Our results provide insight into studying the role of UGT2B variation in hormone-dependent cancers and drug response.

Allele-specific down-regulation of RPTOR expression induced by retinoids contributes to climate adaptations

The mechanistic target of rapamycin (MTOR) pathway regulates cell growth, energy homeostasis, apoptosis, and immune response. The regulatory associated protein of MTOR encoded by the RPTOR gene is a key component of this pathway. A previous survey of candidate genes found that RPTOR contains multiple SNPs with strong correlations between allele frequencies and climate variables, consistent with the action of selective pressures that vary across environments. Using data from a recent genome scan for selection signals, we honed in on a SNP (rs11868112) 26 kb upstream to the transcription start site of RPTOR that exhibits the strongest association with temperature variables. Transcription factor motif scanning and mining of recently mapped transcription factor binding sites identified a binding site for POU class 2 homeobox 1 (POU2F1) spanning the SNP and an adjacent retinoid acid receptor (RAR) binding site. Using expression quantification, chromatin immunoprecipitation (ChIP), and reporter gene assays, we demonstrate that POU2F1 and RARA do bind upstream of the RPTOR gene to regulate its expression in response to retinoids; this regulation is affected by the allele status at rs11868112 with the derived allele resulting in lower expression levels. We propose a model in which the derived allele influences thermogenesis or immune response by altering MTOR pathway activity and thereby increasing fitness in colder climates. Our results show that signatures of genetic adaptations can identify variants with functional effects, consistent with the idea that selection signals may be used for SNP annotation.

Climate has exerted strong selective pressures in human populations during their dispersal, and signatures of these adaptations are still detectable in the geographic distribution of polymorphisms. RPTOR is a key component of the mechanistic target of rapamycin pathway, which regulates cell growth, metabolism, and immune response; and its deregulation is associated with human diseases, including cancer and diabetes. Previous studies showed that variation in RPTOR carry strong signatures of adaptations to different climates. Here, we used evolutionary genetics approaches coupled with transcription factor motif data mining to refine the location of the selection target. We then used functional assays to show that the selected polymorphism resides in a sequence element that regulates gene expression levels in response to retinoids. The derived allele at this SNP, which results in lower expression levels, increases in frequency with decreasing temperatures, consistent with the notion that it confers a selective advantage in colder climates possibly through its effects on energy metabolism or immune response. These results suggest a novel regulatory role for retinoids in MTOR signaling. Moreover, they support the proposal that evolutionary approaches can be informative for SNP functional annotation.

Allelic imbalance (AI) identifies novel tissue-specific cis-regulatory variation for human UGT2B15

Allelic imbalance (AI) is a powerful tool to identify cis-regulatory variation for gene expression. UGT2B15 is an important enzyme involved in the metabolism of multiple endobiotics and xenobiotics. In this study, we measured the relative expression of two alleles at this gene by using SNP rs1902023:G>T. An excess of the G over the T allele was consistently observed in liver (P<0.001), but not in breast (P=0.06) samples, suggesting that SNPs in strong linkage disequilibrium with G253T regulate UGT2B15 expression in liver. Seven such SNPs were identified by resequencing the promoter and exon 1, which define two distinct haplotypes. Reporter gene assays confirmed that one haplotype displayed approximately 20% higher promoter activity compared to the other major haplotype in liver HepG2 (P<0.001), but not in breast MCF-7 (P=0.540) cells. Reporter gene assays with additional constructs pointed to rs34010522:G>T and rs35513228:C>T as the cis-regulatory variants; both SNPs were also evaluated in LNCaP and Caco-2 cells. By ChIP, we showed that the transcription factor Nrf2 binds to the region spanning rs34010522:G>T in all four cell lines. Our results provide a good example for how AI can be used to identify cis-regulatory variation and gain insights into the tissue specific regulation of gene expression.

Adaptive variation regulates the expression of the human SGK1 gene in response to stress

The Serum and Glucocorticoid-regulated Kinase1 (SGK1) gene is a target of the glucocorticoid receptor (GR) and is central to the stress response in many human tissues. Because environmental stress varies across habitats, we hypothesized that natural selection shaped the geographic distribution of genetic variants regulating the level of SGK1 expression following GR activation. By combining population genetics and molecular biology methods, we identified a variant (rs9493857) with marked allele frequency differences between populations of African and European ancestry and with a strong correlation between allele frequency and latitude in worldwide population samples. This SNP is located in a GR-binding region upstream of SGK1 that was identified using a GR ChIP-chip. SNP rs9493857 also lies within a predicted binding site for Oct1, a transcription factor known to cooperate with the GR in the transactivation of target genes. Using ChIP assays, we show that both GR and Oct1 bind to this region and that the ancestral allele at rs9493857 binds the GR-Oct1 complex more efficiently than the derived allele. Finally, using a reporter gene assay, we demonstrate that the ancestral allele is associated with increased glucocorticoid-dependent gene expression when compared to the derived allele. Our results suggest a novel paradigm in which hormonal responsiveness is modulated by sequence variation in the regulatory regions of nuclear receptor target genes. Identifying such functional variants may shed light on the mechanisms underlying inter-individual variation in response to environmental stressors and to hormonal therapy, as well as in the susceptibility to hormone-dependent diseases.

Susceptibility to many common human diseases including hypertension, heart disease, and the metabolic syndrome is associated with increased neuroendocrine signaling in response to environmental stressors. A key component of the human stress response involves increased systemic glucocorticoid secretion that in turn leads to glucocorticoid receptor (GR) activation. As a result, a variety of GR-expressing cell types undergo gene expression changes, thereby providing an integrated physiological response to stress. The SGK1 gene is a well-established GR target that promotes cellular homeostasis in response to stress. Here, we use a combination of population genetics and molecular biology approaches to identify an SNP (rs9493857) in a distant SGK1 GR-binding region with unusually large differences in allele frequency between populations of European and African ancestry. Furthermore, rs9493857 shows a strong correlation between allele frequency and distance from the equator, a pattern consistent with a varying selective advantage across environments. Indeed, the ancestral allele at rs9493857 results in increased GR-binding and glucocorticoid-regulated gene expression, suggesting that an increased stress response (i.e., glucocorticoid responsiveness) was advantageous in ancestral human populations. We speculate that, in modern times, such variation could favor the negative effects of a heightened glucocorticoid response, potentially predisposing individuals to chronic diseases such as metabolic syndrome and hypertension.

Characterization of a novel splicing variant in the RAPTOR gene

The mammalian target of rapamycin (mTOR) plays an essential role in the regulation of cell growth, proliferation and apoptosis. Raptor, the regulatory associated protein of mTOR, is an important member in this signaling pathway. In the present report,we identified and characterized a novel splicing variant of this gene, RAPTOR v2, in which exons 14-17, 474 bp in total, are omitted from the mRNA. This deletion does not change the open reading frame, but causes a nearly complete absence of HEAT repeats, which were shown to be involved in the binding of mTOR substrates. Real time PCR performed on 48 different human tissues demonstrated the ubiquitous presence of this splice variant. Quantification of mRNA levels in lymphoblastoid cell lines (LCL) from 56 unrelated HapMap individuals revealed that the expression of this splicing form is quite variable. One synonymous SNP, rs2289759 in exon 14, was predicted by ESEfinder to cause a significant gain/loss of SRp55 and/or SF2/ASF binding sites, and thus potentially influence splicing. This prediction was confirmed by linear regression analysis between the ratio of RAPTOR v2 to total RAPTOR mRNA levels and the SNP genotype in the above 56 individuals (r=0.281 and P=0.036). Moreover, the functional evaluation indicated that this splicing isoform is expected to retain the ability to bind mTOR, but is unlikely to bind mTOR substrates, hence affecting signal transduction and further cell proliferation.