Application of 'next-generation' sequencing technologies to microbial genetics

Transcriptomic and metabolomic insights into flavor variations in wild and cultivated Agaricus bisporus

Agaricus bisporus is a widely cultivated edible fungus globally. However, the mechanisms underlying the differences in flavor and nutritional traits between wild-type (W) and cultivated-type (C) strains remain unclear, which hinders the artificial breeding of high-quality varieties. This study systematically revealed, for the first time, the molecular and metabolic basis of flavor divergence between wild and cultivated A. bisporus by integrating transcriptomics and metabolomics. A total of 43 strains (23 wild-type and 20 cultivated-type) were analyzed using high-throughput sequencing and liquid chromatography-tandem mass spectrometry (LC-MS/MS) to dissect differences in gene expression and metabolite profiles. Results showed that although total protein and amino acid contents exhibited no significant differences, transcriptomic analysis identified significant upregulation of AGABI2DRAFT_188981 and AGABI2DRAFT_191000 (genes associated with high-affinity methionine permease MUP1) in cultivated strains, suggesting their indirect regulation of flavor formation via methionine metabolism. Metabolomic analysis further revealed a marked increase in uridine levels in cultivated strains (3.2-fold higher than wild-type, p < 0.01), indicating potential medicinal value, while wild strains were enriched with flavor precursors such as fumaric acid and isoleucine (fold change ≥ 2.5). In contrast, cultivated strains accumulated metabolites like 2-hydroxybutyric acid and α-ketoglutarate (VIP > 1.5). This study pioneered the construction of a gene-metabolite correlation network, identifying a strong positive correlation between AGABI2DRAFT_191352 (6-phosphofructokinase) and 2-hydroxybutyric acid (r = 0.82), highlighting the regulatory role of glycolytic flux in flavor metabolism. These findings not only elucidate the impact of artificial cultivation on metabolic reprogramming in A. bisporus but also provide critical molecular targets for targeted breeding of strains with enhanced flavor and nutritional value, offering practical significance for advancing the edible fungi industry.

Novel Enzymes From the Red Sea Brine Pools: Current State and Potential

The Red Sea is a marine environment with unique chemical characteristics and physical topographies. Among the various habitats offered by the Red Sea, the deep-sea brine pools are the most extreme in terms of salinity, temperature and metal contents. Nonetheless, the brine pools host rich polyextremophilic bacterial and archaeal communities. These microbial communities are promising sources for various classes of enzymes adapted to harsh environments - extremozymes. Extremozymes are emerging as novel biocatalysts for biotechnological applications due to their ability to perform catalytic reactions under harsh biophysical conditions, such as those used in many industrial processes. In this review, we provide an overview of the extremozymes from different Red Sea brine pools and discuss the overall biotechnological potential of the Red Sea proteome.

Characterisation of Early Positive mcr-1 Resistance Gene and Plasmidome in Escherichia coli Pathogenic Strains Associated with Variable Phylogroups under Colistin Selection

An antibiotic susceptibility monitoring programme was conducted from 2004 to 2010, resulting in a collection of 143 Escherichia coli cultured from bovine faecal samples (diarrhoea) and milk-aliquots (mastitis). The isolates were subjected to whole-genome sequencing and were distributed in phylogroups A, B1, B2, C, D, E, and G with no correlation for particular genotypes with pathotypes. In fact, the population structure showed that the strains belonging to the different phylogroups matched broadly to ST complexes; however, the isolates are randomly associated with the diseases, highlighting the necessity to investigate the virulence factors more accurately in order to identify the mechanisms by which they cause disease. The antimicrobial resistance was assessed phenotypically, confirming the genomic prediction on three isolates that were resistant to colistin, although one isolate was positive for the presence of the gene mcr-1 but susceptible to colistin. To further characterise the genomic context, the four strains were sequenced by using a single-molecule long read approach. Genetic analyses indicated that these four isolates harboured complex and diverse plasmids encoding not only antibiotic resistant genes (including mcr-1 and bla) but also virulence genes (siderophore, ColV, T4SS). A detailed description of the plasmids of these four E. coli strains, which are linked to bovine mastitis and diarrhoea, is presented for the first time along with the characterisation of the predicted antibiotic resistance genes. The study highlighted the diversity of incompatibility types encoding complex antibiotic resistance elements such as Tn6330, ISEcp1, Tn6029, and IS5075. The mcr-1 resistance determinant was identified in IncHI2 plasmids pCFS3273-1 and pCFS3292-1, thus providing some of the earliest examples of mcr-1 reported in Europe, and these sequences may be a representative of the early mcr-1 plasmidome characterisation in the EU/EEA.

Improvement of Surface-Enhanced Raman Scattering Method for Single Bacterial Cell Analysis

Surface-enhanced Raman scattering (SERS) is a useful tool for label-free analysis of bacteria at the single cell level. However, low reproducibility limits the use of SERS. In this study, for the sake of sensitive and reproducible Raman spectra, we optimized the methods for preparing silver nanoparticles (AgNPs) and depositing AgNPs onto a cell surface. We found that fast dropwise addition of AgNO₃ into the reductant produced smaller and more stable AgNPs, with an average diameter of 45 ± 4 nm. Compared with that observed after simply mixing the bacterial cells with AgNPs, the SERS signal was significantly improved after centrifugation. To optimize the SERS enhancement method, the centrifugal force, method for preparing AgNPs, concentration of AgNPs, ionic strength of the solution used to suspend the cells, and density of the cells were chosen as impact factors and optimized through orthogonal experiments. Finally, the improved method could generate sensitive and reproducible SERS spectra from single Escherichia coli cells, and the SERS signals primarily arose from the cell envelope. We further verified that this optimal method was feasible for the detection of low to 25% incorporation of ¹³C isotopes by the cells and the discrimination of different bacterial species. Our work provides an improved method for generating sensitive and reproducible SERS spectra.

Screening for Candidate Genes Associated with Biocontrol Mechanisms of Bacillus pumilus DX01 Using Tn5 Transposon Mutagenesis and a 2-DE-Based Comparative Proteomic Analysis

A total of 1467 mutants of the biocontrol bacterium Bacillus pumilus DX01 were obtained by Tn5 insertional mutagenesis and subjected to the determination of antagonistic capabilities. Compared with the wild-type strain DX01, the mutant M25 was identified to have the most significant reduction in antagonistic capability against the phytopathogen Bipolaris maydis and extracellular proteinase activity. The integration site of the exogenous T-DNA in the genome of mutant M25 was revealed in the coding region of malony CoA-ACP transacylase (MCAT) gene (mcat), which belongs to a polyketide synthase (PKS) gene cluster, DX01pks of B. pumilus DX01. Furthermore, the whole DX01pks gene cluster was cloned using Illumina Solexa sequencing technology, and it has a modular framework different from the other two gene clusters involved in polyketide synthesis in B. amyloliquefaciens FZB42 (pks1) and B. subtilis 168 (pksX). Finally, in order to gain more insights into the molecular mechanisms of biocontrol of B. pumilus DX01, the changes in the relative level of expression of total proteins between the original strain DX01 and the mutant M25 were detected by 2-DE-based proteomic analysis. A total of twenty differentially expressed proteins were identified upon the mcat gene transposition mutagenesis. Of these proteins, seven proteins were up-regulated in expression level and the other proteins were down-regulated. Taken together, the results in this study showed that Tn5 transposon mutagenesis of B. pumilus DX01 can lead to a significant change of antiphytopathogen ability, and the DX01pks gene cluster possibly play a potential role in the biocontrol processes of this bacterium.

Comprehensive transcriptomic and proteomic analyses of antroquinonol biosynthetic genes and enzymes in Antrodia camphorata

Antroquinonol (AQ) has several remarkable bioactivities in acute myeloid leukaemia and pancreatic cancer, but difficulties in the mass production of AQ hamper its applications. Currently, molecular biotechnology methods, such as gene overexpression, have been widely used to increase the production of metabolites. However, AQ biosynthetic genes and enzymes are poorly understood. In this study, an integrated study coupling RNA-Seq and isobaric tags for relative and absolute quantitation (iTRAQ) were used to identify AQ synthesis-related genes and enzymes in Antrodia camphorata during coenzyme Q0-induced fermentation (FM). The upregulated genes related to acetyl-CoA synthesis indicated that acetyl-CoA enters the mevalonate pathway to form the farnesyl tail precursor of AQ. The metE gene for an enzyme with methyl transfer activity provided sufficient methyl groups for AQ structure formation. The CoQ2 and ubiA genes encode p-hydroxybenzoate polyprenyl transferase, linking coenzyme Q0 and the polyisoprene side chain to form coenzyme Q3. NADH is transformed into NAD+ and releases two electrons, which may be beneficial for the conversion of coenzyme Q3 to AQ. Understanding the biosynthetic genes and enzymes of AQ is important for improving its production by genetic means in the future.

