Following the 'tracks': Tramtrack69 regulates epithelial tube expansion in the Drosophila ovary through Paxillin, Dynamin, and the homeobox protein Mirror.

Epithelial tubes are the infrastructure for organs and tissues, and tube morphogenesis requires precise orchestration of cell signaling, shape, migration, and adhesion. Follicle cells in the Drosophila ovary form a pair of epithelial tubes whose lumens act as molds for the eggshell respiratory filaments, or dorsal appendages (DAs). DA formation is a robust and accessible model for studying the patterning, formation, and expansion of epithelial tubes. Tramtrack69 (TTK69), a transcription factor that exhibits a variable embryonic DNA-binding preference, controls DA lumen volume and shape by promoting tube expansion; the tramtrack mutation twin peaks (ttk(twk)) reduces TTK69 levels late in oogenesis, inhibiting this expansion. Microarray analysis of wild-type and ttk(twk) ovaries, followed by in situ hybridization and RNAi of candidate genes, identified the Phospholipase B-like protein Lamina ancestor (LAMA), the scaffold protein Paxillin, the endocytotic regulator Shibire (Dynamin), and the homeodomain transcription factor Mirror, as TTK69 effectors of DA-tube expansion. These genes displayed enriched expression in DA-tube cells, except lama, which was expressed in all follicle cells. All four genes showed reduced expression in ttk(twk) mutants and exhibited RNAi phenotypes that were enhanced in a ttk(twk)/+ background, indicating ttk(twk) genetic interactions. Although previous studies show that Mirror patterns the follicular epithelium prior to DA tubulogenesis, we show that Mirror has an independent, novel role in tube expansion, involving positive regulation of Paxillin. Thus, characterization of ttk(twk)-differentially expressed genes expands the network of TTK69 effectors, identifies novel epithelial tube-expansion regulators, and significantly advances our understanding of this vital developmental process.

Identification of the Treponema pallidum subsp. pallidum TP0092 (RpoE) regulon and its implications for pathogen persistence in the host and syphilis pathogenesis.

Bacteria often respond to harmful environmental stimuli with the induction of extracytoplasmic function (ECF) sigma (σ) factors that in turn direct RNA polymerase to transcribe specific groups of response genes (or regulons) to minimize cellular damage and favor adaptation to the changed extracellular milieu. In Treponema pallidum subsp. pallidum, the agent of syphilis, the TP0092 gene is predicted to code for the pathogen's only annotated ECF σ factor, homologous to RpoE, known in Escherichia coli to control a key transduction pathway for maintenance of envelope homeostasis in response to external stress and cell growth. Here we have shown that TP0092 is highly transcribed during experimental syphilis. Furthermore, TP0092 transcription levels significantly increase as infection progresses toward immune clearance of the pathogen, suggesting a role for TP0092 in helping T. pallidum respond to harmful stimuli in the host environment. To investigate this hypothesis, we determined the TP0092 regulon at two different time points during infection using chromatin immunoprecipitation followed by high-throughput sequencing. A total of 22 chromosomal regions, all containing putative TP0092-binding sites and corresponding to as many T. pallidum genes, were identified. Noteworthy among them are the genes encoding desulfoferrodoxin and thioredoxin, involved in detoxification of reactive oxygen species (ROS). Because T. pallidum does not possess other enzymes for ROS detoxification, such as superoxide dismutase, catalase, or glutathione peroxidase, our results suggest that the TP0092 regulon is important in protecting the syphilis spirochete from damage caused by ROS produced at the site of infection during the inflammatory response.

Comparative assessment of methods for aligning multiple genome sequences.

