Insertion sequence-caused large-scale rearrangements in the genome of Escherichia coli.

A majority of large-scale bacterial genome rearrangements involve mobile genetic elements such as insertion sequence (IS) elements. Here we report novel insertions and excisions of IS elements and recombination between homologous IS elements identified in a large collection of Escherichia coli mutation accumulation lines by analysis of whole genome shotgun sequencing data. Based on 857 identified events (758 IS insertions, 98 recombinations and 1 excision), we estimate that the rate of IS insertion is 3.5 × 10(-4) insertions per genome per generation and the rate of IS homologous recombination is 4.5 × 10(-5) recombinations per genome per generation. These events are mostly contributed by the IS elements IS1, IS2, IS5 and IS186 Spatial analysis of new insertions suggest that transposition is biased to proximal insertions, and the length spectrum of IS-caused deletions is largely explained by local hopping. For any of the ISs studied there is no region of the circular genome that is favored or disfavored for new insertions but there are notable hotspots for deletions. Some elements have preferences for non-coding sequence or for the beginning and end of coding regions, largely explained by target site motifs. Interestingly, transposition and deletion rates remain constant across the wild-type and 12 mutant E. coli lines, each deficient in a distinct DNA repair pathway. Finally, we characterized the target sites of four IS families, confirming previous results and characterizing a highly specific pattern at IS186 target-sites, 5'-GGGG(N6/N7)CCCC-3'. We also detected 48 long deletions not involving IS elements.

Determinants of spontaneous mutation in the bacterium Escherichia coli as revealed by whole-genome sequencing.

A complete understanding of evolutionary processes requires that factors determining spontaneous mutation rates and spectra be identified and characterized. Using mutation accumulation followed by whole-genome sequencing, we found that the mutation rates of three widely diverged commensal Escherichia coli strains differ only by about 50%, suggesting that a rate of 1-2 × 10(-3) mutations per generation per genome is common for this bacterium. Four major forces are postulated to contribute to spontaneous mutations: intrinsic DNA polymerase errors, endogenously induced DNA damage, DNA damage caused by exogenous agents, and the activities of error-prone polymerases. To determine the relative importance of these factors, we studied 11 strains, each defective for a major DNA repair pathway. The striking result was that only loss of the ability to prevent or repair oxidative DNA damage significantly impacted mutation rates or spectra. These results suggest that, with the exception of oxidative damage, endogenously induced DNA damage does not perturb the overall accuracy of DNA replication in normally growing cells and that repair pathways may exist primarily to defend against exogenously induced DNA damage. The thousands of mutations caused by oxidative damage recovered across the entire genome revealed strong local-sequence biases of these mutations. Specifically, we found that the identity of the 3' base can affect the mutability of a purine by oxidative damage by as much as eightfold.

Interactions and Localization of Escherichia coli Error-Prone DNA Polymerase IV after DNA Damage.

UNLABELLED: Escherichia coli's DNA polymerase IV (Pol IV/DinB), a member of the Y family of error-prone polymerases, is induced during the SOS response to DNA damage and is responsible for translesion bypass and adaptive (stress-induced) mutation. In this study, the localization of Pol IV after DNA damage was followed using fluorescent fusions. After exposure of E. coli to DNA-damaging agents, fluorescently tagged Pol IV localized to the nucleoid as foci. Stepwise photobleaching indicated ∼60% of the foci consisted of three Pol IV molecules, while ∼40% consisted of six Pol IV molecules. Fluorescently tagged Rep, a replication accessory DNA helicase, was recruited to the Pol IV foci after DNA damage, suggesting that the in vitro interaction between Rep and Pol IV reported previously also occurs in vivo. Fluorescently tagged RecA also formed foci after DNA damage, and Pol IV localized to them. To investigate if Pol IV localizes to double-strand breaks (DSBs), an I-SceI endonuclease-mediated DSB was introduced close to a fluorescently labeled LacO array on the chromosome. After DSB induction, Pol IV localized to the DSB site in ∼70% of SOS-induced cells. RecA also formed foci at the DSB sites, and Pol IV localized to the RecA foci. These results suggest that Pol IV interacts with RecA in vivo and is recruited to sites of DSBs to aid in the restoration of DNA replication. IMPORTANCE: DNA polymerase IV (Pol IV/DinB) is an error-prone DNA polymerase capable of bypassing DNA lesions and aiding in the restart of stalled replication forks. In this work, we demonstrate in vivo localization of fluorescently tagged Pol IV to the nucleoid after DNA damage and to DNA double-strand breaks. We show colocalization of Pol IV with two proteins: Rep DNA helicase, which participates in replication, and RecA, which catalyzes recombinational repair of stalled replication forks. Time course experiments suggest that Pol IV recruits Rep and that RecA recruits Pol IV. These findings provide in vivo evidence that Pol IV aids in maintaining genomic stability not only by bypassing DNA lesions but also by participating in the restoration of stalled replication forks.