Dysbiosis is not present in horses with fecal water syndrome when compared to controls in spring and autumn

BACKGROUND

Fecal water syndrome (FWS) is long-standing and common in horses, particularly in central Europe. No large epidemiological data sets exist, and the cause remains elusive. Dysbiosis could play a role in pathogenesis.

OBJECTIVES

To evaluate whether dysbiosis is present in horses with FWS when compared to stable-matched control horses in spring and autumn.

ANIMALS

Fecal samples were collected from horses with FWS (n = 16; 9 mares, 7 geldings) and controls (n = 15; 8 mares, 7 geldings).

METHODS

The bacterial microbiome of samples collected in spring and autumn of 2016 was analyzed using high-throughput sequencing. Differences in relative abundance of bacterial taxa, alpha diversity, and beta diversity indices were assessed between horses with FWS and controls based on season.

RESULTS

Differences in microbial community composition based on time point and health status were not observed on any taxonomic level. Limited differences were seen on linear discriminant analysis effect size analysis. No difference in alpha diversity indices was observed including richness, diversity based on health status, or time point. No effect of health status on microbial community membership structure was observed.

CONCLUSIONS AND CLINICAL IMPORTANCE

Limited differences were found in the bacterial microbiota of horses with and without FWS, regardless of season. Further research is needed to elucidate the role of microbiota in the development of FWS.

dbAMP: an integrated resource for exploring antimicrobial peptides with functional activities and physicochemical properties on transcriptome and proteome data

