Towards Exploring Toxin-Antitoxin Systems in Geobacillus: A Screen for Type II Toxin-Antitoxin System Families in a Thermophilic Genus.

The toxin-antitoxin (TA) systems have been attracting attention due to their role in regulating stress responses in prokaryotes and their biotechnological potential. Much recognition has been given to type II TA system of mesophiles, while thermophiles have received merely limited attention. Here, we are presenting the putative type II TA families encoded on the genomes of four Geobacillus strains. We employed the TA finder tool to mine for TA-coding genes and manually curated the results using protein domain analysis tools. We also used the NCBI BLAST, Operon Mapper, ProOpDB, and sequence alignment tools to reveal the geobacilli TA features. We identified 28 putative TA pairs, distributed over eight TA families. Among the identified TAs, 15 represent putative novel toxins and antitoxins, belonging to the MazEF, MNT-HEPN, ParDE, RelBE, and XRE-COG2856 TA families. We also identified a potentially new TA composite, AbrB-ParE. Furthermore, we are suggesting the Geobacillus acetyltransferase TA (GacTA) family, which potentially represents one of the unique TA families with a reverse gene order. Moreover, we are proposing a hypothesis on the xre-cog2856 gene expression regulation, which seems to involve the c-di-AMP. This study aims for highlighting the significance of studying TAs in Geobacillus and facilitating future experimental research.

Contrasting prevalence of selection and drift in the community structuring of bacteria and microbial eukaryotes.

Whether or not communities of microbial eukaryotes are structured in the same way as bacteria is a general and poorly explored question in ecology. Here, we investigated this question in a set of planktonic lake microbiotas in Eastern Antarctica that represent a natural community ecology experiment. Most of the analysed lakes emerged from the sea during the last 6000 years, giving rise to waterbodies that originally contained marine microbiotas and that subsequently evolved into habitats ranging from freshwater to hypersaline. We show that habitat diversification has promoted selection driven by the salinity gradient in bacterial communities (explaining ∼ 72% of taxa turnover), while microeukaryotic counterparts were predominantly structured by ecological drift (∼72% of the turnover). Nevertheless, we also detected a number of microeukaryotes with specific responses to salinity, indicating that albeit minor, selection has had a role in the structuring of specific members of their communities. In sum, we conclude that microeukaryotes and bacteria inhabiting the same communities can be structured predominantly by different processes. This should be considered in future studies aiming to understand the mechanisms that shape microbial assemblages.

Identification of Putative Novel Class-I Lanthipeptides in Firmicutes: A Combinatorial In Silico Analysis Approach Performed on Genome Sequenced Bacteria and a Close Inspection of Z-Geobacillin Lanthipeptide Biosynthesis Gene Cluster of the Thermophilic Geobacillus sp. Strain ZGt-1.

Lanthipeptides are ribosomally synthesized and post-translationally modified polycyclic peptides. Lanthipeptides that have antimicrobial activity are known as lantibiotics. Accordingly, the discovery of novel lantibiotics constitutes a possible solution for the problem of antibiotic resistance. We utilized the publicly available genome sequences and the bioinformatic tools tailored for the detection of lanthipeptides. We designed our strategy for screening of 252 firmicute genomes and detecting class-I lanthipeptide-coding gene clusters. The designed strategy resulted in identifying 69 class-I lanthipeptide sequences, of which more than 10% were putative novel. The identified putative novel lanthipeptides have not been annotated on the original or the RefSeq genomes, or have been annotated merely as coding for hypothetical proteins. Additionally, we identified bacterial strains that have not been previously recognized as lanthipeptide-producers. Moreover, we suggest corrections for certain firmicute genome annotations, and recommend lanthipeptide records for enriching the bacteriocin genome mining tool (BAGEL) databases. Furthermore, we propose Z-geobacillin, a putative class-I lanthipeptide coded on the genome of the thermophilic strain Geobacillus sp. ZGt-1. We provide lists of putative novel lanthipeptide sequences and of the previously unrecognized lanthipeptide-producing bacterial strains, so they can be prioritized for experimental investigation. Our results are expected to benefit researchers interested in the in vitro production of lanthipeptides.

