A complex intronic enhancer regulates expression of the CFTR gene by direct interaction with the promoter.

Genes can maintain spatiotemporal expression patterns by long-range interactions between cis-acting elements. The cystic fibrosis transmembrane conductance regulator gene (CFTR) is expressed primarily in epithelial cells. An element located within a DNase I-hypersensitive site (DHS) 10 kb into the first intron was previously shown to augment CFTR promoter activity in a tissue-specific manner. Here, we reveal the mechanism by which this element influences CFTR transcription. We employed a high-resolution method of mapping DHS using tiled microarrays to accurately locate the intron 1 DHS. Transfection of promoter-reporter constructs demonstrated that the element displays classical tissue-specific enhancer properties and can independently recruit factors necessary for transcription initiation. In vitro DNase I footprinting analysis identified a protected region that corresponds to a conserved, predicted binding site for hepatocyte nuclear factor 1 (HNF1). We demonstrate by electromobility shift assays (EMSA) and chromatin immunoprecipitation (ChIP) that HNF1 binds to this element both in vitro and in vivo. Moreover, using chromosome conformation capture (3C) analysis, we show that this element interacts with the CFTR promoter in CFTR-expressing cells. These data provide the first insight into the three- dimensional (3D) structure of the CFTR locus and confirm the contribution of intronic cis-acting elements to the regulation of CFTR gene expression.

Co-regulation of nuclear respiratory factor-1 by NFkappaB and CREB links LPS-induced inflammation to mitochondrial biogenesis.

The nuclear respiratory factor-1 (NRF1) gene is activated by lipopolysaccharide (LPS), which might reflect TLR4-mediated mitigation of cellular inflammatory damage via initiation of mitochondrial biogenesis. To test this hypothesis, we examined NRF1 promoter regulation by NFκB, and identified interspecies-conserved κB-responsive promoter and intronic elements in the NRF1 locus. In mice, activation of Nrf1 and its downstream target, Tfam, by Escherichia coli was contingent on NFκB, and in LPS-treated hepatocytes, NFκB served as an NRF1 enhancer element in conjunction with NFκB promoter binding. Unexpectedly, optimal NRF1 promoter activity after LPS also required binding by the energy-state-dependent transcription factor CREB. EMSA and ChIP assays confirmed p65 and CREB binding to the NRF1 promoter and p65 binding to intron 1. Functionality for both transcription factors was validated by gene-knockdown studies. LPS regulation of NRF1 led to mtDNA-encoded gene expression and expansion of mtDNA copy number. In cells expressing plasmid constructs containing the NRF-1 promoter and GFP, LPS-dependent reporter activity was abolished by cis-acting κB-element mutations, and nuclear accumulation of NFκB and CREB demonstrated dependence on mitochondrial H(2)O(2). These findings indicate that TLR4-dependent NFκB and CREB activation co-regulate the NRF1 promoter with NFκB intronic enhancement and redox-regulated nuclear translocation, leading to downstream target-gene expression, and identify NRF-1 as an early-phase component of the host antibacterial defenses.

Genome-wide mRNA expression correlates of viral control in CD4+ T-cells from HIV-1-infected individuals.