Rate and molecular spectrum of spontaneous mutations in the bacterium Escherichia coli as determined by whole-genome sequencing.

Knowledge of the rate and nature of spontaneous mutation is fundamental to understanding evolutionary and molecular processes. In this report, we analyze spontaneous mutations accumulated over thousands of generations by wild-type Escherichia coli and a derivative defective in mismatch repair (MMR), the primary pathway for correcting replication errors. The major conclusions are (i) the mutation rate of a wild-type E. coli strain is ~1 × 10(-3) per genome per generation; (ii) mutations in the wild-type strain have the expected mutational bias for G:C > A:T mutations, but the bias changes to A:T > G:C mutations in the absence of MMR; (iii) during replication, A:T > G:C transitions preferentially occur with A templating the lagging strand and T templating the leading strand, whereas G:C > A:T transitions preferentially occur with C templating the lagging strand and G templating the leading strand; (iv) there is a strong bias for transition mutations to occur at 5'ApC3'/3'TpG5' sites (where bases 5'A and 3'T are mutated) and, to a lesser extent, at 5'GpC3'/3'CpG5' sites (where bases 5'G and 3'C are mutated); (v) although the rate of small (≤4 nt) insertions and deletions is high at repeat sequences, these events occur at only 1/10th the genomic rate of base-pair substitutions. MMR activity is genetically regulated, and bacteria isolated from nature often lack MMR capacity, suggesting that modulation of MMR can be adaptive. Thus, comparing results from the wild-type and MMR-defective strains may lead to a deeper understanding of factors that determine mutation rates and spectra, how these factors may differ among organisms, and how they may be shaped by environmental conditions.

Assessment of Chemical Toxicity in Adult Drosophila Melanogaster.

Human industries generate hundreds of thousands of chemicals, many of which have not been adequately studied for environmental safety or effects on human health. This deficit of chemical safety information is exacerbated by current testing methods in mammals that are expensive, labor-intensive, and time-consuming. Recently, scientists and regulators have been working to develop new approach methodologies (NAMs) for chemical safety testing that are cheaper, more rapid, and reduce animal suffering. One of the key NAMs to emerge is the use of invertebrate organisms as replacements for mammalian models to elucidate conserved chemical modes of action across distantly related species, including humans. To advance these efforts, here, we describe a method that uses the fruit fly, Drosophila melanogaster, to assess chemical safety. The protocol describes a simple, rapid, and inexpensive procedure to measure the viability and feeding behavior of exposed adult flies. In addition, the protocol can be easily adapted to generate samples for genomic and metabolomic approaches. Overall, the protocol represents an important step forward in establishing Drosophila as a standard model for use in precision toxicology.

Identification of novel split-GAL4 drivers for the characterization of enteroendocrine cells in the Drosophila melanogaster midgut.

The Drosophila melanogaster midgut is commonly studied as a model epithelial tissue for many reasons, one of which is the presence of a diverse population of secretory cells called enteroendocrine cells. Subpopulations of these cells secrete various combinations of peptide hormones which have systemic effects on the organism. Many of these hormones are also produced in the Drosophila brain. The split-GAL4 system has been useful for identifying and manipulating discrete groups of cells, but previously characterized split-GAL4 drivers have not driven expression in high proportions of enteroendocrine cells. In this study, we screened candidate split-GAL4 drivers for enteroendocrine cell expression using known reference drivers for this cell type and discovered a new split-GAL4 driver pair that confers expression in a greater number of enteroendocrine cells than previously characterized driver pairs. The new pair demonstrates less brain expression, thereby providing better tools for disentangling the physiological roles of gut- and brain-secreted peptides. We also identified additional split-GAL4 drivers that promote expression in discrete subpopulations of enteroendocrine cells. Overall, the tools reported here will help researchers better target enteroendocrine cell subpopulations.

