Time-series metagenomic analysis reveals robustness of soil microbiome against chemical disturbance

Bacterial Diversity in Old Hydrocarbon Polluted Sediments of Ecuadorian Amazon River Basins

The Ecuadorian Amazon rainforest stands out as one of the world's most biodiverse regions, yet faces significant threats due to oil extraction activities dating back to the 1970s in the northeastern provinces. This research investigates the environmental and societal consequences of prolonged petroleum exploitation and oil spills in Ecuador's Amazon. Conducted in June 2015, the study involved a comprehensive analysis of freshwater sediment samples from 24 locations in the Rio Aguarico and Napo basins. Parameters such as water and air temperature, conductivity, soil pH, and hydrocarbon concentrations were examined. Total petroleum hydrocarbon (TPH) concentrations ranged from 9.4 to 847.4 mg kg-1, with polycyclic aromatic hydrocarbon (PAH) levels varying from 10.15 to 711.1 mg kg-1. The pristane/phytane ratio indicated historic hydrocarbon pollution in 8 of the 15 chemically analyzed sediments. Using non-culturable techniques (Illumina), bacterial analyses identified over 350 ASV, with prominent families including Comamonadaceae, Chitinophagaceae, Anaeromyxobacteraceae, Sphingomonadaceae, and Xanthobacteraceae. Bacterial diversity, assessed in eight samples, exhibited a positive correlation with PAH concentrations. The study provides insights into how microbial communities respond to varying levels of hydrocarbon pollution, shedding light on the enduring impact of oil exploitation in the Amazonian region. Its objective is to deepen our understanding of the environmental and human well-being in the affected area, underscoring the pressing need for remedial actions in the face of ongoing ecological challenges.

Reconstruction of a Soil Microbial Network Induced by Stress Temperature

The microbial community is viewed as a network of diverse microorganisms connected by various interspecific interactions. While the stress gradient hypothesis (SGH) predicts that positive interactions are favored in more stressful environments, the prediction has been less explored in complex microbial communities due to the challenges of identifying interactions. Here, by applying a nonlinear time series analysis to the amplicon-based diversity time series data of the soil microbiota cultured under less stressful (30°C) or more stressful (37°C) temperature conditions, we show how the microbial network responds to temperature stress. While the genera that persisted only under the less stressful condition showed fewer positive effects, the genera that appeared only under the more stressful condition received more positive effects, in agreement with SGH. However, temperature difference also induced reconstruction of the community network, leading to an increased proportion of negative interactions at the whole-community level. The anti-SGH pattern can be explained by the stronger competition caused by increased metabolic rate and population densities. IMPORTANCE By combining amplicon-based diversity survey with recently developed nonlinear analytical tools, we successfully determined the interaction networks of more than 150 natural soil microbial genera under less or more temperature stress and explored the applicability of the stress gradient hypothesis to soil microbiota, shedding new light on the well-known hypothesis.

Soil microbial communities-mediated bioattenuation in simulated aquifer storage and recovery (ASR) condition: Long-term study

This study evaluated the long-term organic removal performance and microbial community shift in simulated aquifer storage and recovery (ASR) conditions. For this purpose, anoxic soil box systems were operated at 15 °C for one year. The results showed that the assimilable organic carbon (AOC) concentration in the anoxic soil box systems was successfully decreased by 79.1%. The dissolved organic carbon (DOC) concentration increased during the initial operational periods; however, it subsequently decreased during long-term operation. Readily biodegradable organic fractions (i.e., low-molecular weight (LMW) neutrals and LMW acids) decreased along with time elapsed, whereas non-biodegradable fraction (i.e., humic substances) increased. Proteobacteria and Acidobacteriota were predominant in the anoxic box systems throughout the operational periods. Firmicutes and Bacteroidota suddenly increased during the initial operational period while Gemmatimonadota slightly increased during prolonged long-term operation. Interestingly, the microbial community structures were significantly shifted with respect to the operational periods while the effects of AOC/NO₃^- addition were negligible. Various bacterial species preferring low temperature or anoxic conditions were detected as predominant bacteria. Some denitrifying (i.e., Noviherbaspirillum denitrificans) and iron reducing bacteria (i.e., Geobacter spp.) appeared during the long-term operation; these bacterial communities also acted as organic degraders in the simulated ASR systems. The findings of this study suggest that the application of natural bioattenuation using indigenous soil microbial communities can be a promising option as an organic carbon management strategy in ASR systems.

