Statistical Approach of Functional Profiling for a Microbial Community

Background

Metagenomics is a relatively new but fast growing field within environmental biology and medical sciences. It enables researchers to understand the diversity of microbes, their functions, cooperation, and evolution in a particular ecosystem. Traditional methods in genomics and microbiology are not efficient in capturing the structure of the microbial community in an environment. Nowadays, high-throughput next-generation sequencing technologies are powerfully driving the metagenomic studies. However, there is an urgent need to develop efficient statistical methods and computational algorithms to rapidly analyze the massive metagenomic short sequencing data and to accurately detect the features/functions present in the microbial community. Although several issues about functions of metagenomes at pathways or subsystems level have been investigated, there is a lack of studies focusing on functional analysis at a low level of a hierarchical functional tree, such as SEED subsystem tree.

Results

A two-step statistical procedure (metaFunction) is proposed to detect all possible functional roles at the low level from a metagenomic sample/community. In the first step a statistical mixture model is proposed at the base of gene codons to estimate the abundances for the candidate functional roles, with sequencing error being considered. As a gene could be involved in multiple biological processes the functional assignment is therefore adjusted by utilizing an error distribution in the second step. The performance of the proposed procedure is evaluated through comprehensive simulation studies. Compared with other existing methods in metagenomic functional analysis the new approach is more accurate in assigning reads to functional roles, and therefore at more general levels. The method is also employed to analyze two real data sets.

Conclusions

metaFunction is a powerful tool in accurate profiling functions in a metagenomic sample.     

Genetic Diversity of Microbial Eukaryotes in Anoxic Sediment of the Saline Meromictic Lake Namako-ike (Japan): On the Detection of Anaerobic or Anoxic-tolerant Lineages of Eukaryotes

Available sequence data on eukaryotic small-subunit ribosomal DNA (SSU rDNA) directly retrieved from various environments have increased recently, and the diversity of microbial eukaryotes (protists) has been shown to be much greater than previously expected. However, the molecular information accumulated to date does still not thoroughly reveal ecological distribution patterns of microbial eukaryotes. In the ongoing challenge to detect anaerobic or anoxic-tolerant lineages of eukaryotes, we directly extracted DNA from the anoxic sediment of a saline meromictic lake, constructed genetic libraries of PCR-amplified SSU rDNA, and performed phylogenetic analyses with the cloned SSU rDNA sequences. Although a few sequences could not be confidently assigned to any major eukaryotic groups in the analyses and are debatable regarding their taxonomic positions, most sequences obtained have affiliations with known major lineages of eukaryotes (Cercozoa, Alveolata, Stramenopiles, and Opisthokonta). Among these sequences, some branched with lineages predominantly composed of uncultured environmental clones retrieved from other anoxic environments, while others were closely related to those of eukaryotic parasites (e.g. Phytomyxea of Cercozoa, Gregarinea of Alveolata, and Ichthyosporea of Opisthokonta).     

High Species Richness of Scinax Treefrogs (Hylidae) in a Threatened Amazonian Landscape Revealed by an Integrative Approach

Rising habitat loss is one of the main drivers of the global amphibian decline. Nevertheless, knowledge of amphibian diversity needed for effective habitat protection is still highly inadequate in remote tropical regions, the greater part of the Amazonia. In this study we integrated molecular, morphological and bioacoustic evidence to evaluate the species richness of the treefrogs genus Scinax over a 1000 km transect across rainforest of the Purus-Madeira interfluve, and along the east bank of the upper Madeira river, Brazilian Amazonia. Analysis revealed that 82% of the regional species richness of Scinax is still undescribed; two nominal species, seven confirmed candidate species, two unconfirmed candidate species, and one deep conspecific lineage were detected in the study area. DNA barcoding based analysis of the 16s rRNA gene indicates possible existence of three discrete species groups within the genus Scinax, in addition to the already-known S. rostratus species Group. Quantifying and characterizing the number of undescribed Scinax taxa on a regional scale, we provide a framework for future taxonomic study in Amazonia. These findings indicate that the level to which Amazonian anura species richness has been underestimated is far greater than expected. Consequently, special attention should be paid both to taxonomic studies and protection of the still-neglected Amazonian Scinax treefrogs.     

Interspecific Comparison of the Transformer Gene of Drosophila Reveals an Unusually High Degree of Evolutionary Divergence