Enteroendocrine cell expression of split-GAL4 drivers bearing regulatory sequences associated with panneuronally expressed genes in Drosophila melanogaster.

In Drosophila melanogaster , hormone-secreting enteroendocrine cells are important for communication from the midgut to other tissues. Many lexA, GAL4, and split-GAL4 drivers that direct gene expression in enteroendocrine cells also confer expression in hormone-secreting cells of the central nervous system. This study examines the midgut expression of selected lexA, GAL4, and split-GAL4 transgenes carrying enhancer fragments previously associated with panneuronal gene expression to assess the experimental usefulness of these drivers for distinguishing the endocrine influences of CNS versus midgut cells on physiological processes.

Axoneme beta-tubulin sequence determines attachment of outer dynein arms.

Axonemes of motile eukaryotic cilia and flagella have a conserved structure of nine doublet microtubules surrounding a central pair of microtubules. Outer and inner dynein arms on the doublets mediate axoneme motility [1]. Outer dynein arms (ODAs) attach to the doublets at specific interfaces [2-5]. However, the molecular contacts of ODA-associated proteins with tubulins of the doublet microtubules are not known. We report here that attachment of ODAs requires glycine 56 in the beta-tubulin internal variable region (IVR). We show that in Drosophila spermatogenesis, a single amino acid change at this position results in sperm axonemes markedly deficient in ODAs. Moreover, we found that axonemal beta-tubulins throughout the phylogeny have invariant glycine 56 and a strongly conserved IVR, whereas nonaxonemal beta-tubulins vary widely in IVR sequences. Our data reveal a deeply conserved physical requirement for assembly of the macromolecular architecture of the motile axoneme. Amino acid 56 projects into the microtubule lumen [6]. Imaging studies of axonemes indicate that several proteins may interact with the doublet-microtubule lumen [3, 4, 7, 8]. This region of beta-tubulin may determine the conformation necessary for correct attachment of ODAs, or there may be sequence-specific interaction between beta-tubulin and a protein involved in ODA attachment or stabilization.

Identification of Split-GAL4 Drivers and Enhancers That Allow Regional Cell Type Manipulations of the Drosophila melanogaster Intestine.

The Drosophila adult midgut is a model epithelial tissue composed of a few major cell types with distinct regional identities. One of the limitations to its analysis is the lack of tools to manipulate gene expression based on these regional identities. To overcome this obstacle, we applied the intersectional split-GAL4 system to the adult midgut and report 653 driver combinations that label cells by region and cell type. We first identified 424 split-GAL4 drivers with midgut expression from ∼7300 drivers screened, and then evaluated the expression patterns of each of these 424 when paired with three reference drivers that report activity specifically in progenitor cells, enteroendocrine cells, or enterocytes. We also evaluated a subset of the drivers expressed in progenitor cells for expression in enteroblasts using another reference driver. We show that driver combinations can define novel cell populations by identifying a driver that marks a distinct subset of enteroendocrine cells expressing genes usually associated with progenitor cells. The regional cell type patterns associated with the entire set of driver combinations are documented in a freely available website, providing information for the design of thousands of additional driver combinations to experimentally manipulate small subsets of intestinal cells. In addition, we show that intestinal enhancers identified with the split-GAL4 system can confer equivalent expression patterns on other transgenic reporters. Altogether, the resource reported here will enable more precisely targeted gene expression for studying intestinal processes, epithelial cell functions, and diseases affecting self-renewing tissues.

Cooperativity between the beta-tubulin carboxy tail and the body of the molecule is required for microtubule function.

