Search

DNA methylation in bacteria frequently serves as a simple immune system, allowing recognition of DNA from foreign sources, such as phages or selfish genetic elements. However, DNA methylation also affects other cell phenotypes in a heritable manner (i.e. epigenetically). While there are several examples of methylation affecting transcription in an epigenetic manner in highly localized contexts, it is not well-established how frequently methylation serves a more general epigenetic function over larger genomic scales. To address this question, here we use Oxford Nanopore sequencing to profile DNA modification marks in three natural isolates of Escherichia coli. We first identify the DNA sequence motifs targeted by the methyltransferases in each strain. We then quantify the frequency of methylation at each of these motifs across the entire genome in different growth conditions. We find that motifs in specific regions of the genome consistently exhibit high or low levels of methylation. Furthermore, we show that there are replicable and consistent differences in methylated regions across different growth conditions. This suggests that during growth, E. coli transiently differentiate into distinct methylation states that depend on the growth state, raising the possibility that measuring DNA methylation alone can be used to infer bacterial growth states without additional information such as transcriptome or proteome data. These results show the utility of using Oxford Nanopore sequencing as an economic means to infer DNA methylation status. They also provide new insights into the dynamics of methylation during bacterial growth and provide evidence of differentiated cell states, a transient analog to what is observed in the differentiation of cell types in multicellular organisms.

Escherichia coli is commonly considered a host-associated bacterium. However, there is evidence that some strains occupy environmental (non-host-associated) niches. Here, we report the complete genomes of 47 Escherichia coli environmental isolates. These will be useful for understanding the dynamics of plasmids, phages, and other repetitive genetic elements.

Background The first step in understanding ecological community diversity and dynamics is quantifying community membership. An increasingly common method for doing so is through metagenomics. Because of the rapidly increasing popularity of this approach, a large number of computational tools and pipelines are available for analysing metagenomic data. However, the majority of these tools have been designed and benchmarked using highly accurate short read data (i.e. Illumina), with few studies benchmarking classification accuracy for long error-prone reads (PacBio or Oxford Nanopore). In addition, few tools have been benchmarked for non-microbial communities. Results Here we compare simulated long reads from Oxford Nanopore and Pacific Biosciences (PacBio) with high accuracy Illumina read sets to systematically investigate the effects of sequence length and taxon type on classification accuracy for metagenomic data from both microbial and non-microbial communities. We show that very generally, classification accuracy is far lower for non-microbial communities, even at low taxonomic resolution (e.g. family rather than genus). We then show that for two popular taxonomic classifiers, long reads can significantly increase classification accuracy, and this is most pronounced for non-microbial communities. Conclusions This work provides insight on the expected accuracy for metagenomic analyses for different taxonomic groups, and establishes the point at which read length becomes more important than error rate for assigning the correct taxon.

Beatrix, Carthage, Daegal, Dulcie, Fancypants, Fenn, Inca, Naira, and Robyn are newly isolated bacteriophages capable of infecting Mycolicibacterium smegmatis mc2 155. We discovered, sequenced, and annotated these New Zealand bacteriophages., Beatrix, Carthage, Daegal, Dulcie, Fancypants, Fenn, Inca, Naira, and Robyn are newly isolated bacteriophages capable of infecting Mycolicibacterium smegmatis mc2 155. We discovered, sequenced, and annotated these New Zealand bacteriophages. These phages illustrate that New Zealand harbors a selection of the highly diverse and distributed mycobacteriophage clusters found globally.

Seven mycobacteriophages from distinct geographical locations were isolated, using Mycobacterium smegmatis mc 2 155 as the host, and then purified and sequenced. All of the genomes are related to cluster A mycobacteriophages, BobSwaget and Lokk in subcluster A2; Fred313, KADY, Stagni, and StepMih in subcluster A3; and MyraDee in subcluster A18, the first phage to be assigned to that subcluster.

