Search

Following 30 years of sequencing, we assessed the phylogenetic diversity (PD) of >1.5 million microbial genomes in public databases, including metagenome-assembled genomes (MAGs) of uncultivated microbes. As compared to the vast diversity uncovered by metagenomic sequences, cultivated taxa account for a modest portion of the overall diversity, 9.73% in bacteria and 6.55% in archaea, while MAGs contribute 48.54% and 57.05%, respectively. Therefore, a substantial fraction of bacterial (41.73%) and archaeal PD (36.39%) still lacks any genomic representation. This unrepresented diversity manifests primarily at lower taxonomic ranks, exemplified by 134,966 species identified in 18,087 metagenomic samples. Our study exposes diversity hotspots in freshwater, marine subsurface, sediment, soil, and other environments, whereas human samples yielded minimal novelty within the context of existing datasets. These results offer a roadmap for future genome recovery efforts, delineating uncaptured taxa in underexplored environments and underscoring the necessity for renewed isolation and sequencing. [ABSTRACT FROM AUTHOR]

We present two standards developed by the Genomic Standards Consortium (GSC) for reporting bacterial and archaeal genome sequences. Both are extensions of the Minimum Information about Any (x) Sequence (MIxS). The standards are the Minimum Information about a Single Amplified Genome (MISAG) and the Minimum Information about a Metagenome-Assembled Genome (MIMAG), including, but not limited to, assembly quality, and estimates of genome completeness and contamination. These standards can be used in combination with other GSC checklists, including the Minimum Information about a Genome Sequence (MIGS), Minimum Information about a Metagenomic Sequence (MIMS), and Minimum Information about a Marker Gene Sequence (MIMARKS). Community-wide adoption of MISAG and MIMAG will facilitate more robust comparative genomic analyses of bacterial and archaeal diversity. [ABSTRACT FROM AUTHOR]

The Genomic Encyclopedia of Bacteria and Archaea (GEBA) project was launched by the JGI in 2007 as a pilot project with the objective of sequencing 250 bacterial and archaeal genomes. The two major goals of that project were (a) to test the hypothesis that there are many benefits to the use the phylogenetic diversity of organisms in the tree of life as a primary criterion for generating their genome sequence and (b) to develop the necessary framework, technology and organization for large-scale sequencing of microbial isolate genomes. While the GEBA pilot project has not yet been entirely completed, both of the original goals have already been successfully accomplished, leading the way for the next phase of the project. Here we propose taking the GEBA project to the next level, by generating high quality draft genomes for 1,000 bacterial and archaeal strains. This represents a combined 16-fold increase in both scale and speed as compared to the GEBA pilot project (250 isolate genomes in 4+ years). We will follow a similar approach for organism selection and sequencing prioritization as was done for the GEBA pilot project (i.e. phylogenetic novelty, availability and growth of cultures of type strains and DNA extraction capability), focusing on type strains as this ensures reproducibility of our results and provides the strongest linkage between genome sequences and other knowledge about each strain. In turn, this project will constitute a pilot phase of a larger effort that will target the genome sequences of all available type strains of the Bacteria and Archaea. [ABSTRACT FROM AUTHOR]

This manuscript calls for an international effort to generate a comprehensive catalog from genome sequences of all the archaeal and bacterial type strains. [ABSTRACT FROM AUTHOR]

Microbiology is experiencing a revolution brought on by recent developments in sequencing technology. The unprecedented volume of microbiome data being generated poses significant challenges that are currently hindering progress in the field. Here, we outline the major bottlenecks and propose a vision to advance microbiome research as a data-driven science. [ABSTRACT FROM AUTHOR]

The article discusses various issues concerning the field of microbial genomics. The author notes that the potential offered by microbial genomics is associated with the need to meet the challenges brought by the rapid development of underlying technology and the growth of data. With the advent of new technologies and the experiences gained in the field of genomics, microbiologists are said to be poised to explore and understand the microbial aspects of the Earth. The article explains the issues concerning the genome project selection bias and the maintenance of biological databases. It also explores the lack of standards in the field of genomics.

The exponential increase of genomic and metagenomic data, fueled in part by recent advancements in sequencing technology, are greatly expanding our understanding of the phylogenetic diversity and metabolic capacity present in the environment. Two of the central challenges that bioinformaticians and ecologists alike must face are the design of bioinformatic resources that facilitate the analysis of genomic and metagenomic data in a comparative context and the efficient capture and organization of the plethora of descriptive information required to usefully describe these data sets. In this commentary, we review three initiatives presented in the “new frontiers” session of the second SCOPE meeting on Microbial Environmental Genomics (MicroEnGen-II, Shanghai, June 12–15, 2006). These are (1) the Integrated Microbial Genomes Resources (IMG), (2) the Genomic Standards Consortium (GSC), and (3) the Natural Environment Research Council (NERC) Environmental Bioinformatics Centre (NEBC). These integrative bioinformatics and data management initiatives underscore the increasingly important role ecologists have to play in the genomic (metagenomic) revolution. [ABSTRACT FROM AUTHOR]

Presents information on a study which quantified the phylogenetic extent of the predicted transcription-associated proteins in Archaea. Materials and methods of the study; Results and discussion.

Focuses on a complex case of the translation initiation factor which leads the initiator transfer ribonucleic acid (tRNA) into the ribosome, while exploring the evolutionary relationships among the components of bacterial initiation factor 2 (IF-2) and eukaryotic IF-2). Materials and methods used; Definition of translation initiation; Results and discussion.

Looks at three processes by which translation is initiated in cells. Background information on these processes; Description of a protein sequencing database in relation to this topic; Discussion on the results from these processes.

The article discusses various reports published within the issue including one on individual bacterial or archaeal phyla in sequenced genomes, one on commentaries about research on microbial genomics or metagenomics and one on cultured microorganisms.

Provides information on the physiological role of tetratrico-peptide-repeat (TPR) proteins in Methanococcus jannaschii organism. Phylogenetic distribution of TPR amino-acid motif; Functional importance and variability; Structural components of TPR proteins.

