Role of natural transformation in the evolution of small cryptic plasmids in Synechocystis sp. PCC 6803

Small cryptic plasmids have no clear effect on the host fitness and their functional repertoire remains obscure. The naturally competent cyanobacterium Synechocystis sp. PCC 6803 harbours several small cryptic plasmids; whether their evolution with this species is supported by horizontal transfer remains understudied. Here, we show that the small cryptic plasmid DNA is transferred in the population exclusively by natural transformation, where the transfer frequency of plasmid-encoded genes is similar to that of chromosome-encoded genes. Establishing a system to follow gene transfer, we compared the transfer frequency of genes encoded in cryptic plasmids pCA2.4 (2378 bp) and pCB2.4 (2345 bp) within and between populations of two Synechocystis sp. PCC 6803 labtypes (termed Kiel and Sevilla). Our results reveal that plasmid gene transfer frequency depends on the recipient labtype. Furthermore, gene transfer via whole plasmid uptake in the Sevilla labtype ranged among the lowest detected transfer rates in our experiments. Our study indicates that horizontal DNA transfer via natural transformation is frequent in the evolution of small cryptic plasmids that reside in naturally competent organisms. Furthermore, we suggest that the contribution of natural transformation to cryptic plasmid persistence in Synechocystis is limited.

SepT, a novel protein specific to multicellular cyanobacteria, influences peptidoglycan growth and septal nanopore formation in Anabaena sp. PCC 7120

IMPORTANCE

Multicellular organization is a requirement for the development of complex organisms, and filamentous cyanobacteria such as Anabaena represent a paradigmatic case of bacterial multicellularity. The Anabaena filament can include hundreds of communicated cells that exchange nutrients and regulators and, depending on environmental conditions, can include different cell types specialized in distinct biological functions. Hence, the specific features of the Anabaena filament and how they are propagated during cell division represent outstanding biological issues. Here, we studied SepT, a novel coiled-coil-rich protein of Anabaena that is located in the intercellular septa and influences the formation of the septal specialized structures that allow communication between neighboring cells along the filament, a fundamental trait for the performance of Anabaena as a multicellular organism.

Worldwide Population Genomics Reveal Long-Term Stability of the Mitochondrial Genome Architecture in a Keystone Marine Plant

Mitochondrial genomes (mitogenomes) of flowering plants are composed of multiple chromosomes. Recombination within and between the mitochondrial chromosomes may generate diverse DNA molecules termed isoforms. The isoform copy number and composition can be dynamic within and among individual plants due to uneven replication and homologous recombination. Nonetheless, despite their functional importance, the level of mitogenome conservation within species remains understudied. Whether the ontogenetic variation translates to evolution of mitogenome composition over generations is currently unknown. Here we show that the mitogenome composition of the seagrass Zostera marina is conserved among worldwide populations that diverged ca. 350,000 years ago. Using long-read sequencing, we characterized the Z. marina mitochondrial genome and inferred the repertoire of recombination-induced configurations. To characterize the mitochondrial genome architecture worldwide and study its evolution, we examined the mitogenome in Z. marina meristematic region sampled in 16 populations from the Pacific and Atlantic oceans. Our results reveal a striking similarity in the isoform relative copy number, indicating a high conservation of the mitogenome composition among distantly related populations and within the plant germline, despite a notable variability during individual ontogenesis. Our study supplies a link between observations of dynamic mitogenomes at the level of plant individuals and long-term mitochondrial evolution.

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.

Ocean current patterns drive the worldwide colonization of eelgrass (Zostera marina)

