Draft genomes of 12 Bifidobacterium isolates from human IBD fecal samples

Twelve Bifidobacterium strains were isolated from fecal samples of inflammatory bowel disease patients and matched "household control" individuals. These include the species Bifidobacterium adolescentis, Bifidobacterium animalis, Bifidobacterium breve, Bifidobacterium catenulatum, Bifidobacterium longum, and Bifidobacterium pseudocatenulatum.

Phables: from fragmented assemblies to high-quality bacteriophage genomes

MOTIVATION

Microbial communities have a profound impact on both human health and various environments. Viruses infecting bacteria, known as bacteriophages or phages, play a key role in modulating bacterial communities within environments. High-quality phage genome sequences are essential for advancing our understanding of phage biology, enabling comparative genomics studies and developing phage-based diagnostic tools. Most available viral identification tools consider individual sequences to determine whether they are of viral origin. As a result of challenges in viral assembly, fragmentation of genomes can occur, and existing tools may recover incomplete genome fragments. Therefore, the identification and characterization of novel phage genomes remain a challenge, leading to the need of improved approaches for phage genome recovery.

RESULTS

We introduce Phables, a new computational method to resolve phage genomes from fragmented viral metagenome assemblies. Phables identifies phage-like components in the assembly graph, models each component as a flow network, and uses graph algorithms and flow decomposition techniques to identify genomic paths. Experimental results of viral metagenomic samples obtained from different environments show that Phables recovers on average over 49% more high-quality phage genomes compared to existing viral identification tools. Furthermore, Phables can resolve variant phage genomes with over 99% average nucleotide identity, a distinction that existing tools are unable to make.

AVAILABILITY AND IMPLEMENTATION

Phables is available on GitHub at https://github.com/Vini2/phables.

Phables: from fragmented assemblies to high-quality bacteriophage genomes

Microbial communities influence both human health and different environments. Viruses infecting bacteria, known as bacteriophages or phages, play a key role in modulating bacterial communities within environments. High-quality phage genome sequences are essential for advancing our understanding of phage biology, enabling comparative genomics studies, and developing phage-based diagnostic tools. Most available viral identification tools consider individual sequences to determine whether they are of viral origin. As a result of the challenges in viral assembly, fragmentation of genomes can occur, leading to the need for new approaches in viral identification. Therefore, the identification and characterisation of novel phages remain a challenge. We introduce Phables, a new computational method to resolve phage genomes from fragmented viral metagenome assemblies. Phables identifies phage-like components in the assembly graph, models each component as a flow network, and uses graph algorithms and flow decomposition techniques to identify genomic paths. Experimental results of viral metagenomic samples obtained from different environments show that Phables recovers on average over 49% more high-quality phage genomes compared to existing viral identification tools. Furthermore, Phables can resolve variant phage genomes with over 99% average nucleotide identity, a distinction that existing tools are unable to make. Phables is available on GitHub at https://github.com/Vini2/phables.

Host interactions of novel Crassvirales species belonging to multiple families infecting bacterial host, Bacteroides cellulosilyticus WH2

Bacteroides, the prominent bacteria in the human gut, play a crucial role in degrading complex polysaccharides. Their abundance is influenced by phages belonging to the Crassvirales order. Despite identifying over 600 Crassvirales genomes computationally, only few have been successfully isolated. Continued efforts in isolation of more Crassvirales genomes can provide insights into phage-host-evolution and infection mechanisms. We focused on wastewater samples, as potential sources of phages infecting various Bacteroides hosts. Sequencing, assembly, and characterization of isolated phages revealed 14 complete genomes belonging to three novel Crassvirales species infecting Bacteroides cellulosilyticus WH2. These species, Kehishuvirus sp. 'tikkala' strain Bc01, Kolpuevirus sp. 'frurule' strain Bc03, and 'Rudgehvirus jaberico' strain Bc11, spanned two families, and three genera, displaying a broad range of virion productions. Upon testing all successfully cultured Crassvirales species and their respective bacterial hosts, we discovered that they do not exhibit co-evolutionary patterns with their bacterial hosts. Furthermore, we observed variations in gene similarity, with greater shared similarity observed within genera. However, despite belonging to different genera, the three novel species shared a unique structural gene that encodes the tail spike protein. When investigating the relationship between this gene and host interaction, we discovered evidence of purifying selection, indicating its functional importance. Moreover, our analysis demonstrated that this tail spike protein binds to the TonB-dependent receptors present on the bacterial host surface. Combining these observations, our findings provide insights into phage-host interactions and present three Crassvirales species as an ideal system for controlled infectivity experiments on one of the most dominant members of the human enteric virome.

Emergent community architecture despite distinct diversity in the global whale shark (Rhincodon typus) epidermal microbiome

Microbiomes confer beneficial physiological traits to their host, but microbial diversity is inherently variable, challenging the relationship between microbes and their contribution to host health. Here, we compare the diversity and architectural complexity of the epidermal microbiome from 74 individual whale sharks (Rhincodon typus) across five aggregations globally to determine if network properties may be more indicative of the microbiome-host relationship. On the premise that microbes are expected to exhibit biogeographic patterns globally and that distantly related microbial groups can perform similar functions, we hypothesized that microbiome co-occurrence patterns would occur independently of diversity trends and that keystone microbes would vary across locations. We found that whale shark aggregation was the most important factor in discriminating taxonomic diversity patterns. Further, microbiome network architecture was similar across all aggregations, with degree distributions matching Erdos-Renyi-type networks. The microbiome-derived networks, however, display modularity indicating a definitive microbiome structure on the epidermis of whale sharks. In addition, whale sharks hosted 35 high-quality metagenome assembled genomes (MAGs) of which 25 were present from all sample locations, termed the abundant 'core'. Two main MAG groups formed, defined here as Ecogroup 1 and 2, based on the number of genes present in metabolic pathways, suggesting there are at least two important metabolic niches within the whale shark microbiome. Therefore, while variability in microbiome diversity is high, network structure and core taxa are inherent characteristics of the epidermal microbiome in whale sharks. We suggest the host-microbiome and microbe-microbe interactions that drive the self-assembly of the microbiome help support a functionally redundant abundant core and that network characteristics should be considered when linking microbiomes with host health.

Host interactions of novel Crassvirales species belonging to multiple families infecting bacterial host, Bacteroides cellulosilyticus WH2

Bacteroides, the prominent bacteria in the human gut, play a crucial role in degrading complex polysaccharides. Their abundance is influenced by phages belonging to the Crassvirales order. Despite identifying over 600 Crassvirales genomes computationally, only few have been successfully isolated. Continued efforts in isolation of more Crassvirales genomes can provide insights into phage-host-evolution and infection mechanisms. We focused on wastewater samples, as potential sources of phages infecting various Bacteroides hosts. Sequencing, assembly, and characterization of isolated phages revealed 14 complete genomes belonging to three novel Crassvirales species infecting Bacteroides cellulosilyticus WH2. These species, Kehishuvirus sp. 'tikkala' strain Bc01, Kolpuevirus sp. 'frurule' strain Bc03, and 'Rudgehvirus jaberico' strain Bc11, spanned two families, and three genera, displaying a broad range of virion productions. Upon testing all successfully cultured Crassvirales species and their respective bacterial hosts, we discovered that they do not exhibit co-evolutionary patterns with their bacterial hosts. Furthermore, we observed variations in gene similarity, with greater shared similarity observed within genera. However, despite belonging to different genera, the three novel species shared a unique structural gene that encodes the tail spike protein. When investigating the relationship between this gene and host interaction, we discovered evidence of purifying selection, indicating its functional importance. Moreover, our analysis demonstrated that this tail spike protein binds to the TonB-dependent receptors present on the bacterial host surface. Combining these observations, our findings provide insights into phage-host interactions and present three Crassvirales species as an ideal system for controlled infectivity experiments on one of the most dominant members of the human enteric virome.

Impact statement

