Bats, Bacteria, and Bat Smell V.2.0: Repeatable Sex-Specific Differences in Scent Organ Microbiota

Reproducibility is a fundamental principle in science, ensuring reliable and valid findings. However, replication studies are scarce, particularly in ecology, due to the emphasis on novelty for publication. We explored the possibility of replicating original findings in the field of microbial and chemical ecology by conducting a conceptual replication of a previous study analysing the sex-specific differences in the microbial communities inhabiting the wing sacs, a scent organ with crucial functions in olfactory communication, of greater sac-winged bat (Saccopteryx bilineata). In the original study, the skin swabs from the antebrachial wing sacs of the males and wing sac rudiments of the females were analysed using culture-dependent methods to test sex-specific differences. The authors demonstrated that males have lower microbial richness and different microbial composition than females. We attempted to reproduce these findings using 16S rRNA sequencing, which offers improved accuracy in pinpointing microbial members than culture-dependent methods because of advanced statistical methods. Our study validated the original study's findings: Males had a lower microbial richness, and the community composition differed between the sexes. Furthermore, in the current study, males had an increased abundance of bacteria that might potentially be involved in odour production and degradation of malodorous substances and antimicrobial production. Our conceptual replication study corroborated that microbes can play a role in shaping their host's olfactory phenotype and consequently influence sexual selection. Furthermore, the current study emphasises the importance of replication efforts and hopefully encourages a culture that values replication studies in scientific practice.

Correction: GABenchToB: A Genome Assembly Benchmark Tuned on Bacteria and Benchtop Sequencers

[This corrects the article DOI: 10.1371/journal.pone.0107014.].

Timing matters: age-dependent impacts of the social environment and host selection on the avian gut microbiota

BACKGROUND

The establishment of the gut microbiota in early life is a critical process that influences the development and fitness of vertebrates. However, the relative influence of transmission from the early social environment and host selection throughout host ontogeny remains understudied, particularly in avian species. We conducted conspecific and heterospecific cross-fostering experiments in zebra finches (Taeniopygia guttata) and Bengalese finches (Lonchura striata domestica) under controlled conditions and repeatedly sampled the faecal microbiota of these birds over the first 3 months of life. We thus documented the development of the gut microbiota and characterised the relative impacts of the early social environment and host selection due to species-specific characteristics and individual genetic backgrounds across ontogeny by using 16S ribosomal RNA gene sequencing.

RESULTS

The taxonomic composition and community structure of the gut microbiota changed across ontogenetic stages; juvenile zebra finches exhibited higher alpha diversity than adults at the post-breeding stage. Furthermore, in early development, the microbial communities of juveniles raised by conspecific and heterospecific foster parents resembled those of their foster family, emphasising the importance of the social environment. In later stages, the social environment continued to influence the gut microbiota, but host selection increased in importance.

CONCLUSIONS

We provided a baseline description of the developmental succession of gut microbiota in zebra finches and Bengalese finches, which is a necessary first step for understanding the impact of the early gut microbiota on host fitness. Furthermore, for the first time in avian species, we showed that the relative strengths of the two forces that shape the establishment and maintenance of the gut microbiota (i.e. host selection and dispersal from the social environment) change during development, with host selection increasing in importance. This finding should be considered when experimentally manipulating the early-life gut microbiota. Our findings also provide new insights into the mechanisms of host selection. Video Abstract.

Changes to the gut microbiota of a wild juvenile passerine in a multidimensional urban mosaic

Urbanisation is a major anthropogenic perturbation presenting novel ecological and evolutionary challenges to wild populations. Symbiotic microorganisms residing in the gastrointestinal tracts (gut) of vertebrates have mutual connections with host physiology and respond quickly to environmental alterations. However, the impact of anthropogenic changes and urbanisation on the gut microbiota remains poorly understood, especially in early development. To address this knowledge gap, we characterised the gut microbiota of juvenile great tits (Parus major) reared in artificial nestboxes and in natural cavities in an urban mosaic, employing two distinct frameworks characterising the urban space. Microbial diversity was influenced by cavity type. Alpha diversity was affected by the amount of impervious surface surrounding the breeding location, and positively correlated with tree cover density. Community composition differed between urban and rural sites: these alterations covaried with sound pollution and distance to the city centre. Overall, the microbial communities reflect and are possibly influenced by the heterogeneous environmental modifications that are typical of the urban space. Strikingly, the choice of framework and environmental variables characterising the urban space can influence the outcomes of such ecological studies. Our results open new perspectives to investigate the impact of microbial symbionts on the adaptive capacity of their hosts.

Multiple Occurrences of a 168-Nucleotide Deletion in SARS-CoV-2 ORF8, Unnoticed by Standard Amplicon Sequencing and Variant Calling Pipelines

Genomic surveillance of the SARS-CoV-2 pandemic is crucial and mainly achieved by amplicon sequencing protocols. Overlapping tiled-amplicons are generated to establish contiguous SARS-CoV-2 genome sequences, which enable the precise resolution of infection chains and outbreaks. We investigated a SARS-CoV-2 outbreak in a local hospital and used nanopore sequencing with a modified ARTIC protocol employing 1200 bp long amplicons. We detected a long deletion of 168 nucleotides in the ORF8 gene in 76 samples from the hospital outbreak. This deletion is difficult to identify with the classical amplicon sequencing procedures since it removes two amplicon primer-binding sites. We analyzed public SARS-CoV-2 sequences and sequencing read data from ENA and identified the same deletion in over 100 genomes belonging to different lineages of SARS-CoV-2, pointing to a mutation hotspot or to positive selection. In almost all cases, the deletion was not represented in the virus genome sequence after consensus building. Additionally, further database searches point to other deletions in the ORF8 coding region that have never been reported by the standard data analysis pipelines. These findings and the fact that ORF8 is especially prone to deletions, make a clear case for the urgent necessity of public availability of the raw data for this and other large deletions that might change the physiology of the virus towards endemism.

