The genome of Polaromonas naphthalenivorans strain CJ2, isolated from coal tar-contaminated sediment, reveals physiological and metabolic versatility and evolution through extensive horizontal gene transfer

Arctic bacterial diversity and connectivity in the coastal margin of the Last Ice Area

Arctic climate change is leading to sea-ice attrition in the Last Ice Area along the northern coast of Canada and Greenland, but less attention has been given to the associated land-based ecosystems. Here we evaluated bacterial community structure in a hydrologically coupled cryo-ecosystem in the region: Thores Glacier, proglacial Thores Lake, and its outlet to the sea. Deep amplicon sequencing revealed that Polaromonas was ubiquitous, but differed genetically among diverse niches. Surface glacier-ice was dominated by Cyanobacteria, while the perennially ice-capped, well-mixed water column of Thores Lake had a unique assemblage of Chloroflexi, Actinobacteriota, and Planctomycetota. Species richness increased downstream, but glacier microbes were little detected in the lake, suggesting strong taxonomic sorting. Ongoing climate change and the retreat of Thores Glacier would lead to complete drainage and loss of the lake microbial ecosystem, indicating the extreme vulnerability of diverse cryohabitats and unique microbiomes in the Last Ice coastal margin.

A global perspective on bacterial diversity in the terrestrial deep subsurface

While recent efforts to catalogue Earth's microbial diversity have focused upon surface and marine habitats, 12-20 % of Earth's biomass is suggested to exist in the terrestrial deep subsurface, compared to ~1.8 % in the deep subseafloor. Metagenomic studies of the terrestrial deep subsurface have yielded a trove of divergent and functionally important microbiomes from a range of localities. However, a wider perspective of microbial diversity and its relationship to environmental conditions within the terrestrial deep subsurface is still required. Our meta-analysis reveals that terrestrial deep subsurface microbiota are dominated by Betaproteobacteria, Gammaproteobacteria and Firmicutes, probably as a function of the diverse metabolic strategies of these taxa. Evidence was also found for a common small consortium of prevalent Betaproteobacteria and Gammaproteobacteria operational taxonomic units across the localities. This implies a core terrestrial deep subsurface community, irrespective of aquifer lithology, depth and other variables, that may play an important role in colonizing and sustaining microbial habitats in the deep terrestrial subsurface. An in silico contamination-aware approach to analysing this dataset underscores the importance of downstream methods for assuring that robust conclusions can be reached from deep subsurface-derived sequencing data. Understanding the global panorama of microbial diversity and ecological dynamics in the deep terrestrial subsurface provides a first step towards understanding the role of microbes in global subsurface element and nutrient cycling.

Micro-synteny conservation analysis revealed the evolutionary history of bacterial biphenyl degradation pathway

Phenolic compounds have been enlisted by the United States Environmental Protection Agency (USEPA) and the European Union (EU) as pollutants of priority concern. The biphenyl degradation pathway plays an essential role in prokaryote polychlorinated biphenyls degradation. Our understanding of prokaryotic pathways and their evolution has dramatically increased in recent years with the advancements in prokaryotic genome sequencing and analysis tools. In this work, we applied bioinformatics tools to study the evolution of the biphenyl degradation pathway focusing on the phylogeny and initiation of four representative species (Burkholderia xenovorans LB400, Polaromonas naphthalenivorans CJ2, Pseudomonas putida F1 and Rhodococcus jostii RHA1). These species contained partial or full concatenated genes from bph gene cluster (i.e. bphRbphA1A2A3A4BCKHJID). The aim was to establish this pathway's origin and development mode in the prokaryotic world. Genomic screening revealed that many bacterial species possess genes for the biphenyl degradation pathway. However, the micro-synteny conservation analysis indicated that massive gene recruitment events might have occurred during the evolution of the biphenyl degradation pathway. Combining with the phylogenetic positions, this work points to the evolutionary process of acquiring the biphenyl degradation pathway by different fragments through horizontal gene transfer in these bacterial groups. This study reports the first-ever evidence of the birth of this pathway in the represented species.

Draft Genome Sequence of Polaromonas eurypsychrophila AER18D-145, Isolated from a Uranium Tailings Management Facility in Northern Saskatchewan, Canada

The 4.8-Mbp draft genome sequence of Polaromonas eurypsychrophila AER18D-145, isolated from a uranium tailings management facility, is reported. The sequence may provide insights into the mechanisms of the hypertolerance of this strain to extreme conditions and help determine its potential for bioremediation applications.

Functional filtering and random processes affect the assembly of microbial communities of snow algae blooms at Maritime Antarctic

The increasing temperatures at the West Antarctic Peninsula (Maritime Antarctic) could lead to a higher occurrence of snow algal blooms which are ubiquitous events that change the snow coloration, reducing albedo and in turn exacerbating melting. However, there is a limited understanding of snow algae blooms biodiversity, composition, and their functional profiles, especially in one of the world's areas most affected by climate change. In this study we used 16S rRNA and 18S rRNA metabarcoding, and shotgun metagenomics to assess the diversity, composition, and functional potential of the snow algae blooms bacterial and eukaryotic communities at three different sites of Maritime Antarctic, between different colors of the algae blooms and between seasonal and semi-permanent snowfields. We tested the hypothesis that the functional potential of snow algae blooms is conserved despite a changing taxonomic composition. Furthermore, we determined taxonomic co-occurrence patterns of bacteria and eukaryotes and assessed the potential for the exchange of metabolites among bacterial taxa. Here, we tested the prediction that there are co-occurring taxa within snow algae whose biotic interactions are marked by the exchange of metabolites. Our results show that the composition of snow algae blooms vary significantly among sites. For instance, a higher abundance of fungi and protists were detected in Fildes Peninsula compared with Doumer Island and O'Higgins. Likewise, the composition varied between snow colors and snow types. However, the functional potential varied only among sampling sites with a higher abundance of genes involved in tolerance to environmental stress at O'Higgins. Co-occurrence patterns of dominant bacterial genera such as Pedobacter, Polaromonas, Flavobacterium and Hymenobacter were recorded, contrasting the absence of co-occurring patterns displayed by Chlamydomonadales algae with other eukaryotes. Finally, genome-scale metabolic models revealed that bacteria within snow algae blooms likely compete for resources instead of forming cooperative communities.

Metagenomic analysis of microbial communities across a transect from low to highly hydrocarbon-contaminated soils in King George Island, Maritime Antarctica

Soil samples from a transect from low to highly hydrocarbon-contaminated soils were collected around the Brazilian Antarctic Station Comandante Ferraz (EACF), located at King George Island, Antarctica. Quantitative PCR (qPCR) analysis of bacterial 16S rRNA genes, 16S rRNA gene (iTag), and shotgun metagenomic sequencing were used to characterize microbial community structure and the potential for petroleum degradation by indigenous microbes. Hydrocarbon contamination did not affect bacterial abundance in EACF soils (bacterial 16S rRNA gene qPCR). However, analysis of 16S rRNA gene sequences revealed a successive change in the microbial community along the pollution gradient. Microbial richness and diversity decreased with the increase of hydrocarbon concentration in EACF soils. The abundance of Cytophaga, Methyloversatilis, Polaromonas, and Williamsia was positively correlated (p-value = <.05) with the concentration of total petroleum hydrocarbons (TPH) and/or polycyclic aromatic hydrocarbons (PAH). Annotation of metagenomic data revealed that the most abundant hydrocarbon degradation pathway in EACF soils was related to alkyl derivative-PAH degradation (mainly methylnaphthalenes) via the CYP450 enzyme family. The abundance of genes related to nitrogen fixation increased in EACF soils as the concentration of hydrocarbons increased. The results obtained here are valuable for the future of bioremediation of petroleum hydrocarbon-contaminated soils in polar environments.

Influence of prokaryotic microorganisms on initial soil formation along a glacier forefield on King George Island, maritime Antarctica

Compared to the 1970s, the edge of the Ecology Glacier on King George Island, maritime Antarctica, is positioned more than 500 m inwards, exposing a large area of new terrain to soil-forming processes and periglacial climate for more than 40 years. To gain information on the state of soil formation and its interplay with microbial activity, three hyperskeletic Cryosols (vegetation cover of 0-80%) deglaciated after 1979 in the foreland of the Ecology Glacier and a Cambic Cryosol (vegetation cover of 100%) distal to the lateral moraine deglaciated before 1956 were investigated by combining soil chemical and microbiological methods. In the upper part of all soils, a decrease in soil pH was observed, but only the Cambic Cryosol showed a clear direction of pedogenic and weathering processes, such as initial silicate weathering indicated by a decreasing Chemical Index of Alteration with depth. Differences in the development of these initial soils could be related to different microbial community compositions and vegetation coverage, despite the short distance among them. We observed-decreasing with depth-the highest bacterial abundances and microbial diversity at vegetated sites. Multiple clusters of abundant amplicon sequence variants were found depending on the site-specific characteristics as well as a distinct shift in the microbial community structure towards more similar communities at soil depths > 10 cm. In the foreland of the Ecology Glacier, the main soil-forming processes on a decadal timescale are acidification and accumulation of soil organic carbon and nitrogen, accompanied by changes in microbial abundances, microbial community compositions, and plant coverage, whereas quantifiable silicate weathering and the formation of pedogenic oxides occur on a centennial to a millennial timescale after deglaciation.