The transformer (tra) gene of Drosophila melanogaster occupies an intermediate position in the regulatory pathway controlling all aspects of somatic sexual differentiation. The female-specific expression of this gene's function is regulated by the Sex lethal (Sxl) gene, through a mechanism involving sex-specific alternative splicing of tra pre-mRNA. The tra gene encodes a protein that is thought to act in conjunction with the transformer-2 (tra-2) gene product to control the sex-specific processing of doublesex (dsx) pre-mRNA. The bifunctional dsx gene carries out opposite functions in the two sexes, repressing female differentiation in males and repressing male differentiation in females. Here we report the results from an evolutionary approach to investigate tra regulation and function, by isolating the tra-homologous genes from selected Drosophila species, and then using the interspecific DNA sequence comparisons to help identify regions of functional significance. The tra-homologous genes from two Sophophoran subgenus species, Drosophila simulans and Drosophila erecta, and two Drosophila subgenus species, Drosophila hydei and Drosophila virilis, were cloned, sequenced and compared to the D. melanogaster tra gene. This comparison reveals an unusually high degree of evolutionary divergence among the tra coding sequences. These studies also highlight a highly conserved sequence within intron one that probably defines a cis-acting regulator of the sex-specific alternative splicing event.     

Microbial Formation of Ethane in Anoxic Estuarine Sediments

Estuarine sediment slurries produced methane and traces of ethane when incubated under hydrogen. Formation of methane occurred over a broad temperature range with an optimum above 65°C. Ethane formation had a temperature optimum at 40°C. Formation of these two gases was inhibited by air, autoclaving, incubation at 4 and 80°C, and by the methanogenic inhibitor, 2-bromoethanesulfonic acid. Ethane production was stimulated by addition of ethylthioethanesulfonic acid, and production from ethylthioethanesulfonic acid was blocked by 2-bromoethanesulfonic acid. A highly purified enrichment culture of a methanogenic bacterium obtained from sediments produced traces of ethane from ethylthioethanesulfonic acid. These results indicate that the small quantities of ethane found in anaerobic sediments can be formed by certain methanogenic bacteria.     

Microbial Formation of Dimethyl Sulfide in Anoxic Sphagnum Peat

Peat bogs dominated by Sphagnum spp. have relatively high areal rates of dimethyl sulfide (DMS) emission to the atmosphere. DMS was produced in anoxic slurries of Sphagnum peat with a linear time course and with an average rate of 40.4 (range, 22.0 to 68.6) nmol per liter of slurry (middot) day(sup-1) observed in nine batches of slurry. Methanethiol (MeSH) was produced at roughly similar rates over the typical 4- to 8-day incubations. DMS and MeSH production in these acidic (pH 4.2 to 4.6) peats were biological, as they were stopped completely by autoclaving and inhibited strongly by addition of antibiotics and 500 (mu)M chloroform. Endogenous DMS production may be due to the degradation of S-methyl-methionine, dimethyl sulfoxide, or methoxyaromatic compounds (e.g., syringic acid), each of which stimulated DMS formation when added at 5 to 10 (mu)M concentrations. However, on the basis of the high rates of thiol (MeSH and ethanethiol) methylation activity that we observed and the availability of endogenous MeSH, we suggest that methylation of MeSH is the major pathway leading to DMS formation in anaerobic peat. Solid-phase adsorption of MeSH plays a key role in its availability for biomethylation reactions. Additions of acetate (1.5 mM) or compounds which could cause acetate to accumulate (e.g., glucose, alanine, and 2-bromoethanesulfonate) suppressed DMS formation. It is likely that acetogenic bacteria are involved in DMS formation, but our data are insufficient to allow firm conclusions about the metabolic pathways or organisms involved. Our observations are the first which point to the methylation of MeSH as the major mechanism for endogenous DMS production in any environment. The rates of net DMS production observed are sufficient to explain the relatively high fluxes of DMS emitted to the atmosphere from Sphagnum sp.-dominated wetlands.     

Validation and Genotyping of Multiple Human Polymorphic Inversions Mediated by Inverted Repeats Reveals a High Degree of Recurrence

