Shane T. Jensen

The Spo0A regulon of Bacillus subtilis

The master regulator for entry into sporulation in Bacillus subtilis is the DNA‐binding protein Spo0A, which has been found to influence, directly or indirectly, the expression of over 500 genes during the early stages of development. To search on a genome‐wide basis for genes under the direct control of Spo0A, we used chromatin immunoprecipitation in combination with gene microarray analysis to identify regions of the chromosome at which an activated form of Spo0A binds in vivo. This information in combination with transcriptional profiling using gene microarrays, gel electrophoretic mobility shift assays, using the DNA‐binding domain of Spo0A, and bioinformatics enabled us to assign 103 genes to the Spo0A regulon in addition to 18 previously known members. Thus, in total, 121 genes, which are organized as 30 single‐gene units and 24 operons, are likely to be under the direct control of Spo0A. Forty of these genes are under the positive control of Spo0A, and 81 are under its negative control. Among newly identified members of the regulon with transcription that was stimulated by Spo0A are genes for metabolic enzymes and genes for efflux pumps. Among  members  with  transcription  that  was in‐hibited by Spo0A are genes encoding components of the DNA replication machinery and genes that govern flagellum biosynthesis and chemotaxis. Also in‐cluded in the regulon are many (25) genes with products that are direct or indirect regulators of gene transcription. Spo0A is a master regulator for sporulation, but many of its effects on the global pattern of gene transcription are likely to be mediated indirectly by regulatory genes under its control.

The Program of Gene Transcription for a Single Differentiating Cell Type during Sporulation in Bacillus subtilis

Asymmetric division during sporulation by Bacillus subtilis generates a mother cell that undergoes a 5-h program of differentiation. The program is governed by a hierarchical cascade consisting of the transcription factors: σE, σK, GerE, GerR, and SpoIIID. The program consists of the activation and repression of 383 genes. The σE factor turns on 262 genes, including those for GerR and SpoIIID. These DNA-binding proteins downregulate almost half of the genes in the σE regulon. In addition, SpoIIID turns on ten genes, including genes involved in the appearance of σK . Next, σK activates 75 additional genes, including that for GerE. This DNA-binding protein, in turn, represses half of the genes that had been activated by σK while switching on a final set of 36 genes. Evidence is presented that repression and activation contribute to proper morphogenesis. The program of gene expression is driven forward by its hierarchical organization and by the repressive effects of the DNA-binding proteins. The logic of the program is that of a linked series of feed-forward loops, which generate successive pulses of gene transcription. Similar regulatory circuits could be a common feature of other systems of cellular differentiation.

The σE regulon and the identification of additional sporulation genes in Bacillus subtilis

We report the identification and characterization on a genome-wide basis of genes under the control of the developmental transcription factor sigma(E) in Bacillus subtilis. The sigma(E) factor governs gene expression in the larger of the two cellular compartments (the mother cell) created by polar division during the developmental process of sporulation. Using transcriptional profiling and bioinformatics we show that 253 genes (organized in 157 operons) appear to be controlled by sigma(E). Among these, 181 genes (organized in 121 operons) had not been previously described as members of this regulon. Promoters for many of the newly identified genes were located by transcription start site mapping. To assess the role of these genes in sporulation, we created null mutations in 98 of the newly identified genes and operons. Of the resulting mutants, 12 (in prkA, ybaN, yhbH, ykvV, ylbJ, ypjB, yqfC, yqfD, ytrH, ytrI, ytvI and yunB) exhibited defects in spore formation. In addition, subcellular localization studies were carried out using in-frame fusions of several of the genes to the coding sequence for GFP. A majority of the fusion proteins localized either to the membrane surrounding the developing spore or to specific layers of the spore coat, although some fusions showed a uniform distribution in the mother cell cytoplasm. Finally, we used comparative genomics to determine that 46 of the sigma(E)-controlled genes in B.subtilis were present in all of the Gram-positive endospore-forming bacteria whose genome has been sequenced, but absent from the genome of the closely related but not endospore-forming bacterium Listeria monocytogenes, thereby defining a core of conserved sporulation genes of probable common ancestral origin. Our findings set the stage for a comprehensive understanding of the contribution of a cell-specific transcription factor to development and morphogenesis.

LEAFY Target Genes Reveal Floral Regulatory Logic, cis Motifs, and a Link to Biotic Stimulus Response