Bacteriophages play a crucial role in shaping microbial communities within the human gut. Among the most dominant bacteriophages in the human gut microbiome are Crassvirales phages, which infect Bacteroides. Despite being widely distributed, only a few Crassvirales genomes have been isolated, leading to a limited understanding of their biology, ecology, and evolution. This study isolated and characterized three novel Crassvirales genomes belonging to two different families, and three genera, but infecting one bacterial host, Bacteroides cellulosilyticus WH2. Notably, the observation confirmed the phages are not co-evolving with their bacterial hosts, rather have a shared ability to exploit similar features in their bacterial host. Additionally, the identification of a critical viral protein undergoing purifying selection and interacting with the bacterial receptors opens doors to targeted therapies against bacterial infections. Given Bacteroides role in polysaccharide degradation in the human gut, our findings advance our understanding of the phage-host interactions and could have important implications for the development of phage-based therapies. These discoveries may hold implications for improving gut health and metabolism to support overall well-being.

Data summary

The genomes used in this research are available on Sequence Read Archive (SRA) within the project, PRJNA737576. Bacteroides cellulosilyticus WH2, Kehishuvirus sp. 'tikkala' strain Bc01, Kolpuevirus sp. ' frurule' strain Bc03, and 'Rudgehvirus jaberico' strain Bc11 are all available on GenBank with accessions NZ_CP072251.1 ( B. cellulosilyticus WH2), QQ198717 (Bc01), QQ198718 (Bc03), and QQ198719 (Bc11), and we are working on making the strains available through ATCC. The 3D protein structures for the three Crassvirales genomes are available to download at doi.org/10.25451/flinders.21946034.

The Promise and Pitfalls of Prophages

Phages dominate every ecosystem on the planet. While virulent phages sculpt the microbiome by killing their bacterial hosts, temperate phages provide unique growth advantages to their hosts through lysogenic conversion. Many prophages benefit their host, and prophages are responsible for genotypic and phenotypic differences that separate individual microbial strains. However, the microbes also endure a cost to maintain those phages: additional DNA to replicate and proteins to transcribe and translate. We have never quantified those benefits and costs. Here, we analysed over two and a half million prophages from over half a million bacterial genome assemblies. Analysis of the whole dataset and a representative subset of taxonomically diverse bacterial genomes demonstrated that the normalised prophage density was uniform across all bacterial genomes above 2 Mbp. We identified a constant carrying capacity of phage DNA per bacterial DNA. We estimated that each prophage provides cellular services equivalent to approximately 2.4 % of the cell's energy or 0.9 ATP per bp per hour. We demonstrate analytical, taxonomic, geographic, and temporal disparities in identifying prophages in bacterial genomes that provide novel targets for identifying new phages. We anticipate that the benefits bacteria accrue from the presence of prophages balance the energetics involved in supporting prophages. Furthermore, our data will provide a new framework for identifying phages in environmental datasets, diverse bacterial phyla, and from different locations.

The Epidermal Microbiome Within an Aggregation of Leopard Sharks (Triakis semifasciata) Has Taxonomic Flexibility with Gene Functional Stability Across Three Time-points

The epidermis of Chondrichthyan fishes consists of dermal denticles with production of minimal but protein-rich mucus that collectively, influence the attachment and biofilm development of microbes, facilitating a unique epidermal microbiome. Here, we use metagenomics to provide the taxonomic and functional characterization of the epidermal microbiome of the Triakis semifasciata (leopard shark) at three time-points collected across 4 years to identify links between microbial groups and host metabolism. Our aims include (1) describing the variation of microbiome taxa over time and identifying recurrent microbiome members (present across all time-points); (2) investigating the relationship between the recurrent and flexible taxa (those which are not found consistently across time-points); (3) describing the functional compositions of the microbiome which may suggest links with the host metabolism; and (4) identifying whether metabolic processes are shared across microbial genera or are unique to specific taxa. Microbial members of the microbiome showed high similarity between all individuals (Bray-Curtis similarity index = 82.7, where 0 = no overlap, 100 = total overlap) with the relative abundance of those members varying across sampling time-points, suggesting flexibility of taxa in the microbiome. One hundred and eighty-eight genera were identified as recurrent, including Pseudomonas, Erythrobacter, Alcanivorax, Marinobacter, and Sphingopxis being consistently abundant across time-points, while Limnobacter and Xyella exhibited switching patterns with high relative abundance in 2013, Sphingobium and Sphingomona in 2015, and Altermonas, Leeuwenhoekiella, Gramella, and Maribacter in 2017. Of the 188 genera identified as recurrent, the top 19 relatively abundant genera formed three recurrent groups. The microbiome also displayed high functional similarity between individuals (Bray-Curtis similarity index = 97.6) with gene function composition remaining consistent across all time-points. These results show that while the presence of microbial genera exhibits consistency across time-points, their abundances do fluctuate. Microbial functions however remain stable across time-points; thus, we suggest the leopard shark microbiomes exhibit functional redundancy. We show coexistence of microbes hosted in elasmobranch microbiomes that encode genes involved in utilizing nitrogen, but not fixing nitrogen, degrading urea, and resistant to heavy metal.

Stingray epidermal microbiomes are species-specific with local adaptations

Marine host-associated microbiomes are affected by a combination of species-specific (e.g., host ancestry, genotype) and habitat-specific features (e.g., environmental physiochemistry and microbial biogeography). The stingray epidermis provides a gradient of characteristics from high dermal denticles coverage with low mucus to reduce dermal denticles and high levels of mucus. Here we investigate the effects of host phylogeny and habitat by comparing the epidermal microbiomes of Myliobatis californica (bat rays) with a mucus rich epidermis, and Urobatis halleri (round rays) with a mucus reduced epidermis from two locations, Los Angeles and San Diego, California (a 150 km distance). We found that host microbiomes are species-specific and distinct from the water column, however composition of M. californica microbiomes showed more variability between individuals compared to U. halleri. The variability in the microbiome of M. californica caused the microbial taxa to be similar across locations, while U. halleri microbiomes were distinct across locations. Despite taxonomic differences, Shannon diversity is the same across the two locations in U. halleri microbiomes suggesting the taxonomic composition are locally adapted, but diversity is maintained by the host. Myliobatis californica and U. halleri microbiomes maintain functional similarity across Los Angeles and San Diego and each ray showed several unique functional genes. Myliobatis californica has a greater relative abundance of RNA Polymerase III-like genes in the microbiome than U. halleri, suggesting specific adaptations to a heavy mucus environment. Construction of Metagenome Assembled Genomes (MAGs) identified novel microbial species within Rhodobacteraceae, Moraxellaceae, Caulobacteraceae, Alcanivoracaceae and Gammaproteobacteria. All MAGs had a high abundance of active RNA processing genes, heavy metal, and antibiotic resistant genes, suggesting the stingray mucus supports high microbial growth rates, which may drive high levels of competition within the microbiomes increasing the antimicrobial properties of the microbes.

Coral and Seawater Metagenomes Reveal Key Microbial Functions to Coral Health and Ecosystem Functioning Shaped at Reef Scale

The coral holobiont is comprised of a highly diverse microbial community that provides key services to corals such as protection against pathogens and nutrient cycling. The coral surface mucus layer (SML) microbiome is very sensitive to external changes, as it constitutes the direct interface between the coral host and the environment. Here, we investigate whether the bacterial taxonomic and functional profiles in the coral SML are shaped by the local reef zone and explore their role in coral health and ecosystem functioning. The analysis was conducted using metagenomes and metagenome-assembled genomes (MAGs) associated with the coral Pseudodiploria strigosa and the water column from two naturally distinct reef environments in Bermuda: inner patch reefs exposed to a fluctuating thermal regime and the more stable outer reefs. The microbial community structure in the coral SML varied according to the local environment, both at taxonomic and functional levels. The coral SML microbiome from inner reefs provides more gene functions that are involved in nutrient cycling (e.g., photosynthesis, phosphorus metabolism, sulfur assimilation) and those that are related to higher levels of microbial activity, competition, and stress response. In contrast, the coral SML microbiome from outer reefs contained genes indicative of a carbohydrate-rich mucus composition found in corals exposed to less stressful temperatures and showed high proportions of microbial gene functions that play a potential role in coral disease, such as degradation of lignin-derived compounds and sulfur oxidation. The fluctuating environment in the inner patch reefs of Bermuda could be driving a more beneficial coral SML microbiome, potentially increasing holobiont resilience to environmental changes and disease.

There Is More than Multiple Choice: Crowd-Sourced Assessment Tips for Online, Hybrid, and Face-to-Face Environments

