Pseudogenes in plasmid genomes reveal past transitions in plasmid mobility.

Evidence for gene non-functionalization due to mutational processes is found in genomes in the form of pseudogenes. Pseudogenes are known to be rare in prokaryote chromosomes, with the exception of lineages that underwent an extreme genome reduction (e.g. obligatory symbionts). Much less is known about the frequency of pseudogenes in prokaryotic plasmids; those are genetic elements that can transfer between cells and may encode beneficial traits for their host. Non-functionalization of plasmid-encoded genes may alter the plasmid characteristics, e.g. mobility, or their effect on the host. Analyzing 10 832 prokaryotic genomes, we find that plasmid genomes are characterized by threefold-higher pseudogene density compared to chromosomes. The majority of plasmid pseudogenes correspond to deteriorated transposable elements. A detailed analysis of enterobacterial plasmids furthermore reveals frequent gene non-functionalization events associated with the loss of plasmid self-transmissibility. Reconstructing the evolution of closely related plasmids reveals that non-functionalization of the conjugation machinery led to the emergence of non-mobilizable plasmid types. Examples are virulence plasmids in Escherichia and Salmonella. Our study highlights non-functionalization of core plasmid mobility functions as one route for the evolution of domesticated plasmids. Pseudogenes in plasmids supply insights into past transitions in plasmid mobility that are akin to transitions in bacterial lifestyle.

Segregational drift hinders the evolution of antibiotic resistance on polyploid replicons.

The emergence of antibiotic resistance under treatment depends on the availability of resistance alleles and their establishment in the population. Novel resistance alleles are encoded either in chromosomal or extrachromosomal genetic elements; both types may be present in multiple copies within the cell. However, the effect of polyploidy on the emergence of antibiotic resistance remains understudied. Here we show that the establishment of resistance alleles in microbial populations depends on the ploidy level. Evolving bacterial populations under selection for antibiotic resistance, we demonstrate that resistance alleles in polyploid elements are lost frequently in comparison to alleles in monoploid elements due to segregational drift. Integrating the experiments with a mathematical model, we find a remarkable agreement between the theoretical and empirical results, confirming our understanding of the allele segregation process. Using the mathematical model, we further show that the effect of polyploidy on the establishment probability of beneficial alleles is strongest for low replicon copy numbers and plateaus for high replicon copy numbers. Our results suggest that the distribution of fitness effects for mutations that are eventually fixed in a population depends on the replicon ploidy level. Our study indicates that the emergence of antibiotic resistance in bacterial pathogens depends on the pathogen ploidy level.

Heteroplasmy is rare in plant mitochondria compared to plastids despite similar mutation rates.

Plant cells harbor two membrane-bound organelles containing their own genetic material -plastids and mitochondria. Although the two organelles co-exist and co-evolve within the same plant cells, they differ in genome copy number, intracellular organization, and mode of segregation. How these attributes affect the time to fixation, or conversely, loss of neutral alleles is currently unresolved. Here we show that mitochondria and plastids share the same mutation rate yet plastid alleles remain in a heteroplasmic state significantly longer compared to mitochondrial alleles. By analyzing genetic variants across populations of the marine flowering plant Zostera marina and simulating organelle allele dynamics, we examine the determinants of allele segregation and allele fixation. Our results suggest that bottlenecks on the cell population, e.g., during branching or seeding, and stratification of the meristematic tissue, are important determinants of mitochondrial allele dynamics. Furthermore, we suggest that the prolonged plastid allele dynamics are due to a yet unknown active plastid partition mechanism. The dissimilarity between plastid and mitochondrial novel allele fixation at different levels of organization may manifest in differences in adaptation processes. Our study uncovers fundamental principles of organelle population genetics that are essential for further investigations of long-term evolution and molecular dating of divergence events.

Denitrification in foraminifera has an ancient origin and is complemented by associated bacteria.