In recent years different types of structural variants (SVs) have been discovered in the human genome and their functional impact has become increasingly clear. Inversions, however, are poorly characterized and more difficult to study, especially those mediated by inverted repeats or segmental duplications. Here, we describe the results of a simple and fast inverse PCR (iPCR) protocol for high-throughput genotyping of a wide variety of inversions using a small amount of DNA. In particular, we analyzed 22 inversions predicted in humans ranging from 5.1 kb to 226 kb and mediated by inverted repeat sequences of 1.6–24 kb. First, we validated 17 of the 22 inversions in a panel of nine HapMap individuals from different populations, and we genotyped them in 68 additional individuals of European origin, with correct genetic transmission in ∼12 mother-father-child trios. Global inversion minor allele frequency varied between 1% and 49% and inversion genotypes were consistent with Hardy-Weinberg equilibrium. By analyzing the nucleotide variation and the haplotypes in these regions, we found that only four inversions have linked tag-SNPs and that in many cases there are multiple shared SNPs between standard and inverted chromosomes, suggesting an unexpected high degree of inversion recurrence during human evolution. iPCR was also used to check 16 of these inversions in four chimpanzees and two gorillas, and 10 showed both orientations either within or between species, providing additional support for their multiple origin. Finally, we have identified several inversions that include genes in the inverted or breakpoint regions, and at least one disrupts a potential coding gene. Thus, these results represent a significant advance in our understanding of inversion polymorphism in human populations and challenge the common view of a single origin of inversions, with important implications for inversion analysis in SNP-based studies.     

Canopy Flow Analysis Reveals the Advantage of Size in the Oldest Communities of Multicellular Eukaryotes

Download : Download video (22MB)     

Crystal Structure of Human Seryl-tRNA Synthetase and Ser-SA Complex Reveals a Molecular Lever Specific to Higher Eukaryotes

1. Download : Download high-res image (197KB)
2. Download : Download full-size image

     

Microbial Dynamics of an Epilithic Mat Community in a High Alpine Stream

Studies were conducted to examine interrelationships between the heterotrophic and phototrophic populations within an epilithic community in the outlet stream of a high alpine lake. Levels of nitrates, phosphates, and total organic compounds in the lake were consistently near the lower limits of detectability. Microscopic examination of the community by phase-contrast light microscopy and scanning electron microscopy revealed diatoms, filamentous algae, and bacteria embedded within a dense gelatinous matrix. Chlorophyll a and primary productivity measurements had peak values in early August, with subsequent declines. Bacterial heterotrophic activity, as measured by V_max, turnover rate, and relative activity, increased significantly as the phototrophic community declined. This trend in heterotrophic activity was not accompanied by an increase in total bacterial numbers as determined by epi-illuminated fluorescence microscopy. These results suggest that the phototrophic community responded to changes in, or interactions among, various chemical and physical factors throughout the study period. The catabolic activity of the sessile bacteria appeared to be positively influenced by changes in the mat environment resulting from the decline of the phototrophic populations.     

Electrochemical Investigation of a Microbial Solar Cell Reveals a Nonphotosynthetic Biocathode Catalyst

Microbial solar cells (MSCs) are microbial fuel cells (MFCs) that generate their own oxidant and/or fuel through photosynthetic reactions. Here, we present electrochemical analyses and biofilm 16S rRNA gene profiling of biocathodes of sediment/seawater-based MSCs inoculated from the biocathode of a previously described sediment/seawater-based MSC. Electrochemical analyses indicate that for these second-generation MSC biocathodes, catalytic activity diminishes over time if illumination is provided during growth, whereas it remains relatively stable if growth occurs in the dark. For both illuminated and dark MSC biocathodes, cyclic voltammetry reveals a catalytic-current–potential dependency consistent with heterogeneous electron transfer mediated by an insoluble microbial redox cofactor, which was conserved following enrichment of the dark MSC biocathode using a three-electrode configuration. 16S rRNA gene profiling showed Gammaproteobacteria, most closely related to Marinobacter spp., predominated in the enriched biocathode. The enriched biocathode biofilm is easily cultured on graphite cathodes, forms a multimicrobe-thick biofilm (up to 8.2 μm), and does not lose catalytic activity after exchanges of the reactor medium. Moreover, the consortium can be grown on cathodes with only inorganic carbon provided as the carbon source, which may be exploited for proposed bioelectrochemical systems for electrosynthesis of organic carbon from carbon dioxide. These results support a scheme where two distinct communities of organisms develop within MSC biocathodes: one that is photosynthetically active and one that catalyzes reduction of O₂ by the cathode, where the former partially inhibits the latter. The relationship between the two communities must be further explored to fully realize the potential for MSC applications.     

Distribution and Diversity of Microbial Eukaryotes in Bathypelagic Waters of the South China Sea