Developing effective assessments of student learning is a challenging task for faculty and even more difficult for those in emerging disciplines that lack readily available resources and standards. With the power of technology-enhanced education and accessible digital learning platforms, instructors are also looking for assessments that work in an online format. This article will be useful for all teachers, but especially for entry-level instructors, in addition to more mature instructors who are looking to become more well versed in assessment, who seek a succinct summary of assessment types to springboard the integration of new forms of assessment of student learning into their courses. In this paper, ten assessment types, all appropriate for face-to-face, blended, and online modalities, are discussed. The assessments are mapped to a set of bioinformatics core competencies with examples of how they have been used to assess student learning. Although bioinformatics is used as the focus of the assessment types, the question types are relevant to many disciplines.

Elasmobranch microbiomes: emerging patterns and implications for host health and ecology

Elasmobranchs (sharks, skates and rays) are of broad ecological, economic, and societal value. These globally important fishes are experiencing sharp population declines as a result of human activity in the oceans. Research to understand elasmobranch ecology and conservation is critical and has now begun to explore the role of body-associated microbiomes in shaping elasmobranch health. Here, we review the burgeoning efforts to understand elasmobranch microbiomes, highlighting microbiome variation among gastrointestinal, oral, skin, and blood-associated niches. We identify major bacterial lineages in the microbiome, challenges to the field, key unanswered questions, and avenues for future work. We argue for prioritizing research to determine how microbiomes interact mechanistically with the unique physiology of elasmobranchs, potentially identifying roles in host immunity, disease, nutrition, and waste processing. Understanding elasmobranch–microbiome interactions is critical for predicting how sharks and rays respond to a changing ocean and for managing healthy populations in managed care.

Mitochondrial genome of the Smoothnose wedgefish Rhynchobatus laevis from the Western Indian Ocean

We present the first mitogenome sequence of the Smoothnose Wedgefish, Rhynchobatus laevis obtained through field sequencing on the MinION handheld sequencer. The mitochondrial genome of R. laevis is 16,560 bp in length and consisted of 13 protein-coding genes (PCGs), 22 tRNA genes, 2 rRNA genes, and a non-coding control region (D-loop). GC content was at 40.1%. The control region was 867 bp in length. Whole mitochondrial genome sequence of R. laevis will enable improved understanding of distribution, abundance, catch and trade rates of the Critically Endangered species.

Mitochondrial genome of the Silvertip shark, Carcharhinus albimarginatus, from the British Indian Ocean Territory

The Chagos archipelago in the British Indian Ocean Territory (BIOT) has been lacking in detailed genetic studies of its chondrichthyan populations. Chondrichthyes in Chagos continue to be endangered through illegal fishing operations, necessitating species distribution and abundance studies to facilitate urgent monitoring and conservation of the species. Here, we present a complete mitochondrial genome of the Silvertip Shark, Carcharhinus albimarginatus sampled in the Chagos archipelago. The mitochondrial genome of C. albimarginatus was 16,706 bp in length and consisted of 13 protein-coding genes, 22 tRNA genes, two rRNA genes, a replication origin and a D-loop region. GC content was at 38.7% and the control region was 1,065 bp in length. We expect that mitogenomes presented here will aid development of molecular assays for species distribution studies. Overall these studies will promote effective conservation of marine ecosystemes in the BIOT.

The skin microbiome of elasmobranchs follows phylosymbiosis, but in teleost fishes, the microbiomes converge

BACKGROUND

The vertebrate clade diverged into Chondrichthyes (sharks, rays, and chimeras) and Osteichthyes fishes (bony fishes) approximately 420 mya, with each group accumulating vast anatomical and physiological differences, including skin properties. The skin of Chondrichthyes fishes is covered in dermal denticles, whereas Osteichthyes fishes are covered in scales and are mucous rich. The divergence time among these two fish groups is hypothesized to result in predictable variation among symbionts. Here, using shotgun metagenomics, we test if patterns of diversity in the skin surface microbiome across the two fish clades match predictions made by phylosymbiosis theory. We hypothesize (1) the skin microbiome will be host and clade-specific, (2) evolutionary difference in elasmobranch and teleost will correspond with a concomitant increase in host-microbiome dissimilarity, and (3) the skin structure of the two groups will affect the taxonomic and functional composition of the microbiomes.

RESULTS

We show that the taxonomic and functional composition of the microbiomes is host-specific. Teleost fish had lower average microbiome within clade similarity compared to among clade comparison, but their composition is not different among clade in a null based model. Elasmobranch's average similarity within clade was not different than across clade and not different in a null based model of comparison. In the comparison of host distance with microbiome distance, we found that the taxonomic composition of the microbiome was related to host distance for the elasmobranchs, but not the teleost fishes. In comparison, the gene function composition was not related to the host-organism distance for elasmobranchs but was negatively correlated with host distance for teleost fishes.

CONCLUSION

Our results show the patterns of phylosymbiosis are not consistent across both fish clades, with the elasmobranchs showing phylosymbiosis, while the teleost fish are not. The discrepancy may be linked to alternative processes underpinning microbiome assemblage, including possible historical host-microbiome evolution of the elasmobranchs and convergent evolution in the teleost which filter specific microbial groups. Our comparison of the microbiomes among fishes represents an investigation into the microbial relationships of the oldest divergence of extant vertebrate hosts and reveals that microbial relationships are not consistent across evolutionary timescales. Video abstract.

Mitochondrial genome of the silky shark Carcharhinus falciformis from the British Indian Ocean Territory Marine Protected Area

We present the first mitochondrial genome of Carcharhinus falciformis from the Chagos Archipelago in the British Indian Ocean Territory (BIOT) Marine Protected Area (MPA). The mitochondrial genome of C. falciformis is 16,701 bp in length and consists of 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and a non-coding control region (D-loop). GC content was at 40.1%. The control region was 1063 bp in length. The complete mitogenome sequence of C. falciformis from the BIOT MPA will enable improved conservation measures of the CITES listed species through studies of species distribution, population abundance, fishing pressure and wildlife trade.

Complete mitochondrial genome of the gray reef shark, Carcharhinus amblyrhynchos (Carcharhiniformes: Carcharhinidae)

We report the first mitochondrial genome sequences for the gray reef shark, Carcharhinus amblyrhynchos. Two specimens from the British Indian Ocean Territory were sequenced independently using two different next generation sequencing methods, namely short read sequencing on the Illumina HiSeq and long read sequencing on the Oxford Nanopore Technologies’ MinION sequencer. The two sequences are 99.9% identical and are 16,705 base pairs (bp) and 16,706 bp in length. The mitogenome contains 22 tRNA genes, two rRNA genes, 13 protein-coding genes and two non-coding regions; the control region and the origin of light-strand replication (OL).

Modeling of the Coral Microbiome: the Influence of Temperature and Microbial Network

Host-associated microbial communities are shaped by extrinsic and intrinsic factors to the holobiont organism. Environmental factors and microbe-microbe interactions act simultaneously on the microbial community structure, making the microbiome dynamics challenging to predict. The coral microbiome is essential to the health of coral reefs and sensitive to environmental changes. Here, we develop a dynamic model to determine the microbial community structure associated with the surface mucus layer (SML) of corals using temperature as an extrinsic factor and microbial network as an intrinsic factor. The model was validated by comparing the predicted relative abundances of microbial taxa to the relative abundances of microbial taxa from the sample data. The SML microbiome from Pseudodiploria strigosa was collected across reef zones in Bermuda, where inner and outer reefs are exposed to distinct thermal profiles. A shotgun metagenomics approach was used to describe the taxonomic composition and the microbial network of the coral SML microbiome. By simulating the annual temperature fluctuations at each reef zone, the model output is statistically identical to the observed data. The model was further applied to six scenarios that combined different profiles of temperature and microbial network to investigate the influence of each of these two factors on the model accuracy. The SML microbiome was best predicted by model scenarios with the temperature profile that was closest to the local thermal environment, regardless of the microbial network profile. Our model shows that the SML microbiome of P. strigosa in Bermuda is primarily structured by seasonal fluctuations in temperature at a reef scale, while the microbial network is a secondary driver.IMPORTANCE Coral microbiome dysbiosis (i.e., shifts in the microbial community structure or complete loss of microbial symbionts) caused by environmental changes is a key player in the decline of coral health worldwide. Multiple factors in the water column and the surrounding biological community influence the dynamics of the coral microbiome. However, by including only temperature as an external factor, our model proved to be successful in describing the microbial community associated with the surface mucus layer (SML) of the coral P. strigosa The dynamic model developed and validated in this study is a potential tool to predict the coral microbiome under different temperature conditions.