Dynamic Changes in the Microbiome of Rice During Shoot and Root Growth Derived From Seeds

Microbes form close associations with host plants including rice as both surface (epiphytes) and internal (endophytes) inhabitants. Yet despite rice being one of the most important cereal crops agriculturally and economically, knowledge of its microbiome, particularly core inhabitants and any functional properties bestowed is limited. In this study, the microbiome in rice seedlings derived directly from seeds was identified, characterized and compared to the microbiome of the seed. Rice seeds were sourced from two different locations in Arkansas, USA of two different rice genotypes (Katy, M202) from two different harvest years (2013, 2014). Seeds were planted in sterile media and bacterial as well as fungal communities were identified through 16S and ITS sequencing, respectively, for four seedling compartments (root surface, root endosphere, shoot surface, shoot endosphere). Overall, 966 bacterial and 280 fungal ASVs were found in seedlings. Greater abundance and diversity were detected for the microbiome associated with roots compared to shoots and with more epiphytes than endophytes. The seedling compartments were the driving factor for microbial community composition rather than other factors such as rice genotype, location and harvest year. Comparison with datasets from seeds revealed that 91 (out of 296) bacterial and 11 (out of 341) fungal ASVs were shared with seedlings with the majority being retained within root tissues. Core bacterial and fungal microbiome shared across seedling samples were identified. Core bacteria genera identified in this study such as Rhizobium, Pantoea, Sphingomonas, and Paenibacillus have been reported as plant growth promoting bacteria while core fungi such as Pleosporales, Alternaria and Occultifur have potential as biocontrol agents.

Probabilistic model based on circular statistics for quantifying coverage depth dynamics originating from DNA replication

BACKGROUND

With the development of DNA sequencing technology, static omics profiling in microbial communities, such as taxonomic and functional gene composition determination, has become possible. Additionally, the recently proposed in situ growth rate estimation method allows the applicable range of current comparative metagenomics to be extended to dynamic profiling. However, with this method, the applicable target range is presently limited. Furthermore, the characteristics of coverage depth during replication have not been sufficiently investigated.

RESULTS

We developed a probabilistic model that mimics coverage depth dynamics. This statistical model explains the bias that occurs in the coverage depth due to DNA replication and errors that arise from coverage depth observation. Although our method requires a complete genome sequence, it involves a stable to low coverage depth (>0.01×). We also evaluated the estimation using real whole-genome sequence datasets and reproduced the growth dynamics observed in previous studies. By utilizing a circular distribution in the model, our method facilitates the quantification of unmeasured coverage depth features, including peakedness, skewness, and degree of density, around the replication origin. When we applied the model to time-series culture samples, the skewness parameter, which indicates the asymmetry, was stable over time; however, the peakedness and degree of density parameters, which indicate the concentration level at the replication origin, changed dynamically. Furthermore, we demonstrated the activity measurement of multiple replication origins in a single chromosome.

CONCLUSIONS

We devised a novel framework for quantifying coverage depth dynamics. Our study is expected to serve as a basis for replication activity estimation from a broader perspective using the statistical model.

Genomic Analysis of γ-Hexachlorocyclohexane-Degrading Sphingopyxis lindanitolerans WS5A3p Strain in the Context of the Pangenome of Sphingopyxis

Sphingopyxis inhabit diverse environmental niches, including marine, freshwater, oceans, soil and anthropogenic sites. The genus includes 20 phylogenetically distinct, valid species, but only a few with a sequenced genome. In this work, we analyzed the nearly complete genome of the newly described species, Sphingopyxislindanitolerans, and compared it to the other available Sphingopyxis genomes. The genome included 4.3 Mbp in total and consists of a circular chromosome, and two putative plasmids. Among the identified set of lin genes responsible for γ-hexachlorocyclohexane pesticide degradation, we discovered a gene coding for a new isoform of the LinA protein. The significant potential of this species in the remediation of contaminated soil is also correlated with the fact that its genome encodes a higher number of enzymes potentially involved in aromatic compound degradation than for most other Sphingopyxis strains. Additional analysis of 44 Sphingopyxis representatives provides insights into the pangenome of Sphingopyxis and revealed a core of 734 protein clusters and between four and 1667 unique proteins per genome.