Antimicrobial peptides (AMPs), naturally encoded from genes and generally contained 10–100 amino acids, are crucial components of the innate immune system and can protect the host from various pathogenic bacteria, as well as viruses. In recent years, the widespread use of antibiotics has inspired the rapid growth of antibiotic-resistant microorganisms that usually induce critical infection and pathogenesis. An increasing interest therefore was motivated to explore natural AMPs that enable the development of new antibiotics. With the potential of AMPs being as new drugs for multidrug-resistant pathogens, we were thus motivated to develop a database (dbAMP, http://csb.cse.yzu.edu.tw/dbAMP/) by accumulating comprehensive AMPs from public domain and manually curating literature. Currently in dbAMP there are 12 389 unique entries, including 4271 experimentally verified AMPs and 8118 putative AMPs along with their functional activities, supported by 1924 research articles. The advent of high-throughput biotechnologies, such as mass spectrometry and next-generation sequencing, has led us to further expand dbAMP as a database-assisted platform for providing comprehensively functional and physicochemical analyses for AMPs based on the large-scale transcriptome and proteome data. Significant improvements available in dbAMP include the information of AMP–protein interactions, antimicrobial potency analysis for ‘cryptic’ region detection, annotations of AMP target species, as well as AMP detection on transcriptome and proteome datasets. Additionally, a Docker container has been developed as a downloadable package for discovering known and novel AMPs on high-throughput omics data. The user-friendly visualization interfaces have been created to facilitate peptide searching, browsing, and sequence alignment against dbAMP entries. All the facilities integrated into dbAMP can promote the functional analyses of AMPs and the discovery of new antimicrobial drugs.

Genome-Wide Profiling of Enterotoxigenic Staphylococcus aureus Strains Used for the Production of Naturally Contaminated Cheeses

Staphylococcus aureus is a major human pathogen and an important cause of livestock infections. More than 20 staphylococcal enterotoxins with emetic activity can be produced by specific strains responsible for staphylococcal food poisoning, one of the most common food-borne diseases. Whole genome sequencing provides a comprehensive view of the genome structure and gene content that have largely been applied in outbreak investigations and genomic comparisons. In this study, six enterotoxigenic S. aureus strains were characterised using a combination of molecular, phenotypical and computational methods. The genomes were analysed for the presence of virulence factors (VFs), where we identified 110 genes and classified them into five categories: adherence (n = 31), exoenzymes (n = 28), genes involved in host immune system evasion (n = 7); iron uptake regulatory system (n = 8); secretion machinery factors and toxins' genes (n = 36), and 39 genes coding for transcriptional regulators related to staphylococcal VFs. Each group of VFs revealed correlations among the six enterotoxigenic strains, and further analysis revealed their accessory genomic content, including mobile genetic elements. The plasmids pLUH02 and pSK67 were detected in the strain ProNaCC1 and ProNaCC7, respectively, carrying out the genes sed, ser, and selj. The genes carried out by prophages were detected in the strain ProNaCC2 (see), ProNaCC4, and ProNaCC7 (both positive for sea). The strain ProNaCC5 resulted positive for the genes seg, sei, sem, sen, seo grouped in an exotoxin gene cluster, and the strain ProNaCC6 resulted positive for seh, a transposon-associated gene. The six strains were used for the production of naturally contaminated cheeses which were tested with the European Screening Method for staphylococcal enterotoxins. The results obtained from the analysis of toxins produced in cheese, combined with the genomic features represent a portrait of the strains that can be used for the production of staphylococcal enterotoxin-positive cheese as reference material.

Comprehensive genomic analysis of an indigenous Pseudomonas pseudoalcaligenes degrading phenolic compounds

Environmental contamination with aromatic compounds is a universal challenge. Aromatic-degrading microorganisms isolated from the same or similar polluted environments seem to be more suitable for bioremediation. Moreover, microorganisms adapted to contaminated environments are able to use toxic compounds as the sole sources of carbon and energy. An indigenous strain of Pseudomonas, isolated from the Mahshahr Petrochemical plant in the Khuzestan province, southwest of Iran, was studied genetically. It was characterized as a novel Gram-negative, aerobic, halotolerant, rod-shaped bacterium designated Pseudomonas YKJ, which was resistant to chloramphenicol and ampicillin. Genome of the strain was completely sequenced using Illumina technology to identify its genetic characteristics. MLST analysis revealed that the YKJ strain belongs to the genus Pseudomonas indicating the highest sequence similarity with Pseudomonas pseudoalcaligenes strain CECT 5344 (99% identity). Core- and pan-genome analysis indicated that P. pseudoalcaligenes contains 1,671 core and 3,935 unique genes for coding DNA sequences. The metabolic and degradation pathways for aromatic pollutants were investigated using the NCBI and KEGG databases. Genomic and experimental analyses showed that the YKJ strain is able to degrade certain aromatic compounds including bisphenol A, phenol, benzoate, styrene, xylene, benzene and chlorobenzene. Moreover, antibiotic resistance and chemotaxis properties of the YKJ strain were found to be controlled by two-component regulatory systems.

Exploring the Hospital Microbiome by High-Resolution 16S rRNA Profiling

The aim of this work was to analyze and compare the bacterial communities of 663 samples from a Brazilian hospital by using high-throughput sequencing of the 16S rRNA gene. To increase taxonomic profiling and specificity of 16S-based identification, a strict sequence quality filtering process was applied for the accurate identification of clinically relevant bacterial taxa. Our results indicate that the hospital environment is predominantly inhabited by closely related species. A massive dominance of a few taxa in all taxonomic levels down to the genera was observed, where the ten most abundant genera in each facility represented 64.4% of all observed taxa, with a major predominance of Acinetobacter and Pseudomonas. The presence of several nosocomial pathogens was revealed. Co-occurrence analysis indicated that the present hospital microbial network had low connectedness, forming a clustered topology, but not structured among groups of nodes (i.e., modules). Furthermore, we were able to detect ecologically relevant relationships between specific microbial taxa, in particular, potential competition between pathogens and non-pathogens. Overall, these results provide new insight into different aspects of a hospital microbiome and indicate that 16S rRNA sequencing may serve as a robust one-step tool for microbiological identification and characterization of a wide range of clinically relevant bacterial taxa in hospital settings with a high resolution.

Comparative transcriptome analysis reveals relationship of three major domesticated varieties of Auricularia auricula-judae

Auricularia auricula-judae is an edible mushroom and a traditional medicine in China as well as the fourth largest cultivated mushroom species in the world. Here for the first time, we present comparative transcriptome analyses of the fruiting bodies of three morphologically distinguishable A. auricula-judae cultivated varieties (Wujin, smooth; Banjin, partially wrinkled; and Quanjin, fully wrinkled) collected from Jilin Province, China. Biological triplicates were performed to determine the expression levels of 13,937 unigenes. Among them, only 13 unigenes were annotated to A. auricula-judae, highlighting the lack of publicly available reference sequences for this economically important species. Principal component analysis (PCA) determined that the gene expression profile of Quanjin was unique when compared to those of Banjin and Wujin. Such relationships were further supported by analyses of annotated and unannotated unigenes, differentially expressed unigenes, gene ontology functions, and the family of peroxidase genes. Using the KEGG database, significant alternations in biological pathways were detected among the three cultivars. This work contributes a large set of A. auricula-judae sequences to public database, establishes the relationships among major cultivars, and provides molecular guidance for breeding and cultivation.

Exploring the Human Microbiome: The Potential Future Role of Next-Generation Sequencing in Disease Diagnosis and Treatment

The interaction between the human microbiome and immune system has an effect on several human metabolic functions and impacts our well-being. Additionally, the interaction between humans and microbes can also play a key role in determining the wellness or disease status of the human body. Dysbiosis is related to a plethora of diseases, including skin, inflammatory, metabolic, and neurological disorders. A better understanding of the host-microbe interaction is essential for determining the diagnosis and appropriate treatment of these ailments. The significance of the microbiome on host health has led to the emergence of new therapeutic approaches focused on the prescribed manipulation of the host microbiome, either by removing harmful taxa or reinstating missing beneficial taxa and the functional roles they perform. Culturing large numbers of microbial taxa in the laboratory is problematic at best, if not impossible. Consequently, this makes it very difficult to comprehensively catalog the individual members comprising a specific microbiome, as well as understanding how microbial communities function and influence host-pathogen interactions. Recent advances in sequencing technologies and computational tools have allowed an increasing number of metagenomic studies to be performed. These studies have provided key insights into the human microbiome and a host of other microbial communities in other environments. In the present review, the role of the microbiome as a therapeutic agent and its significance in human health and disease is discussed. Advances in high-throughput sequencing technologies for surveying host-microbe interactions are also discussed. Additionally, the correlation between the composition of the microbiome and infectious diseases as described in previously reported studies is covered as well. Lastly, recent advances in state-of-the-art bioinformatics software, workflows, and applications for analysing metagenomic data are summarized.

Genome sequencing and assessment of plant growth-promoting properties of a Serratia marcescens strain isolated from vermicompost

Background

Plant-bacteria associations have been extensively studied for their potential in increasing crop productivity in a sustainable manner. Serratia marcescens is a species of Enterobacteriaceae found in a wide range of environments, including soil.

Results

Here we describe the genome sequencing and assessment of plant growth-promoting abilities of S. marcescens UENF-22GI, a strain isolated from mature cattle manure vermicompost. In vitro, S. marcescens UENF-22GI is able to solubilize P and Zn, to produce indole compounds (likely IAA), to colonize hyphae and counter the growth of two phytopathogenic fungi. Inoculation of maize with this strain remarkably increased seedling growth and biomass under greenhouse conditions. The S. marcescens UENF-22GI genome has 5 Mb, assembled in 17 scaffolds comprising 4662 genes (4528 are protein-coding). No plasmids were identified. S. marcescens UENF-22GI is phylogenetically placed within a clade comprised almost exclusively of non-clinical strains. We identified genes and operons that are likely responsible for the interesting plant-growth promoting features that were experimentally described. The S. marcescens UENF-22GI genome harbors a horizontally-transferred genomic island involved in antibiotic production, antibiotic resistance, and anti-phage defense via a novel ADP-ribosyltransferase-like protein and possible modification of DNA by a deazapurine base, which likely contributes to its competitiveness against other bacteria.

Conclusions

Collectively, our results suggest that S. marcescens UENF-22GI is a strong candidate to be used in the enrichment of substrates for plant growth promotion or as part of bioinoculants for agriculture.

Electronic supplementary material

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

The Planktonic Core Microbiome and Core Functions in the Cattle Rumen by Next Generation Sequencing

The cow rumen harbors a great variety of diverse microbes, which form a complex, organized community. Understanding the behavior of this multifarious network is crucial in improving ruminant nutrient use efficiency. The aim of this study was to expand our knowledge by examining 10 Holstein dairy cow rumen fluid fraction whole metagenome and transcriptome datasets. DNA and mRNA sequence data, generated by Ion Torrent, was subjected to quality control and filtering before analysis for core elements. The taxonomic core microbiome consisted of 48 genera belonging to Bacteria (47) and Archaea (1). The genus Prevotella predominated the planktonic core community. Core functional groups were identified using co-occurrence analysis and resulted in 587 genes, from which 62 could be assigned to metabolic functions. Although this was a minimal functional core, it revealed key enzymes participating in various metabolic processes. A diverse and rich collection of enzymes involved in carbohydrate metabolism and other functions were identified. Transcripts coding for enzymes active in methanogenesis made up 1% of the core functions. The genera associated with the core enzyme functions were also identified. Linking genera to functions showed that the main metabolic pathways are primarily provided by Bacteria and several genera may serve as a “back-up” team for the central functions. The key actors in most essential metabolic routes belong to the genus Prevotella. Confirming earlier studies, the genus Methanobrevibacter carries out the overwhelming majority of rumen methanogenesis and therefore methane emission mitigation seems conceivable via targeting the hydrogenotrophic methanogenesis.

DeepSimulator: a deep simulator for Nanopore sequencing

Motivation

Oxford Nanopore sequencing is a rapidly developed sequencing technology in recent years. To keep pace with the explosion of the downstream data analytical tools, a versatile Nanopore sequencing simulator is needed to complement the experimental data as well as to benchmark those newly developed tools. However, all the currently available simulators are based on simple statistics of the produced reads, which have difficulty in capturing the complex nature of the Nanopore sequencing procedure, the main task of which is the generation of raw electrical current signals.

Results

Here we propose a deep learning based simulator, DeepSimulator, to mimic the entire pipeline of Nanopore sequencing. Starting from a given reference genome or assembled contigs, we simulate the electrical current signals by a context-dependent deep learning model, followed by a base-calling procedure to yield simulated reads. This workflow mimics the sequencing procedure more naturally. The thorough experiments performed across four species show that the signals generated by our context-dependent model are more similar to the experimentally obtained signals than the ones generated by the official context-independent pore model. In terms of the simulated reads, we provide a parameter interface to users so that they can obtain the reads with different accuracies ranging from 83 to 97%. The reads generated by the default parameter have almost the same properties as the real data. Two case studies demonstrate the application of DeepSimulator to benefit the development of tools in de novo assembly and in low coverage SNP detection.

Availability and implementation

The software can be accessed freely at: https://github.com/lykaust15/DeepSimulator.

Supplementary information

Supplementary data are available at Bioinformatics online.

Nanopore sequencing and full genome de novo assembly of human cytomegalovirus TB40/E reveals clonal diversity and structural variations

BACKGROUND

Human cytomegalovirus (HCMV) has a double-stranded DNA genome of approximately 235 Kbp that is structurally complex including extended GC-rich repeated regions. Genomic recombination events are frequent in HCMV cultures but have also been observed in vivo. Thus, the assembly of HCMV whole genomes from technologies producing shorter than 500 bp sequences is technically challenging. Here we improved the reconstruction of HCMV full genomes by means of a hybrid, de novo genome-assembly bioinformatics pipeline upon data generated from the recently released MinION MkI B sequencer from Oxford Nanopore Technologies.

RESULTS

The MinION run of the HCMV (strain TB40/E) library resulted in ~ 47,000 reads from a single R9 flowcell and in ~ 100× average read depth across the virus genome. We developed a novel, self-correcting bioinformatics algorithm to assemble the pooled HCMV genomes in three stages. In the first stage of the bioinformatics algorithm, long contigs (N50 = 21,892) of lower accuracy were reconstructed. In the second stage, short contigs (N50 = 5686) of higher accuracy were assembled, while in the final stage the high quality contigs served as template for the correction of the longer contigs resulting in a high-accuracy, full genome assembly (N50 = 41,056). We were able to reconstruct a single representative haplotype without employing any scaffolding steps. The majority (98.8%) of the genomic features from the reference strain were accurately annotated on this full genome construct. Our method also allowed the detection of multiple alternative sub-genomic fragments and non-canonical structures suggesting rearrangement events between the unique (UL /US) and the repeated (T/IRL/S) genomic regions.

CONCLUSIONS

Third generation high-throughput sequencing technologies can accurately reconstruct full-length HCMV genomes including their low-complexity and highly repetitive regions. Full-length HCMV genomes could prove crucial in understanding the genetic determinants and viral evolution underpinning drug resistance, virulence and pathogenesis.

Evaluation of 24-locus MIRU-VNTR genotyping in Mycobacterium tuberculosis cluster investigations in four jurisdictions in the United States, 2006-2010

The U.S. Centers for Disease Control and Prevention (CDC) uses a combination of spacer oligonucleotide typing (spoligotyping) and mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU-VNTR) analyses as part of the National TB Genotyping Service (NTGS). The NTGS expansion from 12-locus MIRU-VNTR (MIRU12) to 24-locus MIRU-VNTR (MIRU24) in 2009 enhanced the ability to discriminate Mycobacterium tuberculosis strains. In the current study, we investigated the MIRU24 concordance among epidemiologic-linked tuberculosis (TB) patients in four U.S. health jurisdictions. We also evaluated the programmatic benefits of combining MIRU24 and spoligotyping with epidemiologic evidence in identifying potential recent TB transmission. We examined 342 TB patients in 42 spoligotype/MIRU12 (PCRType) clusters (equivalent to 46 spoligotype/MIRU24 [GENType] clusters) to identify epidemiologic links among cases. GENType clusters, when compared to PCRType clusters, had 12 times higher odds of epidemiologic links being identified if patients were younger than 25 years and 3 times higher odds if patients resided in the same zip code, or had HIV infection. Sixty (18%) fewer PCRType-clustered patients would need investigations if clusters are defined using GENType instead of PCRType. An important advantage of defining clusters by MIRU24 is resource savings related to the reduced number of clustered cases needing investigation.