The Gut Microbial Composition Is Species-Specific and Individual-Specific in Two Species of Estrildid Finches, the Bengalese Finch and the Zebra Finch

Microbial communities residing in the gastrointestinal tracts of animals have profound impacts on the physiological processes of their hosts. In humans, host-specific and environmental factors likely interact together to shape gut microbial communities, resulting in remarkable inter-individual differences. However, we still lack a full understanding of to what extent microbes are individual-specific and controlled by host-specific factors across different animal taxa. Here, we document the gut microbial characteristics in two estrildid finch species, the Bengalese finch (Lonchura striata domestica) and the zebra finch (Taeniopygia guttata) to investigate between-species and within-species differences. We collected fecal samples from breeding pairs that were housed under strictly controlled environmental and dietary conditions. All individuals were sampled at five different time points over a range of 120 days covering different stages of the reproductive cycle. We found significant species-specific differences in gut microbial assemblages. Over a period of 3 months, individuals exhibited unique, individual-specific microbial profiles. Although we found a strong individual signature in both sexes, within-individual variation in microbial communities was larger in males of both species. Furthermore, breeding pairs had more similar microbial profiles, compared to randomly chosen males and females. Our study conclusively shows that host-specific factors contribute structuring of gut microbiota.

Family matters: skin microbiome reflects the social group and spatial proximity in wild zebra finches

Background

So far, large numbers of studies investigating the microbiome have focused on gut microbiota and less have addressed the microbiome of the skin. Especially in avian taxa our understanding of the ecology and function of these bacteria remains incomplete. The involvement of skin bacteria in intra-specific communication has recently received attention, and has highlighted the need to understand what information is potentially being encoded in bacterial communities. Using next generation sequencing techniques, we characterised the skin microbiome of wild zebra finches, aiming to understand the impact of sex, age and group composition on skin bacteria communities. For this purpose, we sampled skin swabs from both sexes and two age classes (adults and nestlings) of 12 different zebra finch families and analysed the bacterial communities.

Results

Using 16S rRNA sequencing we found no effect of age, sex and family on bacterial diversity (alpha diversity). However, when comparing the composition (beta diversity), we found that animals of social groups (families) harbour highly similar bacterial communities on their skin with respect to community composition. Within families, closely related individuals shared significantly more bacterial taxa than non-related animals. In addition, we found that age (adults vs. nestlings) affected bacterial composition. Finally, we found that spatial proximity of nest sites, and therefore individuals, correlated with the skin microbiota similarity.

Conclusions

Birds harbour very diverse and complex bacterial assemblages on their skin. These bacterial communities are distinguishable and characteristic for intraspecific social groups. Our findings are indicative for a family-specific skin microbiome in wild zebra finches. Genetics and the (social) environment seem to be the influential factors shaping the complex bacterial communities. Bacterial communities associated with the skin have a potential to emit volatiles and therefore these communities may play a role in intraspecific social communication, e.g. via signalling social group membership.

Importance of Defluviitalea raffinosedens for Hydrolytic Biomass Degradation in Co-Culture with Hungateiclostridium thermocellum

Bacterial hydrolysis of polysaccharides is an important step for the production of sustainable energy, for example during the conversion of plant biomass to methane-rich biogas. Previously, Hungateiclostridium thermocellum was identified as cellulolytic key player in thermophilic biogas microbiomes with a great frequency as an accompanying organism. The aim of this study was to physiologically characterize a recently isolated co-culture of H. thermocellum and the saccharolytic bacterium Defluviitalea raffinosedens from a laboratory-scale biogas fermenter. The characterization focused on cellulose breakdown by applying the measurement of cellulose hydrolysis, production of metabolites, and the activity of secreted enzymes. Substrate degradation and the production of volatile metabolites was considerably enhanced when both organisms acted synergistically. The metabolic properties of H. thermocellum have been studied well in the past. To predict the role of D. raffinosedens in this bacterial duet, the genome of D. raffinosedens was sequenced for the first time. Concomitantly, to deduce the prevalence of D. raffinosedens in anaerobic digestion, taxonomic composition and transcriptional activity of different biogas microbiomes were analyzed in detail. Defluviitalea was abundant and metabolically active in reactor operating at highly efficient process conditions, supporting the importance of this organism for the hydrolysis of the raw substrate.

Individual- and Species-Specific Skin Microbiomes in Three Different Estrildid Finch Species Revealed by 16S Amplicon Sequencing

An animals' body is densely populated with bacteria. Although a large number of investigations on physiological microbial colonisation have emerged in recent years, our understanding of the composition, ecology and function of the microbiota remains incomplete. Here, we investigated whether songbirds have an individual-specific skin microbiome that is similar across different body regions. We collected skin microbe samples from three different bird species (Taeniopygia gutatta, Lonchura striata domestica and Stagonopleura gutatta) at two body locations (neck region, preen gland area). To characterise the skin microbes and compare the bacterial composition, we used high-throughput 16S rRNA amplicon sequencing. This method proved suitable for identifying the skin microbiome of birds, even though the bacterial load on the skin appeared to be relatively low. We found that across all species, the two evaluated skin areas of each individual harboured very similar microbial communities, indicative of an individual-specific skin microbiome. Despite experiencing the same environmental conditions and consuming the same diet, significant differences in the skin microbe composition were identified among the three species. The bird species differed both quantitatively and qualitatively regarding the observed bacterial taxa. Although each species harboured its own unique set of skin microbes, we identified a core skin microbiome among the studied species. As microbes are known to influence the host's body odour, our findings of an individual-specific skin microbiome might suggest that the skin microbiome in birds is involved in the odour production and could encode information on the host's genotype.

