Expanding the catalog of cas genes with metagenomes

The CRISPR (clusters of regularly interspaced short palindromic repeats)–Cas adaptive immune system is an important defense system in bacteria, providing targeted defense against invasions of foreign nucleic acids. CRISPR–Cas systems consist of CRISPR loci and cas (CRISPR-associated) genes: sequence segments of invaders are incorporated into host genomes at CRISPR loci to generate specificity, while adjacent cas genes encode proteins that mediate the defense process. We pursued an integrated approach to identifying putative cas genes from genomes and metagenomes, combining similarity searches with genomic neighborhood analysis. Application of our approach to bacterial genomes and human microbiome datasets allowed us to significantly expand the collection of cas genes: the sequence space of the Cas9 family, the key player in the recently engineered RNA-guided platforms for genome editing in eukaryotes, is expanded by at least two-fold with metagenomic datasets. We found genes in cas loci encoding other functions, for example, toxins and antitoxins, confirming the recently discovered potential of coupling between adaptive immunity and the dormancy/suicide systems. We further identified 24 novel Cas families; one novel family contains 20 proteins, all identified from the human microbiome datasets, illustrating the importance of metagenomics projects in expanding the diversity of cas genes.     

A PNPase Dependent CRISPR System in Listeria

The human bacterial pathogen Listeria monocytogenes is emerging as a model organism to study RNA-mediated regulation in pathogenic bacteria. A class of non-coding RNAs called CRISPRs (clustered regularly interspaced short palindromic repeats) has been described to confer bacterial resistance against invading bacteriophages and conjugative plasmids. CRISPR function relies on the activity of CRISPR associated (cas) genes that encode a large family of proteins with nuclease or helicase activities and DNA and RNA binding domains. Here, we characterized a CRISPR element (RliB) that is expressed and processed in the L. monocytogenes strain EGD-e, which is completely devoid of cas genes. Structural probing revealed that RliB has an unexpected secondary structure comprising basepair interactions between the repeats and the adjacent spacers in place of canonical hairpins formed by the palindromic repeats. Moreover, in contrast to other CRISPR-Cas systems identified in Listeria, RliB-CRISPR is ubiquitously present among Listeria genomes at the same genomic locus and is never associated with the cas genes. We showed that RliB-CRISPR is a substrate for the endogenously encoded polynucleotide phosphorylase (PNPase) enzyme. The spacers of the different Listeria RliB-CRISPRs share many sequences with temperate and virulent phages. Furthermore, we show that a cas-less RliB-CRISPR lowers the acquisition frequency of a plasmid carrying the matching protospacer, provided that trans encoded cas genes of a second CRISPR-Cas system are present in the genome. Importantly, we show that PNPase is required for RliB-CRISPR mediated DNA interference. Altogether, our data reveal a yet undescribed CRISPR system whose both processing and activity depend on PNPase, highlighting a new and unexpected function for PNPase in “CRISPRology”.     

SSO1450 - A CAS1 protein from Sulfolobus solfataricus P2 with high affinity for RNA and DNA

Clustered regularly interspaced short palindromic repeats (CRISPR) and their associated protein genes (cas genes) are ubiquitous in archaea and eubacteria. It has been suggested that CRISPR and CAS proteins act as an immune system preventing the invasion of foreign genomic elements at the DNA level. The protein SSO1450 from Sulfolobus solfataricus (Sso) P2 belongs to the CAS1 cluster which is one of the core protein clusters most frequently associated with CRISPR sequences. In this study we show that SSO1450 is a high-affinity nucleic acid binding protein. It binds DNA, RNA and DNA-RNA hybrid apparently sequence non-specific in a multi-site binding mode. Furthermore, SSO1450 promotes the hybridization of complementary nucleic acid strands.     

Molecular identification and characterization of clustered regularly interspaced short palindromic repeats (CRISPRs) in Campylobacter lari

