Skip to main content

Live virus-free or die: coupling of antivirus immunity and programmed suicide or dormancy in prokaryotes

Abstract

Background

The virus-host arms race is a major theater for evolutionary innovation. Archaea and bacteria have evolved diverse, elaborate antivirus defense systems that function on two general principles: i) immune systems that discriminate self DNA from nonself DNA and specifically destroy the foreign, in particular viral, genomes, whereas the host genome is protected, or ii) programmed cell suicide or dormancy induced by infection.

Presentation of the hypothesis

Almost all genomic loci encoding immunity systems such as CRISPR-Cas, restriction-modification and DNA phosphorothioation also encompass suicide genes, in particular those encoding known and predicted toxin nucleases, which do not appear to be directly involved in immunity. In contrast, the immunity systems do not appear to encode antitoxins found in typical toxin-antitoxin systems. This raises the possibility that components of the immunity system themselves act as reversible inhibitors of the associated toxin proteins or domains as has been demonstrated for the Escherichia coli anticodon nuclease PrrC that interacts with the PrrI restriction-modification system. We hypothesize that coupling of diverse immunity and suicide/dormancy systems in prokaryotes evolved under selective pressure to provide robustness to the antivirus response. We further propose that the involvement of suicide/dormancy systems in the coupled antivirus response could take two distinct forms:

1) induction of a dormancy-like state in the infected cell to ‘buy time’ for activation of adaptive immunity; 2) suicide or dormancy as the final recourse to prevent viral spread triggered by the failure of immunity.

Testing the hypothesis

This hypothesis entails many experimentally testable predictions. Specifically, we predict that Cas2 protein present in all cas operons is a mRNA-cleaving nuclease (interferase) that might be activated at an early stage of virus infection to enable incorporation of virus-specific spacers into the CRISPR locus or to trigger cell suicide when the immune function of CRISPR-Cas systems fails. Similarly, toxin-like activity is predicted for components of numerous other defense loci.

Implications of the hypothesis

The hypothesis implies that antivirus response in prokaryotes involves key decision-making steps at which the cell chooses the path to follow by sensing the course of virus infection.

Reviewers

This article was reviewed by Arcady Mushegian, Etienne Joly and Nick Grishin. For complete reviews, go to the Reviewers’ reports section.

Background

Viruses are the most abundant biological entities on earth. In well-characterized habitats such as seawater and soil the number of viral particles exceeds the number of cells by one to two orders of magnitude [13]. Thus, all bacteria and archaea exist in a perennial arms race with the excessively abundant viruses [4, 5]. Consequently, prokaryotes have evolved extremely diverse and elaborate antiviral defense systems that occupy a substantial part of the genome in all free-living prokaryotes [6]. Antivirus defense systems can be classified into two broad categories that differ in their general principles of action. Immune systems function on the self-nonself discrimination principle, i.e. specifically recognize and destroy foreign genomes while protecting the host genome (Figure 1). In addition to targeting viral genomes, these systems are also involved in other inter-genomic conflicts with selfish elements such as plasmids. In contrast, suicide systems execute cell death or dormancy programs that prevent a virus from completing its reproduction in the given infected cell and subsequently infecting other cells [7]. The self-nonself discrimination principle is employed, in particular, by the restriction-modification (RM) systems, which are probably the best characterized defense systems in prokaryotes, to a large extent, because restriction endonucleases are essential experimental tools of molecular biology [810]. Recently, an analogous system of DNA phosphorothioation (known as the DND system) has been characterized [1113]. The RM and DND systems may be considered mechanisms of innate immunity: they do not adapt to a specific infectious agent but simply ensure protection of the self DNA and attack non-self invaders indiscriminately. The CRISPR (clustered regularly interspaced short palindromic repeats)-Cas (CRISPR-associated genes) systems that are encoded in the genomes of the great majority of archaea and many bacteria represent another type of defense machinery that is also based on self-nonself discrimination but fits the definition of adaptive (as opposed to innate) immunity [1418]. Unlike the RM and DND systems that generically distinguish between modified and unmodified recognition sites in DNA, the CRISPR-Cas systems function via adaptation to a specific infectious agent. The CRISPR-Cas first incorporates a unique fragment of the invading DNA into a specific locus in the host genome and then employs the transcripts of these unique spacers to target the cognate sequences in viral or plasmid genomes, achieving extremely high levels of immunity as the result.

Figure 1
figure 1

The defense systems in prokaryotes: innate immunity, adaptive immunity and programmed suicide/dormancy.

The toxin-antitoxin (TA) and abortive infection (ABI) systems that are nearly ubiquitous and highly abundant in bacteria and archaea typically cause cell suicide or dormancy in response to virus infection or other forms of stress [7, 1926]. Under normal conditions, the toxin component of a TA system is kept inactive through the association with the antitoxin but under stress the antitoxin is inactivated and the toxin is unleashed. These systems employ several molecular mechanisms, the most common one being the cleavage of translated mRNA directly on the ribosomes by toxins that possess nuclease (interferase) activity. The evolution of a suicidal response via toxin action might be explained by invoking altruism/kin/group selection: when immune systems fail to prevent virus reproduction, the TA and ABI systems kill the infected cell or render it dormant, hence increasing the survival chance for the surrounding cells [7, 1926]. Such altruistic behavior is particularly likely to be selected for in situations where the neighboring cells are kin emerging from clonal expansion and/or when cells have other factors favoring cooperation, such as aggregation in a biofilm. Indeed, a recent study combining experimental analysis of virus-host interaction with mathematical modeling and simulation has shown that bacterial “Suicidal Defense Against Infection” is efficient in a spatially structured habitat but under well-mixed conditions [27]. The dormancy route not only attains the same end but additionally provides potential for survival of the infected cell itself.

Comparative genomic analysis of the genes encoding defense systems in archaea and bacteria revealed statistically significant genomic clustering of RM, ABI and TA systems (but not CRISPR-Cas) within defense islands [6]. Furthermore, several specific, tight associations were observed. In particular, the abortive infection system AbiU1/AIPR is frequently associated with RM [28]. Two additional architectural associations of the same type have been reported for two previously uncharacterized domains Ymh (Pfam09509) and DUF4145 (Pfam13643), with the latter domain being fused to restriction subunits of RM type I systems.