Elucidation of quantitative structural diversity of remarkable rearrangement regions, shufflons, in IncI2 plasmids

A multiple DNA inversion system, the shufflon, exists in incompatibility (Inc) I1 and I2 plasmids. The shufflon generates variants of the PilV protein, a minor component of the thin pilus. The shufflon is one of the most difficult regions for de novo genome assembly because of its structural diversity even in an isolated bacterial clone. We determined complete genome sequences, including those of IncI2 plasmids carrying mcr-1, of three Escherichia coli strains using single-molecule, real-time (SMRT) sequencing and Illumina sequencing. The sequences assembled using only SMRT sequencing contained misassembled regions in the shufflon. A hybrid analysis using SMRT and Illumina sequencing resolved the misassembled region and revealed that the three IncI2 plasmids, excluding the shufflon region, were highly conserved. Moreover, the abundance ratio of whole-shufflon structures could be determined by quantitative structural variation analysis of the SMRT data, suggesting that a remarkable heterogeneity of whole-shufflon structural variations exists in IncI2 plasmids. These findings indicate that remarkable rearrangement regions should be validated using both long-read and short-read sequencing data and that the structural variation of PilV in the shufflon might be closely related to phenotypic heterogeneity of plasmid-mediated transconjugation involved in horizontal gene transfer even in bacterial clonal populations.

Investigating right ovary degeneration in chick embryos by transcriptome sequencing

In asymmetric chick gonads, the left and right female gonads undergo distinct programs during development, generating a functional ovary on the left side only. Despite some progress being made in recent years, the mechanisms of molecular regulation remain incompletely understood, and little genomic information is available regarding the degeneration of the right ovary in the chick embryo testis. In this study, we performed transcriptome sequencing to investigate differentially expressed genes in the left and right ovaries and gene functions at two critical time points; embryonic days 6 (E6) and 10 (E10). Using high-throughput RNA-sequencing technologies, 539 and 1046 genes were identified as being significantly differentially expressed between 6R-VS-6L and 10R-VS-10L. Gene ontology analysis of the differentially expressed genes revealed enrichment in functional pathways. Among these, candidate genes associated with degeneration of the right ovary in the chick embryo were identified. Identification of a pathway involved in ovarian degeneration provides an important resource for the further study of its molecular mechanisms and functions.

Bioinformatic Analyses of Unique (Orphan) Core Genes of the Genus Acidithiobacillus: Functional Inferences and Use As Molecular Probes for Genomic and Metagenomic/Transcriptomic Interrogation

Using phylogenomic and gene compositional analyses, five highly conserved gene families have been detected in the core genome of the phylogenetically coherent genus Acidithiobacillus of the class Acidithiobacillia. These core gene families are absent in the closest extant genus Thermithiobacillus tepidarius that subtends the Acidithiobacillus genus and roots the deepest in this class. The predicted proteins encoded by these core gene families are not detected by a BLAST search in the NCBI non-redundant database of more than 90 million proteins using a relaxed cut-off of 1.0e⁻⁵. None of the five families has a clear functional prediction. However, bioinformatic scrutiny, using pI prediction, motif/domain searches, cellular location predictions, genomic context analyses, and chromosome topology studies together with previously published transcriptomic and proteomic data, suggests that some may have functions associated with membrane remodeling during cell division perhaps in response to pH stress. Despite the high level of amino acid sequence conservation within each family, there is sufficient nucleotide variation of the respective genes to permit the use of the DNA sequences to distinguish different species of Acidithiobacillus, making them useful additions to the armamentarium of tools for phylogenetic analysis. Since the protein families are unique to the Acidithiobacillus genus, they can also be leveraged as probes to detect the genus in environmental metagenomes and metatranscriptomes, including industrial biomining operations, and acid mine drainage (AMD).

Bg10: A Novel Metagenomics Alcohol-Tolerant and Glucose-Stimulated GH1 ß-Glucosidase Suitable for Lactose-Free Milk Preparation

New ß-glucosidases with product (glucose) or ethanol tolerances are greatly desired to make industrial processes more marketable and efficient. Therefore, this report describes the in silico/vitro characterization of Bg10, a metagenomically derived homodimeric ß-glucosidase that exhibited a Vmax of 10.81 ± 0.43 μM min-1, Kcat of 175.1± 6.91 min-1, and Km of 0.49 ± 0.12 mM at a neutral pH and 37°C when pNP-ß-D-glucopyranoside was used as the substrate, and the enzyme retained greater than 80% activity within the respective pH and temperature ranges of 6.5 to 8.0 and 35 to 40°C. The enzyme was stimulated by its product, glucose; consequently, the Bg10 activity against 50 and 100 mM of glucose were increased by 36.8% and 22%, respectively, while half of the activity was retained at 350 mM. Moreover, the Bg10 was able to hydrolyse 55% (milk sample) and 100% (purified sugar) of the lactose at low (6°C) and optimum (37°C) temperatures, respectively, suggesting the possibility of further optimization of the reaction for lactose-free dairy production. In addition, the enzyme was able to fully hydrolyse 40 mM of cellobiose at one hour and was tolerant to ethanol up to concentrations of 500 mM (86% of activity), while a 1 M concentration still resulted in 41% residual activity, which could be interesting for biofuel production.

Gene Turnover Contributes to the Evolutionary Adaptation of Acidithiobacillus caldus: Insights from Comparative Genomics

Acidithiobacillus caldus is an extremely acidophilic sulfur-oxidizer with specialized characteristics, such as tolerance to low pH and heavy metal resistance. To gain novel insights into its genetic complexity, we chosen six A. caldus strains for comparative survey. All strains analyzed in this study differ in geographic origins as well as in ecological preferences. Based on phylogenomic analysis, we clustered the six A. caldus strains isolated from various ecological niches into two groups: group 1 strains with smaller genomes and group 2 strains with larger genomes. We found no obvious intraspecific divergence with respect to predicted genes that are related to central metabolism and stress management strategies between these two groups. Although numerous highly homogeneous genes were observed, high genetic diversity was also detected. Preliminary inspection provided a first glimpse of the potential correlation between intraspecific diversity at the genome level and environmental variation, especially geochemical conditions. Evolutionary genetic analyses further showed evidence that the difference in environmental conditions might be a crucial factor to drive the divergent evolution of A. caldus species. We identified a diverse pool of mobile genetic elements including insertion sequences and genomic islands, which suggests a high frequency of genetic exchange in these harsh habitats. Comprehensive analysis revealed that gene gains and losses were both dominant evolutionary forces that directed the genomic diversification of A. caldus species. For instance, horizontal gene transfer and gene duplication events in group 2 strains might contribute to an increase in microbial DNA content and novel functions. Moreover, genomes undergo extensive changes in group 1 strains such as removal of potential non-functional DNA, which results in the formation of compact and streamlined genomes. Taken together, the findings presented herein show highly frequent gene turnover of A. caldus species that inhabit extremely acidic environments, and shed new light on the contribution of gene turnover to the evolutionary adaptation of acidophiles.

Amino Acids Hydrolyzed from Animal Carcasses Are a Good Additive for the Production of Bio-organic Fertilizer