PCR amplifications using primers for the clustered regularly interspaced short palindromic repeats (CRISPRs)-associated gene 1 (cas1), cas2, putative (p)-cas and CRISPRs genes generated cas1, cas2, p-cas and CRISPRs genes segments with 9–28 of 28 urease-positive thermophilic Campylobacter (UPTC) isolates, respectively. The p-cas and CRISPRs genes segments were amplified with 10 of 11 and 0 of 11 urease-negative (UN) Campylobacter lari isolates, respectively. When the nucleotide sequences of the CRISPRs consensus sequence repeats of each 33–37 base pairs from the 18 Campylobacter jejuni isolates were aligned, as well as from the four C. jejuni reference and UPTC CF89-12 strains, the repeats were identified as being almost identical. Although a total of all 18 C. jejuni isolates examined gave PCR-positive signals for the CRISPRs genes, it was, interestingly, suggested that many numbers of C. lari and C. jejuni isolates may possibly carry cas but not CRISPRs genes within their CRISPRs loci. In addition, PCR amplification by using a novel primer pair of f-ClCRISPR-ladder and ClCRISPRs-R, which were novel to this study, with the UPTC CF89-12 strain was shown to be useful for the detection of the putative CRISPRs separated by the non-repetitive unique spacer regions, with the electrophoretic ladder DNA profile following 5.0 % polyacrylamide gel electrophoresis. Secondary structure models of the CRISPRs repeats were predicted with UPTC CF89-12 and two C. jejuni strains.     

The Small,Slow and Specialized CRISPR and Anti-CRISPR of Escherichia and Salmonella

Prokaryotes thrive in spite of the vast number and diversity of their viruses. This partly results from the evolution of mechanisms to inactivate or silence the action of exogenous DNA. Among these, Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) are unique in providing adaptive immunity against elements with high local resemblance to genomes of previously infecting agents. Here, we analyze the CRISPR loci of 51 complete genomes of Escherichia and Salmonella. CRISPR are in two pairs of loci in Escherichia, one single pair in Salmonella, each pair showing a similar turnover rate, repeat sequence and putative linkage to a common set of cas genes. Yet, phylogeny shows that CRISPR and associated cas genes have different evolutionary histories, the latter being frequently exchanged or lost. In our set, one CRISPR pair seems specialized in plasmids often matching genes coding for the replication, conjugation and antirestriction machinery. Strikingly, this pair also matches the cognate cas genes in which case these genes are absent. The unexpectedly high conservation of this anti-CRISPR suggests selection to counteract the invasion of mobile elements containing functional CRISPR/cas systems. There are few spacers in most CRISPR, which rarely match genomes of known phages. Furthermore, we found that strains divergent less than 250 thousand years ago show virtually identical CRISPR. The lack of congruence between cas, CRISPR and the species phylogeny and the slow pace of CRISPR change make CRISPR poor epidemiological markers in enterobacteria. All these observations are at odds with the expectedly abundant and dynamic repertoire of spacers in an immune system aiming at protecting bacteria from phages. Since we observe purifying selection for the maintenance of CRISPR these results suggest that alternative evolutionary roles for CRISPR remain to be uncovered.     

Thermoadaptation trait revealed by the genome sequence of thermophilic Geobacillus kaustophilus

We present herein the first complete genome sequence of a thermophilic Bacillus-related species, Geobacillus kaustophilus HTA426, which is composed of a 3.54 Mb chromosome and a 47.9 kb plasmid, along with a comparative analysis with five other mesophilic bacillar genomes. Upon orthologous grouping of the six bacillar sequenced genomes, it was found that 1257 common orthologous groups composed of 1308 genes (37%) are shared by all the bacilli, whereas 839 genes (24%) in the G.kaustophilus genome were found to be unique to that species. We were able to find the first prokaryotic sperm protamine P1 homolog, polyamine synthase, polyamine ABC transporter and RNA methylase in the 839 unique genes; these may contribute to thermophily by stabilizing the nucleic acids. Contrasting results were obtained from the principal component analysis (PCA) of the amino acid composition and synonymous codon usage for highlighting the thermophilic signature of the G.kaustophilus genome. Only in the PCA of the amino acid composition were the Bacillus-related species located near, but were distinguishable from, the borderline distinguishing thermophiles from mesophiles on the second principal axis. Further analysis revealed some asymmetric amino acid substitutions between the thermophiles and the mesophiles, which are possibly associated with the thermoadaptation of the organism.     