Whether the connection between immune and suicidal defense systems involves any form(s) of specific functional cooperation or occurs fortuitously due to non-adaptive clustering of horizontally transferred genes remains unclear [6]. However, in one well-studied model such cooperation has been discovered. This case involves the Escherichia coli anticodon nuclease (ACNase) PrrC which contributes to the T4 phage exclusion mechanism as a component of the RM type Ic system PrrI [29]. The prrC gene co-localizes with the cluster of genes for the 3 subunits of a RM system (hsdMSR), and this genomic association is conserved in diverse bacteria, suggestive of functional coupling. Normally, PrrC persists as an inactive, latent endoribonuclease in association with the Hsd complex but it can be allosterically activated either by unmodified DNA, or by increased levels of dTTP or by the small anti-restriction peptide encoded by the T4-like enterobacteriophages. The activated PrrC ACNase cleaves the anticodon of tRNALys in a GTP-dependent manner, with GTP hydrolysis catalyzed by the N-terminal ABC NTPase domain of PrrC. The cleavage of tRNALys inhibits translation and abrogates the growth of pnk (polynucleotide kinase) or rnl1 (RNA ligase) mutants of T4. In wild type T4 and related phages, Pnl and Rnl1 jointly function as an antidote to PrrC by repairing the cleaved tRNA at the 2′,3′ cyclic phosphate and 5′ OH [30, 31]. The ACNase RloC that is homologous to PrrC, although usually not linked directly to RM systems, seems to function under the same logic, i.e. as an antiviral contingency that acts when DNA restriction is alleviated under genotoxic stress [32, 33]. Although PrrC and RloC are not known to cause bacterial suicide or persistence, they clearly function on the exact same principle, namely inhibiting translation to prevent virus reproduction in situations when an innate anti-phage immunity mechanism fails. Indeed, an analogous but completely independent mechanism of translation inhibition has been demonstrated to be critical in a parallel defense system that is deployed against T4 by certain strains of Escherichia coli. Here, a defective prophage produces a toxin known as Lit that contains a zincin-like metallopeptidase domain and cleaves the elongation factor Tu thereby instigating cell suicide or at least abolishing translation in the bacterial cell infected by T4. This system has been shown to prevent the spread of T4 to the remaining cells in the colony [34, 35].

We were interested in the possibility that echeloning of suicidal/persistence-based defenses and directed attack on viral nucleic acids via genomic coupling could be a general phenomenon in prokaryotes. Here we present evidence of association of TA/ABI systems with RM, DND and CRISPR-Cas systems and outline the hypothesis that all classes of prokaryotic immunity systems function in direct cooperation with suicide or persistence (dormancy) mechanisms.

Presentation of the hypothesis

We first consider the CRISPR-Cas systems of adaptive immunity. The recent experimental breakthroughs seem to have revealed all the genes that are responsible for DNA/RNA targeting [15, 3643] and spacer integration [44] by CRISPR-Cas. However, several genes that are stably associated with CRISPR-Cas loci still do not fit into this scheme [45] (Figure 2A).

Figure 2
figure 2

Organization of the genomic loci that encode prokaryotic immune systems including toxin genes. The core genes of CRISPR-Cas, RM, and DND systems in predicted operons are shown by pink arrows; genes with (predicted) toxin activity are shown by different colors, and the (predicted) toxin domains are denoted with red outline. The Csa3 protein in the Type IA system lacks the HEPN domain. Abbreviations: HEPN - higher eukarytoes and prokaryotes nucleotide-binding domain; Sir2, ParB and REase, DEDD are nucleases from distinct superfamilies. A. CRISPR-Cas. Gene names follow the nomenclature and classification from [18]. B. Restriction-modification. Gene names follow the nomenclature and classification from [46]. C. Phosphorothioation. Gene names follow the nomenclature from [12]. D. Prokaryotic argonaut, pAgo.

The only two genes that are universal among the CRISPR-cas systems are cas1 and cas2 both of which are implicated in the first stage of CRISPR-Cas function, adaptation, or spacer integration [44]. However, despite the genetic evidence of the involvement of both genes, all the enzymatic activities required for adaptation appear to be provided by Cas1 alone [44]. The Cas2 protein has been shown to possess a sequence-specific endoribonuclease activity [47] and is a homolog of the VapDHi toxin of the VapDHi/VapX TA system from Haemophilus influenzae[4750]. Therefore it appears highly likely that like VapDHi, Cas2 is an mRNA interferase that specifically cleaves ribosome-associated mRNAs, an activity that is mechanistically irrelevant for the adaptation stage of the CRISPR-Cas function. Furthermore, certain Cas2 homologs occur independently of an antitoxin in predicted operons or in genomic contexts corresponding to mobile selfish elements that are able to mediate their replication in conjunction with transfer. In these cases, the Cas2-homologs are typically in the neighborhood of a gene encoding a resolvase (e.g. gi: 57504998, from Campylobacter coli RM2228) or a Mob-type relaxase (gi: 313144877 from Helicobacter cinaedi). Conceivably, the endoRNase activity of these Cas2 homologs might foster addiction to the respective mobile elements via toxin-like action upon disruption of the selfish element. Furthermore, the combination of the Cas2-like protein with a DNA resolving/cleaving enzyme is reminiscent of the CRISPR-Cas systems where Cas1, the partner of Cas2, plays the key role in the integration of the acquired spacer DNA. This parallel seems to support the potential origin of the Cas1-Cas2 dyad of the CRISPR-Cas systems from an ancient mobile element [47] similar to the aforementioned elements that combine Cas2-like genes with genes encoding DNA resolving/cleaving enzymes. Therefore, we hypothesize that Cas2 retains its ancestral toxin-like EndoRNase function within the CRISPR-Cas systems but the interferase activity is kept in check through reversible inhibition by the Cas2-Cas1 interaction. This scheme is proposed as a direct analogy to the control mechanism of PrrC which is kept in its latent state via the interaction with the associated RM system.