Using Drosophila spermatogenesis as a model, we show that function of the beta-tubulin C-terminal tail (CTT) is not independent of the body of the molecule. For optimal microtubule function, the beta-tubulin CTT and body must match. beta2 is the only beta-tubulin used in meiosis and spermatid differentiation. beta1-tubulin is used in basal bodies, but beta1 cannot replace beta2. However, when beta1 is co-expressed with beta2, both beta-tubulins are equally incorporated into all microtubules, and males exhibit near wild type fertility. In contrast, co-expression of beta2beta1C and beta1beta2C, two reciprocal chimeric molecules with bodies and tails swapped, results in defects in meiosis, cytoskeletal microtubules, and axonemes; males produce few functional sperm and few or no progeny. In these experiments, all the same beta-tubulin parts are present, but unlike the co-assembled native beta-tubulins, the "trans" configuration of the co-assembled chimeras is poorly functional. Our data thus reveal essential intra-molecular interactions between the CTT and other parts of the beta-tubulin molecule, even though the CTT is a flexible surface feature of tubulin heterodimers and microtubules. In addition, we show that Drosophila sperm tail length depends on the total tubulin pool available for axoneme assembly and spermatid elongation. D. melanogaster and other Drosophila species have extraordinarily long sperm tails, the length of which is remarkably constant in wild type flies. We show that in males of experimental genotypes that express wild type tubulins but have half the amount of the normal tubulin pool size, sperm tails are substantially shorter than wild type.

Determinants of Base-Pair Substitution Patterns Revealed by Whole-Genome Sequencing of DNA Mismatch Repair Defective Escherichia coli.

Mismatch repair (MMR) is a major contributor to replication fidelity, but its impact varies with sequence context and the nature of the mismatch. Mutation accumulation experiments followed by whole-genome sequencing of MMR-defective Escherichia coli strains yielded ≈30,000 base-pair substitutions (BPSs), revealing mutational patterns across the entire chromosome. The BPS spectrum was dominated by A:T to G:C transitions, which occurred predominantly at the center base of 5'NAC3'+5'GTN3' triplets. Surprisingly, growth on minimal medium or at low temperature attenuated these mutations. Mononucleotide runs were also hotspots for BPSs, and the rate at which these occurred increased with run length. Comparison with ≈2000 BPSs accumulated in MMR-proficient strains revealed that both kinds of hotspots appeared in the wild-type spectrum and so are likely to be sites of frequent replication errors. In MMR-defective strains transitions were strand biased, occurring twice as often when A and C rather than T and G were on the lagging-strand template. Loss of nucleotide diphosphate kinase increases the cellular concentration of dCTP, which resulted in increased rates of mutations due to misinsertion of C opposite A and T. In an mmr ndk double mutant strain, these mutations were more frequent when the template A and T were on the leading strand, suggesting that lagging-strand synthesis was more error-prone, or less well corrected by proofreading, than was leading strand synthesis.

Adaptive evolution of bindin in the genus Heliocidaris is correlated with the shift to direct development.

Sea urchins are widely used to study both fertilization and development. In this study we combine the two fields to examine the evolution of reproductive isolation in the genus Heliocidaris. Heliocidaris tuberculata develops indirectly via a feeding larva, whereas the only other species in the genus, H. erythrogramma, has evolved direct development through a nonfeeding larva. We estimated the time of divergence between H. erythrogramma and H. tuberculata from mitochondrial DNA divergence, quantified levels of gametic compatibility between the two species in cross-fertilization assays, and examined the mode of evolution of the sperm protein bindin by sequencing multiple alleles of the two species. Bindin is the major component of the sea urchin sperm acrosomal vesicle, and is involved in sperm-egg attachment and fusion. Based on our analyses, we conclude that: the two species of Heliocidaris diverged less than five million years ago, indicating that direct development can evolve rapidly in sea urchins; since their divergence, the two species have become gametically incompatible; Heliocidaris bindin has evolved under positive selection; and this positive selection is concentrated on the branch leading to H. erythrogramma. Three hypotheses can explain the observed pattern of selection on bindin: (1) it is a correlated response to the evolution of direct development in H. erythrogramma; (2) it is the result of an intraspecific process acting in H. erythrogramma but not in H. tuberculata; or (3) it is the product of reinforcement on the species that invests more energy into each egg to avoid hybridization.

Detection of structural variants involving repetitive regions in the reference genome.