Identifying and characterizing mobile genetic elements in sequencing data is essential for understanding their diversity, ecology, biotechnological applications and impact on public health. Here we introduce geNomad, a classification and annotation framework that combines information from gene content and a deep neural network to identify sequences of plasmids and viruses. geNomad uses a dataset of more than 200,000 marker protein profiles to provide functional gene annotation and taxonomic assignment of viral genomes. Using a conditional random field model, geNomad also detects proviruses integrated into host genomes with high precision. In benchmarks, geNomad achieved high classification performance for diverse plasmids and viruses (Matthews correlation coefficient of 77.8% and 95.3%, respectively), substantially outperforming other tools. Leveraging geNomad's speed and scalability, we processed over 2.7 trillion base pairs of sequencing data, leading to the discovery of millions of viruses and plasmids that are available through the IMG/VR and IMG/PR databases. geNomad is available at https://portal.nersc.gov/genomad. geNomad identifies mobile genetic elements in sequencing data. [ABSTRACT FROM AUTHOR]

Root nodule bacteria (RNB) or 'rhizobia' are a type of plant growth promoting bacteria, typified by their ability to fix nitrogen for their plant host, fixing nearly 65% of the nitrogen currently utilized in sustainable agricultural production of legume crops and pastures. In this study, we sequenced the genomes of 110 RNB from diverse hosts and biogeographical regions, and undertook a global exploration of all available RNB genera with the aim of identifying novel genetic determinants of symbiotic association and plant growth promotion. Specifically, we performed a subtractive comparative analysis with non-RNB genomes, employed relevant transcriptomic data, and leveraged phylogenetic distribution patterns and sequence signatures based on known precepts of symbiotic- and host-microbe interactions. A total of 184 protein families were delineated, including known factors for nodulation and nitrogen fixation, and candidates with previously unexplored functions, for which a role in host-interaction, -regulation, biocontrol, and more, could be posited. These analyses expand our knowledge of the RNB purview and provide novel targets for strain improvement in the ultimate quest to enhance plant productivity and agricultural sustainability. [ABSTRACT FROM AUTHOR]

Presents a comments on the annotation of genome sequences. Comparison of various selected cases of potential genome sequences; Characteristics of the Chlamydia genome; Level of the archaeal proteins.

Delineates a potential relationship between the NusG family, a bacterial transcription elongation factor, and RL24 ribosomal protein family. Report of a putative member of NusG with greater similarity to RL24 ribosomes; Identification of a novel motiff, called KOW after the authors' surname-initials; Structure and biochemical function of the KOW motiff.

The most common practice in studying and cataloguing prokaryotic diversity involves the grouping of sequences into operational taxonomic units (OTUs) at the 97% 16S rRNA gene sequence identity level, often using partial gene sequences, such as PCR-generated amplicons. Due to the high sequence conservation of rRNA genes, organisms belonging to closely related yet distinct species may be grouped under the same OTU. However, it remains unclear how much diversity has been underestimated by this practice. To address this question, we compared the OTUs of genomes defined at the 97% or 98.5% 16S rRNA gene identity level against OTUs of the same genomes defined at the 95% whole-genome average nucleotide identity (ANI), which is a much more accurate proxy for species. Our results show that OTUs resulting from a 98.5% 16S rRNA gene identity cutoff are more accurate than 97% compared to 95% ANI (90.5% versus 89.9% accuracy) but indistinguishable from any other threshold in the 98.29 to 98.78% range. Even with the more stringent thresholds, however, the 16S rRNA gene-based approach commonly underestimates the number of OTUs by ~12%, on average, compared to the ANI-based approach (~14% underestimation when using the 97% identity threshold). More importantly, the degree of underestimation can become 50% or more for certain taxa, such as the genera Pseudomonas, Burkholderia, Escherichia, Campylobacter, and Citrobacter. These results provide a quantitative view of the degree of underestimation of extant prokaryotic diversity by 16S rRNA gene-defined OTUs and suggest that genomic resolution is often necessary. [ABSTRACT FROM AUTHOR]

The type strain of the prospective Bradyrhizobium shewense sp. nov. ERR11T, was isolated from a nodule of the leguminous tree Erythrina brucei native to Ethiopia. The type strain Bradyrhizobium yuanmingense CCBAU 10071T, was isolated from the nodules of Lespedeza cuneata in Beijing, China. The genomes of ERR11T and CCBAU 10071T were sequenced by DOE-JGI and deposited at the DOE-JGI genome portal as well as at the European Nucleotide Archive. The genome of ERR11T is 9,163,226 bp in length and has 102 scaffolds, containing 8548 protein-coding and 86 RNA genes. The CCBAU 10071T genome is arranged in 108 scaffolds and consists of 8,201,522 bp long and 7776 protein-coding and 85 RNA genes. Both genomes contain symbiotic genes, which are homologous to the genes found in the complete genome sequence of Bradyrhizobium diazoefficiens USDA110T. The genes encoding for nodulation and nitrogen fixation in ERR11T showed high sequence similarity with homologous genes found in the draft genome of peanut-nodulating Bradyrhizobium arachidis LMG 26795T. The nodulation genes nolYAnodD2D1YABCSUIJ- nolO-nodZ of ERR11T and CCBAU 10071T are organized in a similar way to the homologous genes identified in the genomes of USDA110T, Bradyrhizobium ottawaense USDA 4 and Bradyrhizobium liaoningense CCBAU 05525. The genomes harbor hupSLCFHK and hypBFDE genes that code the expression of hydrogenase, an enzyme that helps rhizobia to uptake hydrogen released by the N2-fixation process and genes encoding denitrification functions napEDABC and norCBQD for nitrate and nitric oxide reduction, respectively. The genome of ERR11T also contains nosRZDFYLX genes encoding nitrous oxide reductase. Based on multilocus sequence analysis of housekeeping genes, the novel species, which contains eight strains formed a unique group close to the B. ottawaense branch. Genome Average Nucleotide Identity (ANI) calculated between the genome sequences of ERR11T and closely related sequences revealed that strains belonging to B. ottawaense branch (USDA4 and CCBAU15615), were the closest strains to the strain ERR11T with 95.2% ANI. Type strain ERR11T showed the highest DDH predicted value with CCBAU15615 (58.5%), followed by USDA 4 (53.1%). Nevertheless, the ANI and DDH values obtained between ERR11T and CCBAU 15615 or USDA 4 were below the cutoff values (ANI ≥ 96.5%; DDH ≥ 70%) for strains belonging to the same species, suggesting that ERR11T is a new species. Therefore, based on the phylogenetic analysis, ANI and DDH values, we formally propose the creation of B. shewense sp. nov. with strain ERR11T (HAMBI 3532T=LMG 30162T) as the type strain. [ABSTRACT FROM AUTHOR]

