Matthew L. Eaton

Identification of functional elements and regulatory circuits by Drosophila modENCODE

From Genome to Regulatory Networks For biologists, having a genome in hand is only the beginning—much more investigation is still needed to characterize how the genome is used to help to produce a functional organism (see the Perspective by Blaxter). In this vein, Gerstein et al. (p. 1775) summarize for the Caenorhabditis elegans genome, and The modENCODE Consortium (p. 1787) summarize for the Drosophila melanogaster genome, full transcriptome analyses over developmental stages, genome-wide identification of transcription factor binding sites, and high-resolution maps of chromatin organization. Both studies identified regions of the nematode and fly genomes that show highly occupied targets (or HOT) regions where DNA was bound by more than 15 of the transcription factors analyzed and the expression of related genes were characterized. Overall, the studies provide insights into the organization, structure, and function of the two genomes and provide basic information needed to guide and correlate both focused and genome-wide studies. The Drosophila modENCODE project demonstrates the functional regulatory network of flies. To gain insight into how genomic information is translated into cellular and developmental programs, the Drosophila model organism Encyclopedia of DNA Elements (modENCODE) project is comprehensively mapping transcripts, histone modifications, chromosomal proteins, transcription factors, replication proteins and intermediates, and nucleosome properties across a developmental time course and in multiple cell lines. We have generated more than 700 data sets and discovered protein-coding, noncoding, RNA regulatory, replication, and chromatin elements, more than tripling the annotated portion of the Drosophila genome. Correlated activity patterns of these elements reveal a functional regulatory network, which predicts putative new functions for genes, reveals stage- and tissue-specific regulators, and enables gene-expression prediction. Our results provide a foundation for directed experimental and computational studies in Drosophila and related species and also a model for systematic data integration toward comprehensive genomic and functional annotation.

A cis-regulatory map of the Drosophila genome

Systematic annotation of gene regulatory elements is a major challenge in genome science. Direct mapping of chromatin modification marks and transcriptional factor binding sites genome-wide has successfully identified specific subtypes of regulatory elements. In Drosophila several pioneering studies have provided genome-wide identification of Polycomb response elements, chromatin states, transcription factor binding sites, RNA polymerase II regulation and insulator elements; however, comprehensive annotation of the regulatory genome remains a significant challenge. Here we describe results from the modENCODE cis-regulatory annotation project. We produced a map of the Drosophila melanogaster regulatory genome on the basis of more than 300 chromatin immunoprecipitation data sets for eight chromatin features, five histone deacetylases and thirty-eight site-specific transcription factors at different stages of development. Using these data we inferred more than 20,000 candidate regulatory elements and validated a subset of predictions for promoters, enhancers and insulators in vivo. We identified also nearly 2,000 genomic regions of dense transcription factor binding associated with chromatin activity and accessibility. We discovered hundreds of new transcription factor co-binding relationships and defined a transcription factor network with over 800 potential regulatory relationships.

Conserved nucleosome positioning defines replication origins

The origin recognition complex (ORC) specifies replication origin location. The Saccharomyces cerevisiae ORC recognizes the ARS (autonomously replicating sequence) consensus sequence (ACS), but only a subset of potential genomic sites are bound, suggesting other chromosomal features influence ORC binding. Using high-throughput sequencing to map ORC binding and nucleosome positioning, we show that yeast origins are characterized by an asymmetric pattern of positioned nucleosomes flanking the ACS. The origin sequences are sufficient to maintain a nucleosome-free origin; however, ORC is required for the precise positioning of nucleosomes flanking the origin. These findings identify local nucleosomes as an important determinant for origin selection and function.

Chromatin signatures of the Drosophila replication program.