Benthic foraminifera are unicellular eukaryotes that inhabit sediments of aquatic environments. Several foraminifera of the order Rotaliida are known to store and use nitrate for denitrification, a unique energy metabolism among eukaryotes. The rotaliid Globobulimina spp. has been shown to encode an incomplete denitrification pathway of bacterial origin. However, the prevalence of denitrification genes in foraminifera remains unknown, and the missing denitrification pathway components are elusive. Analyzing transcriptomes and metagenomes of 10 foraminiferal species from the Peruvian oxygen minimum zone, we show that denitrification genes are highly conserved in foraminifera. We infer the last common ancestor of denitrifying foraminifera, which enables us to predict the ability to denitrify for additional foraminiferal species. Additionally, an examination of the foraminiferal microbiota reveals evidence for a stable interaction with Desulfobacteraceae, which harbor genes that complement the foraminiferal denitrification pathway. Our results provide evidence that foraminiferal denitrification is complemented by the foraminifera-associated microbiome. The interaction of foraminifera with their resident bacteria is at the basis of foraminiferal adaptation to anaerobic environments that manifested in ecological success in oxygen depleted habitats.

Essential gene acquisition destabilizes plasmid inheritance.

Extra-chromosomal genetic elements are important drivers of evolutionary transformations and ecological adaptations in prokaryotes with their evolutionary success often depending on their 'utility' to the host. Examples are plasmids encoding antibiotic resistance genes, which are known to proliferate in the presence of antibiotics. Plasmids carrying an essential host function are recognized as permanent residents in their host. Essential plasmids have been reported in several taxa where they often encode essential metabolic functions; nonetheless, their evolution remains poorly understood. Here we show that essential genes are rarely encoded on plasmids; evolving essential plasmids in Escherichia coli we further find that acquisition of an essential chromosomal gene by a plasmid can lead to plasmid extinction. A comparative genomics analysis of Escherichia isolates reveals few plasmid-encoded essential genes, yet these are often integrated into plasmid-related functions; an example is the GroEL/GroES chaperonin. Experimental evolution of a chaperonin-encoding plasmid shows that the acquisition of an essential gene reduces plasmid fitness regardless of the stability of plasmid inheritance. Our results suggest that essential plasmid emergence leads to a dose effect caused by gene redundancy. The detrimental effect of essential gene acquisition on plasmid inheritance constitutes a barrier for plasmid-mediated lateral gene transfer and supplies a mechanistic understanding for the rarity of essential genes in extra-chromosomal genetic elements.

Segregational Drift Constrains the Evolutionary Rate of Prokaryotic Plasmids.

Plasmids are extrachromosomal genetic elements in prokaryotes that have been recognized as important drivers of microbial ecology and evolution. Plasmids are found in multiple copies inside their host cell where independent emergence of mutations may lead to intracellular genetic heterogeneity. The intracellular plasmid diversity is thus subject to changes upon cell division. However, the effect of plasmid segregation on plasmid evolution remains understudied. Here, we show that genetic drift during cell division-segregational drift-leads to the rapid extinction of novel plasmid alleles. We established a novel experimental approach to control plasmid allele frequency at the levels of a single cell and the whole population. Following the dynamics of plasmid alleles in an evolution experiment, we find that the mode of plasmid inheritance-random or clustered-is an important determinant of plasmid allele dynamics. Phylogenetic reconstruction of our model plasmid in clinical isolates furthermore reveals a slow evolutionary rate of plasmid-encoded genes in comparison to chromosomal genes. Our study provides empirical evidence that genetic drift in plasmid evolution occurs at multiple levels: the host cell and the population of hosts. Segregational drift has implications for the evolutionary rate heterogeneity of extrachromosomal genetic elements.

Insertion and deletion evolution reflects antibiotics selection pressure in a Mycobacterium tuberculosis outbreak.

