Christopher D. Brown

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.

Quorum-Sensing Genes in Pseudomonas aeruginosa Biofilms: Their Role and Expression Patterns

ABSTRACT Acylated homoserine lactone molecules are used by a number of gram-negative bacteria to regulate cell density-dependent gene expression by a mechanism known as quorum sensing (QS). InPseudomonas aeruginosa, QS or cell-to-cell signaling controls expression of a number of virulence factors, as well as biofilm differentiation. In this study, we investigated the role played by the las and rhl QS systems during the early stages of static biofilm formation when cells are adhering to a surface and forming microcolonies. These studies revealed a marked difference in biofilm formation between the PAO1 parent and the QS mutants when glucose, but not citrate, was used as the sole carbon source. To further elucidate the contribution of lasI andrhlI to biofilm maturation, we utilized fusions to unstable green fluorescent protein in concert with confocal microscopy to perform real-time temporal and spatial studies of these genes in a flowing environment. During the course of 8-day biofilm development,lasI expression was found to progressively decrease over time. Conversely, rhlI expression remained steady throughout biofilm development but occurred in a lower percentage of cells. Spatial analysis revealed that lasI andrhlI were maximally expressed in cells located at the substratum and that expression decreased with increasing biofilm height. Because QS was shown previously to be involved in biofilm differentiation, these findings have important implications for the design of biofilm prevention and eradication strategies.

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.

A Comprehensive Map of Insulator Elements for the Drosophila Genome

Insulators are DNA sequences that control the interactions among genomic regulatory elements and act as chromatin boundaries. A thorough understanding of their location and function is necessary to address the complexities of metazoan gene regulation. We studied by ChIP–chip the genome-wide binding sites of 6 insulator-associated proteins—dCTCF, CP190, BEAF-32, Su(Hw), Mod(mdg4), and GAF—to obtain the first comprehensive map of insulator elements in Drosophila embryos. We identify over 14,000 putative insulators, including all classically defined insulators. We find two major classes of insulators defined by dCTCF/CP190/BEAF-32 and Su(Hw), respectively. Distributional analyses of insulators revealed that particular sub-classes of insulator elements are excluded between cis-regulatory elements and their target promoters; divide differentially expressed, alternative, and divergent promoters; act as chromatin boundaries; are associated with chromosomal breakpoints among species; and are embedded within active chromatin domains. Together, these results provide a map demarcating the boundaries of gene regulatory units and a framework for understanding insulator function during the development and evolution of Drosophila.

Integrative and Comparative Genomic Analysis of HPV-Positive and HPV-Negative Head and Neck Squamous Cell Carcinomas

Purpose: The genetic differences between human papilloma virus (HPV)–positive and –negative head and neck squamous cell carcinomas (HNSCC) remain largely unknown. To identify differential biology and novel therapeutic targets for both entities, we determined mutations and copy-number aberrations in a large cohort of locoregionally advanced HNSCC. Experimental Design: We performed massively parallel sequencing of 617 cancer-associated genes in 120 matched tumor/normal samples (42.5% HPV-positive). Mutations and copy-number aberrations were determined and results validated with a secondary method. Results: The overall mutational burden in HPV-negative and HPV-positive HNSCC was similar with an average of 15.2 versus 14.4 somatic exonic mutations in the targeted cancer-associated genes. HPV-negative tumors showed a mutational spectrum concordant with published lung squamous cell carcinoma analyses with enrichment for mutations in TP53, CDKN2A, MLL2, CUL3, NSD1, PIK3CA, and NOTCH genes. HPV-positive tumors showed unique mutations in DDX3X, FGFR2/3 and aberrations in PIK3CA, KRAS, MLL2/3, and NOTCH1 were enriched in HPV-positive tumors. Currently targetable genomic alterations were identified in FGFR1, DDR2, EGFR, FGFR2/3, EPHA2, and PIK3CA. EGFR, CCND1, and FGFR1 amplifications occurred in HPV-negative tumors, whereas 17.6% of HPV-positive tumors harbored mutations in fibroblast growth factor receptor genes (FGFR2/3), including six recurrent FGFR3 S249C mutations. HPV-positive tumors showed a 5.8% incidence of KRAS mutations, and DNA-repair gene aberrations, including 7.8% BRCA1/2 mutations, were identified. Conclusions: The mutational makeup of HPV-positive and HPV-negative HNSCC differs significantly, including targetable genes. HNSCC harbors multiple therapeutically important genetic aberrations, including frequent aberrations in the FGFR and PI3K pathway genes. Clin Cancer Res; 21(3); 632–41. ©2014 AACR. See related commentary by Krigsfeld and Chung, p. 495

Genetic effects on gene expression across human tissues

Characterization of the molecular function of the human genome and its variation across individuals is essential for identifying the cellular mechanisms that underlie human genetic traits and diseases. The Genotype-Tissue Expression (GTEx) project aims to characterize variation in gene expression levels across individuals and diverse tissues of the human body, many of which are not easily accessible. Here we describe genetic effects on gene expression levels across 44 human tissues. We find that local genetic variation affects gene expression levels for the majority of genes, and we further identify inter-chromosomal genetic effects for 93 genes and 112 loci. On the basis of the identified genetic effects, we characterize patterns of tissue specificity, compare local and distal effects, and evaluate the functional properties of the genetic effects. We also demonstrate that multi-tissue, multi-individual data can be used to identify genes and pathways affected by human disease-associated variation, enabling a mechanistic interpretation of gene regulation and the genetic basis of disease.

