Identifying Structural Variation in Haploid Microbial Genomes from Short-Read Resequencing Data Using Breseq

Mutations that alter chromosomal structure play critical roles in evolution and disease, including in the origin of new lifestyles and pathogenic traits in microbes. Large-scale rearrangements in genomes are often mediated by recombination events involving new or existing copies of mobile genetic elements, recently duplicated genes, or other repetitive sequences. Most current software programs for predicting structural variation from short-read DNA resequencing data are intended primarily for use on human genomes. They typically disregard information in reads mapping to repeat sequences, and significant post-processing and manual examination of their output is often required to rule out false-positive predictions and precisely describe mutational events. Results: We have implemented an algorithm for identifying structural variation from DNA resequencing data as part of the breseq computational pipeline for predicting mutations in haploid microbial genomes. Our method evaluates the support for new sequence junctions present in a clonal sample from split-read alignments to a reference genome, including matches to repeat sequences. Then, it uses a statistical model of read coverage evenness to accept or reject these predictions. Finally, breseq combines predictions of new junctions and deleted chromosomal regions to output biologically relevant descriptions of mutations and their effects on genes. We demonstrate the performance of breseq on simulated Escherichia coli genomes with deletions generating unique breakpoint sequences, new insertions of mobile genetic elements, and deletions mediated by mobile elements. Then, we reanalyze data from an E. coli K-12 mutation accumulation evolution experiment in which structural variation was not previously identified. Transposon insertions and large-scale chromosomal changes detected by breseq account for similar to 25% of spontaneous mutations in this strain. In all cases, we find that breseq is able to reliably predict structural variation with modest read-depth coverage of the reference genome (>40-fold). Conclusions: Using breseq to predict structural variation should be useful for studies of microbial epidemiology, experimental evolution, synthetic biology, and genetics when a reference genome for a closely related strain is available. In these cases, breseq can discover mutations that may be responsible for important or unintended changes in genomes that might otherwise go undetected.U.S. National Institutes of Health R00-GM087550U.S. National Science Foundation (NSF) DEB-0515729NSF BEACON Center for the Study of Evolution in Action DBI-0939454Cancer Prevention & Research Institute of Texas (CPRIT) RP130124University of Texas at Austin startup fundsUniversity of Texas at AustinCPRIT Cancer Research TraineeshipMolecular Bioscience

The genetic architecture of low-temperature adaptation in the wine yeast Saccharomyces cerevisiae

[Background] Low-temperature growth and fermentation of wine yeast can enhance wine aroma and make them highly desirable traits for the industry. Elucidating response to cold in Saccharomyces cerevisiae is, therefore, of paramount importance to select or genetically improve new wine strains. As most enological traits of industrial importance in yeasts, adaptation to low temperature is a polygenic trait regulated by many interacting loci.[Results] In order to unravel the genetic determinants of low-temperature fermentation, we mapped quantitative trait loci (QTLs) by bulk segregant analyses in the F13 offspring of two Saccharomyces cerevisiae industrial strains with divergent performance at low temperature. We detected four genomic regions involved in the adaptation at low temperature, three of them located in the subtelomeric regions (chromosomes XIII, XV and XVI) and one in the chromosome XIV. The QTL analysis revealed that subtelomeric regions play a key role in defining individual variation, which emphasizes the importance of these regions’ adaptive nature.[Conclusions] The reciprocal hemizygosity analysis (RHA), run to validate the genes involved in low-temperature fermentation, showed that genetic variation in mitochondrial proteins, maintenance of correct asymmetry and distribution of phospholipid in the plasma membrane are key determinants of low-temperature adaptation.This work has been financially supported from the Spanish Government through MINECO and FEDER funds (AGL2013-47300-C3-3-R and PCIN-2015-143 grants) and from Generalitat Valenciana through PROMETEOII/2014/042 grant, awarded to JMG. This study has been carried out in the context of the European Project ERA-IB “YeastTempTation” EGR thanks the Spanish government for an FPI grant BES-2011-044498 and MM also thanks the Generalitat Valenciana for a VALi+d ACIF/2015/194 grant. We acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).Peer reviewe

Chance and necessity in the genome evolution of endosymbiotic bacteria of insects

