Human-chimpanzee differences in a FZD8 enhancer alter cell-cycle dynamics in the developing neocortex.

The human neocortex differs from that of other great apes in several notable regards, including altered cell cycle, prolonged corticogenesis, and increased size [1-5]. Although these evolutionary changes most likely contributed to the origin of distinctively human cognitive faculties, their genetic basis remains almost entirely unknown. Highly conserved non-coding regions showing rapid sequence changes along the human lineage are candidate loci for the development and evolution of uniquely human traits. Several studies have identified human-accelerated enhancers [6-14], but none have linked an expression difference to a specific organismal trait. Here we report the discovery of a human-accelerated regulatory enhancer (HARE5) of FZD8, a receptor of the Wnt pathway implicated in brain development and size [15, 16]. Using transgenic mice, we demonstrate dramatic differences in human and chimpanzee HARE5 activity, with human HARE5 driving early and robust expression at the onset of corticogenesis. Similar to HARE5 activity, FZD8 is expressed in neural progenitors of the developing neocortex [17-19]. Chromosome conformation capture assays reveal that HARE5 physically and specifically contacts the core Fzd8 promoter in the mouse embryonic neocortex. To assess the phenotypic consequences of HARE5 activity, we generated transgenic mice in which Fzd8 expression is under control of orthologous enhancers (Pt-HARE5::Fzd8 and Hs-HARE5::Fzd8). In comparison to Pt-HARE5::Fzd8, Hs-HARE5::Fzd8 mice showed marked acceleration of neural progenitor cell cycle and increased brain size. Changes in HARE5 function unique to humans thus alter the cell-cycle dynamics of a critical population of stem cells during corticogenesis and may underlie some distinctive anatomical features of the human brain.

Pol II docking and pausing at growth and stress genes in C. elegans.

Fluctuations in nutrient availability profoundly impact gene expression. Previous work revealed postrecruitment regulation of RNA polymerase II (Pol II) during starvation and recovery in Caenorhabditis elegans, suggesting that promoter-proximal pausing promotes rapid response to feeding. To test this hypothesis, we measured Pol II elongation genome wide by two complementary approaches and analyzed elongation in conjunction with Pol II binding and expression. We confirmed bona fide pausing during starvation and also discovered Pol II docking. Pausing occurs at active stress-response genes that become downregulated in response to feeding. In contrast, "docked" Pol II accumulates without initiating upstream of inactive growth genes that become rapidly upregulated upon feeding. Beyond differences in function and expression, these two sets of genes have different core promoter motifs, suggesting alternative transcriptional machinery. Our work suggests that growth and stress genes are both regulated postrecruitment during starvation but at initiation and elongation, respectively, coordinating gene expression with nutrient availability.

Interactions of chromatin context, binding site sequence content, and sequence evolution in stress-induced p53 occupancy and transactivation.

Cellular stresses activate the tumor suppressor p53 protein leading to selective binding to DNA response elements (REs) and gene transactivation from a large pool of potential p53 REs (p53REs). To elucidate how p53RE sequences and local chromatin context interact to affect p53 binding and gene transactivation, we mapped genome-wide binding localizations of p53 and H3K4me3 in untreated and doxorubicin (DXR)-treated human lymphoblastoid cells. We examined the relationships among p53 occupancy, gene expression, H3K4me3, chromatin accessibility (DNase 1 hypersensitivity, DHS), ENCODE chromatin states, p53RE sequence, and evolutionary conservation. We observed that the inducible expression of p53-regulated genes was associated with the steady-state chromatin status of the cell. Most highly inducible p53-regulated genes were suppressed at baseline and marked by repressive histone modifications or displayed CTCF binding. Comparison of p53RE sequences residing in different chromatin contexts demonstrated that weaker p53REs resided in open promoters, while stronger p53REs were located within enhancers and repressed chromatin. p53 occupancy was strongly correlated with similarity of the target DNA sequences to the p53RE consensus, but surprisingly, inversely correlated with pre-existing nucleosome accessibility (DHS) and evolutionary conservation at the p53RE. Occupancy by p53 of REs that overlapped transposable element (TE) repeats was significantly higher (p<10-7) and correlated with stronger p53RE sequences (p<10-110) relative to nonTE-associated p53REs, particularly for MLT1H, LTR10B, and Mer61 TEs. However, binding at these elements was generally not associated with transactivation of adjacent genes. Occupied p53REs located in L2-like TEs were unique in displaying highly negative PhyloP scores (predicted fast-evolving) and being associated with altered H3K4me3 and DHS levels. These results underscore the systematic interaction between chromatin status and p53RE context in the induced transactivation response. This p53 regulated response appears to have been tuned via evolutionary processes that may have led to repression and/or utilization of p53REs originating from primate-specific transposon elements.

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.