Do we treat our patients or rather periodontal microbes with adjunctive antibiotics in periodontal therapy? A 16S rDNA microbial community analysis

Empiric antibiotics are often used in combination with mechanical debridement to treat patients suffering from periodontitis and to eliminate disease-associated pathogens. Until now, only a few next generation sequencing 16S rDNA amplicon based publications with rather small sample sizes studied the effect of those interventions on the subgingival microbiome. Therefore, we studied subgingival samples of 89 patients with chronic periodontitis (solely non-smokers) before and two months after therapy. Forty-seven patients received mechanical periodontal therapy only, whereas 42 patients additionally received oral administered amoxicillin plus metronidazole (500 and 400 mg, respectively; 3x/day for 7 days). Samples were sequenced with Illumina MiSeq 300 base pairs paired end technology (V3 and V4 hypervariable regions of the 16S rDNA). Inter-group differences before and after therapy of clinical variables (percentage of sites with pocket depth ≥ 5mm, percentage of sites with bleeding on probing) and microbiome variables (diversity, richness, evenness, and dissimilarity) were calculated, a principal coordinate analysis (PCoA) was conducted, and differential abundance of agglomerated ribosomal sequence variants (aRSVs) classified on genus level was calculated using a negative binomial regression model. We found statistically noticeable decreased richness, and increased dissimilarity in the antibiotic, but not in the placebo group after therapy. The PCoA revealed a clear compositional separation of microbiomes after therapy in the antibiotic group, which could not be seen in the group receiving mechanical therapy only. This difference was even more pronounced on aRSV level. Here, adjunctive antibiotics were able to induce a microbiome shift by statistically noticeably reducing aRSVs belonging to genera containing disease-associated species, e.g., Porphyromonas, Tannerella, Treponema, and Aggregatibacter, and by noticeably increasing genera containing health-associated species. Mechanical therapy alone did not statistically noticeably affect any disease-associated taxa. Despite the difference in microbiome modulation both therapies improved the tested clinical parameters after two months. These results cast doubt on the relevance of the elimination and/or reduction of disease-associated taxa as a main goal of periodontal therapy.

Bioinformatics for NGS-based metagenomics and the application to biogas research

Metagenomics has proven to be one of the most important research fields for microbial ecology during the last decade. Starting from 16S rRNA marker gene analysis for the characterization of community compositions to whole metagenome shotgun sequencing which additionally allows for functional analysis, metagenomics has been applied in a wide spectrum of research areas. The cost reduction paired with the increase in the amount of data due to the advent of next-generation sequencing led to a rapidly growing demand for bioinformatic software in metagenomics. By now, a large number of tools that can be used to analyze metagenomic datasets has been developed. The Bielefeld-Gießen center for microbial bioinformatics as part of the German Network for Bioinformatics Infrastructure bundles and imparts expert knowledge in the analysis of metagenomic datasets, especially in research on microbial communities involved in anaerobic digestion residing in biogas reactors. In this review, we give an overview of the field of metagenomics, introduce into important bioinformatic tools and possible workflows, accompanied by application examples of biogas surveys successfully conducted at the Center for Biotechnology of Bielefeld University.

GABenchToB: a genome assembly benchmark tuned on bacteria and benchtop sequencers

De novo genome assembly is the process of reconstructing a complete genomic sequence from countless small sequencing reads. Due to the complexity of this task, numerous genome assemblers have been developed to cope with different requirements and the different kinds of data provided by sequencers within the fast evolving field of next-generation sequencing technologies. In particular, the recently introduced generation of benchtop sequencers, like Illumina's MiSeq and Ion Torrent's Personal Genome Machine (PGM), popularized the easy, fast, and cheap sequencing of bacterial organisms to a broad range of academic and clinical institutions. With a strong pragmatic focus, here, we give a novel insight into the line of assembly evaluation surveys as we benchmark popular de novo genome assemblers based on bacterial data generated by benchtop sequencers. Therefore, single-library assemblies were generated, assembled, and compared to each other by metrics describing assembly contiguity and accuracy, and also by practice-oriented criteria as for instance computing time. In addition, we extensively analyzed the effect of the depth of coverage on the genome assemblies within reasonable ranges and the k-mer optimization problem of de Bruijn Graph assemblers. Our results show that, although both MiSeq and PGM allow for good genome assemblies, they require different approaches. They not only pair with different assembler types, but also affect assemblies differently regarding the depth of coverage where oversampling can become problematic. Assemblies vary greatly with respect to contiguity and accuracy but also by the requirement on the computing power. Consequently, no assembler can be rated best for all preconditions. Instead, the given kind of data, the demands on assembly quality, and the available computing infrastructure determines which assembler suits best. The data sets, scripts and all additional information needed to replicate our results are freely available at ftp://ftp.cebitec.uni-bielefeld.de/pub/GABenchToB.

Metagenome analyses reveal the influence of the inoculant Lactobacillus buchneri CD034 on the microbial community involved in grass ensiling

