Bacterial Genome Evolution Research Articles

MSphere of Influence: How I Learned to Love Bacteria and their Tangled Evolutionary Tree.

ABSTRACTTera Levin works in the fields of evolution, microbiology, and genetics, studying how adaptation shapes the molecular interactions between eukaryotic hosts and bacterial pathogens. In this mSphere of Influence article, she reflects on how the paper “Population genomics of early events in the ecological differentiation of bacteria” by Shapiro et al. (B. J. Shapiro, J. Friedman, O. X. Cordero, S. P. Preheim, et al., Science 336:48–51, 2012, https://doi.org/10.1126/science.1218198) changed the way she thinks about bacterial gene and genome evolution.

Assessing Biosynthetic Gene Cluster Diversity of Specialized Metabolites in the Conserved Gut Symbionts of Herbivorous Turtle Ants.

Cephalotes are herbivorous ants (>115 species) feeding on low-nitrogen food sources, and they rely on gut symbionts to supplement their diet by recycling nitrogen food waste into amino acids. These conserved gut symbionts, which encompass five bacterial orders, have been studied previously for their primary nitrogen metabolism; however, little is known about their ability to biosynthesize specialized metabolites which can play a role in bacterial interactions between communities living in close proximity in the gut. To evaluate the biosynthetic potential of their gut symbionts, we mine 14 cultured isolate genomes and gut metagenomes across 17 Cephalotes species to explore the biodiversity of biosynthetic gene clusters (BGCs) producing specialized metabolites. The diversity of BGCs across Cephalotes phylogeny was analyzed using sequence similarity networking and BGC phylogenetic reconstruction. Our results reveal that the conserved gut symbionts involved in the nutritional symbiosis possess 80% of all the 233 BGCs retrieved in this work. Furthermore, the phylogenetic analysis of BGCs reveals different patterns of distribution, suggesting different mechanisms of conservation. A siderophore BGC shows high similarity in a single symbiont across different ant host species, whereas a BGC encoding the production of non-ribosomal peptides (NRPs) found different symbionts within a single host species. Additionally, BGCs were abundant in four of the five bacterial orders of conserved symbionts co-occurring in the hindgut. However, one major symbiont localized alone in the midgut lack BGCs. Because the spatial isolation prevents direct interaction with other symbionts, this result supports the idea that BGCs are maintained in bacteria living in close proximity but are dispensable for an alone-living symbiont. These findings together pave the way for studying the mechanisms of BGC conservation and evolution in gut bacterial genomes associated with Cephalotes. This work also provides a genetic background for further study, aiming to characterize bacterial specialized metabolites and to understand their functional role in multipartite mutualisms between conserved gut symbionts and Cephalotes turtle ants.

Xenogeneic Silencing and Bacterial Genome Evolution: Mechanisms for DNA Recognition Imply Multifaceted Roles of Xenogeneic Silencers.

Horizontal gene transfer (HGT) is a major driving force for bacterial evolution. To avoid the deleterious effects due to the unregulated expression of newly acquired foreign genes, bacteria have evolved specific proteins named xenogeneic silencers to recognize foreign DNA sequences and suppress their transcription. As there is considerable diversity in genomic base compositions among bacteria, how xenogeneic silencers distinguish self- from nonself DNA in different bacteria remains poorly understood. This review summarizes the progress in studying the DNA binding preferences and the underlying molecular mechanisms of known xenogeneic silencer families, represented by H-NS of Escherichia coli, Lsr2 of Mycobacterium, MvaT of Pseudomonas, and Rok of Bacillus. Comparative analyses of the published data indicate that the differences in DNA recognition mechanisms enable these xenogeneic silencers to have clear characteristics in DNA sequence preferences, which are further correlated with different host genomic features. These correlations provide insights into the mechanisms of how these xenogeneic silencers selectively target foreign DNA in different genomic backgrounds. Furthermore, it is revealed that the genomic AT contents of bacterial species with the same xenogeneic silencer family proteins are distributed in a limited range and are generally lower than those species without any known xenogeneic silencers in the same phylum/class/genus, indicating that xenogeneic silencers have multifaceted roles on bacterial genome evolution. In addition to regulating horizontal gene transfer, xenogeneic silencers also act as a selective force against the GC to AT mutational bias found in bacterial genomes and help the host genomic AT contents maintained at relatively low levels.

Evolutionary Rates are Correlated Between Buchnera Endosymbionts andthe Mitochondrial Genomes of Their Aphid Hosts.

The evolution of bacterial endosymbiont genomes is strongly influenced by host-driven selection. Factors affecting host genome evolution will potentially affect endosymbiont genomes in similar ways. One potential outcome is correlations in molecular rates between the genomes of the symbiotic partners. Recently, we presented the first evidence of such correlations between the mitochondrial genomes of cockroaches and the genomes of their endosymbiont (Blattabacterium cuenoti). Here we investigate whether similar patterns are found in additional host-symbiont partners. We use partial genome data from multiple strains of the bacterial endosymbionts Buchnera aphidicola and Sulcia muelleri, and the mitochondrial genomes of their sap-feeding insect hosts. Both endosymbionts show phylogenetic congruence with the mitochondria of their hosts, a result that is expected due to their identical mode of inheritance. We compared root-to-tip distances and branch lengths of phylogenetically independent species pairs. Both analyses showed a highly significant correlation of molecular rates between thegenomes of Buchnera and themitochondrial genomes of their hosts. A similar correlation was detected between Sulcia and their hosts, but was not statistically significant. Our results indicate that evolutionary rate correlations between hosts and long-term symbionts may be a widespread phenomenon.

Genomic and molecular evolutionary dynamics of transcriptional response regulator genes in bacterial species of the Harveyi clade of Vibrio.

Transcriptional response regulators (TRR) are the most abundant signal transducers in prokaryotic systems that mediate intracellular changes in response to environmental signals. They are involved in a wide range of biological processes that allow bacteria to persist in particular habitats. There is strong evidence that the bacterial habitat and their lifestyle influence the size of their TRR genetic repertoire. Therefore, it would be expected that the evolution of bacterial genomes could be linked to natural selection processes. To test this hypothesis, we explored the evolutionary dynamics of TRR genes of the widely studied Harveyi clade of the genus Vibrio at the molecular and genomic levels. Our results suggest that the TRR genetic repertoire of the species belonging to the Harveyi clade is a product of genomic reduction and expansion. The gene loss and gains that drive their genomic reduction and expansion could be attributed to natural selection and random genetic drift. It seems that natural selection acts to maintain the ancestral state of core TRR genes (shared by all species) by purifying processes and could be driving the loss of some accessory (found in certain species) genes through the diversification of sequences. The neutrality observed in gene gain could be attributed to spontaneous events as horizontal gene transfer driven by stochastic events as occurs in random genetic drift.

Effect of arsenite and growth in biofilm conditions on the evolution of Thiomonas sp. CB2.

Thiomonas bacteria are ubiquitous at acid mine drainage sites and play key roles in the remediation of water at these locations by oxidizing arsenite to arsenate, favouring the sorption of arsenic by iron oxides and their coprecipitation. Understanding the adaptive capacities of these bacteria is crucial to revealing how they persist and remain active in such extreme conditions. Interestingly, it was previously observed that after exposure to arsenite, when grown in a biofilm, some strains of Thiomonas bacteria develop variants that are more resistant to arsenic. Here, we identified the mechanisms involved in the emergence of such variants in biofilms. We found that the percentage of variants generated increased in the presence of high concentrations of arsenite (5.33 mM), especially in the detached cells after growth under biofilm-forming conditions. Analysis of gene expression in the parent strain CB2 revealed that genes involved in DNA repair were upregulated in the conditions where variants were observed. Finally, we assessed the phenotypes and genomes of the subsequent variants generated to evaluate the number of mutations compared to the parent strain. We determined that multiple point mutations accumulated after exposure to arsenite when cells were grown under biofilm conditions. Some of these mutations were found in what is referred to as ICE19, a genomic island (GI) carrying arsenic-resistance genes, also harbouring characteristics of an integrative and conjugative element (ICE). The mutations likely favoured the excision and duplication of this GI. This research aids in understanding how Thiomonas bacteria adapt to highly toxic environments, and, more generally, provides a window to bacterial genome evolution in extreme environments.

The roles of replication-transcription conflict in mutagenesis and evolution of genome organization.

Replication-transcription conflicts promote mutagenesis and give rise to evolutionary signatures, with fundamental importance to genome stability ranging from bacteria to metastatic cancer cells. This review focuses on the interplay between replication-transcription conflicts and the evolution of gene directionality. In most bacteria, the majority of genes are encoded on the leading strand of replication such that their transcription is co-directional with the direction of DNA replication fork movement. This gene strand bias arises primarily due to negative selection against deleterious consequences of head-on replication-transcription conflict. However, many genes remain head-on. Can head-on orientation provide some benefit? We combine insights from both mechanistic and evolutionary studies, review published work, and analyze gene expression data to evaluate an emerging model that head-on genes are temporal targets for adaptive mutagenesis during stress. We highlight the alternative explanation that genes in the head-on orientation may simply be the result of genomic inversions and relaxed selection acting on nonessential genes. We seek to clarify how the mechanisms of replication-transcription conflict, in concert with other mutagenic mechanisms, balanced by natural selection, have shaped bacterial genome evolution.

Palaeomicrobiology: Application of Ancient DNA Sequencing to Better Understand Bacterial Genome Evolution and Adaptation

Next Generation Sequencing (NGS) has unlocked access to the wide range of non-cultivable microorganisms, including those present in the ancient past. The study of microorganisms from ancient sources (palaeomicrobiology) using DNA sequencing now provides a unique opportunity to examine ancient microbial genomic content, explore pathogenicity, and understand microbial evolution in greater detail than ever before. As a result, current studies have focused on reconstructing the evolutionary history of a number of human pathogens involved in ancient and historic pandemic events. These studies have opened the door for a variety of future palaeomicrobiology studies, which can focus on commensal microorganisms, species from non-human hosts, information from host-genomics, and the use of bacteria as proxies for additional information about past human health, behavior, migration, and culture. Here, we describe the origin and the historical and recent advances in the field of palaeomicrobiology, review some of the most notable ancient pathogenic microorganism studies, and provide perspectives on how NGS and whole genome information from ancient microorganisms contribute to our understanding of bacterial evolution on a broader scale. We conclude by exploring the application of newly developed tools in palaeomicrobiology and discussing how future studies can improve our current understanding of the co-evolution relationships between humans and non-pathogenic microbes.

Genome-Wide Analysis Reveals Genetic Potential for Aromatic Compounds Biodegradation of Sphingopyxis.

Members of genus Sphingopyxis are frequently found in diverse eco-environments worldwide and have been traditionally considered to play vital roles in the degradation of aromatic compounds. Over recent decades, many aromatic-degrading Sphingopyxis strains have been isolated and recorded, but little is known about their genetic nature related to aromatic compounds biodegradation. In this study, bacterial genomes of 19 Sphingopyxis strains were used for comparative analyses. Phylogeny showed an ambiguous relatedness between bacterial strains and their habitat specificity, while clustering based on Cluster of Orthologous Groups suggested the potential link of functional profile with substrate-specific traits. Pan-genome analysis revealed that 19 individuals were predicted to share 1,066 orthologous genes, indicating a high genetic homogeneity among Sphingopyxis strains. Notably, KEGG Automatic Annotation Server results suggested that most genes pertaining aromatic compounds biodegradation were predicted to be involved in benzoate, phenylalanine, and aminobenzoate metabolism. Among them, β-ketoadipate biodegradation might be the main pathway in Sphingopyxis strains. Further inspection showed that a number of mobile genetic elements varied in Sphingopyxis genomes, and plasmid-mediated gene transfer coupled with prophage- and transposon-mediated rearrangements might play prominent roles in the evolution of bacterial genomes. Collectively, our findings presented that Sphingopyxis isolates might be the promising candidates for biodegradation of aromatic compounds in pollution sites.

Markov-Modulated Continuous-Time Markov Chains to Identify Site- and Branch-Specific Evolutionary Variation in BEAST

Markov models of character substitution on phylogenies form the foundation of phylogenetic inference frameworks. Early models made the simplifying assumption that the substitution process is homogeneous over time and across sites in the molecular sequence alignment. While standard practice adopts extensions that accommodate heterogeneity of substitution rates across sites, heterogeneity in the process over time in a site-specific manner remains frequently overlooked. This is problematic, as evolutionary processes that act at the molecular level are highly variable, subjecting different sites to different selective constraints over time, impacting their substitution behavior. We propose incorporating time variability through Markov-modulated models (MMMs), which extend covarion-like models and allow the substitution process (including relative character exchange rates as well as the overall substitution rate) at individual sites to vary across lineages. We implement a general MMM framework in BEAST, a popular Bayesian phylogenetic inference software package, allowing researchers to compose a wide range of MMMs through flexible XML specification. Using examples from bacterial, viral, and plastid genome evolution, we show that MMMs impact phylogenetic tree estimation and can substantially improve model fit compared to standard substitution models. Through simulations, we show that marginal likelihood estimation accurately identifies the generative model and does not systematically prefer the more parameter-rich MMMs. To mitigate the increased computational demands associated with MMMs, our implementation exploits recent developments in BEAGLE, a high-performance computational library for phylogenetic inference. [Bayesian inference; BEAGLE; BEAST; covarion, heterotachy; Markov-modulated models; phylogenetics.]

Deciphering the role of insertion sequences in the evolution of bacterial epidemic pathogens with panISa software.

Next-generation sequencing (NGS) is now widely used in microbiology to explore genome evolution and the structure of pathogen outbreaks. Bioinformatics pipelines readily detect single-nucleotide polymorphisms or short indels. However, bacterial genomes also evolve through the action of small transposable elements called insertion sequences (ISs), which are difficult to detect due to their short length and multiple repetitions throughout the genome. We designed panISa software for the ab initio detection of IS insertions in the genomes of prokaryotes. PanISa has been released as open source software (GPL3) available from https://github.com/bvalot/panISa. In this study, we assessed the utility of this software for evolutionary studies, by reanalysing five published datasets for outbreaks of human major pathogens in which ISs had not been specifically investigated. We reanalysed the raw data from each study, by aligning the reads against reference genomes and running panISa on the alignments. Each hit was automatically curated and IS-related events were validated on the basis of nucleotide sequence similarity, by comparison with the ISFinder database. In Acinetobacter baumannii , the panISa pipeline identified ISAba1 or ISAba125 upstream from the ampC gene, which encodes a cephalosporinase in all third-generation cephalosporin-resistant isolates. In the genomes of Vibrio cholerae isolates, we found that early Haitian isolates had the same ISs as Nepalese isolates, confirming the inferred history of the contamination of this island. In Enterococcus faecalis , panISa identified regions of high plasticity, including a pathogenicity island enriched in IS-related events. The overall distribution of ISs deduced with panISa was consistent with SNP-based phylogenic trees, for all species considered. The role of ISs in pathogen evolution has probably been underestimated due to difficulties detecting these transposable elements. We show here that panISa is a useful addition to the bioinformatics toolbox for analyses of the evolution of bacterial genomes. PanISa will facilitate explorations of the functional impact of ISs and improve our understanding of prokaryote evolution.

Horizontally Acquired Homologs of Xenogeneic Silencers: Modulators of Gene Expression Encoded by Plasmids, Phages and Genomic Islands.

Acquisition of mobile elements by horizontal gene transfer can play a major role in bacterial adaptation and genome evolution by providing traits that contribute to bacterial fitness. However, gaining foreign DNA can also impose significant fitness costs to the host bacteria and can even produce detrimental effects. The efficiency of horizontal acquisition of DNA is thought to be improved by the activity of xenogeneic silencers. These molecules are a functionally related group of proteins that possess affinity for the acquired DNA. Binding of xenogeneic silencers suppresses the otherwise uncontrolled expression of genes from the newly acquired nucleic acid, facilitating their integration to the bacterial regulatory networks. Even when the genes encoding for xenogeneic silencers are part of the core genome, homologs encoded by horizontally acquired elements have also been identified and studied. In this article, we discuss the current knowledge about horizontally acquired xenogeneic silencer homologs, focusing on those encoded by genomic islands, highlighting their distribution and the major traits that allow these proteins to become part of the host regulatory networks.

Recording mobile DNA in the gut microbiota using an Escherichia coli CRISPR-Cas spacer acquisition platform

The flow of genetic material between bacteria is central to the adaptation and evolution of bacterial genomes. However, our knowledge about DNA transfer within complex microbiomes is lacking, with most studies of horizontal gene transfer (HGT) relying on bioinformatic analyses of genetic elements maintained on evolutionary timescales or experimental measurements of phenotypically trackable markers. Here, we utilize the CRISPR-Cas spacer acquisition process to detect DNA acquisition events from complex microbiota in real-time and at nucleotide resolution. In this system, an E. coli recording strain is exposed to a microbial sample and spacers are acquired from transferred plasmids and permanently stored in genomic CRISPR arrays. Sequencing and analysis of acquired spacers enables identification of the transferred plasmids. This approach allowed us to identify individual mobile elements without relying on phenotypic markers or post-transfer replication. We found that HGT into the recording strain in human clinical fecal samples can be extensive and is driven by different plasmid types, with the IncX type being the most actively transferred.

Characteristics of two myoviruses induced from the coastal photoheterotrophic bacterium Porphyrobacter sp. YT40.

In this study, we characterized two induced myoviruses from one marine photoheterotrophic bacterium Porphyrobacter sp. YT40 belonging to the Sphingomonadales family in Alphaproteobacteria. The genome sequence of prophage A is ∼36.9kb with an average GC content of 67.1%, and its core or functional genes are homologous to Mu or Mu-like phages. Furthermore, induced viral particles from prophage A show a knob-like neck structure, which is only found in bacteriophage Mu. The genome size of prophage B is ∼36.8kb with an average GC content of 65.3%. Prophage B contains a conserved gene cluster Q-P-O-N-M-L, which is unique in P2 phages. Induced viral particles from prophage B display an icosahedral head with a diameter of ∼55nm and a 130±5nm long contractile tail. To our knowledge, this is the first report that characterizes the induced P2-like phage in marine Alphaproteobacteria. Phylogeny analyses suggest that these two types of prophages are commonly found in sequenced bacteria of the Sphingomonadales family. This study sheds light on the ongoing interaction between marine bacteria and phages, and improves our understanding of bacterial genomic plasticity and evolution.

IslandCafe: Compositional Anomaly and Feature Enrichment Assessment for Delineation of Genomic Islands.

One of the evolutionary forces driving bacterial genome evolution is the acquisition of clusters of genes through horizontal gene transfer (HGT). These genomic islands may confer adaptive advantages to the recipient bacteria, such as, the ability to thwart antibiotics, become virulent or hypervirulent, or acquire novel metabolic traits. Methods for detecting genomic islands either search for markers or features typical of islands or examine anomaly in oligonucleotide composition against the genome background. The former tends to underestimate, missing islands that have the markers either lost or degraded, while the latter tends to overestimate, due to their inability to discriminate compositional atypicality arising because of HGT from those that are a consequence of other biological factors. We propose here a framework that exploits the strengths of both these approaches while bypassing the pitfalls of either. Genomic islands lacking markers are identified by their association with genomic islands with markers. This was made possible by performing marker enrichment and phyletic pattern analyses within an integrated framework of recursive segmentation and clustering. The proposed method, IslandCafe, compared favorably with frequently used methods for genomic island detection on synthetic test datasets and on a test-set of known islands from 15 well-characterized bacterial species. Furthermore, IslandCafe identified novel islands with imprints of likely horizontal acquisition.

Global-level population genomics reveals differential effects of geography and phylogeny on horizontal gene transfer in soil bacteria

Although microorganisms are known to dominate Earth's biospheres and drive biogeochemical cycling, little is known about the geographic distributions of microbial populations or the environmental factors that pattern those distributions. We used a global-level hierarchical sampling scheme to comprehensively characterize the evolutionary relationships and distributional limitations of the nitrogen-fixing bacterial symbionts of the crop chickpea, generating 1,027 draft whole-genome sequences at the level of bacterial populations, including 14 high-quality PacBio genomes from a phylogenetically representative subset. We find that diverse Mesorhizobium taxa perform symbiosis with chickpea and have largely overlapping global distributions. However, sampled locations cluster based on the phylogenetic diversity of Mesorhizobium populations, and diversity clusters correspond to edaphic and environmental factors, primarily soil type and latitude. Despite long-standing evolutionary divergence and geographic isolation, the diverse taxa observed to nodulate chickpea share a set of integrative conjugative elements (ICEs) that encode the major functions of the symbiosis. This symbiosis ICE takes 2 forms in the bacterial chromosome-tripartite and monopartite-with tripartite ICEs confined to a broadly distributed superspecies clade. The pairwise evolutionary relatedness of these elements is controlled as much by geographic distance as by the evolutionary relatedness of the background genome. In contrast, diversity in the broader gene content of Mesorhizobium genomes follows a tight linear relationship with core genome phylogenetic distance, with little detectable effect of geography. These results illustrate how geography and demography can operate differentially on the evolution of bacterial genomes and offer useful insights for the development of improved technologies for sustainable agriculture.

Clustered Core- and Pan-Genome Content on Rhodobacteraceae Chromosomes.

In Bacteria, chromosome replication starts at a single origin of replication and proceeds on both replichores. Due to its asymmetric nature, replication influences chromosome structure and gene organization, mutation rate, and expression. To date, little is known about the distribution of highly conserved genes over the bacterial chromosome. Here, we used a set of 101 fully sequenced Rhodobacteraceae representatives to analyze the relationship between conservation of genes within this family and their distance from the origin of replication. Twenty-two of the analyzed species had core genes clustered significantly closer to the origin of replication with representatives of the genus Celeribacter being the most apparent example. Interestingly, there were also eight species with the opposite organization. In particular, Rhodobaca barguzinensis and Loktanella vestfoldensis showed a significant increase of core genes with distance from the origin of replication. The uneven distribution of low-conserved regions is in particular pronounced for genomes in which the halves of one replichore differ in their conserved gene content. Phage integration and horizontal gene transfer partially explain the scattered nature of Rhodobacteraceae genomes. Our findings lay the foundation for a better understanding of bacterial genome evolution and the role of replication therein.

Complete Sequences of Two Plasmids Found in a Brazilian Bacillus thuringiensis Serovar israelensis Strain.

Plasmids play a crucial role in the evolution of bacterial genomes by mediating horizontal gene transfer. In this work, we sequenced two plasmids found in a Brazilian Bacillus thuringiensis serovar israelensis strain which showed 100% nucleotide identities with Bacillus thuringiensis serovar kurstaki plasmids.

Bacterial 'Grounded' Prophages: Hotspots for Genetic Renovation and Innovation.

Bacterial genomes are highly plastic allowing the generation of variants through mutations and acquisition of genetic information. The fittest variants are then selected by the econiche thereby allowing the bacterial adaptation and colonization of the habitat. Larger genomes, however, may impose metabolic burden and hence bacterial genomes are optimized by the loss of frivolous genetic information. The activity of temperate bacteriophages has acute consequences on the bacterial population as well as the bacterial genome through lytic and lysogenic cycles. Lysogeny is a selective advantage as the prophage provides immunity to the lysogen against secondary phage attack. Since the non-lysogens are eliminated by the lytic phages, lysogens multiply and colonize the habitat. Nevertheless, all lysogens have an imminent risk of lytic cycle activation and cell lysis. However, a mutation in the attachment sites or in the genes that encode the specific recombinase responsible for prophage excision could result in ‘grounding’ of the prophage. Since the lysogens with grounded prophage are immune to respective phage infection as well as dodge the induction of lytic cycle, we hypothesize that the selection of these mutant lysogens is favored relative to their normal lysogenic counterparts. These grounded prophages offer several advantages to the bacterial genome evolution through propensity for genetic variations including inversions, deletions, and insertions via horizontal gene transfer. We propose that the grounded prophages expedite bacterial genome evolution by acting as ‘genetic buffer zones’ thereby increasing the frequency as well as the diversity of variations on which natural selection favors the beneficial variants. The grounded prophages are also hotspots for horizontal gene transfer wherein several ecologically significant genes such as those involved in stress tolerance, antimicrobial resistance, and novel metabolic pathways, are integrated. Moreover, the high frequency of genetic changes within prophages also allows proportionate probability for the de novo genesis of genetic information. Through sequence analyses of well-characterized E. coli prophages we exemplify various roles of grounded prophages in E. coli ecology and evolution. Therefore, the temperate prophages are one of the most significant drivers of bacterial genome evolution and sites of biogenesis of genetic information.

Lamarckian realities: the CRISPR-Cas system and beyond

In his target article, Koonin discusses the insights into the evolution of bacterial genomes provided by the CRISPR-Cas system. This evolved defense system is based on intrinsic processes of genome engineering which, as he argues, enable Lamarckian inheritance. In this commentary I discuss some historical and conceptual issues that pertain to Koonin’s analysis of this aspect of the CRISPR-Cas system, extending and qualifying his discussion.