In genome evolution, genetic variants are the source of diversity, which natural selection acts upon. Treatment of human tuberculosis (TB) induces a strong selection pressure for the emergence of antibiotic resistance-conferring variants in the infecting Mycobacterium tuberculosis (MTB) strains. MTB evolution in response to treatment has been intensively studied and mainly attributed to point substitutions. However, the frequency and contribution of insertions and deletions (indels) to MTB genome evolution remains poorly understood. Here, we analyzed a multi-drug resistant MTB outbreak for the presence of high-quality indels and substitutions. We find that indels are significantly enriched in genes conferring antibiotic resistance. Furthermore, we show that indels are inherited during the outbreak and follow a molecular clock with an evolutionary rate of 5.37e-9 indels/site/year, which is 23 times lower than the substitution rate. Inherited indels may co-occur with substitutions in genes along related biological pathways; examples are iron storage and resistance to second-line antibiotics. This suggests that epistatic interactions between indels and substitutions affect antibiotic resistance and compensatory evolution in MTB.

Metabolic preference of nitrate over oxygen as an electron acceptor in foraminifera from the Peruvian oxygen minimum zone.

Benthic foraminifera populate a diverse range of marine habitats. Their ability to use alternative electron acceptors-nitrate (NO3-) or oxygen (O2)-makes them important mediators of benthic nitrogen cycling. Nevertheless, the metabolic scaling of the two alternative respiration pathways and the environmental determinants of foraminiferal denitrification rates are yet unknown. We measured denitrification and O2 respiration rates for 10 benthic foraminifer species sampled in the Peruvian oxygen minimum zone (OMZ). Denitrification and O2 respiration rates significantly scale sublinearly with the cell volume. The scaling is lower for O2 respiration than for denitrification, indicating that NO3- metabolism during denitrification is more efficient than O2 metabolism during aerobic respiration in foraminifera from the Peruvian OMZ. The negative correlation of the O2 respiration rate with the surface/volume ratio is steeper than for the denitrification rate. This is likely explained by the presence of an intracellular NO3- storage in denitrifying foraminifera. Furthermore, we observe an increasing mean cell volume of the Peruvian foraminifera, under higher NO3- availability. This suggests that the cell size of denitrifying foraminifera is not limited by O2 but rather by NO3- availability. Based on our findings, we develop a mathematical formulation of foraminiferal cell volume as a predictor of respiration and denitrification rates, which can further constrain foraminiferal biogeochemical cycling in biogeochemical models. Our findings show that NO3- is the preferred electron acceptor in foraminifera from the OMZ, where the foraminiferal contribution to denitrification is governed by the ratio between NO3- and O2.

A novel septal protein of multicellular heterocystous cyanobacteria is associated with the divisome.

Cyanobacteria are unique among the eubacteria as they possess a hybrid Gram phenotype, having an outer membrane but also a comparably thick peptidoglycan sheet. Furthermore, the cyanobacterial divisome includes proteins specific for both the Gram types as well as cyanobacteria-specific proteins. Cells in multicellular cyanobacteria share a continuous periplasm and their cytoplasms are connected by septal junctions that enable communication between cells in the filament. The localization of septal junction proteins depends on interaction with the divisome, however additional yet unknown proteins may be involved in this process. Here, we characterized Alr3364 (termed SepI), a novel septal protein that interacts with the divisome in the multicellular heterocystous cyanobacterium Anabaena sp. strain PCC 7120. SepI localized to the Z-ring and the intercellular septa but did not interact with FtsZ. Instead, SepI interacted with the divisome proteins ZipN, SepF and FtsI and with the septal protein SepJ. The inactivation of sepI led to a defect in cell filament integrity, colony and cell morphology, septum size, nanopore formation and peptidoglycan biogenesis, and inability to differentiate heterocysts. Our results show that SepI plays a role in intercellular communication and furthermore indicate that SepI functions in the coordination of septal junction localization during cell division.

Colonization dynamics of Pantoea agglomerans in the wheat root habitat.