Multiple sequence alignment is a difficult computational problem. There have been compelling pleas for methods to assess whole-genome multiple sequence alignments and compare the alignments produced by different tools. We assess the four ENCODE alignments, each of which aligns 28 vertebrates on 554 Mbp of total input sequence. We measure the level of agreement among the alignments and compare their coverage and accuracy. We find a disturbing lack of agreement among the alignments not only in species distant from human, but even in mouse, a well-studied model organism. Overall, the assessment shows that Pecan produces the most accurate or nearly most accurate alignment in all species and genomic location categories, while still providing coverage comparable to or better than that of the other alignments in the placental mammals. Our assessment reveals that constructing accurate whole-genome multiple sequence alignments remains a significant challenge, particularly for noncoding regions and distantly related species.

Algorithms for locating extremely conserved elements in multiple sequence alignments.

BACKGROUND: In 2004, Bejerano et al. announced the startling discovery of hundreds of "ultraconserved elements", long genomic sequences perfectly conserved across human, mouse, and rat. Their announcement stimulated a flurry of subsequent research. RESULTS: We generalize the notion of ultraconserved element in a natural way from extraordinary human-rodent conservation to extraordinary conservation over an arbitrary set of species. We call these "Extremely Conserved Elements". There is a linear time algorithm to find all such Extremely Conserved Elements in any multiple sequence alignment, provided that the conservation is required to be across all the aligned species. For the general case of conservation across an arbitrary subset of the aligned species, we show that the question of whether there exists an Extremely Conserved Element is NP-complete. We illustrate the linear time algorithm by cataloguing all 177 Extremely Conserved Elements in the currently available 44-vertebrate whole-genome alignment, and point out some of the characteristics of these elements. CONCLUSIONS: The NP-completeness in the case of conservation across an arbitrary subset of the aligned species implies that it is unlikely an efficient algorithm exists for this general case. Despite this fact, for the interesting case of conservation across all or most of the aligned species, our algorithm is efficient enough to be practical. The 177 Extremely Conserved Elements that we catalog demonstrate many of the characteristics of the original ultraconserved elements of Bejerano et al.

Meta-analysis of inter-species liver co-expression networks elucidates traits associated with common human diseases.

Co-expression networks are routinely used to study human diseases like obesity and diabetes. Systematic comparison of these networks between species has the potential to elucidate common mechanisms that are conserved between human and rodent species, as well as those that are species-specific characterizing evolutionary plasticity. We developed a semi-parametric meta-analysis approach for combining gene-gene co-expression relationships across expression profile datasets from multiple species. The simulation results showed that the semi-parametric method is robust against noise. When applied to human, mouse, and rat liver co-expression networks, our method out-performed existing methods in identifying gene pairs with coherent biological functions. We identified a network conserved across species that highlighted cell-cell signaling, cell-adhesion and sterol biosynthesis as main biological processes represented in genome-wide association study candidate gene sets for blood lipid levels. We further developed a heterogeneity statistic to test for network differences among multiple datasets, and demonstrated that genes with species-specific interactions tend to be under positive selection throughout evolution. Finally, we identified a human-specific sub-network regulated by RXRG, which has been validated to play a different role in hyperlipidemia and Type 2 diabetes between human and mouse. Taken together, our approach represents a novel step forward in integrating gene co-expression networks from multiple large scale datasets to leverage not only common information but also differences that are dataset-specific.

Assessing the discordance of multiple sequence alignments.

Multiple sequence alignments have wide applicability in many areas of computational biology, including comparative genomics, functional annotation of proteins, gene finding, and modeling evolutionary processes. Because of the computational difficulty of multiple sequence alignment and the availability of numerous tools, it is critical to be able to assess the reliability of multiple alignments. We present a tool called StatSigMA to assess whether multiple alignments of nucleotide or amino acid sequences are contaminated with one or more unrelated sequences. There are numerous applications for which StatSigMA can be used. Two such applications are to distinguish homologous sequences from nonhomologous ones and to compare alignments produced by various multiple alignment tools. We present examples of both types of applications.

TP0262 is a modulator of promoter activity of tpr Subfamily II genes of Treponema pallidum ssp. pallidum.