Alkylphenols and Chlorophenols Remediation in Vertical Flow Constructed Wetlands: Removal Efficiency and Microbial Community Response

This study aims to investigate the effect of two different groups of phenolic compounds (the alkylphenols nonylphenol (NP) and octylphenol (OP), and the chlorophenol pentachlorophenol (PCP)) on constructed wetlands (CWs) performance, including on organic matter, nutrients and contaminants removal efficiency, and on microbial community structure in the plant bed substrate. CWs were assembled at lab scale simulating a vertical flow configuration and irrigated along eight weeks with Ribeira de Joane (an urban stream) water not doped (control) or doped with a mixture of NP and OP or with PCP (at a 100 μg·L−1 concentration each). The presence of the phenolic contaminants did not interfere in the removal of organic matter or nutrients in CWs in the long term. Removals of NP and OP were >99%, whereas PCP removals varied between 87% and 98%, mainly due to biodegradation. Microbial richness, diversity and dominance in CWs substrate were generally not affected by phenolic compounds, with only PCP decreasing diversity. Microbial community structure, however, showed that there was an adaptation of the microbial community to the presence of each contaminant, with several specialist genera being enriched following exposure. The three more abundant specialist genera were Methylotenera and Methylophilus (methylophilaceae family) and Hyphomicrobium (hyphomicrobiaceae family) when the systems were exposed to a mixture of NP and OP. When exposed to PCP, the three more abundant genera were Denitromonas (Rhodocyclaceae family), Xenococcus_PCC_7305 (Xenococcaceae family) and Rhodocyclaceae_uncultured (Rhodocyclaceae family). To increase CWs efficiency in the elimination of phenolic compounds, namely PCP which was not totally removed, strategies to stimulate (namely biostimulation) or increase (namely bioaugmentation) the presence of these bacteria should be explore. This study clearly shows the potential of vertical flow CWs for the removal of phenolic compounds, a still little explored subject, contributing to promote the use of CWs as nature-based solutions to remediate water contaminated with different families of persistent and/or emergent contaminants.

VV Reduction by Polaromonas spp. in Vanadium Mine Tailings

Vanadium (V) is an important metal with critical industrial and medical applications. Elevated V contamination, however, can be a threat to the environment and human health. Microorganisms can reduce the more toxic and mobile V^V to the less toxic and immobile V^IV, which could be a detoxification and energy metabolism strategy adopted by V-reducing bacteria (VRB). The limited understanding of microbial responses to V contamination and the mechanisms for V^V reduction, however, hamper our capability to attenuate V contamination. This study focused on determining the microbial responses to elevated V concentration and the mechanisms of V^V reduction in V tailings. The bacterial communities were characterized and compared between the V tailings and the less contaminated adjacent mineral soils. Further, V^V-reducing enrichments indicated that bacteria associated with Polaromonas, a genus belonging to the family Burkholderiaceae, were potentially responsible for V^V reduction. Retrieved metagenome-assembled genomes (MAGs) suggested that the Polaromonas spp. encoded genes (cymA, omcA, and narG) were responsible for V^V reduction. Additionally, Polaromonas spp. was metabolically versatile and could use both organic and inorganic electron donors. The metabolic versatility of Polaromonas spp. may be important for its ability to flourish in the V tailings.

Molecular Characterization and Comparative Genomics of IncQ-3 Plasmids Conferring Resistance to Various Antibiotics Isolated from a Wastewater Treatment Plant in Warsaw (Poland)

As small, mobilizable replicons with a broad host range, IncQ plasmids are widely distributed among clinical and environmental bacteria. They carry antibiotic resistance genes, and it has been shown that they confer resistance to β-lactams, fluoroquinolones, aminoglycosides, trimethoprim, sulphonamides, and tetracycline. The previously proposed classification system divides the plasmid group into four subgroups, i.e., IncQ-1, IncQ-2, IncQ-3, and IncQ-4. The last two subgroups have been poorly described so far. The aim of this study was to analyze five newly identified IncQ-3 plasmids isolated from a wastewater treatment plant in Poland and to compare them with all known plasmids belonging to the IncQ-3 subgroup whose sequences were retrieved from the NCBI database. The complete nucleotide sequences of the novel plasmids were annotated and bioinformatic analyses were performed, including identification of core genes and auxiliary genetic load. Furthermore, functional experiments testing plasmid mobility were carried out. Phylogenetic analysis based on three core genes (repA, mobA/repB, and mobC) revealed the presence of three main clusters of IncQ-3 replicons. Apart from having a highly conserved core, the analyzed IncQ-3 plasmids were vectors of antibiotic resistance genes, including (I) the qnrS2 gene that encodes fluoroquinolone resistance and (II) β-lactam, trimethoprim, and aminoglycoside resistance genes within integron cassettes.

Complete genome of Polaromonas vacuolata KCTC 22033T isolated from beneath Antarctic Sea ice

Polaromonas vacuolata KCTC 22033^T is an obligate aerobic, Gram-negative, psychrophilic and rod-shaped bacterium isolated from beneath the sea ice off the coast of the Palmer Peninsula, Anvers Islands, Antarctica. P. vacuolata is the type species of Polaromonas genus and the first example of gas vacuolate Betaproteobacteria isolated from marine habitats. Here, we report a complete genome of P. vacuolata KCTC 22033^T, which consists of 3,837,686 bp (G + C content of 52.07%) with a single chromosome, 3461 protein-coding genes, 56 tRNAs and 6 rRNA operons. Genomic analysis revealed the presence of genes involved in bacterial adaptation under saline conditions, cold adaptation via the production of gas vesicles and cell adhesion proteins, and a photoheterotrophic lifestyle when challenged by starvation. Intriguingly, several of these genes were likely acquired from species outside the Polaromonas genus. The genomic information therefore describes the unique evolution and adaptation of P. vacuolata to its extraordinary habitat, i.e., beneath the Antarctic sea ice.

Functional Analysis of a Polluted River Microbiome Reveals a Metabolic Potential for Bioremediation

The objective of this study is to understand the functional and metabolic potential of the microbial communities along the Apatlaco River and highlight activities related to bioremediation and its relationship with the Apatlaco’s pollutants, to enhance future design of more accurate bioremediation processes. Water samples were collected at four sampling sites along the Apatlaco River (S1–S4) and a whole metagenome shotgun sequencing was performed to survey and understand the microbial metabolic functions with potential for bioremediation. A HMMER search was used to detect sequence homologs related to polyethylene terephthalate (PET) and polystyrene biodegradation, along with bacterial metal tolerance in Apatlaco River metagenomes. Our results suggest that pollution is a selective pressure which enriches microorganisms at polluted sites, displaying metabolic capacities to tolerate and transform the contamination. According to KEGG annotation, all sites along the river have bacteria with genes related to xenobiotic biodegradation. In particular, functions such as environmental processing, xenobiotic biodegradation and glycan biosynthesis are over-represented in polluted samples, in comparison to those in the clean water site. This suggests a functional specialization in the communities that inhabit each perturbated point. Our results can contribute to the determination of the partition in a metabolic niche among different Apatlaco River prokaryotic communities, that help to contend with and understand the effect of anthropogenic contamination.

Novel copper-containing membrane monooxygenases (CuMMOs) encoded by alkane-utilizing Betaproteobacteria

Copper-containing membrane monooxygenases (CuMMOs) are encoded by xmoCAB(D) gene clusters and catalyze the oxidation of methane, ammonia, or some short-chain alkanes and alkenes. In a metagenome constructed from an oilsands tailings pond we detected an xmoCABD gene cluster with <59% derived protein sequence identity to genes from known bacteria. Stable isotope probing experiments combined with a specific xmoA qPCR assay demonstrated that the bacteria possessing these genes were incapable of methane assimilation, but did grow on ethane and propane. Single-cell amplified genomes (SAGs) from propane-enriched samples were screened with the specific PCR assay to identify bacteria possessing the target gene cluster. Multiple SAGs of Betaproteobacteria belonging to the genera Rhodoferax and Polaromonas possessed homologues of the metagenomic xmoCABD gene cluster. Unexpectedly, each of these two genera also possessed other xmoCABD paralogs, representing two additional lineages in phylogenetic analyses. Metabolic reconstructions from SAGs predicted that neither bacterium encoded enzymes with the potential to support catabolic methane or ammonia oxidation, but that both were capable of higher n-alkane degradation. The involvement of the encoded CuMMOs in alkane oxidation was further suggested by reverse transcription PCR analyses, which detected elevated transcription of the xmoA genes upon enrichment of water samples with propane as the sole energy source. Enrichments, isotope incorporation studies, genome reconstructions, and gene expression studies therefore all agreed that the unknown xmoCABD operons did not encode methane or ammonia monooxygenases, but rather n-alkane monooxygenases. This study broadens the known diversity of CuMMOs and identifies these enzymes in non-nitrifying Betaproteobacteria.