Plants are colonized by microbial communities that have diverse implications for plant development and health. The establishment of a stable plant-bacteria interaction depends on a continuous coexistence over generations. Transmission via the seed is considered as the main route for vertical inheritance of plant-associated bacteria. Nonetheless, the ecological principles that govern the plant colonization by seed endophytes remain understudied. Here we quantify the contribution of arrival time and colonization history to bacterial colonization of the wheat root. Establishing a common seed endophyte, Pantoea agglomerans, and wheat as a model system enabled us to document bacterial colonization of the plant roots during the early stages of germination. Using our system, we estimate the carrying capacity of the wheat roots as 108 cells g-1 , which is robust among individual plants and over time. Competitions in planta reveal a significant advantage of early incoming colonizers over late-incoming colonizers. Priming for the wheat environment had little effect on the colonizer success. Our experiments thus provide empirical data on the root colonization dynamics of a seed endophyte. The persistence of seed endophyte bacteria with the plant population over generations may contribute to the stable transmission that is one route for the evolution of a stable host-associated lifestyle.

Discovery of multi-operon colinear syntenic blocks in microbial genomes.

MOTIVATION: An important task in comparative genomics is to detect functional units by analyzing gene-context patterns. Colinear syntenic blocks (CSBs) are groups of genes that are consistently encoded in the same neighborhood and in the same order across a wide range of taxa. Such CSBs are likely essential for the regulation of gene expression in prokaryotes. Recent results indicate that colinearity can be conserved across multiple operons, thus motivating the discovery of multi-operon CSBs. This computational task raises scalability challenges in large datasets. RESULTS: We propose an efficient algorithm for the discovery of cross-strand multi-operon CSBs in large genomic datasets. The proposed algorithm uses match-point arithmetic, which is scalable for large datasets of microbial genomes in terms of running time and space requirements. The algorithm is implemented and incorporated into a tool with a graphical user interface, called CSBFinder-S. We applied CSBFinder-S to data mine 1485 prokaryotic genomes and analyzed the identified cross-strand CSBs. Our results indicate that most of the syntenic blocks are exclusively colinear. Additional results indicate that transcriptional regulation by overlapping transcriptional genes is abundant in bacteria. We demonstrate the utility of CSBFinder-S to identify common function of the gene-pair PulEF in multiple contexts, including Type 2 Secretion System, Type 4 Pilus System and DNA uptake machinery. AVAILABILITY AND IMPLEMENTATION: CSBFinder-S software and code are publicly available at https://github.com/dinasv/CSBFinder. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Segregational Drift and the Interplay between Plasmid Copy Number and Evolvability.

The ubiquity of plasmids in all prokaryotic phyla and habitats and their ability to transfer between cells marks them as prominent constituents of prokaryotic genomes. Many plasmids are found in their host cell in multiple copies. This leads to an increased mutational supply of plasmid-encoded genes and genetically heterogeneous plasmid genomes. Nonetheless, the segregation of plasmid copies into daughter cells during cell division is considered to occur in the absence of selection on the plasmid alleles. We investigate the implications of random genetic drift of multicopy plasmids during cell division-termed here "segregational drift"-to plasmid evolution. Performing experimental evolution of low- and high-copy non-mobile plasmids in Escherichia coli, we find that the evolutionary rate of multicopy plasmids does not reflect the increased mutational supply expected according to their copy number. In addition, simulated evolution of multicopy plasmid alleles demonstrates that segregational drift leads to increased loss frequency and extended fixation time of plasmid mutations in comparison to haploid chromosomes. Furthermore, an examination of the experimentally evolved hosts reveals a significant impact of the plasmid type on the host chromosome evolution. Our study demonstrates that segregational drift of multicopy plasmids interferes with the retention and fixation of novel plasmid variants. Depending on the selection pressure on newly emerging variants, plasmid genomes may evolve slower than haploid chromosomes, regardless of their higher mutational supply. We suggest that plasmid copy number is an important determinant of plasmid evolvability due to the manifestation of segregational drift.

CSBFinder: discovery of colinear syntenic blocks across thousands of prokaryotic genomes.