High-quality bio-organic fertilizers (BIOs) cannot be produced without the addition of some proteins. In this study, compound liquid amino acids (CLAA) from animal carcasses were utilized as additives into matured composts to create novel BIOs containing plant growth-promoting rhizobacteria (PGPR). The results showed that adding CLAA and inoculating bacteria meanwhile resulted in failed solid-state fermentation (SSF) due to the higher H(+) contents. While after pre-compost for 4 days before PGPR inoculation, treatments of matured chicken or pig manure added with 0.2 ml g(-1) of CLAA resulted in a maximum biomass of functional strains. Illumine-MiSeq sequencing and Real-Time PCR results showed that the CLAA addition decreased the bacterial abundance and richness, altered the bacterial community structure and changed the relative abundance of some microbial groups. This study offers a high value-added utilization of waste protein resources for producing economical, high-quality BIO.

Role of Escherichia coli in Biofuel Production

Increased energy consumption coupled with depleting petroleum reserves and increased greenhouse gas emissions have renewed our interest in generating fuels from renewable energy sources via microbial fermentation. Central to this problem is the choice of microorganism that catalyzes the production of fuels at high volumetric productivity and yield from cheap and abundantly available renewable energy sources. Microorganisms that are metabolically engineered to redirect renewable carbon sources into desired fuel products are contemplated as best choices to obtain high volumetric productivity and yield. Considering the availability of vast knowledge in genomic and metabolic fronts, Escherichia coli is regarded as a primary choice for the production of biofuels. Here, we reviewed the microbial production of liquid biofuels that have the potential to be used either alone or in combination with the present-day fuels. We specifically highlighted the metabolic engineering and synthetic biology approaches used to improve the production of biofuels from E. coli over the past few years. We also discussed the challenges that still exist for the biofuel production from E. coli and their possible solutions.

De Novo Assembly of Human Herpes Virus Type 1 (HHV-1) Genome, Mining of Non-Canonical Structures and Detection of Novel Drug-Resistance Mutations Using Short- and Long-Read Next Generation Sequencing Technologies

Human herpesvirus type 1 (HHV-1) has a large double-stranded DNA genome of approximately 152 kbp that is structurally complex and GC-rich. This makes the assembly of HHV-1 whole genomes from short-read sequencing data technically challenging. To improve the assembly of HHV-1 genomes we have employed a hybrid genome assembly protocol using data from two sequencing technologies: the short-read Roche 454 and the long-read Oxford Nanopore MinION sequencers. We sequenced 18 HHV-1 cell culture-isolated clinical specimens collected from immunocompromised patients undergoing antiviral therapy. The susceptibility of the samples to several antivirals was determined by plaque reduction assay. Hybrid genome assembly resulted in a decrease in the number of contigs in 6 out of 7 samples and an increase in N(G)50 and N(G)75 of all 7 samples sequenced by both technologies. The approach also enhanced the detection of non-canonical contigs including a rearrangement between the unique (UL) and repeat (T/IRL) sequence regions of one sample that was not detectable by assembly of 454 reads alone. We detected several known and novel resistance-associated mutations in UL23 and UL30 genes. Genome-wide genetic variability ranged from <1% to 53% of amino acids in each gene exhibiting at least one substitution within the pool of samples. The UL23 gene had one of the highest genetic variabilities at 35.2% in keeping with its role in development of drug resistance. The assembly of accurate, full-length HHV-1 genomes will be useful in determining genetic determinants of drug resistance, virulence, pathogenesis and viral evolution. The numerous, complex repeat regions of the HHV-1 genome currently remain a barrier towards this goal.

Detection and characterization of two co-infection variant strains of avian orthoreovirus (ARV) in young layer chickens using next-generation sequencing (NGS)

Using next-generation sequencing (NGS) for full genomic characterization studies of the newly emerging avian orthoreovirus (ARV) field strains isolated in Pennsylvania poultry, we identified two co-infection ARV variant strains from one ARV isolate obtained from ARV-affected young layer chickens. The de novo assembly of the ARV reads generated 19 contigs of two different ARV variant strains according to 10 genome segments of each ARV strain. The two variants had the same M2 segment. The complete genomes of each of the two variant strains were 23,493 bp in length, and 10 dsRNA segments ranged from 1192 bp (S4) to 3958 bp (L1), encoding 12 viral proteins. Sequence comparison of nucleotide (nt) and amino acid (aa) sequences of all 10 genome segments revealed 58.1–100% and 51.4–100% aa identity between the two variant strains, and 54.3–89.4% and 49.5–98.1% aa identity between the two variants and classic vaccine strains. Phylogenetic analysis revealed a moderate to significant nt sequence divergence between the two variant and ARV reference strains. These findings have demonstrated the first naturally occurring co-infection of two ARV variants in commercial young layer chickens, providing scientific evidence that multiple ARV strains can be simultaneously present in one host species of chickens.

Recent Developments in Using Advanced Sequencing Technologies for the Genomic Studies of Lignin and Cellulose Degrading Microorganisms

Lignin is a complex polyphenyl aromatic compound which exists in tight associations with cellulose and hemicellulose to form plant primary and secondary cell wall. Lignocellulose is an abundant renewable biomaterial present on the earth. It has gained much attention in the scientific community in recent years because of its potential applications in bio-based industries. Microbial degradation of lignocellulose polymers was well studied in wood decaying fungi. Based on the plant materials they degrade these fungi were classified as white rot, brown rot and soft rot. However, some groups of bacteria belonging to the actinomycetes, α-proteobacteria and β-proteobacteria were also found to be efficient in degrading lignocellulosic biomass but not well understood unlike the fungi. In this review we focus on recent advancements deployed for finding and understanding the lignocellulose degradation by microorganisms. Conventional molecular methods like sequencing 16s rRNA and Inter Transcribed Spacer (ITS) regions were used for identification and classification of microbes. Recent progression in genomics mainly next generation sequencing technologies made the whole genome sequencing of microbes possible in a great ease. The whole genome sequence studies reveals high quality information about genes and canonical pathways involved in the lignin and other cell wall components degradation.

Opening the black box of Anaplasma phagocytophilum diversity: current situation and future perspectives

Anaplasma phagocytophilum is a zoonotic obligate intracellular bacterium known to be transmitted by ticks belonging to the Ixodes persulcatus complex. This bacterium can infect several mammalian species, and is known to cause diseases with variable symptoms in many domestic animals. Specifically, it is the causative agent of tick-borne fever (TBF), a disease of important economic impact in European domestic ruminants, and human granulocytic anaplasmosis (HGA), an emerging zoonotic disease in Asia, USA and Europe. A. phagocytophilum epidemiological cycles are complex and involve different ecotypes, vectors, and mammalian host species. Moreover, the epidemiology of A. phagocytophilum infection differs greatly between Europe and the USA. These different epidemiological contexts are associated with considerable variations in bacterial strains. Until recently, few A. phagocytophilum molecular typing tools were available, generating difficulties in completely elucidating the epidemiological cycles of this bacterium. Over the last few years, many A. phagocytophilum typing techniques have been developed, permitting in-depth epidemiological exploration. Here, we review the current knowledge and future perspectives regarding A. phagocytophilum epidemiology and phylogeny, and then focus on the molecular typing tools available for studying A. phagocytophilum genetic diversity.

Genetic markers, genotyping methods & next generation sequencing in Mycobacterium tuberculosis

Molecular epidemiology (ME) is one of the main areas in tuberculosis research which is widely used to study the transmission epidemics and outbreaks of tubercle bacilli. It exploits the presence of various polymorphisms in the genome of the bacteria that can be widely used as genetic markers. Many DNA typing methods apply these genetic markers to differentiate various strains and to study the evolutionary relationships between them. The three widely used genotyping tools to differentiate Mycobacterium tuberculosis strains are IS6110 restriction fragment length polymorphism (RFLP), spacer oligotyping (Spoligotyping), and mycobacterial interspersed repeat units - variable number of tandem repeats (MIRU-VNTR). A new prospect towards ME was introduced with the development of whole genome sequencing (WGS) and the next generation sequencing (NGS) methods, where the entire genome is sequenced that not only helps in pointing out minute differences between the various sequences but also saves time and the cost. NGS is also found to be useful in identifying single nucleotide polymorphisms (SNPs), comparative genomics and also various aspects about transmission dynamics. These techniques enable the identification of mycobacterial strains and also facilitate the study of their phylogenetic and evolutionary traits.

VirAmp: a galaxy-based viral genome assembly pipeline

Background

Advances in next generation sequencing make it possible to obtain high-coverage sequence data for large numbers of viral strains in a short time. However, since most bioinformatics tools are developed for command line use, the selection and accessibility of computational tools for genome assembly and variation analysis limits the ability of individual labs to perform further bioinformatics analysis.

Findings