A PAH-degrading bacterial community enriched with contaminated agricultural soil and its utility for microbial bioremediation

A bacterial community was enriched with polycyclic aromatic hydrocarbons (PAHs) polluted soil to better study PAH degradation by indigenous soil bacteria. The consortium degraded more than 52% of low molecular weight and 35% of high molecular weight (HMW) PAHs during 16 days in a soil leachate medium. 16S rRNA gene high-throughput sequencing and quantitative polymerase chain reaction analyses for alpha subunit genes of ring-hydroxylating-dioxygenase (RHDα) suggested that Proteobacteria and Actinobacteria at the phylum level, Pseudomonas, Methylobacillus, Nocardioides, Methylophilaceae, Achromobacter, Pseudoxanthomonas, and Caulobacter at the generic level were involved in PAH degradation and might have the ability to carry RHDα genes (nidA and nahAc). The community was selected and collected according to biomass and RHDα gene contents, and added back to the PAH-polluted soil. The 16 EPA priority PAHs decreased from 95.23 to 23.41 mg kg^-1 over 35 days. Compared with soil without the introduction of this bacterial community, adding the community with RHDα genes significantly decreased soil PAH contents, particularly HMW PAHs. The metabolic rate of PAHs in soil was positively correlated with nidA and nahAc gene contents. These results indicate that adding an indigenous bacterial consortium containing RHDα genes to contaminated soil may be a feasible and environmentally friendly method to clean up PAHs in agricultural soil.

Plasmids carrying antimicrobial resistance genes in Enterobacteriaceae

Bacterial antimicrobial resistance (AMR) is constantly evolving and horizontal gene transfer through plasmids plays a major role. The identification of plasmid characteristics and their association with different bacterial hosts provides crucial knowledge that is essential to understand the contribution of plasmids to the transmission of AMR determinants. Molecular identification of plasmid and strain genotypes elicits a distinction between spread of AMR genes by plasmids and dissemination of these genes by spread of bacterial clones. For this reason several methods are used to type the plasmids, e.g. PCR-based replicon typing (PBRT) or relaxase typing. Currently, there are 28 known plasmid types in Enterobacteriaceae distinguished by PBRT. Frequently reported plasmids [IncF, IncI, IncA/C, IncL (previously designated IncL/M), IncN and IncH] are the ones that bear the greatest variety of resistance genes. The purpose of this review is to provide an overview of all known AMR-related plasmid families in Enterobacteriaceae, the resistance genes they carry and their geographical distribution.

Insight Into Metabolic Versatility of an Aromatic Compounds-Degrading Arthrobacter sp. YC-RL1

The genus Arthrobacter is ubiquitously distributed in different natural environments. Many xenobiotic-degrading Arthrobacter strains have been isolated and described; however, few have been systematically characterized with regard to multiple interrelated metabolic pathways and the genes that encode them. In this study, the biodegradability of seven aromatic compounds by Arthrobacter sp. YC-RL1 was investigated. Strain YC-RL1 could efficiently degrade p-xylene (PX), naphthalene, phenanthrene, biphenyl, p-nitrophenol (PNP), and bisphenol A (BPA) under both separated and mixed conditions. Based on the detected metabolic intermediates, metabolic pathways of naphthalene, biphenyl, PNP, and BPA were proposed, which indicated that strain YC-RL1 harbors systematic metabolic pathways toward aromatic compounds. Further, genomic analysis uncovered part of genes involved in the proposed pathways. Both intradiol and extradiol ring-cleavage dioxygenase genes were identified in the genome of strain YC-RL1. Meanwhile, gene clusters predicted to encode the degradation of biphenyl (bph), para-substituted phenols (npd) and protocatechuate (pca) were identified, and bphA1A2BCD was proposed to be a novel biphenyl-degrading gene cluster. The complete metabolic pathway of biphenyl was deduced via intermediates and functional gene analysis (bph and pca gene clusters). One of the these genes encoding ring-cleavage dioxygenase in bph gene cluster, a predicted 2,3-dihydroxybiphenyl 1,2-dioxygenase (BphC) gene, was cloned and its activity was confirmed by heterologous expression. This work systematically illuminated the metabolic versatility of aromatic compounds in strain YC-RL1 via the combination of metabolites identification, genomics analysis and laboratory experiments. These results suggested that strain YC-RL1 might be a promising candidate for the bioremediation of aromatic compounds pollution sites.

Microbial connectivity and sorting in a High Arctic watershed

Aquatic ecosystems in the High Arctic are facing unprecedented changes as a result of global warming effects on the cryosphere. Snow pack is a central feature of northern landscapes, but the snow microbiome and its microbial connectivity to adjacent and downstream habitats have been little explored. To evaluate these aspects, we sampled along a hydrologic continuum at Ward Hunt Lake (latitude 83°N) in the Canadian High Arctic, from snow banks, water tracks in the permafrost catchment, the upper and lower strata of the lake, and the lake outlet and its coastal marine mixing zone. The microbial communities were analyzed by high-throughput sequencing of 16 and 18S rRNA to determine the composition of potentially active Bacteria, Archaea and microbial Eukarya. Each habitat had distinct microbial assemblages, with highest species richness in the subsurface water tracks that connected the melting snow to the lake. However, up to 30% of phylotypes were shared along the hydrologic continuum, showing that many taxa originating from the snow can remain in the active fraction of downstream microbiomes. The results imply that changes in snowfall associated with climate warming will affect microbial community structure throughout all spatially connected habitats within snow-fed polar ecosystems.

Metataxonomy and functionality of wood-tar degrading microbial consortia

Wood-tar is a liquid material obtained by wood gasification process, and comprises several polycyclic aromatic hydrocarbons (PAH). Tar biodegradation is a very challenging task, due to its toxicity and to its complex chemistry. The 'microbial resource management' concerns the use of environmental microbial communities potentially able to provide us services. We applied this concept in tar biodegradation. Tar composed by several PAH (including phenanthrene, acenaphthylene and fluorene) was subjected to a biodegradation process in triplicate microcosms spiked with a microbial community collected from PAH-rich soils. In 20 days, 98.9% of tar was mineralized or adsorbed to floccules, while negative controls showed poor PAH reduction. The dynamics of fungal and bacterial communities was assessed through Automated Ribosomal Intergenic Spacer Analysis (ARISA), 454 pyrosequencing of the fungal ITS and of the bacterial 16S rRNA. Quantification of the degrading bacterial communities was performed via quantitative Real Time PCR of the 16S rRNA genes and of the cathecol 2,3-dioxygenase genes. Results showed the importance of fungal tar-degrading populations in the first period of incubation, followed by a complex bacterial dynamical growth ruled by co-feeding behaviors.

Plasmids of Psychrotolerant Polaromonas spp. Isolated From Arctic and Antarctic Glaciers - Diversity and Role in Adaptation to Polar Environments

Cold-active bacteria of the genus Polaromonas (class Betaproteobacteria) are important components of glacial microbiomes. In this study, extrachromosomal replicons of 26 psychrotolerant Polaromonas strains, isolated from Arctic and Antarctic glaciers, were identified, sequenced, and characterized. The plasmidome of these strains consists of 13 replicons, ranging in size from 3,378 to 101,077 bp. In silico sequence analyses identified the conserved backbones of these plasmids, composed of genes required for plasmid replication, stable maintenance, and conjugal transfer. Host range analysis revealed that all of the identified plasmids are narrow-host-range replicons, only able to replicate in bacteria of closely related genera (Polaromonas and Variovorax) of the Comamonadaceae family. Special attention was paid to the identification of plasmid auxiliary genetic information, which may contribute to the adaptation of bacteria to environmental conditions occurring in glaciers. Detailed analysis revealed the presence of genes encoding proteins potentially involved in (i) protection against reactive oxygen species, ultraviolet radiation, and low temperatures; (ii) transport and metabolism of organic compounds; (iii) transport of metal ions; and (iv) resistance to heavy metals. Some of the plasmids also carry genes required for the molecular assembly of iron-sulfur [Fe-S] clusters. Functional analysis of the predicted heavy metal resistance determinants demonstrated that their activity varies, depending on the host strain. This study provides the first molecular insight into the mobile DNA of Polaromonas spp. inhabiting polar glaciers. It has generated valuable data on the structure and properties of a pool of plasmids and highlighted their role in the biology of psychrotolerant Polaromonas strains and their adaptation to the environmental conditions of Arctic and Antarctic glaciers.