Comparative Analysis of the Orphan CRISPR2 Locus in 242 Enterococcus faecalis Strains

Clustered, Regularly Interspaced Short Palindromic Repeats and their associated Cas proteins (CRISPR-Cas) provide prokaryotes with a mechanism for defense against mobile genetic elements (MGEs). A CRISPR locus is a molecular memory of MGE encounters. It contains an array of short sequences, called spacers, that generally have sequence identity to MGEs. Three different CRISPR loci have been identified among strains of the opportunistic pathogen Enterococcus faecalis. CRISPR1 and CRISPR3 are associated with the cas genes necessary for blocking MGEs, but these loci are present in only a subset of E. faecalis strains. The orphan CRISPR2 lacks cas genes and is ubiquitous in E. faecalis, although its spacer content varies from strain to strain. Because CRISPR2 is a variable locus occurring in all E. faecalis, comparative analysis of CRISPR2 sequences may provide information about the clonality of E. faecalis strains. We examined CRISPR2 sequences from 228 E. faecalis genomes in relationship to subspecies phylogenetic lineages (sequence types; STs) determined by multilocus sequence typing (MLST), and to a genome phylogeny generated for a representative 71 genomes. We found that specific CRISPR2 sequences are associated with specific STs and with specific branches on the genome tree. To explore possible applications of CRISPR2 analysis, we evaluated 14 E. faecalis bloodstream isolates using CRISPR2 analysis and MLST. CRISPR2 analysis identified two groups of clonal strains among the 14 isolates, an assessment that was confirmed by MLST. CRISPR2 analysis was also used to accurately predict the ST of a subset of isolates. We conclude that CRISPR2 analysis, while not a replacement for MLST, is an inexpensive method to assess clonality among E. faecalis isolates, and can be used in conjunction with MLST to identify recombination events occurring between STs.     

Germ Warfare in a Microbial Mat Community: CRISPRs Provide Insights into the Co-Evolution of Host and Viral Genomes

CRISPR arrays and associated cas genes are widespread in bacteria and archaea and confer acquired resistance to viruses. To examine viral immunity in the context of naturally evolving microbial populations we analyzed genomic data from two thermophilic Synechococcus isolates (Syn OS-A and Syn OS-B′) as well as a prokaryotic metagenome and viral metagenome derived from microbial mats in hotsprings at Yellowstone National Park. Two distinct CRISPR types, distinguished by the repeat sequence, are found in both the Syn OS-A and Syn OS-B′ genomes. The genome of Syn OS-A contains a third CRISPR type with a distinct repeat sequence, which is not found in Syn OS-B′, but appears to be shared with other microorganisms that inhabit the mat. The CRISPR repeats identified in the microbial metagenome are highly conserved, while the spacer sequences (hereafter referred to as “viritopes” to emphasize their critical role in viral immunity) were mostly unique and had no high identity matches when searched against GenBank. Searching the viritopes against the viral metagenome, however, yielded several matches with high similarity some of which were within a gene identified as a likely viral lysozyme/lysin protein. Analysis of viral metagenome sequences corresponding to this lysozyme/lysin protein revealed several mutations all of which translate into silent or conservative mutations which are unlikely to affect protein function, but may help the virus evade the host CRISPR resistance mechanism. These results demonstrate the varied challenges presented by a natural virus population, and support the notion that the CRISPR/viritope system must be able to adapt quickly to provide host immunity. The ability of metagenomics to track population-level variation in viritope sequences allows for a culture-independent method for evaluating the fast co-evolution of host and viral genomes and its consequence on the structuring of complex microbial communities.     

Functional innovations of three chronological mesohexaploid Brassica rapa genomes

Background

The Brassicaceae family is an exemplary model for studying plant polyploidy. The Brassicaceae knowledge-base includes the well-annotated Arabidopsis thaliana reference sequence; well-established evidence for three rounds of whole genome duplication (WGD); and the conservation of genomic structure, with 24 conserved genomic blocks (GBs). The recently released Brassica rapa draft genome provides an ideal opportunity to update our knowledge of the conserved genomic structures in Brassica, and to study evolutionary innovations of the mesohexaploid plant, B. rapa.