Suppression of substrate inhibition in phenanthrene-degrading Mycobacterium by co-cultivation with a non-degrading Burkholderia strain

In natural environments contaminated by recalcitrant organic pollutants, efficient biodegradation of such pollutants has been suggested to occur through the cooperation of different bacterial species. A phenanthrene-degrading bacterial consortium, MixEPa4, from polluted soil was previously shown to include a phenanthrene-degrading strain, Mycobacterium sp. EPa45, and a non-polycyclic aromatic hydrocarbon (PAH)-degrading strain, Burkholderia sp. Bcrs1W. In this study, we show that addition of phenanthrene to rich liquid medium resulted in the transient growth arrest of EPa45 during its degradation of phenanthrene. RNA-sequencing analysis of the growth-arrested cells showed the phenanthrene-dependent induction of genes that were predicted to be involved in the catabolism of this compound, and many other cell systems, such as a ferric iron-uptake, were up-regulated, implying iron deficiency of the cells. This negative effect of phenanthrene became much more apparent when using phenanthrene-containing minimal agar medium; colony formation of EPa45 on such agar was significantly inhibited in the presence of phenanthrene and its intermediate degradation products. However, growth inhibition was suppressed by the co-residence of viable Bcrs1W cells. Various Gram-negative bacterial strains, including the three other strains from MixEPa4, also exhibited varying degrees of suppression of the growth inhibition effect on EPa45, strongly suggesting that this effect is not strain-specific. Growth inhibition of EPa45 was also observed by other PAHs, biphenyl and naphthalene, and these two compounds and phenanthrene also inhibited the growth of another mycobacterial strain, M. vanbaalenii PYR-1, that can use them as carbon sources. These phenomena of growth inhibition were also suppressed by Bcrs1W. Our findings suggest that, in natural environments, various non-PAH-degrading bacterial strains play potentially important roles in the facilitation of PAH degradation by the co-residing mycobacteria.

Biodegradation of Phenanthrene and Heavy Metal Removal by Acid-Tolerant Burkholderia fungorum FM-2

Phenanthrene (PHE) is a common pollutant of acidic and non-acidic environments that is recalcitrant to biodegradation. Herein, Burkholderia fungorum FM-2 (GenBank accession no. KM263605) was isolated from oil-contaminated soil in Xinjiang and characterized morphologically, physiologically, and phylogenetically. Environmental parameters including PHE concentration, pH, temperature, and salinity were optimized, and heavy metal tolerance was investigated. The MIC of strain FM-2 tolerant to Pb(II) and Cd(II) was 50 and 400 mg L⁻¹, respectively, while the MIC of Zn(II) was >1,200 mg L⁻¹. Atypically for a B. fungorum strain, FM-2 utilized PHE (300 mg L⁻¹) as a sole carbon source over a wide pH range (between pH 3 and 9). PHE and heavy metal metabolism were assessed using gas chromatography (GC), inductively coupled plasma optical emission spectroscopy (ICP-OES), scanning electron microscopy-energy dispersive X-ray spectroscopy (SEM-EDS), Fourier-transform infrared (FTIR) spectroscopy and ultraviolet (UV) absorption spectrometry. The effects of heavy metals on the bioremediation of PHE in soil were investigated, and the findings suggest that FM-2 has potential for combined bioremediation of soils co-contaminated with PHE and heavy metals.

The Relationship Between Microbial Community Structures and Environmental Parameters Revealed by Metagenomic Analysis of Hot Spring Water in the Kirishima Area, Japan

Diverse microorganisms specifically inhabit extreme environments, such as hot springs and deep-sea hydrothermal vents. To test the hypothesis that the microbial community structure is predictable based on environmental factors characteristic of such extreme environments, we conducted correlation analyses of microbial taxa/functions and environmental factors using metagenomic and 61 types of physicochemical data of water samples from nine hot springs in the Kirishima area (Kyusyu, Japan), where hot springs with diverse chemical properties are distributed in a relatively narrow area. Our metagenomic analysis revealed that the samples can be classified into two major types dominated by either phylum Crenarchaeota or phylum Aquificae. The correlation analysis showed that Crenarchaeota dominated in nutrient-rich environments with high concentrations of ions and total carbons, whereas Aquificae dominated in nutrient-poor environments with low ion concentrations. These environmental factors were also important explanatory variables in the generalized linear models constructed to predict the abundances of Crenarchaeota or Aquificae. Functional enrichment analysis of genes also revealed that the separation of the two major types is primarily attributable to genes involved in autotrophic carbon fixation, sulfate metabolism and nitrate reduction. Our results suggested that Aquificae and Crenarchaeota play a vital role in the Kirishima hot spring water ecosystem through their metabolic pathways adapted to each environment. Our findings provide a basis to predict microbial community structures in hot springs from environmental parameters, and also provide clues for the exploration of biological resources in extreme environments.