Silage is green fodder conserved by lactic acid fermentation performed by epiphytic lactic acid bacteria under anaerobic conditions. To improve the ensiling process and the quality of the resulting silage, starter cultures are added to the fresh forage. A detailed analysis of the microbial community playing a role in grass ensiling has been carried out by high throughput sequencing technologies. Moreover, the influence of the inoculant Lactobacillus buchneri CD034 on the microbial community composition was studied. For this purpose, grass was ensiled untreated or inoculated with L. buchneri CD034. The fresh forage as well as silages after 14 and 58 days of fermentation were characterized physico-chemically. Characteristic silage conditions such as increased titers of lactic acid bacteria and higher concentrations of acetic acid were observed in the inoculated silage in comparison to the untreated samples. Taxonomic community profiles deduced from 16S rDNA amplicon sequences indicated that the relative abundance of Lactococci diminished in the course of fermentations and that the proportion of bacteria belonging to the phyla Proteobacteria and Bacteroidetes increased during the fermentation of untreated silage. In the inoculated silage, members of these phyla were repressed due to an increased abundance of Lactobacilli. In addition, metagenome analyses of silage samples confirmed taxonomic profiles based on 16S rDNA amplicons. Moreover, Lactobacillus plantarum, Lactobacillus brevis and Lactococcus lactis were found to be dominant species within silages as analyzed by means of fragment recruitments of metagenomic sequence reads on complete reference genome sequences. Fragment recruitments also provided clear evidence for the competitiveness of the inoculant strain L. buchneri CD034 during the fermentation of the inoculated silage. The inoculation strain was able to outcompete other community members and also affected physico-chemical characteristics of the silage.

Changes in root bacterial communities associated to two different development stages of canola (Brassica napus L. var oleifera) evaluated through next-generation sequencing technology

Crop production may benefit from plant growth-promoting bacteria. The knowledge on bacterial communities is indispensable in agricultural systems that intend to apply beneficial bacteria to improve plant health and production of crops such as canola. In this work, the diversity of root bacterial communities associated to two different developmental phases of canola (Brassica napus L.) plants was assessed through the application of new generation sequencing technology. Total bacterial DNA was extracted from root samples from two different growth states of canola (rosette and flowering). It could be shown how bacterial communities inside the roots changed with the growing stage of the canola plants. There were differences in the abundance of the genera, family, and even the phyla identified for each sample. While in both root samples Proteobacteria was the most common phylum, at the rosette stage, the most common bacteria belonged to the family Pseudomonadaceae and the genus Pseudomonas, and in the flowering stage, the Xanthomonadaceae family and the genus Xanthomonas dominated the community. This implies in a switch in the predominant bacteria in the different developmental stages of the plant, suggesting that the plant itself interferes with the associated microbial community.

Bacterial community shift in treated periodontitis patients revealed by ion torrent 16S rRNA gene amplicon sequencing

Periodontitis, one of the most common diseases in the world, is caused by a mixture of pathogenic bacteria and inflammatory host responses and often treated by antimicrobials as an adjunct to scaling and root planing (SRP). Our study aims to elucidate explorative and descriptive temporal shifts in bacterial communities between patients treated by SRP alone versus SRP plus antibiotics. This is the first metagenomic study using an Ion Torrent Personal Genome Machine (PGM). Eight subgingival plaque samples from four patients with chronic periodontitis, taken before and two months after intervention were analyzed. Amplicons from the V6 hypervariable region of the 16S rRNA gene were generated and sequenced each on a 314 chip. Sequencing reads were clustered into operational taxonomic units (OTUs, 3% distance), described by community metrics, and taxonomically classified. Reads ranging from 599,933 to 650,416 per sample were clustered into 1,648 to 2,659 non-singleton OTUs, respectively. Increased diversity (Shannon and Simpson) in all samples after therapy was observed regardless of the treatment type whereas richness (ACE) showed no correlation. Taxonomic analysis revealed different microbial shifts between both therapy approaches at all taxonomic levels. Most remarkably, the genera Porphyromonas, Tannerella, Treponema, and Filifactor all harboring periodontal pathogenic species were removed almost only in the group treated with SPR and antibiotics. For the species T. forsythia and P. gingivalis results were corroborated by real-time PCR analysis. In the future, hypothesis free metagenomic analysis could be the key in understanding polymicrobial diseases and be used for therapy monitoring. Therefore, as read length continues to increase and cost to decrease, rapid benchtop sequencers like the PGM might finally be used in routine diagnostic.

Importance of rare taxa for bacterial diversity in the rhizosphere of Bt- and conventional maize varieties

Ribosomal 16S rRNA gene pyrosequencing was used to explore whether the genetically modified (GM) Bt-maize hybrid MON 89034 × MON 88017, expressing three insecticidal recombinant Cry proteins of Bacillus thuringiensis, would alter the rhizosphere bacterial community. Fine roots of field cultivated Bt-maize and three conventional maize varieties were analyzed together with coarse roots of the Bt-maize. A total of 547,000 sequences were obtained. Library coverage was 100% at the phylum and 99.8% at the genus rank. Although cluster analyses based on relative abundances indicated no differences at higher taxonomic ranks, genera abundances pointed to variety specific differences. Genera-based clustering depended solely on the 49 most dominant genera while the remaining 461 rare genera followed a different selection. A total of 91 genera responded significantly to the different root environments. As a benefit of pyrosequencing, 79 responsive genera were identified that might have been overlooked with conventional cloning sequencing approaches owing to their rareness. There was no indication of bacterial alterations in the rhizosphere of the Bt-maize beyond differences found between conventional varieties. B. thuringiensis-like phylotypes were present at low abundance (0.1% of Bacteria) suggesting possible occurrence of natural Cry proteins in the rhizospheres. Although some genera indicated potential phytopathogenic bacteria in the rhizosphere, their abundances were not significantly different between conventional varieties and Bt-maize. With an unprecedented sensitivity this study indicates that the rhizosphere bacterial community of a GM maize did not respond abnormally to the presence of three insecticidal proteins in the root tissue.

