Phylosymbiosis and metabolomics resolve phenotypically plastic and cryptic sponge species in the genus Agelas across the Caribbean basin

Fundamental to holobiont biology is recognising how variation in microbial composition and function relates to host phenotypic variation. Sponges often exhibit considerable phenotypic plasticity and also harbour dense microbial communities that function to protect and nourish hosts. One of the most prominent sponge genera on Caribbean coral reefs is Agelas. Using a comprehensive set of morphological (growth form, spicule), chemical and molecular data on 13 recognised species of Agelas in the Caribbean basin, we were able to define only five species (=clades) and found that many morphospecies designations were incongruent with phylogenomic and population genetic analyses. Microbial communities were also strongly differentiated between phylogenetic species, showing little evidence of cryptic divergence and relatively low correlation with morphospecies assignment. Metagenomic analyses also showed strong correspondence to phylogenetic species, and to a lesser extent, geographical and morphological characters. Surprisingly, the variation in secondary metabolites produced by sponge holobionts was explained by geography and morphospecies assignment, in addition to phylogenetic species, and covaried significantly with a subset of microbial symbionts. Spicule characteristics were highly plastic, under greater impact from geographical location than phylogeny. Our results suggest that while phenotypic plasticity is rampant in Agelas, morphological differences within phylogenetic species affect functionally important ecological traits, including the composition of the symbiotic microbial communities and metabolomic profiles.

Limited expression of Nrf2 in neurons across the central nervous system

Nrf2, encoded by the gene Nfe2l2, is a broadly expressed transcription factor that regulates gene expression in response to reactive oxygen species (ROS) and oxidative stress. It is commonly referred to as a ubiquitous pathway, but this generalization overlooks work indicating that Nrf2 is essentially unexpressed in some neuronal populations. To explore whether this pattern extends throughout the central nervous system (CNS), we quantified Nfe2l2 expression and chromatin accessibility at the Nfe2l2 locus across multiple single cell datasets. In both the mouse and human CNS, Nfe2l2 was repressed in almost all mature neurons, but highly expressed in non-neuronal support cells, and this pattern was robust across multiple human CNS diseases. A subset of key Nrf2 target genes, like Slc7a11, also remained low in neurons. Thus, these data suggest that while most cells express Nfe2l2, with activity determined by ROS levels, neurons actively avoid Nrf2 activity by keeping Nfe2l2 expression low.

An antioxidant response element regulates the HIF1α axis in breast cancer cells

The redox sensitive transcription factor NRF2 is a central regulator of the transcriptional response to reactive oxygen species (ROS). NRF2 is widely recognized for its ROS-responsive upregulation of antioxidant genes that are essential for mitigating the damaging effects of oxidative stress. However, multiple genome-wide approaches have suggested that NRF2's regulatory reach extends well beyond the canonical antioxidant genes, with the potential to regulate many noncanonical target genes. Recent work from our lab and others suggests HIF1A, which encodes the hypoxia-responsive transcription factor HIF1α, is one such noncanonical NRF2 target. These studies found that NRF2 activity is associated with high HIF1A expression in multiple cellular contexts, HIF1A expression is partially dependent on NRF2, and there is a putative NRF2 binding site (antioxidant response element, or ARE) approximately 30 kilobases upstream of HIF1A. These findings all support a model in which HIF1A is a direct target of NRF2, but did not confirm the functional importance of the upstream ARE in HIF1A expression. Here we use CRISPR/Cas9 genome editing to mutate this ARE in its genomic context and test the impact on HIF1A expression. We find that mutation of this ARE in a breast cancer cell line (MDA-MB-231) eliminates NRF2 binding and decreases HIF1A expression at the transcript and protein levels, and disrupts HIF1α target genes as well as phenotypes driven by these HIF1α targets. Taken together, these results indicate that this NRF2 targeted upstream ARE plays an important role in the expression of HIF1A and activity of the HIF1α axis in MDA-MB-231 cells.

Limited Expression of Nrf2 in Neurons Across the Central Nervous System

Nrf2 is a broadly expressed transcription factor that regulates gene expression in response to reactive oxygen species (ROS) and oxidative stress. It is commonly referred to as a ubiquitous pathway, but this generalization overlooks work indicating that Nrf2 is essentially unexpressed in some neuronal populations. To explore whether this pattern extends throughout the central nervous system (CNS), we quantified Nrf2 expression and chromatin accessibility at the Nrf2 locus across multiple single cell datasets. In both the mouse and human CNS, Nrf2 was repressed in almost all mature neurons, but highly expressed in non-neuronal support cells, and this pattern was robust across multiple human CNS diseases. A subset of key Nrf2 target genes, like Slc7a11 , also remained low in neurons. Thus, these data suggest that while most cells express Nrf2, with activity determined by ROS levels, neurons actively avoid Nrf2 activity by keeping Nrf2 expression low.

Ribbon boosts ribosomal protein gene expression to coordinate organ form and function

Cell growth is well defined for late (postembryonic) stages of development, but evidence for early (embryonic) cell growth during postmitotic morphogenesis is limited. Here, we report early cell growth as a key characteristic of tubulogenesis in the Drosophila embryonic salivary gland (SG) and trachea. A BTB/POZ domain nuclear factor, Ribbon (Rib), mediates this early cell growth. Rib binds the transcription start site of nearly every SG-expressed ribosomal protein gene (RPG) and is required for full expression of all RPGs tested. Rib binding to RPG promoters in vitro is weak and not sequence specific, suggesting that specificity is achieved through cofactor interactions. Accordingly, we demonstrate Rib’s ability to physically interact with each of the three known regulators of RPG transcription. Surprisingly, Rib-dependent early cell growth in another tubular organ, the embryonic trachea, is not mediated by direct RPG transcription. These findings support a model of early cell growth customized by transcriptional regulatory networks to coordinate organ form and function.

Transcriptional reprogramming by oxidative stress occurs within a predefined chromatin accessibility landscape

