RAREFAN: A webservice to identify REPINs and RAYTs in bacterial genomes

Compared to eukaryotes, repetitive sequences are rare in bacterial genomes and usually do not persist for long. Yet, there is at least one class of persistent prokaryotic mobile genetic elements: REPINs. REPINs are non-autonomous transposable elements replicated by single-copy transposases called RAYTs. REPIN-RAYT systems are mostly vertically inherited and have persisted in individual bacterial lineages for millions of years. Discovering and analyzing REPIN populations and their corresponding RAYT transposases in bacterial species can be rather laborious, hampering progress in understanding REPIN-RAYT biology and evolution. Here we present RAREFAN, a webservice that identifies REPIN populations and their corresponding RAYT transposase in a given set of bacterial genomes. We demonstrate RAREFAN’s capabilities by analyzing a set of 49 Stenotrophomonas maltophilia genomes, containing nine different REPIN-RAYT systems. We guide the reader through the process of identifying and analyzing REPIN-RAYT systems across S. maltophilia, highlighting erroneous associations between REPIN and RAYTs, and providing solutions on how to find correct associations. RAREFAN enables rapid, large-scale detection of REPINs and RAYTs, and provides insight into the fascinating world of intragenomic sequence populations in bacterial genomes. RAREFAN is available at http://rarefan.evolbio.mpg.de

Bistability in a Metabolic Network Underpins the De Novo Evolution of Colony Switching in Pseudomonas fluorescens

Phenotype switching is commonly observed in nature. This prevalence has allowed the elucidation of a number of underlying molecular mechanisms. However, little is known about how phenotypic switches arise and function in their early evolutionary stages. The first opportunity to provide empirical insight was delivered by an experiment in which populations of the bacterium Pseudomonas fluorescens SBW25 evolved, de novo, the ability to switch between two colony phenotypes. Here we unravel the molecular mechanism behind colony switching, revealing how a single nucleotide change in a gene enmeshed in central metabolism (carB) generates such a striking phenotype. We show that colony switching is underpinned by ON/OFF expression of capsules consisting of a colanic acid-like polymer. We use molecular genetics, biochemical analyses, and experimental evolution to establish that capsule switching results from perturbation of the pyrimidine biosynthetic pathway. Of central importance is a bifurcation point at which uracil triphosphate is partitioned towards either nucleotide metabolism or polymer production. This bifurcation marks a cell-fate decision point whereby cells with relatively high pyrimidine levels favour nucleotide metabolism (capsule OFF), while cells with lower pyrimidine levels divert resources towards polymer biosynthesis (capsule ON). This decision point is present and functional in the wild-type strain. Finally, we present a simple mathematical model demonstrating that the molecular components of the decision point are capable of producing switching. Despite its simple mutational cause, the connection between genotype and phenotype is complex and multidimensional, offering a rare glimpse of how noise in regulatory networks can provide opportunity for evolution

Dissecting HIV Virulence: Heritability of Setpoint Viral Load, CD4+ T-Cell Decline, and Per-Parasite Pathogenicity.

Pathogen strains may differ in virulence because they attain different loads in their hosts, or because they induce different disease-causing mechanisms independent of their load. In evolutionary ecology, the latter is referred to as "per-parasite pathogenicity". Using viral load and CD4+ T-cell measures from 2014 HIV-1 subtype B-infected individuals enrolled in the Swiss HIV Cohort Study, we investigated if virulence-measured as the rate of decline of CD4+ T cells-and per-parasite pathogenicity are heritable from donor to recipient. We estimated heritability by donor-recipient regressions applied to 196 previously identified transmission pairs, and by phylogenetic mixed models applied to a phylogenetic tree inferred from HIV pol sequences. Regressing the CD4+ T-cell declines and per-parasite pathogenicities of the transmission pairs did not yield heritability estimates significantly different from zero. With the phylogenetic mixed model, however, our best estimate for the heritability of the CD4+ T-cell decline is 17% (5-30%), and that of the per-parasite pathogenicity is 17% (4-29%). Further, we confirm that the set-point viral load is heritable, and estimate a heritability of 29% (12-46%). Interestingly, the pattern of evolution of all these traits differs significantly from neutrality, and is most consistent with stabilizing selection for the set-point viral load, and with directional selection for the CD4+ T-cell decline and the per-parasite pathogenicity. Our analysis shows that the viral genotype affects virulence mainly by modulating the per-parasite pathogenicity, while the indirect effect via the set-point viral load is minor

Within-Genome Evolution of REPINs: a New Family of Miniature Mobile DNA in Bacteria

Repetitive sequences are a conserved feature of many bacterial genomes. While first reported almost thirty years ago, and frequently exploited for genotyping purposes, little is known about their origin, maintenance, or processes affecting the dynamics of within-genome evolution. Here, beginning with analysis of the diversity and abundance of short oligonucleotide sequences in the genome of Pseudomonas fluorescens SBW25, we show that over-represented short sequences define three distinct groups (GI, GII, and GIII) of repetitive extragenic palindromic (REP) sequences. Patterns of REP distribution suggest that closely linked REP sequences form a functional replicative unit: REP doublets are over-represented, randomly distributed in extragenic space, and more highly conserved than singlets. In addition, doublets are organized as inverted repeats, which together with intervening spacer sequences are predicted to form hairpin structures in ssDNA or mRNA. We refer to these newly defined entities as REPINs (REP doublets forming hairpins) and identify short reads from population sequencing that reveal putative transposition intermediates. The proximal relationship between GI, GII, and GIII REPINs and specific REP-associated tyrosine transposases (RAYTs), combined with features of the putative transposition intermediate, suggests a mechanism for within-genome dissemination. Analysis of the distribution of REPs in a range of RAYT–containing bacterial genomes, including Escherichia coli K-12 and Nostoc punctiforme, show that REPINs are a widely distributed, but hitherto unrecognized, family of miniature non-autonomous mobile DNA

