Metagenomics reveals that detoxification systems are underrepresented in marine bacterial communities

Background

Environmental shotgun sequencing (metagenomics) provides a new way to study communities in microbial ecology. We here use sequence data from the Global Ocean Sampling (GOS) expedition to investigate toxicant selection pressures revealed by the presence of detoxification genes in marine bacteria. To capture a broad range of potential toxicants we selected detoxification protein families representing systems protecting microorganisms from a variety of stressors, such as metals, organic compounds, antibiotics and oxygen radicals.

Results

Using a bioinformatics procedure based on comparative analysis to finished bacterial genomes we found that the amount of detoxification genes present in marine microorganisms seems surprisingly small. The underrepresentation is particularly evident for toxicant transporters and proteins involved in detoxifying metals. Exceptions are enzymes involved in oxidative stress defense where peroxidase enzymes are more abundant in marine bacteria compared to bacteria in general. In contrast, catalases are almost completely absent from the open ocean environment, suggesting that peroxidases and peroxiredoxins constitute a core line of defense against reactive oxygen species (ROS) in the marine milieu.

Conclusions

We found no indication that detoxification systems would be generally more abundant close to the coast compared to the open ocean. On the contrary, for several of the protein families that displayed a significant geographical distribution, like peroxidase, penicillin binding transpeptidase and divalent ion transport protein, the open ocean samples showed the highest abundance. Along the same lines, the abundance of most detoxification proteins did not increase with estimated pollution. The low level of detoxification systems in marine bacteria indicate that the majority of marine bacteria have a low capacity to adapt to increased pollution. Our study exemplifies the use of metagenomics data in ecotoxicology, and in particular how anthropogenic consequences on life in the sea can be examined.

Electronic supplementary material

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

Sequence heterogeneity of TT virus and closely related viruses

TT virus (TTV) is a recently discovered infectious agent originally obtained from transfusion-related hepatitis. However, the causative link between the TTV infection and liver disease remains uncertain. Recent studies demonstrated that genome sequences of different TTV strains are significantly divergent. To assess genetic heterogeneity of the TTV genome in more detail, a sequence analysis of PCR fragments (271 bp) amplified from open reading frame 1 (ORF1) was performed. PCR fragments were amplified from 5 to 40% of serum specimens obtained from patients with different forms of hepatitis who reside in different countries (e.g., China, Egypt, Vietnam, and the United States) and from normal human specimens obtained from U.S. residents. A total of 170 PCR fragments were sequenced and compared to sequences derived from the corresponding TTV genome region deposited in GenBank. Genotypes 2 and 3 were found to be significantly more genetically related than any other TTV genotype. Moreover, three sequences were shown to be almost equally related to both genotypes 2 and 3. These observations suggest a merger of genotypes 2 and 3 into one genotype, 2/3. Additionally, five new groups of TTV sequences were identified. One group represents a new genotype, whereas the other four groups were shown to be more evolutionary distant from all known TTV sequences. The evolutionary distances between these four groups were also shown to be greater than between TTV genotypes. The phylogenetic analysis suggested that these four new genetic groups represent closely related yet different viral species. Thus, TTV exists as a "swarm" of at least five closely related but different viruses. These observations suggest a high degree of genetic complexity within the TTV population. The finding of the additional TTV-related species should be taken into consideration when the association between TTV infections and human diseases of unknown etiology is studied.     

Genetic diversity in marine algal virus communities as revealed by sequence analysis of DNA polymerase genes.

Algal-virus-specific PCR primers were used to amplify DNA polymerase gene (pol) fragments (683 to 689 bp) from the virus-sized fraction (0.02 to 0.2 microns) concentrated from inshore and offshore water samples collected from the Gulf of Mexico. Algal-virus-like DNA pol genes were detected in five samples collected from the surface and deep chlorophyll maximum. PCR products from an offshore station were cloned, and the genetic diversity of 33 fragments was examined by restriction fragment length polymorphism and sequence analysis. The five different genotypes or operational taxonomic units (OTUs) that were identified on the basis of restriction fragment length polymorphism banding patterns were present in different relative abundances (9 to 34%). One clone from each OTU was sequenced, and phylogenetic analysis showed that all of the OTUs fell within the family Phycodnaviridae. Four of the OTUs fell within a group of viruses (MpV) which infect the photosynthetic picoplankter Micromonas pusilla. The genetic diversity among these genotypes was as large as that previously found for MpV isolates from different oceans. The remaining genotype formed its own clade between viruses which infect M. pusilla and Chrysochromulina brevifilum. These results imply that marine virus communities contain a diverse assemblage of MpV-like viruses, as well as other unknown members of the Phycodnaviridae.     