Reactive oxygen species (ROS) are important signaling molecules in many physiological processes, yet excess ROS leads to cell damage and can lead to pathology. Accordingly, cells need to maintain tight regulation of ROS levels, and ROS-responsive transcriptional reprogramming is central to this process. Although it has long been recognized that oxidative stress leads to rapid, significant changes in gene expression, the impact of oxidative stress on the underlying chromatin accessibility landscape remained unclear. Here, we asked whether ROS-responsive transcriptional reprogramming is accompanied by reprogramming of the chromatin environment in MCF7 human breast cancer cells. Using a time-course exposure to multiple inducers of oxidative stress, we determined that the widespread ROS-responsive changes in gene expression induced by ROS occur with minimal changes to the chromatin environment. While we did observe changes in chromatin accessibility, these changes were: (1) far less numerous than gene expression changes after oxidative stress, and (2) occur within pre-existing regions of accessible chromatin. Transcription factor (TF) footprinting analysis of our ATAC-seq experiments identified 5 TFs or TF families with evidence for ROS-responsive changes in DNA binding: NRF2, AP-1, p53, NFY, and SP/KLF. Importantly, several of these (AP-1, NF-Y, and SP/KLF factors) have not been previously implicated as widespread regulators in the response to ROS. In summary, we have characterized genome-wide changes in gene expression and chromatin accessibility in response to ROS treatment of MCF7 cells, and we have found that regulation of the large-scale transcriptional response to excess ROS is primarily constrained by the cell's pre-existing chromatin landscape.

CrebA increases secretory capacity through direct transcriptional regulation of the secretory machinery, a subset of secretory cargo, and other key regulators

Specialization of many cells, including the acinar cells of the salivary glands and pancreas, milk-producing cells of mammary glands, mucus-secreting goblet cells, antibody-producing plasma cells, and cells that generate the dense extracellular matrices of bone and cartilage, requires scaling up both secretory machinery and cell-type specific secretory cargo. Using tissue-specific genome-scale analyses, we determine how increases in secretory capacity are coordinated with increases in secretory load in the Drosophila salivary gland (SG), an ideal model for gaining mechanistic insight into the functional specialization of secretory organs. Our findings show that CrebA, a bZIP transcription factor, directly binds genes encoding the core secretory machinery, including protein components of the signal recognition particle and receptor, ER cargo translocators, Cop I and Cop II vesicles, as well as the structural proteins and enzymes of these organelles. CrebA directly binds a subset of SG cargo genes and CrebA binds and boosts expression of Sage, a SG-specific transcription factor essential for cargo expression. To further enhance secretory output, CrebA binds and activates Xbp1 and Tudor-SN. Thus, CrebA directly upregulates the machinery of secretion and additional factors to increase overall secretory capacity in professional secretory cells; concomitant increases in cargo are achieved both directly and indirectly.

New regulators of Drosophila eye development identified from temporal transcriptome changes

In the last larval instar, uncommitted progenitor cells in the Drosophila eye primordium start to adopt individual retinal cell fates, arrest their growth and proliferation, and initiate terminal differentiation into photoreceptor neurons and other retinal cell types. To explore the regulation of these processes, we have performed mRNA-Seq studies of the larval eye and antennal primordial at multiple developmental stages. A total of 10,893 fly genes were expressed during these stages and could be adaptively clustered into gene groups, some of whose expression increases or decreases in parallel with the cessation of proliferation and onset of differentiation. Using in situ hybridization of a sample of 98 genes to verify spatial and temporal expression patterns, we estimate that 534 genes or more are transcriptionally upregulated during retinal differentiation, and 1367 or more downregulated as progenitor cells differentiate. Each group of co-expressed genes is enriched for regulatory motifs recognized by co-expressed transcription factors, suggesting that they represent coherent transcriptional regulatory programs. Using available mutant strains, we describe novel roles for the transcription factors SoxNeuro (SoxN), H6-like homeobox (Hmx), CG10253, without children (woc), Structure specific recognition protein (Ssrp), and multisex combs (mxc).

The Drosophila MLR COMPASS complex is essential for programming cis-regulatory information and maintaining epigenetic memory during development

The MLR COMPASS complex monomethylates H3K4 that serves to epigenetically mark transcriptional enhancers to drive proper gene expression during animal development. Chromatin enrichment analyses of the Drosophila MLR complex reveals dynamic association with promoters and enhancers in embryos with late stage enrichments biased toward both active and poised enhancers. RNAi depletion of the Cmi (also known as Lpt) subunit that contains the chromatin binding PHD finger domains attenuates enhancer functions, but unexpectedly results in inappropriate enhancer activation during stages when hormone responsive enhancers are poised, revealing critical epigenetic roles involved in both the activation and repression of enhancers depending on developmental context. Cmi is necessary for robust H3K4 monomethylation and H3K27 acetylation that mark active enhancers, but not for the chromatin binding of Trr, the MLR methyltransferase. Our data reveal two likely major regulatory modes of MLR function, contributions to enhancer commissioning in early embryogenesis and bookmarking enhancers to enable rapid transcriptional re-activation at subsequent developmental stages.

A distinct class of antioxidant response elements is consistently activated in tumors with NRF2 mutations

NRF2 is a redox-responsive transcription factor that regulates expression of cytoprotective genes via its interaction with DNA sequences known as antioxidant response elements (AREs). NRF2 activity is induced by oxidative stress, but oxidative stress is not the only context in which NRF2 can be activated. Mutations that disrupt the interaction between NRF2 and KEAP1, an inhibitor of NRF2, lead to NRF2 hyperactivation and promote oncogenesis. The mechanisms underlying NRF2's oncogenic properties remain unclear, but likely involve aberrant expression of select NRF2 target genes. We tested this possibility using an integrative genomics approach to get a precise view of the direct NRF2 target genes dysregulated in tumors with NRF2 hyperactivating mutations. This approach revealed a core set of 32 direct NRF2 targets that are consistently upregulated in NRF2 hyperactivated tumors. This set of NRF2 "cancer target genes" includes canonical redox-related NRF2 targets, as well as target genes that have not been previously linked to NRF2 activation. Importantly, NRF2-driven upregulation of this gene set is largely independent of the organ system where the tumor developed. One key distinguishing feature of these NRF2 cancer target genes is that they are regulated by high affinity AREs that fall within genomic regions possessing a ubiquitously permissive chromatin signature. This implies that these NRF2 cancer target genes are responsive to oncogenic NRF2 in most tissues because they lack the regulatory constraints that restrict expression of most other NRF2 target genes. This NRF2 cancer target gene set also serves as a reliable proxy for NRF2 activity, and high NRF2 activity is associated with significant decreases in survival in multiple cancer types. Overall, the pervasive upregulation of these NRF2 cancer targets across multiple cancers, and their association with negative outcomes, suggests that these will be central to dissecting the functional implications of NRF2 hyperactivation in several cancer contexts.