The dissipation and microbial ecotoxicity of tebuconazole and its transformation products in soil under standard laboratory and simulated winter conditions

Tebuconazole (TBZ) is a widely used triazole fungicide at EU level on cereals and vines. It is relatively persistent in soil where it is transformed to various transformation products (TPs) which might be environmentally relevant. We assessed the dissipation of TBZ in soil under contrasting incubation conditions (standard vs winter simulated) that are relevant to its application scheme, determined its transformation pathway using advanced analytical tools and ¹⁴C-labeled TBZ and assessed its soil microbial toxicity. Mineralization of ¹⁴C-triazole-ring-labeled TBZ was negligible but up to 11% of ¹⁴C-penyl-ring-labeled TBZ evolved as ¹⁴CO₂ within 150 days of incubation. TBZ persistence increased at higher dose rates (×10 compared to the recommended agronomical dose ×1) and under winter simulated conditions compared to standard incubation conditions (at ×1 dose rate DT₅₀ of 202 and 88 days, respectively). Non-target suspect screening enabled the detection of 22 TPs of TBZ, among which 17 were unknown. Mass spectrometry analysis led to the identification of 1-(4-chlorophenyl) ethanone, a novel TP of TBZ, the formation of which and decay in soil was determined by gas chromatography mass spectrometry. Three hypothetical transformation pathways of TBZ, all converging to 1H-1,2,4-triazole are proposed based on suspect screening. The ecotoxicological effect of TBZ and of its TPs was assessed by measuring by qPCR the abundance of the total bacteria and the relative abundance of 11 prokaryotic taxa and 4 functional groups. A transient impact of TBZ on the relative abundance of all prokaryotic taxa (except α-proteobacteria and Bacteroidetes) and one functional microbial group (pcaH-carrying microorganisms) was observed. However the direction of the effect (positive or negative) varied, and in certain cases, depended on the incubation conditions. Proteobacteria was the most responsive phylum to TBZ with recovery observed 20 days after treatment. The ecotoxicological effects on the soil microorganisms were not correlated with 1-(4-chlorophenyl) ethanone.

Latent environment allocation of microbial community data

As data for microbial community structures found in various environments has increased, studies have examined the relationship between environmental labels given to retrieved microbial samples and their community structures. However, because environments continuously change over time and space, mixed states of some environments and its effects on community formation should be considered, instead of evaluating effects of discrete environmental categories. Here we applied a hierarchical Bayesian model to paired datasets containing more than 30,000 samples of microbial community structures and sample description documents. From the training results, we extracted latent environmental topics that associate co-occurring microbes with co-occurring word sets among samples. Topics are the core elements of environmental mixtures and the visualization of topic-based samples clarifies the connections of various environments. Based on the model training results, we developed a web application, LEA (Latent Environment Allocation), which provides the way to evaluate typicality and heterogeneity of microbial communities in newly obtained samples without confining environmental categories to be compared. Because topics link words and microbes, LEA also enables to search samples semantically related to the query out of 30,000 microbiome samples.

In the past decade, microbiomes from various natural and human symbiotic environments have been thoroughly studied. However, our knowledge is limited as to what types of environments affect the structure of a microbial community. In the first place, how can we define “environments”, in particular, the environmental entities that are often continuously varying and difficult to discretely categorize? We assumed that environments could be represented from microbiome data because the structure of microbial communities reflect the state of the environment. We applied a probabilistic topic model to a dataset containing taxonomic composition data and natural language sample descriptions of >30,000 microbiome samples and extracted “latent environments” of the microbial communities, which are core elements of environmental mixtures. Integrating the training results of the model, we developed a web application to explore the microbiome universe and to place new metagenomic data on this universe like a global positioning system. Our tool shows what kinds of the environment naturally exist and are similar to each other on the perspective of the structural patterns of microbiome, and provides the way to evaluate the commonality and the heterogeneity of users’ microbiome samples.