MOTIVATION: Identification of conserved syntenic blocks across microbial genomes is important for several problems in comparative genomics such as gene annotation, study of genome organization and evolution and prediction of gene interactions. Current tools for syntenic block discovery do not scale up to the large quantity of prokaryotic genomes available today. RESULTS: We present a novel methodology for the discovery, ranking and taxonomic distribution analysis of colinear syntenic blocks (CSBs)-groups of genes that are consistently located close to each other, in the same order, across a wide range of taxa. We present an efficient algorithm that identifies CSBs in large genomic datasets. The algorithm is implemented and incorporated in a novel tool with a graphical user interface, denoted CSBFinder, that ranks the discovered CSBs according to a probabilistic score and clusters them to families according to their gene content similarity. We apply CSBFinder to data mine 1487 prokaryotic genomes including chromosomes and plasmids. For post-processing analysis, we generate heatmaps for visualizing the distribution of CSB family members across various taxa. We exemplify the utility of CSBFinder in operon prediction, in deciphering unknown gene function and in taxonomic analysis of colinear syntenic blocks. AVAILABILITY AND IMPLEMENTATION: CSBFinder software and code are publicly available at https://github.com/dinasv/CSBFinder. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Rates of Mutation and Recombination in Siphoviridae Phage Genome Evolution over Three Decades.

The evolution of asexual organisms is driven not only by the inheritance of genetic modification but also by the acquisition of foreign DNA. The contribution of vertical and horizontal processes to genome evolution depends on their rates per year and is quantified by the ratio of recombination to mutation. These rates have been estimated for bacteria; however, no estimates have been reported for phages. Here, we delineate the contribution of mutation and recombination to dsDNA phage genome evolution. We analyzed 34 isolates of the 936 group of Siphoviridae phages using a Lactococcus lactis strain from a single dairy over 29 years. We estimate a constant substitution rate of 1.9 × 10-4 substitutions per site per year due to mutation that is within the range of estimates for eukaryotic RNA and DNA viruses. The reconstruction of recombination events reveals a constant rate of five recombination events per year and 4.5 × 10-3 nucleotide alterations due to recombination per site per year. Thus, the recombination rate exceeds the substitution rate, resulting in a relative effect of recombination to mutation (r/m) of ∼24 that is homogenous over time. Especially in the early transcriptional region, we detect frequent gene loss and regain due to recombination with phages of the 936 group, demonstrating the role of the 936 group pangenome as a reservoir of genetic variation. The observed substitution rate homogeneity conforms to the neutral theory of evolution; hence, the neutral theory can be applied to phage genome evolution and also to genetic variation brought about by recombination.

Plasticity first: molecular signatures of a complex morphological trait in filamentous cyanobacteria.

BACKGROUND: Filamentous cyanobacteria that differentiate multiple cell types are considered the peak of prokaryotic complexity and their evolution has been studied in the context of multicellularity origins. Species that form true-branching filaments exemplify the most complex cyanobacteria. However, the mechanisms underlying the true-branching morphology remain poorly understood despite of several investigations that focused on the identification of novel genes or pathways. An alternative route for the evolution of novel traits is based on existing phenotypic plasticity. According to that scenario - termed genetic assimilation - the fixation of a novel phenotype precedes the fixation of the genotype. RESULTS: Here we show that the evolution of transcriptional regulatory elements constitutes a major mechanism for the evolution of new traits. We found that supplementation with sucrose reconstitutes the ancestral branchless phenotype of two true-branching Fischerella species and compared the transcription start sites (TSSs) between the two phenotypic states. Our analysis uncovers several orthologous TSSs whose transcription level is correlated with the true-branching phenotype. These TSSs are found in genes that encode components of the septosome and elongasome (e.g., fraC and mreB). CONCLUSIONS: The concept of genetic assimilation supplies a tenable explanation for the evolution of novel traits but testing its feasibility is hindered by the inability to recreate and study the evolution of present-day traits. We present a novel approach to examine transcription data for the plasticity first route and provide evidence for its occurrence during the evolution of complex colony morphology in true-branching cyanobacteria. Our results reveal a route for evolution of the true-branching phenotype in cyanobacteria via modification of the transcription level of pre-existing genes. Our study supplies evidence for the 'plasticity-first' hypothesis and highlights the importance of transcriptional regulation in the evolution of novel traits.