Identification of a functional antioxidant response element at the HIF1A locus

Reactive oxygen species (ROS), which are a byproduct of oxidative metabolism, serve as signaling molecules in a number of physiological settings. However, if their levels are not tightly maintained, excess ROS lead to potentially cytotoxic oxidative stress. Accordingly, several transcriptional regulatory networks have evolved to include components that are highly ROS-responsive. Depending on the context, these regulatory networks can leverage ROS to respond to nutrient conditions, metabolism, or other physiological signals, or to respond to oxidative stress. However, ROS signaling is complex, so regulatory interactions between various ROS-responsive transcription factors are still being mapped out. Here we show that the transcription factor NRF2, a key regulator of the adaptive response to oxidative stress, directly regulates expression of HIF1A, which encodes HIF1α, a key transcriptional regulator of the adaptive response to hypoxia. We used an integrative genomics approach to identify HIF1A as a ROS-responsive transcript and we found an NRF2-bound antioxidant response element (ARE) approximately 30 kilobases upstream of HIF1A. This ARE sequence is deeply conserved, and we verified that it is directly bound and activated by NRF2. In addition, we found that HIF1A is upregulated in breast and bladder tumors with high NRF2 activity. Taken together, our results demonstrate that NRF2 targets a functional ARE at the HIF1A locus, and reveal a direct regulatory connection between two important oxygen responsive transcription factors.

cis-regulatory architecture of a short-range EGFR organizing center in the Drosophila melanogaster leg

We characterized the establishment of an Epidermal Growth Factor Receptor (EGFR) organizing center (EOC) during leg development in Drosophila melanogaster. Initial EGFR activation occurs in the center of leg discs by expression of the EGFR ligand Vn and the EGFR ligand-processing protease Rho, each through single enhancers, vnE and rhoE, that integrate inputs from Wg, Dpp, Dll and Sp1. Deletion of vnE and rhoE eliminates vn and rho expression in the center of the leg imaginal discs, respectively. Animals with deletions of both vnE and rhoE (but not individually) show distal but not medial leg truncations, suggesting that the distal source of EGFR ligands acts at short-range to only specify distal-most fates, and that multiple additional ‘ring’ enhancers are responsible for medial fates. Further, based on the cis-regulatory logic of vnE and rhoE we identified many additional leg enhancers, suggesting that this logic is broadly used by many genes during Drosophila limb development.

The EGFR signaling pathway plays a major role in innumerable developmental processes in all animals and its deregulation leads to different types of cancer, as well as many other developmental diseases in humans. Here we explored the integration of inputs from the Wnt- and TGF-beta signaling pathways and the leg-specifying transcription factors Distal-less and Sp1 at enhancer elements of EGFR ligands. These enhancers trigger a specific EGFR-dependent developmental output in the fly leg that is limited to specifying distal-most fates. Our findings suggest that activation of the EGFR pathway during fly leg development occurs through the activation of multiple EGFR ligand enhancers that are active at different positions along the proximo-distal axis. Similar enhancer elements are likely to control EGFR activation in humans as well. Such DNA elements might be ‘hot spots’ that cause formation of EGFR-dependent tumors if mutations in them occur. Thus, understanding the molecular characteristics of such DNA elements could facilitate the detection and treatment of cancer.

Members of the Candidate Phyla Radiation are functionally differentiated by carbon- and nitrogen-cycling capabilities

BACKGROUND

The Candidate Phyla Radiation (CPR) is a recently described expansion of the tree of life that represents more than 15% of all bacterial diversity and potentially contains over 70 different phyla. Despite this broad phylogenetic variation, these microorganisms appear to feature little functional diversity, with members generally characterized as obligate fermenters. Additionally, much of the data describing CPR phyla has been generated from a limited number of environments, constraining our knowledge of their functional roles and biogeographical distribution. To expand our understanding of subsurface CPR microorganisms, we sampled four separate groundwater wells over 2 years across three Ohio counties.

RESULTS

Samples were analyzed using 16S rRNA gene amplicon and shotgun metagenomic sequencing. Amplicon results indicated that CPR members comprised between 2 and 20% of the microbial communities with relative abundances stable through time in Athens and Greene samples but dynamic in Licking groundwater. Shotgun metagenomic analyses generated 71 putative CPR genomes, representing roughly 32 known phyla and 2 putative new phyla, Candidatus Brownbacteria and Candidatus Hugbacteria. While these genomes largely mirrored metabolic characteristics of known CPR members, some features were previously uncharacterized. For instance, nitrite reductase, encoded by nirK, was found in four of our Parcubacteria genomes and multiple CPR genomes from other studies, indicating a potentially undescribed role for these microorganisms in denitrification. Additionally, glycoside hydrolase (GH) family profiles for our 71 genomes and over 2000 other CPR genomes were analyzed to characterize their carbon-processing potential. Although common trends were present throughout the radiation, differences highlighted potential mechanisms that could allow microorganisms across the CPR to occupy various subsurface niches. For example, members of the Microgenomates superphylum appear to potentially degrade a wider range of carbon substrates than other CPR phyla.

CONCLUSIONS

CPR members are present across a range of environments and often constitute a significant fraction of the microbial population in groundwater systems, particularly. Further sampling of such environments will resolve this portion of the tree of life at finer taxonomic levels, which is essential to solidify functional differences between members that populate this phylogenetically broad region of the tree of life.

A Polymorphic Antioxidant Response Element Links NRF2/sMAF Binding to Enhanced MAPT Expression and Reduced Risk of Parkinsonian Disorders

The NRF2/sMAF protein complex regulates the oxidative stress response by occupying cis-acting enhancers containing an antioxidant response element (ARE). Integrating genome-wide maps of NRF2/sMAF occupancy with disease-susceptibility loci, we discovered eight polymorphic AREs linked to 14 highly ranked disease-risk SNPs in individuals of European ancestry. Among these SNPs was rs242561, located within a regulatory region of the MAPT gene (encoding microtubule-associated protein Tau). It was consistently occupied by NRF2/sMAF in multiple experiments and its strong-binding allele associated with higher mRNA levels in cell lines and human brain tissue. Induction of MAPT transcription by NRF2 was confirmed using a human neuroblastoma cell line and a Nrf2-deficient mouse model. Most importantly, rs242561 displayed complete linkage disequilibrium with a highly protective allele identified in multiple GWASs of progressive supranuclear palsy, Parkinson's disease, and corticobasal degeneration. These observations suggest a potential role for NRF2/sMAF in tauopathies and a possible role for NRF2 pathway activators in disease prevention.