Currents are unique drivers of oceanic phylogeography and thus determine the distribution of marine coastal species, along with past glaciations and sea-level changes. Here we reconstruct the worldwide colonization history of eelgrass (Zostera marina L.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Two west-to-east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival of Z. marina in the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present there for ~243 ky (thousand years). Mediterranean populations were founded ~44 kya, while extant distributions along western and eastern Atlantic shores were founded at the end of the Last Glacial Maximum (~19 kya), with at least one major refuge being the North Carolina region. The recent colonization and five- to sevenfold lower genomic diversity of the Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.

Phylogenomic Testing of Root Hypotheses

The determination of the last common ancestor (LCA) of a group of species plays a vital role in evolutionary theory. Traditionally, an LCA is inferred by the rooting of a fully resolved species tree. From a theoretical perspective, however, inference of the LCA amounts to the reconstruction of just one branch - the root branch - of the true species tree, and should therefore be a much easier task than the full resolution of the species tree. Discarding the reliance on a hypothesised species tree and its rooting leads us to re-evaluate what phylogenetic signal is directly relevant to LCA inference, and to recast the task as that of sampling the total evidence from all gene families at the genomic scope. Here we reformulate LCA and root inference in the framework of statistical hypothesis testing and outline an analytical procedure to formally test competing a-priori LCA hypotheses and to infer confidence sets for the earliest speciation events in the history of a group of species. Applying our methods to two demonstrative datasets we show that our inference of the opisthokonta LCA is well in agreement with the common knowledge. Inference of the proteobacteria LCA shows that it is most closely related to modern Epsilonproteobacteria, raising the possibility that it may have been characterized by a chemolithoautotrophic and anaerobic life-style. Our inference is based on data comprising between 43% (opisthokonta) and 86% (proteobacteria) of all gene families. Approaching LCA inference within a statistical framework renders the phylogenomic inference powerful and robust.

Out of the 'host' box: extreme off-host conditions alter the infectivity and virulence of a parasitic bacterium

Disease agents play an important role in the ecology and life history of wild and cultivated populations and communities. While most studies focus on the adaptation of parasites to their hosts, the adaptation of free-living parasite stages to their external (off-host) environment may tell us a lot about the factors that shape the distribution of parasites. Pasteuria ramosa is an endoparasitic bacterium of the water flea Daphnia with a wide geographical distribution. Its transmission stages rest outside of the host and thus experience varying environmental regimes. We examined the life history of P. ramosa populations from four environmental conditions (i.e. groups of habitats): the factorial combinations of summer-dry water bodies or not, and winter-freeze water bodies or not. Our goal was to examine how the combination of winter temperature and summer dryness affects the parasite's ability to attach to its host and to infect it. We subjected samples of the four groups of habitats to temperatures of 20, 33, 46 and 60°C in dry and wet conditions, and exposed a susceptible clone of Daphnia magna to the treated spores. We found that spores which had undergone desiccation endured higher temperatures better than spores kept wet, both regarding attachment and subsequent infection. Furthermore, spores treated with heightened temperatures were much less infective and virulent. Even under high temperatures (60°C), exposed spores from all populations were able to attach to the host cuticle, albeit they were unable to establish infection. Our work highlights the sensitivity of a host-free resting stage of a bacterial parasite to the external environment. Long heatwaves and harsh summers, which are becoming more frequent owing to recent climate changes, may therefore pose a problem for parasite survival. This article is part of the theme issue 'Infectious disease ecology and evolution in a changing world'.

Assessing Functional Cognition and Health During COVID-19 Pandemic: Gender Differences Among Community-Dwelling Adults

INTRODUCTION/OBJECTIVES

The COVID-19 pandemic has long-term implications for adult health and function, whether or not people were infected with the disease. Although cognitive disruptions are among the major symptoms of COVID-19, most research focused on managing medical symptoms, such as respiratory symptoms or pain. Thus, less is known about the pandemic's long-term implications for assessing functional cognition. This study aimed to examine COVID-19's effects on community-dwelling adults' functional cognition and health, comparing gender differences.

METHODS

This cross-sectional study divided 118 community-dwelling adults (25 previously infected with COVID-19) into gender groups. Primary outcome measures included the Daily Living Questionnaire (DLQ) and short form health status survey, SF-12.

RESULTS

No significant differences were found in functional cognition or health between participants who had contracted COVID-19 and those who remained healthy, but men had better functional cognition and health measures in comparison with women.

CONCLUSIONS

Gender differences in functional cognition and health state may relate to gender-based family roles. It is essential to assess functional cognition of young adults who were exposed to a pandemic, such as COVID-19, because it may significantly affect their health and functional status. The DLQ is a reliable, valid assessment of functional cognition that may suit individuals who previously contracted COVID-19.

Fixation dynamics of beneficial alleles in prokaryotic polyploid chromosomes and plasmids

Theoretical population genetics has been mostly developed for sexually reproducing diploid and for monoploid (haploid) organisms, focusing on eukaryotes. The evolution of bacteria and archaea is often studied by models for the allele dynamics in monoploid populations. However, many prokaryotic organisms harbor multicopy replicons—chromosomes and plasmids—and theory for the allele dynamics in populations of polyploid prokaryotes remains lacking. Here, we present a population genetics model for replicons with multiple copies in the cell. Using this model, we characterize the fixation process of a dominant beneficial mutation at 2 levels: the phenotype and the genotype. Our results show that depending on the mode of replication and segregation, the fixation of the mutant phenotype may precede genotypic fixation by many generations; we term this time interval the heterozygosity window. We furthermore derive concise analytical expressions for the occurrence and length of the heterozygosity window, showing that it emerges if the copy number is high and selection strong. Within the heterozygosity window, the population is phenotypically adapted, while both alleles persist in the population. Replicon ploidy thus allows for the maintenance of genetic variation following phenotypic adaptation and consequently for reversibility in adaptation to fluctuating environmental conditions.

Pangenome evolution in environmentally transmitted symbionts of deep-sea mussels is governed by vertical inheritance

Microbial pangenomes vary across species; their size and structure are determined by genetic diversity within the population and by gene loss and horizontal gene transfer (HGT). Many bacteria are associated with eukaryotic hosts where the host colonization dynamics may impact bacterial genome evolution. Host-associated lifestyle has been recognized as a barrier to HGT in parentally transmitted bacteria. However, pangenome evolution of environmentally acquired symbionts remains understudied, often due to limitations in symbiont cultivation. Using high-resolution metagenomics, here we study pangenome evolution of two co-occurring endosymbionts inhabiting Bathymodiolus brooksi mussels from a single cold seep. The symbionts, sulfur-oxidizing (SOX) and methane-oxidizing (MOX) gamma-proteobacteria, are environmentally acquired at an early developmental stage and individual mussels may harbor multiple strains of each symbiont species. We found differences in the accessory gene content of both symbionts across individual mussels, which are reflected by differences in symbiont strain composition. Compared to core genes, accessory genes are enriched in genome plasticity functions. We found no evidence for recent horizontal gene transfer between both symbionts. A comparison between the symbiont pangenomes revealed that the MOX population is less diverged and contains fewer accessory genes, supporting that the MOX association with B. brooksi is more recent in comparison to that of SOX. Our results show that the pangenomes of both symbionts evolved mainly by vertical inheritance. We conclude that genome evolution of environmentally transmitted symbionts that associate with individual hosts over their lifetime is affected by a narrow symbiosis where the frequency of HGT is constrained.​.

Gene sharing among plasmids and chromosomes reveals barriers for antibiotic resistance gene transfer

The emergence of antibiotic resistant bacteria is a major threat to modern medicine. Rapid adaptation to antibiotics is often mediated by the acquisition of plasmids carrying antibiotic resistance (ABR) genes. Nonetheless, the determinants of plasmid-mediated ABR gene transfer remain debated. Here, we show that the propensity of ABR gene transfer via plasmids is higher for accessory chromosomal ABR genes in comparison with core chromosomal ABR genes, regardless of the resistance mechanism. Analysing the pattern of ABR gene occurrence in the genomes of 2635 Enterobacteriaceae isolates, we find that 33% of the 416 ABR genes are shared between chromosomes and plasmids. Phylogenetic reconstruction of ABR genes occurring on both plasmids and chromosomes supports their evolution by lateral gene transfer. Furthermore, accessory ABR genes (encoded in less than 10% of the chromosomes) occur more abundantly in plasmids in comparison with core ABR genes (encoded in greater than or equal to 90% of the chromosomes). The pattern of ABR gene occurrence in plasmids and chromosomes is similar to that in the total Escherichia genome. Our results thus indicate that the previously recognized barriers for gene acquisition by lateral gene transfer apply also to ABR genes. We propose that the functional complexity of the underlying ABR mechanism is an important determinant of ABR gene transferability. This article is part of the theme issue 'The secret lives of microbial mobile 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.

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.

Mobile genetic elements have been extensively studied due to their role as agents of genetic innovation and rapid adaptation in prokaryotes. Specifically, prokaryotic plasmids have been the focus of investigation in the context of bacterial survival under growth limiting conditions with the prime example of resistance to antibiotics and heavy metals. In contrast, plasmids that encode for functions that are essential to their host viability are rarely described. We investigate the evolution of plasmids that encode for genes previously identified as essential for bacterial life. Our analysis of Escherichia isolates reveals only few plasmid-encoded essential genes, which likely function in the plasmid rather than the host life cycle. Following the evolution of plasmids encoding an essential gene in Escherichia coli in real time, we further find that the acquisition of a chromosomal essential gene may lead to plasmid loss. Our study supplies data and a mechanistic understanding on the rarity of essential genes in mobile genetic elements. We conclude that prokaryotic plasmids are rarely essential for their bacterial host.

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 10⁸ 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.

Darwinian individuality of extrachromosomal genetic elements calls for population genetics tinkering

Plasmids are extrachromosomal genetic elements commonly found in prokaryotic (and sometime eukaryotic) cells. Small plasmids are often considered cryptic and their effect on the host is elusive, while large plasmids may encode functions that are essential for the host lifestyle and attain a secondary chromosome status. Plasmids are thus an important source of raw material for microbial genome evolution outside the mainstream of bacterial chromosomes. The discovery of plasmid-mediated antibiotic resistance led to extensive research on the contribution of plasmids to the environmental dimensions of antibiotic resistance and the evolution of plasmid-host interactions following the acquisition of plasmids encoding for antibiotic resistance. Recent experimental studies revealed the importance of intracellular plasmid diversity for plasmid-host interactions. Here we describe the evolutionary forces at play during plasmid evolution in a top-down approach: this includes the effect of processes at the level of the host population and the consideration of plasmids as Darwinian individuals within the host cell.

Two novel heteropolymer-forming proteins maintain the multicellular shape of the cyanobacterium Anabaena sp. PCC 7120

Polymerizing and filament-forming proteins are instrumental for numerous cellular processes such as cell division and growth. Their function in stabilization and localization of protein complexes and replicons is achieved by a filamentous structure. Known filamentous proteins assemble into homopolymers consisting of single subunits - for example, MreB and FtsZ in bacteria - or heteropolymers that are composed of two subunits, for example, keratin and α/β tubulin in eukaryotes. Here, we describe two novel coiled-coil-rich proteins (CCRPs) in the filament-forming cyanobacterium Anabaena sp. PCC 7120 (hereafter Anabaena) that assemble into a heteropolymer and function in the maintenance of the Anabaena multicellular shape (termed trichome). The two CCRPs - Alr4504 and Alr4505 (named ZicK and ZacK) - are strictly interdependent for the assembly of protein filaments in vivo and polymerize nucleotide independently in vitro, similar to known intermediate filament (IF) proteins. A ΔzicKΔzacK double mutant is characterized by a zigzagged cell arrangement and hence a loss of the typical linear Anabaena trichome shape. ZicK and ZacK interact with themselves, with each other, with the elongasome protein MreB, the septal junction protein SepJ and the divisome associate septal protein SepI. Our results suggest that ZicK and ZacK function in cooperation with SepJ and MreB to stabilize the Anabaena trichome and are likely essential for the manifestation of the multicellular shape in Anabaena. Our study reveals the presence of filament-forming IF-like proteins whose function is achieved through the formation of heteropolymers in cyanobacteria.

The Order of Trait Emergence in the Evolution of Cyanobacterial Multicellularity

The transition from unicellular to multicellular organisms is one of the most significant events in the history of life. Key to this process is the emergence of Darwinian individuality at the higher level: Groups must become single entities capable of reproduction for selection to shape their evolution. Evolutionary transitions in individuality are characterized by cooperation between the lower level entities and by division of labor. Theory suggests that division of labor may drive the transition to multicellularity by eliminating the trade off between two incompatible processes that cannot be performed simultaneously in one cell. Here, we examine the evolution of the most ancient multicellular transition known today, that of cyanobacteria, where we reconstruct the sequence of ecological and phenotypic trait evolution. Our results show that the prime driver of multicellularity in cyanobacteria was the expansion in metabolic capacity offered by nitrogen fixation, which was accompanied by the emergence of the filamentous morphology and succeeded by a reproductive life cycle. This was followed by the progression of multicellularity into higher complexity in the form of differentiated cells and patterned multicellularity.

Plasmid evolution

Plasmids are genetic elements that colonize and replicate in prokaryotic cells (Box 1). They are considered a major driving force of prokaryote evolution, as they can migrate between populations, making them potent agents of lateral DNA transfer and microbial warfare. The importance of plasmids goes beyond microbial evolution, as they are widely used as vectors for genetic engineering in basic research (e.g., random mutagenesis) as well as applications in biotechnology (e.g., insulin production), synthetic biology, agriculture (e.g., genetic engineering of crops) and medicine (e.g., biopharmaceuticals).

Intracellular Competitions Reveal Determinants of Plasmid Evolutionary Success

Plasmids are autonomously replicating genetic elements that are ubiquitous in all taxa and habitats where they constitute an integral part of microbial genomes. The stable inheritance of plasmids depends on their segregation during cell division and their long-term persistence in a host population is thought to largely depend on their impact on the host fitness. Nonetheless, many plasmids found in nature are lacking a clear trait that is advantageous to their host; the determinants of plasmid evolutionary success in the absence of plasmid benefit to the host remain understudied. Here we show that stable plasmid inheritance is an important determinant of plasmid evolutionary success. Borrowing terminology from evolutionary biology of cellular living forms, we hypothesize that Darwinian fitness is key for the plasmid evolutionary success. Performing intracellular plasmid competitions between non-mobile plasmids enables us to compare the evolutionary success of plasmid genotypes within the host, i.e., the plasmid fitness. Intracellular head-to-head competitions between stable and unstable variants of the same model plasmid revealed that the stable plasmid variant has a higher fitness in comparison to the unstable plasmid. Preemptive plasmid competitions reveal that plasmid fitness may depend on the order of plasmid arrival in the host. Competitions between plasmids characterized by similar stability of inheritance reveal plasmid fitness differences depending on the plasmid-encoded trait. Our results further reveal that competing plasmids can be maintained in coexistence following plasmid fusions that maintain unstable plasmid variants over time. Plasmids are not only useful accessory genetic elements to their host but they are also evolving and replicating entities, similarly to cellular living forms. There is a clear link between plasmid genetics and plasmid evolutionary success – hence plasmids are evolving entities whose fitness is quantifiable.

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.

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.

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.

Somatic genetic drift and multilevel selection in a clonal seagrass

All multicellular organisms are genetic mosaics owing to somatic mutations. The accumulation of somatic genetic variation in clonal species undergoing asexual (or clonal) reproduction may lead to phenotypic heterogeneity among autonomous modules (termed ramets). However, the abundance and dynamics of somatic genetic variation under clonal reproduction remain poorly understood. Here we show that branching events in a seagrass (Zostera marina) clone or genet lead to population bottlenecks of tissue that result in the evolution of genetically differentiated ramets in a process of somatic genetic drift. By studying inter-ramet somatic genetic variation, we uncovered thousands of single nucleotide polymorphisms that segregated among ramets. Ultra-deep resequencing of single ramets revealed that the strength of purifying selection on mosaic genetic variation was greater within than among ramets. Our study provides evidence for multiple levels of selection during the evolution of seagrass genets. Somatic genetic drift during clonal propagation leads to the emergence of genetically unique modules that constitute an elementary level of selection and individuality in long-lived clonal species.

Identification and characterization of novel filament-forming proteins in cyanobacteria

Filament-forming proteins in bacteria function in stabilization and localization of proteinaceous complexes and replicons; hence they are instrumental for myriad cellular processes such as cell division and growth. Here we present two novel filament-forming proteins in cyanobacteria. Surveying cyanobacterial genomes for coiled-coil-rich proteins (CCRPs) that are predicted as putative filament-forming proteins, we observed a higher proportion of CCRPs in filamentous cyanobacteria in comparison to unicellular cyanobacteria. Using our predictions, we identified nine protein families with putative intermediate filament (IF) properties. Polymerization assays revealed four proteins that formed polymers in vitro and three proteins that formed polymers in vivo. Fm7001 from Fischerella muscicola PCC 7414 polymerized in vitro and formed filaments in vivo in several organisms. Additionally, we identified a tetratricopeptide repeat protein - All4981 - in Anabaena sp. PCC 7120 that polymerized into filaments in vitro and in vivo. All4981 interacts with known cytoskeletal proteins and is indispensable for Anabaena viability. Although it did not form filaments in vitro, Syc2039 from Synechococcus elongatus PCC 7942 assembled into filaments in vivo and a Δsyc2039 mutant was characterized by an impaired cytokinesis. Our results expand the repertoire of known prokaryotic filament-forming CCRPs and demonstrate that cyanobacterial CCRPs are involved in cell morphology, motility, cytokinesis and colony integrity.

The Effect of Population Bottleneck Size and Selective Regime on Genetic Diversity and Evolvability in Bacteria

Population bottlenecks leading to a drastic reduction of the population size are common in the evolutionary dynamics of natural populations; their occurrence is known to have implications for genome evolution due to genetic drift, the consequent reduction in genetic diversity, and the rate of adaptation. Nevertheless, an empirical characterization of the effect of population bottleneck size on evolutionary dynamics of bacteria is currently lacking. In this study, we show that selective conditions have a stronger effect on the evolutionary history of bacteria in comparison to population bottlenecks. We evolved Escherichia coli populations under three different population bottleneck sizes (small, medium, and large) in two temperature regimes (37 °C and 20 °C). We find a high genetic diversity in the large in comparison to the small bottleneck size. Nonetheless, the cold temperature led to reduced genetic diversity regardless the bottleneck size; hence, the temperature has a stronger effect on the genetic diversity in comparison to the bottleneck size. A comparison of the fitness gain among the evolved populations reveals a similar pattern where the temperature has a significant effect on the fitness. Our study demonstrates that population bottlenecks are an important determinant of bacterial evolvability; their consequences depend on the selective conditions and are best understood via their effect on the standing genetic variation.

Quantification of Plasmid-Mediated Antibiotic Resistance in an Experimental Evolution Approach

Plasmids play a major role in microbial ecology and evolution as vehicles of lateral gene transfer and reservoirs of accessory gene functions in microbial populations. This is especially the case under rapidly changing environments such as fluctuating antibiotics exposure. We recently showed that plasmids maintain antibiotic resistance genes in Escherichia coli without positive selection for the plasmid presence. Here we describe an experimental system that allows following both the plasmid genotype and phenotype in long-term evolution experiments. We use molecular techniques to design a model plasmid that is subsequently introduced to an experimental evolution batch system approach in an E. coli host. We follow the plasmid frequency over time by applying replica plating of the E. coli populations while quantifying the antibiotic resistance persistence. In addition, we monitor the conformation of plasmids in host cells by analyzing the extent of plasmid multimer formation by plasmid nicking and agarose gel electrophoresis. Such an approach allows us to visualize not only the genome size of evolving plasmids but also their topological conformation-a factor highly important for plasmid inheritance. Our system combines molecular strategies with traditional microbiology approaches and provides a set-up to follow plasmids in bacterial populations over a long time. The presented approach can be applied to study a wide range of mobile genetic elements in the future.

Horizontally transmitted symbiont populations in deep-sea mussels are genetically isolated

Eukaryotes are habitats for bacterial organisms where the host colonization and dispersal among individual hosts have consequences for the bacterial ecology and evolution. Vertical symbiont transmission leads to geographic isolation of the microbial population and consequently to genetic isolation of microbiotas from individual hosts. In contrast, the extent of geographic and genetic isolation of horizontally transmitted microbiota is poorly characterized. Here we show that chemosynthetic symbionts of individual Bathymodiolus brooksi mussels constitute genetically isolated subpopulations. The reconstruction of core genome-wide strains from high-resolution metagenomes revealed distinct phylogenetic clades. Nucleotide diversity and strain composition vary along the mussel life span and individual hosts show a high degree of genetic isolation. Our results suggest that the uptake of environmental bacteria is a restricted process in B. brooksi, where self-infection of the gill tissue results in serial founder effects during symbiont evolution. We conclude that bacterial colonization dynamics over the host life cycle is thus an important determinant of population structure and genome evolution of horizontally transmitted symbionts.

Comparative analysis of amplicon and metagenomic sequencing methods reveals key features in the evolution of animal metaorganisms

BACKGROUND

The interplay between hosts and their associated microbiome is now recognized as a fundamental basis of the ecology, evolution, and development of both players. These interdependencies inspired a new view of multicellular organisms as "metaorganisms." The goal of the Collaborative Research Center "Origin and Function of Metaorganisms" is to understand why and how microbial communities form long-term associations with hosts from diverse taxonomic groups, ranging from sponges to humans in addition to plants.

METHODS

In order to optimize the choice of analysis procedures, which may differ according to the host organism and question at hand, we systematically compared the two main technical approaches for profiling microbial communities, 16S rRNA gene amplicon and metagenomic shotgun sequencing across our panel of ten host taxa. This includes two commonly used 16S rRNA gene regions and two amplification procedures, thus totaling five different microbial profiles per host sample.

CONCLUSION

While 16S rRNA gene-based analyses are subject to much skepticism, we demonstrate that many aspects of bacterial community characterization are consistent across methods. The resulting insight facilitates the selection of appropriate methods across a wide range of host taxa. Overall, we recommend single- over multi-step amplification procedures, and although exceptions and trade-offs exist, the V3 V4 over the V1 V2 region of the 16S rRNA gene. Finally, by contrasting taxonomic and functional profiles and performing phylogenetic analysis, we provide important and novel insight into broad evolutionary patterns among metaorganisms, whereby the transition of animals from an aquatic to a terrestrial habitat marks a major event in the evolution of host-associated microbial composition.

Interactions and Coadaptation in Plant Metaorganisms

Plants associate with a wide diversity of microorganisms. Some microorganisms engage in intimate associations with the plant host, collectively forming a metaorganism. Such close coexistence with plants requires specific adaptations that allow microorganisms to overcome plant defenses and inhabit plant tissues during growth and reproduction. New data suggest that the plant immune system has a broader role beyond pathogen recognition and also plays an important role in the community assembly of the associated microorganism. We propose that core microorganisms undergo coadaptation with their plant host, notably in response to the plant immune system allowing them to persist and propagate in their host. Microorganisms, which are vertically transmitted from generation to generation via plant seeds, putatively compose highly adapted species and may have plant-beneficial functions. The extent to which plant domestication has impacted the underlying genetics of plant-microbe associations remains poorly understood. We propose that the ability of domesticated plants to select and maintain advantageous microbial partners may have been affected. In this review, we discuss factors that impact plant metaorganism assembly and function. We underline the importance of microbe-microbe interactions in plant tissues, as they are still poorly studied but may have a great impact on plant health.

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.

Evolthon: A community endeavor to evolve lab evolution

In experimental evolution, scientists evolve organisms in the lab, typically by challenging them to new environmental conditions. How best to evolve a desired trait? Should the challenge be applied abruptly, gradually, periodically, sporadically? Should one apply chemical mutagenesis, and do strains with high innate mutation rate evolve faster? What are ideal population sizes of evolving populations? There are endless strategies, beyond those that can be exposed by individual labs. We therefore arranged a community challenge, Evolthon, in which students and scientists from different labs were asked to evolve Escherichia coli or Saccharomyces cerevisiae for an abiotic stress—low temperature. About 30 participants from around the world explored diverse environmental and genetic regimes of evolution. After a period of evolution in each lab, all strains of each species were competed with one another. In yeast, the most successful strategies were those that used mating, underscoring the importance of sex in evolution. In bacteria, the fittest strain used a strategy based on exploration of different mutation rates. Different strategies displayed variable levels of performance and stability across additional challenges and conditions. This study therefore uncovers principles of effective experimental evolutionary regimens and might prove useful also for biotechnological developments of new strains and for understanding natural strategies in evolutionary arms races between species. Evolthon constitutes a model for community-based scientific exploration that encourages creativity and cooperation.

This Community Page article describes Evolthon; a first-of-its-kind community-based effort, involving about 30 participant labs around the world, aiming to explore the best strategy for evolving microorganisms to cope with an environmental challenge.

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⁻⁴ 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⁻³ 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.

A Novel Eukaryotic Denitrification Pathway in Foraminifera

Benthic foraminifera are unicellular eukaryotes inhabiting sediments of aquatic environments. Several species were shown to store and use nitrate for complete denitrification, a unique energy metabolism among eukaryotes. The population of benthic foraminifera reaches high densities in oxygen-depleted marine habitats, where they play a key role in the marine nitrogen cycle. However, the mechanisms of denitrification in foraminifera are still unknown, and the possibility of a contribution of associated bacteria is debated. Here, we present evidence for a novel eukaryotic denitrification pathway that is encoded in foraminiferal genomes. Large-scale genome and transcriptomes analyses reveal the presence of a denitrification pathway in foraminifera species of the genus Globobulimina. This includes the enzymes nitrite reductase (NirK) and nitric oxide reductase (Nor) as well as a wide range of nitrate transporters (Nrt). A phylogenetic reconstruction of the enzymes' evolutionary history uncovers evidence for an ancient acquisition of the foraminiferal denitrification pathway from prokaryotes. We propose a model for denitrification in foraminifera, where a common electron transport chain is used for anaerobic and aerobic respiration. The evolution of hybrid respiration in foraminifera likely contributed to their ecological success, which is well documented in palaeontological records since the Cambrian period.

Recombination Signal in Mycobacterium tuberculosis Stems from Reference-guided Assemblies and Alignment Artefacts

DNA acquisition via genetic recombination is considered advantageous as it has the potential to bring together beneficial mutations that emerge independently within a population. Furthermore, recombination is considered to contribute to the maintenance of genome stability by purging slightly deleterious mutations. The prevalence of recombination differs among prokaryotic species and depends on the accessibility of DNA transfer mechanisms. An exceptional example is the human pathogen Mycobacterium tuberculosis (MTB) where no clear transfer mechanisms have been so far characterized and the presence of recombination is questioned. Here, we analyze completely assembled MTB genomes in search for evidence of recombination. We find that putative recombination events are enriched in strains reconstructed by reference-guided assembly and in regions with unreliable alignments. In addition, assembly and alignment artefacts introduce phylogenetic signals that are conflicting the established MTB phylogeny. Our results reveal that the so far reported recombination events in MTB are likely to stem from methodological artefacts. We conclude that no reliable signal of recombination is observed in the currently available MTB genomes. Moreover, our study demonstrates the limitations of reference-guided genome assembly for phylogenetic reconstructions. Rigorously de novo assembled genomes of high quality are mandatory in order to distinguish true evolutionary signal from noise, in particular for low diversity species such as MTB.

Carrying Capacity and Colonization Dynamics of Curvibacter in the Hydra Host Habitat

Most eukaryotic species are colonized by a microbial community – the microbiota – that is acquired during early life stages and is critical to host development and health. Much research has focused on the microbiota biodiversity during the host life, however, empirical data on the basic ecological principles that govern microbiota assembly is lacking. Here we quantify the contribution of colonizer order, arrival time and colonization history to microbiota assembly on a host. We established the freshwater polyp Hydra vulgaris and its dominant colonizer Curvibacter as a model system that enables the visualization and quantification of colonizer population size at the single cell resolution, in vivo, in real time. We estimate the carrying capacity of a single Hydra polyp as 2 × 10⁵Curvibacter cells, which is robust among individuals and time. Colonization experiments reveal a clear priority effect of first colonizers that depends on arrival time and colonization history. First arriving colonizers achieve a numerical advantage over secondary colonizers within a short time lag of 24 h. Furthermore, colonizers primed for the Hydra habitat achieve a numerical advantage in the absence of a time lag. These results follow the theoretical expectations for any bacterial habitat with a finite carrying capacity. Thus, Hydra colonization and succession processes are largely determined by the habitat occupancy over time and Curvibacter colonization history. Our experiments provide empirical data on the basic steps of host-associated microbiota establishment – the colonization stage. The presented approach supplies a framework for studying habitat characteristics and colonization dynamics within the host–microbe setting.

Evolution of Chaperonin Gene Duplication in Stigonematalean Cyanobacteria (Subsection V)

Chaperonins promote protein folding and are known to play a role in the maintenance of cellular stability under stress conditions. The group I bacterial chaperonin complex comprises GroEL, that forms a barrel-like oligomer, and GroES that forms the lid. In most eubacteria the GroES/GroEL chaperonin is encoded by a single-copy bicistronic operon, whereas in cyanobacteria up to three groES/groEL paralogs have been documented. Here we study the evolution and functional diversification of chaperonin paralogs in the heterocystous, multi-seriate filament forming cyanobacterium Chlorogloeopsis fritschii PCC 6912. The genome of C. fritschii encodes two groES/groEL operons (groESL1, groESL1.2) and a monocistronic groEL gene (groEL2). A phylogenetic reconstruction reveals that the groEL2 duplication is as ancient as cyanobacteria, whereas the groESL1.2 duplication occurred at the ancestor of heterocystous cyanobacteria. A comparison of the groEL paralogs transcription levels under different growth conditions shows that they have adapted distinct transcriptional regulation. Our results reveal that groEL1 and groEL1.2 are upregulated during diazotrophic conditions and the localization of their promoter activity points towards a role in heterocyst differentiation. Furthermore, protein–protein interaction assays suggest that paralogs encoded in the two operons assemble into hybrid complexes. The monocistronic encoded GroEL2 is not forming oligomers nor does it interact with the co-chaperonins. Interaction between GroES1.2 and GroEL1.2 could not be documented, suggesting that the groESL1.2 operon does not encode a functional chaperonin complex. Functional complementation experiments in Escherichia coli show that only GroES1/GroEL1 and GroES1/GroEL1.2 can substitute the native operon. In summary, the evolutionary consequences of chaperonin duplication in cyanobacteria include the retention of groESL1 as a housekeeping gene, subfunctionalization of groESL1.2 and neofunctionalization of the monocistronic groEL2 paralog.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1053-5) contains supplementary material, which is available to authorized users.