Ectomycorrhizal fungi decompose soil organic matter using oxidative mechanisms adapted from saprotrophic ancestors.

Ectomycorrhizal fungi are thought to have a key role in mobilizing organic nitrogen that is trapped in soil organic matter (SOM). However, the extent to which ectomycorrhizal fungi decompose SOM and the mechanism by which they do so remain unclear, considering that they have lost many genes encoding lignocellulose-degrading enzymes that are present in their saprotrophic ancestors. Spectroscopic analyses and transcriptome profiling were used to examine the mechanisms by which five species of ectomycorrhizal fungi, representing at least four origins of symbiosis, decompose SOM extracted from forest soils. In the presence of glucose and when acquiring nitrogen, all species converted the organic matter in the SOM extract using oxidative mechanisms. The transcriptome expressed during oxidative decomposition has diverged over evolutionary time. Each species expressed a different set of transcripts encoding proteins associated with oxidation of lignocellulose by saprotrophic fungi. The decomposition 'toolbox' has diverged through differences in the regulation of orthologous genes, the formation of new genes by gene duplications, and the recruitment of genes from diverse but functionally similar enzyme families. The capacity to oxidize SOM appears to be common among ectomycorrhizal fungi. We propose that the ancestral decay mechanisms used primarily to obtain carbon have been adapted in symbiosis to scavenge nutrients instead.

Antimicrobial Protein Candidates from the Thermophilic Geobacillus sp. Strain ZGt-1: Production, Proteomics, and Bioinformatics Analysis.

A thermophilic bacterial strain, Geobacillus sp. ZGt-1, isolated from Zara hot spring in Jordan, was capable of inhibiting the growth of the thermophilic G. stearothermophilus and the mesophilic Bacillus subtilis and Salmonella typhimurium on a solid cultivation medium. Antibacterial activity was not observed when ZGt-1 was cultivated in a liquid medium; however, immobilization of the cells in agar beads that were subjected to sequential batch cultivation in the liquid medium at 60 °C showed increasing antibacterial activity up to 14 cycles. The antibacterial activity was lost on protease treatment of the culture supernatant. Concentration of the protein fraction by ammonium sulphate precipitation followed by denaturing polyacrylamide gel electrophoresis separation and analysis of the gel for antibacterial activity against G. stearothermophilus showed a distinct inhibition zone in 15-20 kDa range, suggesting that the active molecule(s) are resistant to denaturation by SDS. Mass spectrometric analysis of the protein bands around the active region resulted in identification of 22 proteins with molecular weight in the range of interest, three of which were new and are here proposed as potential antimicrobial protein candidates by in silico analysis of their amino acid sequences. Mass spectrometric analysis also indicated the presence of partial sequences of antimicrobial enzymes, amidase, and dd-carboxypeptidase.

Species and gene divergence in Littorina snails detected by array comparative genomic hybridization.

BACKGROUND: Array comparative genomic hybridization (aCGH) is commonly used to screen different types of genetic variation in humans and model species. Here, we performed aCGH using an oligonucleotide gene-expression array for a non-model species, the intertidal snail Littorina saxatilis. First, we tested what types of genetic variation can be detected by this method using direct re-sequencing and comparison to the Littorina genome draft. Secondly, we performed a genome-wide comparison of four closely related Littorina species: L. fabalis, L. compressa, L. arcana and L. saxatilis and of populations of L. saxatilis found in Spain, Britain and Sweden. Finally, we tested whether we could identify genetic variation underlying "Crab" and "Wave" ecotypes of L. saxatilis. RESULTS: We could reliably detect copy number variations, deletions and high sequence divergence (i.e. above 3%), but not single nucleotide polymorphisms. The overall hybridization pattern and number of significantly diverged genes were in close agreement with earlier phylogenetic reconstructions based on single genes. The trichotomy of L. arcana, L. compressa and L. saxatilis could not be resolved and we argue that these divergence events have occurred recently and very close in time. We found evidence for high levels of segmental duplication in the Littorina genome (10% of the transcripts represented on the array and up to 23% of the analyzed genomic fragments); duplicated genes and regions were mostly the same in all analyzed species. Finally, this method discriminated geographically distant populations of L. saxatilis, but we did not detect any significant genome divergence associated with ecotypes of L. saxatilis. CONCLUSIONS: The present study provides new information on the sensitivity and the potential use of oligonucleotide arrays for genotyping of non-model organisms. Applying this method to Littorina species yields insights into genome evolution following the recent species radiation and supports earlier single-gene based phylogenies. Genetic differentiation of L. saxatilis ecotypes was not detected in this study, despite pronounced innate phenotypic differences. The reason may be that these differences are due to single-nucleotide polymorphisms.