Little is known about the biodiversity of microbial eukaryotes in the South China Sea, especially in waters at bathyal depths. Here, we employed SSU rDNA gene sequencing to reveal the diversity and community structure across depth and distance gradients in the South China Sea. Vertically, the highest alpha diversity was found at 75‐m depth. The communities of microbial eukaryotes were clustered into shallow‐, middle‐, and deep‐water groups according to the depth from which they were collected, indicating a depth‐related diversity and distribution pattern. Rhizaria sequences dominated the microeukaryote community and occurred in all samples except those from less than 50‐m deep, being most abundant near the sea floor where they contributed ca. 64–97% and 40–74% of the total sequences and OTUs recovered, respectively. A large portion of rhizarian OTUs has neither a nearest named neighbor nor a nearest neighbor in the GenBank database which indicated the presence of new phylotypes in the South China Sea. Given their overwhelming abundance and richness, further phylogenetic analysis of rhizarians were performed and three new genetic clusters were revealed containing sequences retrieved from the deep waters of the South China Sea. Our results shed light on the diversity and community structure of microbial eukaryotes in this not yet fully explored area.     

Diversity and Distribution of Marine Microbial Eukaryotes in the Arctic Ocean and Adjacent Seas

We analyzed microbial eukaryote diversity in perennially cold arctic marine waters by using 18S rRNA gene clone libraries. Samples were collected during concurrent oceanographic missions to opposite sides of the Arctic Ocean Basin and encompassed five distinct water masses. Two deep water Arctic Ocean sites and the convergence of the Greenland, Norwegian, and Barents Seas were sampled from 28 August to 2 September 2002. An additional sample was obtained from the Beaufort Sea (Canada) in early October 2002. The ribotypes were diverse, with different communities among sites and between the upper mixed layer and just below the halocline. Eukaryotes from the remote Canada Basin contained new phylotypes belonging to the radiolarian orders Acantharea, Polycystinea, and Taxopodida. A novel group within the photosynthetic stramenopiles was also identified. One sample closest to the interior of the Canada Basin yielded only four major taxa, and all but two of the sequences recovered belonged to the polar diatom Fragilariopsis and a radiolarian. Overall, 42% of the sequences were <98% similar to any sequences in GenBank. Moreover, 15% of these were <95% similar to previously recovered sequences, which is indicative of endemic or undersampled taxa in the North Polar environment. The cold, stable Arctic Ocean is a threatened environment, and climate change could result in significant loss of global microbial biodiversity.     

Molecular Analyses of the Microbial Community Composition of an Anoxic Basin of a Municipal Wastewater Treatment Plant Reveal a Novel Lineage of Proteobacteria

A culture-independent molecular phylogenetic approach was used to study prokaryotic diversity in an anoxic activated sludge from a municipal wastewater treatment plant. Two 16S rRNA gene libraries were constructed using total genomic DNA and amplified by polymerase chain reaction using primers specific for archaeal or bacterial domains. Phylogenetic analysis of 132 and 249 almost full-length 16S rRNA genes for Archaea and Bacteria, respectively, was done using the ARB software package. Phylogenetic groups affiliated with the Archaea belong to Euryarchaeota (93.8% of the operational taxonomic units [OTUs]) and Crenarchaeota (6.2% of the OTUs). Within the bacterial library, 84.8% of the OTUs represent novel putative phylotypes never described before and affiliated with ten divisions. The Proteobacteria phylum is the most abundant and diversified phylogenetic group representing 60.4% of the OTUs, followed by Bacteroidetes (22.1%) and gram-positives (6.1%). Interestingly, we detected a novel Proteobacteria monophyletic group distinct from the five known subclasses, which we named New Lineage of Proteobacteria (NLP) lineage, and it is composed of eight clones representing 4.6% of the Proteobacteria. A new 16S rRNA-targeted hybridization probe was designed and fluorescent in situ hybridization analyses shows representatives of NLP as cocci-shaped microorganisms. The Chloroflexi, Acidobacterium, and Nitrospira phyla and TM7 candidate division are each represented by ≤3% of clone sequences. A comprehensive set of eight 16S and 23S rRNA-targeted oligonucleotide probes was used to quantify these major groups by dot blot hybridization within 12 samples. The Proteobacteria accounted for 82.5 ± 4.9%, representing the most abundant phyla. The Bacteroidetes and Planctomycetales groups accounted for 4.9 ± 1.3% and 4 ± 1.7%, respectively. Firmicutes and Actinobacteria together accounted for only 1.9 ± 0.5%. The set of probes covers 93.4 ± 14% of the total bacterial population rRNA within the anoxic basin.     

Taxonomic Identity in Microbial Eukaryotes: A Practical Approach Using the Testate Amoeba Centropyxis to Resolve Conflicts Between Old and New Taxonomic Descriptions