We have developed a multi-step viral genome assembly pipeline named VirAmp, which combines existing tools and techniques and presents them to end users via a web-enabled Galaxy interface. Our pipeline allows users to assemble, analyze, and interpret high coverage viral sequencing data with an ease and efficiency that was not possible previously. Our software makes a large number of genome assembly and related tools available to life scientists and automates the currently recommended best practices into a single, easy to use interface. We tested our pipeline with three different datasets from human herpes simplex virus (HSV).

Conclusions

VirAmp provides a user-friendly interface and a complete pipeline for viral genome analysis. We make our software available via an Amazon Elastic Cloud disk image that can be easily launched by anyone with an Amazon web service account. A fully functional demonstration instance of our system can be found at http://viramp.com/. We also maintain detailed documentation on each tool and methodology at http://docs.viramp.com.

Electronic supplementary material

The online version of this article (doi:10.1186/s13742-015-0060-y) contains supplementary material, which is available to authorized users.

MicroRNA-like small RNAs prediction in the development of Antrodia cinnamomea

Antrodia cinnamomea, a precious, host-specific brown-rot fungus that has been used as a folk medicine in Taiwan for centuries is known to have diverse bioactive compounds with potent pharmaceutical activity. In this study, different fermentation states of A. cinnamomea (wild-type fruiting bodies and liquid cultured mycelium) were sequenced using the next-generation sequencing (NGS) technique. A 45.58 Mb genome encoding 6,522 predicted genes was obtained. High quality reads were assembled into a total of 13,109 unigenes. Using a previously constructed pipeline to search for microRNAs (miRNAs), we then identified 4 predicted conserved miRNA and 63 novel predicted miRNA-like small RNA (milRNA) candidates. Target prediction revealed several interesting proteins involved in tri-terpenoid synthesis, mating type recognition, chemical or physical sensory protein and transporters predicted to be regulated by the miRNAs and milRNAs.

Comparative whole-genome analysis of clinical isolates reveals characteristic architecture of Mycobacterium tuberculosis pangenome

The tubercle complex consists of closely related mycobacterium species which appear to be variants of a single species. Comparative genome analysis of different strains could provide useful clues and insights into the genetic diversity of the species. We integrated genome assemblies of 96 strains from Mycobacterium tuberculosis complex (MTBC), which included 8 Indian clinical isolates sequenced and assembled in this study, to understand its pangenome architecture. We predicted genes for all the 96 strains and clustered their respective CDSs into homologous gene clusters (HGCs) to reveal a hard-core, soft-core and accessory genome component of MTBC. The hard-core (HGCs shared amongst 100% of the strains) was comprised of 2,066 gene clusters whereas the soft-core (HGCs shared amongst at least 95% of the strains) comprised of 3,374 gene clusters. The change in the core and accessory genome components when observed as a function of their size revealed that MTBC has an open pangenome. We identified 74 HGCs that were absent from reference strains H37Rv and H37Ra but were present in most of clinical isolates. We report PCR validation on 9 candidate genes depicting 7 genes completely absent from H37Rv and H37Ra whereas 2 genes shared partial homology with them accounting to probable insertion and deletion events. The pangenome approach is a promising tool for studying strain specific genetic differences occurring within species. We also suggest that since selecting appropriate target genes for typing purposes requires the expected target gene be present in all isolates being typed, therefore estimating the core-component of the species becomes a subject of prime importance.

Large-scale contamination of microbial isolate genomes by Illumina PhiX control

With the rapid growth and development of sequencing technologies, genomes have become the new go-to for exploring solutions to some of the world’s biggest challenges such as searching for alternative energy sources and exploration of genomic dark matter. However, progress in sequencing has been accompanied by its share of errors that can occur during template or library preparation, sequencing, imaging or data analysis. In this study we screened over 18,000 publicly available microbial isolate genome sequences in the Integrated Microbial Genomes database and identified more than 1000 genomes that are contaminated with PhiX, a control frequently used during Illumina sequencing runs. Approximately 10% of these genomes have been published in literature and 129 contaminated genomes were sequenced under the Human Microbiome Project. Raw sequence reads are prone to contamination from various sources and are usually eliminated during downstream quality control steps. Detection of PhiX contaminated genomes indicates a lapse in either the application or effectiveness of proper quality control measures. The presence of PhiX contamination in several publicly available isolate genomes can result in additional errors when such data are used in comparative genomics analyses. Such contamination of public databases have far-reaching consequences in the form of erroneous data interpretation and analyses, and necessitates better measures to proofread raw sequences before releasing them to the broader scientific community.

Sequencing and characterizing the genome of Estrella lausannensis as an undergraduate project: training students and biological insights

With the widespread availability of high-throughput sequencing technologies, sequencing projects have become pervasive in the molecular life sciences. The huge bulk of data generated daily must be analyzed further by biologists with skills in bioinformatics and by "embedded bioinformaticians," i.e., bioinformaticians integrated in wet lab research groups. Thus, students interested in molecular life sciences must be trained in the main steps of genomics: sequencing, assembly, annotation and analysis. To reach that goal, a practical course has been set up for master students at the University of Lausanne: the "Sequence a genome" class. At the beginning of the academic year, a few bacterial species whose genome is unknown are provided to the students, who sequence and assemble the genome(s) and perform manual annotation. Here, we report the progress of the first class from September 2010 to June 2011 and the results obtained by seven master students who specifically assembled and annotated the genome of Estrella lausannensis, an obligate intracellular bacterium related to Chlamydia. The draft genome of Estrella is composed of 29 scaffolds encompassing 2,819,825 bp that encode for 2233 putative proteins. Estrella also possesses a 9136 bp plasmid that encodes for 14 genes, among which we found an integrase and a toxin/antitoxin module. Like all other members of the Chlamydiales order, Estrella possesses a highly conserved type III secretion system, considered as a key virulence factor. The annotation of the Estrella genome also allowed the characterization of the metabolic abilities of this strictly intracellular bacterium. Altogether, the students provided the scientific community with the Estrella genome sequence and a preliminary understanding of the biology of this recently-discovered bacterial genus, while learning to use cutting-edge technologies for sequencing and to perform bioinformatics analyses.

Bioinformatics approaches for viral metagenomics in plants using short RNAs: model case of study and application to a Cicer arietinum population

Over the past years deep sequencing experiments have opened novel doors to reconstruct viral populations in a high-throughput and cost-effective manner. Currently a substantial number of studies have been performed which employ next generation sequencing techniques to either analyze known viruses by means of a reference-guided approach or to discover novel viruses using a de novo-based strategy. Taking advantage of the well-known Cymbidium ringspot virus we have carried out a comparison of different bioinformatics tools to reconstruct the viral genome based on 21–27 nt short (s)RNA sequencing with the aim to identify the most efficient pipeline. The same approach was applied to a population of plants constituting an ancient variety of Cicer arietinum with red seeds. Among the discovered viruses, we describe the presence of a Tobamovirus referring to the Tomato mottle mosaic virus (NC_022230), which was not yet observed on C. arietinum nor revealed in Europe and a viroid referring to Hop stunt viroid (NC_001351.1) never reported in chickpea. Notably, a reference sequence guided approach appeared the most efficient in such kind of investigation. Instead, the de novo assembly reached a non-appreciable coverage although the most prominent viral species could still be identified. Advantages and limitations of viral metagenomics analysis using sRNAs are discussed.

Development and validation of an rDNA operon based primer walking strategy applicable to de novo bacterial genome finishing

Advances in sequencing technology have drastically increased the depth and feasibility of bacterial genome sequencing. However, little information is available that details the specific techniques and procedures employed during genome sequencing despite the large numbers of published genomes. Shotgun approaches employed by second-generation sequencing platforms has necessitated the development of robust bioinformatics tools for in silico assembly, and complete assembly is limited by the presence of repetitive DNA sequences and multi-copy operons. Typically, re-sequencing with multiple platforms and laborious, targeted Sanger sequencing are employed to finish a draft bacterial genome. Here we describe a novel strategy based on the identification and targeted sequencing of repetitive rDNA operons to expedite bacterial genome assembly and finishing. Our strategy was validated by finishing the genome of Paenibacillus polymyxa strain CR1, a bacterium with potential in sustainable agriculture and bio-based processes. An analysis of the 38 contigs contained in the P. polymyxa strain CR1 draft genome revealed 12 repetitive rDNA operons with varied intragenic and flanking regions of variable length, unanimously located at contig boundaries and within contig gaps. These highly similar but not identical rDNA operons were experimentally verified and sequenced simultaneously with multiple, specially designed primer sets. This approach also identified and corrected significant sequence rearrangement generated during the initial in silico assembly of sequencing reads. Our approach reduces the required effort associated with blind primer walking for contig assembly, increasing both the speed and feasibility of genome finishing. Our study further reinforces the notion that repetitive DNA elements are major limiting factors for genome finishing. Moreover, we provided a step-by-step workflow for genome finishing, which may guide future bacterial genome finishing projects.