There is great interindividual variability in HIV-1 viral setpoint after seroconversion, some of which is known to be due to genetic differences among infected individuals. Here, our focus is on determining, genome-wide, the contribution of variable gene expression to viral control, and to relate it to genomic DNA polymorphism. RNA was extracted from purified CD4+ T-cells from 137 HIV-1 seroconverters, 16 elite controllers, and 3 healthy blood donors. Expression levels of more than 48,000 mRNA transcripts were assessed by the Human-6 v3 Expression BeadChips (Illumina). Genome-wide SNP data was generated from genomic DNA using the HumanHap550 Genotyping BeadChip (Illumina). We observed two distinct profiles with 260 genes differentially expressed depending on HIV-1 viral load. There was significant upregulation of expression of interferon stimulated genes with increasing viral load, including genes of the intrinsic antiretroviral defense. Upon successful antiretroviral treatment, the transcriptome profile of previously viremic individuals reverted to a pattern comparable to that of elite controllers and of uninfected individuals. Genome-wide evaluation of cis-acting SNPs identified genetic variants modulating expression of 190 genes. Those were compared to the genes whose expression was found associated with viral load: expression of one interferon stimulated gene, OAS1, was found to be regulated by a SNP (rs3177979, p = 4.9E-12); however, we could not detect an independent association of the SNP with viral setpoint. Thus, this study represents an attempt to integrate genome-wide SNP signals with genome-wide expression profiles in the search for biological correlates of HIV-1 control. It underscores the paradox of the association between increasing levels of viral load and greater expression of antiviral defense pathways. It also shows that elite controllers do not have a fully distinctive mRNA expression pattern in CD4+ T cells. Overall, changes in global RNA expression reflect responses to viral replication rather than a mechanism that might explain viral control.

Novel associations of UDP-glucuronosyltransferase 2B gene variants with prostate cancer risk in a multiethnic study.

BACKGROUND: We have previously shown that a functional polymorphism of the UGT2B15 gene (rs1902023) was associated with increased risk of prostate cancer (PC). Novel functional polymorphisms of the UGT2B17 and UGT2B15 genes have been recently characterized by in vitro assays but have not been evaluated in epidemiologic studies. METHODS: Fifteen functional SNPs of the UGT2B17 and UGT2B15 genes, including cis-acting UGT2B gene SNPs, were genotyped in African American and Caucasian men (233 PC cases and 342 controls). Regression models were used to analyze the association between SNPs and PC risk. RESULTS: After adjusting for race, age and BMI, we found that six UGT2B15 SNPs (rs4148269, rs3100, rs9994887, rs13112099, rs7686914 and rs7696472) were associated with an increased risk of PC in log-additive models (p < 0.05). A SNP cis-acting on UGT2B17 and UGT2B15 expression (rs17147338) was also associated with increased risk of prostate cancer (OR = 1.65, 95% CI = 1.00-2.70); while a stronger association among men with high Gleason sum was observed for SNPs rs4148269 and rs3100. CONCLUSIONS: Although small sample size limits inference, we report novel associations between UGT2B15 and UGT2B17 variants and PC risk. These associations with PC risk in men with high Gleason sum, more frequently found in African American men, support the relevance of genetic differences in the androgen metabolism pathway, which could explain, in part, the high incidence of PC among African American men. Larger studies are required.

New tools provide a second look at HDV ribozyme structure, dynamics and cleavage.

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme essential for processing viral transcripts during rolling circle viral replication. The first crystal structure of the cleaved ribozyme was solved in 1998, followed by structures of uncleaved, mutant-inhibited and ion-complexed forms. Recently, methods have been developed that make the task of modeling RNA structure and dynamics significantly easier and more reliable. We have used ERRASER and PHENIX to rebuild and re-refine the cleaved and cis-acting C75U-inhibited structures of the HDV ribozyme. The results correct local conformations and identify alternates for RNA residues, many in functionally important regions, leading to improved R values and model validation statistics for both structures. We compare the rebuilt structures to a higher resolution, trans-acting deoxy-inhibited structure of the ribozyme, and conclude that although both inhibited structures are consistent with the currently accepted hammerhead-like mechanism of cleavage, they do not add direct structural evidence to the biochemical and modeling data. However, the rebuilt structures (PDBs: 4PR6, 4PRF) provide a more robust starting point for research on the dynamics and catalytic mechanism of the HDV ribozyme and demonstrate the power of new techniques to make significant improvements in RNA structures that impact biologically relevant conclusions.

Punctuated evolution and transitional hybrid network in an ancestral cell cycle of fungi.