Under this hypothesis, when CRISPR-Cas fails to contain virus growth due to viral counter-attack and/or the level of genotoxic stress increases due to the accumulation of nucleotide metabolites synthesized by viral enzymes, Cas2 is activated (possibly through degradation of Cas1) and abrogates translation, probably leading to cell suicide or dormancy (Figure 3). The requirement of Cas2 for spacer integration might stem from regulation or stabilization of Cas1 via Cas1-Cas2 complex formation that also reversibly inactivates Cas2. Perhaps even more important, Cas2 might act at the initial stage of the CRISPR-Cas response by rendering the infected cell dormant and hence ‘buying time’ to allow the CRISPR-Cas system to integrate virus-specific spacers for effective future use in antivirus response. The strong sequence conservation of Cas2 across the CRISPR-Cas systems, together with its sequence-specific endoribonuclease activity, appear to be better consistent with “self” RNAs being targeted rather than non-self RNAs which would likely select for greater diversity of Cas2 [51]. Of interest in this regard is the fusion of Cas2 with a 3′-5′ exonuclease RNaseH fold domain (e.g. LSEI_0356 from Lactobacillus casei ATCC 334), which is present in several Type I-E CRISPR-Cas loci. In these fusion proteins, Cas2 appears to be inactivated as a result of substitution of the catalytic residues so that the ribonuclease activity is most likely supplied by the 3′-5′exonuclease-like domain [52].

Figure 3
figure 3

The immunity-dormancy/suicide coupling hypothesis: Route 1 – toxin action before immunity activation. The coupling between dormancy-suicide and immunity is specifically illustrated by the CRISPR-Cas system that is hypothesized to adapt by inserting a cognate spacer during a dormancy-like phase induced by the action of the toxin (typically, Cas2).