Next-generation sequencing techniques are now commonly used to characterize structural variations (SVs) in population genomics and elucidate their associations with phenotypes. Many of the computational tools developed for detecting structural variations work by mapping paired-end reads to a reference genome and identifying the discordant read-pairs whose mapped loci in the reference genome deviate from the expected insert size and orientation. However, repetitive regions in the reference genome represent a major challenge in SV detection, because the paired-end reads from these regions may be mapped to multiple loci in the reference genome, resulting in spuriously discordant read-pairs. To address this issue, we have developed an algorithmic approach for read mapping and SV detection based on the framework of A-Bruijn graphs. Instead of mapping reads to a linear sequence of the reference genome, we propose to map reads onto the A-Bruijn graph constructed from the reference genome in which all instances of the same repeat are collapsed into a single edge. As a result, any given read, either from repetitive regions or not, will be mapped to a unique location in the A-Bruijn graph, and each discordant read-pair in the A-Bruijn graph indicates a potentially true SV event. We also developed a simple clustering algorithm to derive valid clusters of these discordant read-pairs, each supporting a different SV event. Finally, we demonstrate the performance of this approach, compared to existing approaches, by identifying transposition events of insertion sequence (IS) elements, a class of simple mobile genetic elements (MGEs), in E. coli by using simulated and real paired-end sequence data acquired from E. coli mutation accumulation lines.

On the mutational topology of the bacterial genome.

By sequencing the genomes of 34 mutation accumulation lines of a mismatch-repair defective strain of Escherichia coli that had undergone a total of 12,750 generations, we identified 1625 spontaneous base-pair substitutions spread across the E. coli genome. These mutations are not distributed at random but, instead, fall into a wave-like spatial pattern that is repeated almost exactly in mirror image in the two separately replicated halves of the bacterial chromosome. The pattern is correlated to genomic features, with mutation densities greatest in regions predicted to have high superhelicity. Superimposed upon this pattern are regional hotspots, some of which are located where replication forks may collide or be blocked. These results suggest that, as they traverse the chromosome, the two replication forks encounter parallel structural features that change the fidelity of DNA replication.

A molecularly defined duplication set for the X chromosome of Drosophila melanogaster.

We describe a molecularly defined duplication kit for the X chromosome of Drosophila melanogaster. A set of 408 overlapping P[acman] BAC clones was used to create small duplications (average length 88 kb) covering the 22-Mb sequenced portion of the chromosome. The BAC clones were inserted into an attP docking site on chromosome 3L using ΦC31 integrase, allowing direct comparison of different transgenes. The insertions complement 92% of the essential and viable mutations and deletions tested, demonstrating that almost all Drosophila genes are compact and that the current annotations of the genome are reasonably accurate. Moreover, almost all genes are tolerated at twice the normal dosage. Finally, we more precisely mapped two regions at which duplications cause diplo-lethality in males. This collection comprises the first molecularly defined duplication set to cover a whole chromosome in a multicellular organism. The work presented removes a long-standing barrier to genetic analysis of the Drosophila X chromosome, will greatly facilitate functional assays of X-linked genes in vivo, and provides a model for functional analyses of entire chromosomes in other species.

Axoneme specialization embedded in a "generalist" beta-tubulin.

The relationship between the primary structure of the beta-tubulin C-terminal tail (CTT) and axoneme structure and function is explored using the spermatogenesis-specific beta2-tubulin of Drosophila. We previously showed that all beta-tubulins used for motile 9 + 2 axonemes contain a conserved sequence motif in the proximal part of the CTT, the beta-tubulin axoneme motif. The differential ability of tubulin isoforms and abilities of beta2-tubulin C-terminal truncations to form axonemes led us to hypothesize that the axoneme motif is essential for axoneme formation and the distal half of the CTT was less important. The studies we report here indicate that it is not that simple. Unexpectedly, some changes in the core sequence of the axoneme motif did not disrupt formation of motile axonemes. And, while deletion of the distal CTT did not disrupt the ability to produce functional sperm [Popodi et al., Cell Motil Cytoskeleton 2005;62:48-64], changing the amino acid sequence in this region can. Thus both regions are important. The deep conservation of the axoneme motif in all eukaryotic groups implies that the presence of the sequence motif confers a functional advantage. The central pair is the axoneme structure most sensitive to perturbations in tubulin molecules; we hypothesize central pair assembly is facilitated by the presence of this motif. Our data reveal that beta2-tubulin has robust properties for axoneme assembly, and that axonemal specializations are embedded in both the CTT and the body of the beta2 molecule.