Results

Three chronological B. rapa genomes (recent, young, and old) were reconstructed with sequence divergences, revealing a trace of recursive WGD events. A total of 636 fast evolving genes were unevenly distributed throughout the recent and young genomes. The representative Gene Ontology (GO) terms for these genes were ‘stress response’ and ‘development’ both through a change in protein modification or signaling, rather than by enhancing signal recognition. In retention patterns analysis, 98% of B. rapa genes were retained as collinear gene pairs; 77% of those were singly-retained in recent or young genomes resulting from death of the ancestral copies, while others were multi-retained as long retention genes. GO enrichments indicated that single retention genes mainly function in the interpretation of genetic information, whereas, multi-retention genes were biased toward signal response, especially regarding development and defense. In the recent genome, 13,302, 5,790, and 20 gene pairs were multi-retained following Brassica whole genome triplication (WGT) events with 2, 3, and 4 homoeologous copies, respectively. Enriched GO-slim terms from B. rapa homomoelogues imply that a major effect of the B. rapa WGT may have been to acquire environmental adaptability or to change the course of development. These homoeologues seem to more frequently undergo subfunctionalization with spatial expression patterns compared with other possible events including nonfunctionalization and neofunctionalization.

Conclusion

We refined Brassicaceae GB information using the latest genomic resources, and distinguished three chronologically ordered B. rapa genomes. B. rapa genes were categorized into fast evolving, single- and multi-retention genes, and long retention genes by their substitution rates and retention patterns. Representative functions of the categorized genes were elucidated, providing better understanding of B. rapa evolution and the Brassica genus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-606) contains supplementary material, which is available to authorized users.     

Genome Sequence of Lactobacillus pentosus KCA1: Vaginal Isolate from a Healthy Premenopausal Woman

The vaginal microbiota, in particular Lactobacillus species, play an important role in female health through modulation of immunity, countering pathogens and maintaining a pH below 4.7. We report the isolation and genome sequence of Lactobacillus pentosus strain KCA1 (formally known as L. plantarum) from the vagina of a healthy Nigerian woman. The genome was sequenced using Illumina GA II technology. The resulting 16,920,226 paired-end reads were assembled with the Velvet tool. Contigs were annotated using the RAST server, and manually curated. A comparative analysis with the available genomes of L. pentosus IG1 and L. plantarum WCFS1 showed that over 15% of the predicted functional activities are found only in this strain. The strain has a chromosome sequence of 3,418,159 bp with a G+C content of 46.4%, and is devoid of plasmids. Novel gene clusters or variants of known genes relative to the reference genomes were found. In particular, the strain has loci encoding additional putative mannose phosphotransferase systems. Clusters of genes include those for utilization of hydantoin, isopropylmalate, malonate, rhamnosides, and genes for assimilation of polyglycans, suggesting the metabolic versatility of L. pentosus KCA1. Loci encoding putative phage defense systems were also found including clustered regularly interspaced short palindromic repeats (CRISPRs), abortive infection (Abi) systems and toxin-antitoxin systems (TA). A putative cluster of genes for biosynthesis of a cyclic bacteriocin precursor, here designated as pentocin KCA1 (penA) were identified. These findings add crucial information for understanding the genomic and geographic diversity of vaginal lactobacilli.     

Bioinformatics analyses of <Emphasis Type="Italic">Shigella</Emphasis> CRISPR structure and spacer classification