Evidence from GC-TRFLP that bacterial communities in soil are lognormally distributed

The Species Abundance Distribution (SAD) is a fundamental property of ecological communities and the form and formation of SADs have been examined for a wide range of communities including those of microorganisms. Progress in understanding microbial SADs, however, has been limited by the remarkable diversity and vast size of microbial communities. As a result, few microbial systems have been sampled with sufficient depth to generate reliable estimates of the community SAD. We have used a novel approach to characterize the SAD of bacterial communities by coupling genomic DNA fractionation with analysis of terminal restriction fragment length polymorphisms (GC-TRFLP). Examination of a soil microbial community through GC-TRFLP revealed 731 bacterial operational taxonomic units (OTUs) that followed a lognormal distribution. To recover the same 731 OTUs through analysis of DNA sequence data is estimated to require analysis of 86,264 16S rRNA sequences. The approach is examined and validated through construction and analysis of simulated microbial communities in silico. Additional simulations performed to assess the potential effects of PCR bias show that biased amplification can cause a community whose distribution follows a power-law function to appear lognormally distributed. We also show that TRFLP analysis, in contrast to GC-TRFLP, is not able to effectively distinguish between competing SAD models. Our analysis supports use of the lognormal as the null distribution for studying the SAD of bacterial communities as for plant and animal communities.     

Nearly identical bacteriophage structural gene sequences are widely distributed in both marine and freshwater environments

Primers were designed to amplify a 592-bp region within a conserved structural gene (g20) found in some cyanophages. The goal was to use this gene as a proxy to infer genetic richness in natural cyanophage communities and to determine if sequences were more similar in similar environments. Gene products were amplified from samples from the Gulf of Mexico, the Arctic, Southern, and Northeast and Southeast Pacific Oceans, an Arctic cyanobacterial mat, a catfish production pond, lakes in Canada and Germany, and a depth of ca. 3,246 m in the Chuckchi Sea. Amplicons were separated by denaturing gradient gel electrophoresis, and selected bands were sequenced. Phylogenetic analysis revealed four previously unknown groups of g20 clusters, two of which were entirely found in freshwater. Also, sequences with >99% identities were recovered from environments that differed greatly in temperature and salinity. For example, nearly identical sequences were recovered from the Gulf of Mexico, the Southern Pacific Ocean, an Arctic freshwater cyanobacterial mat, and Lake Constance, Germany. These results imply that closely related hosts and the viruses infecting them are distributed widely across environments or that horizontal gene exchange occurs among phage communities from very different environments. Moreover, the amplification of g20 products from deep in the cyanobacterium-sparse Chuckchi Sea suggests that this primer set targets bacteriophages other than those infecting cyanobacteria.     

Spatial analysis of ectomycorrhizal fungi reveals that root tip communities are structured by competitive interactions

Microbial ecology has made large advances over the last decade, mostly because of improvements in molecular analysis techniques that have enabled the detection and identification of progressively larger numbers of microbial species. However, determining the ecological patterns and processes taking place in communities of microbes remains a significant challenge. Are communities randomly assembled through dispersal and priority effects, or do species interact with each other leading to positive and negative associations? For mycorrhizal fungi, evidence is accumulating that stochastic and competitive interactions between species may both have a role in shaping community structure. Could the methodological approach, which is often incidence based, impact the outcomes detected? Here, we applied an incidence‐based Terminal Restriction Fragment Length Polymorphism (T‐RFLP) database approach to examine species diversity and ecological interactions within a community of ectomycorrhizal (ECM) fungi. Co‐occurrence analysis revealed that the ECM community colonizing root tips was strongly structured by competitive interactions, or ecological processes generating a similar spatial pattern, rather than neutral processes. Analysis of β‐diversity indicated that community structure was significantly more similar (spatially autocorrelated) at distances equal to or <3.41 m. The eight most frequently encountered species in the root tip community of ECM fungi displayed significant competitive interactions with at least one other species, showing that the incidence‐based approach was capable of detecting this sort of ecological information.     