Inactivation of the von Hippel-Lindau tumor suppressor leads to selective expression of a human endogenous retrovirus in kidney cancer.

A human endogenous retrovirus type E (HERV-E) was recently found to be selectively expressed in most renal cell carcinomas (RCCs). Importantly, antigens derived from this provirus are immunogenic, stimulating cytotoxic T cells that kill RCC cells in vitro and in vivo. Here, we show HERV-E expression is restricted to the clear cell subtype of RCC (ccRCC) characterized by an inactivation of the von Hippel-Lindau (VHL) tumor-suppressor gene with subsequent stabilization of hypoxia-inducible transcription factors (HIFs)-1α and -2α. HERV-E expression in ccRCC linearly correlated with HIF-2α levels and could be silenced in tumor cells by either transfection of normal VHL or small interfering RNA inhibition of HIF-2α. Using chromatin immunoprecipitation, we demonstrated that HIF-2α can serve as transcriptional factor for HERV-E by binding with HIF response element (HRE) localized in the proviral 5' long terminal repeat (LTR). Remarkably, the LTR was found to be hypomethylated only in HERV-E-expressing ccRCC while other tumors and normal tissues possessed a hypermethylated LTR preventing proviral expression. Taken altogether, these findings provide the first evidence that inactivation of a tumor suppressor gene can result in aberrant proviral expression in a human tumor and give insights needed for translational research aimed at boosting human immunity against antigenic components of this HERV-E.

High-coverage sequencing and annotated assemblies of the budgerigar genome.

BACKGROUND: Parrots belong to a group of behaviorally advanced vertebrates and have an advanced ability of vocal learning relative to other vocal-learning birds. They can imitate human speech, synchronize their body movements to a rhythmic beat, and understand complex concepts of referential meaning to sounds. However, little is known about the genetics of these traits. Elucidating the genetic bases would require whole genome sequencing and a robust assembly of a parrot genome. FINDINGS: We present a genomic resource for the budgerigar, an Australian Parakeet (Melopsittacus undulatus) -- the most widely studied parrot species in neuroscience and behavior. We present genomic sequence data that includes over 300× raw read coverage from multiple sequencing technologies and chromosome optical maps from a single male animal. The reads and optical maps were used to create three hybrid assemblies representing some of the largest genomic scaffolds to date for a bird; two of which were annotated based on similarities to reference sets of non-redundant human, zebra finch and chicken proteins, and budgerigar transcriptome sequence assemblies. The sequence reads for this project were in part generated and used for both the Assemblathon 2 competition and the first de novo assembly of a giga-scale vertebrate genome utilizing PacBio single-molecule sequencing. CONCLUSIONS: Across several quality metrics, these budgerigar assemblies are comparable to or better than the chicken and zebra finch genome assemblies built from traditional Sanger sequencing reads, and are sufficient to analyze regions that are difficult to sequence and assemble, including those not yet assembled in prior bird genomes, and promoter regions of genes differentially regulated in vocal learning brain regions. This work provides valuable data and material for genome technology development and for investigating the genomics of complex behavioral traits.

HDAC inhibitors cause site-specific chromatin remodeling at PU.1-bound enhancers in K562 cells