Transcriptional regulation in Treponema pallidum ssp. pallidum is poorly understood, primarily because this organism cannot be cultivated in vitro or genetically manipulated. We have recently shown a phase variation mechanism controlling transcription initiation of Subfamily II tpr (T. pallidumrepeat) genes (tprE, tprG and tprJ), a group of virulence factor candidates. Furthermore, the same study suggested that additional mechanisms might influence the level of transcription of these tprs. The T. pallidum genome sequence has revealed a few open reading frames with similarity to known bacterial transcription factors, including four catabolite activator protein homologues. In this work, sequences matching the Escherichia coli cAMP receptor protein (CRP) binding motif were identified in silico upstream of tprE, tprG and tprJ. Using elecrophoretic mobility shift assay and DNaseI footprinting assay, recombinant TP0262, a T. pallidum CRP homologue, was shown to bind specifically to amplicons obtained from the tpr promoters containing putative CRP binding motifs. Using a heterologous reporter system, binding of TP0262 to these promoters was shown to either increase (tprE and tprJ) or decrease (tprG) tpr promoter activity. This is the first characterization of a T. pallidum transcriptional modulator that influences tpr promoter activity.

How accurately is ncRNA aligned within whole-genome multiple alignments?

BACKGROUND: Multiple alignment of homologous DNA sequences is of great interest to biologists since it provides a window into evolutionary processes. At present, the accuracy of whole-genome multiple alignments, particularly in noncoding regions, has not been thoroughly evaluated. RESULTS: We evaluate the alignment accuracy of certain noncoding regions using noncoding RNA alignments from Rfam as a reference. We inspect the MULTIZ 17-vertebrate alignment from the UCSC Genome Browser for all the human sequences in the Rfam seed alignments. In particular, we find 638 instances of chimeric and partial alignments to human noncoding RNA elements, of which at least 225 can be improved by straightforward means. As a byproduct of our procedure, we predict many novel instances of known ncRNA families that are suggested by the alignment. CONCLUSION: MULTIZ does a fairly accurate job of aligning these genomes in these difficult regions. However, our experiments indicate that better alignments exist in some regions.

Identification of 22 candidate structured RNAs in bacteria using the CMfinder comparative genomics pipeline.

We applied a computational pipeline based on comparative genomics to bacteria, and identified 22 novel candidate RNA motifs. We predicted six to be riboswitches, which are mRNA elements that regulate gene expression on binding a specific metabolite. In separate studies, we confirmed that two of these are novel riboswitches. Three other riboswitch candidates are upstream of either a putative transporter gene in the order Lactobacillales, citric acid cycle genes in Burkholderiales or molybdenum cofactor biosynthesis genes in several phyla. The remaining riboswitch candidate, the widespread Genes for the Environment, for Membranes and for Motility (GEMM) motif, is associated with genes important for natural competence in Vibrio cholerae and the use of metal ions as electron acceptors in Geobacter sulfurreducens. Among the other motifs, one has a genetic distribution similar to a previously published candidate riboswitch, ykkC/yxkD, but has a different structure. We identified possible non-coding RNAs in five phyla, and several additional cis-regulatory RNAs, including one in epsilon-proteobacteria (upstream of purD, involved in purine biosynthesis), and one in Cyanobacteria (within an ATP synthase operon). These candidate RNAs add to the growing list of RNA motifs involved in multiple cellular processes, and suggest that many additional RNAs remain to be discovered.

A computational pipeline for high- throughput discovery of cis-regulatory noncoding RNA in prokaryotes.

Noncoding RNAs (ncRNAs) are important functional RNAs that do not code for proteins. We present a highly efficient computational pipeline for discovering cis-regulatory ncRNA motifs de novo. The pipeline differs from previous methods in that it is structure-oriented, does not require a multiple-sequence alignment as input, and is capable of detecting RNA motifs with low sequence conservation. We also integrate RNA motif prediction with RNA homolog search, which improves the quality of the RNA motifs significantly. Here, we report the results of applying this pipeline to Firmicute bacteria. Our top-ranking motifs include most known Firmicute elements found in the RNA family database (Rfam). Comparing our motif models with Rfam's hand-curated motif models, we achieve high accuracy in both membership prediction and base-pair-level secondary structure prediction (at least 75% average sensitivity and specificity on both tasks). Of the ncRNA candidates not in Rfam, we find compelling evidence that some of them are functional, and analyze several potential ribosomal protein leaders in depth.