Integration of two ancestral chaperone systems into one: the evolution of eukaryotic molecular chaperones in light of eukaryogenesis.

Eukaryotic genomes are mosaics of genes acquired from their prokaryotic ancestors, the eubacterial endosymbiont that gave rise to the mitochondrion and its archaebacterial host. Genomic footprints of the prokaryotic merger at the origin of eukaryotes are still discernable in eukaryotic genomes, where gene expression and function correlate with their prokaryotic ancestry. Molecular chaperones are essential in all domains of life as they assist the functional folding of their substrate proteins and protect the cell against the cytotoxic effects of protein misfolding. Eubacteria and archaebacteria code for slightly different chaperones, comprising distinct protein folding pathways. Here we study the evolution of the eukaryotic protein folding pathways following the endosymbiosis event. A phylogenetic analysis of all 64 chaperones encoded in the Saccharomyces cerevisiae genome revealed 25 chaperones of eubacterial ancestry, 11 of archaebacterial ancestry, 10 of ambiguous prokaryotic ancestry, and 18 that may represent eukaryotic innovations. Several chaperone families (e.g., Hsp90 and Prefoldin) trace their ancestry to only one prokaryote group, while others, such as Hsp40 and Hsp70, are of mixed ancestry, with members contributed from both prokaryotic ancestors. Analysis of the yeast chaperone-substrate interaction network revealed no preference for interaction between chaperones and substrates of the same origin. Our results suggest that the archaebacterial and eubacterial protein folding pathways have been reorganized and integrated into the present eukaryotic pathway. The highly integrated chaperone system of yeast is a manifestation of the central role of chaperone-mediated folding in maintaining cellular fitness. Most likely, both archaebacterial and eubacterial chaperone systems were essential at the very early stages of eukaryogenesis, and the retention of both may have offered new opportunities for expanding the scope of chaperone-mediated folding.

Acquisition of 1,000 eubacterial genes physiologically transformed a methanogen at the origin of Haloarchaea.

Archaebacterial halophiles (Haloarchaea) are oxygen-respiring heterotrophs that derive from methanogens--strictly anaerobic, hydrogen-dependent autotrophs. Haloarchaeal genomes are known to have acquired, via lateral gene transfer (LGT), several genes from eubacteria, but it is yet unknown how many genes the Haloarchaea acquired in total and, more importantly, whether independent haloarchaeal lineages acquired their genes in parallel, or as a single acquisition at the origin of the group. Here we have studied 10 haloarchaeal and 1,143 reference genomes and have identified 1,089 haloarchaeal gene families that were acquired by a methanogenic recipient from eubacteria. The data suggest that these genes were acquired in the haloarchaeal common ancestor, not in parallel in independent haloarchaeal lineages, nor in the common ancestor of haloarchaeans and methanosarcinales. The 1,089 acquisitions include genes for catabolic carbon metabolism, membrane transporters, menaquinone biosynthesis, and complexes I-IV of the eubacterial respiratory chain that functions in the haloarchaeal membrane consisting of diphytanyl isoprene ether lipids. LGT on a massive scale transformed a strictly anaerobic, chemolithoautotrophic methanogen into the heterotrophic, oxygen-respiring, and bacteriorhodopsin-photosynthetic haloarchaeal common ancestor.

Directed networks reveal genomic barriers and DNA repair bypasses to lateral gene transfer among prokaryotes.