Although cell cycle control is an ancient, conserved, and essential process, some core animal and fungal cell cycle regulators share no more sequence identity than non-homologous proteins. Here, we show that evolution along the fungal lineage was punctuated by the early acquisition and entrainment of the SBF transcription factor through horizontal gene transfer. Cell cycle evolution in the fungal ancestor then proceeded through a hybrid network containing both SBF and its ancestral animal counterpart E2F, which is still maintained in many basal fungi. We hypothesize that a virally-derived SBF may have initially hijacked cell cycle control by activating transcription via the cis-regulatory elements targeted by the ancestral cell cycle regulator E2F, much like extant viral oncogenes. Consistent with this hypothesis, we show that SBF can regulate promoters with E2F binding sites in budding yeast.

Punctuated evolution and transitional hybrid network in an ancestral cell cycle of fungi

© Medina et al.Although cell cycle control is an ancient, conserved, and essential process, some core animal and fungal cell cycle regulators share no more sequence identity than non-homologous proteins. Here, we show that evolution along the fungal lineage was punctuated by the early acquisition and entrainment of the SBF transcription factor through horizontal gene transfer. Cell cycle evolution in the fungal ancestor then proceeded through a hybrid network containing both SBF and its ancestral animal counterpart E2F, which is still maintained in many basal fungi. We hypothesize that a virally-derived SBF may have initially hijacked cell cycle control by activating transcription via the cis-regulatory elements targeted by the ancestral cell cycle regulator E2F, much like extant viral oncogenes. Consistent with this hypothesis, we show that SBF can regulate promoters with E2F binding sites in budding yeast.

Modeling the evolution of regulatory elements by simultaneous detection and alignment with phylogenetic pair HMMs.

The computational detection of regulatory elements in DNA is a difficult but important problem impacting our progress in understanding the complex nature of eukaryotic gene regulation. Attempts to utilize cross-species conservation for this task have been hampered both by evolutionary changes of functional sites and poor performance of general-purpose alignment programs when applied to non-coding sequence. We describe a new and flexible framework for modeling binding site evolution in multiple related genomes, based on phylogenetic pair hidden Markov models which explicitly model the gain and loss of binding sites along a phylogeny. We demonstrate the value of this framework for both the alignment of regulatory regions and the inference of precise binding-site locations within those regions. As the underlying formalism is a stochastic, generative model, it can also be used to simulate the evolution of regulatory elements. Our implementation is scalable in terms of numbers of species and sequence lengths and can produce alignments and binding-site predictions with accuracy rivaling or exceeding current systems that specialize in only alignment or only binding-site prediction. We demonstrate the validity and power of various model components on extensive simulations of realistic sequence data and apply a specific model to study Drosophila enhancers in as many as ten related genomes and in the presence of gain and loss of binding sites. Different models and modeling assumptions can be easily specified, thus providing an invaluable tool for the exploration of biological hypotheses that can drive improvements in our understanding of the mechanisms and evolution of gene regulation.

Four-decade responses of soil trace elements to an aggrading old-field forest: B, Mn, Zn, Cu, and Fe.