Clustered regularly interspaced short palindromic repeats (CRISPR) are inheritable genetic elements of a variety of archaea and bacteria and indicative of the bacterial ecological adaptation, conferring acquired immunity against invading foreign nucleic acids. Shigella is an important pathogen for anthroponosis. This study aimed to analyze the features of Shigella CRISPR structure and classify the spacers through bioinformatics approach. Among 107 Shigella, 434 CRISPR structure loci were identified with two to seven loci in different strains. CRISPR-Q1, CRISPR-Q4 and CRISPR-Q5 were widely distributed in Shigella strains. Comparison of the first and last repeats of CRISPR1, CRISPR2 and CRISPR3 revealed several base variants and different stem-loop structures. A total of 259 cas genes were found among these 107 Shigella strains. The cas gene deletions were discovered in 88 strains. However, there is one strain that does not contain cas gene. Intact clusters of cas genes were found in 19 strains. From comprehensive analysis of sequence signature and BLAST and CRISPRTarget score, the 708 spacers were classified into three subtypes: Type I, Type II and Type III. Of them, Type I spacer referred to those linked with one gene segment, Type II spacer linked with two or more different gene segments, and Type III spacer undefined. This study examined the diversity of CRISPR/cas system in Shigella strains, demonstrated the main features of CRISPR structure and spacer classification, which provided critical information for elucidation of the mechanisms of spacer formation and exploration of the role the spacers play in the function of the CRISPR/cas system.     

Strain-level genomic variation in natural populations of Lebetimonas from an erupting deep-sea volcano

Chemolithoautotrophic Epsilonproteobacteria are ubiquitous in sulfidic, oxygen-poor habitats, including hydrothermal vents, marine oxygen minimum zones, marine sediments and sulfidic caves and have a significant role in cycling carbon, hydrogen, nitrogen and sulfur in these environments. The isolation of diverse strains of Epsilonproteobacteria and the sequencing of their genomes have revealed that this group has the metabolic potential to occupy a wide range of niches, particularly at dynamic deep-sea hydrothermal vents. We expand on this body of work by examining the population genomics of six strains of Lebetimonas, a vent-endemic, thermophilic, hydrogen-oxidizing Epsilonproteobacterium, from a single seamount in the Mariana Arc. Using Lebetimonas as a model for anaerobic, moderately thermophilic organisms in the warm, anoxic subseafloor environment, we show that genomic content is highly conserved and that recombination is limited between closely related strains. The Lebetimonas genomes are shaped by mobile genetic elements and gene loss as well as the acquisition of novel functional genes by horizontal gene transfer, which provide the potential for adaptation and microbial speciation in the deep sea. In addition, these Lebetimonas genomes contain two operons of nitrogenase genes with different evolutionary origins. Lebetimonas expressed nifH during growth with nitrogen gas as the sole nitrogen source, thus providing the first evidence of nitrogen fixation in any Epsilonproteobacteria from deep-sea hydrothermal vents. In this study, we provide a comparative overview of the genomic potential within the Nautiliaceae as well as among more distantly related hydrothermal vent Epsilonproteobacteria to broaden our understanding of microbial adaptation and diversity in the deep sea.     

Characterization of the endonuclease SSO2001 from Sulfolobus solfataricus P2

Clustered regularly interspaced short palindromic repeats (CRISPR) and their associated protein genes (cas genes) have been suggested to act as an immune system in archaea and bacteria mimicking the eukaryotic RNA interference (RNAi) system. We have investigated the properties of the protein SSO2001 from Sulfolobus solfataricus (Sso) P2, which is part of the cas gene cluster. This study shows that SSO2001 is an endonuclease specifically digesting double-stranded oligonucleotides and preferably cleaving at G:C pairs. Point mutations identify both highly conserved aspartate and glutamate residues as being crucial for the nuclease activity. The catalytic activity shows an optimum at neutral pH and pH 3.     

Diversity, activity, and evolution of CRISPR loci in Streptococcus thermophilus

