Meiofauna metabarcoding in Lima estuary (Portugal) suggests high taxon replacement within a background of network stability

Meiofauna organisms play an important role in ecological and sedimentary processes in estuarine ecosystems. Recently, the application of environmental DNA (eDNA) for investigating meiofauna in different environments, improved the accessibility to its diversity and composition in a scale, frequency and depth previously unattainable. Nevertheless, little attention has been given to the description of baseline patterns of coupled spatial and temporal dynamics of meiobenthic communities. In an earlier study conducted in Lima estuary (NW Portugal), using eDNA metabarcoding of sediment samples, high levels of meiofauna Operational Taxonomic Units (OTUs) turnover were recorded, between sampling points only a few metres apart, and among sampling stations along the estuary. In order to verify the consistency of these patterns, in the current study we re-assessed Lima estuary's meiofauna communities approximately 1 year after, applying the same methodological approach (targeting segments of the COI and 18S rRNA genes), and expanding HTS-data analyses through the use of association networks. A high degree of spatial turnover was found both within and between sampling stations and this was consistent for both markers and years. As a consequence, most of the beta-diversity was accounted by OTU replacement with only a minor contribution from OTU richness. Despite the high levels of OTU replacement, relatively stable network properties were found in meiofaunal communities, irrespective of the sampled year. Network properties appear to shift sharply from the downstream/high salinity area of the estuary to the mesohaline medium-upstream areas, suggesting high resilience and redundancy of meiofaunal communities along the estuarine gradient. The recognition of meiofauna's networks features may improve the understanding of the ecology and dynamics of these communities that apparently hold large portions of variable elements, thereby making difficult their analyses solely based on the OTU/species composition.This study was funded by the project "The NextSea: Next generation monitoring of coastal ecosystems in a scenario of global change'' (operacao NORTE-01-0145-FEDER-000032), supported by Norte Portugal Regional Operational Programme (NORTE 2020), under the PORTUGAL 2020 Partnership Agreement, through the European Regional Development Fund (ERDF). Maria Fais and Sofia Duarte were supported, respectively by a Ph.D. (SFRH/BD/113547/2015) and a post-doc fellowship (SFRH/BPD/109842/2015) from FCT. The authors would like to thank Prof. Jesus Troncoso (University of Vigo) and Prof. Pedro Gomes (University of Minho) for their availability and hospitality for the lab processing; Dr. Bruno Sampieri and Barbara Leite (University of Minho) for the great support during the implementation of the practical stage of the research. The authors would like also to thank the anonymous reviewers for critically reading the manuscript and suggesting substantial improvements

Identification of Prophages in Bacterial Genomes by Dinucleotide Relative Abundance Difference

BACKGROUND: Prophages are integrated viral forms in bacterial genomes that have been found to contribute to interstrain genetic variability. Many virulence-associated genes are reported to be prophage encoded. Present computational methods to detect prophages are either by identifying possible essential proteins such as integrases or by an extension of this technique, which involves identifying a region containing proteins similar to those occurring in prophages. These methods suffer due to the problem of low sequence similarity at the protein level, which suggests that a nucleotide based approach could be useful. METHODOLOGY: Earlier dinucleotide relative abundance (DRA) have been used to identify regions, which deviate from the neighborhood areas, in genomes. We have used the difference in the dinucleotide relative abundance (DRAD) between the bacterial and prophage DNA to aid location of DNA stretches that could be of prophage origin in bacterial genomes. Prophage sequences which deviate from bacterial regions in their dinucleotide frequencies are detected by scanning bacterial genome sequences. The method was validated using a subset of genomes with prophage data from literature reports. A web interface for prophage scan based on this method is available at http://bicmku.in:8082/prophagedb/dra.html. Two hundred bacterial genomes which do not have annotated prophages have been scanned for prophage regions using this method. CONCLUSIONS: The relative dinucleotide distribution difference helps detect prophage regions in genome sequences. The usefulness of this method is seen in the identification of 461 highly probable loci pertaining to prophages which have not been annotated so earlier. This work emphasizes the need to extend the efforts to detect and annotate prophage elements in genome sequences

Acquisition of the Sda1-encoding bacteriophage does not enhance virulence of the serotype M1 Streptococcus pyogenes strain SF370

The resurgence of invasive disease caused by Streptococcus pyogenes (group A Streptococcus [GAS]) in the past 30 years has paralleled the emergence and global dissemination of the highly virulent M1T1 clone. The GAS M1T1 clone has diverged from the ancestral M1 serotype by horizontal acquisition of two unique bacteriophages, encoding the potent DNase Sda1/SdaD2 and the superantigen SpeA, respectively. The phage-encoded DNase promotes escape from neutrophil extracellular traps and is linked to enhanced virulence of the M1T1 clone. In this study, we successfully used in vitro lysogenic conversion to transfer the Sda1-encoding phage from the M1T1 clonal strain 5448 to the nonclonal M1 isolate SF370 and determined the impact of this horizontal gene transfer event on virulence. Although Sda1 was expressed in SF370 lysogens, no capacity of the phage-converted strain to survive human neutrophil killing, switch to a hyperinvasive covRS mutant form, or cause invasive lethal infection in a humanized plasminogen mouse model was observed. This work suggests that the hypervirulence of the M1T1 clone is due to the unique synergic effect of the M1T1 clone bacteriophage-specific virulence factor Sda1 acting in concert with the M1T1 clone-specific genetic scaffold