The evolution of selfish genetic elements within bacterial genomes : a thesis submitted in partial fulfilment of the requirements for the degree of Ph.D. in Molecular Evolution at Massey University, Auckland, New Zealand

Content removed due to copyright restrictions: Bertels, F., & Rainey, P.B. (2011). Curiosities of REPINs and RAYTs. Mobile Genetic Elements, 1(4), 1-7. doi.:10.4161/mge.1.4.18610.Genes that increase their copy number relative to that of the host genome are termed selfish. Selfish genes are found ubiquitously in bacterial genomes. Within genomes they can often be identified due to their repetitive nature. Short repetitive sequences such as repetitive extragenic palindromic (REP) sequences have been proposed to be selfish genetic elements. However, evidence for the selfishness of REPs is scarce due to the lack of knowledge about their origin, evolution and mechanisms of dispersal. Here, REPs are studied in the model bacterium Pseudomonas fluorescens SBW25. The evidence provided suggests that REPs are part of a greater mobile genetic element, which is termed REP doublet forming hairpins (REPINs). Subsequently, I investigate the cause of REPIN dispersal: a putative transposase. The transposase, named REP-associated tyrosine transposase (RAYT) shares essential motifs with the IS200 family of insertion sequences. However, unlike insertion sequences, RAYTs are found only as single copy genes. This indicates that RAYTs may not be entirely selfish; instead they may have been co-opted by the host to perform a beneficial function. Finally, two more repetitive sequence classes are studied in the SBW25 genome. Interestingly, both sequence classes consist of a protein coding sequence and a sequence that forms a stable secondary structure in single stranded DNA or RNA. This arrangement is reminiscent of bacterial toxin-antitoxin (TA) systems. Evidence from sequence analyses suggests that the repetitive nature of these elements in SBW25 may be the result of cooperation between REPINs or other replicative elements and the TA systems. The presented analyses show that despite the streamlined nature of bacterial genomes selfish genetic elements frequently arise, replicate and probably increase theirpersistence and spread through cooperation with addictive and duplicative elements respectively. persistence and spread through cooperation with addictive and duplicative elements respectively

Simulating bacterial evolution forward-in-time

A recommendation – based on reviews by three anonymous reviewers – of the article: Cury J, Haller BC, Achaz G, Jay F (2021) Simulation of bacterial populations with SLiM. bioRxiv, 2020.09.28.316869, version 5 peer-reviewed and recommended by Peer community in Evolutionary Biology. https://doi.org/10.1101/2020.09.28.31686

Discovering Complete Quasispecies in Bacterial Genomes

Mobile genetic elements can be found in almost all genomes. Possibly the most common nonautonomous mobile genetic elements in bacteria are repetitive extragenic palindromic doublets forming hairpins (REPINs) that can occur hundreds of times within a genome. The sum of all REPINs in a genome can be viewed as an evolving population because REPINs replicate and mutate. In contrast to most other biological populations, we know the exact composition of the REPIN population and the sequence of each member of the population. Here, we model the evolution of REPINs as quasispecies. We fit our quasispecies model to 10 different REPIN populations from 10 different bacterial strains and estimate effective duplication rates. Our estimated duplication rates range from ~5 3 1029 to 15 3 1029 duplications per bacterial generation per REPIN. The small range and the low level of the REPIN duplication rates suggest a universal trade-off between the survival of the REPIN population and the reduction of the mutational load for the host genome. The REPIN populations we investigated also possess features typical of other natural populations. One population shows hallmarks of a population that is going extinct, another population seems to be growing in size, and we also see an example of competition between two REPIN populations. © 2017 by the Genetics Society of America

Estimating the mutational fitness effects distribution during early HIV infection

The evolution of HIV during acute infection is often considered a neutral process. Recent analysis of sequencing data from this stage of infection, however, showed high levels of shared mutations between independent viral populations. This suggests that selection might play a role in the early stages of HIV infection. We adapted an existing model for random evolution during acute HIV-infection to include selection. Simulations of this model were used to fit a global mutational fitness effects distribution to previously published sequencing data of the env gene of individuals with acute HIV infection. Measures of sharing between viral populations were used as summary statistics to compare the data to the simulations. We confirm that evolution during acute infection is significantly different from neutral. The distribution of mutational fitness effects is best fit by a distribution with a low, but significant fraction of beneficial mutations and a high fraction of deleterious mutations. While most mutations are neutral or deleterious in this model, about 5% of mutations are beneficial. These beneficial mutations will, on average, result in a small but significant increase in fitness. When assuming no epistasis, this indicates that, at the moment of transmission, HIV is near, but not on the fitness peak for early infection