Diversity and Assembling Processes of Bacterial Communities in Cryoconite Holes of a Karakoram Glacier

Cryoconite holes are small ponds that form on the surface of glaciers that contain a dark debris, the cryoconite, at the bottom and host active ecological communities. Differences in the structure of bacterial communities have been documented among Arctic and mountain glaciers, and among glaciers in different areas of the world. In this study, we investigated the structure of bacterial communities of cryoconite holes of Baltoro Glacier, a large (62 km in length and 524 km² of surface) glacier of the Karakoram, by high-throughput sequencing of the V5-V6 hypervariable regions of the 16S rRNA gene. We found that Betaproteobacteria dominated bacterial communities, with large abundance of genera Polaromonas, probably thanks to its highly versatile metabolism, and Limnohabitans, which may have been favoured by the presence of supraglacial lakes in the area where cryoconite holes were sampled. Variation in bacterial communities among different sampling areas of the glacier could be explained by divergent selective processes driven by variation in environmental conditions, particularly pH, which was the only environmental variable that significantly affected the structure of bacterial communities. This variability may be due to both temporal and spatial patterns of variation in environmental conditions.

The ICEXTD of Azoarcus sp. CIB, an integrative and conjugative element with aerobic and anaerobic catabolic properties

Integrative and conjugative elements (ICE) play a major role in aerobic degradation of aromatic compounds, but they have not yet been shown to be involved in anaerobic degradation. We have characterized here the ICE_XTD element which endows to the beta-proteobacterium Azoarcus sp. CIB with the ability to utilize aromatic hydrocarbons. The core region of ICE_XTD , which shows a remarkable synteny with that of ICEclc-like elements, allows its own intracellular and intercellular mobility. ICE_XTD integrates at the tRNAGly of the host chromosome, but it can also excise to produce a ready to transfer circular form. The adaptation modules of ICE_XTD represent a unique combination of gene clusters for aerobic (tod genes) and anaerobic (bss-bbs and mbd genes) degradation of certain aromatic hydrocarbons, e.g., toluene, m-xylene and cumene. Transfer of ICE_XTD to other Azoarcus strains, e.g., A. evansii, confers them the ability to degrade aromatic hydrocarbons both aerobically and anaerobically. Interestingly, ICE_XTD allows Cupriavidus pinatubonensis, a bacterium unable to degrade anaerobically aromatic compounds, to grow with m-xylene under anoxic conditions. Thus, ICE_XTD constitutes the first mobile genetic element able to expand the catabolic abilities of certain bacteria for the removal of aromatic hydrocarbons either in the presence or absence of oxygen.

Comparative Genomics and Metabolic Analysis Reveals Peculiar Characteristics of Rhodococcus opacus Strain M213 Particularly for Naphthalene Degradation

The genome of Rhodococcus opacus strain M213, isolated from a fuel-oil contaminated soil, was sequenced and annotated which revealed a genome size of 9,194,165 bp encoding 8680 putative genes and a G+C content of 66.72%. Among the protein coding genes, 71.77% were annotated as clusters of orthologous groups of proteins (COGs); 55% of the COGs were present as paralog clusters. Pulsed field gel electrophoresis (PFGE) analysis of M213 revealed the presence of three different sized replicons- a circular chromosome and two megaplasmids (pNUO1 and pNUO2) estimated to be of 750Kb 350Kb in size, respectively. Conversely, using an alternative approach of optical mapping, the plasmid replicons appeared as a circular ~1.2 Mb megaplasmid and a linear, ~0.7 Mb megaplasmid. Genome-wide comparative analysis of M213 with a cohort of sequenced Rhodococcus species revealed low syntenic affiliation with other R. opacus species including strains B4 and PD630. Conversely, a closer affiliation of M213, at the functional (COG) level, was observed with the catabolically versatile R. jostii strain RHA1 and other Rhodococcii such as R. wratislaviensis strain IFP 2016, R. imtechensis strain RKJ300, Rhodococcus sp. strain JVH1, and Rhodococcus sp. strain DK17, respectively. An in-depth, genome-wide comparison between these functional relatives revealed 971 unique genes in M213 representing 11% of its total genome; many associating with catabolic functions. Of major interest was the identification of as many as 154 genomic islands (GEIs), many with duplicated catabolic genes, in particular for PAHs; a trait that was confirmed by PCR-based identification of naphthalene dioxygenase (NDO) as a representative gene, across PFGE-resolved replicons of strain M213. Interestingly, several plasmid/GEI-encoded genes, that likely participate in degrading naphthalene (NAP) via a peculiar pathway, were also identified in strain M213 using a combination of bioinformatics, metabolic analysis and gene expression measurements of selected catabolic genes by RT-PCR. Taken together, this study provides a comprehensive understanding of the genome plasticity and ecological competitiveness of strain M213 likely facilitated by horizontal gene transfer (HGT), bacteriophage attacks and genomic reshuffling- aspects that continue to be understudied and thus poorly understood, in particular for the soil-borne Rhodococcii.

Linking Microbial Community and Catabolic Gene Structures during the Adaptation of Three Contaminated Soils under Continuous Long-Term Pollutant Stress

Three types of contaminated soil from three geographically different areas were subjected to a constant supply of benzene or benzene/toluene/ethylbenzene/xylenes (BTEX) for a period of 3 months. Different from the soil from Brazil (BRA) and Switzerland (SUI), the Czech Republic (CZE) soil which was previously subjected to intensive in situ bioremediation displayed only negligible changes in community structure. BRA and SUI soil samples showed a clear succession of phylotypes. A rapid response to benzene stress was observed, whereas the response to BTEX pollution was significantly slower. After extended incubation, actinobacterial phylotypes increased in relative abundance, indicating their superior fitness to pollution stress. Commonalities but also differences in the phylotypes were observed. Catabolic gene surveys confirmed the enrichment of actinobacteria by identifying the increase of actinobacterial genes involved in the degradation of pollutants. Proteobacterial phylotypes increased in relative abundance in SUI microcosms after short-term stress with benzene, and catabolic gene surveys indicated enriched metabolic routes. Interestingly, CZE soil, despite staying constant in community structure, showed a change in the catabolic gene structure. This indicates that a highly adapted community, which had to adjust its gene pool to meet novel challenges, has been enriched.

Identification of benzo[a]pyrene-metabolizing bacteria in forest soils by using DNA-based stable-isotope probing

DNA-based stable-isotope probing (DNA-SIP) was used in this study to investigate the uncultivated bacteria with benzo[a]pyrene (BaP) metabolism capacities in two Chinese forest soils (Mt. Maoer in Heilongjiang Province and Mt. Baicaowa in Hubei Province). We characterized three different phylotypes with responsibility for BaP degradation, none of which were previously reported as BaP-degrading microorganisms by SIP. In Mt. Maoer soil microcosms, the putative BaP degraders were classified as belonging to the genus Terrimonas (family Chitinophagaceae, order Sphingobacteriales), whereas Burkholderia spp. were the key BaP degraders in Mt. Baicaowa soils. The addition of metabolic salicylate significantly increased BaP degradation efficiency in Mt. Maoer soils, and the BaP-metabolizing bacteria shifted to the microorganisms in the family Oxalobacteraceae (genus unclassified). Meanwhile, salicylate addition did not change either BaP degradation or putative BaP degraders in Mt. Baicaowa. Polycyclic aromatic hydrocarbon ring-hydroxylating dioxygenase (PAH-RHD) genes were amplified, sequenced, and quantified in the DNA-SIP (13)C heavy fraction to further confirm the BaP metabolism. By illuminating the microbial diversity and salicylate additive effects on BaP degradation across different soils, the results increased our understanding of BaP natural attenuation and provided a possible approach to enhance the bioremediation of BaP-contaminated soils.

A novel family of integrases associated with prophages and genomic islands integrated within the tRNA-dihydrouridine synthase A (dusA) gene

Genomic islands play a key role in prokaryotic genome plasticity. Genomic islands integrate into chromosomal loci such as transfer RNA genes and protein coding genes, whilst retaining various cargo genes that potentially bestow novel functions on the host organism. A gene encoding a putative integrase was identified at a single site within the 5′ end of the dusA gene in the genomes of over 200 bacteria. This integrase was discovered to be a component of numerous genomic islands, which appear to share a target site within the dusA gene. dusA encodes the tRNA-dihydrouridine synthase A enzyme, which catalyses the post-transcriptional reduction of uridine to dihydrouridine in tRNA. Genomic islands encoding homologous dusA-associated integrases were found at a much lower frequency within the related dusB and dusC genes, and non-dus genes. Excision of these dusA-associated islands from the chromosome as circularized intermediates was confirmed by polymerase chain reaction. Analysis of the dusA-associated islands indicated that they were highly diverse, with the integrase gene representing the only universal common feature.