Horizontal DNA transfer mechanisms of bacteria as weapons of intragenomic conflict

Horizontal DNA transfer (HDT) is a pervasive mechanism of diversification in many microbial species, but its primary evolutionary role remains controversial. Much recent research has emphasised the adaptive benefit of acquiring novel DNA, but here we argue instead that intragenomic conflict provides a coherent framework for understanding the evolutionary origins of HDT. To test this hypothesis, we developed a mathematical model of a clonally descended bacterial population undergoing HDT through transmission of mobile genetic elements (MGEs) and genetic transformation. Including the known bias of transformation toward the acquisition of shorter alleles into the model suggested it could be an effective means of counteracting the spread of MGEs. Both constitutive and transient competence for transformation were found to provide an effective defence against parasitic MGEs; transient competence could also be effective at permitting the selective spread of MGEs conferring a benefit on their host bacterium. The coordination of transient competence with cell-cell killing, observed in multiple species, was found to result in synergistic blocking of MGE transmission through releasing genomic DNA for homologous recombination while simultaneously reducing horizontal MGE spread by lowering the local cell density. To evaluate the feasibility of the functions suggested by the modelling analysis, we analysed genomic data from longitudinal sampling of individuals carrying Streptococcus pneumoniae. This revealed the frequent within-host coexistence of clonally descended cells that differed in their MGE infection status, a necessary condition for the proposed mechanism to operate. Additionally, we found multiple examples of MGEs inhibiting transformation through integrative disruption of genes encoding the competence machinery across many species, providing evidence of an ongoing "arms race." Reduced rates of transformation have also been observed in cells infected by MGEs that reduce the concentration of extracellular DNA through secretion of DNases. Simulations predicted that either mechanism of limiting transformation would benefit individual MGEs, but also that this tactic's effectiveness was limited by competition with other MGEs coinfecting the same cell. A further observed behaviour we hypothesised to reduce elimination by transformation was MGE activation when cells become competent. Our model predicted that this response was effective at counteracting transformation independently of competing MGEs. Therefore, this framework is able to explain both common properties of MGEs, and the seemingly paradoxical bacterial behaviours of transformation and cell-cell killing within clonally related populations, as the consequences of intragenomic conflict between self-replicating chromosomes and parasitic MGEs. The antagonistic nature of the different mechanisms of HDT over short timescales means their contribution to bacterial evolution is likely to be substantially greater than previously appreciated

The Defective Prophage Pool of Escherichia coli O157: Prophage–Prophage Interactions Potentiate Horizontal Transfer of Virulence Determinants

Bacteriophages are major genetic factors promoting horizontal gene transfer (HGT) between bacteria. Their roles in dynamic bacterial genome evolution have been increasingly highlighted by the fact that many sequenced bacterial genomes contain multiple prophages carrying a wide range of genes. Enterohemorrhagic Escherichia coli O157 is the most striking case. A sequenced strain (O157 Sakai) possesses 18 prophages (Sp1–Sp18) that encode numerous genes related to O157 virulence, including those for two potent cytotoxins, Shiga toxins (Stx) 1 and 2. However, most of these prophages appeared to contain multiple genetic defects. To understand whether these defective prophages have the potential to act as mobile genetic elements to spread virulence determinants, we looked closely at the Sp1–Sp18 sequences, defined the genetic defects of each Sp, and then systematically analyzed all Sps for their biological activities. We show that many of the defective prophages, including the Stx1 phage, are inducible and released from O157 cells as particulate DNA. In fact, some prophages can even be transferred to other E. coli strains. We also show that new Stx1 phages are generated by recombination between the Stx1 and Stx2 phage genomes. The results indicate that these defective prophages are not simply genetic remnants generated in the course of O157 evolution, but rather genetic elements with a high potential for disseminating virulence-related genes and other genetic traits to other bacteria. We speculate that recombination and various other types of inter-prophage interactions in the O157 prophage pool potentiate such activities. Our data provide new insights into the potential activities of the defective prophages embedded in bacterial genomes and lead to the formulation of a novel concept of inter-prophage interactions in defective prophage communities

Host Specific Diversity in Lactobacillus johnsonii as Evidenced by a Major Chromosomal Inversion and Phage Resistance Mechanisms