Expansion of the redox-sensitive proteome coincides with the plastid endosymbiosis

The redox-sensitive proteome (RSP) consists of protein thiols that undergo redox reactions, playing an important role in coordinating cellular processes. Here, we applied a large-scale phylogenomic reconstruction approach in the model diatom Phaeodactylum tricornutum to map the evolutionary origins of the eukaryotic RSP. The majority of P. tricornutum redox-sensitive cysteines (76%) is specific to eukaryotes, yet these are encoded in genes that are mostly of a prokaryotic origin (57%). Furthermore, we find a threefold enrichment in redox-sensitive cysteines in genes that were gained by endosymbiotic gene transfer during the primary plastid acquisition. The secondary endosymbiosis event coincides with frequent introduction of reactive cysteines into existing proteins. While the plastid acquisition imposed an increase in the production of reactive oxygen species, our results suggest that it was accompanied by significant expansion of the RSP, providing redox regulatory networks the ability to cope with fluctuating environmental conditions.

DnaK-Dependent Accelerated Evolutionary Rate in Prokaryotes

Many proteins depend on an interaction with molecular chaperones in order to fold into a functional tertiary structure. Previous studies showed that protein interaction with the GroEL/GroES chaperonine and Hsp90 chaperone can buffer the impact of slightly deleterious mutations in the protein sequence. This capacity of GroEL/GroES to prevent protein misfolding has been shown to accelerate the evolution of its client proteins. Whether other bacterial chaperones have a similar effect on their client proteins is currently unknown. Here, we study the impact of DnaK (Hsp70) chaperone on the evolution of its client proteins. Evolutionary parameters were derived from comparison of the Escherichia coli proteome to 1,808,565 orthologous proteins in 1,149 proteobacterial genomes. Our analysis reveals a significant positive correlation between protein binding frequency with DnaK and evolutionary rate. Proteins with high binding affinity to DnaK evolve on average 4.3-fold faster than proteins in the lowest binding affinity class at the genus resolution. Differences in evolutionary rates of DnaK interactor classes are still significant after adjusting for possible effects caused by protein expression level. Furthermore, we observe an additive effect of DnaK and GroEL chaperones on the evolutionary rates of their common interactors. Finally, we found pronounced similarities in the physicochemical profiles that characterize proteins belonging to DnaK and GroEL interactomes. Our results thus implicate DnaK-mediated folding as a major component in shaping protein evolutionary dynamics in bacteria and supply further evidence for the long-term manifestation of chaperone-mediated folding on genome evolution.