Characterisation of a transcriptome to find sequence differences between two differentially migrating subspecies of the willow warbler Phylloscopus trochilus.

BACKGROUND: Animal migration requires adaptations in morphological, physiological and behavioural traits. Several of these traits have been shown to possess a strong heritable component in birds, but little is known about their genetic architecture. Here we used 454 sequencing of brain-derived transcriptomes from two differentially migrating subspecies of the willow warbler Phylloscopus trochilus to detect genes potentially underlying traits associated with migration. RESULTS: The transcriptome sequencing resulted in 1.8 million reads following filtering steps. Most of the reads (84%) were successfully mapped to the genome of the zebra finch Taeniopygia gutatta. The mapped reads were situated within at least 12,101 predicted zebra finch genes, with the greatest sequencing depth in exons. Reads that were mapped to intergenic regions were generally located close to predicted genes and possibly located in uncharacterized untranslated regions (UTRs). Out of 85,000 single nucleotide polymorphisms (SNPs) with a minimum sequencing depth of eight reads from each of two subspecies-specific pools, only 55 showed high differentiation, confirming previous studies showing that most of the genetic variation is shared between the subspecies. Validation of a subset of the most highly differentiated SNPs using Sanger sequencing demonstrated that several of them also were differentiated between an independent set of individuals of each subspecies. These SNPs were clustered in two chromosome regions that are likely to be influenced by divergent selection between the subspecies and that could potentially be associated with adaptations to their different migratory strategies. CONCLUSIONS: Our study represents the first large-scale sequencing analysis aiming at detecting genes underlying migratory phenotypes in birds and provides new candidates for genes potentially involved in migration.

The molecular components of the extracellular protein-degradation pathways of the ectomycorrhizal fungus Paxillus involutus.

Proteins contribute to a major part of the organic nitrogen (N) in forest soils. This N is mobilized and becomes available to trees as a result of the depolymerizing activities of symbiotic ectomycorrhizal fungi. The mechanisms by which these fungi depolymerize proteins and assimilate the released N are poorly characterized. Biochemical analysis and transcriptome profiling were performed to examine the proteolytic machinery and the uptake system of the ectomycorrhizal basidiomycete Paxillus involutus during the assimilation of organic N from various protein sources and extracts of organic matter. All substrates induced secretion of peptidase activity with an acidic pH optimum, mostly contributed by aspartic peptidases. The peptidase activity was transiently repressed by ammonium. Transcriptional analysis revealed a large number of extracellular endo- and exopeptidases. The expression levels of these peptidases were regulated in parallel with transporters and enzymes involved in the assimilation and metabolism of the released peptides and amino acids. For the first time the molecular components of the protein degradation pathways of an ectomycorrhizal fungus are described. The data suggest that the transcripts encoding these components are regulated in response to the chemical properties and the availability of the protein substrates.

Insights into evolution of multicellular fungi from the assembled chromosomes of the mushroom Coprinopsis cinerea (Coprinus cinereus).

The mushroom Coprinopsis cinerea is a classic experimental model for multicellular development in fungi because it grows on defined media, completes its life cycle in 2 weeks, produces some 10(8) synchronized meiocytes, and can be manipulated at all stages in development by mutation and transformation. The 37-megabase genome of C. cinerea was sequenced and assembled into 13 chromosomes. Meiotic recombination rates vary greatly along the chromosomes, and retrotransposons are absent in large regions of the genome with low levels of meiotic recombination. Single-copy genes with identifiable orthologs in other basidiomycetes are predominant in low-recombination regions of the chromosome. In contrast, paralogous multicopy genes are found in the highly recombining regions, including a large family of protein kinases (FunK1) unique to multicellular fungi. Analyses of P450 and hydrophobin gene families confirmed that local gene duplications drive the expansions of paralogous copies and the expansions occur in independent lineages of Agaricomycotina fungi. Gene-expression patterns from microarrays were used to dissect the transcriptional program of dikaryon formation (mating). Several members of the FunK1 kinase family are differentially regulated during sexual morphogenesis, and coordinate regulation of adjacent duplications is rare. The genomes of C. cinerea and Laccaria bicolor, a symbiotic basidiomycete, share extensive regions of synteny. The largest syntenic blocks occur in regions with low meiotic recombination rates, no transposable elements, and tight gene spacing, where orthologous single-copy genes are overrepresented. The chromosome assembly of C. cinerea is an essential resource in understanding the evolution of multicellularity in the fungi.