A new era in palaeomicrobiology: prospects for ancient dental calculus as a long-term record of the human oral microbiome

The field of palaeomicrobiology is dramatically expanding thanks to recent advances in high-throughput biomolecular sequencing, which allows unprecedented access to the evolutionary history and ecology of human-associated and environmental microbes. Recently, human dental calculus has been shown to be an abundant, nearly ubiquitous, and long-term reservoir of the ancient oral microbiome, preserving not only microbial and host biomolecules but also dietary and environmental debris. Modern investigations of native human microbiota have demonstrated that the human microbiome plays a central role in health and chronic disease, raising questions about changes in microbial ecology, diversity and function through time. This paper explores the current state of ancient oral microbiome research and discusses successful applications, methodological challenges and future possibilities in elucidating the intimate evolutionary relationship between humans and their microbes.

Detailed investigation of the microbial community in foaming activated sludge reveals novel foam formers

Foaming of activated sludge (AS) causes adverse impacts on wastewater treatment operation and hygiene. In this study, we investigated the microbial communities of foam, foaming AS and non-foaming AS in a sewage treatment plant via deep-sequencing of the taxonomic marker genes 16S rRNA and mycobacterial rpoB and a metagenomic approach. In addition to Actinobacteria, many genera (e.g., Clostridium XI, Arcobacter, Flavobacterium) were more abundant in the foam than in the AS. On the other hand, deep-sequencing of rpoB did not detect any obligate pathogenic mycobacteria in the foam. We found that unknown factors other than the abundance of Gordonia sp. could determine the foaming process, because abundance of the same species was stable before and after a foaming event over six months. More interestingly, although the dominant Gordonia foam former was the closest with G. amarae, it was identified as an undescribed Gordonia species by referring to the 16S rRNA gene, gyrB and, most convincingly, the reconstructed draft genome from metagenomic reads. Our results, based on metagenomics and deep sequencing, reveal that foams are derived from diverse taxa, which expands previous understanding and provides new insight into the underlying complications of the foaming phenomenon in AS.

Microbial experimental evolution as a novel research approach in the Vibrionaceae and squid-Vibrio symbiosis

The Vibrionaceae are a genetically and metabolically diverse family living in aquatic habitats with a great propensity toward developing interactions with eukaryotic microbial and multicellular hosts (as either commensals, pathogens, and mutualists). The Vibrionaceae frequently possess a life history cycle where bacteria are attached to a host in one phase and then another where they are free from their host as either part of the bacterioplankton or adhered to solid substrates such as marine sediment, riverbeds, lakebeds, or floating particulate debris. These two stages in their life history exert quite distinct and separate selection pressures. When bound to solid substrates or to host cells, the Vibrionaceae can also exist as complex biofilms. The association between bioluminescent Vibrio spp. and sepiolid squids (Cephalopoda: Sepiolidae) is an experimentally tractable model to study bacteria and animal host interactions, since the symbionts and squid hosts can be maintained in the laboratory independently of one another. The bacteria can be grown in pure culture and the squid hosts raised gnotobiotically with sterile light organs. The partnership between free-living Vibrio symbionts and axenic squid hatchlings emerging from eggs must be renewed every generation of the cephalopod host. Thus, symbiotic bacteria and animal host can each be studied alone and together in union. Despite virtues provided by the Vibrionaceae and sepiolid squid-Vibrio symbiosis, these assets to evolutionary biology have yet to be fully utilized for microbial experimental evolution. Experimental evolution studies already completed are reviewed, along with exploratory topics for future study.

Comparative genomics of first available bovine Anaplasma phagocytophilum genome obtained with targeted sequence capture

BACKGROUND

Anaplasma phagocytophilum is a zoonotic and obligate intracellular bacterium transmitted by ticks. In domestic ruminants, it is the causative agent of tick-borne fever, which causes significant economic losses in Europe. As A. phagocytophilum is difficult to isolate and cultivate, only nine genome sequences have been published to date, none of which originate from a bovine strain.Our goals were to; 1/ develop a sequencing methodology which efficiently circumvents the difficulties associated with A. phagocytophilum isolation and culture; 2/ describe the first genome of a bovine strain; and 3/ compare it with available genomes, in order to both explore key genomic features at the species level, and to identify candidate genes that could be specific to bovine strains.

RESULTS

DNA was extracted from a bovine blood sample infected by A. phagocytophilum. Following a whole genome capture approach, A. phagocytophilum DNA was enriched 197-fold in the sample and then sequenced using Illumina technology. In total, 58.9% of obtained reads corresponded to the A. phagocytophilum genome, covering 85.3% of the HZ genome. Then by performing comparisons with nine previously-sequenced A. phagocytophilum genomes, we determined the core genome of these ten strains. Following analysis, 1281 coding DNA sequences, including 1001 complete sequences, were detected in the A. phagocytophilum bovine genome, of which four appeared to be unique to the bovine isolate. These four coding DNA sequences coded for "hypothetical proteins of unknown function" and require further analysis. We also identified nine proteins common to both European domestic ruminants tested.

CONCLUSION

Using a whole genome capture approach, we have sequenced the first A. phagocytophilum genome isolated from a cow. To the best of our knowledge, this is the first time that this method has been used to selectively enrich pathogenic bacterial DNA from samples also containing host DNA. The four proteins unique to the A. phagocytophilum bovine genome could be involved in host tropism, therefore their functions need to be explored.

Seasonal variability in airborne biotic contaminants in swine confinement buildings

Little is known about the seasonal dynamics of biotic contaminants in swine confinement buildings (SCBs). The biotic contaminants of seven SCBs were monitored during one visit in the winter and one during the summer. Paired-end Illumina sequencing of the 16S rRNA gene, V3 region, was used to examine seasonal shifts in bacterial community composition and diversity. The abundances of 16S rRNA genes and six tetracycline resistance genes (tetB, tetH, tetZ, tetO, tetQ, and tetW) were also quantified using real-time PCR. Bacterial abundances, community composition and diversity all showed strong seasonal patterns defined by winter peaks in abundance and diversity. Microclimatic variables of SCBs, particularly air speed, PM2.5 and total suspended particles (TSP) were found significantly correlated to abundances, community composition, and diversity of bacterial bioaerosols. Seasonal fluctuations were also observed for four tetracycline resistance genes, tetH, tetO, tetQ, and tetW. The frequency of occurrences of these resistance genes were significantly higher in samples collected during winter and was also significantly correlated with air speed, PM2.5 and TSP. Overall, our results indicate that biotic contaminants in SCBs exhibit seasonal trends, and these could be associated with the microclimatic variables of SCBs. The correlations established in the current study could be helpful in establishing better management strategies to minimize the potential health impacts on both livestock and humans working in this environment.

Platinum: a database of experimentally measured effects of mutations on structurally defined protein-ligand complexes