Genomic and ecological attributes of marine bacteriophages encoding bacterial virulence genes

Background

Bacteriophages encode genes that modify bacterial functions during infection. The acquisition of phage-encoded virulence genes is a major mechanism for the rise of bacterial pathogens. In coral reefs, high bacterial density and lysogeny has been proposed to exacerbate reef decline through the transfer of phage-encoded virulence genes. However, the functions and distribution of these genes in phage virions on the reef remain unknown.

Results

Here, over 28,000 assembled viral genomes from the free viral community in Atlantic and Pacific Ocean coral reefs were queried against a curated database of virulence genes. The diversity of virulence genes encoded in the viral genomes was tested for relationships with host taxonomy and bacterial density in the environment. These analyses showed that bacterial density predicted the profile of virulence genes encoded by phages. The Shannon diversity of virulence-encoding phages was negatively related with bacterial density, leading to dominance of fewer genes at high bacterial abundances. A statistical learning analysis showed that reefs with high microbial density were enriched in viruses encoding genes enabling bacterial recognition and invasion of metazoan epithelium. Over 60% of phages could not have their hosts identified due to limitations of host prediction tools; for those which hosts were identified, host taxonomy was not an indicator of the presence of virulence genes.

Conclusions

This study described bacterial virulence factors encoded in the genomes of bacteriophages at the community level. The results showed that the increase in microbial densities that occurs during coral reef degradation is associated with a change in the genomic repertoire of bacteriophages, specifically in the diversity and distribution of bacterial virulence genes. This suggests that phages are implicated in the rise of pathogens in disturbed marine ecosystems.

Microbial abundance, composition, and function in nectar are shaped by flower visitor identity

Microbial dispersal is essential for establishment in new habitats, but the role of vector identity is poorly understood in community assembly and function. Here, we compared microbial assembly and function in floral nectar visited by legitimate pollinators (hummingbirds) and nectar robbers (carpenter bees). We assessed effects of visitation on the abundance and composition of culturable bacteria and fungi and their taxonomy and function using shotgun metagenomics and nectar chemistry. We also compared metagenome-assembled genomes (MAGs) of Acinetobacter, a common and highly abundant nectar bacterium, among visitor treatments. Visitation increased microbial abundance, but robbing resulted in 10× higher microbial abundance than pollination. Microbial communities differed among visitor treatments: robbed flowers were characterized by predominant nectar specialists within Acetobacteraceae and Metschnikowiaceae, with a concurrent loss of rare taxa, and these resulting communities harbored genes relating to osmotic stress, saccharide metabolism and specialized transporters. Gene differences were mirrored in function: robbed nectar contained a higher percentage of monosaccharides. Draft genomes of Acinetobacter revealed distinct amino acid and saccharide utilization pathways in strains isolated from robbed versus pollinated flowers. Our results suggest an unrecognized cost of nectar robbing for pollination and distinct effects of visitor type on interactions between plants and pollinators. Overall, these results suggest vector identity is an underappreciated factor structuring microbial community assembly and function.

PHANOTATE: a novel approach to gene identification in phage genomes

MOTIVATION

Currently there are no tools specifically designed for annotating genes in phages. Several tools are available that have been adapted to run on phage genomes, but due to their underlying design, they are unable to capture the full complexity of phage genomes. Phages have adapted their genomes to be extremely compact, having adjacent genes that overlap and genes completely inside of other longer genes. This non-delineated genome structure makes it difficult for gene prediction using the currently available gene annotators. Here we present PHANOTATE, a novel method for gene calling specifically designed for phage genomes. Although the compact nature of genes in phages is a problem for current gene annotators, we exploit this property by treating a phage genome as a network of paths: where open reading frames are favorable, and overlaps and gaps are less favorable, but still possible. We represent this network of connections as a weighted graph, and use dynamic programing to find the optimal path.

RESULTS

We compare PHANOTATE to other gene callers by annotating a set of 2133 complete phage genomes from GenBank, using PHANOTATE and the three most popular gene callers. We found that the four programs agree on 82% of the total predicted genes, with PHANOTATE predicting more genes than the other three. We searched for these extra genes in both GenBank's non-redundant protein database and all of the metagenomes in the sequence read archive, and found that they are present at levels that suggest that these are functional protein-coding genes.

AVAILABILITY AND IMPLEMENTATION

https://github.com/deprekate/PHANOTATE.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Global phylogeography and ancient evolution of the widespread human gut virus crAssphage

Microbiomes are vast communities of microorganisms and viruses that populate all natural ecosystems. Viruses have been considered to be the most variable component of microbiomes, as supported by virome surveys and examples of high genomic mosaicism. However, recent evidence suggests that the human gut virome is remarkably stable compared with that of other environments. Here, we investigate the origin, evolution and epidemiology of crAssphage, a widespread human gut virus. Through a global collaboration, we obtained DNA sequences of crAssphage from more than one-third of the world's countries and showed that the phylogeography of crAssphage is locally clustered within countries, cities and individuals. We also found fully colinear crAssphage-like genomes in both Old-World and New-World primates, suggesting that the association of crAssphage with primates may be millions of years old. Finally, by exploiting a large cohort of more than 1,000 individuals, we tested whether crAssphage is associated with bacterial taxonomic groups of the gut microbiome, diverse human health parameters and a wide range of dietary factors. We identified strong correlations with different clades of bacteria that are related to Bacteroidetes and weak associations with several diet categories, but no significant association with health or disease. We conclude that crAssphage is a benign cosmopolitan virus that may have coevolved with the human lineage and is an integral part of the normal human gut virome.

Global phylogeography and ancient evolution of the widespread human gut virus crAssphage

Microbiomes are vast communities of microorganisms and viruses that populate all natural ecosystems. Viruses have been considered to be the most variable component of microbiomes, as supported by virome surveys and examples of high genomic mosaicism. However, recent evidence suggests that the human gut virome is remarkably stable compared with that of other environments. Here, we investigate the origin, evolution and epidemiology of crAssphage, a widespread human gut virus. Through a global collaboration, we obtained DNA sequences of crAssphage from more than one-third of the world's countries and showed that the phylogeography of crAssphage is locally clustered within countries, cities and individuals. We also found fully colinear crAssphage-like genomes in both Old-World and New-World primates, suggesting that the association of crAssphage with primates may be millions of years old. Finally, by exploiting a large cohort of more than 1,000 individuals, we tested whether crAssphage is associated with bacterial taxonomic groups of the gut microbiome, diverse human health parameters and a wide range of dietary factors. We identified strong correlations with different clades of bacteria that are related to Bacteroidetes and weak associations with several diet categories, but no significant association with health or disease. We conclude that crAssphage is a benign cosmopolitan virus that may have coevolved with the human lineage and is an integral part of the normal human gut virome.

'Genome skimming' with the MinION hand-held sequencer identifies CITES-listed shark species in India's exports market

Chondrichthyes - sharks, rays, skates, and chimeras, are among the most threatened and data deficient vertebrate species. Global demand for shark and ray derived products, drives unregulated and exploitative fishing practices, which are in turn facilitated by the lack of ecological data required for effective conservation of these species. Here, we describe a Next Generation Sequencing method (using the MinION, a hand-held portable sequencing device from Oxford Nanopore Technologies), and analyses pipeline for molecular ecological studies in Chondrichthyes. Using this method, the complete mitochondrial genome and nuclear intergenic and protein-coding sequences were obtained by direct sequencing of genomic DNA obtained from shark fin tissue. Recovered loci include mitochondrial barcode sequences- Cytochrome oxidase I, NADH2, 16S rRNA and 12S rRNA- and nuclear genetic loci such as 5.8S rRNA, Internal Transcribed Spacer 2, and 28S rRNA regions, which are commonly used for taxonomic identification. Other loci recovered were the nuclear protein-coding genes for antithrombin or SerpinC, Immunoglobulin lambda light chain, Preprogehrelin, selenium binding protein 1(SBP1), Interleukin-1 beta (IL-1β) and Recombination-Activating Gene 1 (RAG1). The median coverage across all genetic loci was 20x and sequence accuracy was ≥99.8% compared to reference sequences. Analyses of the nuclear ITS2 region and the mitochondrial protein-encoding loci allowed accurate taxonomic identification of the shark specimen as Carcharhinus falciformis, a CITES Appendix II species. MinION sequencing provided 1,152,211 bp of new shark genome, increasing the number of sequenced shark genomes to five. Phylogenetic analyses using both mitochondrial and nuclear loci provided evidence that Prionace glauca is nested within Carcharhinus, suggesting the need for taxonomic reassignment of P. glauca. We increased genomic information about a shark species for ecological and population genetic studies, enabled accurate identification of the shark tissue for biodiversity indexing and resolved phylogenetic relationships among multiple taxa. The method was independent of amplification bias, and adaptable for field assessments of other Chondrichthyes and wildlife species in the future.