Measuring the accuracy of genome-size multiple alignments.

Whole-genome alignments are invaluable for comparative genomics. Before doing any comparative analysis on a region of interest, one must have confidence in that region's alignment. We provide a methodology to measure the accuracy of arbitrary regions of these alignments, and apply it to the UCSC Genome Browser's 17-vertebrate alignment. We identify 9.7% (21 Mbp) of the human chromosome 1 alignment as suspiciously aligned. We present independent evidence that many of these suspicious regions represent misalignments.

Mis-translation of a computationally designed protein yields an exceptionally stable homodimer: implications for protein engineering and evolution.

We recently used computational protein design to create an extremely stable, globular protein, Top7, with a sequence and fold not observed previously in nature. Since Top7 was created in the absence of genetic selection, it provides a rare opportunity to investigate aspects of the cellular protein production and surveillance machinery that are subject to natural selection. Here we show that a portion of the Top7 protein corresponding to the final 49 C-terminal residues is efficiently mis-translated and accumulates at high levels in Escherichia coli. We used circular dichroism, size-exclusion chromatography, small-angle X-ray scattering, analytical ultra-centrifugation, and NMR spectroscopy to show that the resulting C-terminal fragment (CFr) protein adopts a compact, extremely stable, homo-dimeric structure. Based on the solution structure, we engineered an even more stable variant of CFr by disulfide-induced covalent circularisation that should be an excellent platform for design of novel functions. The accumulation of high levels of CFr exposes the high error rate of the protein translation machinery. The rarity of correspondingly stable fragments in natural proteins coupled with the observation that high quality ribosome binding sites are found to occur within E. coli protein-coding regions significantly less often than expected by random chance implies a stringent evolutionary pressure against protein sub-fragments that can independently fold into stable structures. The symmetric self-association between two identical mis-translated CFr sub-domains to generate an extremely stable structure parallels a mechanism for natural protein-fold evolution by modular recombination of protein sub-structures.

MicroFootPrinter: a tool for phylogenetic footprinting in prokaryotic genomes.

Phylogenetic footprinting is a method for the discovery of regulatory elements in a set of homologous regulatory regions, usually collected from multiple species. It does so by identifying the most conserved motifs in those homologous regions. This note describes web software that has been designed specifically for this purpose in prokaryotic genomes, making use of the phylogenetic relationships among the homologous sequences in order to make more accurate predictions. The software is called MicroFootPrinter and is available at http://bio.cs.washington.edu/software.html.

Analysis of computational approaches for motif discovery.

Recently, we performed an assessment of 13 popular computational tools for discovery of transcription factor binding sites (M. Tompa, N. Li, et al., "Assessing Computational Tools for the Discovery of Transcription Factor Binding Sites", Nature Biotechnology, Jan. 2005). This paper contains follow-up analysis of the assessment results, and raises and discusses some important issues concerning the state of the art in motif discovery methods: 1. We categorize the objective functions used by existing tools, and design experiments to evaluate whether any of these objective functions is the right one to optimize. 2. We examine various features of the data sets that were used in the assessment, such as sequence length and motif degeneracy, and identify which features make data sets hard for current motif discovery tools. 3. We identify an important feature that has not yet been used by existing tools and propose a new objective function that incorporates this feature.

Discovery of regulatory elements in vertebrates through comparative genomics.