Drug resistance is a major challenge for the treatment of many diseases and a significant concern throughout the drug development process. The ability to understand and predict the effects of mutations on protein–ligand affinities and their roles in the emergence of resistance would significantly aid treatment and drug design strategies. In order to study and understand the impacts of missense mutations on the interaction of ligands with the proteome, we have developed Platinum (http://structure.bioc.cam.ac.uk/platinum). This manually curated, literature-derived database, comprising over 1000 mutations, associates for the first time experimental information on changes in affinity with three-dimensional structures of protein–ligand complexes. To minimize differences arising from experimental techniques and to directly compare binding affinities, Platinum considers only changes measured by the same group and with the same amino-acid sequence used for structure determination, providing a direct link between protein structure, how a ligand binds and how mutations alter the affinity of the ligand of the protein. We believe Platinum will be an invaluable resource for understanding the effects of mutations that give rise to drug resistance, a major problem emerging in pandemics including those caused by the influenza virus, in infectious diseases such as tuberculosis, in cancer and in many other life-threatening illnesses.

Whole genome sequencing and analysis of plant growth promoting bacteria isolated from the rhizosphere of plantation crops coconut, cocoa and arecanut

Coconut, cocoa and arecanut are commercial plantation crops that play a vital role in the Indian economy while sustaining the livelihood of more than 10 million Indians. According to 2012 Food and Agricultural organization's report, India is the third largest producer of coconut and it dominates the production of arecanut worldwide. In this study, three Plant Growth Promoting Rhizobacteria (PGPR) from coconut (CPCRI-1), cocoa (CPCRI-2) and arecanut (CPCRI-3) characterized for the PGP activities have been sequenced. The draft genome sizes were 4.7 Mb (56% GC), 5.9 Mb (63.6% GC) and 5.1 Mb (54.8% GB) for CPCRI-1, CPCRI-2, CPCRI-3, respectively. These genomes encoded 4056 (CPCRI-1), 4637 (CPCRI-2) and 4286 (CPCRI-3) protein-coding genes. Phylogenetic analysis revealed that both CPCRI-1 and CPCRI-3 belonged to Enterobacteriaceae family, while, CPCRI-2 was a Pseudomonadaceae family member. Functional annotation of the genes predicted that all three bacteria encoded genes needed for mineral phosphate solubilization, siderophores, acetoin, butanediol, 1-aminocyclopropane-1-carboxylate (ACC) deaminase, chitinase, phenazine, 4-hydroxybenzoate, trehalose and quorum sensing molecules supportive of the plant growth promoting traits observed in the course of their isolation and characterization. Additionally, in all the three CPCRI PGPRs, we identified genes involved in synthesis of hydrogen sulfide (H2S), which recently has been proposed to aid plant growth. The PGPRs also carried genes for central carbohydrate metabolism indicating that the bacteria can efficiently utilize the root exudates and other organic materials as energy source. Genes for production of peroxidases, catalases and superoxide dismutases that confer resistance to oxidative stresses in plants were identified. Besides these, genes for heat shock tolerance, cold shock tolerance and glycine-betaine production that enable bacteria to survive abiotic stress were also identified.

Infectious disease, shifting climates, and opportunistic predators: cumulative factors potentially impacting wild salmon declines

Emerging diseases are impacting animals under high-density culture, yet few studies assess their importance to wild populations. Microparasites selected for enhanced virulence in culture settings should be less successful maintaining infectivity in wild populations, as once the host dies, there are limited opportunities to infect new individuals. Instead, moderately virulent microparasites persisting for long periods across multiple environments are of greatest concern. Evolved resistance to endemic microparasites may reduce susceptibilities, but as barriers to microparasite distributions are weakened, and environments become more stressful, unexposed populations may be impacted and pathogenicity enhanced. We provide an overview of the evolutionary and ecological impacts of infectious diseases in wild salmon and suggest ways in which modern technologies can elucidate the microparasites of greatest potential import. We present four case studies that resolve microparasite impacts on adult salmon migration success, impact of river warming on microparasite replication, and infection status on susceptibility to predation. Future health of wild salmon must be considered in a holistic context that includes the cumulative or synergistic impacts of multiple stressors. These approaches will identify populations at greatest risk, critically needed to manage and potentially ameliorate the shifts in current or future trajectories of wild populations.

Validation of high throughput sequencing and microbial forensics applications

High throughput sequencing (HTS) generates large amounts of high quality sequence data for microbial genomics. The value of HTS for microbial forensics is the speed at which evidence can be collected and the power to characterize microbial-related evidence to solve biocrimes and bioterrorist events. As HTS technologies continue to improve, they provide increasingly powerful sets of tools to support the entire field of microbial forensics. Accurate, credible results allow analysis and interpretation, significantly influencing the course and/or focus of an investigation, and can impact the response of the government to an attack having individual, political, economic or military consequences. Interpretation of the results of microbial forensic analyses relies on understanding the performance and limitations of HTS methods, including analytical processes, assays and data interpretation. The utility of HTS must be defined carefully within established operating conditions and tolerances. Validation is essential in the development and implementation of microbial forensics methods used for formulating investigative leads attribution. HTS strategies vary, requiring guiding principles for HTS system validation. Three initial aspects of HTS, irrespective of chemistry, instrumentation or software are: 1) sample preparation, 2) sequencing, and 3) data analysis. Criteria that should be considered for HTS validation for microbial forensics are presented here. Validation should be defined in terms of specific application and the criteria described here comprise a foundation for investigators to establish, validate and implement HTS as a tool in microbial forensics, enhancing public safety and national security.

Whole-genome sequencing of clarithromycin resistant Helicobacter pylori characterizes unidentified variants of multidrug resistant efflux pump genes

BACKGROUND

Clarithromycin (CLR) is the key drug in eradication therapy of Helicobacter pylori (H. pylori) infection, and widespread use of CLR has led to an increase in primary CLR-resistant H. pylori. The known mechanism of CLR resistance has been established in A2146G and A2147G mutations in the 23S rRNA gene, but evidence of the involvement of other genetic mechanisms is lacking. Using the MiSeq platform, whole-genome sequencing of the 19 clinical strains and the reference strain ATCC26695 was performed to identify single nucleotide variants (SNVs) of multi-drug resistant efflux pump genes in the CLR-resistant phenotype.

RESULTS

Based on sequencing data of ATCC26695, over one million sequencing reads with over 50-fold coverage were sufficient to detect SNVs, but not indels in the bacterial genome. Sequencing reads of the clinical isolates ranged from 1.82 to 10.8 million, and average coverage ranged from 90.9- to 686.3-fold, which were acceptable criteria for detecting SNVs. Utilizing the conventional approach of allele-specific PCR, point mutations in the 23S rRNA gene were detected in 12 clinical resistant isolates, but not in 7 clinical susceptible isolates. All sequencing reads of CLR-resistant strains had a G mutation in an identical position of the 23S rRNA gene. In addition, genetic variants of four gene clusters (hp0605-hp0607, hp0971-hp0969, hp1327-hp1329, and hp1489-hp1487) of TolC homologues, which have been implicated in multi-drug resistance, were examined. Specific SNVs were dominantly found in resistant strains.

CONCLUSIONS

Gene clusters of TolC homologues are involved in CLR susceptibility profiles in individual H. pylori strains. Whole-genome sequencing has yielded novel understanding of genotype-phenotype relationships.

Mycobacterial phylogenomics: an enhanced method for gene turnover analysis reveals uneven levels of gene gain and loss among species and gene families

Species of the genus Mycobacterium differ in several features, from geographic ranges, and degree of pathogenicity, to ecological and host preferences. The recent availability of several fully sequenced genomes for a number of these species enabled the comparative study of the genetic determinants of this wide lifestyle diversity. Here, we applied two complementary phylogenetic-based approaches using information from 19 Mycobacterium genomes to obtain a more comprehensive view of the evolution of this genus. First, we inferred the phylogenetic relationships using two new approaches, one based on a Mycobacterium-specific amino acid substitution matrix and the other on a gene content dissimilarity matrix. Then, we utilized our recently developed gain-and-death stochastic models to study gene turnover dynamics in this genus in a maximum-likelihood framework. We uncovered a scenario that differs markedly from traditional 16S rRNA data and improves upon recent phylogenomic approaches. We also found that the rates of gene gain and death are high and unevenly distributed both across species and across gene families, further supporting the utility of the new models of rate heterogeneity applied in a phylogenetic context. Finally, the functional annotation of the most expanded or contracted gene families revealed that the transposable elements and the fatty acid metabolism-related gene families are the most important drivers of gene content evolution in Mycobacterium.

Deep mRNA sequencing reveals stage-specific transcriptome alterations during microsclerotia development in the smoke tree vascular wilt pathogen, Verticillium dahliae

Background

Verticillium dahliae is a soil-borne fungus that causes vascular wilt diseases in a wide range of plant hosts. V. dahliae produces multicelled, melanized resting bodies, also known as microsclerotia (MS) that can survive for years in the soil. The MS are the primary source of infection of the Verticillium disease cycle. Thus, MS formation marks an important event in the disease cycle of V. dahliae.

Results

In this study, next generation sequencing technology of RNA-Seq was employed to investigate the global transcriptomic dynamics of MS development to identify differential gene expression at several stages of MS formation in strain XS11 of V. dahliae, isolated from smoke tree. We observed large-scale changes in gene expression during MS formation, such as increased expression of genes involved in protein metabolism and carbohydrate metabolism. Genes involved in glycolytic pathway and melanin biosynthesis were dramatically up-regulated in MS. Cluster analyses revealed increased expression of genes encoding products involved in primary metabolism and stress responses throughout MS development. Differential expression of ubiquitin-dependent protein catabolism and cell death-associated genes during MS development were revealed. Homologs of genes located in the lineage-specific (LS) regions of V. dahliae strain VdLs.17, were either not expressed or showed low expression. Furthermore, alternative splicing (AS) events were analyzed, revealing that over 95.0% AS events involve retention of introns (RI).

Conclusions

These data reveal the dynamics of transcriptional regulation during MS formation and were used to construct a comprehensive high-resolution gene expression map. This map provides a key resource for understanding the biology and molecular basis of MS development of V. dahliae.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-324) contains supplementary material, which is available to authorized users.