Rhodoliths holobionts in a changing ocean: host-microbes interactions mediate coralline algae resilience under ocean acidification

Background

Life in the ocean will increasingly have to contend with a complex matrix of concurrent shifts in environmental properties that impact their physiology and control their life histories. Rhodoliths are coralline red algae (Corallinales, Rhodophyta) that are photosynthesizers, calcifiers, and ecosystem engineers and therefore represent important targets for ocean acidification (OA) research. Here, we exposed live rhodoliths to near-future OA conditions to investigate responses in their photosynthetic capacity, calcium carbonate production, and associated microbiome using carbon uptake, decalcification assays, and whole genome shotgun sequencing metagenomic analysis, respectively. The results from our live rhodolith assays were compared to similar manipulations on dead rhodolith (calcareous skeleton) biofilms and water column microbial communities, thereby enabling the assessment of host-microbiome interaction under climate-driven environmental perturbations.

Results

Under high pCO₂ conditions, live rhodoliths exhibited positive physiological responses, i.e. increased photosynthetic activity, and no calcium carbonate biomass loss over time. Further, whereas the microbiome associated with live rhodoliths remained stable and resembled a healthy holobiont, the microbial community associated with the water column changed after exposure to elevated pCO₂.

Conclusions

Our results suggest that a tightly regulated microbial-host interaction, as evidenced by the stability of the rhodolith microbiome recorded here under OA-like conditions, is important for host resilience to environmental stress. This study extends the scarce comprehension of microbes associated with rhodolith beds and their reaction to increased pCO₂, providing a more comprehensive approach to OA studies by assessing the host holobiont.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5064-4) contains supplementary material, which is available to authorized users.

Elevated temperature drives kelp microbiome dysbiosis, while elevated carbon dioxide induces water microbiome disruption

Global climate change includes rising temperatures and increased pCO2 concentrations in the ocean, with potential deleterious impacts on marine organisms. In this case study we conducted a four-week climate change incubation experiment, and tested the independent and combined effects of increased temperature and partial pressure of carbon dioxide (pCO2), on the microbiomes of a foundation species, the giant kelp Macrocystis pyrifera, and the surrounding water column. The water and kelp microbiome responded differently to each of the climate stressors. In the water microbiome, each condition caused an increase in a distinct microbial order, whereas the kelp microbiome exhibited a reduction in the dominant kelp-associated order, Alteromondales. The water column microbiomes were most disrupted by elevated pCO2, with a 7.3 fold increase in Rhizobiales. The kelp microbiome was most influenced by elevated temperature and elevated temperature in combination with elevated pCO2. Kelp growth was negatively associated with elevated temperature, and the kelp microbiome showed a 5.3 fold increase Flavobacteriales and a 2.2 fold increase alginate degrading enzymes and sulfated polysaccharides. In contrast, kelp growth was positively associated with the combination of high temperature and high pCO2 'future conditions', with a 12.5 fold increase in Planctomycetales and 4.8 fold increase in Rhodobacteriales. Therefore, the water and kelp microbiomes acted as distinct communities, where the kelp was stabilizing the microbiome under changing pCO2 conditions, but lost control at high temperature. Under future conditions, a new equilibrium between the kelp and the microbiome was potentially reached, where the kelp grew rapidly and the commensal microbes responded to an increase in mucus production.

Optimizing and evaluating the reconstruction of Metagenome-assembled microbial genomes

Background

Microbiome/host interactions describe characteristics that affect the host's health. Shotgun metagenomics includes sequencing a random subset of the microbiome to analyze its taxonomic and metabolic potential. Reconstruction of DNA fragments into genomes from metagenomes (called metagenome-assembled genomes) assigns unknown fragments to taxa/function and facilitates discovery of novel organisms. Genome reconstruction incorporates sequence assembly and sorting of assembled sequences into bins, characteristic of a genome. However, the microbial community composition, including taxonomic and phylogenetic diversity may influence genome reconstruction. We determine the optimal reconstruction method for four microbiome projects that had variable sequencing platforms (IonTorrent and Illumina), diversity (high or low), and environment (coral reefs and kelp forests), using a set of parameters to select for optimal assembly and binning tools.

Methods

We tested the effects of the assembly and binning processes on population genome reconstruction using 105 marine metagenomes from 4 projects. Reconstructed genomes were obtained from each project using 3 assemblers (IDBA, MetaVelvet, and SPAdes) and 2 binning tools (GroopM and MetaBat). We assessed the efficiency of assemblers using statistics that including contig continuity and contig chimerism and the effectiveness of binning tools using genome completeness and taxonomic identification.

Results

We concluded that SPAdes, assembled more contigs (143,718 ± 124 contigs) of longer length (N50 = 1632 ± 108 bp), and incorporated the most sequences (sequences-assembled = 19.65%). The microbial richness and evenness were maintained across the assembly, suggesting low contig chimeras. SPAdes assembly was responsive to the biological and technological variations within the project, compared with other assemblers. Among binning tools, we conclude that MetaBat produced bins with less variation in GC content (average standard deviation: 1.49), low species richness (4.91 ± 0.66), and higher genome completeness (40.92 ± 1.75) across all projects. MetaBat extracted 115 bins from the 4 projects of which 66 bins were identified as reconstructed metagenome-assembled genomes with sequences belonging to a specific genus. We identified 13 novel genomes, some of which were 100% complete, but show low similarity to genomes within databases.

Conclusions

In conclusion, we present a set of biologically relevant parameters for evaluation to select for optimal assembly and binning tools. For the tools we tested, SPAdes assembler and MetaBat binning tools reconstructed quality metagenome-assembled genomes for the four projects. We also conclude that metagenomes from microbial communities that have high coverage of phylogenetically distinct, and low taxonomic diversity results in highest quality metagenome-assembled genomes.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4294-1) contains supplementary material, which is available to authorized users.

Microbial processes driving coral reef organic carbon flow

Coral reefs are one of the most productive ecosystems on the planet, with primary production rates compared to that of rain forests. Benthic organisms release 10-50% of their gross organic production as mucus that stimulates heterotrophic microbial metabolism in the water column. As a result, coral reef microbes grow up to 50 times faster than open ocean communities. Anthropogenic disturbances cause once coral-dominated reefs to become dominated by fleshy organisms, with several outcomes for trophic relationships. Here we review microbial processes implicated in organic carbon flux in coral reefs displaying species phase shifts. The first section presents microbial players and interactions within the coral holobiont that contribute to reef carbon flow. In the second section, we identify four ecosystem-level microbial features that directly respond to benthic species phase shifts: community composition, biomass, metabolism and viral predation. The third section discusses the significance of microbial consumption of benthic organic matter to reef trophic relationships. In the fourth section, we propose that the 'microbial phase shifts' discussed here are conducive to lower resilience, facilitating the transition to new degradation states in coral reefs.

Aura-biomes are present in the water layer above coral reef benthic macro-organisms

As coral reef habitats decline worldwide, some reefs are transitioning from coral- to algal-dominated benthos with the exact cause for this shift remaining elusive. Increases in the abundance of microbes in the water column has been correlated with an increase in coral disease and reduction in coral cover. Here we investigated how multiple reef organisms influence microbial communities in the surrounding water column. Our study consisted of a field assessment of microbial communities above replicate patches dominated by a single macro-organism. Metagenomes were constructed from 20 L of water above distinct macro-organisms, including (1) the coral Mussismilia braziliensis, (2) fleshy macroalgae (Stypopodium, Dictota and Canistrocarpus), (3) turf algae, and (4) the zoanthid Palythoa caribaeorum and were compared to the water microbes collected 3 m above the reef. Microbial genera and functional potential were annotated using MG-RAST and showed that the dominant benthic macro-organisms influence the taxa and functions of microbes in the water column surrounding them, developing a specific “aura-biome”. The coral aura-biome reflected the open water column, and was associated with Synechococcus and functions suggesting oligotrophic growth, while the fleshy macroalgae aura-biome was associated with Ruegeria, Pseudomonas, and microbial functions suggesting low oxygen conditions. The turf algae aura-biome was associated with Vibrio, Flavobacterium, and functions suggesting pathogenic activity, while zoanthids were associated with Alteromonas and functions suggesting a stressful environment. Because each benthic organism has a distinct aura-biome, a change in benthic cover will change the microbial community of the water, which may lead to either the stimulation or suppression of the recruitment of benthic organisms.