GABA(B) receptors assemble from principle and auxiliary subunits. The principle subunits GABA(B1) and GABA(B2) form functional heteromeric GABA(B(1,2)) receptors that associate with homotetramers of auxiliary KCTD8, -12, -12b, or -16 (named after their K(+) channel tetramerization domain) subunits. These auxiliary subunits constitute receptor subtypes with distinct functional properties. KCTD12 and -12b generate desensitizing receptor responses while KCTD8 and -16 generate largely non-desensitizing receptor responses. The structural elements of the KCTDs underlying these differences in desensitization are unknown. KCTDs are modular proteins comprising a T1 tetramerization domain, which binds to GABA(B2), and a H1 homology domain. KCTD8 and -16 contain an additional C-terminal H2 homology domain that is not sequence-related to the H1 domains. No functions are known for the H1 and H2 domains. Here we addressed which domains and sequence motifs in KCTD proteins regulate desensitization of the receptor response. We found that the H1 domains in KCTD12 and -12b mediate desensitization through a particular sequence motif, T/NFLEQ, which is not present in the H1 domains of KCTD8 and -16. In addition, the H2 domains in KCTD8 and -16 inhibit desensitization when expressed C-terminal to the H1 domains but not when expressed as a separate protein in trans. Intriguingly, the inhibitory effect of the H2 domain is sequence-independent, suggesting that the H2 domain sterically hinders desensitization by the H1 domain. Evolutionary analysis supports that KCTD12 and -12b evolved desensitizing properties by liberating their H1 domains from antagonistic H2 domains and acquisition of the T/NFLEQ motif.

Studies of microbial evolutionary dynamics are being transformed by the availability of affordable high-throughput sequencing technologies, which allow whole-genome sequencing of hundreds of related taxa in a single study. Reconstructing a phylogenetic tree of these taxa is generally a crucial step in any evolutionary analysis. Instead of constructing genome assemblies for all taxa, annotating these assemblies, and aligning orthologous genes, many recent studies 1) directly map raw sequencing reads to a single reference sequence, 2) extract single nucleotide polymorphisms (SNPs), and 3) infer the phylogenetic tree using maximum likelihood methods from the aligned SNP positions. However, here we show that, when using such methods to reconstruct phylogenies from sets of simulated sequences, both the exclusion of nonpolymorphic positions and the alignment to a single reference genome, introduce systematic biases and errors in phylogeny reconstruction. To address these problems, we developed a new method that combines alignments from mappings to multiple reference sequences and show that this successfully removes biases from the reconstructed phylogenies. We implemented this method as a web server named REALPHY (Reference sequence Alignment-based Phylogeny builder), which fully automates phylogenetic reconstruction from raw sequencing reads., Molecular Biology and Evolution, 31 (5), ISSN:0737-4038, ISSN:1537-1719

X-linked lymphoproliferative disease (XLP) is a primary immunodeficiency characterized by extreme susceptibility to Epstein-Barr virus. The XLP disease gene product SH2D1A (SAP) interacts via its SH2 domain with a motif (TIYXXV) present in the cytoplasmic tail of the cell-surface receptors CD150/SLAM, CD84, CD229/Ly-9, and CD244/2B4. Characteristically, the SH2D1A three-pronged interaction with Tyr(281) of CD150 can occur in absence of phosphorylation. Here we analyze the effect of SH2D1A protein missense mutations identified in 10 XLP families. Two sets of mutants were found: (i) mutants with a marked decreased protein half-life (e.g. Y7C, S28R, Q99P, P101L, V102G, and X129R) and (ii) mutants with structural changes that differently affect the interaction with the four receptors. In the second group, mutations that disrupt the interaction between the SH2D1A hydrophobic cleft and Val +3 of its binding motif (e.g. T68I) and mutations that interfere with the SH2D1A phosphotyrosine-binding pocket (e.g. C42W) abrogated SH2D1A binding to all four receptors. Surprisingly, a mutation in SH2D1A able to interfere with Thr -2 of the CD150 binding motif (mutant T53I) severely impaired non-phosphotyrosine interactions while preserving unaffected the binding of SH2D1A to phosphorylated CD150. Mutant T53I, however, did not bind to CD229 and CD224, suggesting that SH2D1A controls several critical signaling pathways in T and natural killer cells. Because no correlation is present between identified types of mutations and XLP patient clinical presentation, additional unidentified genetic or environmental factors must play a strong role in XLP disease manifestations.