The ectomycorrhizal fungus Paxillus involutus converts organic matter in plant litter using a trimmed brown-rot mechanism involving Fenton chemistry.

Soils in boreal forests contain large stocks of carbon. Plants are the main source of this carbon through tissue residues and root exudates. A major part of the exudates are allocated to symbiotic ectomycorrhizal fungi. In return, the plant receives nutrients, in particular nitrogen from the mycorrhizal fungi. To capture the nitrogen, the fungi must at least partly disrupt the recalcitrant organic matter-protein complexes within which the nitrogen is embedded. This disruption process is poorly characterized. We used spectroscopic analyses and transcriptome profiling to examine the mechanism by which the ectomycorrhizal fungus Paxillus involutus degrades organic matter when acquiring nitrogen from plant litter. The fungus partially degraded polysaccharides and modified the structure of polyphenols. The observed chemical changes were consistent with a hydroxyl radical attack, involving Fenton chemistry similar to that of brown-rot fungi. The set of enzymes expressed by Pa. involutus during the degradation of the organic matter was similar to the set of enzymes involved in the oxidative degradation of wood by brown-rot fungi. However, Pa. involutus lacked transcripts encoding extracellular enzymes needed for metabolizing the released carbon. The saprotrophic activity has been reduced to a radical-based biodegradation system that can efficiently disrupt the organic matter-protein complexes and thereby mobilize the entrapped nutrients. We suggest that the released carbon then becomes available for further degradation and assimilation by commensal microbes, and that these activities have been lost in ectomycorrhizal fungi as an adaptation to symbiotic growth on host photosynthate. The interdependence of ectomycorrhizal symbionts and saprophytic microbes would provide a key link in the turnover of nutrients and carbon in forest ecosystems.

Major histocompatibility complex class II compatibility, but not class I, predicts mate choice in a bird with highly developed olfaction.

Mate choice for major histocompatibility complex (MHC) compatibility has been found in several taxa, although rarely in birds. MHC is a crucial component in adaptive immunity and by choosing an MHC-dissimilar partner, heterozygosity and potentially broad pathogen resistance is maximized in the offspring. The MHC genotype influences odour cues and preferences in mammals and fish and hence olfactory-based mate choice can occur. We tested whether blue petrels, Halobaena caerulea, choose partners based on MHC compatibility. This bird is long-lived, monogamous and can discriminate between individual odours using olfaction, which makes it exceptionally well suited for this analysis. We screened MHC class I and II B alleles in blue petrels using 454-pyrosequencing and quantified the phylogenetic, functional and allele-sharing similarity between individuals. Partners were functionally more dissimilar at the MHC class II B loci than expected from random mating (p = 0.033), whereas there was no such difference at the MHC class I loci. Phylogenetic and non-sequence-based MHC allele-sharing measures detected no MHC dissimilarity between partners for either MHC class I or II B. Our study provides evidence of mate choice for MHC compatibility in a bird with a high dependency on odour cues, suggesting that MHC odour-mediated mate choice occurs in birds.

Recruitment of members from the rare biosphere of marine bacterioplankton communities after an environmental disturbance.