Clustered regularly interspaced short palindromic repeats (CRISPR) are hypervariable loci widely distributed in prokaryotes that provide acquired immunity against foreign genetic elements. Here, we characterize a novel Streptococcus thermophilus locus, CRISPR3, and experimentally demonstrate its ability to integrate novel spacers in response to bacteriophage. Also, we analyze CRISPR diversity and activity across three distinct CRISPR loci in several S. thermophilus strains. We show that both CRISPR repeats and cas genes are locus specific and functionally coupled. A total of 124 strains were studied, and 109 unique spacer arrangements were observed across the three CRISPR loci. Overall, 3,626 spacers were analyzed, including 2,829 for CRISPR1 (782 unique), 173 for CRISPR2 (16 unique), and 624 for CRISPR3 (154 unique). Sequence analysis of the spacers revealed homology and identity to phage sequences (77%), plasmid sequences (16%), and S. thermophilus chromosomal sequences (7%). Polymorphisms were observed for the CRISPR repeats, CRISPR spacers, cas genes, CRISPR motif, locus architecture, and specific sequence content. Interestingly, CRISPR loci evolved both via polarized addition of novel spacers after exposure to foreign genetic elements and via internal deletion of spacers. We hypothesize that the level of diversity is correlated with relative CRISPR activity and propose that the activity is highest for CRISPR1, followed by CRISPR3, while CRISPR2 may be degenerate. Globally, the dynamic nature of CRISPR loci might prove valuable for typing and comparative analyses of strains and microbial populations. Also, CRISPRs provide critical insights into the relationships between prokaryotes and their environments, notably the coevolution of host and viral genomes.     

Genomic repeats, genome plasticity and the dynamics of Mycoplasma evolution

Mycoplasmas evolved by a drastic reduction in genome size, but their genomes contain numerous repeated sequences with important roles in their evolution. We have established a bioinformatic strategy to detect the major recombination hot-spots in the genomes of Mycoplasma pneumoniae, Mycoplasma genitalium, Ureaplasma urealyticum and Mycoplasma pulmonis. This allowed the identification of large numbers of potentially variable regions, as well as a comparison of the relative recombination potentials of different genomic regions. Different trends are perceptible among mycoplasmas, probably due to different functional and structural constraints. The largest potential for illegitimate recombination in M.pulmonis is found at the vsa locus and its comparison in two different strains reveals numerous changes since divergence. On the other hand, the main M.pneumoniae and M.genitalium adhesins rely on large distant repeats and, hence, homologous recombination for variation. However, the relation between the existence of repeats and antigenic variation is not necessarily straightforward, since repeats of P1 adhesin were found to be anti-correlated with epitopes recognized by patient antibodies. These different strategies have important consequences for the structures of genomes, since large distant repeats correlate well with the major chromosomal rearrangements. Probably to avoid such events, mycoplasmas strongly avoid inverse repeats, in comparison to co-oriented repeats.     

Core and accessory genome architecture in a group of Pseudomonas aeruginosa Mu-like phages

Background

Bacteriophages that infect the opportunistic pathogen Pseudomonas aeruginosa have been classified into several groups. One of them, which includes temperate phage particles with icosahedral heads and long flexible tails, bears genomes whose architecture and replication mechanism, but not their nucleotide sequences, are like those of coliphage Mu. By comparing the genomic sequences of this group of P. aeruginosa phages one could draw conclusions about their ontogeny and evolution.

Results

Two newly isolated Mu-like phages of P. aeruginosa are described and their genomes sequenced and compared with those available in the public data banks. The genome sequences of the two phages are similar to each other and to those of a group of P. aeruginosa transposable phages. Comparing twelve of these genomes revealed a common genomic architecture in the group. Each phage genome had numerous genes with homologues in all the other genomes and a set of variable genes specific for each genome. The first group, which comprised most of the genes with assigned functions, was named “core genome”, and the second group, containing mostly short ORFs without assigned functions was called “accessory genome”. Like in other phage groups, variable genes are confined to specific regions in the genome.

Conclusion

Based on the known and inferred functions for some of the variable genes of the phages analyzed here, they appear to confer selective advantages for the phage survival under particular host conditions. We speculate that phages have developed a mechanism for horizontally acquiring genes to incorporate them at specific loci in the genome that help phage adaptation to the selective pressures imposed by the host.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1146) contains supplementary material, which is available to authorized users.     

Ectopic Repo suppresses expression of castor gene in Drosophila central nervous system