Profiling of the metabolically active community from a production-scale biogas plant by means of high-throughput metatranscriptome sequencing

Structural composition and gene content of a biogas-producing microbial community from a production-scale biogas plant fed with renewable primary products was recently analyzed by means of a metagenome sequencing approach. To determine the transcriptionally active part of the same biogas community and to identify key transcripts for the biogas production process, the metatranscriptome of the microorganisms was sequenced for the first time. The metatranscriptome sequence dataset generated on the Genome Sequencer FLX platform is represented by 484,920 sequence reads. Taxonomic profiling of the active part of the community by classification of 16S ribosomal sequence tags revealed that members of the Euryarchaeota and Firmicutes account for the dominant phyla. Only smaller fractions of the 16S ribosomal sequence tags were assigned to the phyla Bacteroidetes, Actinobacteria and Synergistetes. Among the mRNA-derived sequence tags from the metatranscriptome dataset, transcripts encoding enzymes involved in substrate hydrolysis, acidogenesis, acetate formation and methanogenesis could be identified. Transcripts for enzymes functioning in methanogenesis are among the most abundant mRNA tags indicating that the corresponding pathway is very active in the methanogenic sub-community. As a frame of reference for evaluation of metatranscriptome sequence data, the 16S rDNA-based taxonomic profile of the community was analyzed by means of high-throughput 16S rDNA amplicon sequencing. Processing of the obtained amplicon reads resulted in 18,598 high-quality 16S rDNA sequences covering the V3-V4 hypervariable region of the 16S rRNA gene. Comparison of the taxonomic profiles deduced from 16S rDNA amplicon sequences and the metatranscriptome dataset indicates a high transcriptional activity of archaeal species. Overall, it was shown that the most abundant species dominating the community also contributed the majority of the transcripts. In the future, key transcripts for the biogas production process will provide valuable markers for evaluation of the performance of biogas-producing microbial communities with the objective to optimize the biotechnology of this process.

Comparative and joint analysis of two metagenomic datasets from a biogas fermenter obtained by 454-pyrosequencing

Biogas production from renewable resources is attracting increased attention as an alternative energy source due to the limited availability of traditional fossil fuels. Many countries are promoting the use of alternative energy sources for sustainable energy production. In this study, a metagenome from a production-scale biogas fermenter was analysed employing Roche's GS FLX Titanium technology and compared to a previous dataset obtained from the same community DNA sample that was sequenced on the GS FLX platform. Taxonomic profiling based on 16S rRNA-specific sequences and an Environmental Gene Tag (EGT) analysis employing CARMA demonstrated that both approaches benefit from the longer read lengths obtained on the Titanium platform. Results confirmed Clostridia as the most prevalent taxonomic class, whereas species of the order Methanomicrobiales are dominant among methanogenic Archaea. However, the analyses also identified additional taxa that were missed by the previous study, including members of the genera Streptococcus, Acetivibrio, Garciella, Tissierella, and Gelria, which might also play a role in the fermentation process leading to the formation of methane. Taking advantage of the CARMA feature to correlate taxonomic information of sequences with their assigned functions, it appeared that Firmicutes, followed by Bacteroidetes and Proteobacteria, dominate within the functional context of polysaccharide degradation whereas Methanomicrobiales represent the most abundant taxonomic group responsible for methane production. Clostridia is the most important class involved in the reductive CoA pathway (Wood-Ljungdahl pathway) that is characteristic for acetogenesis. Based on binning of 16S rRNA-specific sequences allocated to the dominant genus Methanoculleus, it could be shown that this genus is represented by several different species. Phylogenetic analysis of these sequences placed them in close proximity to the hydrogenotrophic methanogen Methanoculleus bourgensis. While rarefaction analyses still indicate incomplete coverage, examination of the GS FLX Titanium dataset resulted in the identification of additional genera and functional elements, providing a far more complete coverage of the community involved in anaerobic fermentative pathways leading to methane formation.

WebCARMA: a web application for the functional and taxonomic classification of unassembled metagenomic reads

Background

Metagenomics is a new field of research on natural microbial communities. High-throughput sequencing techniques like 454 or Solexa-Illumina promise new possibilities as they are able to produce huge amounts of data in much shorter time and with less efforts and costs than the traditional Sanger technique. But the data produced comes in even shorter reads (35-100 basepairs with Illumina, 100-500 basepairs with 454-sequencing). CARMA is a new software pipeline for the characterisation of species composition and the genetic potential of microbial samples using short, unassembled reads.

Results

In this paper, we introduce WebCARMA, a refined version of CARMA available as a web application for the taxonomic and functional classification of unassembled (ultra-)short reads from metagenomic communities. In addition, we have analysed the applicability of ultra-short reads in metagenomics.

Conclusions

We show that unassembled reads as short as 35 bp can be used for the taxonomic classification of a metagenome. The web application is freely available at http://webcarma.cebitec.uni-bielefeld.de.

EMMA 2--a MAGE-compliant system for the collaborative analysis and integration of microarray data

Background

Understanding transcriptional regulation by genome-wide microarray studies can contribute to unravel complex relationships between genes. Attempts to standardize the annotation of microarray data include the Minimum Information About a Microarray Experiment (MIAME) recommendations, the MAGE-ML format for data interchange, and the use of controlled vocabularies or ontologies. The existing software systems for microarray data analysis implement the mentioned standards only partially and are often hard to use and extend. Integration of genomic annotation data and other sources of external knowledge using open standards is therefore a key requirement for future integrated analysis systems.