The transition from vegetative growth to flower formation is critical for the survival of flowering plants. The plant-specific transcription factor LEAFY (LFY) has central, evolutionarily conserved roles in this process, both in the formation of the first flower and later in floral patterning. We performed genome-wide binding and expression studies to elucidate the molecular mechanisms by which LFY executes these roles. Our study reveals that LFY directs an elaborate regulatory network in control of floral homeotic gene expression. LFY also controls the expression of genes that regulate the response to external stimuli in Arabidopsis. Thus, our findings support a key role for LFY in the coordination of reproductive stage development and disease response programs in plants that may ensure optimal allocation of plant resources for reproductive fitness. Finally, motif analyses reveal a possible mechanism for stage-specific LFY recruitment and suggest a role for LFY in overcoming polycomb repression.

BioOptimizer: a Bayesian scoring function approach to motif discovery

MOTIVATION Transcription factors (TFs) bind directly to short segments on the genome, often within hundreds to thousands of base pairs upstream of gene transcription start sites, to regulate gene expression. The experimental determination of TFs binding sites is expensive and time-consuming. Many motif-finding programs have been developed, but no program is clearly superior in all situations. Practitioners often find it difficult to judge which of the motifs predicted by these algorithms are more likely to be biologically relevant. RESULTS We derive a comprehensive scoring function based on a full Bayesian model that can handle unknown site abundance, unknown motif width and two-block motifs with variable-length gaps. An algorithm called BioOptimizer is proposed to optimize this scoring function so as to reduce noise in the motif signal found by any motif-finding program. The accuracy of BioOptimizer, which can be used in conjunction with several existing programs, is shown to be superior to using any of these motif-finding programs alone when evaluated by both simulation studies and application to sets of co-regulated genes in bacteria. In addition, this scoring function formulation enables us to compare objectively different predicted motifs and select the optimal ones, effectively combining the strengths of existing programs. AVAILABILITY BioOptimizer is available for download at www.fas.harvard.edu/~junliu/BioOptimizer/

Potential Mechanisms for Cancer Resistance in Elephants and Comparative Cellular Response to DNA Damage in Humans

IMPORTANCE Evolutionary medicine may provide insights into human physiology and pathophysiology, including tumor biology. OBJECTIVE To identify mechanisms for cancer resistance in elephants and compare cellular response to DNA damage among elephants, healthy human controls, and cancer-prone patients with Li-Fraumeni syndrome (LFS). DESIGN, SETTING, AND PARTICIPANTS A comprehensive survey of necropsy data was performed across 36 mammalian species to validate cancer resistance in large and long-lived organisms, including elephants (n = 644). The African and Asian elephant genomes were analyzed for potential mechanisms of cancer resistance. Peripheral blood lymphocytes from elephants, healthy human controls, and patients with LFS were tested in vitro in the laboratory for DNA damage response. The study included African and Asian elephants (n = 8), patients with LFS (n = 10), and age-matched human controls (n = 11). Human samples were collected at the University of Utah between June 2014 and July 2015. EXPOSURES Ionizing radiation and doxorubicin. MAIN OUTCOMES AND MEASURES Cancer mortality across species was calculated and compared by body size and life span. The elephant genome was investigated for alterations in cancer-related genes. DNA repair and apoptosis were compared in elephant vs human peripheral blood lymphocytes. RESULTS Across mammals, cancer mortality did not increase with body size and/or maximum life span (eg, for rock hyrax, 1% [95% CI, 0%-5%]; African wild dog, 8% [95% CI, 0%-16%]; lion, 2% [95% CI, 0%-7%]). Despite their large body size and long life span, elephants remain cancer resistant, with an estimated cancer mortality of 4.81% (95% CI, 3.14%-6.49%), compared with humans, who have 11% to 25% cancer mortality. While humans have 1 copy (2 alleles) of TP53, African elephants have at least 20 copies (40 alleles), including 19 retrogenes (38 alleles) with evidence of transcriptional activity measured by reverse transcription polymerase chain reaction. In response to DNA damage, elephant lymphocytes underwent p53-mediated apoptosis at higher rates than human lymphocytes proportional to TP53 status (ionizing radiation exposure: patients with LFS, 2.71% [95% CI, 1.93%-3.48%] vs human controls, 7.17% [95% CI, 5.91%-8.44%] vs elephants, 14.64% [95% CI, 10.91%-18.37%]; P < .001; doxorubicin exposure: human controls, 8.10% [95% CI, 6.55%-9.66%] vs elephants, 24.77% [95% CI, 23.0%-26.53%]; P < .001). CONCLUSIONS AND RELEVANCE Compared with other mammalian species, elephants appeared to have a lower-than-expected rate of cancer, potentially related to multiple copies of TP53. Compared with human cells, elephant cells demonstrated increased apoptotic response following DNA damage. These findings, if replicated, could represent an evolutionary-based approach for understanding mechanisms related to cancer suppression.