Rhizobium aethiopicum sp. nov. is a newly proposed species within the genus Rhizobium. This species includes six rhizobial strains; which were isolated from root nodules of the legume plant Phaseolus vulgaris growing in soils of Ethiopia. The species fixes nitrogen effectively in symbiosis with the host plant P. vulgaris, and is composed of aerobic, Gram-negative staining, rod-shaped bacteria. The genome of type strain HBR26T of R. aethiopicum sp. nov. was one of the rhizobial genomes sequenced as a part of the DOE JGI 2014 Genomic Encyclopedia project designed for soil and plant-associated and newly described type strains. The genome sequence is arranged in 62 scaffolds and consists of 6,557,588 bp length, with a 61% G + C content and 6221 protein-coding and 86 RNAs genes. The genome of HBR26T contains repABC genes (plasmid replication genes) homologous to the genes found in five different Rhizobium etli CFN42T plasmids, suggesting that HBR26T may have five additional replicons other than the chromosome. In the genome of HBR26T, the nodulation genes nodB, nodC, nodS, nodI, nodJ and nodD are located in the same module, and organized in a similar way as nod genes found in the genome of other known common bean-nodulating rhizobial species. nodA gene is found in a different scaffold, but it is also very similar to nodA genes of other bean-nodulating rhizobial strains. Though HBR26T is distinct on the phylogenetic tree and based on ANI analysis (the highest value 90.2% ANI with CFN42T) from other bean-nodulating species, these nod genes and most nitrogen-fixing genes found in the genome of HBR26T share high identity with the corresponding genes of known bean-nodulating rhizobial species (96-100% identity). This suggests that symbiotic genes might be shared between bean-nodulating rhizobia through horizontal gene transfer. R. aethiopicum sp. nov. was grouped into the genus Rhizobium but was distinct from all recognized species of that genus by phylogenetic analyses of combined sequences of the housekeeping genes recA and glnII. The closest reference type strains for HBR26T were R. etli CFN42T (94% similarity of the combined recA and glnII sequences) and Rhizobium bangladeshense BLR175T (93%). Genomic ANI calculation based on protein-coding genes also revealed that the closest reference strains were R. bangladeshense BLR175T and R. etli CFN42T with ANI values 91.8 and 90.2%, respectively. Nevertheless, the ANI values between HBR26T and BLR175T or CFN42T are far lower than the cutoff value of ANI (> = 96%) between strains in the same species, confirming that HBR26T belongs to a novel species. Thus, on the basis of phylogenetic, comparative genomic analyses and ANI results, we formally propose the creation of R. aethiopicum sp. nov. with strain HBR26T (=HAMBI 3550T=LMG 29711T) as the type strain. The genome assembly and annotation data is deposited in the DOE JGI portal and also available at European Nucleotide Archive under accession numbers FMAJ01000001-FMAJ01000062. [ABSTRACT FROM AUTHOR]

Members of the order Methanomicrobiales are abundant, and sometimes dominant, hydrogenotrophic (H2-CO2 utilizing) methanoarchaea in a broad range of anoxic habitats. Despite their key roles in greenhouse gas emissions and waste conversion to methane, little is known about the physiological and genomic bases for their widespread distribution and abundance. In this study, we compared the genomes of nine diverse Methanomicrobiales strains, examined their pangenomes, reconstructed gene flow and identified genes putatively mediating their success across different habitats. Most strains slowly increased gene content whereas one, Methanocorpusculum labreanum, evidenced genome downsizing. Peat-dwelling Methanomicrobiales showed adaptations centered on improved transport of scarce inorganic nutrients and likely use H+ rather than Na+ transmembrane chemiosmotic gradients during energy conservation. In contrast, other Methanomicrobiales show the potential to concurrently use Na+ and H+ chemiosmotic gradients. Analyses also revealed that the Methanomicrobiales lack a canonical electron bifurcation system (MvhABGD) known to produce low potential electrons in other orders of hydrogenotrophic methanogens. Additional putative differences in anabolic metabolism suggest that the dynamics of interspecies electron transfer from Methanomicrobiales syntrophic partners can also differ considerably. Altogether, these findings suggest profound differences in electron trafficking in the Methanomicrobiales compared with other hydrogenotrophs, and warrant further functional evaluations. [ABSTRACT FROM AUTHOR]

The substrates of the Brazilian campos rupestres, a grassland ecosystem, have extremely low concentrations of phosphorus and nitrogen, imposing restrictions to plant growth. Despite that, this ecosystem harbors almost 15% of the Brazilian plant diversity, raising the question of how plants acquire nutrients in such a harsh environment. Here, we set out to uncover the taxonomic profile, the compositional and functional differences and similarities, and the nutrient turnover potential of microbial communities associated with two plant species of the campos rupestres-dominant family Velloziaceae that grow over distinct substrates (soil and rock). Using amplicon sequencing data, we show that, despite the pronounced composition differentiation, the plant-associated soil and rock communities share a core of highly efficient colonizers that tend to be highly abundant and is enriched in 21 bacterial families. Functional investigation of metagenomes and 522 metagenome-assembled genomes revealed that the microorganisms found associated to plant roots are enriched in genes involved in organic compound intake, and phosphorus and nitrogen turnover. We show that potential for phosphorus transport, mineralization, and solubilization are mostly found within bacterial families of the shared microbiome, such as Xanthobacteraceae and Bryobacteraceae. We also detected the full repertoire of nitrogen cycle-related genes and discovered a lineage of Isosphaeraceae that acquired nitrogen-fixing potential via horizontal gene transfer and might be also involved in nitrification via a metabolic handoff association with Binataceae. We highlight that plant-associated microbial populations in the campos rupestres harbor a genetic repertoire with potential to increase nutrient availability and that the microbiomes of biodiversity hotspots can reveal novel mechanisms of nutrient turnover. [ABSTRACT FROM AUTHOR]