BACKGROUND: Small molecule inhibitors of histone deacetylases (HDACi) hold promise as anticancer agents for particular malignancies. However, clinical use is often confounded by toxicity, perhaps due to indiscriminate hyperacetylation of cellular proteins. Therefore, elucidating the mechanisms by which HDACi trigger differentiation, cell cycle arrest, or apoptosis of cancer cells could inform development of more targeted therapies. We used the myelogenous leukemia line K562 as a model of HDACi-induced differentiation to investigate chromatin accessibility (DNase-seq) and expression (RNA-seq) changes associated with this process. RESULTS: We identified several thousand specific regulatory elements [~10 % of total DNase I-hypersensitive (DHS) sites] that become significantly more or less accessible with sodium butyrate or suberanilohydroxamic acid treatment. Most of the differential DHS sites display hallmarks of enhancers, including being enriched for non-promoter regions, associating with nearby gene expression changes, and increasing luciferase reporter expression in K562 cells. Differential DHS sites were enriched for key hematopoietic lineage transcription factor motifs, including SPI1 (PU.1), a known pioneer factor. We found PU.1 increases binding at opened DHS sites with HDACi treatment by ChIP-seq, but PU.1 knockdown by shRNA fails to block the chromatin accessibility and expression changes. A machine-learning approach indicates H3K27me3 initially marks PU.1-bound sites that open with HDACi treatment, suggesting these sites are epigenetically poised. CONCLUSIONS: We find HDACi treatment of K562 cells results in site-specific chromatin remodeling at epigenetically poised regulatory elements. PU.1 shows evidence of a pioneer role in this process by marking poised enhancers but is not required for transcriptional activation.

Expression of the Ets transcription factor Erm is regulated through a conventional PKC signaling pathway in the Molt4 lymphoblastic cell line.

Erm, a member of the PEA3 group within the Ets family of transcription factors, is expressed in murine and human lymphocytes. Here, we show that in the human Molt4 lymphoblastic cell line, the erm gene expression is regulated by the conventional PKC (cPKC) pathway. To better characterize the molecular mechanism by which cPKC regulates Erm transcription in Molt4 cells, we tested proximal promoter deletions of the human gene, and identified a specific cPKC-regulated region between positions -420 and -115 upstream of the first exon.

Sequence variation in the CHAT locus shows no association with late-onset Alzheimer's disease

There is substantial evidence for a susceptibility gene for late-onset Alzheimer's disease (AD) on chromosome 10. One of the characteristic features of AD is the degeneration and dysfunction of the cholinergic system. The genes encoding choline acetyltransferase (ChAT) and its vesicular transporter (VAChT), CHAT and SLC18A3 respectively, map to the linked region of chromosome 10 and are therefore both positional and obvious functional candidate genes for late-onset AD. We have screened both genes for sequence variants and investigated each for association with late-onset AD in up to 500 late-onset AD cases and 500 control DNAs collected in the UK. We detected a total of 17 sequence variants. Of these, 14 were in CHAT, comprising three non-synonymous variants (D7N in the S exon, A120T in exon 5 and L243F in exon 8), one synonymous change (H547H), nine single-nucleotide polymorphisms in intronic, untranslated or promoter regions, and a variable number of tandem repeats in intron 7. Three non-coding SNPs were detected in SLC18A3. None demonstrated any reproducible association with late-onset AD in our samples. Levels of linkage disequilibrium were generally low across the CHAT locus but two of the coding variants, D7N and A120T, proved to be in complete linkage disequilibrium.

Cytokine gene polymorphisms in ischaemic heart disease:Investigation using family based tests of association