Results

The EMMA 2 software has been designed to resolve shortcomings with respect to full MAGE-ML and ontology support and makes use of modern data integration techniques. We present a software system that features comprehensive data analysis functions for spotted arrays, and for the most common synthesized oligo arrays such as Agilent, Affymetrix and NimbleGen. The system is based on the full MAGE object model. Analysis functionality is based on R and Bioconductor packages and can make use of a compute cluster for distributed services.

Conclusion

Our model-driven approach for automatically implementing a full MAGE object model provides high flexibility and compatibility. Data integration via SOAP-based web-services is advantageous in a distributed client-server environment as the collaborative analysis of microarray data is gaining more and more relevance in international research consortia. The adequacy of the EMMA 2 software design and implementation has been proven by its application in many distributed functional genomics projects. Its scalability makes the current architecture suited for extensions towards future transcriptomics methods based on high-throughput sequencing approaches which have much higher computational requirements than microarrays.

Protonation reactions in relation to the coupling mechanism of bovine cytochrome c oxidase

Identification of the locations of protonatable sites in cytochrome c oxidase that are influenced by reactions in the binuclear centre is critical to assessment of proposed coupling mechanisms, and to controversies on where the pumping steps occur. One such protonation site is that which governs interconversion of the isoelectronic 607 nm 'P(M)' and 580 nm 'F' forms of the two-electron-reduced oxygen intermediate. Low pH favours protonation of a site that is close to an electron paramagnetic resonance (EPR)-silent radical species in P(M), and this induces a partial electronic redistribution to form an EPR-detectable tryptophan radical in F. A further protonatable group that must be close to the binuclear centre has been detected in bacterial oxidases by Fourier transform infrared spectroscopy from pH-dependent changes in the haem-bound CO vibration frequency at low temperatures. However, in bovine cytochrome c oxidase under similar conditions of measurement, haem-bound CO remains predominantly in a single 1963 cm(-1) form between pH 6.5 and 8.5, indicating that this group is not present. Lack of pH dependence extends to the protein region of the CO photolysis spectra and suggests that both the reduced and the reduced/CO states do not have titratable groups that affect the binuclear centre strongly in the pH range 6.5-8.5. This includes the conserved glutamic acid residue E242 whose pK appears to be above 8.5 even in the fully oxidised enzyme. The results are discussed in relation to recent ideas on coupling mechanism.

Reaction of bovine cytochrome c oxidase with hydrogen peroxide produces a tryptophan cation radical and a porphyrin cation radical

Oxidized bovine cytochrome c oxidase reacts with hydrogen peroxide to generate two electron paramagnetic resonance (EPR) free radical signals (Fabian, M., and Palmer, G. (1995) Biochemistry 34, 13802-13810). These radicals are associated with the binuclear center and give rise to two overlapped EPR signals, one signal being narrower in line width (DeltaHptp = 12 G) than the other (DeltaHptp = 45 G). We have used electron nuclear double resonance (ENDOR) spectrometry to identify the two different chemical species giving rise to these two EPR signals. Comparison of the ENDOR spectrum associated with the narrow signal with that of compound I of horseradish peroxidase (formed by reaction of that enzyme with hydrogen peroxide) demonstrates that the two species are virtually identical. The chemical species giving rise to the narrow signal is therefore identified as an exchange-coupled porphyrin cation radical similar to that formed in horseradish peroxidase compound I. Comparison of the ENDOR spectrum of compound ES (formed by the reaction of hydrogen peroxide with cytochrome c peroxidase) with that of the broad signal indicates that the chemical species giving rise to the broad EPR signal in cytochrome c oxidase is probably an exchange coupled tryptophan cation radical. This is substantiated using H(2)O/D(2)O solvent exchange experiments where the ENDOR difference spectrum of the broad EPR signal of cytochrome c oxidase shows a feature consistent with hyperfine coupling to the exchangeable N(1) proton of a tryptophan cation radical.

The reactions of hydrogen peroxide with bovine cytochrome c oxidase

Oxidised cytochrome c oxidase is known to react with two molecules of hydrogen peroxide to form consecutively 607 nm 'Peroxy' and 580-nm 'Ferryl' species. These are widely used as model compounds for the equivalent P and F intermediates of the catalytic cycle. However, kinetic analysis of the reaction with H(2)O(2) in the pH range 6.0-9.0 reveals a more complex situation. In particular, as the pH is lowered, a 580-nm compound can be formed by reaction with a single H(2)O(2). This species, termed F(&z.rad;), is spectrally similar, but not identical, to F. The reactions are equivalent to those previously reported for the bo type quinol oxidase from Escherichia coli (T. Brittain, R.H. Little, C. Greenwood, N.J. Watmough, FEBS Lett. 399 (1996) 21-25) where it was proposed that F(&z.rad;) is produced directly from P. However, in the bovine oxidase F(&z.rad;) does not appear in samples of the 607-nm form, P(M), produced by CO/O(2) treatment, even at low pH, although this form is shown to be identical to the H(2)O(2)-derived P state, P(H), on the basis of spectral characteristics and kinetics of reaction with H(2)O(2). Furthermore, lowering the pH of a sample of P(M) or P(H) generated at high pH results in F(&z.rad;) formation only on a minutes time scale. It is concluded that P and F(&z.rad;) are not in a rapid, pH-dependent equilibrium, but instead are formed by distinct pathways and cannot interconvert in a simple manner, and that the crucial difference between them lies in their patterns of protonation.

Effects of mutation of the conserved glutamic acid-286 in subunit I of cytochrome c oxidase from Rhodobacter sphaeroides