DNA replication initiates from thousands of start sites throughout the Drosophila genome and must be coordinated with other ongoing nuclear processes such as transcription to ensure genetic and epigenetic inheritance. Considerable progress has been made toward understanding how chromatin modifications regulate the transcription program; in contrast, we know relatively little about the role of the chromatin landscape in defining how start sites of DNA replication are selected and regulated. Here, we describe the Drosophila replication program in the context of the chromatin and transcription landscape for multiple cell lines using data generated by the modENCODE consortium. We find that while the cell lines exhibit similar replication programs, there are numerous cell line-specific differences that correlate with changes in the chromatin architecture. We identify chromatin features that are associated with replication timing, early origin usage, and ORC binding. Primary sequence, activating chromatin marks, and DNA-binding proteins (including chromatin remodelers) contribute in an additive manner to specify ORC-binding sites. We also generate accurate and predictive models from the chromatin data to describe origin usage and strength between cell lines. Multiple activating chromatin modifications contribute to the function and relative strength of replication origins, suggesting that the chromatin environment does not regulate origins of replication as a simple binary switch, but rather acts as a tunable rheostat to regulate replication initiation events.

WNT11 expression is induced by ERRα and β-catenin and acts in an autocrine manner to increase cancer cell migration

Elevated expression of the orphan nuclear receptor estrogen-related receptor α (ERRα) has been associated with a negative outcome in several cancers, although the mechanism(s) by which this receptor influences the pathophysiology of this disease and how its activity is regulated remain unknown. Using a chemical biology approach, it was determined that compounds, previously shown to inhibit canonical Wnt signaling, also inhibited the transcriptional activity of ERRα. The significance of this association was revealed in a series of biochemical and genetic experiments that show that (a) ERRα, β-catenin (β-cat), and lymphoid enhancer-binding factor-1 form macromolecular complexes in cells, (b) ERRα transcriptional activity is enhanced by β-cat expression and vice versa, and (c) there is a high level of overlap among genes previously shown to be regulated by ERRα or β-cat. Furthermore, silencing of ERRα and β-cat expression individually or together dramatically reduced the migratory capacity of breast, prostate, and colon cancer cells in vitro. This increased migration could be attributed to the ERRα/β-cat-dependent induction of WNT11. Specifically, using (a) conditioned medium from cells overexpressing recombinant WNT11 or (b) WNT11 neutralizing antibodies, we were able to show that this protein was the key mediator of the promigratory activities of ERRα/β-cat. Together, these data provide evidence for an autocrine regulatory loop involving transcriptional upregulation of WNT11 by ERRα and β-cat that influences the migratory capacity of cancer cells.

Developmental control of gene copy number by repression of replication initiation and fork progression

Precise DNA replication is crucial for genome maintenance, yet this process has been inherently difficult to study on a genome-wide level in untransformed differentiated metazoan cells. To determine how metazoan DNA replication can be repressed, we examined regions selectively under-replicated in Drosophila polytene salivary glands, and found they are transcriptionally silent and enriched for the repressive H3K27me3 mark. In the first genome-wide analysis of binding of the origin recognition complex (ORC) in a differentiated metazoan tissue, we find that ORC binding is dramatically reduced within these large domains, suggesting reduced initiation as one mechanism leading to under-replication. Inhibition of replication fork progression by the chromatin protein SUUR is an additional repression mechanism to reduce copy number. Although repressive histone marks are removed when SUUR is mutated and copy number restored, neither transcription nor ORC binding is reinstated. Tethering of the SUUR protein to a specific site is insufficient to block replication, however. These results establish that developmental control of DNA replication, at both the initiation and elongation stages, is a mechanism to change gene copy number during differentiation.

The conserved bromo-adjacent homology domain of yeast Orc1 functions in the selection of DNA replication origins within chromatin

The origin recognition complex (ORC) binds to the specific positions on chromosomes that serve as DNA replication origins. Although ORC is conserved from yeast to humans, the DNA sequence elements that specify ORC binding are not. In particular, metazoan ORC shows no obvious DNA sequence specificity, whereas yeast ORC binds to a specific DNA sequence within all yeast origins. Thus, whereas chromatin must play an important role in metazoan ORCs ability to recognize origins, it is unclear whether chromatin plays a role in yeast ORCs recognition of origins. This study focused on the role of the conserved N-terminal bromo-adjacent homology domain of yeast Orc1 (Orc1BAH). Recent studies indicate that BAH domains are chromatin-binding modules. We show that the Orc1BAH domain was necessary for ORCs stable association with yeast chromosomes, and was physiologically relevant to DNA replication in vivo. This replication role was separable from the Orc1BAH domains previously defined role in transcriptional silencing. Genome-wide analyses of ORC binding in ORC1 and orc1bahDelta cells revealed that the Orc1BAH domain contributed to ORCs association with most yeast origins, including a class of origins highly dependent on the Orc1BAH domain for ORC association (orc1bahDelta-sensitive origins). Orc1bahDelta-sensitive origins required the Orc1BAH domain for normal activity on chromosomes and plasmids, and were associated with a distinct local nucleosome structure. These data provide molecular insights into how the Orc1BAH domain contributes to ORCs selection of replication origins, as well as new tools for examining conserved mechanisms governing ORCs selection of origins within eukaryotic chromosomes.