Clades of Adeno-associated viruses are widely disseminated in human tissues

The potential for using Adeno-associated virus (AAV) as a vector for human gene therapy has stimulated interest in the Dependovirus genus. Serologic data suggest that AAV infections are prevalent in humans, although analyses of viruses and viral sequences from clinical samples are extremely limited. Molecular techniques were used in this study to successfully detect endogenous AAV sequences in 18% of all human tissues screened, with the liver and bone marrow being the most predominant sites. Sequence characterization of rescued AAV DNAs indicated a diverse array of molecular forms which segregate into clades whose members share functional and serologic similarities. One of the most predominant human clades is a hybrid of two previously described AAV serotypes, while another clade was found in humans and several species of nonhuman primates, suggesting a cross-species transmission of this virus. These data provide important information regarding the biology of parvoviruses in humans and their use as gene therapy vectors.     

Assemblage of pollinator communities in four widely isolated nature reserves of southern Punjab,Pakistan

Pollinators are key components of the global biodiversity, providing vital ecosystem services of pollination to crops and wild plants. Understanding pollinator assemblage and how it changes spatially is important in order to device effective ecosystem management planning. This is the first empirical study to explore pollinator community structure in four widely isolated nature reserves of southern Punjab, Pakistan, i.e., the central agricultural plain (Pirowal Forest), the Indus delta plain (Ghazi Ghat Wetland), hilly mountains (Fort Munro), and a sandy desert (Lal Suhanra Forest). Six assemblage parameters, i.e., abundance, richness, dominance, evenness, Simpson index, and Shannon-Wiener index, were studied over seven consecutive months (February to August 2012). Maximum abundance and richness of the pollinators were recorded at Ghazi Ghat (1354 individuals and 84 species), followed by Pirowal (1331 individuals and 80 species), Lal Suhanra (1197 individuals and 77 species), and Fort Munro (808 individuals and 65 species). The highest values of Simpson and evenness indexes were recorded at Fort Munro, while the Shannon-Wiener and dominance indexes were highest at Ghazi Ghat and Pirowal. Species richness was highest for bees, while it was lowest for flies, yet the latter were the highest in abundance among the four pollinator groups. The current study is a first account of pollinator assemblage structure in four widely isolated forest reserves of Punjab, Pakistan. This will support ecosystem management planning, and opens doors for further research in order to reveal complexities in pollinator assemblage structure especially in relation with plant communities.     

Viruses occurring in East Africa that are related to peanut mottle virus

Viruses occurring in Cassia bicapsularis in Northern Tanzania, in Voandzeia subterranea in north western and eastern Tanzania, and in Phaseolus lunatus in the Kenya highlands, were all serologically related to peanut mottle virus. Their host ranges, and the symptoms they induced in test plants, were very similar, and they differed only in degree of virulence in some host species. The Voandzeia isolate did not infect groundnut, and only the Phaseolus isolate infected two species in the Cucurbitaceae. All the isolates infected Chenopodium amaranticolor, a species which formerly was reported as being immune to peanut mottle and thus considered of diagnostic value. In Africa, variation in peanut mottle virus isolates seems to be associated with host species and ecology, and there is at present no evidence for naturally occurring variants within a host species as occurs in groundnut in America. Three of the four isolates were purified by homogenising together infected leaf tissue, chloroform and 0.5 M sodium citrate buffer containing 1% 2-mercapto-ethanol at pH 8, in the proportion 1: 1:2 respectively, and precipitating the virus from the clarified homogenate with 5% w/v polyethylene glycol. When centrifuged in sucrose density gradients such preparations gave a single, bright specific light scat-tering zone with no haze.     

Two paralogous families of a two-gene subtilisin operon are widely distributed in oral treponemes