The skin microbiome of the common thresher shark (Alopias vulpinus) has low taxonomic and gene function β-diversity

The health of sharks, like all organisms, is linked to their microbiome. At the skin interface, sharks have dermal denticles that protrude above the mucus, which may affect the types of microbes that occur here. We characterized the microbiome from the skin of the common thresher shark (Alopias vulpinus) to investigate the structure and composition of the skin microbiome. On average 618 812 (80.9% ± S.D. 0.44%) reads per metagenomic library contained open reading frames; of those, between 7.6% and 12.8% matched known protein sequences. Genera distinguishing the A. vulpinus microbiome from the water column included, Pseudoalteromonas (12.8% ± 4.7 of sequences), Erythrobacter (5. 3% ± 0.5) and Idiomarina (4.2% ± 1.2) and distinguishing gene pathways included, cobalt, zinc and cadmium resistance (2.2% ± 0.1); iron acquisition (1.2% ± 0.1) and ton/tol transport (1.3% ± 0.08). Taxonomic community overlap (100 - dissimilarity index) was greater in the skin microbiome (77.6), relative to the water column microbiome (70.6) and a reference host-associated microbiome (algae: 71.5). We conclude the A. vulpinus skin microbiome is influenced by filtering processes, including biochemical and biophysical components of the shark skin and result in a structured microbiome.

Bacterial Community Associated with the Reef Coral Mussismilia braziliensis's Momentum Boundary Layer over a Diel Cycle

Corals display circadian physiological cycles, changing from autotrophy during the day to heterotrophy during the night. Such physiological transition offers distinct environments to the microbial community associated with corals: an oxygen-rich environment during daylight hours and an oxygen-depleted environment during the night. Most studies of coral reef microbes have been performed on samples taken during the day, representing a bias in the understanding of the composition and function of these communities. We hypothesized that coral circadian physiology alters the composition and function of microbial communities in reef boundary layers. Here, we analyzed microbial communities associated with the momentum boundary layer (MBL) of the Brazilian endemic reef coral Mussismilia braziliensis during a diurnal cycle, and compared them to the water column. We determined microbial abundance and nutrient concentration in samples taken within a few centimeters of the coral's surface every 6 h for 48 h, and sequenced microbial metagenomes from a subset of the samples. We found that dominant taxa and functions in the coral MBL community were stable over the time scale of our sampling, with no significant shifts between night and day samples. Interestingly, the two water column metagenomes sampled 1 m above the corals were also very similar to the MBL metagenomes. When all samples were analyzed together, nutrient concentration significantly explained 40% of the taxonomic dissimilarity among dominant genera in the community. Functional profiles were highly homogenous and not significantly predicted by any environmental variables measured. Our data indicated that water flow may overrule the effects of coral physiology in the MBL bacterial community, at the scale of centimeters, and suggested that sampling resolution at the scale of millimeters may be necessary to address diurnal variation in community composition.

Nearshore Pelagic Microbial Community Abundance Affects Recruitment Success of Giant Kelp, Macrocystis pyrifera

Marine microbes mediate key ecological processes in kelp forest ecosystems and interact with macroalgae. Pelagic and biofilm-associated microbes interact with macroalgal propagules at multiple stages of recruitment, yet these interactions have not been described for Macrocystis pyrifera. Here we investigate the influence of microbes from coastal environments on recruitment of giant kelp, M. pyrifera. Through repeated laboratory experiments, we tested the effects of altered pelagic microbial abundance on the settlement and development of the microscopic propagules of M. pyrifera during recruitment. M. pyrifera zoospores were reared in laboratory microcosms exposed to environmental microbial communities from seawater during the complete haploid stages of the kelp recruitment cycle, including zoospore release, followed by zoospore settlement, to gametophyte germination and development. We altered the microbial abundance states differentially in three independent experiments with repeated trials, where microbes were (a) present or absent in seawater, (b) altered in community composition, and (c) altered in abundance. Within the third experiment, we also tested the effect of nearshore versus offshore microbial communities on the macroalgal propagules. Distinct pelagic microbial communities were collected from two southern California temperate environments reflecting contrasting intensity of human influence, the nearshore Point Loma kelp forest and the offshore Santa Catalina Island kelp forest. The Point Loma kelp forest is a high impacted coastal region adjacent to the populous San Diego Bay; whereas the kelp forest at Catalina Island is a low impacted region of the Channel Islands, 40 km offshore the southern California coast, and is adjacent to a marine protected area. Kelp gametophytes reared with nearshore Point Loma microbes showed lower survival, growth, and deteriorated morphology compared to gametophytes with the offshore Catalina Island microbial community, and these effects were magnified under high microbial abundances. Reducing abundance of Point Loma microbes restored M. pyrifera propagule success. Yet an intermediate microbial abundance was optimal for kelp propagules reared with Catalina Island microbes, suggesting that microbes also have a beneficial influence on kelp. Our study shows that pelagic microbes from nearshore and offshore environments are differentially influencing kelp propagule success, which has significant implications for kelp recruitment and kelp forest ecosystem health.

Genomic Comparison of the Closely-Related Salmonella enterica Serovars Enteritidis, Dublin and Gallinarum

The Salmonella enterica serovars Enteritidis, Dublin, and Gallinarum are closely related but differ in virulence and host range. To identify the genetic elements responsible for these differences and to better understand how these serovars are evolving, we sequenced the genomes of Enteritidis strain LK5 and Dublin strain SARB12 and compared these genomes to the publicly available Enteritidis P125109, Dublin CT 02021853 and Dublin SD3246 genome sequences. We also compared the publicly available Gallinarum genome sequences from biotype Gallinarum 287/91 and Pullorum RKS5078. Using bioinformatic approaches, we identified single nucleotide polymorphisms, insertions, deletions, and differences in prophage and pseudogene content between strains belonging to the same serovar. Through our analysis we also identified several prophage cargo genes and pseudogenes that affect virulence and may contribute to a host-specific, systemic lifestyle. These results strongly argue that the Enteritidis, Dublin and Gallinarum serovars of Salmonella enterica evolve by acquiring new genes through horizontal gene transfer, followed by the formation of pseudogenes. The loss of genes necessary for a gastrointestinal lifestyle ultimately leads to a systemic lifestyle and niche exclusion in the host-specific serovars.

Sequencing at sea: challenges and experiences in Ion Torrent PGM sequencing during the 2013 Southern Line Islands Research Expedition

Genomics and metagenomics have revolutionized our understanding of marine microbial ecology and the importance of microbes in global geochemical cycles. However, the process of DNA sequencing has always been an abstract extension of the research expedition, completed once the samples were returned to the laboratory. During the 2013 Southern Line Islands Research Expedition, we started the first effort to bring next generation sequencing to some of the most remote locations on our planet. We successfully sequenced twenty six marine microbial genomes, and two marine microbial metagenomes using the Ion Torrent PGM platform on the Merchant Yacht Hanse Explorer. Onboard sequence assembly, annotation, and analysis enabled us to investigate the role of the microbes in the coral reef ecology of these islands and atolls. This analysis identified phosphonate as an important phosphorous source for microbes growing in the Line Islands and reinforced the importance of L-serine in marine microbial ecosystems. Sequencing in the field allowed us to propose hypotheses and conduct experiments and further sampling based on the sequences generated. By eliminating the delay between sampling and sequencing, we enhanced the productivity of the research expedition. By overcoming the hurdles associated with sequencing on a boat in the middle of the Pacific Ocean we proved the flexibility of the sequencing, annotation, and analysis pipelines.

A highly abundant bacteriophage discovered in the unknown sequences of human faecal metagenomes