We have analyzed issues of reliability in studies in which comparative genomic approaches have been applied to the discovery of regulatory elements at a genome-wide level in vertebrates. We point out some potential problems with such studies, including difficulties in accurately identifying orthologous promoter regions. Many of these subtle analytical problems have become apparent only when studying the more complex vertebrate genomes. By determining motif reliability, we compared existing tools when applied to the discovery of vertebrate regulatory elements. We then used a statistical clustering method to produce a computational catalog of high quality putative regulatory elements from vertebrates, some of which are widely conserved among vertebrates and many of which are novel regulatory elements. The results provide a glimpse into the wealth of information that comparative genomics can yield and suggest the need for further improvement of genome-wide comparative computational techniques.

Statistics of local multiple alignments.

SUMMARY: BLAST statistics have been shown to be extremely useful for searching for significant similarity hits, for amino acid and nucleotide sequences. Although these statistics are well understood for pairwise comparisons, there has been little success developing statistical scores for multiple alignments. In particular, there is no score for multiple alignment that is well founded and treated as a standard. We extend the BLAST theory to multiple alignments. Following some simple assumptions, we present and justify a significance score for multiple segments of a local multiple alignment. We demonstrate its usefulness in distinguishing high and moderate quality multiple alignments from low quality ones, with supporting experiments on orthologous vertebrate promoter sequences.

Assessing computational tools for the discovery of transcription factor binding sites.

The prediction of regulatory elements is a problem where computational methods offer great hope. Over the past few years, numerous tools have become available for this task. The purpose of the current assessment is twofold: to provide some guidance to users regarding the accuracy of currently available tools in various settings, and to provide a benchmark of data sets for assessing future tools.

PhyME: a probabilistic algorithm for finding motifs in sets of orthologous sequences.

BACKGROUND: This paper addresses the problem of discovering transcription factor binding sites in heterogeneous sequence data, which includes regulatory sequences of one or more genes, as well as their orthologs in other species. RESULTS: We propose an algorithm that integrates two important aspects of a motif's significance - overrepresentation and cross-species conservation - into one probabilistic score. The algorithm allows the input orthologous sequences to be related by any user-specified phylogenetic tree. It is based on the Expectation-Maximization technique, and scales well with the number of species and the length of input sequences. We evaluate the algorithm on synthetic data, and also present results for data sets from yeast, fly, and human. CONCLUSIONS: The results demonstrate that the new approach improves motif discovery by exploiting multiple species information.

Evolutionarily conserved sequence elements that positively regulate IFN-gamma expression in T cells.

Our understanding of mechanisms by which the expression of IFN-gamma is regulated is limited. Herein, we identify two evolutionarily conserved noncoding sequence elements (IFNgCNS1 and IFNg CNS2) located approximately 5 kb upstream and approximately 18 kb downstream of the initiation codon of the murine Ifng gene. When linked to the murine Ifng gene (-3.4 to +5.6 kb) and transiently transfected into EL-4 cells, these elements clearly enhanced IFN-gamma expression in response to ionomycin and phorbol 12-myristate 13-acetate and weakly enhanced expression in response to T-bet. A DNase I hypersensitive site and extragenic transcripts at IFNgCNS2 correlated positively with the capacity of primary T cell subsets to produce IFN-gamma. Transcriptionally favorable histone modifications in the Ifng promoter, intronic regions, IFNgCNS2, and, although less pronounced, IFNgCNS1 increased as naïve T cells differentiated into IFN-gamma-producing effector CD8+ and T helper (TH) 1 T cells, but not into TH2 T cells. Like IFN-gamma expression, these histone modifications were T-bet-dependent in CD4+ cells, but not CD8+ T cells. These findings define two distal regulatory elements associated with T cell subset-specific IFN-gamma expression.

YMF: A program for discovery of novel transcription factor binding sites by statistical overrepresentation.

A fundamental challenge facing biologists is to identify DNA binding sites for unknown regulatory factors, given a collection of genes believed to be coregulated. The program YMF identifies good candidates for such binding sites by searching for statistically overrepresented motifs. More specifically, YMF enumerates all motifs in the search space and is guaranteed to produce those motifs with greatest z-scores. This note describes the YMF web software, available at http://bio.cs.washington.edu/software.html.