A statin-dependent QTL for GATM expression is associated with statin-induced myopathy.

Statins are prescribed widely to lower plasma low-density lipoprotein (LDL) concentrations and cardiovascular disease risk and have been shown to have beneficial effects in a broad range of patients. However, statins are associated with an increased risk, albeit small, of clinical myopathy and type 2 diabetes. Despite evidence for substantial genetic influence on LDL concentrations, pharmacogenomic trials have failed to identify genetic variations with large effects on either statin efficacy or toxicity, and have produced little information regarding mechanisms that modulate statin response. Here we identify a downstream target of statin treatment by screening for the effects of in vitro statin exposure on genetic associations with gene expression levels in lymphoblastoid cell lines derived from 480 participants of a clinical trial of simvastatin treatment. This analysis identified six expression quantitative trait loci (eQTLs) that interacted with simvastatin exposure, including rs9806699, a cis-eQTL for the gene glycine amidinotransferase (GATM) that encodes the rate-limiting enzyme in creatine synthesis. We found this locus to be associated with incidence of statin-induced myotoxicity in two separate populations (meta-analysis odds ratio = 0.60). Furthermore, we found that GATM knockdown in hepatocyte-derived cell lines attenuated transcriptional response to sterol depletion, demonstrating that GATM may act as a functional link between statin-mediated lowering of cholesterol and susceptibility to statin-induced myopathy.

Integrative Analysis of Head and Neck Cancer Identifies Two Biologically Distinct HPV and Three Non-HPV Subtypes

Purpose: Current classification of head and neck squamous cell carcinomas (HNSCC) based on anatomic site and stage fails to capture biologic heterogeneity or adequately inform treatment. Experimental Design: Here, we use gene expression-based consensus clustering, copy number profiling, and human papillomavirus (HPV) status on a clinically homogenous cohort of 134 locoregionally advanced HNSCCs with 44% HPV+ tumors together with additional cohorts, which in total comprise 938 tumors, to identify HNSCC subtypes and discover several subtype-specific, translationally relevant characteristics. Results: We identified five subtypes of HNSCC, including two biologically distinct HPV subtypes. One HPV+ and one HPV− subtype show a prominent immune and mesenchymal phenotype. Prominent tumor infiltration with CD8+ lymphocytes characterizes this inflamed/mesenchymal subtype, independent of HPV status. Compared with other subtypes, the two HPV subtypes show low expression and no copy number events for EGFR/HER ligands. In contrast, the basal subtype is uniquely characterized by a prominent EGFR/HER signaling phenotype, negative HPV-status, as well as strong hypoxic differentiation not seen in other subtypes. Conclusion: Our five-subtype classification provides a comprehensive overview of HPV+ as well as HPV− HNSCC biology with significant translational implications for biomarker development and personalized care for patients with HNSCC. Clin Cancer Res; 21(4); 870–81. ©2014 AACR.

Functional Architecture and Evolution of Transcriptional Elements That Drive Gene Coexpression

Transcriptional coexpression of interacting gene products is required for complex molecular processes; however, the function and evolution of cis-regulatory elements that orchestrate coexpression remain largely unexplored. We mutagenized 19 regulatory elements that drive coexpression of Ciona muscle genes and obtained quantitative estimates of the cis-regulatory activity of the 77 motifs that comprise these elements. We found that individual motif activity ranges broadly within and among elements, and among different instantiations of the same motif type. The activity of orthologous motifs is strongly constrained, although motif arrangement, type, and activity vary greatly among the elements of different co-regulated genes. Thus, the syntactical rules governing this regulatory function are flexible but become highly constrained evolutionarily once they are established in a particular element.

Rapid growth of a hepatocellular carcinoma and the driving mutations revealed by cell-population genetic analysis of whole-genome data

We present the analysis of the evolution of tumors in a case of hepatocellular carcinoma. This case is particularly informative about cancer growth dynamics and the underlying driving mutations. We sampled nine different sections from three tumors and seven more sections from the adjacent nontumor tissues. Selected sections were subjected to exon as well as whole-genome sequencing. Putative somatic mutations were then individually validated across all 9 tumor and 7 nontumor sections. Among the mutations validated, 24 were amino acid changes; in addition, 22 large indels/copy number variants (>1 Mb) were detected. These somatic mutations define four evolutionary lineages among tumor cells. Separate evolution and expansion of these lineages were recent and rapid, each apparently having only one lineage-specific protein-coding mutation. Hence, by using a cell-population genetic definition, this approach identified three coding changes (CCNG1, P62, and an indel/fusion gene) as tumor driver mutations. These three mutations, affecting cell cycle control and apoptosis, are functionally distinct from mutations that accumulated earlier, many of which are involved in inflammation/immunity or cell anchoring. These distinct functions of mutations at different stages may reflect the genetic interactions underlying tumor growth.