We have studied the effects of mutations, E286Q and E286D, of the conserved glutamate in subunit I of cytochrome c oxidase from Rhodobacter sphaeroides with a view to evaluating the role of this residue in redox-linked proton translocation. The mutation E286D did not have any dramatic effects on enzyme properties and retained 50% of wild-type catalytic activity. For E286Q a fraction of the binuclear center was trapped in an unreactive, spectrally distinct form which is most likely due to misfolded protein, but the majority of E286Q reacted normally with formate and cyanide in the oxidized state, and with carbon monoxide and cyanide in the dithionite-reduced form. The mutation also had little effect on the pH-dependent redox properties of haem a in the reactive fraction. However, formation of the P state from oxidized enzyme with hydrogen peroxide or by aerobic incubation with carbon monoxide was inhibited. In particular, only an F-type product was obtained, at less than 25% yield, in the reaction with hydrogen peroxide. The aerobic steady state in the presence of ferrous cytochrome c was characterized by essentially fully reduced haem a and ferric haem a3, suggesting that the mutation hinders electron transfer from haem a to the binuclear center. Under these conditions or after reoxidation, on a seconds time scale, of haem a3 following anaerobiosis, there was no indication of accumulation of significant amounts of P state. We propose that the glutamate is implicated in several steps in the catalytic cycle, O --> R, P --> F, and, possibly, F --> O. The results are discussed in relation to the "glutamate trap" model for proton translocation.

On the mechanism of quinol oxidation in the bc1 complex

The question of whether significant levels of a semiquinone can be generated in the Qo site of the bc1 complex under conditions of oxidant-induced reduction is relevant to the mechanism of bifurcation of electron transfer in this site. It has already been reported that beef heart submitochondrial particles under such conditions exhibit an EPR-detectable semiquinone, which is distinct from Q-i and which was attributed to a semiquinone in the Qo site (de Vries, S., Albracht, S. P. J., Berden, J. A., and Slater, E. C. (1981) J. Biol. Chem. 256, 11996-11998). However, we show here that this signal, which can be generated to a level of around 0.1 per bc1 monomer, is insensitive to the Qo site inhibitors myxothiazol, E-beta-methoxyacrylate-stilbene, and stigmatellin, indicating that it does not arise from a Q-o species. Based on sensitivities to inhibitors of other Q sites, up to 60% of the signal may arise from semiquinones of complexes I and II. We further show that the iron-sulfur center remains EPR silent under oxidant-induced reduction conditions. Overall, the results indicate that, under conditions of oxidant-induced reduction, the Qo site is occupied primarily by quinol with the iron-sulfur center oxidized, or, possibly, by an antiferromagnetically coupled semiquinone/reduced iron-sulfur center pair, which are EPR silent. This is discussed in relation to proposed mechanisms of quinol oxidation in the Qo site, and we describe a minimal intermediate-controlled bifurcation model based on rate constants by which bifurcated electron transfer at the Qo site might occur.

Protonmotive mechanism of heme-copper oxidases

The mechanism of coupling of proton and electron transfer in oxidases is reviewed and related to the structural information that is now available. A "glutamate trap" mechanism for proton/electron coupling is described.

Effects of mutation of the conserved lysine-362 in cytochrome c oxidase from Rhodobacter sphaeroides

We describe the effects of a mutation, K362M, of the conserved lysine in cytochrome c oxidase from Rhodobacter sphaeroides, a residue located in a putative proton channel that may convey substrate protons to the binuclear center. Spectra of the "as prepared", ferricyanide-oxidized, and dithionite-reduced forms of the mutant protein confirm that the redox centers remain intact. Ligand binding kinetics of the ferricyanide-oxidized enzyme and of the dithionite-reducible fraction are similar to those of the wild type, indicating that the K channel is not the major route for CO, cyanide, formate, or peroxide entry into the structure. Protonation of the lysine residue is not redox-linked to heme a or CuB as judged from the essentially unaltered midpoint potentials of these centers in the cyanide-ligated enzyme. A difficulty in electron transfer from heme a to the binuclear center is indicated by the slow and only partial reduction of heme a3 by dithionite or ferrocytochrome c and by the presence of some reduced heme a in the as prepared mutant enzyme and under steady-state conditions. Further characterization of the K362M enzyme in the steady state shows that up to one electron, but not two, can enter the binuclear center easily. It is this inability to form the two-electron-reduced, oxygen-reactive R state that prevents activity. A model is proposed where the K channel serves as a dielectric well of high dielectric strength and low proton conductivity, rather than as a pathway for proton entry to the binuclear center. The function of this structure would be to decrease the cost of introducing a transiently uncompensated charge into the binuclear center prior to formation of a stable, charge-compensated R-state.

Effects of decyl-aurachin D and reversed electron transfer in cytochrome bd

Decyl-aurachin D is a near-stoichiometric inhibitor of cytochrome bd from Azotobacter vinelandii. Interaction of decyl-aurachin D with the oxidase induces a redshift of the alpha-band and Soret band of a b-type cytochrome, probably b-558, suggesting close proximity of the inhibitor binding site to this haem and hence to the proposed quinol binding domain. The compound does not affect the oxygen binding site directly as judged from unchanged CO recombination kinetics to haem d in dithionite-reduced enzyme. Although in the presence of ubiquinol-1 a decyl-aurachin D containing sample generates levels of haem reduction and catalytic intermediates similar to the control, the approach to this steady state is severely inhibited. In addition to the spectral effect on b-558, decyl-aurachin D raises the midpoint potential of haem b-558, but also lowers that of haem b-595. Consistent with the shift in midpoint potentials, electron backflow from haem d to the b-type haems can be observed in decyl-aurachin D inhibited samples following photolysis of the mixed-valence CO-ligated form of the enzyme. The data show that decyl-aurachin D acts on the donor side of haem b-558 without substantially affecting internal electron transfer rates or the oxygen reduction site.