Certain oral treponemes express a highly proteolytic phenotype and have been associated with periodontal diseases. The periodontal pathogen Treponema denticola produces dentilisin, a serine protease of the subtilisin family. The two-gene operon prcA-prtP is required for expression of active dentilisin (PrtP), a putative lipoprotein attached to the treponeme's outer membrane or sheath. The purpose of this study was to examine the diversity and structure of treponemal subtilisin-like proteases in order to better understand their distribution and function. The complete sequences of five prcA-prtP operons were determined for Treponema lecithinolyticum, "Treponema vincentii," and two canine species. Partial operon sequences were obtained for T. socranskii subsp. 04 as well as 450- to 1,000-base fragments of prtP genes from four additional treponeme strains. Phylogenetic analysis demonstrated that the sequences fall into two paralogous families. The first family includes the sequence from T. denticola. Treponemes possessing this operon family express chymotrypsin-like protease activity and can cleave the substrate N-succinyl-alanyl-alanyl-prolyl-phenylalanine-p-nitroanilide (SAAPFNA). Treponemes possessing the second paralog family do not possess chymotrypsin-like activity or cleave SAAPFNA. Despite examination of a range of protein and peptide substrates, the specificity of the second protease family remains unknown. Each of the fully sequenced prcA and prtP genes contains a 5' hydrophobic leader sequence with a treponeme lipobox. The two paralogous families of treponeme subtilisins represent a new subgroup within the subtilisin family of proteases and are the only subtilisin lipoprotein family. The present study demonstrated that the subtilisin paralogs comprising a two-gene operon are widely distributed among treponemes.     

Sequence analysis reveals new membrane anchor of reaction centre-bound cytochromes possibly related to PufX

Most of the bacterial photosynthetic reaction centres known to date contain a cytochrome subunit with four covalently bound haem groups. In the case of Blastochloris viridis, this reaction centre subunit is anchored in the membrane by a lipid molecule covalently attached to the cysteine which forms the N-terminus of the mature protein after processing by a signal peptidase. We show that posttranslational N-terminal cleavage of the cytochrome subunit does not occur in the aerobic photosynthetic bacterium Roseobacter denitrificans. From sequence analysis of the resulting elongated N-terminus it follows that a transmembrane helix is anchoring the reaction centre-bound cytochrome in the membrane. Comparative sequence analysis strongly suggests that all cytochrome subunits lacking the lipid coupling cysteine share this structural feature. Comparison of the N-terminal segment of the cytochrome subunit of Roseobacter denitrificans with the sequences of the PufX proteins from Rhodobacter sphaeroides and Rhodobacter capsulatus suggests a phylogenetic relation.     

Statistical analysis of vertebrate sequences reveals that long genes are scarce in GC-rich isochores

We compared the exon/intron organization of vertebrate genes belonging to different isochore classes, as predicted by their GC content at third codon position. Two main features have emerged from the analysis of sequences published in GenBank: (1) genes coding for long proteins (i.e., 500 aa) are almost two times more frequent in GC-poor than in GC-rich isochores; (2) intervening sequences (=sum of introns) are on average three times longer in GC-poor than in GC-rich isochores. These patterns are observed among human, mouse, rat, cow, and even chicken genes and are therefore likely to be common to all warm-blooded vertebrates. Analysis of Xenopus sequences suggests that the same patterns exist in cold-blooded vertebrates. It could be argued that such results do not reflect the reality because sequence databases are not representative of entire genomes. However, analysis of biases in GenBank revealed that the observed discrepancies between GC-rich and GC-poor isochores are not artifactual, and are probably largely underestimated. We investigated the distribution of microsatellites and interspersed repeats in introns of human and mouse genes from different isochores. This analysis confirmed previous studies showing that Ll repeats are almost absent from GC-rich isochores. Microsatellites and SINES (Alu, B1, B2) are found at roughly equal frequencies in introns from all isochore classes. Globally, the presence of repeated sequences does not account for the increased intron length in GC-poor isochores. The relationships between gene structure and global genome organization and evolution are discussed.     

Functional analysis of a DNA-shuffled movement protein reveals that microtubules are dispensable for the cell-to-cell movement of tobacco mosaic virus

Microtubules interact strongly with the viral movement protein (MP) of Tobacco mosaic virus (TMV) and are thought to transport the viral genome between plant cells. We describe a functionally enhanced DNA-shuffled movement protein (MP(R3)) that remained bound to the vertices of the cortical endoplasmic reticulum, showing limited affinity for microtubules. A single amino acid change was shown to confer the MP(R3) phenotype. Disruption of the microtubule cytoskeleton in situ with pharmacological agents, or by silencing of the alpha-tubulin gene, had no significant effect on the spread of TMV vectors expressing wild-type MP (MP(WT)) and did not prevent the accumulation of MP(WT) in plasmodesmata. Thus, cell-to-cell trafficking of TMV can occur independently of microtubules. The MP(R3) phenotype was reproduced when infection sites expressing MP(WT) were treated with a specific proteasome inhibitor, indicating that the degradation of MP(R3) is impaired. We suggest that the improved viral transport functions of MP(R3) arise from evasion of a host degradation pathway.     