[EN] An open question in evolutionary biology is how does the selection¿drift balance determine the fates of biological interactions. We searched for signatures of selection and drift in genomes of five endosymbiotic bacterial groups known to evolve under strong genetic drift. Although most genes in endosymbiotic bacteria showed evidence of relaxed purifying selection, many genes in these bacteria exhibited stronger selective constraints than their orthologs in free-living bacterial relatives. Remarkably, most of these highly constrained genes had no role in the host¿symbiont interactions but were involved in either buffering the deleterious consequences of drift or other host-unrelated functions, suggesting that they have either acquired new roles or their role became more central in endosymbiotic bacteria. Experimental evolution of Escherichia coli under strong genetic drift revealed remarkable similarities in the mutational spectrum, genome reduction patterns and gene losses to endosymbiotic bacteria of insects. Interestingly, the transcriptome of the experimentally evolved lines showed a generalized deregulation of the genome that affected genes encoding proteins involved in mutational buffering, regulation and amino acid biosynthesis, patterns identical to those found in endosymbiotic bacteria. Our results indicate that drift has shaped endosymbiotic associations through a change in the functional landscape of bacterial genes and that the host had only a small role in such a shiftThis work was supported by Science Foundation Ireland (12/IP/1637) and grants from the Spanish Ministerio de Economia y Competitividad (MINECO-FEDER; BFU2012-36346 and BFU2015-66073-P) to MAF. DAP and CT were supported by Juan de la Cierva fellowships from MINECO (references: JCI-2011-11089 and JCA-2012-14056, respectively). DAP is supported by funds from the University of Nevada, Reno, NV, USA.Sabater-Muñoz, B.; Toft, C.; Alvarez-Ponce, D.; Fares Riaño, MA. (2017). Chance and necessity in the genome evolution of endosymbiotic bacteria of insects. The ISME Journal. 11(6):1291-1304. https://doi.org/10.1038/ismej.2017.18S12911304116Aguilar-Rodriguez J, Sabater-Munoz B, Montagud-Martinez R, Berlanga V, Alvarez-Ponce D, Wagner A et al. (2016). The molecular chaperone DnaK is a source of mutational robustness. Genome Biol Evol 8: 2979–2991.Alvarez-Ponce D, Sabater-Munoz B, Toft C, Ruiz-Gonzalez MX, Fares MA . (2016). Essentiality is a strong determinant of protein rates of evolution during mutation accumulation experiments in Escherichia coli. Genome Biol Evol 8: 2914–2927.Anders S, Huber W . (2010). Differential expression analysis for sequence count data. Genome Biol 11: R106.Archibald J . (2014) One Plus One Equals One: Symbiosis and the Evolution of Complex Life. Oxford University Press: Oxford, UK.Aussel L, Loiseau L, Hajj Chehade M, Pocachard B, Fontecave M, Pierrel F et al. (2014). ubiJ, a new gene required for aerobic growth and proliferation in macrophage, is involved in coenzyme Q biosynthesis in Escherichia coli and Salmonella enterica serovar Typhimurium. J Bacteriol 196: 70–79.Baumann P, Baumann L, Clark MA . (1996). Levels of Buchnera aphidicola chaperonin groEL during growth of the aphid Schizaphis graminum. Curr Microbiol 32: 7.Benjamini Y, Yekutieli Y . (2005). False discovery rate controlling confidence intervals for selected parameters. J Am Stat Assoc 100: 10.Bennett GM, Moran NA . (2015). Heritable symbiosis: the advantages and perils of an evolutionary rabbit hole. Proc Natl Acad Sci USA 112: 10169–10176.Bermingham J, Rabatel A, Calevro F, Vinuelas J, Febvay G, Charles H et al. (2009). Impact of host developmental age on the transcriptome of the symbiotic bacterium Buchnera aphidicola in the pea aphid (Acyrthosiphon pisum. Appl Environ Microbiol 75: 7294–7297.Bogumil D, Dagan T . (2010). Chaperonin-dependent accelerated substitution rates in prokaryotes. Genome Biol Evol 2: 602–608.Carbon S, Ireland A, Mungall CJ, Shu S, Marshall B, Lewis S et al. (2009). AmiGO: online access to ontology and annotation data. Bioinformatics 25: 288–289.Chen Z, Wang Y, Li Y, Li Y, Fu N, Ye J et al. (2012). Esre: a novel essential non-coding RNA in Escherichia coli. FEBS Lett 586: 1195–1200.Clark JW, Hossain S, Burnside CA, Kambhampati S . (2001). Coevolution between a cockroach and its bacterial endosymbiont: a biogeographical perspective. Proc Biol Sci 268: 393–398.Dale C, Wang B, Moran N, Ochman H . (2003). Loss of DNA recombinational repair enzymes in the initial stages of genome degeneration. Mol Biol Evol 20: 1188–1194.Deatherage DE, Barrick JE . (2014). Identification of mutations in laboratory-evolved microbes from next-generation sequencing data using breseq. Methods Mol Biol 1151: 165–188.Douglas AE . (2003). The nutritional physiology of aphids. Adv Insect Physiol 31: 68.Fares MA, Barrio E, Sabater-Munoz B, Moya A . (2002a). The evolution of the heat-shock protein GroEL from Buchnera, the primary endosymbiont of aphids, is governed by positive selection. Mol Biol Evol 19: 1162–1170.Fares MA, Ruiz-Gonzalez MX, Moya A, Elena SF, Barrio E . (2002b). Endosymbiotic bacteria: groEL buffers against deleterious mutations. Nature 417: 398.Gancedo C, Flores CL, Gancedo JM . (2016). The expanding landscape of moonlighting proteins in yeasts. Microbiol Mol Biol Rev 80: 765–777.Gerardo NM, Altincicek B, Anselme C, Atamian H, Barribeau SM, de Vos M et al. (2010). Immunity and other defenses in pea aphids, Acyrthosiphon pisum. Genome Biol 11: R21.Gomez-Valero L, Latorre A, Silva FJ . (2004). The evolutionary fate of nonfunctional DNA in the bacterial endosymbiont Buchnera aphidicola. Mol Biol Evol 21: 2172–2181.Gomez-Valero L, Silva FJ, Christophe Simon J, Latorre A . (2007). Genome reduction of the aphid endosymbiont Buchnera aphidicola in a recent evolutionary time scale. Gene 389: 87–95.Gonzalez-Domenech CM, Belda E, Patino-Navarrete R, Moya A, Pereto J, Latorre A . (2012). Metabolic stasis in an ancient symbiosis: genome-scale metabolic networks from two Blattabacterium cuenoti strains, primary endosymbionts of cockroaches. BMC Microbiol 12 (Suppl 1): S5.Hansen AK, Moran NA . (2011). Aphid genome expression reveals host-symbiont cooperation in the production of amino acids. Proc Natl Acad Sci USA 108: 2849–2854.Hansen AK, Moran NA . (2014). The impact of microbial symbionts on host plant utilization by herbivorous insects. Mol Ecol 23: 1473–1496.Henderson B, Fares MA, Lund PA . (2013). Chaperonin 60: a paradoxical, evolutionarily conserved protein family with multiple moonlighting functions. Biol Rev Camb Philos Soc 88: 955–987.Humphreys NJ, Douglas AE . (1997). Partitioning of symbiotic bacteria between generations of an insect: a quantitative study of a Buchnera sp. in the pea aphid (Acyrthosiphon pisum reared at different temperatures. Appl Environ Microbiol 63: 3294–3296.International Aphid Genomics Consortium. (2010). Genome sequence of the pea aphid Acyrthosiphon pisum. PLoS Biol 8: e1000313.Kadibalban AS, Bogumil D, Landan G, Dagan T . (2016). DnaK-dependent accelerated evolutionary rate in prokaryotes. Genome Biol Evol 8: 1590–1599.Katoh K, Standley DM . (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol 30: 772–780.Kelkar YD, Ochman H . (2013). Genome reduction promotes increase in protein functional complexity in bacteria. Genetics 193: 303–307.Koga R, Meng XY, Tsuchida T, Fukatsu T . (2012). Cellular mechanism for selective vertical transmission of an obligate insect symbiont at the bacteriocyte-embryo interface. Proc Natl Acad Sci USA 109: E1230–E1237.Kuo CH, Moran NA, Ochman H . (2009). The consequences of genetic drift for bacterial genome complexity. Genome Res 19: 1450–1454.Kuo CH, Ochman H . (2009). Deletional bias across the three domains of life. Genome Biol Evol 1: 145–152.Law R, Lewis DH . (1983). Biotic environments and the maintenance of sex-some evidence from mutualistic symbioses. Biol J Linnean Soc 20: 28.Liu XD, Xie L, Wei Y, Zhou X, Jia B, Liu J et al. (2014). Abiotic stress resistance, a novel moonlighting function of ribosomal protein RPL44 in the halophilic fungus Aspergillus glaucus. Appl Environ Microbiol 80: 4294–4300.Lohse M, Bolger AM, Nagel A, Fernie AR, Lunn JE, Stitt M et al. (2012). RobiNA: a user-friendly, integrated software solution for RNA-Seq-based transcriptomics. Nucleic Acids Res 40: W622–W627.Macdonald SJ, Lin GG, Russell CW, Thomas GH, Douglas AE . (2012). The central role of the host cell in symbiotic nitrogen metabolism. Proc Biol Sci 279: 2965–2973.McClure R, Balasubramanian D, Sun Y, Bobrovskyy M, Sumby P, Genco CA et al. (2013). Computational analysis of bacterial RNA-Seq data. Nucleic Acids Res 41: e140.McCutcheon JP, Moran NA . (2012). Extreme genome reduction in symbiotic bacteria. Nat Rev Microbiol 10: 13–26.McFall-Ngai M, Hadfield MG, Bosch TC, Carey HV, Domazet-Loso T, Douglas AE et al. (2013). Animals in a bacterial world, a new imperative for the life sciences. Proc Natl Acad Sci USA 110: 3229–3236.Mira A, Ochman H, Moran NA . (2001). Deletional bias and the evolution of bacterial genomes. Trends Genet 17: 589–596.Moran NA . (1996). Accelerated evolution and Muller's rachet in endosymbiotic bacteria. Proc Natl Acad Sci USA 93: 2873–2878.Moran NA, Dunbar HE, Wilcox JL . (2005). Regulation of transcription in a reduced bacterial genome: nutrient-provisioning genes of the obligate symbiont Buchnera aphidicola. J Bacteriol 187: 4229–4237.Moran NA, McCutcheon JP, Nakabachi A . (2008). Genomics and evolution of heritable bacterial symbionts. Annu Rev Genet 42: 165–190.Moran NA, McLaughlin HJ, Sorek R . (2009). The dynamics and time scale of ongoing genomic erosion in symbiotic bacteria. Science 323: 379–382.Nakabachi A, Ishida K, Hongoh Y, Ohkuma M, Miyagishima SY . (2014). Aphid gene of bacterial origin encodes a protein transported to an obligate endosymbiont. Curr Biol 24: R640–R641.Nilsson AI, Koskiniemi S, Eriksson S, Kugelberg E, Hinton JC, Andersson DI . (2005). Bacterial genome size reduction by experimental evolution. Proc Natl Acad Sci USA 102: 12112–12116.Patino-Navarrete R, Moya A, Latorre A, Pereto J . (2013). Comparative genomics of Blattabacterium cuenoti: the frozen legacy of an ancient endosymbiont genome. Genome Biol Evol 5: 351–361.Pettersson ME, Berg OG . (2007). Muller's ratchet in symbiont populations. Genetica 130: 199–211.Price DR, Feng H, Baker JD, Bavan S, Luetje CW, Wilson AC . (2014). Aphid amino acid transporter regulates glutamine supply to intracellular bacterial symbionts. Proc Natl Acad Sci USA 111: 320–325.Reyes-Prieto M, Vargas-Chavez C, Latorre A, Moya A . (2015). SymbioGenomesDB: a database for the integration and access to knowledge on host-symbiont relationships. Database 2015: bav109 (1–8).Robinson MD, McCarthy DJ, Smyth GK . (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26: 139–140.Sabater-Muñoz B, Prats-Escriche M, Montagud-Martinez R, Lopez-Cerdan A, Toft C, Aguilar-Rodriguez J et al. (2015). Fitness trade-offs determine the role of the molecular chaperonin groel in buffering mutations. Mol Biol Evol 32: 2681–2693.Schlicker A, Domingues FS, Rahnenfuhrer J, Lengauer T . (2006). A new measure for functional similarity of gene products based on Gene Ontology. BMC Bioinformatics 7: 302.Shigenobu S, Watanabe H, Hattori M, Sakaki Y, Ishikawa H . (2000). Genome sequence of the endocellular bacterial symbiont of aphids Buchnera sp. APS. Nature 407: 81–86.Supek F, Bosnjak M, Skunca N, Smuc T . (2011). REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS ONE 6: e21800.Tamas I, Klasson L, Canback B, Naslund AK, Eriksson AS, Wernegreen JJ et al. (2002). 50 million years of genomic stasis in endosymbiotic bacteria. Science 296: 2376–2379.Toft C, Fares MA . (2008). The evolution of the flagellar assembly pathway in endosymbiotic bacterial genomes. Mol Biol Evol 25: 2069–2076.van Ham RC, Kamerbeek J, Palacios C, Rausell C, Abascal F, Bastolla U et al. (2003). Reductive genome evolution in Buchnera aphidicola. Proc Natl Acad Sci USA 100: 581–586.Wernegreen JJ . (2002). Genome evolution in bacterial endosymbionts of insects. Nat Rev Genet 3: 850–861.Wernegreen JJ . (2011). Reduced selective constraint in endosymbionts: elevation in radical amino acid replacements occurs genome-wide. PLoS One 6: e28905.Williams TA, Fares MA . (2010). The effect of chaperonin buffering on protein evolution. Genome Biol Evol 2: 609–619.Yang Z . (2007). PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol 24: 1586–1591