Genetic diversity and genomic rearrangements are a driving force in bacterial evolution and niche adaptation. We sequenced and annotated the genome of Lactobacillus johnsonii DPC6026, a strain isolated from the porcine intestinal tract. Although the genome of DPC6026 is similar in size (1.97mbp) and GC content (34.8%) to the sequenced human isolate L. johnsonii NCC 533, a large symmetrical inversion of approximately 750 kb differentiated the two strains. Comparative analysis among 12 other strains of L. johnsonii including 8 porcine, 3 human and 1 poultry isolate indicated that the genome architecture found in DPC6026 is more common within the species than that of NCC 533. Furthermore a number of unique features were annotated in DPC6026, some of which are likely to have been acquired by horizontal gene transfer (HGT) and contribute to protection against phage infection. A putative type III restriction-modification system was identified, as were novel Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) elements. Interestingly, these particular elements are not widely distributed among L. johnsonii strains. Taken together these data suggest intra-species genomic rearrangements and significant genetic diversity within the L. johnsonii species and indicate towards a host-specific divergence of L. johnsonii strains with respect to genome inversion and phage exposure

Alternative Sigma Factor σH Modulates Prophage Integration and Excision in Staphylococcus aureus

The prophage is one of the most important components of variable regions in bacterial genomes. Some prophages carry additional genes that may enhance the toxicity and survival ability of their host bacteria. This phenomenon is predominant in Staphylococcus aureus, a very common human pathogen. Bioinformatics analysis of several staphylococcal prophages revealed a highly conserved 40-bp untranslated region upstream of the int gene. A small transcript encoding phage integrase was identified to be initiated from the region, demonstrating that the untranslated region contained a promoter for int. No typical recognition sequence for either σA or σB was identified in the 40-bp region. Experiments both in vitro and in vivo demonstrated that σH recognized the promoter and directed transcription. Genetic deletion of sigH altered the int expression, and subsequently, the excision proportion of prophage DNAs. Phage assays further showed that sigH affected the ability of spontaneous lysis and lysogenization in S. aureus, suggesting that sigH plays a role in stabilizing the lysogenic state. These findings revealed a novel mechanism of prophage integration specifically regulated by a host-source alternative sigma factor. This mechanism suggests a co-evolution strategy of staphylococcal prophages and their host bacteria

Cryptic prophages help bacteria cope with adverse environments

Phages are the most abundant entity in the biosphere and outnumber bacteria by a factor of 10. Phage DNA may also constitute 20% of bacterial genomes; however, its role is ill defined. Here, we explore the impact of cryptic prophages on cell physiology by precisely deleting all nine prophage elements (166 kbp) using Escherichia coli. We find that cryptic prophages contribute significantly to resistance to sub-lethal concentrations of quinolone and β-lactam antibiotics primarily through proteins that inhibit cell division (for example, KilR of rac and DicB of Qin). Moreover, the prophages are beneficial for withstanding osmotic, oxidative and acid stresses, for increasing growth, and for influencing biofilm formation. Prophage CPS-53 proteins YfdK, YfdO and YfdS enhanced resistance to oxidative stress, prophages e14, CPS-53 and CP4-57 increased resistance to acid, and e14 and rac proteins increased early biofilm formation. Therefore, cryptic prophages provide multiple benefits to the host for surviving adverse environmental conditions

The complete genome sequence and genetic analysis of ΦCA82 a novel uncultured microphage from the turkey gastrointestinal system

The genomic DNA sequence of a novel enteric uncultured microphage, ΦCA82 from a turkey gastrointestinal system was determined utilizing metagenomics techniques. The entire circular, single-stranded nucleotide sequence of the genome was 5,514 nucleotides. The ΦCA82 genome is quite different from other microviruses as indicated by comparisons of nucleotide similarity, predicted protein similarity, and functional classifications. Only three genes showed significant similarity to microviral proteins as determined by local alignments using BLAST analysis. ORF1 encoded a predicted phage F capsid protein that was phylogenetically most similar to the Microviridae ΦMH2K member's major coat protein. The ΦCA82 genome also encoded a predicted minor capsid protein (ORF2) and putative replication initiation protein (ORF3) most similar to the microviral bacteriophage SpV4. The distant evolutionary relationship of ΦCA82 suggests that the divergence of this novel turkey microvirus from other microviruses may reflect unique evolutionary pressures encountered within the turkey gastrointestinal system

Assessing the Diversity and Specificity of Two Freshwater Viral Communities through Metagenomics

Transitions between saline and fresh waters have been shown to be infrequent for microorganisms. Based on host-specific interactions, the presence of specific clades among hosts suggests the existence of freshwater-specific viral clades. Yet, little is known about the composition and diversity of the temperate freshwater viral communities, and even if freshwater lakes and marine waters harbor distinct clades for particular viral sub-families, this distinction remains to be demonstrated on a community scale