A bacterial community may be resistant to environmental disturbances if some of its species show metabolic flexibility and physiological tolerance to the changing conditions. Alternatively, disturbances can change the composition of the community and thereby potentially affect ecosystem processes. The impact of disturbance on the composition of bacterioplankton communities was examined in continuous seawater cultures. Bacterial assemblages from geographically closely connected areas, the Baltic Sea (salinity 7 and high dissolved organic carbon [DOC]) and Skagerrak (salinity 28 and low DOC), were exposed to gradual opposing changes in salinity and DOC over a 3-week period such that the Baltic community was exposed to Skagerrak salinity and DOC and vice versa. Denaturing gradient gel electrophoresis and clone libraries of PCR-amplified 16S rRNA genes showed that the composition of the transplanted communities differed significantly from those held at constant salinity. Despite this, the growth yields (number of cells ml(-1)) were similar, which suggests similar levels of substrate utilization. Deep 454 pyrosequencing of 16S rRNA genes showed that the composition of the disturbed communities had changed due to the recruitment of phylotypes present in the rare biosphere of the original community. The study shows that members of the rare biosphere can become abundant in a bacterioplankton community after disturbance and that those bacteria can have important roles in maintaining ecosystem processes.

Insight into trade-off between wood decay and parasitism from the genome of a fungal forest pathogen.

Parasitism and saprotrophic wood decay are two fungal strategies fundamental for succession and nutrient cycling in forest ecosystems. An opportunity to assess the trade-off between these strategies is provided by the forest pathogen and wood decayer Heterobasidion annosum sensu lato. We report the annotated genome sequence and transcript profiling, as well as the quantitative trait loci mapping, of one member of the species complex: H. irregulare. Quantitative trait loci critical for pathogenicity, and rich in transposable elements, orphan and secreted genes, were identified. A wide range of cellulose-degrading enzymes are expressed during wood decay. By contrast, pathogenic interaction between H. irregulare and pine engages fewer carbohydrate-active enzymes, but involves an increase in pectinolytic enzymes, transcription modules for oxidative stress and secondary metabolite production. Our results show a trade-off in terms of constrained carbohydrate decomposition and membrane transport capacity during interaction with living hosts. Our findings establish that saprotrophic wood decay and necrotrophic parasitism involve two distinct, yet overlapping, processes.

Phylogenetic analysis suggests that habitat filtering is structuring marine bacterial communities across the globe.

The phylogenetic structure and community composition were analysed in an existing data set of marine bacterioplankton communities to elucidate the evolutionary and ecological processes dictating the assembly. The communities were sampled from coastal waters at nine locations distributed worldwide and were examined through the use of comprehensive clone libraries of 16S ribosomal RNA genes. The analyses show that the local communities are phylogenetically different from each other and that a majority of them are phylogenetically clustered, i.e. the species (operational taxonomic units) were more related to each other than expected by chance. Accordingly, the local communities were assembled non-randomly from the global pool of available bacterioplankton. Further, the phylogenetic structures of the communities were related to the water temperature at the locations. In agreement with similar studies, including both macroorganisms and bacteria, these results suggest that marine bacterial communities are structured by "habitat filtering", i.e. through non-random colonization and invasion determined by environmental characteristics. Different bacterial types seem to have different ecological niches that dictate their survival in different habitats. Other eco-evolutionary processes that may contribute to the observed phylogenetic patterns are discussed. The results also imply a mapping between phenotype and phylogenetic relatedness which facilitates the use of community phylogenetic structure analysis to infer ecological and evolutionary assembly processes.

RAMI: a tool for identification and characterization of phylogenetic clusters in microbial communities.

MOTIVATION: The most common approach to estimate microbial diversity is based on the analysis of DNA sequences of specific target genes including ribosomal genes. Commonly, the sequences are grouped into operational taxonomic units based on genetic distance (sequence similarity) instead of genetic change (patristic distance). This method may fail to adequately identify clusters of evolutionary related sequences and it provides no information on the phylogenetic structure of the community. An ease-of-use web application for this purpose has been missing. RESULTS: We have developed RAMI, which clusters related nodes in a phylogenetic tree based on the patristic distance. RAMI also produces indices of cluster properties and other indices used in population and community studies on-the-fly. AVAILABILITY: RAMI is licensed under GNU GPL and can be run or downloaded from http://www.acgt.se/online.html. SUPPLEMENTARY INFORMATION: http://www.acgt.se/RAMI/SuppInfo.

ARWEN: a program to detect tRNA genes in metazoan mitochondrial nucleotide sequences.