The power of next-generation sequencing has resulted in an explosive growth in the number of projects aiming to understand the metagenomic diversity of complex microbial environments. The interdisciplinary nature of this microbiome research community, along with the absence of reporting standards for microbiome data and samples, poses a significant challenge for follow-up studies. Commonly used names of metagenomes and metatranscriptomes in public databases currently lack the essential information necessary to accurately describe and classify the underlying samples, which makes a comparative analysis difficult to conduct and often results in misclassified sequences in data repositories. The Genomes OnLine Database (GOLD) (https:// gold.jgi.doe.gov/) at the Department of Energy Joint Genome Institute has been at the forefront of addressing this challenge by developing a standardized nomenclature system for naming microbiome samples. GOLD, currently in its twenty-fifth anniversary, continues to enrich the research community with hundreds of thousands of metagenomes and metatranscriptomes with well-curated and easy-to-understand names. Through this manuscript, we describe the overall naming process that can be easily adopted by researchers worldwide. Additionally, we propose the use of this naming system as a best practice for the scientific community to facilitate better interoperability and reusability of microbiome data. [ABSTRACT FROM AUTHOR]

Serratia plymuthica AS13 is a plant-associated Gammaproteobacteria, isolated from rapeseed roots. It is of special interest because of its ability to inhibit fungal pathogens of rapeseed and to promote plant growth. The complete genome of S. plymuthica AS13 consists of a 5,442,549 bp circular chromosome. The chromosome contains 4,951 protein-coding genes, 87 tRNA genes and 7 rRNA operons. This genome was sequenced as part of the project enti-tled "Genomics of four rapeseed plant growth promoting bacteria with antagonistic effect on plant pathogens" within the 2010 DOE-JGI Community Sequencing Program (CSP2010). [ABSTRACT FROM AUTHOR]

Assembling individual genomes from complex community metagenomic data remains a challenging issue for environmental studies. We evaluated the quality of genome assemblies from community short read data (Illumina 100 bp pair-ended sequences) using datasets recovered from freshwater and soil microbial communities as well as in silico simulations. Our analyses revealed that the genome of a single genotype (or species) can be accurately assembled from a complex metagenome when it shows at least about 20 × coverage. At lower coverage, however, the derived assemblies contained a substantial fraction of non-target sequences (chimeras), which explains, at least in part, the higher number of hypothetical genes recovered in metagenomic relative to genomic projects. We also provide examples of how to detect intrapopulation structure in metagenomic datasets and estimate the type and frequency of errors in assembled genes and contigs from datasets of varied species complexity. [ABSTRACT FROM AUTHOR]

Next-generation sequencing (NGS) is commonly used in metagenomic studies of complex microbial communities but whether or not different NGS platforms recover the same diversity from a sample and their assembled sequences are of comparable quality remain unclear. We compared the two most frequently used platforms, the Roche 454 FLX Titanium and the Illumina Genome Analyzer (GA) II, on the same DNA sample obtained from a complex freshwater planktonic community. Despite the substantial differences in read length and sequencing protocols, the platforms provided a comparable view of the community sampled. For instance, derived assemblies overlapped in ∼90% of their total sequences and in situ abundances of genes and genotypes (estimated based on sequence coverage) correlated highly between the two platforms (R2.0.9). Evaluation of base-call error, frameshift frequency, and contig length suggested that Illumina offered equivalent, if not better, assemblies than Roche 454. The results from metagenomic samples were further validated against DNA samples of eighteen isolate genomes, which showed a range of genome sizes and G+C% content. We also provide quantitative estimates of the errors in gene and contig sequences assembled from datasets characterized by different levels of complexity and G+C% content. For instance, we noted that homopolymer-associated, single-base errors affected ∼1% of the protein sequences recovered in Illumina contigs of 10 × coverage and 50% G+C; this frequency increased to ∼3% when non-homopolymer errors were also considered. Collectively, our results should serve as a useful practical guide for choosing proper sampling strategies and data possessing protocols for future metagenomic studies. [ABSTRACT FROM AUTHOR]

The production of low-cost biofuels in engineered microorganisms is of great interest due to the continual increase in the world's energy demands. Biodiesel is a renewable fuel that can potentially be produced in microbes cost-effectively. Fatty acid methyl esters (FAMEs) are a common component of biodiesel and can be synthesized from either triacylglycerol or free fatty acids (FFAs). Here we report the identification of a novel bacterial fatty acid methyltransferase (FAMT) that catalyzes the formation of FAMEs and 3-hydroxyl fatty acid methyl esters (3-OH-FAMEs) from the respective free acids and S-adenosylmethionine (AdoMet). FAMT exhibits a higher specificity toward 3-hydroxy free fatty acids (3-OH-FFAs) than FFAs, synthesizing 3-hydroxy fatty acid methyl esters (3-OH-FAMEs) in vivo. We have also identified bacterial members of the fatty acyl-acyl carrier protein (ACP) thioesterase (FAT) enzyme family with distinct acyl chain specificities. These bacterial FATs exhibit increased specificity toward 3-hydroxyacyl-ACP, generating 3-OH-FFAs, which can subsequently be utilized by FAMTs to produce 3-OH-FAMEs. PhaG (3-hydroxyacyl ACP:coenzyme A [CoA] transacylase) constitutes an alternative route to 3-OH-FFA synthesis; the coexpression of PhaG with FAMT led to the highest level of accumulation of 3-OH-FAMEs and FAMEs. The availability of AdoMet, the second substrate for FAMT, is an important factor regulating the amount of methyl esters produced by bacterial cells. Our results indicate that the deletion of the global methionine regulator metJ and the overexpression of methionine adenosyltransferase result in increased methyl ester synthesis. [ABSTRACT FROM AUTHOR]

With advances in DNA sequencing and miniaturized molecular biology workflows, rapid and affordable sequencing of single-cell genomes has become a reality. Compared to 16S rRNA gene surveys and shotgun metagenomics, large-scale application of single-cell genomics to whole microbial communities provides an integrated snapshot of community composition and function, directly links mobile elements to their hosts, and enables analysis of population heterogeneity of the dominant community members. To that end, we sequenced nearly 500 single-cell genomes from a low diversity hot spring sediment sample from Dewar Creek, British Columbia, and compared this approach to 16S rRNA gene amplicon and shotgun metagenomics applied to the same sample. We found that the broad taxonomic profiles were similar across the three sequencing approaches, though several lineages were missing from the 16S rRNA gene amplicon dataset, likely the result of primer mismatches. At the functional level, we detected a large array of mobile genetic elements present in the single-cell genomes but absent from the corresponding same species metagenome-assembled genomes. Moreover, we performed a single-cell population genomic analysis of the three most abundant community members, revealing differences in population structure based on mutation and recombination profiles. While the average pairwise nucleotide identities were similar across the dominant species-level lineages, we observed differences in the extent of recombination between these dominant populations. Most intriguingly, the creek's Hydrogenobacter sp. population appeared to be so recombinogenic that it more closely resembled a sexual species than a clonally evolving microbe. Together, this work demonstrates that a randomized single-cell approach can be useful for the exploration of previously uncultivated microbes from community composition to population structure. [ABSTRACT FROM AUTHOR]

Purpose: Studying the legume nodule microbiome is important for understanding the development and nutrition of the plants inhabited by the various microbes within and upon them. We analyzed the microbiomes of these underground organs from both an important crop plant (Medicago sativa) and a related legume (M. polymorpha) using metagenomic and culture-based techniques to identify the main cultivatable contributors to plant growth enhancement. Methods: Using high-throughput sequencing, culturing, and in planta techniques, we identified and analyzed a broad population of the bacterial taxa within Medicago nodules and the surrounding soil. Results: Fifty-one distinct bacterial strains were isolated and characterized from nodules of both Medicago species and their growth-promoting activities were studied. Sequencing of 16S rRNA gene amplicons showed that in addition to Ensifer, the dominant genus, a large number of Gram-positive bacteria belonging to the Firmicutes and Actinobacteria were also present. After performing ecological and plant growth-promoting trait analyses, selecting the most promising strains, and then performing in planta assays, we found that strains of Bacillus and Micromonospora among others could play important roles in supporting the growth, health, and productivity of the host plant. Conclusion: To our knowledge, the comparison of the biodiversity of the microbiota of undomesticated vs. cultivated Medicago roots and nodules is novel and shows the range of potential Plant Growth-Promoting Bacteria that could be used for plants of agricultural interest. These and other nodule-isolated microbes could also serve as inoculants with rhizobia with the goal of replacing synthetic fertilizers and pesticides for sustainable agriculture. [ABSTRACT FROM AUTHOR]

Reports a study that introduces a method to identifies gene-fusion events in complete genomes, solely based on sequence comparison. Why other biological, biochemical and genetic methods are tedious, labor intensive and inaccurate; Benefits of predicting that two proteins will interact based on structure or sequence alone; How researchers predicted functional associations of proteins in 64 unique fusion events.

Motivation Two key steps in the analysis of uncultured viruses recovered from metagenomes are the taxonomic classification of the viral sequences and the identification of putative host(s). Both steps rely mainly on the assignment of viral proteins to orthologs in cultivated viruses. Viral Protein Families (VPFs) can be used for the robust identification of new viral sequences in large metagenomics datasets. Despite the importance of VPF information for viral discovery, VPFs have not yet been explored for determining viral taxonomy and host targets. Results In this work, we classified the set of VPFs from the IMG/VR database and developed VPF-Class. VPF-Class is a tool that automates the taxonomic classification and host prediction of viral contigs based on the assignment of their proteins to a set of classified VPFs. Applying VPF-Class on 731K uncultivated virus contigs from the IMG/VR database, we were able to classify 363K contigs at the genus level and predict the host of over 461K contigs. In the RefSeq database, VPF-class reported an accuracy of nearly 100% to classify dsDNA, ssDNA and retroviruses, at the genus level, considering a membership ratio and a confidence score of 0.2. The accuracy in host prediction was 86.4%, also at the genus level, considering a membership ratio of 0.3 and a confidence score of 0.5. And, in the prophages dataset, the accuracy in host prediction was 86% considering a membership ratio of 0.6 and a confidence score of 0.8. Moreover, from the Global Ocean Virome dataset, over 817K viral contigs out of 1 million were classified. Availability and implementation The implementation of VPF-Class can be downloaded from https://github.com/biocom-uib/vpf-tools. Supplementary information Supplementary data are available at Bioinformatics online. [ABSTRACT FROM AUTHOR]

Millions of new viral sequences have been identified from metagenomes, but the quality and completeness of these sequences vary considerably. Here we present CheckV, an automated pipeline for identifying closed viral genomes, estimating the completeness of genome fragments and removing flanking host regions from integrated proviruses. CheckV estimates completeness by comparing sequences with a large database of complete viral genomes, including 76,262 identified from a systematic search of publicly available metagenomes, metatranscriptomes and metaviromes. After validation on mock datasets and comparison to existing methods, we applied CheckV to large and diverse collections of metagenome-assembled viral sequences, including IMG/VR and the Global Ocean Virome. This revealed 44,652 high-quality viral genomes (that is, >90% complete), although the vast majority of sequences were small fragments, which highlights the challenge of assembling viral genomes from short-read metagenomes. Additionally, we found that removal of host contamination substantially improved the accurate identification of auxiliary metabolic genes and interpretation of viral-encoded functions. The quality of viral genomes assembled from metagenome data is assessed by CheckV. [ABSTRACT FROM AUTHOR]

The article offers information on Synergistetes, a bacterial phylum. It notes that this phylum has 11 genera and one uncultured genus and is comprised of 300 16S ribosomal RNA gene sequences. It notes that the Synergistetes has compact genomes and these genomes share 35% and 50% of their genes. It adds that the phylum has been found in many anaerobic habitats, such as animal and insect gastrointestinal tracts and wastewater treatment systems, and is associated to periodontal disease. It notes that despite the habitat diversity of Synergistetes, its isolates demonstrate physiological resemblance to one another and are anaerobic. It mentions that although its genomes have the genes that could be used for lipopolysaccharide biosynthesis, they lack genes for the TolAQR-Pal complex.

Comprehensive, high-quality reference genomes are required for functional characterization and taxonomic assignment of the human gut microbiota. We present the Unified Human Gastrointestinal Genome (UHGG) collection, comprising 204,938 nonredundant genomes from 4,644 gut prokaryotes. These genomes encode >170 million protein sequences, which we collated in the Unified Human Gastrointestinal Protein (UHGP) catalog. The UHGP more than doubles the number of gut proteins in comparison to those present in the Integrated Gene Catalog. More than 70% of the UHGG species lack cultured representatives, and 40% of the UHGP lack functional annotations. Intraspecies genomic variation analyses revealed a large reservoir of accessory genes and single-nucleotide variants, many of which are specific to individual human populations. The UHGG and UHGP collections will enable studies linking genotypes to phenotypes in the human gut microbiome. More than 200,000 gut prokaryotic reference genomes and the proteins they encode are collated, providing comprehensive resources for microbiome researchers. [ABSTRACT FROM AUTHOR]

Our current knowledge of host–virus interactions in biofilms is limited to computational predictions based on laboratory experiments with a small number of cultured bacteria. However, natural biofilms are diverse and chiefly composed of uncultured bacteria and archaea with no viral infection patterns and lifestyle predictions described to date. Herein, we predict the first DNA sequence-based host–virus interactions in a natural biofilm. Using single-cell genomics and metagenomics applied to a hot spring mat of the Cone Pool in Mono County, California, we provide insights into virus–host range, lifestyle and distribution across different mat layers. Thirty-four out of 130 single cells contained at least one viral contig (26%), which, together with the metagenome-assembled genomes, resulted in detection of 59 viruses linked to 34 host species. Analysis of single-cell amplification kinetics revealed a lack of active viral replication on the single-cell level. These findings were further supported by mapping metagenomic reads from different mat layers to the obtained host–virus pairs, which indicated a low copy number of viral genomes compared to their hosts. Lastly, the metagenomic data revealed high layer specificity of viruses, suggesting limited diffusion to other mat layers. Taken together, these observations indicate that in low mobility environments with high microbial abundance, lysogeny is the predominant viral lifestyle, in line with the previously proposed "Piggyback-the-Winner" theory. [ABSTRACT FROM AUTHOR]

Bacteria and archaea possessing the hgcAB gene pair methylate inorganic mercury (Hg) to form highly toxic methylmercury. HgcA consists of a corrinoid binding domain and a transmembrane domain, and HgcB is a dicluster ferredoxin. However, their detailed structure and function have not been thoroughly characterized. We modeled the HgcAB complex by combining metagenome sequence data mining, coevolution analysis, and Rosetta structure calculations. In addition, we overexpressed HgcA and HgcB in Escherichia coli, confirmed spectroscopically that they bind cobalamin and [4Fe-4S] clusters, respectively, and incorporated these cofactors into the structural model. Surprisingly, the two domains of HgcA do not interact with each other, but HgcB forms extensive contacts with both domains. The model suggests that conserved cysteines in HgcB are involved in shuttling HgII, methylmercury, or both. These findings refine our understanding of the mechanism of Hg methylation and expand the known repertoire of corrinoid methyltransferases in nature. Connor J. Cooper et al. expressed HgcA and HgcB in Escherichia coli and modeled the structure of the HgcAB complex by combining metagenome sequence data, coevolution analysis, and ab initio structure calculations. This study provides insights into the biochemical mechanism of mercury (Hg) methylation. [ABSTRACT FROM AUTHOR]

Kiloniella laminariae is a true marine bacterium and the first member of the family and order, the Kiloniellaceae and Kiloniellales. K. laminariae LD81T (= DSM 19542T) was isolated from the marine macroalga Saccharina latissima and is a mesophilic, typical marine chemoheterotrophic aerobic bacterium with antifungal activity. Phylogenetic analysis of 16S rRNA gene sequence revealed the similarity of K. laminariae LD81T not only with three validly described species of the genus Kiloniella, but also with undescribed isolates and clone sequences from marine samples in the range of 93.6–96.7%. We report on the analysis of the draft genome of this alphaproteobacterium and describe some selected features. The 4.4 Mb genome has a G + C content of 51.4%, contains 4213 coding sequences including 51 RNA genes as well as 4162 protein-coding genes, and is a part of the Genomic Encyclopaedia of Bacteria and Archaea (GEBA) project. The genome provides insights into a number of metabolic properties, such as carbon and sulfur metabolism, and indicates the potential for denitrification and the biosynthesis of secondary metabolites. Comparative genome analysis was performed with K. laminariae LD81T and the animal-associated species Kiloniella majae M56.1T from a spider crab, Kiloniella spongiae MEBiC09566T from a sponge as well as Kiloniella litopenai P1-1 from a white shrimp, which all inhabit quite different marine habitats. The analysis revealed that the K. laminariae LD81T contains 1397 unique genes, more than twice the amount of the other species. Unique among others is a mixed PKS/NRPS biosynthetic gene cluster with similarity to the biosynthetic gene cluster responsible for the production of syringomycin. [ABSTRACT FROM AUTHOR]

Background: Cupriavidus strain STM 6070 was isolated from nickel-rich soil collected near Koniambo massif, New Caledonia, using the invasive legume trap host Mimosa pudica. STM 6070 is a heavy metal-tolerant strain that is highly effective at fixing nitrogen with M. pudica. Here we have provided an updated taxonomy for STM 6070 and described salient features of the annotated genome, focusing on heavy metal resistance (HMR) loci and heavy metal efflux (HME) systems. Results: The 6,771,773 bp high-quality-draft genome consists of 107 scaffolds containing 6118 protein-coding genes. ANI values show that STM 6070 is a new species of Cupriavidus. The STM 6070 symbiotic region was syntenic with that of the M. pudica-nodulating Cupriavidus taiwanensis LMG 19424T. In contrast to the nickel and zinc sensitivity of C. taiwanensis strains, STM 6070 grew at high Ni2+ and Zn2+ concentrations. The STM 6070 genome contains 55 genes, located in 12 clusters, that encode HMR structural proteins belonging to the RND, MFS, CHR, ARC3, CDF and P-ATPase protein superfamilies. These HMR molecular determinants are putatively involved in arsenic (ars), chromium (chr), cobalt-zinc-cadmium (czc), copper (cop, cup), nickel (nie and nre), and silver and/or copper (sil) resistance. Seven of these HMR clusters were common to symbiotic and non-symbiotic Cupriavidus species, while four clusters were specific to STM 6070, with three of these being associated with insertion sequences. Within the specific STM 6070 HMR clusters, three novel HME-RND systems (nieIC cep nieBA, czcC2B2A2, and hmxB zneAC zneR hmxS) were identified, which constitute new candidate genes for nickel and zinc resistance. Conclusions: STM 6070 belongs to a new Cupriavidus species, for which we have proposed the name Cupriavidus neocaledonicus sp. nov.. STM6070 harbours a pSym with a high degree of gene conservation to the pSyms of M. pudica-nodulating C. taiwanensis strains, probably as a result of recent horizontal transfer. The presence of specific HMR clusters, associated with transposase genes, suggests that the selection pressure of the New Caledonian ultramafic soils has driven the specific adaptation of STM 6070 to heavy-metal-rich soils via horizontal gene transfer. [ABSTRACT FROM AUTHOR]

Background: The genus Trichococcus currently contains nine species: T. flocculiformis, T. pasteurii, T. palustris, T. collinsii, T. patagoniensis, T. ilyis, T. paludicola, T. alkaliphilus, and T. shcherbakoviae. In general, Trichococcus species can degrade a wide range of carbohydrates. However, only T. pasteurii and a non-characterized strain of Trichococcus, strain ES5, have the capacity of converting glycerol to mainly 1,3-propanediol. Comparative genomic analysis of Trichococcus species provides the opportunity to further explore the physiological potential and uncover novel properties of this genus. Results: In this study, a genotype-phenotype comparative analysis of Trichococcus strains was performed. The genome of Trichococcus strain ES5 was sequenced and included in the comparison with the other nine type strains. Genes encoding functions related to e.g. the utilization of different carbon sources (glycerol, arabinan and alginate), antibiotic resistance, tolerance to low temperature and osmoregulation could be identified in all the sequences analysed. T. pasteurii and Trichococcus strain ES5 contain a operon with genes encoding necessary enzymes for 1,3-PDO production from glycerol. All the analysed genomes comprise genes encoding for cold shock domains, but only five of the Trichococcus species can grow at 0 °C. Protein domains associated to osmoregulation mechanisms are encoded in the genomes of all Trichococcus species, except in T. palustris, which had a lower resistance to salinity than the other nine studied Trichococcus strains. Conclusions: Genome analysis and comparison of ten Trichococcus strains allowed the identification of physiological traits related to substrate utilization and environmental stress resistance (e.g. to cold and salinity). Some substrates were used by single species, e.g. alginate by T. collinsii and arabinan by T. alkaliphilus. Strain ES5 may represent a subspecies of Trichococcus flocculiformis and contrary to the type strain (DSM 2094T), is able to grow on glycerol with the production of 1,3-propanediol. [ABSTRACT FROM AUTHOR]

Burkholderia cenocepacia TAtl-371 was isolated from the rhizosphere of a tomato plant growing in Atlatlahucan, Morelos, Mexico. This strain exhibited a broad antimicrobial spectrum against bacteria, yeast, and fungi. Here, we report and describe the improved, high-quality permanent draft genome of B. cenocepacia TAtl-371, which was sequenced using a combination of PacBio RS and PacBio RS II sequencing methods. The 7,496,106 bp genome of the TAtl-371 strain is arranged in three scaffolds, contains 6722 protein-coding genes, and 99 RNA only-encoding genes. Genome analysis revealed genes related to biosynthesis of antimicrobials such as non-ribosomal peptides, siderophores, chitinases, and bacteriocins. Moreover, analysis of bacterial growth on different carbon and nitrogen sources shows that the strain retains its antimicrobial ability. [ABSTRACT FROM AUTHOR]

CRISPR-Cas systems provide microbes with adaptive immunity to infectious nucleic acids and are widely employed as genome editing tools. These tools use RNA-guided Cas proteins whose large size (950 to 1400 amino acids) has been considered essential to their specific DNA- or RNA-targeting activities. Here we present a set of CRISPR-Cas systems from uncultivated archaea that contain Cas14, a family of exceptionally compact RNAguided nucleases (400 to 700 amino acids). Despite their small size, Cas14 proteins are capable of targeted single-stranded DNA (ssDNA) cleavage without restrictive sequence requirements. Moreover, target recognition by Cas14 triggers nonspecific cutting of ssDNA molecules, an activity that enables high-fidelity single-nucleotide polymorphism genotyping (Cas14-DETECTR). Metagenomic data show that multiple CRISPR-Cas14 systems evolved independently and suggest a potential evolutionary origin of singleeffector CRISPR-based adaptive immunity. [ABSTRACT FROM AUTHOR]

Summary: Pseudomonas is a large and diverse genus of Gammaproteobacteria. To provide a framework for discovery of evolutionary and taxonomic relationships of these bacteria, we compared the genomes of type strains of 163 species and 3 additional subspecies of Pseudomonas, including 118 genomes sequenced herein. A maximum likelihood phylogeny of the 166 type strains based on protein sequences of 100 single‐copy orthologous genes revealed thirteen groups of Pseudomonas, composed of two to sixty three species each. Pairwise average nucleotide identities and alignment fractions were calculated for the data set of the 166 type strains and 1224 genomes of Pseudomonas available in public databases. Results revealed that 394 of the 1224 genomes were distinct from any type strain, suggesting that the type strains represent only a fraction of the genomic diversity of the genus. The core genome of Pseudomonas was determined to contain 794 genes conferring primarily housekeeping functions. The results of this study provide a phylogenetic framework for future studies aiming to resolve the classification and phylogenetic relationships, identify new gene functions and phenotypes, and explore the ecological and metabolic potential of the Pseudomonas spp. [ABSTRACT FROM AUTHOR]

Productivity of ruminant livestock depends on the rumen microbiota, which ferment indigestible plant polysaccharides into nutrients used for growth. Understanding the functions carried out by the rumen microbiota is important for reducing greenhouse gas production by ruminants and for developing biofuels from lignocellulose. We present 410 cultured bacteria and archaea, together with their reference genomes, representing every cultivated rumen-associated archaeal and bacterial family. We evaluate polysaccharide degradation, short-chain fatty acid production and methanogenesis pathways, and assign specific taxa to functions. A total of 336 organisms were present in available rumen metagenomic data sets, and 134 were present in human gut microbiome data sets. Comparison with the human microbiome revealed rumen-specific enrichment for genes encoding de novo synthesis of vitamin B12, ongoing evolution by gene loss and potential vertical inheritance of the rumen microbiome based on underrepresentation of markers of environmental stress. We estimate that our Hungate genome resource represents ∼75% of the genus-level bacterial and archaeal taxa present in the rumen. [ABSTRACT FROM AUTHOR]

A bacterium belonging to the phylum Synergistetes, genus Dethiosulfovibrio was isolated in 2007 from a saline spring in Colombia. Dethiosulfovibrio salsuginis USBA 82T (DSM 21565T = KCTC 5659T) is a mesophilic, strictly anaerobic, slightly halophilic, Gram negative bacterium with a diderm cell envelope. The strain ferments peptides, amino acids and a few organic acids. Here we present the description of the complete genome sequencing and annotation of the type species Dethiosulfovibrio salsuginis USBA 82T. The genome consisted of 2.68 Mbp with a 53. 7% G + C. A total of 2609 genes were predicted and of those, 2543 were protein coding genes and 66 were RNA genes. We detected in USBA 82T genome six Synergistetes conserved signature indels (CSIs), specific for Jonquetella, Pyramidobacter and Dethiosulfovibrio. The genome of D. salsuginis contained, as expected, genes related to amino acid transport, amino acid metabolism and thiosulfate reduction. These genes represent the major gene groups of Synergistetes, related with their phenotypic traits, and interestingly, 11.8% of the genes in the genome belonged to the amino acid fermentation COG category. In addition, we identified in the genome some ammonification genes such as nitrate reductase genes. The presence of proline operon genes could be related to de novo synthesis of proline to protect the cell in response to high osmolarity. Our bioinformatics workflow included antiSMASH and BAGEL3 which allowed us to identify bacteriocins genes in the genome. [ABSTRACT FROM AUTHOR]

Paraburkholderia caballeronis is a plant-associated bacterium. Strain TNe-841T was isolated from the rhizosphere of tomato (Solanum lycopersicum L. var. lycopersicum) growing in Nepantla Mexico State. Initially this bacterium was found to effectively nodulate Phaseolus vulgaris L. However, from an analysis of the genome of strain TNe-841T and from repeat inoculation experiments, we found that this strain did not nodulate bean and also lacked nodulation genes, suggesting that the genes were lost. The genome consists of 7,115,141 bp with a G + C content of 67.01%. The sequence includes 6251 protein-coding genes and 87 RNA genes. [ABSTRACT FROM AUTHOR]

Here we present the physiological features of Pseudomonas extremaustralis strain USBA-GBX-515 (CMPUJU 515), isolated from soils in Superparamo ecosystems, > 4000 m.a.s.l, in the northern Andes of South America, as well as the thorough analysis of the draft genome. Strain USBA-GBX-515 is a Gram-negative rod shaped bacterium of 1.0-3. 0 μm x 0.5-1 μm, motile and unable to form spores, it grows aerobically and cells show one single flagellum. Several genetic indices, the phylogenetic analysis of the 16S rRNA gene sequence and the phenotypic characterization confirmed that USBA-GBX-515 is a member of Pseudomonas genus and, the similarity of the 16S rDNA sequence was 100% with P. extremaustralis strain CT14-3T. The draft genome of P. extremaustralis strain USBA-GBX-515 consisted of 6,143,638 Mb with a G + C content of 60.9 mol%. A total of 5665 genes were predicted and of those, 5544 were protein coding genes and 121 were RNA genes. The distribution of genes into COG functional categories showed that most genes were classified in the category of amino acid transport and metabolism (10.5%) followed by transcription (8.4%) and signal transduction mechanisms (7.3%). We performed experimental analyses of the lipolytic activity and results showed activity mainly on short chain fatty acids. The genome analysis demonstrated the existence of two genes, lip515A and est515A, related to a triacylglycerol lipase and carboxylesterase, respectively. Ammonification genes were also observed, mainly nitrate reductase genes. Genes related with synthesis of poly-hydroxyalkanoates (PHAs), especially poly-hydroxybutyrates (PHBs), were detected. The phaABC and phbABC operons also appeared complete in the genome. P. extremaustralis strain USBA-GBX-515 conserves the same gene organization of the type strain CT14-3T. We also thoroughly analyzed the potential for production of secondary metabolites finding close to 400 genes in 32 biosynthetic gene clusters involved in their production. [ABSTRACT FROM AUTHOR]

The permanent draft genome sequence of Actinotignum schaalii DSM 15541T is presented. The annotated genome includes 2,130,987 bp, with 1777 protein-coding and 58 rRNA-coding genes. Genome sequence analysis revealed absence of genes encoding for: components of the PTS systems, enzymes of the TCA cycle, glyoxylate shunt and gluconeogensis. Genomic data revealed that A. schaalii is able to oxidize carbohydrates via glycolysis, the nonoxidative pentose phosphate and the Entner-Doudoroff pathways. Besides, the genome harbors genes encoding for enzymes involved in the conversion of pyruvate to lactate, acetate and ethanol, which are found to be the end products of carbohydrate fermentation. The genome contained the gene encoding Type I fatty acid synthase required for de novo FAS biosynthesis. The plsY and plsX genes encoding the acyltransferases necessary for phosphatidic acid biosynthesis were absent from the genome. The genome harbors genes encoding enzymes responsible for isoprene biosynthesis via the mevalonate (MVA) pathway. Genes encoding enzymes that confer resistance to reactive oxygen species (ROS) were identified. In addition, A. schaalii harbors genes that protect the genome against viral infections. These include restriction-modification (RM) systems, type II toxin-antitoxin (TA), CRISPR-Cas and abortive infection system. A. schaalii genome also encodes several virulence factors that contribute to adhesion and internalization of this pathogen such as the tad genes encoding proteins required for pili assembly, the nanI gene encoding exo-alpha-sialidase, genes encoding heat shock proteins and genes encoding type VII secretion system. These features are consistent with anaerobic and pathogenic lifestyles. Finally, resistance to ciprofloxacin occurs by mutation in chromosomal genes that encode the subunits of DNA-gyrase (GyrA) and topisomerase IV (ParC) enzymes, while resistant to metronidazole was due to the frxA gene, which encodes NADPH-flavin oxidoreductase. [ABSTRACT FROM AUTHOR]

Dinitrogen (N2)-fixation by cyanobacteria in symbiosis with feathermosses is the primary pathway of biological nitrogen (N) input into boreal forests. Despite its significance, little is known about the cyanobacterial gene repertoire and regulatory rewiring needed for the establishment and maintenance of the symbiosis. To determine gene acquisitions and regulatory changes allowing cyanobacteria to form and maintain this symbiosis, we compared genomically closely related symbiotic-competent and -incompetent Nostoc strains using a proteogenomics approach and an experimental set up allowing for controlled chemical and physical contact between partners. Thirty-two gene families were found only in the genomes of symbiotic strains, including some never before associated with cyanobacterial symbiosis. We identified conserved orthologs that were differentially expressed in symbiotic strains, including protein families involved in chemotaxis and motility, NO regulation, sulfate/phosphate transport, and glycosyl-modifying and oxidative stress-mediating exoenzymes. The physical moss-cyanobacteria epiphytic symbiosis is distinct from other cyanobacteria-plant symbioses, with Nostoc retaining motility, and lacking modulation of N2-fixation, photosynthesis, GS-GOGAT cycle and heterocyst formation. The results expand our knowledge base of plant-cyanobacterial symbioses, provide a model of information and material exchange in this ecologically significant symbiosis, and suggest new currencies, namely nitric oxide and aliphatic sulfonates, may be involved in establishing and maintaining the cyanobacteria-feathermoss symbiosis. [ABSTRACT FROM AUTHOR]