Genomics of microbial plasmids: classification and identification based on replication and transfer systems and host taxonomy

Plasmids are important "vehicles" for the communication of genetic information between bacteria. The exchange of plasmids transmits pathogenically and environmentally relevant traits to the host bacteria, promoting their rapid evolution and adaptation to various environments. Over the past six decades, a large number of plasmids have been identified and isolated from different microbes. With the revolution of sequencing technology, more than 4600 complete sequences of plasmids found in bacteria, archaea, and eukaryotes have been determined. The classification of a wide variety of plasmids is not only important to understand their features, host ranges, and microbial evolution but is also necessary to effectively use them as genetic tools for microbial engineering. This review summarizes the current situation of the classification of fully sequenced plasmids based on their host taxonomy and their features of replication and conjugative transfer. The majority of the fully sequenced plasmids are found in bacteria in the Proteobacteria, Firmicutes, Spirochaetes, Actinobacteria, Cyanobacteria and Euryarcheota phyla, and key features of each phylum are included. Recent advances in the identification of novel types of plasmids and plasmid transfer by culture-independent methods using samples from natural environments are also discussed.

High correlation between genotypes and phenotypes of environmental bacteria Comamonas testosteroni strains

BACKGROUND

Members of Comamonas testosteroni are environmental microorganisms that are usually found in polluted environment samples. They utilize steroids and aromatic compounds but rarely sugars, and show resistance to multiple heavy metals and multiple drugs. However, comprehensive genomic analysis among the C. testosteroni strains is lacked.

RESULTS

To understand the genome bases of the features of C. testosteroni, we sequenced 10 strains of this species and analyzed them together with other related published genome sequences. The results revealed that: 1) the strains of C. testosteroni have genome sizes ranging from 5.1 to 6.0 Mb and G + C contents ranging from 61.1% to 61.8%. The pan-genome contained 10,165 gene families and the core genome contained 3,599 gene families. Heap's law analysis indicated that the pan-genome of C. testosteroni may be open (α = 0.639); 2) by analyzing 31 phenotypes of 11 available C. testosteroni strains, 99.4% of the genotypes (putative genes) were found to be correlated to the phenotypes, indicating a high correlation between phenotypes and genotypes; 3) gene clusters for nitrate reduction, steroids degradation and metal and multi-drug resistance were found and were highly conserved among all the genomes of this species; 4) the genome similarity of C. testosteroni may be related to the geographical distances.

CONCLUSIONS

This work provided an overview on the genomes of C. testosteroni and new genome resources that would accelerate the further investigations of this species. Importantly, this work focused on the analysis of potential genetic determinants for the typical characters and found high correlation between the phenotypes and their corresponding genotypes.

Genomic and metabolic analysis of fluoranthene degradation pathway in Celeribacter indicus P73T

Celeribacter indicus P73(T), isolated from deep-sea sediment from the Indian Ocean, is capable of degrading a wide range of polycyclic aromatic hydrocarbons (PAHs) and is the first fluoranthene-degrading bacterium within the family Rhodobacteraceae. Here, the complete genome sequence of strain P73(T) is presented and analyzed. Besides a 4.5-Mb circular chromosome, strain P73(T) carries five plasmids, and encodes 4827 predicted protein-coding sequences. One hundred and thirty-eight genes, including 14 dioxygenase genes, were predicted to be involved in the degradation of aromatic compounds, and most of these genes are clustered in four regions. P73_0346 is the first fluoranthene 7,8-dioxygenase to be discovered and the first fluoranthene dioxygenase within the toluene/biphenyl family. The degradative genes in regions B and D in P73(T) are absent in Celeribacter baekdonensis B30, which cannot degrade PAHs. Four intermediate metabolites [acenaphthylene-1(2H)-one, acenaphthenequinone, 1,2-dihydroxyacenaphthylene, and 1,8-naphthalic anhydride] of fluoranthene degradation by strain P73(T) were detected as the main intermediates, indicating that the degradation of fluoranthene in P73(T) was initiated by dioxygenation at the C-7,8 positions. Based on the genomic and metabolitic results, we propose a C-7,8 dioxygenation pathway in which fluoranthene is mineralized to TCA cycle intermediates.

Draft Genome Sequence of a Metabolically Diverse Antarctic Supraglacial Stream Organism, Polaromonas sp. Strain CG9_12, Determined Using Pacific Biosciences Single-Molecule Real-Time Sequencing Technology

Polaromonas species are found in a diversity of environments and are particularly common in icy ecosystems. Polaromonas sp. strain CG9_12 is an aerobic, Gram-negative, catalase-positive, white-pigmented bacterium of the Proteobacteria phylum. Here, we present the draft genome sequence of Polaromonas sp. strain CG9_12, isolated from an Antarctic supraglacial stream.

Vertical distribution of bacterial communities in high arsenic sediments of Hetao Plain, Inner Mongolia

Vertical distribution of bacterial communities was detected in high arsenic (As) sediments in a representative high As area in Inner Mongolia. Nineteen sediment samples were collected from a 30 m borehole and detected by geochemistry and molecular ecological approaches including polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE), 16S rRNA gene clone library and 454 pyrosequencing. As contents ranged from 42.1 to 111.3 mg kg(-1) which fluctuated with different depth and significantly high in clay and mild clay sediment samples at depth of 8, 20, 25 and 28 m respectively. The ratios of As(III) to total As generally increased with depth but As(V) dominated in all sediment samples. High concentrations of total As, sulfur, iron and total organic carbon were generally found in clay and low in sand samples. Both DGGE patterns and 454 pyrosequencing results indicated that bacterial communities dynamically diversified with increasing depth and were dominated by Firmicutes, Bacteroidetes, Proteobacteria and Chloroflexi. Most of the sediment samples were dominated by populations including Sporosarcina, Acinetobacter, Pseudomonas, Halomonas, Polaromonas, Paenibacillus and Flavobacterium. These populations were found with high similarities with those microbes capable of denitrification, sulfur oxidation, organic matter degradation and As resistance and reduction. These results implied that microbes might play an important role in As mobilization in the shallow aquifers of Hetao Plain, Inner Mongolia.

Ecology drives the distribution of specialized tyrosine metabolism modules in fungi