MOTIVATION: Mitochondrial genomes encode their own transfer RNAs (tRNAs). These are often degenerate in sequence and structure compared to tRNAs in their bacterial ancestors. This is one of the reasons why current tRNA gene predictor programs perform poorly identifying mitochondrial tRNA genes. As a consequence there is a need for a new program with the specific aim of predicting these tRNAs. RESULTS: In this study, we present the software ARWEN that identifies tRNA genes in metazoan mitochondrial nucleotide sequences. ARWEN detects close to 100% of previously annotated genes. AVAILABILITY: An online version, software for download and test results are available at www.acgt.se/online.html

ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences.

A computer program, ARAGORN, identifies tRNA and tmRNA genes. The program employs heuristic algorithms to predict tRNA secondary structure, based on homology with recognized tRNA consensus sequences and ability to form a base-paired cloverleaf. tmRNA genes are identified using a modified version of the BRUCE program. ARAGORN achieves a detection sensitivity of 99% from a set of 1290 eubacterial, eukaryotic and archaeal tRNA genes and detects all complete tmRNA sequences in the tmRNA database, improving on the performance of the BRUCE program. Recently discovered tmRNA genes in the chloroplasts of two species from the 'green' algae lineage are detected. The output of the program reports the proposed tRNA secondary structure and, for tmRNA genes, the secondary structure of the tRNA domain, the tmRNA gene sequence, the tag peptide and a list of organisms with matching tmRNA peptide tags.

BRUCE: a program for the detection of transfer-messenger RNA genes in nucleotide sequences.

A computer program, BRUCE, was developed for the identification of transfer-messenger RNA (tmRNA) genes. The program employs heuristic algorithms to search for a tRNA(Ala)-like secondary structure surrounding a short sequence encoding the tag peptide. In the 57 completely sequenced bacterial genomes where tmRNA genes have been reported previously, BRUCE identified all with no false positives. In addition, BRUCE found 99 of the 100 tmRNAs identified previously in other bacteria, red chloroplasts and cyanelles. The output of the program reports the proposed tRNA secondary structure, the tmRNA gene sequence and the tag peptide.

The Genome of Haemoproteus tartakovskyi and Its Relationship to Human Malaria Parasites.

The phylogenetic relationships among hemosporidian parasites, including the origin of Plasmodium falciparum, the most virulent malaria parasite of humans, have been heavily debated for decades. Studies based on multiple-gene sequences have helped settle many of these controversial phylogenetic issues. However, denser taxon sampling and genome-wide analyses are needed to confidently resolve the evolutionay relationships among hemosporidian parasites. Genome sequences of several Plasmodium parasites are available but only for species infecting primates and rodents. To root the phylogenetic tree of Plasmodium, genomic data from related parasites of birds or reptiles are required. Here, we use a novel approach to isolate parasite DNA from microgametes and describe the first genome of a bird parasite in the sister genus to Plasmodium, Haemoproteus tartakovskyi Similar to Plasmodium parasites, H. tartakovskyi has a small genome (23.2 Mb, 5,990 genes) and a GC content (25.4%) closer to P. falciparum (19.3%) than to Plasmodium vivax (42.3%). Combined with novel transcriptome sequences of the bird parasite Plasmodium ashfordi, our phylogenomic analyses of 1,302 orthologous genes demonstrate that mammalian-infecting malaria parasites are monophyletic, thus rejecting the repeatedly proposed hypothesis that the ancestor of Laverania parasites originated from a secondary host shift from birds to humans. Genes and genomic features previously found to be shared between P. falciparum and bird malaria parasites, but absent in other mammal malaria parasites, are therefore signatures of maintained ancestral states. We foresee that the genome of H. tartakovskyi will open new directions for comparative evolutionary analyses of malarial adaptive traits.

Genome Sequence of Geobacillus sp. Strain ZGt-1, an Antibacterial Peptide-Producing Bacterium from Hot Springs in Jordan.

This paper reports the draft genome sequence of the firmicute Geobacillus sp. strain ZGt-1, an antibacterial peptide producer isolated from the Zara hot spring in Jordan. This study is the first report on genomic data from a thermophilic bacterial strain isolated in Jordan.