The proximal region of the beta-tubulin C-terminal tail is sufficient for axoneme assembly.

We have used Drosophila testis-specific beta2-tubulin to determine sequence requirements for different microtubules. The beta2-tubulin C-terminal tail has unique sperm-specific functions [Dev Biol 158:267-286 (2003)] and is also important for forming stable heterodimers with alpha-tubulin, a general function common to all microtubules [Mol Biol Cell 12(7):2185-2194 (2001)]. beta-tubulins utilized in motile 9 + 2 axonemes contain a C-terminal sequence "axoneme motif" [Science 275 (1997) 70-73]. C-terminal truncated beta2-tubulin cannot form the sperm tail axoneme. Here we show that a partially truncated beta2-tubulin (beta2Delta7) containing only the proximal portion of the C-terminal tail, including the axoneme motif, can support production of functional motile sperm. We conclude that these proximal eight amino acids specify the binding site for protein(s) essential to support assembly of the motile axoneme. Males that express beta2Delta7, although they are fertile, produce fewer sperm than wild type males. Beta2Delta7 causes a slightly increased error rate in spermatogenesis attributable to loss of stabilizing properties intrinsic to the full-length C-terminal tail. Therefore, beta2Delta7 males would be at a selective disadvantage and it is likely that the full-length C-terminus would be essential in the wild and in evolution.

Patterns of gene expression in the developing adult sea urchin central nervous system reveal multiple domains and deep-seated neural pentamery.

The adult sea urchin central nervous system (CNS) is composed of five radial nerve cords connected to a circular nerve ring. Although much is known about the molecular mechanisms underlying the development and function of the nervous systems of many invertebrate and vertebrate species, virtually nothing is known about these processes in echinoderms. We have isolated a set of clones from a size-selected cDNA library prepared from the nervous system of the sea urchin Heliocidaris erythrogramma for use as probes. A total of 117 expressed sequence clones were used to search the GenBank database. Identified messages include genes that encode signaling proteins, cytoskeletal elements, cell surface proteins and receptors, cell proliferation and differentiation factors, transport and channel proteins, and a RNA DEAD box helicase. Expression was analyzed by RNA gel blot hybridization to document expression through development. Many of the genes have apparently neural limited expression and function, but some have been co-opted into new roles, notably associated with exocytotic events at fertilization. Localization of gene expression by whole-mount in situ hybridization shows that the morphologically simple sea urchin radial CNS exhibits complex organization into localized transcriptional domains. The transcription patterns reflect the morphological pentamery of the echinoderm CNS and provide no indication of an underlying functional bilateral symmetry in the CNS.

Regulatory punctuated equilibrium and convergence in the evolution of developmental pathways in direct-developing sea urchins.

We made hybrid crosses between closely and distantly related sea urchin species to test two hypotheses about the evolution of gene regulatory systems in the evolution of ontogenetic pathways and larval form. The first hypothesis is that gene regulatory systems governing development evolve in a punctuational manner during periods of rapid morphological evolution but are relatively stable over long periods of slow morphological evolution. We compared hybrids between direct and indirect developers from closely and distantly related families. Hybrids between eggs of the direct developer Heliocidaris erythrogramma and sperm of the 4-million year distant species H. tuberculata, an indirect developer, restored feeding larval structures and paternal gene expression that were lost in the evolution of the direct-developing maternal parent. Hybrids resulting from the cross between eggs of H. erythrogramma and sperm of the 40-million year distant indirect-developer Pseudoboletia maculata are strikingly similar to hybrids between the congeneric hybrids. The marked similarities in ontogenetic trajectory and morphological outcome in crosses of involving either closely or distantly related indirect developing species indicates that their regulatory mechanisms interact with those of H. erythrogramma in the same way, supporting remarkable conservation of molecular control pathways among indirect developers. Second, we tested the hypothesis that convergent developmental pathways in independently evolved direct developers reflect convergence of the underlying regulatory systems. Crosses between two independently evolved direct-developing species from two 70-million year distant families, H. erythrogramma and Holopneustes purpurescens, produced harmoniously developing hybrid larvae that maintained the direct mode of development and did not exhibit any obvious restoration of indirect-developing features. These results are consistent with parallel evolution of direct-developing features in these two lineages.