Lateral gene transfer (LGT) plays a major role in prokaryote evolution with only a few genes that are resistant to it; yet the nature and magnitude of barriers to lateral transfer are still debated. Here, we implement directed networks to investigate donor-recipient events of recent lateral gene transfer among 657 sequenced prokaryote genomes. For 2,129,548 genes investigated, we detected 446,854 recent lateral gene transfer events through nucleotide pattern analysis. Among these, donor-recipient relationships could be specified through phylogenetic reconstruction for 7% of the pairs, yielding 32,028 polarized recent gene acquisition events, which constitute the edges of our directed networks. We find that the frequency of recent LGT is linearly correlated both with genome sequence similarity and with proteome similarity of donor-recipient pairs. Genome sequence similarity accounts for 25% of the variation in gene-transfer frequency, with proteome similarity adding only 1% to the variability explained. The range of donor-recipient GC content similarity within the network is extremely narrow, with 86% of the LGTs occurring between donor-recipient pairs having ≤5% difference in GC content. Hence, genome sequence similarity and GC content similarity are strong barriers to LGT in prokaryotes. But they are not insurmountable, as we detected 1530 recent transfers between distantly related genomes. The directed network revealed that recipient genomes of distant transfers encode proteins of nonhomologous end-joining (NHEJ; a DNA repair mechanism) far more frequently than the recipient lacking that mechanism. This implicates NHEJ in genes spread across distantly related prokaryotes through bypassing the donor-recipient sequence similarity barrier.

The genome of the obligate intracellular parasite Trachipleistophora hominis: new insights into microsporidian genome dynamics and reductive evolution.

The dynamics of reductive genome evolution for eukaryotes living inside other eukaryotic cells are poorly understood compared to well-studied model systems involving obligate intracellular bacteria. Here we present 8.5 Mb of sequence from the genome of the microsporidian Trachipleistophora hominis, isolated from an HIV/AIDS patient, which is an outgroup to the smaller compacted-genome species that primarily inform ideas of evolutionary mode for these enormously successful obligate intracellular parasites. Our data provide detailed information on the gene content, genome architecture and intergenic regions of a larger microsporidian genome, while comparative analyses allowed us to infer genomic features and metabolism of the common ancestor of the species investigated. Gene length reduction and massive loss of metabolic capacity in the common ancestor was accompanied by the evolution of novel microsporidian-specific protein families, whose conservation among microsporidians, against a background of reductive evolution, suggests they may have important functions in their parasitic lifestyle. The ancestor had already lost many metabolic pathways but retained glycolysis and the pentose phosphate pathway to provide cytosolic ATP and reduced coenzymes, and it had a minimal mitochondrion (mitosome) making Fe-S clusters but not ATP. It possessed bacterial-like nucleotide transport proteins as a key innovation for stealing host-generated ATP, the machinery for RNAi, key elements of the early secretory pathway, canonical eukaryotic as well as microsporidian-specific regulatory elements, a diversity of repetitive and transposable elements, and relatively low average gene density. Microsporidian genome evolution thus appears to have proceeded in at least two major steps: an ancestral remodelling of the proteome upon transition to intracellular parasitism that involved reduction but also selective expansion, followed by a secondary compaction of genome architecture in some, but not all, lineages.

The evolution of antibiotic resistance islands occurs within the framework of plasmid lineages.

Bacterial pathogens carrying multidrug resistance (MDR) plasmids are a major threat to human health. The acquisition of antibiotic resistance genes (ARGs) in plasmids is often facilitated by mobile genetic elements that copy or translocate ARGs between DNA molecules. The agglomeration of mobile elements in plasmids generates resistance islands comprising multiple ARGs. However, whether the emergence of resistance islands is restricted to specific MDR plasmid lineages remains understudied. Here we show that the agglomeration of ARGs in resistance islands is biased towards specific large plasmid lineages. Analyzing 6784 plasmids in 2441 Escherichia, Salmonella, and Klebsiella isolates, we quantify that 84% of the ARGs in MDR plasmids are found in resistance islands. We furthermore observe rapid evolution of ARG combinations in resistance islands. Most regions identified as resistance islands are shared among closely related plasmids but rarely among distantly related plasmids. Our results suggest the presence of barriers for the dissemination of ARGs between plasmid lineages, which are related to plasmid genetic properties, host range and the plasmid evolutionary history. The agglomeration of ARGs in plasmids is attributed to the workings of mobile genetic elements that operate within the framework of existing plasmid lineages.