Metagenomics, or sequencing of the genetic material from a complete microbial community, is a promising tool to discover novel microbes and viruses. Viral metagenomes typically contain many unknown sequences. Here we describe the discovery of a previously unidentified bacteriophage present in the majority of published human faecal metagenomes, which we refer to as crAssphage. Its ~97 kbp genome is six times more abundant in publicly available metagenomes than all other known phages together; it comprises up to 90% and 22% of all reads in virus-like particle (VLP)-derived metagenomes and total community metagenomes, respectively; and it totals 1.68% of all human faecal metagenomic sequencing reads in the public databases. The majority of crAssphage-encoded proteins match no known sequences in the database, which is why it was not detected before. Using a new co-occurrence profiling approach, we predict a Bacteroides host for this phage, consistent with Bacteroides-related protein homologues and a unique carbohydrate-binding domain encoded in the phage genome.

Metagenomic studies of microbial communities often report DNA sequences from unidentified viruses. Here, Dutilh et al. analyse metagenomic data to reveal the complete genome of an abundant, ubiquitous virus from human faeces, and predict that the virus infects bacteria of the Bacteroides group.

Combining de novo and reference-guided assembly with scaffold_builder

Genome sequencing has become routine, however genome assembly still remains a challenge despite the computational advances in the last decade. In particular, the abundance of repeat elements in genomes makes it difficult to assemble them into a single complete sequence. Identical repeats shorter than the average read length can generally be assembled without issue. However, longer repeats such as ribosomal RNA operons cannot be accurately assembled using existing tools. The application Scaffold_builder was designed to generate scaffolds – super contigs of sequences joined by N-bases – based on the similarity to a closely related reference sequence. This is independent of mate-pair information and can be used complementarily for genome assembly, e.g. when mate-pairs are not available or have already been exploited. Scaffold_builder was evaluated using simulated pyrosequencing reads of the bacterial genomes Escherichia coli 042, Lactobacillus salivarius UCC118 and Salmonella enterica subsp. enterica serovar Typhi str. P-stx-12. Moreover, we sequenced two genomes from Salmonella enterica serovar Typhimurium LT2 G455 and Salmonella enterica serovar Typhimurium SDT1291 and show that Scaffold_builder decreases the number of contig sequences by 53% while more than doubling their average length. Scaffold_builder is written in Python and is available at http://edwards.sdsu.edu/scaffold_builder. A web-based implementation is additionally provided to allow users to submit a reference genome and a set of contigs to be scaffolded.

Multivariate analysis of functional metagenomes

Metagenomics is a primary tool for the description of microbial and viral communities. The sheer magnitude of the data generated in each metagenome makes identifying key differences in the function and taxonomy between communities difficult to elucidate. Here we discuss the application of seven different data mining and statistical analyses by comparing and contrasting the metabolic functions of 212 microbial metagenomes within and between 10 environments. Not all approaches are appropriate for all questions, and researchers should decide which approach addresses their questions. This work demonstrated the use of each approach: for example, random forests provided a robust and enlightening description of both the clustering of metagenomes and the metabolic processes that were important in separating microbial communities from different environments. All analyses identified that the presence of phage genes within the microbial community was a predictor of whether the microbial community was host-associated or free-living. Several analyses identified the subtle differences that occur with environments, such as those seen in different regions of the marine environment.

Oxygen minimum zones harbour novel viral communities with low diversity

Oxygen minimum zones (OMZs) are oceanographic features that affect ocean productivity and biodiversity, and contribute to ocean nitrogen loss and greenhouse gas emissions. Here we describe the viral communities associated with the Eastern Tropical South Pacific (ETSP) OMZ off Iquique, Chile for the first time through abundance estimates and viral metagenomic analysis. The viral-to-microbial ratio (VMR) in the ETSP OMZ fluctuated in the oxycline and declined in the anoxic core to below one on several occasions. The number of viral genotypes (unique genomes as defined by sequence assembly) ranged from 2040 at the surface to 98 in the oxycline, which is the lowest viral diversity recorded to date in the ocean. Within the ETSP OMZ viromes, only 4.95% of genotypes were shared between surface and anoxic core viromes using reciprocal BLASTn sequence comparison. ETSP virome comparison with surface marine viromes (Sargasso Sea, Gulf of Mexico, Kingman Reef, Chesapeake Bay) revealed a dissimilarity of ETSP OMZ viruses to those from other oceanic regions. From the 1.4 million non-redundant DNA sequences sampled within the altered oxygen conditions of the ETSP OMZ, more than 97.8% were novel. Of the average 3.2% of sequences that showed similarity to the SEED non-redundant database, phage sequences dominated the surface viromes, eukaryotic virus sequences dominated the oxycline viromes, and phage sequences dominated the anoxic core viromes. The viral community of the ETSP OMZ was characterized by fluctuations in abundance, taxa and diversity across the oxygen gradient. The ecological significance of these changes was difficult to predict; however, it appears that the reduction in oxygen coincides with an increased shedding of eukaryotic viruses in the oxycline, and a shift to unique viral genotypes in the anoxic core.

Assembly rules of reef corals are flexible along a steep climatic gradient

Coral reefs, one of the world's most complex and vulnerable ecosystems, face an uncertain future in coming decades as they continue to respond to anthropogenic climate change, overfishing, pollution, and other human impacts [1, 2]. Traditionally, marine macroecology is based on presence/absence data from taxonomic checklists or geographic ranges, providing a qualitative overview of spatial shifts in species richness that treats rare and common species equally [3, 4]. As a consequence, regional and long-term shifts in relative abundances of individual taxa are poorly understood. Here we apply a more rigorous quantitative approach to examine large-scale spatial variation in the species composition and abundance of corals on midshelf reefs along the length of Australia's Great Barrier Reef, a biogeographic region where species richness is high and relatively homogeneous [5]. We demonstrate that important functional components of coral assemblages "sample" space differently at 132 sites separated by up to 1740 km, leading to complex latitudinal shifts in patterns of absolute and relative abundance. The flexibility in community composition that we document along latitudinal environmental gradients indicates that climate change is likely to result in a reassortment of coral reef taxa rather than wholesale loss of entire reef ecosystems.

Substrate type determines metagenomic profiles from diverse chemical habitats

Environmental parameters drive phenotypic and genotypic frequency variations in microbial communities and thus control the extent and structure of microbial diversity. We tested the extent to which microbial community composition changes are controlled by shifting physiochemical properties within a hypersaline lagoon. We sequenced four sediment metagenomes from the Coorong, South Australia from samples which varied in salinity by 99 Practical Salinity Units (PSU), an order of magnitude in ammonia concentration and two orders of magnitude in microbial abundance. Despite the marked divergence in environmental parameters observed between samples, hierarchical clustering of taxonomic and metabolic profiles of these metagenomes showed striking similarity between the samples (>89%). Comparison of these profiles to those derived from a wide variety of publically available datasets demonstrated that the Coorong sediment metagenomes were similar to other sediment, soil, biofilm and microbial mat samples regardless of salinity (>85% similarity). Overall, clustering of solid substrate and water metagenomes into discrete similarity groups based on functional potential indicated that the dichotomy between water and solid matrices is a fundamental determinant of community microbial metabolism that is not masked by salinity, nutrient concentration or microbial abundance.

The lagoon at Caroline/Millennium atoll, Republic of Kiribati: natural history of a nearly pristine ecosystem

A series of surveys were carried out to characterize the physical and biological parameters of the Millennium Atoll lagoon during a research expedition in April of 2009. Millennium is a remote coral atoll in the Central Pacific belonging to the Republic of Kiribati, and a member of the Southern Line Islands chain. The atoll is among the few remaining coral reef ecosystems that are relatively pristine. The lagoon is highly enclosed, and was characterized by reticulate patch and line reefs throughout the center of the lagoon as well as perimeter reefs around the rim of the atoll. The depth reached a maximum of 33.3 m in the central region of the lagoon, and averaged between 8.8 and 13.7 m in most of the pools. The deepest areas were found to harbor large platforms of Favia matthaii, which presumably provided a base upon which the dominant corals (Acropora spp.) grew to form the reticulate reef structure. The benthic algal communities consisted mainly of crustose coralline algae (CCA), microfilamentous turf algae and isolated patches of Halimeda spp. and Caulerpa spp. Fish species richness in the lagoon was half of that observed on the adjacent fore reef. The lagoon is likely an important nursery habitat for a number of important fisheries species including the blacktip reef shark and Napoleon wrasse, which are heavily exploited elsewhere around the world but were common in the lagoon at Millennium. The lagoon also supports an abundance of giant clams (Tridacna maxima). Millennium lagoon provides an excellent reference of a relatively undisturbed coral atoll. As with most coral reefs around the world, the lagoon communities of Millennium may be threatened by climate change and associated warming, acidification and sea level rise, as well as sporadic local resource exploitation which is difficult to monitor and enforce because of the atoll's remote location. While the remote nature of Millennium has allowed it to remain one of the few nearly pristine coral reef ecosystems in the world, it is imperative that this ecosystem receives protection so that it may survive for future generations.