Figure S6 is incorrectly a duplicate of Figure S7. The correct Figure S6 can be viewed here: Click here for additional data file.(244K, pdf) [^]

Gene expression is subject to random perturbations that lead to fluctuations in the rate of protein production. As a consequence, for any given protein, genetically identical organisms living in a constant environment will contain different amounts of that particular protein, resulting in different phenotypes. This phenomenon is known as “phenotypic noise.” In bacterial systems, previous studies have shown that, for specific genes, both transcriptional and translational processes affect phenotypic noise. Here, we focus on how the promoter regions of genes affect noise and ask whether levels of promoter-mediated noise are correlated with genes' functional attributes, using data for over 60% of all promoters in Escherichia coli. We find that essential genes and genes with a high degree of evolutionary conservation have promoters that confer low levels of noise. We also find that the level of noise cannot be attributed to the evolutionary time that different genes have spent in the genome of E. coli. In contrast to previous results in eukaryotes, we find no association between promoter-mediated noise and gene expression plasticity. These results are consistent with the hypothesis that, in bacteria, natural selection can act to reduce gene expression noise and that some of this noise is controlled through the sequence of the promoter region alone., PLoS Genetics, 8 (1), ISSN:1553-7390, ISSN:1553-7404

Genetically identical populations of unicellular organisms often show marked variation in some phenotypic traits. To investigate the molecular causes and possible biological functions of this phenotypic noise, it would be useful to have a method to identify genes whose expression varies stochastically on a certain time scale. Here, we developed such a method and used it for identifying genes with high levels of phenotypic noise in Salmonella enterica ssp. I serovar Typhimurium (S. Typhimurium). We created a genomic plasmid library fused to a green fluorescent protein (GFP) reporter and subjected replicate populations harboring this library to fluctuating selection for GFP expression using fluorescent-activated cell sorting (FACS). After seven rounds of fluctuating selection, the populations were strongly enriched for promoters that showed a high amount of noise in gene expression. Our results indicate that the activity of some promoters of S. Typhimurium varies on such a short time scale that these promoters can absorb rapid fluctuations in the direction of selection, as imposed during our experiment. The genomic fragments that conferred the highest levels of phenotypic variation were promoters controlling the synthesis of flagella, which are associated with virulence and host–pathogen interactions. This confirms earlier reports that phenotypic noise may play a role in pathogenesis and indicates that these promoters have among the highest levels of noise in the S. Typhimurium genome. This approach can be applied to many other bacterial and eukaryotic systems as a simple method for identifying genes with noisy expression., PLoS Genetics, 4, ISSN:1553-7390, ISSN:1553-7404

Growth cones are guided to their final destination by intermediate targets. Here, we identify intermediate targets and signaling components acting on zebrafish habenula commissural axons. Live imaging establishes that axons pause at the medial habenula before and after crossing the roof plate. esrom mutants axons fail to advance beyond the ipsilateral medial habenula. Tsc2 function is reduced in mutant axons, indicating cell autonomous defects in signaling. Consistent with signaling properties changing outside the roof plate, EphB is surface localized on axon segments within a zone demarcated by the medial habenula. wnt4a is expressed in the medial habenula and morpholino knockdown causes loss of the commissure. Electroporation of truncated Ryk causes axons to reenter the midline after reaching the contralateral habenula. These data identify Esrom as a mediator of growth cone navigation at an intermediate target and underscore the importance of midline boundaries as signaling centers for commissure formation.