DNA-binding sequence specificity of DUX4

Background

Misexpression of the double homeodomain transcription factor DUX4 results in facioscapulohumeral muscular dystrophy (FSHD). A DNA-binding consensus with two tandem TAAT motifs based on chromatin IP peaks has been discovered; however, the consensus has multiple variations (flavors) of unknown relative activity. In addition, not all peaks have this consensus, and the Pitx1 promoter, the first DUX4 target sequence mooted, has a different TAAT-rich sequence. Furthermore, it is not known whether and to what extent deviations from the consensus affect DNA-binding affinity and transcriptional activation potential.

Results

Here, we take both unbiased and consensus sequence-driven approaches to determine the DNA-binding specificity of DUX4 and its tolerance to mismatches at each site within its consensus sequence. We discover that the best binding and the greatest transcriptional activation are observed when the two TAAT motifs are separated by a C residue. The second TAAT motif in the consensus sequence is actually (T/C)AAT. We find that a T is preferred here. DUX4 has no transcriptional activity on “half-sites”, i.e., those bearing only a single TAAT motif. We further find that DUX4 does not bind to the TAATTA motif in the Pitx1 promoter, that Pitx1 sequences have no competitive band shift activity, and that the Pitx1 sequence is transcriptionally inactive, calling into question PITX1 as a DUX4 target gene. Finally, by multimerizing binding sites, we find that DUX4 transcriptional activation demonstrates tremendous synergy and that at low DNA concentrations, at least two motifs are necessary to detect a transcriptional response.

Conclusions

These studies illuminate the DNA-binding sequence preferences of DUX4.

Electronic supplementary material

The online version of this article (doi:10.1186/s13395-016-0080-z) contains supplementary material, which is available to authorized users.

Beyond antioxidant genes in the ancient Nrf2 regulatory network

Nrf2, a basic leucine zipper transcription factor encoded by the gene NFE2L2, is a master regulator of the transcriptional response to oxidative stress. Nrf2 is structurally and functionally conserved from insects to humans, and it heterodimerizes with the small MAF transcription factors to bind a consensus DNA sequence (the antioxidant response element, or ARE) and regulate gene expression. We have used genome-wide chromatin immunoprecipitation and gene expression data to identify direct Nrf2 target genes in Drosophila and humans. These data have allowed us to construct the deeply conserved ancient Nrf2 regulatory network-target genes that are conserved from Drosophila to human. The ancient network consists of canonical antioxidant genes, as well as genes related to proteasomal pathways and metabolism and a number of less expected genes. We have also used enhancer reporter assays and electrophoretic mobility-shift assays to confirm Nrf2-mediated regulation of ARE activity at a number of these novel target genes. Interestingly, the ancient network also highlights a prominent negative feedback loop; this, combined with the finding that Nrf2-mediated regulatory output is tightly linked to the quality of the ARE it is targeting, suggests that precise regulation of nuclear Nrf2 concentration is necessary to achieve proper quantitative regulation of distinct gene sets. Together, these findings highlight the importance of balance in the Nrf2-ARE pathway and indicate that Nrf2-mediated regulation of xenobiotic metabolism, glucose metabolism, and proteostasis has been central to this pathway since its inception.

Ribbon regulates morphogenesis of the Drosophila embryonic salivary gland through transcriptional activation and repression

Transcription factors affect spatiotemporal patterns of gene expression often regulating multiple aspects of tissue morphogenesis, including cell-type specification, cell proliferation, cell death, cell polarity, cell shape, cell arrangement and cell migration. In this work, we describe a distinct role for Ribbon (Rib) in controlling cell shape/volume increases during elongation of the Drosophila salivary gland (SG). Notably, the morphogenetic changes in rib mutants occurred without effects on general SG cell attributes such as specification, proliferation and apoptosis. Moreover, the changes in cell shape/volume in rib mutants occurred without compromising epithelial-specific morphological attributes such as apicobasal polarity and junctional integrity. To identify the genes regulated by Rib, we performed ChIP-seq analysis in embryos driving expression of GFP-tagged Rib specifically in the SGs. To learn if the Rib binding sites identified in the ChIP-seq analysis were linked to changes in gene expression, we performed microarray analysis comparing RNA samples from age-matched wild-type and rib null embryos. From the superposed ChIP-seq and microarray gene expression data, we identified 60 genomic sites bound by Rib likely to regulate SG-specific gene expression. We confirmed several of the identified Rib targets by qRT-pCR and/or in situ hybridization. Our results indicate that Rib regulates cell growth and tissue shape in the Drosophila salivary gland via a diverse array of targets through both transcriptional activation and repression. Furthermore, our results suggest that autoregulation of rib expression may be a key component of the SG morphogenetic gene network.

Deconvolving the recognition of DNA shape from sequence

Protein-DNA binding is mediated by the recognition of the chemical signatures of the DNA bases and the 3D shape of the DNA molecule. Because DNA shape is a consequence of sequence, it is difficult to dissociate these modes of recognition. Here, we tease them apart in the context of Hox-DNA binding by mutating residues that, in a co-crystal structure, only recognize DNA shape. Complexes made with these mutants lose the preference to bind sequences with specific DNA shape features. Introducing shape-recognizing residues from one Hox protein to another swapped binding specificities in vitro and gene regulation in vivo. Statistical machine learning revealed that the accuracy of binding specificity predictions improves by adding shape features to a model that only depends on sequence, and feature selection identified shape features important for recognition. Thus, shape readout is a direct and independent component of binding site selection by Hox proteins.

Evolution of H3K27me3-marked chromatin is linked to gene expression evolution and to patterns of gene duplication and diversification

Histone modifications are critical for the regulation of gene expression, cell type specification, and differentiation. However, evolutionary patterns of key modifications that regulate gene expression in differentiating organisms have not been examined. Here we mapped the genomic locations of the repressive mark histone 3 lysine 27 trimethylation (H3K27me3) in four species of Drosophila, and compared these patterns to those in C. elegans. We found that patterns of H3K27me3 are highly conserved across species, but conservation is substantially weaker among duplicated genes. We further discovered that retropositions are associated with greater evolutionary changes in H3K27me3 and gene expression than tandem duplications, indicating that local chromatin constraints influence duplicated gene evolution. These changes are also associated with concomitant evolution of gene expression. Our findings reveal the strong conservation of genomic architecture governed by an epigenetic mark across distantly related species and the importance of gene duplication in generating novel H3K27me3 profiles.