Bioaugmentation of thiabendazole-contaminated soils from a wastewater disposal site: Factors driving the efficacy of this strategy and the diversity of the indigenous soil bacterial community

The application of the fungicide thiabendazole (TBZ) in fruit packaging plants (FPP) results in the production of effluents which are often disposed in adjacent field sites. These require remediation to prevent further environmental dispersal of TBZ. We assessed the bioaugmentation potential of a newly isolated TBZ-degrading bacterial consortium in a naturally contaminated soil (NCS) exhibiting a natural gradient of TBZ levels (12000, 400, 250 and 12 mg kg^-1). The effect of aging on bioaugmentation efficacy was comparatively tested in a soil with similar physicochemical properties and soil microbiota, which was artificially, contaminated with the same TBZ levels (ACS). The impact of bioaugmentation and TBZ on the bacterial diversity in the NCS was explored via amplicon sequencing. Bioaugmentation effectively removed TBZ from both soils at levels up to 400 mg kg^-1 but failed at the highest contamination level (12000 mg kg^-1). Dissipation of TBZ in bioaugmented samples showed a concentration-dependent pattern, while aging of TBZ had a slight effect on bioaugmentation efficiency. Bioaugmentation had no impact on the soil bacterial diversity, in contrast to TBZ contamination. Soils from the hotspots of TBZ contamination (12000 mg kg^-1) showed a drastically lower α-diversity driven by the dominance of β- and γ-proteobacteria at the expense of all other bacterial phyla, especially Actinobacteria. Overall, bioaugmentation with specialized microbial inocula could be an effective solution for the recovery of disposal sites contaminated with persistent chemicals like TBZ.

'TIME': A Web Application for Obtaining Insights into Microbial Ecology Using Longitudinal Microbiome Data