Computational Discovery of Gene Regulatory Binding Motifs: A Bayesian Perspective

The Bayesian approach together with Markov chain Monte Carlo techniques has provided an attractive solution to many important bioinformatics problems such as multiple sequence alignment, microarray analysis and the discovery of gene regulatory binding motifs. The employment of such methods and, more broadly, explicit statistical modeling, has revolutionized the field of computational biology. After reviewing several heuristicsbased computational methods, this article presents a systematic account of Bayesian formulations and solutions to the motif discovery problem. Generalizations are made to further enhance the Bayesian approach. Motivated by the need of a speedy algorithm, we also provide a perspective of the problem from the viewpoint of optimizing a scoring function. We observe that scoring functions resulting from proper posterior distributions, or approximations to such distributions, showed the best performance and can be used to improve upon existing motif-finding programs. Simulation analyses and a real-data example are used to support our observation.

Mapping structural landmarks, ligand binding sites, and missense mutations to the collagen IV heterotrimers predicts major functional domains, novel interactions, and variation in phenotypes in inherited diseases affecting basement membranes.

Collagen IV is the major protein found in basement membranes. It comprises three heterotrimers (α1α1α2, α3α4α5, and α5α5α6) that form distinct networks, and are responsible for membrane strength and integrity. We constructed linear maps of the collagen IV heterotrimers (“interactomes”) that indicated major structural landmarks, known and predicted ligand‐binding sites, and missense mutations, in order to identify functional and disease‐associated domains, potential interactions between ligands, and genotype–phenotype relationships. The maps documented more than 30 known ligand‐binding sites as well as motifs for integrins, heparin, von Willebrand factor (VWF), decorin, and bone morphogenetic protein (BMP). They predicted functional domains for angiogenesis and haemostasis, and disease domains for autoimmunity, tumor growth and inhibition, infection, and glycation. Cooperative ligand interactions were indicated by binding site proximity, for example, between integrins, matrix metalloproteinases, and heparin. The maps indicated that mutations affecting major ligand‐binding sites, for example, for Von Hippel Lindau (VHL) protein in the α1 chain or integrins in the α5 chain, resulted in distinctive phenotypes (Hereditary Angiopathy, Nephropathy, Aneurysms, and muscle Cramps [HANAC] syndrome, and early‐onset Alport syndrome, respectively). These maps further our understanding of basement membrane biology and disease, and suggest novel membrane interactions, functions, and therapeutic targets. Hum Mutat 32:127–143, 2011.

Clustering of genes into regulons using integrated modeling-COGRIM

We present a Bayesian hierarchical model and Gibbs Sampling implementation that integrates gene expression, ChIP binding, and transcription factor motif data in a principled and robust fashion. COGRIM was applied to both unicellular and mammalian organisms under different scenarios of available data. In these applications, we demonstrate the ability to predict gene-transcription factor interactions with reduced numbers of false-positive findings and to make predictions beyond what is obtained when single types of data are considered.

Single-Cell Analysis of Proxy Reporter Allele-Marked Epithelial Cells Establishes Intestinal Stem Cell Hierarchy

Summary The recent development of targeted murine reporter alleles as proxies for intestinal stem cell activity has led to significant advances in our understanding of somatic stem cell hierarchies and dynamics. Analysis of these reporters has led to a model in which an indispensable reserve stem cell at the top of the hierarchy (marked by Bmi1 and Hopx reporters) gives rise to active intestinal stem cells (marked by an Lgr5 reporter). Despite these advances, controversy exists regarding the specificity and fidelity with which these alleles distinguish intestinal stem cell populations. Here, we undertake a comprehensive comparison of widely used proxy reporters including both CreERT2 and EGFP cassettes targeted to the Lgr5, Bmi1, and Hopx loci. Single-cell transcriptional profiling of these populations and their progeny reveals that reserve and active intestinal stem cells are molecularly and functionally distinct, supporting a two-stem-cell model for intestinal self-renewal.