Gene clusters encoding accessory or environmentally specialized metabolic pathways likely play a significant role in the evolution of fungal genomes. Two such gene clusters encoding enzymes associated with the tyrosine metabolism pathway (KEGG #00350) have been identified in the filamentous fungus Aspergillus fumigatus. The l-tyrosine degradation (TD) gene cluster encodes a functional module that facilitates breakdown of the phenolic amino acid, l-tyrosine through a homogentisate intermediate, but is also involved in the production of pyomelanin, a fungal pathogenicity factor. The gentisate catabolism (GC) gene cluster encodes a functional module likely involved in phenolic compound degradation, which may enable metabolism of biphenolic stilbenes in multiple lineages. Our investigation of the evolution of the TD and GC gene clusters in 214 fungal genomes revealed spotty distributions partially shaped by gene cluster loss and horizontal gene transfer (HGT). Specifically, a TD gene cluster shows evidence of HGT between the extremophilic, melanized fungi Exophiala dermatitidis and Baudoinia compniacensis, and a GC gene cluster shows evidence of HGT between Sordariomycete and Dothideomycete grass pathogens. These results suggest that the distribution of specialized tyrosine metabolism modules is influenced by both the ecology and phylogeny of fungal species.

Correction of sequence-dependent ambiguous bases (Ns) from the 454 pyrosequencing system

Pyrosequencing of the 16S ribosomal RNA gene (16S) has become one of the most popular methods to assess microbial diversity. Pyrosequencing reads containing ambiguous bases (Ns) are generally discarded based on the assumptions of their non-sequence-dependent formation and high error rates. However, taxonomic composition differed by removal of reads with Ns. We determined whether Ns from pyrosequencing occur in a sequence-dependent manner. Our reads and the corresponding flow value data revealed occurrence of sequence-specific N errors with a common sequential pattern (a homopolymer + a few nucleotides with bases other than the homopolymer + N) and revealed that the nucleotide base of the homopolymer is the true base for the following N. Using an algorithm reflecting this sequence-dependent pattern, we corrected the Ns in the 16S (86.54%), bphD (81.37%) and nifH (81.55%) amplicon reads from a mock community with high precisions of 95.4, 96.9 and 100%, respectively. The new N correction method was applicable for determining most of Ns in amplicon reads from a soil sample, resulting in reducing taxonomic biases associated with N errors and in shotgun sequencing reads from public metagenome data. The method improves the accuracy and precision of microbial community analysis and genome sequencing using 454 pyrosequencing.

Architecture and gene repertoire of the flexible genome of the extreme acidophile Acidithiobacillus caldus

BACKGROUND

Acidithiobacillus caldus is a sulfur oxidizing extreme acidophile and the only known mesothermophile within the Acidithiobacillales. As such, it is one of the preferred microbes for mineral bioprocessing at moderately high temperatures. In this study, we explore the genomic diversity of A. caldus strains using a combination of bioinformatic and experimental techniques, thus contributing first insights into the elucidation of the species pangenome.

PRINCIPAL FINDINGS

Comparative sequence analysis of A. caldus ATCC 51756 and SM-1 indicate that, despite sharing a conserved and highly syntenic genomic core, both strains have unique gene complements encompassing nearly 20% of their respective genomes. The differential gene complement of each strain is distributed between the chromosomal compartment, one megaplasmid and a variable number of smaller plasmids, and is directly associated to a diverse pool of mobile genetic elements (MGE). These include integrative conjugative and mobilizable elements, genomic islands and insertion sequences. Some of the accessory functions associated to these MGEs have been linked previously to the flexible gene pool in microorganisms inhabiting completely different econiches. Yet, others had not been unambiguously mapped to the flexible gene pool prior to this report and clearly reflect strain-specific adaption to local environmental conditions.

SIGNIFICANCE

For many years, and because of DNA instability at low pH and recurrent failure to genetically transform acidophilic bacteria, gene transfer in acidic environments was considered negligible. Findings presented herein imply that a more or less conserved pool of actively excising MGEs occurs in the A. caldus population and point to a greater frequency of gene exchange in this econiche than previously recognized. Also, the data suggest that these elements endow the species with capacities to withstand the diverse abiotic and biotic stresses of natural environments, in particular those associated with its extreme econiche.

Analysis of metabolic evolution in bacteria using whole-genome metabolic models

Recent advances in the automation of metabolic model reconstruction have led to the availability of draft-quality metabolic models (predicted reaction complements) for multiple bacterial species. These reaction complements can be considered as trait representations and can be used for ancestral state reconstruction to infer the most likely metabolic complements of common ancestors of all bacteria with generated metabolic models. We present here an ancestral state reconstruction for 141 extant bacteria and analyze the reaction gains and losses for these bacteria with respect to their lifestyles and pathogenic nature. A simulated annealing approach is used to look at coordinated metabolic gains and losses in two bacteria. The main losses of Onion yellows phytoplasma OY-M, an obligate intracellular pathogen, are shown (as expected) to be in cell wall biosynthesis. The metabolic gains made by Clostridium difficile CD196 in adapting to its current habitat in the human colon is also analyzed. Our analysis shows that the capability to utilize N-Acetyl-neuraminic acid as a carbon source has been gained, rather than having been present in the Clostridium ancestor, as has the capability to synthesize phthiocerol dimycocerosate, which could potentially aid the evasion of the host immune response. We have shown that the availability of large numbers of metabolic models, along with conventional approaches, has enabled a systematic method to analyze metabolic evolution in the bacterial domain.

Novel tripartite aromatic acid transporter essential for terephthalate uptake in Comamonas sp. strain E6

It has been suggested that a novel type of aromatic acid transporter, which is similar to the tripartite tricarboxylate transporter (TTT), is involved in terephthalate (TPA) uptake by Comamonas sp. strain E6. This suggestion was based on the presence of the putative TPA-binding protein gene, tphC, in the TPA catabolic operon. The tphC gene is essential for growth on TPA and is similar to the genes encoding TTT-like substrate-binding proteins. Here we identified two sets of E6 genes, tctBA and tpiBA, which encode TTT-like cytoplasmic transmembrane proteins. Disruption of tctA showed no influence on TPA uptake but resulted in a complete loss of the uptake of citrate. This loss suggests that tctA is involved in citrate uptake. On the other hand, disruption of tpiA or tpiB demonstrated that both genes are essential for TPA uptake. Only when both tphC and tpiBA were introduced with the TPA catabolic genes into cells of a non-TPA-degrading Pseudomonas strain did the resting cells of the transformant acquire the ability to convert TPA. From all these results, it was concluded that the TPA uptake system consists of the TpiA-TpiB membrane components and TPA-binding TphC. Interestingly, not only was the tpiA mutant of E6 unable to grow on TPA or isophthalate, it also showed significant growth delays on o-phthalate and protocatechuate. These results suggested that the TpiA-TpiB membrane components are able to interact with multiple substrate-binding proteins. The tpiBA genes were constitutively transcribed as a single operon in E6 cells, whereas the transcription of tphC was positively regulated by TphR. TPA uptake by E6 cells was completely inhibited by a protonophore, carbonyl cyanide m-chlorophenyl hydrazone, indicating that the TPA uptake system requires a proton motive force.

Bacterial communities in tetrachloroethene-polluted groundwaters: a case study

The compositions of bacterial groundwater communities of three sites contaminated with chlorinated ethenes were analyzed by pyrosequencing their 16S rRNA genes. For each location, the entire and the active bacterial populations were characterized by independent molecular analysis of the community DNA and RNA. The sites were selected to cover a broad range of different environmental conditions and contamination levels, with tetrachloroethene (PCE) and trichloroethene (TCE) being the primary contaminants. Before sampling the biomass, a long-term monitoring of the polluted locations revealed high concentrations of cis-1,2-dichloroethene (cDCE) and vinyl chloride (VC), which are toxic by-products of the incomplete bacterial degradation of PCE and TCE. The applied pyrosequencing technique enabled known dechlorinators to be identified at a very low detection level (<0.25%) without compromising the detailed analysis of the entire bacterial community of these sites. The study revealed that only a few species dominated the bacterial communities, with Albidiferax ferrireducens being the only highly prominent member found at all three sites. Only a limited number of OTUs with abundances of up to 1% and high sequence identities to known dechlorinating microorganisms were retrieved from the RNA pools of the two highly contaminated sites. The dechlorinating consortium was likely to be comprised of cDCE-assimilating bacteria (Polaromonas spp.), anaerobic organohalide respirers (mainly Geobacter spp.), and Burkholderia spp. involved in cometabolic dechlorination processes, together with methylotrophs (Methylobacter spp.). The deep sequencing results suggest that the indigenous dechlorinating consortia present at the investigated sites can be used as a starting point for future bioremediation activities by stimulating their anaerobic and aerobic chloroethene degradation capacities (i.e. reductive dechlorination, and metabolic and cometabolic oxidation).

Denitrifying bacterial communities affect current production and nitrous oxide accumulation in a microbial fuel cell

The biocathodic reduction of nitrate in Microbial Fuel Cells (MFCs) is an alternative to remove nitrogen in low carbon to nitrogen wastewater and relies entirely on microbial activity. In this paper the community composition of denitrifiers in the cathode of a MFC is analysed in relation to added electron acceptors (nitrate and nitrite) and organic matter in the cathode. Nitrate reducers and nitrite reducers were highly affected by the operational conditions and displayed high diversity. The number of retrieved species-level Operational Taxonomic Units (OTUs) for narG, napA, nirS and nirK genes was 11, 10, 31 and 22, respectively. In contrast, nitrous oxide reducers remained virtually unchanged at all conditions. About 90% of the retrieved nosZ sequences grouped in a single OTU with a high similarity with Oligotropha carboxidovorans nosZ gene. nirS-containing denitrifiers were dominant at all conditions and accounted for a significant amount of the total bacterial density. Current production decreased from 15.0 A·m⁻³ NCC (Net Cathodic Compartment), when nitrate was used as an electron acceptor, to 14.1 A·m⁻³ NCC in the case of nitrite. Contrarily, nitrous oxide (N₂O) accumulation in the MFC was higher when nitrite was used as the main electron acceptor and accounted for 70% of gaseous nitrogen. Relative abundance of nitrite to nitrous oxide reducers, calculated as (qnirS+qnirK)/qnosZ, correlated positively with N₂O emissions. Collectively, data indicate that bacteria catalysing the initial denitrification steps in a MFC are highly influenced by main electron acceptors and have a major influence on current production and N₂O accumulation.

The dynamic bacterial communities of a melting High Arctic glacier snowpack

Snow environments can occupy over a third of land surface area, but little is known about the dynamics of snowpack bacteria. The effect of snow melt on bacterial community structure and diversity of surface environments of a Svalbard glacier was examined using analyses of 16S rRNA genes via T-RFLP, qPCR and 454 pyrosequencing. Distinct community structures were found in different habitat types, with changes over 1 week apparent, in particular for the dominant bacterial class present, Betaproteobacteria. The differences observed were consistent with influences from depositional mode (snowfall vs aeolian dusts), contrasting snow with dust-rich snow layers and near-surface ice. Contrary to that, slush as the decompositional product of snow harboured distinct lineages of bacteria, further implying post-depositional changes in community structure. Taxa affiliated to the betaproteobacterial genus Polaromonas were particularly dynamic, and evidence for the presence of betaproteobacterial ammonia-oxidizing bacteria was uncovered, inviting the prospect that the dynamic bacterial communities associated with snowpacks may be active in supraglacial nitrogen cycling and capable of rapid responses to changes induced by snowmelt. Furthermore the potential of supraglacial snowpack ecosystems to respond to transient yet spatially extensive melting episodes such as that observed across most of Greenland's ice sheet in 2012 merits further investigation.

Genomic and functional analyses of the gentisate and protocatechuate ring-cleavage pathways and related 3-hydroxybenzoate and 4-hydroxybenzoate peripheral pathways in Burkholderia xenovorans LB400

In this study, the gentisate and protocatechuate pathways in Burkholderia xenovorans LB400 were analyzed by genomic and functional approaches, and their role in 3-hydroxybenzoate (3-HBA) and 4-hydroxybenzoate (4-HBA) degradation was proposed. The LB400 genome possesses two identical mhbRTDHI gene clusters encoding the gentisate pathway and one mhbM gene encoding a 3-HBA 6-hydroxylase that converts 3-HBA into gentisate. The pca genes encoding the protocatechuate pathway and the pobA gene encoding the 4-HBA 3-monooxygenase that oxidizes 4-HBA into protocatechuate are arranged in gene clusters and single genes mainly at the minor chromosome, but also at the major chromosome and the megaplasmid. Strain LB400 was able to grow on gentisate, protocatechuate, 3-HBA and 4-HBA. Transcriptional analyses showed that the mhbD gene encoding the gentisate 1,2-dioxygenase was expressed during growth on 3-HBA, 4-HBA and gentisate, whereas the pcaG gene encoding the protocatechuate 3,4-dioxygenase was expressed only during growth on 4-HBA and protocatechuate. The mhbM gene encoding the 3-HBA 6-hydroxylase was transcribed in strain LB400 during growth on HBAs, gentisate, protocatechuate and glucose. The pobA gene encoding the 4-HBA 3-monooxygenase was expressed during growth on HBAs and glucose. 3-HBA- and 4-HBA-grown LB400 cells showed gentisate 1,2-dioxygenase activity, whereas protocatechuate 3,4-dioxygenase activity was observed only in 4-HBA-grown cells. The mhbR gene encoding a MarR-type transcriptional regulator that probably regulates the expression of the MhbT transporter, and the pcaQ and pcaR genes encoding LysR-type transcriptional regulators that regulate pcaHG and pcaIJBDC genes, respectively, were transcribed during growth on both HBAs, gentisate, protocatechuate and glucose, suggesting a basal constitutive expression. The results indicate active gentisate, protocatechuate, 3-HBA and 4-HBA catabolic pathways in B. xenovorans LB400 and suggest that 3-HBA is channeled exclusively through the gentisate route, whereas 4-HBA is funneled into the protocatechuate central pathway and potentially into the gentisate pathway.

Dynamics of bacterial communities in two unpolluted soils after spiking with phenanthrene: soil type specific and common responders

Considering their key role for ecosystem processes, it is important to understand the response of microbial communities in unpolluted soils to pollution with polycyclic aromatic hydrocarbons (PAH). Phenanthrene, a model compound for PAH, was spiked to a Cambisol and a Luvisol soil. Total community DNA from phenanthrene-spiked and control soils collected on days 0, 21, and 63 were analyzed based on PCR-amplified 16S rRNA gene fragments. Denaturing gradient gel electrophoresis (DGGE) fingerprints of bacterial communities increasingly deviated with time between spiked and control soils. In taxon specific DGGE, significant responses of Alphaproteobacteria and Actinobacteria became only detectable after 63 days, while significant effects on Betaproteobacteria were detectable in both soils after 21 days. Comparison of the taxonomic distribution of bacteria in spiked and control soils on day 63 as revealed by pyrosequencing indicated soil type specific negative effects of phenanthrene on several taxa, many of them belonging to the Gamma-, Beta-, or Deltaproteobacteria. Bacterial richness and evenness decreased in spiked soils. Despite the significant differences in the bacterial community structure between both soils on day 0, similar genera increased in relative abundance after PAH spiking, especially Sphingomonas and Polaromonas. However, this did not result in an increased overall similarity of the bacterial communities in both soils.

Distribution of nitrogen fixation and nitrogenase-like sequences amongst microbial genomes

Background

The metabolic capacity for nitrogen fixation is known to be present in several prokaryotic species scattered across taxonomic groups. Experimental detection of nitrogen fixation in microbes requires species-specific conditions, making it difficult to obtain a comprehensive census of this trait. The recent and rapid increase in the availability of microbial genome sequences affords novel opportunities to re-examine the occurrence and distribution of nitrogen fixation genes. The current practice for computational prediction of nitrogen fixation is to use the presence of the nifH and/or nifD genes.

Results

Based on a careful comparison of the repertoire of nitrogen fixation genes in known diazotroph species we propose a new criterion for computational prediction of nitrogen fixation: the presence of a minimum set of six genes coding for structural and biosynthetic components, namely NifHDK and NifENB. Using this criterion, we conducted a comprehensive search in fully sequenced genomes and identified 149 diazotrophic species, including 82 known diazotrophs and 67 species not known to fix nitrogen. The taxonomic distribution of nitrogen fixation in Archaea was limited to the Euryarchaeota phylum; within the Bacteria domain we predict that nitrogen fixation occurs in 13 different phyla. Of these, seven phyla had not hitherto been known to contain species capable of nitrogen fixation. Our analyses also identified protein sequences that are similar to nitrogenase in organisms that do not meet the minimum-gene-set criteria. The existence of nitrogenase-like proteins lacking conserved co-factor ligands in both diazotrophs and non-diazotrophs suggests their potential for performing other, as yet unidentified, metabolic functions.

Conclusions

Our predictions expand the known phylogenetic diversity of nitrogen fixation, and suggest that this trait may be much more common in nature than it is currently thought. The diverse phylogenetic distribution of nitrogenase-like proteins indicates potential new roles for anciently duplicated and divergent members of this group of enzymes.

Comparative genomics reveals adaptation by Alteromonas sp. SN2 to marine tidal-flat conditions: cold tolerance and aromatic hydrocarbon metabolism

Alteromonas species are globally distributed copiotrophic bacteria in marine habitats. Among these, sea-tidal flats are distinctive: undergoing seasonal temperature and oxygen-tension changes, plus periodic exposure to petroleum hydrocarbons. Strain SN2 of the genus Alteromonas was isolated from hydrocarbon-contaminated sea-tidal flat sediment and has been shown to metabolize aromatic hydrocarbons there. Strain SN2's genomic features were analyzed bioinformatically and compared to those of Alteromonas macleodii ecotypes: AltDE and ATCC 27126. Strain SN2's genome differs from that of the other two strains in: size, average nucleotide identity value, tRNA genes, noncoding RNAs, dioxygenase gene content, signal transduction genes, and the degree to which genes collected during the Global Ocean Sampling project are represented. Patterns in genetic characteristics (e.g., GC content, GC skew, Karlin signature, CRISPR gene homology) indicate that strain SN2's genome architecture has been altered via horizontal gene transfer (HGT). Experiments proved that strain SN2 was far more cold tolerant, especially at 5°C, than the other two strains. Consistent with the HGT hypothesis, a total of 15 genomic islands in strain SN2 likely confer ecological fitness traits (especially membrane transport, aromatic hydrocarbon metabolism, and fatty acid biosynthesis) specific to the adaptation of strain SN2 to its seasonally cold sea-tidal flat habitat.

Role of nitrogen fixation in the autecology of Polaromonas naphthalenivorans in contaminated sediments

Polaromonas naphthalenivorans strain CJ2 is a Gram-negative betaproteobacterium that was identified, using stable isotope probing in 2003, as a dominant in situ degrader of naphthalene in coal tar-contaminated sediments. The sequenced genome of strain CJ2 revealed several genes conferring nitrogen fixation within a 65.6 kb region of strain CJ2's chromosome that is absent in the genome of its closest sequenced relative Polaromonas sp. strain JS666. Laboratory growth and nitrogenase assays verified that these genes are functional, providing an alternative source of nitrogen in N-free media when using naphthalene or pyruvate as carbon sources. Knowing this, we investigated if nitrogen-fixation activity could be detected in microcosms containing sediments from the field site where strain CJ2 was isolated. Inducing nitrogen limitation with the addition of glucose or naphthalene stimulated nitrogenase activity in amended sediments, as detected using the acetylene reduction assay. With the use of fluorescence microscopy, we screened the microcosm sediments for the presence of active strain CJ2 cells using a dual-labelling approach. When we examined the carbon-amended microcosm sediments stained with both a strain CJ2-specific fluorescent in situ hybridization probe and a polyclonal fluorescently tagged antibody, we were able to detect dual-labelled active cells. In contrast, in sediments that received no carbon addition (showing no nitrogenase activity), no dual-labelled cells were detected. Furthermore, the naphthalene amendment enhanced the proportion of active strain CJ2 cells in the sediment relative to a glucose amendment. Field experiments performed in sediments where strain CJ2 was isolated showed nitrogenase activity in response to dosing with naphthalene. Dual-label fluorescence staining of these sediments showed a fivefold increase in active strain CJ2 in the sediments dosed with naphthalene over those dosed with deionized water. These experiments show that nitrogen fixation may play an important role in naphthalene biodegradation by strain CJ2 and contribute to its ecological success.

The phn Island: A New Genomic Island Encoding Catabolism of Polynuclear Aromatic Hydrocarbons

Bacteria are key in the biodegradation of polycyclic aromatic hydrocarbons (PAH), which are widespread environmental pollutants. At least six genotypes of PAH degraders are distinguishable via phylogenies of the ring-hydroxylating dioxygenase (RHD) that initiates bacterial PAH metabolism. A given RHD genotype can be possessed by a variety of bacterial genera, suggesting horizontal gene transfer (HGT) is an important process for dissemination of PAH-degrading genes. But, mechanisms of HGT for most RHD genotypes are unknown. Here, we report in silico and functional analyses of the phenanthrene-degrading bacterium Delftia sp. Cs1-4, a representative of the phn_AFK2 RHD group. The phn_AFK2 genotype predominates PAH degrader communities in some soils and sediments, but, until now, their genomic biology has not been explored. In the present study, genes for the entire phenanthrene catabolic pathway were discovered on a novel ca. 232 kb genomic island (GEI), now termed the phn island. This GEI had characteristics of an integrative and conjugative element with a mobilization/stabilization system similar to that of SXT/R391-type GEI. But, it could not be grouped with any known GEI, and was the first member of a new GEI class. The island also carried genes predicted to encode: synthesis of quorum sensing signal molecules, fatty acid/polyhydroxyalkanoate biosynthesis, a type IV secretory system, a PRTRC system, DNA mobilization functions and >50 hypothetical proteins. The 50% G + C content of the phn gene cluster differed significantly from the 66.7% G + C level of the island as a whole and the strain Cs1-4 chromosome, indicating a divergent phylogenetic origin for the phn genes. Collectively, these studies added new insights into the genetic elements affecting the PAH biodegradation capacity of microbial communities specifically, and the potential vehicles of HGT in general.

Horizontal transfer of PAH catabolism genes in Mycobacterium: evidence from comparative genomics and isolated pyrene-degrading bacteria

Biodegradation of high molecular weight polycyclic aromatic hydrocarbons (PAHs), such as pyrene and benzo[a]pyrene, has only been observed in a few genera, namely fast-growing Mycobacterium and Rhodococcus. In M. vanbaalenii PYR-1, multiple aromatic ring hydroxylating dioxygenase (ARHDOs) genes including pyrene dioxygenases nidAB and nidA3B3 are localized in one genomic region. Here we examine the homologous genomic regions in four other PAH-degrading Mycobacterium (strains JLS, KMS, and MCS, and M. gilvum PYR-GCK), presenting evidence for past horizontal gene transfer events. Seven distinct types of ARHDO genes are present in all five genomes, and display conserved syntenic architecture with respect to gene order, orientation, and association with other genes. Duplications and putative integrase and transposase genes suggest past gene shuffling. To corroborate these observations, pyrene-degrading strains were isolated from two PAH-contaminated sediments: Chattanooga Creek (Tennessee) and Lake Erie (western basin). Some were related to fast-growing Mycobacterium spp. and carried both nidA and nidA3 genes. Other isolates belonged to Microbacteriaceae and Intrasporangiaceae presenting the first evidence of pyrene degradation in these families. These isolates had nidA (and some, nidA3) genes that were homologous to Mycobacterial ARHDO genes, suggesting that horizontal gene transfer events have occurred.

Polaromonas glacialis sp. nov. and Polaromonas cryoconiti sp. nov., isolated from alpine glacier cryoconite

The taxonomic positions of two Gram-staining-negative, psychrophilic bacteria, which were isolated from alpine glacier cryoconite and designated strains Cr4-12(T) and Cr4-35(T), were investigated using a polyphasic approach. Both novel strains contained ubiquinone Q-8 as the sole quinone, summed feature 3 (C(16:1)ω7c and/or C(16:1)ω6c) and C(16:0) as the dominant cellular fatty acids, putrescine and 2-hydroxyputrescine as the major polyamines, and diphosphatidylglycerol, phosphatidylglycerol and phosphatidylethanolamine as the major polar lipids. The genomic DNA G+C contents of strains Cr4-12(T) and Cr4-35(T) were 61.3 mol% and 60.7 mol%, respectively. Phylogenetic analysis based on 16S rRNA gene sequences indicated that the two strains belonged to the genus Polaromonas. Although the 16S rRNA gene sequences of strains Cr4-12(T) and Cr4-35(T) were very similar (98.7% sequence similarity), hybridizations indicated a DNA-DNA relatedness value of only 26.9% between the two novel strains. In pairwise comparisons with the type strains of recognized Polaromonas species, strains Cr4-12(T) and Cr4-35(T) showed 16S rRNA gene sequence similarities of 96.4-98.5% and 96.5-98.4%, respectively. Based on the phenotypic and phylogenetic evidence and DNA-DNA relatedness data, strains Cr4-12(T) and Cr4-35(T) represent two novel species within the genus Polaromonas, for which the names Polaromonas glacialis sp. nov. and Polaromonas cryoconiti sp. nov., respectively, are proposed. The type strain of Polaromonas glacialis sp. nov. is Cr4-12(T) (=DSM 24062(T) =LMG 26049(T) =KACC 15089(T)) and that of Polaromonas cryoconiti sp. nov. is Cr4-35(T) (=DSM 24248(T) =LMG 26050(T) =KACC 15090(T)).

Diversity, biology and evolution of IncQ-family plasmids

Plasmids of IncQ-family are distinguished by having a unique strand-displacement mechanism of replication that is capable of functioning in a wide variety of bacterial hosts. In addition, these plasmids are highly mobilizable and therefore very promiscuous. Common features of the replicons have been used to identify IncQ-family plasmids in DNA sequence databases and in this way several unstudied plasmids have been compared to more well-studied IncQ plasmids. We propose that IncQ plasmids can be divided into four subgroups based on a number of mutually supportive criteria. The most important of these are the amino acid sequences of their three essential replication proteins and the observation that the replicon of each subgroup has become fused to four different lineages of mobilization genes. This review of IncQ-family plasmid diversity has highlighted several events in the evolution of these plasmids and raised several questions for further research.

Gentisate 1,2-dioxygenase, in the third naphthalene catabolic gene cluster of Polaromonas naphthalenivorans CJ2, has a role in naphthalene degradation

Polaromonas naphthalenivorans strain CJ2 metabolizes naphthalene via the gentisate pathway and has recently been shown to carry a third copy of gentisate 1,2-dioxygenase (GDO), encoded by nagI3, within a previously uncharacterized naphthalene catabolic gene cluster. The role of this cluster (especially nagI3) in naphthalene metabolism of strain CJ2 was investigated by documenting patterns in regulation, transcription and enzyme activity. Transcriptional analysis of wild-type cells showed the third cluster to be polycistronic and that nagI3 was expressed at a relatively high level. Individual knockout mutants of all three nagI genes were constructed and their influence on both GDO activity and cell growth was evaluated. Of the three knockout strains, CJ2ΔnagI3 showed severely diminished GDO activity and grew slowest on aromatic substrates. These observations are consistent with the hypothesis that nagI3 may prevent toxic intracellular levels of gentisate from accumulating in CJ2 cells. All three nagI genes from strain CJ2 were cloned into Escherichia coli: the nagI2 and nagI3 genes were successfully overexpressed. The subunit mass of the GDOs were ~36–39 kDa, and their structures were deduced to be dimeric. The K m values of NagI2 and NagI3 were 31 and 10 µM, respectively, indicating that the higher affinity of NagI3 for gentisate may protect the wild-type cells from gentisate toxicity. These results provide clues for explaining why the third gene cluster, particularly the nagI3 gene, is important in strain CJ2. The organization of genes in the third gene cluster matched that of clusters in Polaromonas sp. JS666 and Leptothrix cholodnii SP-6. While horizontal gene transfer (HGT) is one hypothesis for explaining this genetic motif, gene duplication within the ancestral lineage is equally valid. The HGT hypothesis was discounted by noting that the nagI3 allele of strain CJ2 did not share high sequence identity with its homologues in Polaromonas sp. JS666 and L. cholodnii SP-6.