DNA copy number evolution in Drosophila cell lines

BackgroundStructural rearrangements of the genome resulting in genic imbalance due to copy number change are often deleterious at the organismal level, but are common in immortalized cell lines and tumors, where they may be an advantage to cells. In order to explore the biological consequences of copy number changes in the Drosophila genome, we resequenced the genomes of 19 tissue-culture cell lines and generated RNA-Seq profiles.ResultsOur work revealed dramatic duplications and deletions in all cell lines. We found three lines of evidence indicating that copy number changes were due to selection during tissue culture. First, we found that copy numbers correlated to maintain stoichiometric balance in protein complexes and biochemical pathways, consistent with the gene balance hypothesis. Second, while most copy number changes were cell line-specific, we identified some copy number changes shared by many of the independent cell lines. These included dramatic recurrence of increased copy number of the PDGF/VEGF receptor, which is also over-expressed in many cancer cells, and of bantam, an anti-apoptosis miRNA. Third, even when copy number changes seemed distinct between lines, there was strong evidence that they supported a common phenotypic outcome. For example, we found that proto-oncogenes were over-represented in one cell line (S2-DRSC), whereas tumor suppressor genes were under-represented in another (Kc167).ConclusionOur study illustrates how genome structure changes may contribute to selection of cell lines in vitro. This has implications for other cell-level natural selection progressions, including tumorigenesis.

Multimodal Regulation of E2F1 Gene Expression by Progestins

ABSTRACT An analysis of mRNA expression in T47D breast cancer cells treated with the synthetic progestin R5020 revealed a subset of progesterone receptor (PR) target genes that are enriched for E2F binding sites. Following up on this observation, we determined that PR-B acts in both direct and indirect manners to positively upregulate E2F1 expression in T47D cells. The direct effects of PR on E2F1 expression were confirmed by chromatin immunoprecipitation (ChIP) analysis, which indicated that the agonist-bound receptor was recruited to several enhancer elements proximal to the E2F1 transcript. However, we also noted that cycloheximide partially inhibits R5020 induction of E2F1 expression, indicating that the ligand-dependent actions of PR on this gene may involve additional indirect regulatory pathways. In support of this hypothesis, we demonstrated that treatment with R5020 significantly increases both hyperphosphorylation of Rb and recruitment of E2F1 to its own promoter, thus activating a positive feedback loop that further amplifies its transcription. Furthermore, we established that PR-mediated induction of Krüppel-like factor 15 (KLF15), which can bind to GC-rich DNA within the E2F1 promoter, is required for maximal induction of E2F1 expression by progestins. Taken together, these results suggest a new paradigm for multimodal regulation of target gene expression by PR.

Challenging the spliceosome machine

BackgroundUsing cDNA copies of transcripts and corresponding genomic sequences from the Berkeley Drosophila Genome Project, a set of 24,753 donor and acceptor splice sites were computed with a scanning algorithm that tested for single nucleotide insertion, deletion and substitution polymorphisms. Using this dataset, we developed a progressive partitioning approach to examining the effects of challenging the spliceosome system.ResultsOur analysis shows that information content increases near splice sites flanking progressively longer introns and exons, suggesting that longer splice elements require stronger binding of spliceosome components. Information also increases at splice sites near very short introns and exons, suggesting that short splice elements have crowding problems. We observe that the information found at individual splice sites depends upon a balance of splice element lengths in the vicinity, including both flanking and non-adjacent introns and exons.ConclusionThese results suggest an interdependence of multiple splicing events along the pre-mRNA, which may have implications for how the macromolecular spliceosome machine processes sets of neighboring splice sites.