Geographical features are the predominant driver of molecular diversification in widely distributed North American whipsnakes

Allopatric divergence following the formation of geographical features has been implicated as a major driver of evolutionary diversification. Widespread species complexes provide opportunities to examine allopatric divergence across varying degrees of isolation in both time and space. In North America, several geographical features may play such a role in diversification, including the Mississippi River, Pecos River, Rocky Mountains, Cochise Filter Barrier, Gulf of California and Isthmus of Tehuantepec. We used thousands of nuclear single nucleotide polymorphisms (SNPs) and mitochondrial DNA from several species of whipsnakes (genera Masticophis and Coluber) distributed across North and Central America to investigate the role that these geographical features have played on lineage divergence. We hypothesize that these features restrict gene flow and separate whipsnakes into diagnosable genomic clusters. We performed genomic clustering and phylogenetic reconstructions at the species and population levels using Bayesian and likelihood analyses and quantified migration levels across geographical features to assess the degree of genetic isolation due to allopatry. Our analyses suggest that (i) major genetic divisions are often consistent with isolation by geographical features, (ii) migration rates between clusters are asymmetrical across major geographical features, and (iii) areas that receive proportionally more migrants possess higher levels of genetic diversity. Collectively, our findings suggest that multiple features of the North American landscape contributed to allopatric divergence in this widely distributed snake group.     

Phosphoenolpyruvate:glycose phosphotransferase system in species of Vibrio, a widely distributed marine bacterial genus.

The genus Vibrio is one of the most common and widely distributed groups of marine bacteria. Studies on the physiology of marine Vibrio species were initiated by examining 15 species for the bacterial phosphoenolpyruvate:glycose phosphotransferase system (PTS). All species tested contained a PTS analogous to the glucose-specific (IIGlc) system in enteric bacteria. Crude extracts of the cells showed immunological cross-reactivity with antibodies to enzyme I, HPr, and IIIGlc from Salmonella typhimurium when assayed by the rocket-line method. Toluene-permeabilized cells of 11 species were tested and were active in phosphorylating methyl alpha-D-glucoside with phosphoenolpyruvate but not ATP as the phosphoryl donor. Membranes from 10 species were assayed, and they phosphorylated methyl alpha-D-glucoside when supplemented with a phospho-IIIGlc-generating system composed of homogeneous proteins from enteric bacteria. Toluene-permeabilized cells and membranes of seven species were assayed, as were phosphorylated fructose and 2-deoxyglucose. IIIGlc was isolated from Vibrio fluvialis and was active in phosphorylating methyl alpha-D-glucoside when supplemented with a phospho-HPr-generating system composed of homogeneous proteins from Escherichia coli and membranes from either E. coli or V. fluvialis. These results show that the bacterial PTS is widely distributed in the marine environment and that it is likely to have a significant role in marine bacterial physiology and in the marine ecosystem.     

Metabolic profiling reveals that time related physiological changes in mammalian cell perfusion cultures are bioreactor scale independent

Metabolic profiling was used to characterize the time course of cell physiology both in laboratory- and manufacturing-scale mammalian cell perfusion cultures. Two independent experiments were performed involving three vials from the same BHK cell bank, used to inoculate three laboratory-scale bioreactors, from which four manufacturing-scale cultures were initiated. It was shown that metabolomic analysis can indeed enhance the prime variable dataset for the monitoring of perfusion cultures by providing a higher resolution view of the metabolic state. Metabolic profiles could capture physiological state shifts over the course of the perfusion cultures and indicated a metabolic “signature” of the phase transitions, which was not observable from prime variable data. Specifically, the vast majority of metabolites had lower concentrations in the middle compared to the other two phases. Notably, metabolomics provided orthogonal (to prime variables) evidence that all cultures followed this same metabolic state shift with cell age, independently of bioreactor scale.