A suggested mechanism for the catalytic cycle of cytochrome bd terminal oxidase based on kinetic analysis

The apparent oxygen affinity of cytochrome bd from Escherichia coli and Azotobacter vinelandii has been measured using oxymyoglobin as a sensitive monitor of oxygen concentration. In membrane preparations, the Km(O2) and respiratory rate varied with the nature of the primary substrate used (malate, lactate, reduced nicotinamide adenine dinucleotide (NADH), or ubiquinol-1). At maximum respiratory rates, the Km(O2) for cytochrome bd from A. vinelandii was 4.1 microM, approximately 2 times higher than the corresponding value for the E. coli enzyme. There were no significant differences between the Km(O2) values for membrane-bound and purified cytochrome bd from A. vinelandii when ubiquinol-1 was used as primary substrate. The kinetic parameters Km(O2) and Vmax provide a value of 2.8 x 10(8) M-1 s-1 for the bimolecular rate constant for oxygen reaction with the enzyme, suggesting that this reaction is diffusion-controlled. Kinetic analysis indicates a mechanism involving a ternary complex. A scheme for the reaction mechanism of cytochrome bd is proposed.

Cytochrome bd oxidase from Azotobacter vinelandii. Purification and quantitation of ligand binding to the oxygen reduction site

Cytochrome bd has been purified from Azotobacter vinelandii by a new simplified procedure. The heme and total iron content has been measured, as has the number of high affinity CO and NO binding sites. Spectral changes indicate high affinity binding of CO and NO to heme d only, with a stoichiometry of 1 molecule of gas per 2 molecules of heme b or per 3 atoms of iron. The results clearly define a stoichiometry of one heme d per complex. Low affinity binding of CO and NO to heme b595 also occurs at higher ligand concentrations. EPR heme-nitrosyl signals are seen with NO bound to both hemes b595 and d but with no indication of spin exchange coupling. Exposure of the air-oxidized complex to alkaline pH results in removal of molecular oxygen from heme d and a change in line shape of the high spin region of the EPR spectrum. Cyanide binds to both heme d and heme b595 in the air-oxidized complex, displacing molecular oxygen from heme d. The rate of cyanide binding to heme d as assessed by spectral changes at 650 nm does not correlate with the rate of binding to heme b595 as assessed by the loss of the high spin EPR signal. In addition, the cyanide binding rate in the presence of reductant is only 3 times that of the rate of binding to the air-oxidized enzyme, in contrast to the copper-containing oxidases where strong redox cooperativity makes these two rates differ by a factor of at least 10(6). The results do not support the idea of the presence of two strongly interacting hemes in a binuclear center.

Antimycin inhibition of the cytochrome bd complex from Azotobacter vinelandii indicates the presence of a branched electron transfer pathway for the oxidation of ubiquinol

Antimycin A and UHBDT inhibit the activity of the purified cytochrome bd complex from Azotobacter vinelandii. Inhibition of activity is non-competitive and antimycin A binding induces a shift to the red in the spectrum of a b-type haem. No inhibitory effects were seen with myxothiazol. Steady-state experiments indicate that the site of inhibition for antimycin A lies on the low-potential side of haem b558. In the presence of antimycin A at concentrations sufficient to inhibit respiration, some direct electron transfer from ubiquinol-1 to haem b595 and haem d still occurs. The results are consistent with a branched electron transfer pathway from ubiquinol to the oxygen reduction site.

Slow ('resting') forms of mitochondrial cytochrome c oxidase consist of two kinetically distinct conformations of the binuclear CuB/a3 centre--relevance to the mechanism of proton translocation

We have purified slow ('resting') cytochrome oxidase from bovine heart, free of contamination with fast ('pulsed') enzyme. This form of the enzyme shows two kinetic phases of reduction of haem a3 by dithionite (k = 0.020 +/- 0.005 s-1 and k = 0.005 +/- 0.002 s-1). The presence of ligands that bind to the oxidized or reduced binuclear centre (formate or carbon monoxide respectively) has no effect on these rates. Varying the dithionite concentration also has no effect on either phase, although at low dithionite concentrations a lag phase is observed as the rate of haem a reduction is slower. The results are consistent with a model for reduction of the slow enzyme where the rate of electron transfer to the binuclear centre is the limiting step, rather than an equilibrium model where the haem a3 redox potential is low. Increasing the pH decreases the rate of the slower phase of dithionite reduction, but has no effect on the faster phase. EPR studies show that the slow phase (only) correlates with the disappearance of the g' = 12/g' = 2.95 signals, with the same pH dependence; again the presence of formate has no effect on these results. Deconvolution of the oxidized optical spectra shows that the enzyme reduced in the slow phase has a blue-shifted Soret band, relative to that reduced in the faster phase. Incubation of the oxidized enzyme at high pH causes a line-broadening of both the g' = 12 and g' = 2.95 EPR signals with no obvious effect on the amount of signal. The results are interpreted in a model where the presence of a carboxylate bridge between haem a3 and CuB defines the slow enzyme. It is suggested that the two rates of dithionite reduction are the result of different ligation to CuB--where water is the ligand the binuclear centre is FeIV/CuI (EPR-silent) and where hydroxide is the ligand the binuclear centre is FeIII/CuII (g' = 12/g' = 2.95 EPR signals).