In the ancient and acidic Ultisol soils of the Southern Piedmont, USA, we studied changes in trace element biogeochemistry over four decades, a period during which formerly cultivated cotton fields were planted with pine seedlings that grew into mature forest stands. In 16 permanent plots, we estimated 40-year accumulations of trace elements in forest biomass and O horizons (between 1957 and 1997), and changes in bioavailable soil fractions indexed by extractions of 0.05 mol/L HCl and 0.2 mol/L acid ammonium oxalate (AAO). Element accumulations in 40-year tree biomass plus O horizons totaled 0.9, 2.9, 4.8, 49.6, and 501.3 kg/ha for Cu, B, Zn, Mn, and Fe, respectively. In response to this forest development, samples of the upper 0.6-m of mineral soil archived in 1962 and 1997 followed one of three patterns. (1) Extractable B and Mn were significantly depleted, by -4.1 and -57.7 kg/ha with AAO, depletions comparable to accumulations in biomass plus O horizons, 2.9 and 49.6 kg/ha, respectively. Tree uptake of B and Mn from mineral soil greatly outpaced resupplies from atmospheric deposition, mineral weathering, and deep-root uptake. (2) Extractable Zn and Cu changed little during forest growth, indicating that nutrient resupplies kept pace with accumulations by the aggrading forest. (3) Oxalate-extractable Fe increased substantially during forest growth, by 275.8 kg/ha, about 10-fold more than accumulations in tree biomass (28.7 kg/ha). The large increases in AAO-extractable Fe in surficial 0.35-m mineral soils were accompanied by substantial accretions of Fe in the forest's O horizon, by 473 kg/ha, amounts that dwarfed inputs via litterfall and canopy throughfall, indicating that forest Fe cycling is qualitatively different from that of other macro- and micronutrients. Bioturbation of surficial forest soil layers cannot account for these fractions and transformations of Fe, and we hypothesize that the secondary forest's large inputs of organic additions over four decades has fundamentally altered soil Fe oxides, potentially altering the bioavailability and retention of macro- and micronutrients, contaminants, and organic matter itself. The wide range of responses among the ecosystem's trace elements illustrates the great dynamics of the soil system over time scales of decades.

Chromatin accessibility mapping identifies mediators of basal transcription and retinoid-induced repression of OTX2 in medulloblastoma.

Despite an emerging understanding of the genetic alterations giving rise to various tumors, the mechanisms whereby most oncogenes are overexpressed remain unclear. Here we have utilized an integrated approach of genomewide regulatory element mapping via DNase-seq followed by conventional reporter assays and transcription factor binding site discovery to characterize the transcriptional regulation of the medulloblastoma oncogene Orthodenticle Homeobox 2 (OTX2). Through these studies we have revealed that OTX2 is differentially regulated in medulloblastoma at the level of chromatin accessibility, which is in part mediated by DNA methylation. In cell lines exhibiting chromatin accessibility of OTX2 regulatory regions, we found that autoregulation maintains OTX2 expression. Comparison of medulloblastoma regulatory elements with those of the developing brain reveals that these tumors engage a developmental regulatory program to drive OTX2 transcription. Finally, we have identified a transcriptional regulatory element mediating retinoid-induced OTX2 repression in these tumors. This work characterizes for the first time the mechanisms of OTX2 overexpression in medulloblastoma. Furthermore, this study establishes proof of principle for applying ENCODE datasets towards the characterization of upstream trans-acting factors mediating expression of individual genes.

The upstream enhancer elements of the G6PC promoter are critical for optimal G6PC expression in murine glycogen storage disease type Ia.

Glycogen storage disease type-Ia (GSD-Ia) patients deficient in glucose-6-phosphatase-α (G6Pase-α or G6PC) manifest impaired glucose homeostasis characterized by fasting hypoglycemia, growth retardation, hepatomegaly, nephromegaly, hyperlipidemia, hyperuricemia, and lactic acidemia. Two efficacious recombinant adeno-associated virus pseudotype 2/8 (rAAV8) vectors expressing human G6Pase-α have been independently developed. One is a single-stranded vector containing a 2864-bp of the G6PC promoter/enhancer (rAAV8-GPE) and the other is a double-stranded vector containing a shorter 382-bp minimal G6PC promoter/enhancer (rAAV8-miGPE). To identify the best construct, a direct comparison of the rAAV8-GPE and the rAAV8-miGPE vectors was initiated to determine the best vector to take forward into clinical trials. We show that the rAAV8-GPE vector directed significantly higher levels of hepatic G6Pase-α expression, achieved greater reduction in hepatic glycogen accumulation, and led to a better toleration of fasting in GSD-Ia mice than the rAAV8-miGPE vector. Our results indicated that additional control elements in the rAAV8-GPE vector outweigh the gains from the double-stranded rAAV8-miGPE transduction efficiency, and that the rAAV8-GPE vector is the current choice for clinical translation in human GSD-Ia.