The vertebrate muscle Z-disc: sarcomere anchor for structure and signalling

The Z-disc, appearing as a fine dense line forming sarcomere boundaries in striated muscles, when studied in detail reveals crosslinked filament arrays that transmit tension and house myriads of proteins with diverse functions. At the Z-disc the barbed ends of the antiparallel actin filaments from adjoining sarcomeres interdigitate and are crosslinked primarily by layers of α-actinin. The Z-disc is therefore the site of polarity reversal of the actin filaments, as needed to interact with the bipolar myosin filaments in successive sarcomeres. The layers of α-actinin determine the Z-disc width: fast fibres have narrow (~30–50 nm) Z-discs and slow and cardiac fibres have wide (~100 nm) Z-discs. Comprehensive reviews on the roles of the numerous proteins located at the Z-disc in signalling and disease have been published; the aim here is different, namely to review the advances in structural aspects of the Z-disc

Extracellular vesicles are independent metabolic units with asparaginase activity.

Extracellular vesicles (EVs) are membrane particles involved in the exchange of a broad range of bioactive molecules between cells and the microenvironment. Although it has been shown that cells can traffic metabolic enzymes via EVs, much remains to be elucidated with regard to their intrinsic metabolic activity. Accordingly, herein we assessed the ability of neural stem/progenitor cell (NSC)-derived EVs to consume and produce metabolites. Our metabolomics and functional analyses both revealed that EVs harbor L-asparaginase activity, catalyzed by the enzyme asparaginase-like protein 1 (Asrgl1). Critically, we show that Asrgl1 activity is selective for asparagine and is devoid of glutaminase activity. We found that mouse and human NSC EVs traffic Asrgl1. Our results demonstrate, for the first time, that NSC EVs function as independent metabolic units that are able to modify the concentrations of critical nutrients, with the potential to affect the physiology of their microenvironment.This work has received support from the Italian Multiple Sclerosis Association (AISM, grant 2010/R/31 and grant 2014/PMS/4 to SP), the Italian Ministry of Health (GR08-7 to SP), the European Research Council (ERC) under the ERC-2010-StG Grant agreement n° 260511-SEM_SEM, the Medical Research Council, the Engineering and Physical Sciences Research Council, and the Biotechnology and Biological Sciences Research Council UK Regenerative Medicine Platform Hub “Acellular Approaches for Therapeutic Delivery” (MR/K026682/1 to SP), The Evelyn Trust (RG 69865 to SP), The Bascule Charitable Trust (RG 75149 to SP) and core support grant from the Wellcome Trust and Medical Research Council to the Wellcome Trust – MRC Cambridge Stem Cell Institute. N.I. was supported by a FEBS long-term fellowship. C.F., A.S.H., and E.G. were funded by the Medical Research Council, Core Fund SKAG006

Critique of the review of 'Water fluoridation for the prevention of dental caries' published by the Cochrane Collaboration in 2015

The Cochrane Review on water fluoridation for the prevention of dental caries was published in 2015 and attracted considerable interest and comment, especially in countries with extensive water fluoridation programmes. The Review had two objectives: (i) to evaluate the effects of water fluoridation (artificial or natural) on the prevention of dental caries, and (ii) to evaluate the effects of water fluoridation (artificial or natural) on dental fluorosis. The authors concluded, inter alia, that there was very little contemporary evidence, meeting the Review's inclusion criteria, that evaluated the effectiveness of water fluoridation for the prevention of dental caries. The purpose of this critique is to examine the conduct of the above Review, and to put it into context in the wider body of evidence regarding the effectiveness of water fluoridation. While the overall conclusion that water fluoridation is effective in caries prevention agrees with previous reviews, many important public health questions could not be answered by the Review because of the restrictive criteria used to judge adequacy of study design and risk of bias. The potential benefits of using wider criteria in order to achieve a fuller understanding of the effectiveness of water fluoridation are discussed