ABSTRACT. The present work focuses on 12 taxa of the genus Centropyxis Stein, 1857 to explore the conflict between traditional and contemporary taxonomic practices. We examined the morphology, biometry, and ecology of 2,120 Centropyxis individuals collected from Tiete River, Sao Paulo, Brazil; with these new data we studied the consistency of previously described species, varieties, and forms. We encountered transitional forms of test morphology that undermine specific and varietal distinctions for three species and nine varieties. Biometrical analyses made comparing the organisms at the species level suggest a lack of separation between Centropyxis aculeata and Centropyxis discoides, and a possible distinction for Centropyxis ecornis based on spine characteristics. However, incongruence between recent and previous surveys makes taking any taxonomic–nomenclatural actions inadvisable, as they would only add to the confusion. We suggest an explicit and objective taxonomic practice in order to enhance our taxonomic and species concepts for microbial eukaryotes. This will allow more precise inferences of taxon identity for studies in other areas.     

Strong Seasonality of Marine Microbial Eukaryotes in a High-Arctic Fjord (Isfjorden,in West Spitsbergen,Norway)

The Adventfjorden time series station (IsA) in Isfjorden, West Spitsbergen, Norway, was sampled frequently from December 2011 to December 2012. The community composition of microbial eukaryotes (size, 0.45 to 10 μm) from a depth of 25 m was determined using 454 sequencing of the 18S V4 region amplified from both DNA and RNA. The compositional changes throughout the year were assessed in relation to in situ fjord environmental conditions. Size fractionation analyses of chlorophyll a showed that the photosynthetic biomass was dominated by small cells (<10 μm) most of the year but that larger cells dominated during the spring and summer. The winter and early-spring communities were more diverse than the spring and summer/autumn communities. Dinophyceae were predominant throughout the year. The Arctic Micromonas ecotype was abundant mostly in the early-bloom and fall periods, whereas heterotrophs, such as marine stramenopiles (MASTs), Picozoa, and the parasitoid marine alveolates (MALVs), displayed higher relative abundance in the winter than in other seasons. Our results emphasize the extreme seasonality of Arctic microbial eukaryotic communities driven by the light regime and nutrient availability but point to the necessity of a thorough knowledge of hydrography for full understanding of their succession and variability.     

Comparative Genomic Study Reveals a Transition from TA Richness in Invertebrates to GC Richness in Vertebrates at CpG Flanking Sites： An Indication for Context-Dependent Mutagenicity of Methylated CpG Sites

Vertebrate genomes are characterized with CpG deficiency, particularly for GCpoor regions. The GC content-related CpG deficiency is probably caused by context-dependent deamination of methylated CpG sites. This hypothesis was examined in this study by comparing nucleotide frequencies at CpG flanking positions among invertebrate and vertebrate genomes. The finding is a transition of nucleotide preference of 5＇ T to 5＇ A at the invertebrate-vertebrate boundary, indicating that a large number of CpG sites with 5＇ Ts were depleted because of global DNA methylation developed in vertebrates. At genome level, we investigated CpG observed/expected （obs/exp） values in 500 bp fragments, and found that higher CpG obs/exp value is shown in GC-poor regions of invertebrate genomes （except sea urchin） but in GC-rich sequences of vertebrate genomes. We next compared GC content at CpG flanking positions with genomic average, showing that the GC content is lower than the average in invertebrate genomes, but higher than that in vertebrate genomes. These results indicate that although 5＇ T and 5＇ A are different in inducing deamination of methylated CpG sites, GC content is even more important in affecting the deamination rate. In all the tests, the results of sea urchin are similar to vertebrates perhaps due to its fractional DNA methylation. CpG deficiency is therefore suggested to be mainly a result of high mutation rates of methylated CpG sites in GC-poor regions.     

Relationship between Degree of Dominance and Species Richness in Grass Communities with Different Productivities

The aim of this study is to test the assumption that the relationship between the degree of dominance and local species richness may be different in grass communities with different productivities. Alpine, subalpine, and low-mountain grasslands, as well as subalpine mires, alpine communities of low-snow habitats and those with long-term snow cover, steppe communities, and the grass layer of low-mountain forest communities of the Western Caucasus and Ciscaucasia, are used as objects of research. The data on the phytomass of 419 plots with an area of 0.25 m² are studied. The results show that, the higher the mean productivity of communities is, the closer the relationship between the degree of dominance and species richness is, and the closest relationship is observed in meadow communities. Possible causes of these relationships are considered. It is reasonably suggested that this may be due to the features of the organization of plant communities with high and low productivity (in particular, high or low intensity of interspecific competition).