The ventral nerve cord (VNC) of the Drosophila embryo is derived from neuroblasts (NBs). NBs divide in a stem cell lineage to generate a series of ganglion mother cells (GMCs), each of which divides once to produce a pair of neurons or glial cells. One of the NB genes, castor (cas), is expressed in a subset of NBs and has never been identified in neurons and the peripheral nervous system; cas plays a role in axonogenesis. But its limited expression along the dorsal-ventral axis within the central nervous system has not been investigated yet. In the present study, we examined the expression patterns of both genes using confocal microscopy to determine the effects of repo mutation on cas expression. Cas was mainly expressed in layers different from repo-expressed layers during early embryogenesis: repo was expressed mostly from deep to mid layers, while cas, from mid to superficial layers. Loss-of-function of repo did not result in an ectopic expression of cas, but rather, a scattering of cas-expressing cells. However, repo gain-of-function mutation caused repression of cas. In addition, repo-expressing cells seemed to block the migration of cas-expressing cells.     

Degradation of biological macromolecules supports uncultured microbial populations in Guaymas Basin hydrothermal sediments

Hydrothermal sediments contain large numbers of uncultured heterotrophic microbial lineages. Here, we amended Guaymas Basin sediments with proteins, polysaccharides, nucleic acids or lipids under different redox conditions and cultivated heterotrophic thermophiles with the genomic potential for macromolecule degradation. We reconstructed 20 metagenome-assembled genomes (MAGs) of uncultured lineages affiliating with known archaeal and bacterial phyla, including endospore-forming Bacilli and candidate phylum Marinisomatota. One Marinisomatota MAG had 35 different glycoside hydrolases often in multiple copies, seven extracellular CAZymes, six polysaccharide lyases, and multiple sugar transporters. This population has the potential to degrade a broad spectrum of polysaccharides including chitin, cellulose, pectin, alginate, chondroitin, and carrageenan. We also describe thermophiles affiliating with the genera Thermosyntropha, Thermovirga, and Kosmotoga with the capability to make a living on nucleic acids, lipids, or multiple macromolecule classes, respectively. Several populations seemed to lack extracellular enzyme machinery and thus likely scavenged oligo- or monomers (e.g., MAGs affiliating with Archaeoglobus) or metabolic products like hydrogen (e.g., MAGs affiliating with Thermodesulfobacterium or Desulforudaceae). The growth of methanogens or the production of methane was not observed in any condition, indicating that the tested macromolecules are not degraded into substrates for methanogenesis in hydrothermal sediments. We provide new insights into the niches, and genomes of microorganisms that actively degrade abundant necromass macromolecules under oxic, sulfate-reducing, and fermentative thermophilic conditions. These findings improve our understanding of the carbon flow across trophic levels and indicate how primary produced biomass sustains complex and productive ecosystems.Subject terms: Water microbiology, Environmental sciences     

Comparative genomics of two newly isolated Dehalococcoides strains and an enrichment using a genus microarray

Comparative genomics of Dehalococcoides strains and an enrichment were performed using a microarray targeting genes from all available sequenced genomes of the Dehalococcoides genus. The microarray was designed with 4305 probe sets to target 98.6% of the open-reading frames from strains 195, CBDB1, BAV1 and VS. The microarrays were validated and applied to query the genomes of two recently isolated Dehalococcoides strains, ANAS1 and ANAS2, and their enrichment source (ANAS) to understand the genome–physiology relationships. Strains ANAS1 and ANAS2 can both couple the reduction of trichloroethene, cis-dichloroethene (DCE) and 1,1-DCE, but not tetrachloroethene and trans-DCE with growth, whereas only strain ANAS2 couples vinyl chloride reduction to growth. Comparative genomic analysis showed that the genomes of both strains are similar to each other and to strain 195, except for genes that are within the previously defined integrated elements or high-plasticity regions. Combined results of the two isolates closely matched the results obtained using genomic DNA of the ANAS enrichment. The genome similarities, together with the distinct chlorinated ethene usage of strains ANAS1, ANAS2 and 195 demonstrate that closely phylogenetically related strains can be physiologically different. This incongruence between physiology and core genome phylogeny seems to be related to the presence of distinct reductive dehalogenase-encoding genes with assigned chlorinated ethene functions (pceA, tceA in strain 195; tceA in strain ANAS1; vcrA in strain ANAS2). Overall, the microarrays are a valuable high-throughput tool for comparative genomics of unsequenced Dehalococcoides-containing samples to provide insights into their gene content and dechlorination functions.