The Contribution of Genetic Recombination to CRISPR Array Evolution

CRISPR (clustered regularly interspaced short palindromic repeats) is a microbial immune system against foreign DNA. Recognition sequences (spacers) encoded within the CRISPR array mediate the immune reaction in a sequence-specific manner. The known mechanisms for the evolution of CRISPR arrays include spacer acquisition from foreign DNA elements at the time of invasion and array erosion through spacer deletion. Here, we consider the contribution of genetic recombination between homologous CRISPR arrays to the evolution of spacer repertoire. Acquisition of spacers from exogenic arrays via recombination may confer the recipient with immunity against unencountered antagonists. For this purpose, we develop a novel method for the detection of recombination in CRISPR arrays by modeling the spacer order in arrays from multiple strains from the same species. Because the evolutionary signal of spacer recombination may be similar to that of pervasive spacer deletions or independent spacer acquisition, our method entails a robustness analysis of the recombination inference by a statistical comparison to resampled and perturbed data sets. We analyze CRISPR data sets from four bacterial species: two Gammaproteobacteria species harboring CRISPR type I and two Streptococcus species harboring CRISPR type II loci. We find that CRISPR array evolution in Escherichia coli and Streptococcus agalactiae can be explained solely by vertical inheritance and differential spacer deletion. In Pseudomonas aeruginosa, we find an excess of single spacers potentially incorporated into the CRISPR locus during independent acquisition events. In Streptococcus thermophilus, evidence for spacer acquisition by recombination is present in 5 out of 70 strains. Genetic recombination has been proposed to accelerate adaptation by combining beneficial mutations that arose in independent lineages. However, for most species under study, we find that CRISPR evolution is shaped mainly by spacer acquisition and loss rather than recombination. Since the evolution of spacer content is characterized by a rapid turnover, it is likely that recombination is not beneficial for improving phage resistance in the strains under study, or that it cannot be detected in the resolution of intraspecies comparisons.

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.

Genomes of Stigonematalean cyanobacteria (subsection V) and the evolution of oxygenic photosynthesis from prokaryotes to plastids

Cyanobacteria forged two major evolutionary transitions with the invention of oxygenic photosynthesis and the bestowal of photosynthetic lifestyle upon eukaryotes through endosymbiosis. Information germane to understanding those transitions is imprinted in cyanobacterial genomes, but deciphering it is complicated by lateral gene transfer (LGT). Here, we report genome sequences for the morphologically most complex true-branching cyanobacteria, and for Scytonema hofmanni PCC 7110, which with 12,356 proteins is the most gene-rich prokaryote currently known. We investigated components of cyanobacterial evolution that have been vertically inherited, horizontally transferred, and donated to eukaryotes at plastid origin. The vertical component indicates a freshwater origin for water-splitting photosynthesis. Networks of the horizontal component reveal that 60% of cyanobacterial gene families have been affected by LGT. Plant nuclear genes acquired from cyanobacteria define a lower bound frequency of 611 multigene families that, in turn, specify diazotrophic cyanobacterial lineages as having a gene collection most similar to that possessed by the plastid ancestor.