Atherosclerosis has an inflammatory basis, with cytokines, cellular adhesion molecules and pro-inflammatory cells having important roles in the initiation and progression of this process. Interleukin (IL) 6, IL-10 and transforming growth factor (TGF) β have been proposed as important modulators of the atherosclerotic process, with IL-6 having a pro-inflammatory, atherogenic effect and IL-10 and TGF-β having anti-inflammatory, protective roles. The possible role of functional polymorphisms in the promoter regions of the IL-6, IL-10 and TGF-β genes in the susceptibility to ischaemic heart disease (IHD) was investigated in a well-defined Irish population using two recently described family-based tests of association. We genotyped 1,012 individuals from 386 families with at least one member prematurely affected with IHD. Using the combined transmission disequilibrium test (TDT)/sib-TDT and the pedigree disequilibrium test, no association between any of the IL-6 -174G/C, IL-10 -1082G/A and TGF-β -509C/T polymorphisms and IHD was found. Our data demonstrate that, in an Irish population, these polymorphisms are not associated with IHD. © Springer-Verlag 2004.

Comparative analysis of dna methylation profiles in peripheral blood leukocytes versus lymphoblastoid cell lines

Previous reports have shown that DNA methylation profiles within primary human malignant tissues are altered when these cells are transformed into cancer cell lines. However, it is unclear if similar differences in DNA methylation profiles exist between DNA derived from peripheral blood leukocytes (PBLs) and corresponding Epstein-Barr Virus transformed lymphoblastoid cell lines (LCLs). To assess the utility of LCLs as a resource for methylation studies we have compared DNA methylation profiles in promoter and 5' regions of 318 genes in PBL and LCL sample pairs from patients with type 1 diabetes with or without nephropathy. We identified a total of 27 (similar to 8%) genes that revealed different DNA methylation profiles in PBL compared with LCL-derived DNA samples. In conclusion, although the profiles for most promoter regions were similar between PBL-LCL pairs, our results indicate that LCL-derived DNA may not be suitable for DNA methylation studies at least in diabetic nephropathy.

DNA methylation profiling in cell models of diabetic nephropathy

We have previously identified differentially expressed genes in cell models of diabetic nephropathy and renal biopsies. Here we have performed quantitative DNA methylation profiling in cell models of diabetic nephropathy. Over 3,000 CpG units in the promoter regions of 192 candidate genes were assessed in unstimulated human mesangial cells (HMCs) and proximal tubular epithelial cells (PTCs) compared to HMCs or PTCs exposed to appropriate stimuli. A total of 301 CpG units across 38 genes (similar to 20%) were identified as differentially methylated in unstimulated HMCs versus PTCs. Analysis of amplicon methylation values in unstimulated versus stimulated cell models failed to demonstrate a >20% difference between amplicons. In conclusion, our results demonstrate that specific DNA methylation signatures are present in HMCs and PTCs, and standard protocols for exposure of renal cells to stimuli that alter gene expression may be insufficient to replicate possible alterations in DNA methylation profiles in diabetic nephropathy.

BRCA1 transcriptionally regulates genes associated with the basal breast cancer phenotype

Background Ten to twenty per cent of breast tumours exhibit a basallike genetic profile and these tumours carry a poor prognosis. Breast tumours which contain germline mutations for BRCA1 commonly exhibit a molecular profile similar to basal breast tumours. BRCA1 is a tumour suppressor gene which is mutated in up to 5–10% of breast cancer cases and is involved in multiple cellular processes including DNA damage control, cell cycle checkpoint control, apoptosis, ubiquitination and transcriptional regulation.

Methods Microarray-based profiling was carried out using the HCC1937EV and HCC1937BR breast cancer cell lines. Basal gene and protein expression levels were analysed by qRT-PCR and western blotting. ChIP analyses were performed and demonstrated that BRCA1 regulates basal gene expression through a transcriptional mechanism involving c-myc.

Results We have previously carried out microarray-based expression profiling to examine differences in gene expression when BRCA1 is reconstituted in BRCA1 mutated HCC1937 breast cancer cells. We observed that p-cadherin and the cytokeratin 5 and cytokeratin 17 genes, which are strongly correlated with the basal phenotype, are differentially expressed when BRCA1 is reconstituted. In addition, qRT-PCR and ChIP analysis of BRCA1 reconstituted cells show that BRCA1 represses the expression of these basal genes by a transcriptional mechanism. Furthermore, abrogation of endogenous BRCA1 protein in the T47D cell line using siRNA results in reexpression of these basal genes, suggesting that BRCA1 expression levels may be important in basal gene expression. We have also demonstrated that BRCA1 is physically associated with the promoter regions of basal genes through an association with c-myc. Consequently, we have confirmed that siRNA inhibition of c-myc in T47D cells results in re-expression of these genes.