Most of the Type I and Type III CRISPR-Cas systems encompass additional domains, besides Cas2, which could potentially function similarly to toxins (Figure 2A). Typically, these predicted toxins are either fused to or encoded in the same operon with proteins of COG1517 (also known as Csm6 and Csx1) that are common in CRISPR-Cas systems (Figure 2A). Most COG1517 proteins contain a distinct Rossmann-fold domain [5, 49] but many in addition contain a third, usually C-terminal domain. The domains fused to COG1517 proteins include REase fold DNases (Pfam Clan: PD-(DE)xK) [49], as confirmed by a recently solved crystal structure (pdb code 1XMX), and HEPN domains ([53] and (KSM, VA, EVK, LA, unpublished). The HEPN domain has been predicted to comprise a component of a distinct TA jointly with the minimal nucleotidyltransferase domain [25]; recently, HEPN has been shown to function as a toxin [54]. Our recent analysis suggests that the PrrC/RloC ACNase domain and the KEN domain of eukaryotic unfolded protein response/ antiviral RNAses (Ire1/RNAse L) also are distinct versions of the HEPN fold and accordingly that several HEPN domains are RNases with a toxin-type activity (KSM, VA, EVK, LA, unpublished). Additional notable architectures of this CRISPR-CAS associated Rossmann-like domain include fusions to distinct RNase domains with potential toxin activity such as the RelE (e.g. sll7062 from Synechocystis sp. PCC 6803) and PIN (e.g. APE_2119.1, Aeropyrum pernix K1) [55]. Both these toxins have been shown to possess mRNA interferase activity [24]. Most of the predicted toxin-like genes are located within or in a close proximity of predicted Type III CRISPR-Cas operons [18].

Moreover, the widespread Cas4 protein, a REase fold nuclease that is part of most Type I CRISPR-Cas systems but for which no specific function in the CRISPR-Cas response so far has been established also might exhibit a toxin-like activity. In this respect three observations seem to be relevant: first, the aforementioned fusion of the REase fold nuclease domain with the Rossmannoid domain; second, the only CRISPR-Cas Type I-A locus so far detected in which the COG1517 protein (Csa3) does not contain a potential toxin domain encompasses Cas4 (KSM, unpublished); third, the fusion of Cas4 with Cas1 [18]. The Cas4-Cas1 fusion may be considered in parallel with the fusion of Cas1 with the reverse transcriptase (RT) domain that was described previously [49] and now appears to be widespread in a variety of bacteria (e.g. alr1468 from Nostoc sp. PCC 7120, Franean1_1369 from Frankia sp. EAN1pec, VVA1544 from Vibrio vulnificus YJ016, etc.). The RT domain fused with Cas1 is related to the RT-like proteins of the AbiA and AbiK families that participate in the abortive infection response [56]. Recently, it has been shown that AbiK is not a typical RT but apparently catalyzes non-templated synthesis of random sequence DNA that remains covalently attached to the protein and contributes to abortive infection via a still uncharacterized mechanism [57]. We propose that the Cas1-associated RT functions in a similar fashion. Thus, within the CRISPR-Cas systems, multiple cases of replacement of one type of toxin by another seem to occur.

Similarly to CRISPR-Cas, defense systems that are best classified as innate immunity are also associated with suicide/dormancy genes. In addition to the PrrC-PrrI link discussed above, many other Type I and Type III RM systems encompass domains with potential toxin-like roles including different families of HEPN domains as well as Sir2, ParB and REase fold nucleases that are homologous to abortive infection system subunits and predicted toxins of TA systems [6] (Figure 2B). Thus, in the RM system the same phenomenon of repeated toxin displacement as in the CRISPR-Cas systems seems to take place.

Although the specific functions of the proteins involved in DNA phosphorothioation and those that are involved in the restriction by this system are poorly understood, the sets of genes involved in both processes are relatively well-defined [1113] (Figure 2C). One of the genes (dndB) that is part of the phosphorothioation operon is a negative regulator of the process and might play a role in discriminating between the modification of self and non-self DNA [13, 58]. The N-terminal domain of DndB belongs to the AbiU1/AIPR family of abortive infection proteins that are also linked to Type III RM systems [28]. The AIPR proteins contain an N-terminal HEPN domain that is a toxin and a predicted RNase (see above). Another HEPN domain is present in the DndF protein that is implicated in the restriction process [1113].

Other, poorly characterized prokaryotic defense systems also might contain toxin domains (Figure 2D). In particular, it has been predicted that prokaryotic Argonaut protein (pAgo) is a key component of a distinct defense system that can be encoded by a stand-alone gene or within putative operons together with several nucleases [59]. Notably, one of these nucleases is a close homolog of Cas4 that might function as a toxin in conjunction with CRISPR-Cas systems [59]. In addition, other REase fold DNases as well as a predicted Sir2 superfamily enzyme are often found in association with pAgo [59]. Thus, notably, (predicted) Sir2 superfamily enzymes are distinct components of the ABI system AbiH, and are linked both to Type I RM systems and to the putative defense system centered around pAgo (Figure 2A and 2D). These enzymes might either function as mono-ADP-ribosylating toxins targeting proteins or DNA [60] or as nucleases as has been proposed for the Sir2 enzymes found in systems containing HerA-FtsK-like proteins [61].

Most likely, there are additional, still uncharacterized innate and/or adaptive immunity defense systems in prokaryotes that are also coupled with toxins. One recently studied example is the phage and chemical stress resistance systems centered on the ter genes [62]. Several operons of these systems combine genes encoding nucleases and helicases similar to those found in RM systems with genes encoding COG1517 family proteins as well as putative toxins related to those in typical TA systems. The latter proteins include Doc domains, which modify proteins by addition of NMP moieties, and predicted RNases homologous to RelE, barnase or the PIN domain. It seems likely that these systems also engender higher level functional cooperation between the restriction DNAses that directly target invading DNA and toxins that promote dormancy or suicide. These operons also share components with the Pgl system of reverse restriction-modification that encompasses many genes whose functions remain unclear although some of these encode domains suggestive of toxin-like activities [6] (VA, L.M. Iyer, LA, unpublished).

We are now in a position to formulate the immunity-dormancy/suicide coupling hypothesis according to which immune systems and cell dormancy/suicide systems such as TA are functionally coupled and directly regulate each other’s functions (Figure 3). Stand-alone suicide/dormancy modules, such as typical TA systems, lack a direct, demonstrable role in immunity and might even decrease the fitness of the host, especially when they exhibit their toxic activity in the context of plasmid addiction or restriction attack upon disruption of RM systems. However, the comparative genomics data summarized here indicate that all major antivirus immunity systems identified to date in prokaryotes are linked to proteins that are homologous to toxin subunits of typical suicide/dormancy modules and are predicted to possess toxin-like activities. These toxin-like proteins might not be mechanistically involved in the immune function but rather could act similarly to their toxin counterparts from stand-alone suicide/dormancy modules.

In most cases, the characterized and predicted toxin proteins in immune systems are highly variable both in terms of (potential) targets (i.e. DNA, tRNA, rRNA or mRNA) and mechanism of action (i.e. protein modification or nucleic acid degradation). However, if our prediction of the toxin-like interferase function for Cas2 is valid, the coupling between immunity and suicidal/dormancy response is an intrinsic feature of all CRISPR-Cas systems. Indeed, such coupling might be critical for the effectiveness of the adaptive CRISPR-Cas systems. Given that antivirus defense by these adaptive immunity systems depends on preliminary infection to integrate virus-specific spacer DNA, a dormancy response through the action of Cas2 is likely to help the host to ‘buy time’ to prime the immunity response. Under the current hypothesis, this is the first route of involvement of dormancy/suicide modules in coupled defense whereby the dormancy-like response precedes and enables the immune response (Figure 3). The second route of immunity-suicide coupling follows an even more straightforward biological rationale: when an immunity system fails and/or the level of genotoxic stress increases through the activity of metabolic enzymes encoded by the infecting virus, the cell employs the associated toxins for drastic measures that involve abrogation of key cell processes, typically translation, leading to cell death or dormancy (Figure 4).

Figure 4
figure 4

The immunity-dormancy/suicide coupling hypothesis: Route 2 – toxin action after immunity failure. The coupling between dormancy-suicide and immunity is illustrated by the CRISPR-Cas system as in Figure 3.

Immunity systems do not encode typical antitoxins. Therefore (and by analogy with the PrrC model), we hypothesize that components of the immunity systems themselves (e.g. Cas1) function as reversible inhibitors of the toxic components, and that they are inactivated by genotoxic stress or viral counter-defenses neutralizing the immunity system and thus unleashing the toxin (Figures 3 and 4).

Testing the hypothesis

Clearly, the immunity-suicide coupling hypothesis generates numerous predictions that can be tested experimentally on diverse prokaryotic defense systems. Such tests would involve biochemical experiments to detect the activity of predicted toxin-like proteins, such as Cas2, as well as computational searches for novel toxin-like proteins among products of completely uncharacterized genes within defense loci. Given that many of the toxins are highly diverged, small proteins [25], the latter approach is expected to yield numerous new findings. In addition, in trans interactions between immunity systems and TA modules cannot be ruled out. Beyond direct characterization of enzymatic activities, physical interactions between (predicted) toxins and components of immune systems (e.g. Cas1-Cas2) represent important targets for experimental testing of the hypothesis.

Implications of the hypothesis

A general implication of the immunity-suicide coupling hypothesis is that infected bacterial or archaeal cells routinely make decisions on the nature of their response to virus infection: induction of a dormancy-like state to allow priming of immunity; immediate immune response; or dormancy/suicide in the face of immune system failure. Apparently, these decisions are made through sensing the course of the virus infection. Such decision-making could be highly complex and might potentially involve diverse signal transduction systems similarly to the sporulation-competence decision in Bacillus subtilis that has been mathematically modeled and experimentally characterized [63, 64]. Mathematical modeling of the immunity-dormancy networks has the potential to reveal unexpected properties and functional regimes of prokaryotic defense systems.

Reviewers’ reports

Reviewer 1: Dr. Arcady Mushegian, Stowers Institute for Medical Research, Kansas City

The study by Makarova et al. notes that most of the operons encoding bacterial immunity systems (i.e., restriction-modification, DNA phosphorothioation, and CRISP-Cas) are associated with the known or predicted suicide/dormancy systems (i.e., toxin-antitoxin and abortive infection). The hypothesis is that bacteria are playing zone defense, which apparently can take two different forms: first, if immunity system fails to mount an efficient response to non-self DNA, then a suicide system is activated to kill the cell and prevent the completion of virus infection; second, if immunity system is slow to mount a response to non-self DNA, then a dormancy system may shift cell metabolism in low gear, perhaps slowing down virus reproduction and ‘buying time’ for the immunity system to catch up. Somehow, these two versions of the echeloning are not presented as two distinct responses, despite their opposite outcome, but are discussed interminglingly; perhaps it would be better to state them as the alternatives.

Authors’ response: We agree that the two routes of dormancy-immunity coupling, one where dormancy is proposed to enable immunity and the other one where the transition to dormancy or actual suicide follows immunity failure, are best clearly distinguished. This is how we treat the subject in the revised manuscript including the Abstract and the modified Figure 3and new Figure 4.

The idea, nonetheless, is of a considerable interest. It comes with many specific predictions of novel functions and domain architectures of the involved proteins. In addition to all this, it is hypothesized that Cas system itself has a built-in Kingston valve, consisting of the ubiquitous Cas2 component that acts as a suicidal toxin when released from its inhibitory interaction with Cas1, which has an antitoxin function. The main open question in my mind is whether the tight control of the zone defense of the suicidal variety is possible. Failure or slowness of the cell to

mount the immunity response has to be sensed neither too early (otherwise, futile destruction of the cell that was mounting said response just fine) nor too late (otherwise, destruction only releases already-formed virus progeny). This has to be recognized as an intrinsic difficulty; perhaps dormancy response is more flexible in this regard.

Authors’ response: This comment is very much along the lines of the Implications section in the revised manuscript. As we note in this section, mathematical modeling of the coupled immunity-dormancy response is likely to reveal different regimes and perhaps provide at least tentative answers to the questions posed in this comment that could then become experimentally tractable.

Minor comments:

“[defense island formation] occurs fortuitously due to non-adaptive clustering of horizontally transferred genes” --- perhaps horizontally transferred genes may also cluster adaptively, but not because of their molecular functions but because, for example of the advantage in avoiding the established genome neighborhoods?

Authors’ response: indeed, this is conceivable as discussed previously[6].

“Cas2 .. is a homolog of the VapD toxin of the VapD/VapX TA system [4144]. Therefore it appears highly likely that analogously to VapD, Cas2 is an mRNA interferase…” --- two homologs have the same function; somehow, ‘analogously’ does not belong here.

Authors’ response: perhaps, this was indeed awkward, modified.

is it out of the realm of possibility that Cas-associated reverse transcriptases

could act on virus RNA, rare as the RNA bacteriophages are?

Authors’ response: the data discussed in the article seem to point in a difference direction; yet, we agree, the involvement of the RT in the CRISPR-Cas response to RNA viruses is not entirely inconceivable

Reviewer 2: Dr. Etienne Joly, Université Paul Sabatier, Toulouse

This manuscript presents an interesting hypothesis whereby kin selection/altruism may play a very significant role in the way bacteria combat against viruses, by committing suicide (or dormancy) rather than letting the virus complete its infectious cycle. The main ground for this hypothesis is that many, if not most, antiviral operons contain proteins that appear more likely to play a role in causing the suicide (or dormancy) of the host than against viruses. My main comment is that I wish to congratulate the authors for the remarkable improvement of their manuscript compared to the initial version. The manuscript is now much easier to read and comprehend for someone like me who is not directly involved in this field: it is now much easier to see what was already known, what the speculations and hypotheses presented are, and to see how these could be confirmed or infirmed by the predictions made.

Authors’ response: We greatly appreciate these comments.

Reviewer 3: Dr. Nick Grishin, University of Texas Southwestern Medical Center, Dallas

In this short manuscript an interesting hypothesis is suggested that immunity and suicide in bacteria may be linked by location of immunity genes and suicide nuclease genes within the same operon. Moreover, due to the lack of known antitoxins encoded in such operons, immunity proteins may serve as reversible inhibitors of these nucleases that would cause cell death when immune response fails. The hypothesis is supported by the following evidence: abundance of nucleases without clear functional roles in the operons of all major bacterial immunity systems; extensive functional characterization of CRISPR-cas system failed to pinpoint the function of essential Cas2 suggesting that Cas1 enzymatic activity is sufficient; Cas2 is homologous to the VapD toxin (nuclease); deteriorated Cas2 proteins may be fused with other, possibly active nuclease domains. All these small pieces of incidental evidence and more are not likely to be fortuitous and indeed the hypothesis suggested by the authors seems very attractive and worthy of experimental exploration. However, it is also possible that some of these nucleases do have other functions in immunity and such functions have not been discovered yet. In any case, experimental study of an uncharacterized protein is best to start from a testable functional hypothesis and this work provides a highly plausible one. In addition, introduction and parts of the results section read like an excellent review and a primer in bacterial immunity for readers who are not particularly familiar with the subject. In these regards, I think that the following recent papers might be relevant to the discussion [4, 27]

Authors’ response: We appreciate these comments and cited the suggested, indeed highly relevant references in the revised manuscript.

References

  1. Suttle CA: Marine viruses–major players in the global ecosystem. Nat Rev Microbiol. 2007, 5 (10): 801-812. 10.1038/nrmicro1750.

    Article  PubMed  CAS  Google Scholar 

  2. Suttle CA: Viruses in the sea. Nature. 2005, 437 (7057): 356-361. 10.1038/nature04160.

    Article  PubMed  CAS  Google Scholar 

  3. Edwards RA, Rohwer F: Viral metagenomics. Nat Rev Microbiol. 2005, 3 (6): 504-510. 10.1038/nrmicro1163.

    Article  PubMed  CAS  Google Scholar 

  4. Stern A, Sorek R: The phage-host arms race: shaping the evolution of microbes. Bioessays. 2011, 33 (1): 43-51. 10.1002/bies.201000071.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  5. Buckling A, Brockhurst M: Bacteria-virus coevolution. Adv Exp Med Biol. 2012, 751: 347-370. 10.1007/978-1-4614-3567-9_16.

    Article  PubMed  CAS  Google Scholar 

  6. Makarova KS, Wolf YI, Snir S, Koonin EV: Defense islands in bacterial and archaeal genomes and prediction of novel defense systems. J Bacteriol. 2011, 193 (21): 6039-6056. 10.1128/JB.05535-11.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  7. Labrie SJ, Samson JE, Moineau S: Bacteriophage resistance mechanisms. Nat Rev Microbiol. 2010, 8 (5): 317-327. 10.1038/nrmicro2315.

    Article  PubMed  CAS  Google Scholar 

  8. Kovall RA, Matthews BW: Type II restriction endonucleases: structural, functional and evolutionary relationships. Curr Opin Chem Biol. 1999, 3 (5): 578-583. 10.1016/S1367-5931(99)00012-5.

    Article  PubMed  CAS  Google Scholar 

  9. Williams RJ: Restriction endonucleases: classification, properties, and applications. Mol Biotechnol. 2003, 23 (3): 225-243. 10.1385/MB:23:3:225.

    Article  PubMed  CAS  Google Scholar 

  10. Orlowski J, Bujnicki JM: Structural and evolutionary classification of Type II restriction enzymes based on theoretical and experimental analyses. Nucleic Acids Res. 2008, 36 (11): 3552-3569. 10.1093/nar/gkn175.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  11. Xu T, Yao F, Zhou X, Deng Z, You D: A novel host-specific restriction system associated with DNA backbone S-modification in Salmonella. Nucleic Acids Res. 2010, 38 (20): 7133-7141. 10.1093/nar/gkq610.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  12. He X, Ou HY, Yu Q, Zhou X, Wu J, Liang J, Zhang W, Rajakumar K, Deng Z: Analysis of a genomic island housing genes for DNA S-modification system in Streptomyces lividans 66 and its counterparts in other distantly related bacteria. Mol Microbiol. 2007, 65 (4): 1034-1048. 10.1111/j.1365-2958.2007.05846.x.

    Article  PubMed  CAS  Google Scholar 

  13. Liang J, Wang Z, He X, Li J, Zhou X, Deng Z: DNA modification by sulfur: analysis of the sequence recognition specificity surrounding the modification sites. Nucleic Acids Res. 2007, 35 (9): 2944-2954. 10.1093/nar/gkm176.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  14. Barrangou R, Horvath P: CRISPR: new horizons in phage resistance and strain identification. Annu Rev Food Sci Technol. 2012, 3: 143-162. 10.1146/annurev-food-022811-101134.

    Article  PubMed  CAS  Google Scholar 

  15. Wiedenheft B, Sternberg SH, Doudna JA: RNA-guided genetic silencing systems in bacteria and archaea. Nature. 2012, 482 (7385): 331-338. 10.1038/nature10886.

    Article  PubMed  CAS  Google Scholar 

  16. Karginov FV, Hannon GJ: The CRISPR system: small RNA-guided defense in bacteria and archaea. Mol Cell. 2010, 37 (1): 7-19. 10.1016/j.molcel.2009.12.033.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  17. van der Oost J, Jore MM, Westra ER, Lundgren M, Brouns SJ: CRISPR-based adaptive and heritable immunity in prokaryotes. Trends Biochem Sci. 2009, 34 (8): 401-407. 10.1016/j.tibs.2009.05.002.

    Article  PubMed  CAS  Google Scholar 

  18. Makarova KS, Haft DH, Barrangou R, Brouns SJ, Charpentier E, Horvath P, Moineau S, Mojica FJ, Wolf YI, Yakunin AF, et al: Evolution and classification of the CRISPR-Cas systems. Nat Rev Microbiol. 2011, 9 (6): 467-477. 10.1038/nrmicro2577.

    Article  PubMed  CAS  Google Scholar 

  19. Yamaguchi Y, Inouye M: Regulation of growth and death in Escherichia coli by toxin-antitoxin systems. Nat Rev Microbiol. 2011, 9 (11): 779-790. 10.1038/nrmicro2651.

    Article  PubMed  CAS  Google Scholar 

  20. Chopin MC, Chopin A, Bidnenko E: Phage abortive infection in lactococci: variations on a theme. Curr Opin Microbiol. 2005, 8 (4): 473-479. 10.1016/j.mib.2005.06.006.

    Article  PubMed  CAS  Google Scholar 

  21. Maisonneuve E, Shakespeare LJ, Jorgensen MG, Gerdes K: Bacterial persistence by RNA endonucleases. Proc Natl Acad Sci U S A. 2011, 108 (32): 13206-13211. 10.1073/pnas.1100186108.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  22. Tashiro Y, Kawata K, Taniuchi A, Kakinuma K, May T, Okabe S: RelE-mediated dormancy is enhanced at high cell density in Escherichia coli. J Bacteriol. 2012, 194 (5): 1169-1176. 10.1128/JB.06628-11.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  23. Gerdes K, Christensen SK, Lobner-Olesen A: Prokaryotic toxin-antitoxin stress response loci. Nat Rev Microbiol. 2005, 3 (5): 371-382. 10.1038/nrmicro1147.

    Article  PubMed  CAS  Google Scholar 

  24. Hayes F, Van Melderen L: Toxins-antitoxins: diversity, evolution and function. Crit Rev Biochem Mol Biol. 2011, 46 (5): 386-408. 10.3109/10409238.2011.600437.

    Article  PubMed  CAS  Google Scholar 

  25. Makarova KS, Wolf YI, Koonin EV: Comprehensive comparative-genomic analysis of type 2 toxin-antitoxin systems and related mobile stress response systems in prokaryotes. Biol Direct. 2009, 4: 19-10.1186/1745-6150-4-19.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Blower TR, Salmond GP, Luisi BF: Balancing at survival’s edge: the structure and adaptive benefits of prokaryotic toxin-antitoxin partners. Curr Opin Struct Biol. 2011, 21 (1): 109-118. 10.1016/j.sbi.2010.10.009.

    Article  PubMed  CAS  Google Scholar 

  27. Fukuyo M, Sasaki A, Kobayashi I: Success of a suicidal defense strategy against infection in a structured habitat. Sci Rep. 2012, 2: 238-

    Article  PubMed  PubMed Central  Google Scholar 

  28. Iyer LM, Abhiman S, Aravind L: MutL homologs in restriction-modification systems and the origin of eukaryotic MORC ATPases. Biol Direct. 2008, 3: 8-10.1186/1745-6150-3-8.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Penner M, Morad I, Snyder L, Kaufmann G: Phage T4-coded Stp: double-edged effector of coupled DNA and tRNA-restriction systems. J Mol Biol. 1995, 249 (5): 857-868. 10.1006/jmbi.1995.0343.

    Article  PubMed  CAS  Google Scholar 

  30. Kaufmann G: Anticodon nucleases. Trends Biochem Sci. 2000, 25 (2): 70-74. 10.1016/S0968-0004(99)01525-X.

    Article  PubMed  CAS  Google Scholar 

  31. Uzan M, Miller ES: Post-transcriptional control by bacteriophage T4: mRNA decay and inhibition of translation initiation. Virol J. 2010, 7 (1): 360-10.1186/1743-422X-7-360.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  32. Klaiman D, Steinfels-Kohn E, Krutkina E, Davidov E, Kaufmann G: The wobble nucleotide-excising anticodon nuclease RloC is governed by the zinc-hook and DNA-dependent ATPase of its Rad50-like region. Nucleic Acids Res. 2012, 40 (17): 8568-8578. 10.1093/nar/gks593.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  33. Davidov E, Kaufmann G: RloC: a wobble nucleotide-excising and zinc-responsive bacterial tRNase. Mol Microbiol. 2008, 69 (6): 1560-1574. 10.1111/j.1365-2958.2008.06387.x.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  34. Snyder L: Phage-exclusion enzymes: a bonanza of biochemical and cell biology reagents?. Mol Microbiol. 1995, 15 (3): 415-420. 10.1111/j.1365-2958.1995.tb02255.x.

    Article  PubMed  CAS  Google Scholar 

  35. Aravind L, Anantharaman V, Zhang D, de Souza RF, Iyer LM: Gene flow and biological conflict systems in the origin and evolution of eukaryotes. Front Cell Infect Microbiol. 2012, 2: 89-

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  36. Sashital DG, Wiedenheft B, Doudna JA: Mechanism of foreign DNA selection in a bacterial adaptive immune system. Mol Cell. 2012, 46 (5): 606-615. 10.1016/j.molcel.2012.03.020.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  37. Jinek M, Chylinski K, Fonfara I, Hauer M, Doudna JA, Charpentier E: A programmable dual-RNA-guided DNA endonuclease in adaptive bacterial immunity. Science. 2012, 337 (6096): 816-821. 10.1126/science.1225829.

    Article  PubMed  CAS  Google Scholar 

  38. Datsenko KA, Pougach K, Tikhonov A, Wanner BL, Severinov K, Semenova E: Molecular memory of prior infections activates the CRISPR/Cas adaptive bacterial immunity system. Nat Commun. 2012, 3: 945-

    Article  PubMed  Google Scholar 

  39. Wiedenheft B, van Duijn E, Bultema J, Waghmare S, Zhou K, Barendregt A, Westphal W, Heck A, Boekema E, Dickman M, et al: RNA-guided complex from a bacterial immune system enhances target recognition through seed sequence interactions. Proc Natl Acad Sci U S A. 2011, 108 (25): 10092-10097. 10.1073/pnas.1102716108.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  40. Wiedenheft B, Lander GC, Zhou K, Jore MM, Brouns SJ, van der Oost J, Doudna JA, Nogales E: Structures of the RNA-guided surveillance complex from a bacterial immune system. Nature. 2011, 477 (7365): 486-489. 10.1038/nature10402.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  41. Terns MP, Terns RM: CRISPR-based adaptive immune systems. Curr Opin Microbiol. 2011, 14 (3): 321-327. 10.1016/j.mib.2011.03.005.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  42. Sinkunas T, Gasiunas G, Fremaux C, Barrangou R, Horvath P, Siksnys V: Cas3 is a single-stranded DNA nuclease and ATP-dependent helicase in the CRISPR/Cas immune system. EMBO J. 2011, 30 (7): 1335-1342. 10.1038/emboj.2011.41.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  43. Semenova E, Jore MM, Datsenko KA, Semenova A, Westra ER, Wanner B, van der Oost J, Brouns SJ, Severinov K: Interference by clustered regularly interspaced short palindromic repeat (CRISPR) RNA is governed by a seed sequence. Proc Natl Acad Sci U S A. 2011, 108 (25): 10098-10103. 10.1073/pnas.1104144108.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  44. Yosef I, Goren MG, Qimron U: Proteins and DNA elements essential for the CRISPR adaptation process in Escherichia coli. Nucleic Acids Res. 2012, 40 (12): 5569-5576. 10.1093/nar/gks216.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  45. Makarova KS, Haft DH, Barrangou R, Brouns SJJ, Charpentier E, Horvath P, Moineau S, Mojica FJM, Wolf YI, Yakunin AF, et al: Evolution and classification of the CRISPR/Cas systems. Nature Rev Microbiol. 2011, 9 (6): 467-477. 10.1038/nrmicro2577.

    Article  CAS  Google Scholar 

  46. Roberts RJ, Vincze T, Posfai J, Macelis D: REBASE--enzymes and genes for DNA restriction and modification. Nucleic Acids Res. 2007, 35 (Database issue): D269-D270.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  47. Beloglazova N, Brown G, Zimmerman MD, Proudfoot M, Makarova KS, Kudritska M, Kochinyan S, Wang S, Chruszcz M, Minor W, et al: A novel family of sequence-specific endoribonucleases associated with the clustered regularly interspaced short palindromic repeats. J Biol Chem. 2008, 283 (29): 20361-20371. 10.1074/jbc.M803225200.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  48. Makarova KS, Aravind L, Wolf YI, Koonin EV: Unification of Cas protein families and a simple scenario for the origin and evolution of CRISPR-Cas systems. Biol Direct. 2011, 6 (1): 38-10.1186/1745-6150-6-38.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  49. Makarova KS, Grishin NV, Shabalina SA, Wolf YI, Koonin EV: A putative RNA-interference-based immune system in prokaryotes: computational analysis of the predicted enzymatic machinery, functional analogies with eukaryotic RNAi, and hypothetical mechanisms of action. Biol Direct. 2006, 1 (1): 7-10.1186/1745-6150-1-7.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Kwon AR, Kim JH, Park SJ, Lee KY, Min YH, Im H, Lee I, Lee BJ: Structural and biochemical characterization of HP0315 from Helicobacter pylori as a VapD protein with an endoribonuclease activity. Nucleic Acids Res. 2012, 40 (9): 4216-4228. 10.1093/nar/gkr1305.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  51. Takeuchi N, Wolf YI, Makarova KS, Koonin EV: Nature and intensity of selection pressure on CRISPR-associated genes. J Bacteriol. 2012, 194 (5): 1216-1225. 10.1128/JB.06521-11.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  52. Ulferts R, Ziebuhr J: Nidovirus ribonucleases: structures and functions in viral replication. RNA Biol. 2011, 8 (2): 295-304. 10.4161/rna.8.2.15196.

    Article  PubMed  CAS  Google Scholar 

  53. Grynberg M, Erlandsen H, Godzik A: HEPN: a common domain in bacterial drug resistance and human neurodegenerative proteins. Trends Biochem Sci. 2003, 28 (5): 224-226. 10.1016/S0968-0004(03)00060-4.

    Article  PubMed  CAS  Google Scholar 

  54. Kimelman A, Levy A, Sberro H, Kidron S, Leavitt A, Amitai G, Yoder-Himes DR, Wurtzel O, Zhu Y, Rubin EM, et al: A vast collection of microbial genes that are toxic to bacteria. Genome Res. 2012, 22 (4): 802-809. 10.1101/gr.133850.111.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  55. Anantharaman V, Aravind L: New connections in the prokaryotic toxin-antitoxin network: relationship with the eukaryotic nonsense-mediated RNA decay system. Genome Biol. 2003, 4 (12): R81-10.1186/gb-2003-4-12-r81.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Kojima KK, Kanehisa M: Systematic survey for novel types of prokaryotic retroelements based on gene neighborhood and protein architecture. Mol Biol Evol. 2008, 25 (7): 1395-1404. 10.1093/molbev/msn081.

    Article  PubMed  CAS  Google Scholar 

  57. Wang C, Villion M, Semper C, Coros C, Moineau S, Zimmerly S: A reverse transcriptase-related protein mediates phage resistance and polymerizes untemplated DNA in vitro. Nucleic Acids Res. 2011, 39 (17): 7620-7629. 10.1093/nar/gkr397.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  58. Zhang D, de Souza RF, Anantharaman V, Iyer LM, Aravind L: Polymorphic toxin systems: comprehensive characterization of trafficking modes, processing, mechanisms of action, immunity and ecology using comparative genomics. Biol Direct. 2012, 7 (1): 18-10.1186/1745-6150-7-18.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  59. Makarova KS, Wolf YI, van der Oost J, Koonin EV: Prokaryotic homologs of Argonaute proteins are predicted to function as key components of a novel system of defense against mobile genetic elements. Biol Direct. 2009, 4: 29-10.1186/1745-6150-4-29.

    Article  PubMed  PubMed Central  Google Scholar 

  60. de Souza RF, Aravind L: Identification of novel components of NAD-utilizing metabolic pathways and prediction of their biochemical functions. Mol Biosyst. 2012, 8 (6): 1661-1677. 10.1039/c2mb05487f.

    Article  PubMed  Google Scholar 

  61. Iyer LM, Makarova KS, Koonin EV, Aravind L: Comparative genomics of the FtsK-HerA superfamily of pumping ATPases: implications for the origins of chromosome segregation, cell division and viral capsid packaging. Nucleic Acids Res. 2004, 32 (17): 5260-5279. 10.1093/nar/gkh828.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  62. Anantharaman V, Iyer LM, Aravind L: Ter-dependent stress response systems: novel pathways related to metal sensing, production of a nucleoside-like metabolite, and DNA-processing. Mol Biosyst. 2012, in press

    Google Scholar 

  63. Schultz D, Wolynes PG, Ben Jacob E, Onuchic JN: Deciding fate in adverse times: sporulation and competence in Bacillus subtilis. Proc Natl Acad Sci U S A. 2009, 106 (50): 21027-21034. 10.1073/pnas.0912185106.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  64. Jabbari S, Heap JT, King JR: Mathematical modelling of the sporulation-initiation network in Bacillus subtilis revealing the dual role of the putative quorum-sensing signal molecule PhrA. Bull Math Biol. 2011, 73 (1): 181-211. 10.1007/s11538-010-9530-7.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Yuri Wolf and Tatiana Senkevich for critical reading of the manuscript and useful discussions. The authors’ research is supported by the US Department of Health and Human Services intramural funds (to National Library of Medicine).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Eugene V Koonin.

Additional information

Competing interest

The authors declare no competing interest.

Authors’ contributions

KSM and VA collected and analyzed the data; KSM, LA and EVK formulated the hypothesis and wrote the manuscript that was read and approved by all authors.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Makarova, K.S., Anantharaman, V., Aravind, L. et al. Live virus-free or die: coupling of antivirus immunity and programmed suicide or dormancy in prokaryotes. Biol Direct 7, 40 (2012). https://doi.org/10.1186/1745-6150-7-40

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1745-6150-7-40

Keywords