Realization of the importance of microbiome studies, coupled with the decreasing sequencing cost, has led to the exponential growth of microbiome data. A number of these microbiome studies have focused on understanding changes in the microbial community over time. Such longitudinal microbiome studies have the potential to offer unique insights pertaining to the microbial social networks as well as their responses to perturbations. In this communication, we introduce a web based framework called 'TIME' (Temporal Insights into Microbial Ecology'), developed specifically to obtain meaningful insights from microbiome time series data. The TIME web-server is designed to accept a wide range of popular formats as input with options to preprocess and filter the data. Multiple samples, defined by a series of longitudinal time points along with their metadata information, can be compared in order to interactively visualize the temporal variations. In addition to standard microbiome data analytics, the web server implements popular time series analysis methods like Dynamic time warping, Granger causality and Dickey Fuller test to generate interactive layouts for facilitating easy biological inferences. Apart from this, a new metric for comparing metagenomic time series data has been introduced to effectively visualize the similarities/differences in the trends of the resident microbial groups. Augmenting the visualizations with the stationarity information pertaining to the microbial groups is utilized to predict the microbial competition as well as community structure. Additionally, the 'causality graph analysis' module incorporated in TIME allows predicting taxa that might have a higher influence on community structure in different conditions. TIME also allows users to easily identify potential taxonomic markers from a longitudinal microbiome analysis. We illustrate the utility of the web-server features on a few published time series microbiome data and demonstrate the ease with which it can be used to perform complex analysis.

Differentiating current and past PCB and PCDD/F sources: The role of a large contaminated soil site in an industrialized city area

Cities and contaminated areas can be primary or secondary sources of polychlorinated biphenyls (PCBs), dibenzo-p-dioxins (PCDDs), dibenzofurans (PCDFs), and other chemicals, into air and soil and can influence the regional level of some of these pollutants. In a contaminated site, the evaluation of such emissions can be crucial in the choice of the remediation technology to be adopted. In the city of Brescia (Northern Italy), more than 100 ha of agricultural areas were contaminated with PCBs, PCDD/Fs and heavy metals, originating from the activities of a former PCB factory. In order to evaluate the current emissions of PCBs and PCDD/Fs from the contaminated site, in a location where other current sources are present, we compared measured and predicted air concentrations, resulting from chemical volatilization from soils as well as fingerprints of Brescia soils and of soils contaminated by specific sources. The results confirm that the contaminated area is still a current and important secondary source of PCBs to the air, and to a lesser extent of PCDFs (especially the more volatile), but not for PCDDs. PCBs in soils have fingerprints similar to highly chlorinated mixtures, indicating contamination by these mixtures and/or a long weathering process. PCB 209 is also present at important levels. PCDD fingerprints in soil cannot be related to current emission sources, while PCDFs are compatible to industrial and municipal waste incineration, although weathering and/or natural attenuation may have played a role in modifying such soil fingerprints. Finally, we combined chemical and microbiological analyses to provide an integrated approach to evaluate soil fingerprints and their variation in a wider perspective, which accounts for the mutual effects between contamination and soil microbiota, a pivotal hint for addressing in situ bioremediation activities.

Prokaryotic Virus Orthologous Groups (pVOGs): a resource for comparative genomics and protein family annotation

Viruses are the most abundant and diverse biological entities on earth, and while most of this diversity remains completely unexplored, advances in genome sequencing have provided unprecedented glimpses into the virosphere. The Prokaryotic Virus Orthologous Groups (pVOGs, formerly called Phage Orthologous Groups, POGs) resource has aided in this task over the past decade by using automated methods to keep pace with the rapid increase in genomic data. The uses of pVOGs include functional annotation of viral proteins, identification of genes and viruses in uncharacterized DNA samples, phylogenetic analysis, large-scale comparative genomics projects, and more. The pVOGs database represents a comprehensive set of orthologous gene families shared across multiple complete genomes of viruses that infect bacterial or archaeal hosts (viruses of eukaryotes will be added at a future date). The pVOGs are constructed within the Clusters of Orthologous Groups (COGs) framework that is widely used for orthology identification in prokaryotes. Since the previous release of the POGs, the size has tripled to nearly 3000 genomes and 300 000 proteins, and the number of conserved orthologous groups doubled to 9518. User-friendly webpages are available, including multiple sequence alignments and HMM profiles for each VOG. These changes provide major improvements to the pVOGs database, at a time of rapid advances in virus genomics. The pVOGs database is hosted jointly at the University of Iowa at http://dmk-brain.ecn.uiowa.edu/pVOGs and the NCBI at ftp://ftp.ncbi.nlm.nih.gov/pub/kristensen/pVOGs/home.html.

Reflection on Molecular Approaches Influencing State-of-the-Art Bioremediation Design: Culturing to Microbial Community Fingerprinting to Omics

Bioremediation is generally viewed as a cost effective and sustainable technology because it relies on microbes to transform pollutants into benign compounds. Advances in molecular biological analyses allow unprecedented microbial detection and are increasingly incorporated into bioremediation. Throughout history, state-of-the-art techniques have informed bioremediation strategies. However, the insights those techniques provided were not as in depth as those provided by recently developed omics tools. Advances in next generation sequencing (NGS) have now placed metagenomics and metatranscriptomics within reach of environmental engineers. As NGS costs decrease, metagenomics and metatranscriptomics have become increasingly feasible options to rapidly scan sites for specific degradative functions and identify microorganisms important in pollutant degradation. These omic techniques are capable of revolutionizing biological treatment in environmental engineering by allowing highly sensitive characterization of previously uncultured microorganisms. Omics enables the discovery of novel microorganisms for use in bioaugmentation and supports systematic optimization of biostimulation strategies. This review describes the omics journey from roots in biology and medicine to its current status in environmental engineering including potential future directions in commercial application.

Spatial and Temporal Dynamics of Pacific Oyster Hemolymph Microbiota across Multiple Scales

Unveiling the factors and processes that shape the dynamics of host associated microbial communities (microbiota) under natural conditions is an important part of understanding and predicting an organism's response to a changing environment. The microbiota is shaped by host (i.e., genetic) factors as well as by the biotic and abiotic environment. Studying natural variation of microbial community composition in multiple host genetic backgrounds across spatial as well as temporal scales represents a means to untangle this complex interplay. Here, we combined a spatially-stratified with a longitudinal sampling scheme within differentiated host genetic backgrounds by reciprocally transplanting Pacific oysters between two sites in the Wadden Sea (Sylt and Texel). To further differentiate contingent site from host genetic effects, we repeatedly sampled the same individuals over a summer season to examine structure, diversity and dynamics of individual hemolymph microbiota following experimental removal of resident microbiota by antibiotic treatment. While a large proportion of microbiome variation could be attributed to immediate environmental conditions, we observed persistent effects of antibiotic treatment and translocation suggesting that hemolymph microbial community dynamics is subject to within-microbiome interactions and host population specific factors. In addition, the analysis of spatial variation revealed that the within-site microenvironmental heterogeneity resulted in high small-scale variability, as opposed to large-scale (between-site) stability. Similarly, considerable within-individual temporal variability was in contrast with the overall temporal stability at the site level. Overall, our longitudinal, spatially-stratified sampling design revealed that variation in hemolymph microbiota is strongly influenced by site and immediate environmental conditions, whereas internal microbiome dynamics and oyster-related factors add to their long-term stability. The combination of small and large scale resolution of spatial and temporal observations therefore represents a crucial but underused tool to study host-associated microbiome dynamics.

Comparison of the complete genome sequences of four γ-hexachlorocyclohexane-degrading bacterial strains: insights into the evolution of bacteria able to degrade a recalcitrant man-made pesticide

γ-Hexachlorocyclohexane (γ-HCH) is a recalcitrant man-made chlorinated pesticide. Here, the complete genome sequences of four γ-HCH-degrading sphingomonad strains, which are most unlikely to have been derived from one ancestral γ-HCH degrader, were compared. Together with several experimental data, we showed that (i) all the four strains carry almost identical linA to linE genes for the conversion of γ-HCH to maleylacetate (designated “specific” lin genes), (ii) considerably different genes are used for the metabolism of maleylacetate in one of the four strains, and (iii) the linKLMN genes for the putative ABC transporter necessary for γ-HCH utilization exhibit structural divergence, which reflects the phylogenetic relationship of their hosts. Replicon organization and location of the lin genes in the four genomes are significantly different with one another, and that most of the specific lin genes are located on multiple sphingomonad-unique plasmids. Copies of IS6100, the most abundant insertion sequence in the four strains, are often located in close proximity to the specific lin genes. Analysis of the footprints of target duplication upon IS6100 transposition and the experimental detection of IS6100 transposition strongly suggested that the IS6100 transposition has caused dynamic genome rearrangements and the diversification of lin-flanking regions in the four strains.

Establishment of a multi-species biofilm model and metatranscriptomic analysis of biofilm and planktonic cell communities

We collected several biofilm samples from Japanese rivers and established a reproducible multi-species biofilm model that can be analyzed in laboratories. Bacterial abundance at the generic level was highly similar between the planktonic and biofilm communities, whereas comparative metatranscriptomic analysis revealed many upregulated and downregulated genes in the biofilm. Many genes involved in iron-sulfur metabolism, stress response, and cell envelope function were upregulated; biofilm formation is mediated by an iron-dependent signaling mechanism and the signal is relayed to stress-responsive and cell envelope function genes. Flagella-related gene expression was regulated depending upon the growth phase, indicating different roles of flagella during the adherence, maturation, and dispersal steps of biofilm formation. Downregulation of DNA repair genes was observed, indicating that spontaneous mutation frequency would be elevated within the biofilm and that the biofilm is a cradle for generating novel genetic traits. Although the significance remains unclear, genes for rRNA methyltransferase, chromosome partitioning, aminoacyl-tRNA synthase, and cysteine, methionine, leucine, thiamine, nucleotide, and fatty acid metabolism were found to be differentially regulated. These results indicate that planktonic and biofilm communities are in different dynamic states. Studies on biofilm and sessile cells, which have received less attention, are important for understanding microbial ecology and for designing tailor-made anti-biofilm drugs.

Complete Genome Sequence of a Phenanthrene Degrader, Burkholderia sp. HB-1 (NBRC 110738)

The phenanthrene-degrading Burkholderia sp. HB-1 was isolated from a phenanthrene-enrichment culture seeded with a pristine farm soil sample. We report the complete genome sequence of HB-1, which has been deposited to the stock culture (NBRC 110738) at Biological Resource Center, National Institute of Technology and Evaluation (NITE), Tokyo, Japan. The genome of strain HB-1 comprises two circular chromosomes of 4.1 Mb and 3.1 Mb. The finishing was facilitated by the computational tools GenoFinisher, AceFileViewer, and ShortReadManager.