Conclusions Our results suggest that BRCA1 is involved in the transcriptional regulation of genes associated with the basal phenotype and that BRCA1 controls basal gene expression through a transcriptional mechanism involving c-myc. Further work is now concentrating on defining the relationship between BRCA1 and basal gene expression and how this may affect clinical responses to breast cancer chemotherapy.

Identification of a methylation hotspot in the death receptor Fas/CD95 in bladder cancer

We characterized Fas immunoreactivity, functionality and its role in the response to mitomycin-C (MMC) chemotherapy in vitro in cell lines and in vivo in bladder washings from 23 transitional cell carcinoma of the bladder (TCCB) patients, harvested prior to and during MMC intravesical treatment. Having established the importance of functional Fas, we investigated the methylation and exon 9 mutation as mechanisms of Fas silencing in TCCB. For the first time, we report p53 up-regulation in 9/14 and Fas up-regulation in 7/9 TCCB patients during intravesical MMC treatment. Fas immunoreactivity was strong in the TCCB cell line T24 and in 17/20 (85%) tumor samples from patients with advanced TCCB. T24 and HT1376 cells were resistant to MMC and recombinant Fas ligand, whilst RT4 cells were responsive to Fas ligand and MMC. Using RT4 cells as a model, siRNA targeting p53 significantly reduced MMC-induced p53 and Fas up-regulation and stable DN-FADD transfection decreased MMC-induced apoptosis, suggesting that functional Fas enhances chemotherapy responses in a p53-dependent manner. In HT1376 cells, 5-aza-2-deoxycytidine (12 µM) induced Fas immunoreactivity and reversed methylation at CpG site -548 within the Fas promoter. This site was methylated in 13/24 (54%) TCCB patient samples assessed using Methylation-Specific Polymerase Chain Reaction. There was no methylation at either the p53 enhancer region within the first intron or at the SP-1 binding region in the promoter and no mutation within exon 9 in tumor DNA extracted from 38 patients. Methylation at CpG site -548 is a potential target for demethylating drugs.

Genetic analysis of the dTDP-rhamnose biosynthesis region of the Escherichia coli VW187 (O7:K1) rfb gene cluster:identification of functional homologs of rfbB and rfbA in the rff cluster and correct location of the rffE gene

The O-repeating unit of the Escherichia coli O7-specific lipopolysaccharide is made of galactose, mannose, rhamnose, 4-acetamido-4,6-dideoxyglucose, and N-acetyglucosamine. We have recently characterized the genes involved in the biosynthesis of the sugar precursor GDP-mannose occurring in the E. coli O7:K1 strain VW187 (C. L. Marolda and M. A. Valvano, J. Bacteriol. 175:148-158, 1993). In the present study, we identified and sequenced the rfbBDAC genes encoding the enzymes for the biosynthesis of another precursor, dTDP-rhamnose. These genes are localized on the upstream end of the rfbEcO7 region, and they are strongly conserved compared with similar genes found in various enteric and nonenteric bacteria. Upstream of rfbB we identified a DNA segment containing the rfb promoter and a highly conserved untranslated leader sequence also present in the promoter regions of other surface polysaccharide gene clusters. Also, we have determined that rfbB and rfbA have homologs, rffG (o355) and rffH (o292), respectively, located on the rff cluster, which is involved in the synthesis of enterobacterial common antigen. We provide biochemical evidence that rffG and rffH encode dTDP-glucose dehydratase and glucose-1-phosphate thymidylyltransferase activities, respectively, and we also show that rffG complemented the rfbB defect in the O7+ cosmid pJHCV32. We also demonstrate that rffG is distinct from rffE and map the rffE gene to the second gene of the rff cluster.