The GAAS metagenomic tool and its estimations of viral and microbial average genome size in four major biomes

Metagenomic studies characterize both the composition and diversity of uncultured viral and microbial communities. BLAST-based comparisons have typically been used for such analyses; however, sampling biases, high percentages of unknown sequences, and the use of arbitrary thresholds to find significant similarities can decrease the accuracy and validity of estimates. Here, we present Genome relative Abundance and Average Size (GAAS), a complete software package that provides improved estimates of community composition and average genome length for metagenomes in both textual and graphical formats. GAAS implements a novel methodology to control for sampling bias via length normalization, to adjust for multiple BLAST similarities by similarity weighting, and to select significant similarities using relative alignment lengths. In benchmark tests, the GAAS method was robust to both high percentages of unknown sequences and to variations in metagenomic sequence read lengths. Re-analysis of the Sargasso Sea virome using GAAS indicated that standard methodologies for metagenomic analysis may dramatically underestimate the abundance and importance of organisms with small genomes in environmental systems. Using GAAS, we conducted a meta-analysis of microbial and viral average genome lengths in over 150 metagenomes from four biomes to determine whether genome lengths vary consistently between and within biomes, and between microbial and viral communities from the same environment. Significant differences between biomes and within aquatic sub-biomes (oceans, hypersaline systems, freshwater, and microbialites) suggested that average genome length is a fundamental property of environments driven by factors at the sub-biome level. The behavior of paired viral and microbial metagenomes from the same environment indicated that microbial and viral average genome sizes are independent of each other, but indicative of community responses to stressors and environmental conditions.

Comparative metagenomics reveals host specific metavirulomes and horizontal gene transfer elements in the chicken cecum microbiome

BACKGROUND

The complex microbiome of the ceca of chickens plays an important role in nutrient utilization, growth and well-being of these animals. Since we have a very limited understanding of the capabilities of most species present in the cecum, we investigated the role of the microbiome by comparative analyses of both the microbial community structure and functional gene content using random sample pyrosequencing. The overall goal of this study was to characterize the chicken cecal microbiome using a pathogen-free chicken and one that had been challenged with Campylobacter jejuni.

METHODOLOGY/PRINCIPAL FINDINGS

Comparative metagenomic pyrosequencing was used to generate 55,364,266 bases of random sampled pyrosequence data from two chicken cecal samples. SSU rDNA gene tags and environmental gene tags (EGTs) were identified using SEED subsystems-based annotations. The distribution of phylotypes and EGTs detected within each cecal sample were primarily from the Firmicutes, Bacteroidetes and Proteobacteria, consistent with previous SSU rDNA libraries of the chicken cecum. Carbohydrate metabolism and virulence genes are major components of the EGT content of both of these microbiomes. A comparison of the twelve major pathways in the SEED Virulence Subsystem (metavirulome) represented in the chicken cecum, mouse cecum and human fecal microbiomes showed that the metavirulomes differed between these microbiomes and the metavirulomes clustered by host environment. The chicken cecum microbiomes had the broadest range of EGTs within the SEED Conjugative Transposon Subsystem, however the mouse cecum microbiomes showed a greater abundance of EGTs in this subsystem. Gene assemblies (32 contigs) from one microbiome sample were predominately from the Bacteroidetes, and seven of these showed sequence similarity to transposases, whereas the remaining sequences were most similar to those from catabolic gene families.

CONCLUSION/SIGNIFICANCE

This analysis has demonstrated that mobile DNA elements are a major functional component of cecal microbiomes, thus contributing to horizontal gene transfer and functional microbiome evolution. Moreover, the metavirulomes of these microbiomes appear to associate by host environment. These data have implications for defining core and variable microbiome content in a host species. Furthermore, this suggests that the evolution of host specific metavirulomes is a contributing factor in disease resistance to zoonotic pathogens.

Microbial ecology of four coral atolls in the Northern Line Islands

Microbes are key players in both healthy and degraded coral reefs. A combination of metagenomics, microscopy, culturing, and water chemistry were used to characterize microbial communities on four coral atolls in the Northern Line Islands, central Pacific. Kingman, a small uninhabited atoll which lies most northerly in the chain, had microbial and water chemistry characteristic of an open ocean ecosystem. On this atoll the microbial community was equally divided between autotrophs (mostly Prochlorococcus spp.) and heterotrophs. In contrast, Kiritimati, a large and populated ( approximately 5500 people) atoll, which is most southerly in the chain, had microbial and water chemistry characteristic of a near-shore environment. On Kiritimati, there were 10 times more microbial cells and virus-like particles in the water column and these microbes were dominated by heterotrophs, including a large percentage of potential pathogens. Culturable Vibrios were common only on Kiritimati. The benthic community on Kiritimati had the highest prevalence of coral disease and lowest coral cover. The middle atolls, Palmyra and Tabuaeran, had intermediate densities of microbes and viruses and higher percentages of autotrophic microbes than either Kingman or Kiritimati. The differences in microbial communities across atolls could reflect variation in 1) oceaonographic and/or hydrographic conditions or 2) human impacts associated with land-use and fishing. The fact that historically Kingman and Kiritimati did not differ strongly in their fish or benthic communities (both had large numbers of sharks and high coral cover) suggest an anthropogenic component in the differences in the microbial communities. Kingman is one of the world's most pristine coral reefs, and this dataset should serve as a baseline for future studies of coral reef microbes. Obtaining the microbial data set, from atolls is particularly important given the association of microbes in the ongoing degradation of coral reef ecosystems worldwide.

Baselines and degradation of coral reefs in the Northern Line Islands

Effective conservation requires rigorous baselines of pristine conditions to assess the impacts of human activities and to evaluate the efficacy of management. Most coral reefs are moderately to severely degraded by local human activities such as fishing and pollution as well as global change, hence it is difficult to separate local from global effects. To this end, we surveyed coral reefs on uninhabited atolls in the northern Line Islands to provide a baseline of reef community structure, and on increasingly populated atolls to document changes associated with human activities. We found that top predators and reef-building organisms dominated unpopulated Kingman and Palmyra, while small planktivorous fishes and fleshy algae dominated the populated atolls of Tabuaeran and Kiritimati. Sharks and other top predators overwhelmed the fish assemblages on Kingman and Palmyra so that the biomass pyramid was inverted (top-heavy). In contrast, the biomass pyramid at Tabuaeran and Kiritimati exhibited the typical bottom-heavy pattern. Reefs without people exhibited less coral disease and greater coral recruitment relative to more inhabited reefs. Thus, protection from overfishing and pollution appears to increase the resilience of reef ecosystems to the effects of global warming.

Assessing anchor damage on coral reefs: a case study in selection of environmental indicators

Because environmental conservation can remove scarce natural resources from competing uses, it is important to gain support for conservation programs by demonstrating that management actions have been effective in achieving their goals. One way to do this is to show that selected significant environmental variables (indicators) vary between managed and unmanaged areas or change over time following implementation of a management regime. However, identifying indicators that reflect environmental conditions relevant to management practices has proven difficult. This paper focuses on developing a framework for choosing indicators in a coral reef habitat. The framework consisted of three phases: (1) information gathering to identify candidate variables; (2) field-testing candidate variables at sites that differ in intensity of human activity, thus identifying potential indicators; and (3) evaluating potential indicators against a set of feasibility criteria to identify the most useful indicators. To identify indicators suitable to measure the success of a management strategy to reduce anchor damage to a coral reef, 24 candidate variables were identified and evaluated at sites with different intensities of anchoring. In this study, measures that reflected injuries to coral colonies were generally more efficient than traditional measures of coral cover in describing the effects of anchoring. The number of overturned colonies was identified as the single most useful indicator of coral reef condition associated with anchoring intensities. The indicator selection framework developed here has the advantages of being transparent, cost efficient, and readily transferable to other types of human activities and management strategies.