Diverse patterns of genomic targeting by transcriptional regulators in Drosophila melanogaster

Annotation of regulatory elements and identification of the transcription-related factors (TRFs) targeting these elements are key steps in understanding how cells interpret their genetic blueprint and their environment during development, and how that process goes awry in the case of disease. One goal of the modENCODE (model organism ENCyclopedia of DNA Elements) Project is to survey a diverse sampling of TRFs, both DNA-binding and non-DNA-binding factors, to provide a framework for the subsequent study of the mechanisms by which transcriptional regulators target the genome. Here we provide an updated map of the Drosophila melanogaster regulatory genome based on the location of 84 TRFs at various stages of development. This regulatory map reveals a variety of genomic targeting patterns, including factors with strong preferences toward proximal promoter binding, factors that target intergenic and intronic DNA, and factors with distinct chromatin state preferences. The data also highlight the stringency of the Polycomb regulatory network, and show association of the Trithorax-like (Trl) protein with hotspots of DNA binding throughout development. Furthermore, the data identify more than 5800 instances in which TRFs target DNA regions with demonstrated enhancer activity. Regions of high TRF co-occupancy are more likely to be associated with open enhancers used across cell types, while lower TRF occupancy regions are associated with complex enhancers that are also regulated at the epigenetic level. Together these data serve as a resource for the research community in the continued effort to dissect transcriptional regulatory mechanisms directing Drosophila development.

Transcription factor networks in Drosophila melanogaster

Specific cellular fates and functions depend on differential gene expression, which occurs primarily at the transcriptional level and is controlled by complex regulatory networks of transcription factors (TFs). TFs act through combinatorial interactions with other TFs, cofactors, and chromatin-remodeling proteins. Here, we define protein-protein interactions using a coaffinity purification/mass spectrometry method and study 459 Drosophila melanogaster transcription-related factors, representing approximately half of the established catalog of TFs. We probe this network in vivo, demonstrating functional interactions for many interacting proteins, and test the predictive value of our data set. Building on these analyses, we combine regulatory network inference models with physical interactions to define an integrated network that connects combinatorial TF protein interactions to the transcriptional regulatory network of the cell. We use this integrated network as a tool to connect the functional network of genetic modifiers related to mastermind, a transcriptional cofactor of the Notch pathway.

Absence of a simple code: how transcription factors read the genome

Transcription factors (TFs) influence cell fate by interpreting the regulatory DNA within a genome. TFs recognize DNA in a specific manner; the mechanisms underlying this specificity have been identified for many TFs based on 3D structures of protein-DNA complexes. More recently, structural views have been complemented with data from high-throughput in vitro and in vivo explorations of the DNA-binding preferences of many TFs. Together, these approaches have greatly expanded our understanding of TF-DNA interactions. However, the mechanisms by which TFs select in vivo binding sites and alter gene expression remain unclear. Recent work has highlighted the many variables that influence TF-DNA binding, while demonstrating that a biophysical understanding of these many factors will be central to understanding TF function.

Yorkie promotes transcription by recruiting a histone methyltransferase complex

Hippo signaling limits organ growth by inhibiting the transcriptional coactivator Yorkie. Despite the key role of Yorkie in both normal and oncogenic growth, the mechanism by which it activates transcription has not been defined. We report that Yorkie binding to chromatin correlates with histone H3K4 methylation and is sufficient to locally increase it. We show that Yorkie can recruit a histone methyltransferase complex through binding between WW domains of Yorkie and PPxY sequence motifs of NcoA6, a subunit of the Trithorax-related (Trr) methyltransferase complex. Cell culture and in vivo assays establish that this recruitment of NcoA6 contributes to Yorkie's ability to activate transcription. Mammalian NcoA6, a subunit of Trr-homologous methyltransferase complexes, can similarly interact with Yorkie's mammalian homolog YAP. Our results implicate direct recruitment of a histone methyltransferase complex as central to transcriptional activation by Yorkie, linking the control of cell proliferation by Hippo signaling to chromatin modification.

Divergent transcriptional regulatory logic at the intersection of tissue growth and developmental patterning

The Yorkie/Yap transcriptional coactivator is a well-known regulator of cellular proliferation in both invertebrates and mammals. As a coactivator, Yorkie (Yki) lacks a DNA binding domain and must partner with sequence-specific DNA binding proteins in the nucleus to regulate gene expression; in Drosophila, the developmental regulators Scalloped (Sd) and Homothorax (Hth) are two such partners. To determine the range of target genes regulated by these three transcription factors, we performed genome-wide chromatin immunoprecipitation experiments for each factor in both the wing and eye-antenna imaginal discs. Strong, tissue-specific binding patterns are observed for Sd and Hth, while Yki binding is remarkably similar across both tissues. Binding events common to the eye and wing are also present for Sd and Hth; these are associated with genes regulating cell proliferation and “housekeeping” functions, and account for the majority of Yki binding. In contrast, tissue-specific binding events for Sd and Hth significantly overlap enhancers that are active in the given tissue, are enriched in Sd and Hth DNA binding sites, respectively, and are associated with genes that are consistent with each factor's previously established tissue-specific functions. Tissue-specific binding events are also significantly associated with Polycomb targeted chromatin domains. To provide mechanistic insights into tissue-specific regulation, we identify and characterize eye and wing enhancers of the Yki-targeted bantam microRNA gene and demonstrate that they are dependent on direct binding by Hth and Sd, respectively. Overall these results suggest that both Sd and Hth use distinct strategies – one shared between tissues and associated with Yki, the other tissue-specific, generally Yki-independent and associated with developmental patterning – to regulate distinct gene sets during development.

The Hippo tumor suppressor pathway controls proliferation in a tissue-nonspecific fashion in Drosophila epithelial progenitor tissues via the transcriptional coactivator Yorkie (Yki). However, despite the tissue-nonspecific role that Yki plays in tissue growth, the transcription factors that recruit Yki to DNA, most notably Scalloped (Sd) and Homothorax (Hth), are important regulators of developmental patterning with many tissue-specific functions. Thus, these three transcriptional regulators – Yki, Sd, and Hth – provide a model for exploring the properties of protein-DNA interactions that regulate both tissue-shared and tissue-specific functions. With this goal in mind, we identified the positions in the fly genome that are bound by Yki, Sd, and Hth in the progenitors of the wing and eye-antenna structures of the fly. These data not only provide a global view of the Yki gene regulatory network, they reveal an unusual amount of tissue specificity in the genomic regions targeted by Sd and Hth, but not Yki. The data also reveal that tissue-specific binding is very likely to overlap tissue-specific enhancer regions, provide important clues for how tissue-specific Sd and Hth binding occurs, and support the idea that gene regulatory networks are plastic, with spatial differences in binding significantly impacting network structures.

Interpreting the regulatory genome: the genomics of transcription factor function in Drosophila melanogaster

Researchers have now had access to the fully sequenced Drosophila melanogaster genome for over a decade, and the sequenced genomes of 11 additional Drosophila species have been available for almost 5 years, with more species' genomes becoming available every year [Adams MD, Celniker SE, Holt RA, et al. The genome sequence of Drosophila melanogaster. Science 2000;287:2185-95; Clark AG, Eisen MB, Smith DR, et al. Evolution of genes and genomes on the Drosophila phylogeny. Nature 2007;450:203-18]. Although the best studied of the D. melanogaster transcription factors (TFs) were cloned before sequencing of the genome, the availability of sequence data promised to transform our understanding of TFs and gene regulatory networks. Sequenced genomes have allowed researchers to generate tools for high-throughput characterization of gene expression levels, genome-wide TF localization and analyses of evolutionary constraints on DNA elements across multiple species. With an estimated 700 DNA-binding proteins in the Drosophila genome, it will be many years before each potential sequence-specific TF is studied in detail, yet the last decade of functional genomics research has already impacted our view of gene regulatory networks and TF DNA recognition.

Genome-wide association of Yorkie with chromatin and chromatin-remodeling complexes

The Hippo pathway regulates growth through the transcriptional coactivator Yorkie, but how Yorkie promotes transcription remains poorly understood. We address this by characterizing Yorkie's association with chromatin and by identifying nuclear partners that effect transcriptional activation. Coimmunoprecipitation and mass spectrometry identify GAGA factor (GAF), the Brahma complex, and the Mediator complex as Yorkie-associated nuclear protein complexes. All three are required for Yorkie's transcriptional activation of downstream genes, and GAF and the Brahma complex subunit Moira interact directly with Yorkie. Genome-wide chromatin-binding experiments identify thousands of Yorkie sites, most of which are associated with elevated transcription, based on genome-wide analysis of messenger RNA and histone H3K4Me3 modification. Chromatin binding also supports extensive functional overlap between Yorkie and GAF. Our studies suggest a widespread role for Yorkie as a regulator of transcription and identify recruitment of the chromatin-modifying GAF protein and BRM complex as a molecular mechanism for transcriptional activation by Yorkie.

ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia

Chromatin immunoprecipitation (ChIP) followed by high-throughput DNA sequencing (ChIP-seq) has become a valuable and widely used approach for mapping the genomic location of transcription-factor binding and histone modifications in living cells. Despite its widespread use, there are considerable differences in how these experiments are conducted, how the results are scored and evaluated for quality, and how the data and metadata are archived for public use. These practices affect the quality and utility of any global ChIP experiment. Through our experience in performing ChIP-seq experiments, the ENCODE and modENCODE consortia have developed a set of working standards and guidelines for ChIP experiments that are updated routinely. The current guidelines address antibody validation, experimental replication, sequencing depth, data and metadata reporting, and data quality assessment. We discuss how ChIP quality, assessed in these ways, affects different uses of ChIP-seq data. All data sets used in the analysis have been deposited for public viewing and downloading at the ENCODE (http://encodeproject.org/ENCODE/) and modENCODE (http://www.modencode.org/) portals.

Disentangling the many layers of eukaryotic transcriptional regulation

Regulation of gene expression in eukaryotes is an extremely complex process. In this review, we break down several critical steps, emphasizing new data and techniques that have expanded current gene regulatory models. We begin at the level of DNA sequence where cis-regulatory modules (CRMs) provide important regulatory information in the form of transcription factor (TF) binding sites. In this respect, CRMs function as instructional platforms for the assembly of gene regulatory complexes. We discuss multiple mechanisms controlling complex assembly, including cooperative DNA binding, combinatorial codes, and CRM architecture. The second section of this review places CRM assembly in the context of nucleosomes and condensed chromatin. We discuss how DNA accessibility and histone modifications contribute to TF function. Lastly, new advances in chromosomal mapping techniques have provided increased understanding of intra- and interchromosomal interactions. We discuss how these topological maps influence gene regulatory models.

Control of mitochondrial structure and function by the Yorkie/YAP oncogenic pathway

Mitochondrial structure and function are highly dynamic, but the potential roles for cell signaling pathways in influencing these properties are not fully understood. Reduced mitochondrial function has been shown to cause cell cycle arrest, and a direct role of signaling pathways in controlling mitochondrial function during development and disease is an active area of investigation. Here, we show that the conserved Yorkie/YAP signaling pathway implicated in the control of organ size also functions in the regulation of mitochondria in Drosophila as well as human cells. In Drosophila, activation of Yorkie causes direct transcriptional up-regulation of genes that regulate mitochondrial fusion, such as opa1-like (opa1) and mitochondria assembly regulatory factor (Marf), and results in fused mitochondria with dramatic reduction in reactive oxygen species (ROS) levels. When mitochondrial fusion is genetically attenuated, the Yorkie-induced cell proliferation and tissue overgrowth are significantly suppressed. The function of Yorkie is conserved across evolution, as activation of YAP2 in human cell lines causes increased mitochondrial fusion. Thus, mitochondrial fusion is an essential and direct target of the Yorkie/YAP pathway in the regulation of organ size control during development and could play a similar role in the genesis of cancer.

Cofactor binding evokes latent differences in DNA binding specificity between Hox proteins

Members of transcription factor families typically have similar DNA binding specificities yet execute unique functions in vivo. Transcription factors often bind DNA as multiprotein complexes, raising the possibility that complex formation might modify their DNA binding specificities. To test this hypothesis, we developed an experimental and computational platform, SELEX-seq, that can be used to determine the relative affinities to any DNA sequence for any transcription factor complex. Applying this method to all eight Drosophila Hox proteins, we show that they obtain novel recognition properties when they bind DNA with the dimeric cofactor Extradenticle-Homothorax (Exd). Exd-Hox specificities group into three main classes that obey Hox gene collinearity rules and DNA structure predictions suggest that anterior and posterior Hox proteins prefer DNA sequences with distinct minor groove topographies. Together, these data suggest that emergent DNA recognition properties revealed by interactions with cofactors contribute to transcription factor specificities in vivo.

Genome-wide tissue-specific occupancy of the Hox protein Ultrabithorax and Hox cofactor Homothorax in Drosophila

The Hox genes are responsible for generating morphological diversity along the anterior-posterior axis during animal development. The Drosophila Hox gene Ultrabithorax (Ubx), for example, is required for specifying the identity of the third thoracic (T3) segment of the adult, which includes the dorsal haltere, an appendage required for flight, and the ventral T3 leg. Ubx mutants show homeotic transformations of the T3 leg towards the identity of the T2 leg and the haltere towards the wing. All Hox genes, including Ubx, encode homeodomain containing transcription factors, raising the question of what target genes Ubx regulates to generate these adult structures. To address this question, we carried out whole genome ChIP-chip studies to identify all of the Ubx bound regions in the haltere and T3 leg imaginal discs, which are the precursors to these adult structures. In addition, we used ChIP-chip to identify the sites bound by the Hox cofactor, Homothorax (Hth). In contrast to previous ChIP-chip studies carried out in Drosophila embryos, these binding studies reveal that there is a remarkable amount of tissue- and transcription factor-specific binding. Analyses of the putative target genes bound and regulated by these factors suggest that Ubx regulates many downstream transcription factors and developmental pathways in the haltere and T3 leg. Finally, we discovered additional DNA sequence motifs that in some cases are specific for individual data sets, arguing that Ubx and/or Hth work together with many regionally expressed transcription factors to execute their functions. Together, these data provide the first whole-genome analysis of the binding sites and target genes regulated by Ubx to specify the morphologies of the adult T3 segment of the fly.

Transcription factor choice in the Hippo signaling pathway: homothorax and yorkie regulation of the microRNA bantam in the progenitor domain of the Drosophila eye imaginal disc

The accurate control of cell proliferation and survival is critical for animal development. The Hippo tumor suppressor pathway regulates both of these parameters by controlling the nuclear availability of the transcriptional coactivator Yorkie (Yki), which regulates downstream target genes together with Scalloped (Sd), a DNA-binding protein. Here we provide evidence that Yki can also regulate target genes in conjunction with Homothorax (Hth) and Teashirt (Tsh), two DNA-binding transcription factors expressed in the uncommitted progenitor cells of the Drosophila eye imaginal disc. Clonal analyses demonstrate that Hth and Tsh promote cell proliferation and protect eye progenitor cells from apoptosis. Genetic epistasis experiments suggest that Hth and Tsh execute these functions with Yki, in part by up-regulating the microRNA bantam. A physical interaction between Hth and Yki can be detected in cell culture, and we show that Hth and Yki are bound to a DNA sequence approximately 14 kb upstream of the bantam hairpin in eye imaginal disc cells, arguing that this regulation is direct. These data suggest that the Hippo pathway uses different DNA-binding transcription factors depending on the cellular context. In the eye disc, Hth and Tsh provide spatial information to this pathway, promoting cell proliferation and survival in the progenitor domain.

Species differences and effects of soft coral extracts from Sinnularia maximus on the expression of cytochrome P4501A and 2N in butterflyfishes (Chaetodon spp.)

Cytochrome P450 (CYP) has been shown to confer resistance in numerous terrestrial insects that consume potentially toxic secondary metabolites in plants, but fewer studies have examined the role of critical biotransformation enzymes in allowing marine organisms to consume chemically defended foods. This study examined the expression of CYP1A and CYP2N mRNAs in several butterflyfish species, which can feed on numerous chemically defended soft and hard corals. In addition, the effect of an extract from a soft coral (Sinnularia maxima) on expression of hepatic CYP1A and CYP2 mRNAs was also examined. Fish were fed extracts on days 1, 3 and 5, and expression was examined on day 5. Phylogenetic analyses of the CYP1A cDNA from 12 species of butterflyfish (DNA, amino acid) indicate well-separated groupings according to their feeding strategies. The non-coralline feeding Chaetodon xanthurus exhibited a 7-fold higher basal expression of CYP2N8 relative to the other species studied. Although induction of CYP2N7 expression was observed in C. punctatofasciatus, CYP1A and CYP2N was largely unaffected or diminished by extract treatment in the other species of butterflyfish. These results indicated groupings of feeding strategy with CYP1A phylogeny in Chaetodon, but generally unaltered expression of CYP1A and CYP2N following dietary treatment with an extract from a chemically defended soft coral suggesting an inconclusive role of these isoforms in the detoxification of chemicals in these extracts.

A lexicon for homeodomain-DNA recognition

Decoding the cis-regulatory logic of eukaryotic genomes requires knowledge of the DNA-binding specificities of all transcription factors. New work (Berger et al., 2008; Noyes et al., 2008) provides individual specificities for nearly all Drosophila and mouse homeodomains, key DNA-binding domains in many transcription factors. The data underscore the complexity of determining target specificities in vivo.

Chemical inhibition of the TFIIH-associated kinase Cdk7/Kin28 does not impair global mRNA synthesis

The process of gene transcription requires the recruitment of a hypophosphorylated form of RNA polymerase II (Pol II) to a gene promoter. The TFIIH-associated kinase Cdk7/Kin28 hyperphosphorylates the promoter-bound polymerase; this event is thought to play a crucial role in transcription initiation and promoter clearance. Studies using temperature-sensitive mutants of Kin28 have provided the most compelling evidence for an essential role of its kinase activity in global mRNA synthesis. In contrast, using a small molecule inhibitor that specifically inhibits Kin28 in vivo, we find that the kinase activity is not essential for global transcription. Unlike the temperature-sensitive alleles, the small-molecule inhibitor does not perturb protein-protein interactions nor does it provoke the disassociation of TFIIH from gene promoters. These results lead us to conclude that other functions of TFIIH, rather than the kinase activity, are critical for global gene transcription.

Aplysina red band syndrome: a new threat to Caribbean sponges

A substantial and increasing number of reports have documented dramatic changes and continuing declines in Caribbean coral reef communities over the past 2 decades. To date, the majority of disease reports have focused on scleractinian corals, whereas sponge diseases have been less frequently documented. In this study, we describe Aplysina red band syndrome (ARBS) affecting Caribbean rope sponges of the genus Aplysina observed on shallow reefs in the Bahamas. Visible signs of disease presence included 1 or more rust-colored leading edges, with or without a trailing area of necrotic tissue, such that the lesion forms a contiguous band around part or all of the sponge branch. Microscopic examination of the leading edge of the disease margin indicated that a cyanobacterium was consistently responsible for the coloration. Although the presence of this distinctive coloration was used to characterize the diseased state, it is not yet known whether this cyanobacterium is directly responsible for disease causation. The prevalence of ARBS declined significantly from July to October 2004 before increasing above July levels in January 2005. Transmission studies in the laboratory demonstrated that contact with the leading edge of an active lesion was sufficient to spread ARBS to a previously healthy sponge, suggesting that the etiologic agent, currently undescribed, is contagious. Studies to elucidate the etiologic agent of ARBS are ongoing. Sponges are an essential component of coral reef communities and emerging sponge diseases clearly have the potential to impact benthic community structure on coral reefs.

Hormone replacement therapy and colorectal adenoma recurrence among women in the Polyp Prevention Trial

BACKGROUND

Epidemiologic studies have suggested that estrogen may protect against the development of colorectal cancers and adenomatous polyps. We conducted a prospective study to evaluate the association between hormone replacement therapy (HRT) and adenoma recurrence among perimenopausal and postmenopausal women participating in the Polyp Prevention Trial, a randomized dietary intervention study of individuals with colorectal adenomas.

METHODS

We used a questionnaire and interviews to collect detailed information, at baseline and at each of four annual study visits, from 620 women regarding hormone use, menopausal status, diet, alcohol consumption, and other risk factors. Adenoma recurrence was ascertained by complete colonoscopy at baseline and after 1 and 4 years. Logistic regression models were used to evaluate the association between hormone use and adenoma recurrence after adjusting for intervention group and for age and body mass index at baseline. All statistical tests were two-sided.

RESULTS

Adenomas recurred in 200 women. There was no overall association between adenoma recurrence and either overall hormone use (odds ratio [OR] = 1.01; 95% confidence interval [CI] = 0.70 to 1.45), combined estrogen and progestin use (OR = 0.94; 95% CI = 0.57 to 1.56), or unopposed estrogen use (OR = 1.04; 95% CI = 0.68 to 1.59). HRT use was associated with a reduction in risk for recurrence of distal adenomas (OR = 0.56; 95% CI = 0.32 to 1.00) and a statistically nonsignificant increase in risk for recurrence of proximal adenomas (OR = 1.39; 95% CI = 0.85 to 2.26). We observed a statistically significant interaction between the HRT-adenoma recurrence association and age (P =.02). HRT was associated with a 40% reduced risk of adenoma recurrence among women older than 62 years (OR = 0.58; 95% CI = 0.35 to 0.97) but with an increased risk among women younger than 62 years (OR = 1.99; 95% CI = 1.11 to 3.55).

CONCLUSIONS

HRT was not associated with a reduced risk for overall adenoma recurrence in this trial cohort, although there was a suggestion of an age interaction. The effect of age on the association needs to be confirmed in other adenoma recurrence trials.

Implementation of a 4-y, high-fiber, high-fruit-and-vegetable, low-fat dietary intervention: results of dietary changes in the Polyp Prevention Trial

BACKGROUND

The Polyp Prevention Trial (PPT) was a multicenter randomized clinical trial designed to determine the effects of a high-fiber (4.30 g/MJ), high-fruit-and-vegetable (0.84 servings/MJ), low-fat (20% of energy from fat) diet on the recurrence of adenomatous polyps in the large bowel.

OBJECTIVE

Our goal was to determine whether the PPT intervention plan could effect change in 3 dietary goals and to examine the intervention's effect on the intake of other food groups and nutrients.

DESIGN

Participants with large-bowel adenomatous polyps diagnosed in the past 6 mo were randomly assigned to either the intervention (n = 1037) or the control (n = 1042) group and remained in the trial for 4 y. Three dietary assessment instruments were used to measure dietary change: food-frequency questionnaires (in 100% of the sample), 4-d food records (in a 20% random cohort), and 24-h dietary recalls (in a 10% random sample).

RESULTS

Intervention participants made and sustained significant changes in all PPT goals as measured by the dietary assessment instruments; the control participants' intakes remained essentially the same throughout the trial. The absolute differences between the intervention and control groups over the 4-y period were 9.7% of energy from fat (95% CI: 9.0%, 10.3%), 1.65 g dietary fiber/MJ (95% CI: 1.53, 1.74), and 0.27 servings of fruit and vegetables/MJ (95% CI: 0.25, 0.29). Intervention participants also reported significant changes in the intake of other nutrients and food groups. The intervention group also had significantly higher serum carotenoid concentrations and lower body weights than did the control group.

CONCLUSION

Motivated, free-living individuals, given appropriate support, can make and sustain major dietary changes over a 4-y period.

The efficacy of functional foot orthoses in the control of pain in ankle joint disintegration in hemophilia

Sixteen subjects with hemophilia A of levels 1-5 stage of joint damage were tested over a 6-week period to evaluate the efficacy of functional foot orthoses. The level of ankle bleeds (hemarthrosis) before and after the intervention with functional foot orthoses was determined by evaluating pain, disability, and activity levels. All subjects reported a significant reduction of ankle bleeds coinciding with the intervention of functional foot orthoses. The use of a foot-pain disability measure clearly showed significant reduction in the level of pain experienced by the subjects and in their overall index score. However, the disability and activity index scores showed no significant improvement after the intervention with orthoses. This finding would support the use of functional foot orthoses to treat patients with hemophilia A, as significant reduction in pain levels appears to greatly improve the lives of the patients.

Competition-Mediated Antibiotic Induction in the Marine Bacterium Streptomyces tenjimariensis

Microbial competition for limiting natural resources within a community is thought to be the selective force that promotes biosynthesis of antimicrobial compounds. The marine bacterium Streptomyces tenjimariensis produces the antibiotics istamycin A and B under select laboratory culture conditions; presumably these compounds serve an ecological role under natural conditions. Here we report results of a novel marine microbial competition experiment that examined the impact of co-culture of marine bacteria on istamycin production by S. tenjimariensis. Twelve of the 53 bacterial species tested (i.e., 22.6%) induced Istamycin production; this antibiotic also inhibited growth of the competitor colonies. These results suggest that marine bacterial metabolites serve an ecological role in countering competitive species.