ContigScape: a Cytoscape plugin facilitating microbial genome gap closing

Complete genome sequence of a novel bovine hepacivirus from Yunnan, China

We detected a novel bovine hepacivirus N (HNV) subtype, IME_BovHep_01, in the serum of cattle in Tengchong, Yunnan, China, by high-throughput sequencing. The complete genome of IME_BovHep_01, was sequenced using an Illumina MiSeq sequencer and found to be 8850 nt in length, encoding one hypothetical protein. BLASTn analysis showed that the genome sequence shared similarity with the bovine hepacivirus isolate BovHepV_209/Ger/2014, with 88% query coverage and 70.8% identity. However, the highest similarity was to bovine hepacivirus N strain BRBovHep_RS963, for which only a partial genome sequence is available, with 68% query coverage and 81.5% identity. Sequence comparisons and phylogenetic analysis suggested that IME_BovHep_01 is a novel HNV subtype. Importantly, IME_BovHep_01 is the first member of this new genotype for which the complete genome sequence was determined.

Integrated identification of key genes and pathways in Alzheimer's disease via comprehensive bioinformatical analyses

Background

Alzheimer's disease (AD) is known to be caused by multiple factors, meanwhile the pathogenic mechanism and development of AD associate closely with genetic factors. Existing understanding of the molecular mechanisms underlying AD remains incomplete.

Methods

Gene expression data (GSE48350) derived from post-modern brain was extracted from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) were derived from hippocampus and entorhinal cortex regions between AD patients and healthy controls and detected via Morpheus. Functional enrichment analyses, including Gene Ontology (GO) and pathway analyses of DEGs, were performed via Cytoscape and followed by the construction of protein-protein interaction (PPI) network. Hub proteins were screened using the criteria: nodes degree≥10 (for hippocampus tissues) and ≥ 8 (for entorhinal cortex tissues). Molecular Complex Detection (MCODE) was used to filtrate the important clusters. University of California Santa Cruz (UCSC) and the database of RNA-binding protein specificities (RBPDB) were employed to identify the RNA-binding proteins of the long non-coding RNA (lncRNA).

Results

251 & 74 genes were identified as DEGs, which consisted of 56 & 16 up-regulated genes and 195 & 58 down-regulated genes in hippocampus and entorhinal cortex, respectively. Biological analyses demonstrated that the biological processes and pathways related to memory, transmembrane transport, synaptic transmission, neuron survival, drug metabolism, ion homeostasis and signal transduction were enriched in these genes. 11 genes were identified as hub genes in hippocampus and entorhinal cortex, and 3 hub genes were identified as the novel candidates involved in the pathology of AD. Furthermore, 3 lncRNAs were filtrated, whose binding proteins were closely associated with AD.

Conclusions

Through GO enrichment analyses, pathway analyses and PPI analyses, we showed a comprehensive interpretation of the pathogenesis of AD at a systematic biology level, and 3 novel candidate genes and 3 lncRNAs were identified as novel and potential candidates participating in the pathology of AD. The results of this study could supply integrated insights for understanding the pathogenic mechanism underlying AD and potential novel therapeutic targets.

Genomic Characterization Provides New Insights for Detailed Phage- Resistant Mechanism for Brucella abortus

As the causative agent of cattle brucellosis, Brucella abortus commonly exhibits smooth phenotype (by virtue of colony morphology) that is characteristically sensitive to specific Brucella phages, playing until recently a major role in taxonomical classification of the Brucella species by the phage typing approach. We previously reported the discrepancy between traditional phenotypic typing and MLVA results of a smooth phage-resistant (SPR) strain Bab8416 isolated from a 45-year-old custodial worker with brucellosis in a cattle farm. Here, we performed whole genome sequencing and further obtained a complete genome sequence of strain Bab8416 by a combination of multiple NGS technologies and routine PCR sequencing. The detailed genetic differences between B. abortus SPR Bab8416 and large smooth phage-sensitive (SPS) strains were investigated in a comprehensively comparative genomic study. The large indels between B. abortus SPS strains and Bab8416 showed possible divergence between two evolutionary branches at a far phylogenetic node. Compared to B. abortus SPS strain 9-941 (Bab9-941), the specific re-arrangement event in Bab8416 displaying a closer linear relationship with B. melitensis 16M than other B. abortus strains resulted in the truncation of c-di-GMP synthesis, and 3 c-di-GMP-metabolizing genes, were present in Bab8416 and B. melitensis 16M, but absent in Bab9-941 and other B. abortus strains, indicating potential SPR-associated key determinants and novel molecular mechanisms. Moreover, despite almost completely intact smooth LPS related genes, only one mutated OmpA family protein of Bab8416, functionally related to flagellar and efflux pump, was newly identified. Several point mutations were identified to be Bab8416 specific while a majority of them were verified to be B. abortus ST2 characteristic. In conclusion, our study therefore identifies new SPR-associated factors that could play a role in refining and updating Brucella taxonomic schemes and provides resources for further detailed analysis of mechanism for Brucella phage resistance.

Distribution and characteristics of SGI1/PGI2 genomic island from Proteus strains in China

The emergence of multidrug-resistant Salmonella genomic island 1 (SGI1) and Proteus genomic island (PGI) bearing P. mirabilis present a serious threat to public health. In this study, we screened 288 Proteus isolates recovered from seven provinces in China. Fourteen strains (4.9%) all belonged to P. mirabilis were positive for SGI1/PGI2, including twelve from clinical samples (5.3%) and two from food (3.3%). A Blastn search against GenBank and phylogenetic analyses identified eight different SGI1 variants and one PGI2 variant from the fourteen SGI1/PGI2 variants. All SGI1 variants shared a common backbone and harbored different resistance gene(s), except the sul1 gene at its multidrug-resistant (MDR) region. Among the variants, three novel SGI1 variants, designated as SGI1-PmCA11, SGI1-PmCA14 and SGI1-PmCA46, contained different gene cassettes, which were similar to sequences in plasmids or class 1 integrons of Klebsiella pneumoniae, P. mirabilis, Escherichia coli and Salmonella. Moreover, one novel PGI2, designated as PGI2-PmCA72, had an identical gene cassette to the first class 1 integron from PGI2 (GenBank accession no. MG201402.1) in P. mirabilis, but varied due to missing, replaced, inserted and inverted gene clusters. The four novel SGI1/PGI2 variants contained the cmlA5, dfrA14, bla_OXA-10, aadA15, bla_OXA-1, catB3 and dfrA16 resistance genes, which have never been reported in SGI1/PGI2 variants. Phenotypically, all fourteen SGI1/PGI2-containing strains showed multidrug resistance. All except four strains were resistant to the first, or the second and/or-third generation cephalosporins. Considering the increasing number and the emergence of new SGI1/PGI2 variants, further surveillance is needed to prevent the spreading of the MDR genomic islands among Proteus isolates from human and food.

Comprehensive Comparative Analysis of Cholesterol Catabolic Genes/Proteins in Mycobacterial Species

In dealing with Mycobacterium tuberculosis, the causative agent of the deadliest human disease—tuberculosis (TB)—utilization of cholesterol as a carbon source indicates the possibility of using cholesterol catabolic genes/proteins as novel drug targets. However, studies on cholesterol catabolism in mycobacterial species are scarce, and the number of mycobacterial species utilizing cholesterol as a carbon source is unknown. The availability of a large number of mycobacterial species’ genomic data affords an opportunity to explore and predict mycobacterial species’ ability to utilize cholesterol employing in silico methods. In this study, comprehensive comparative analysis of cholesterol catabolic genes/proteins in 93 mycobacterial species was achieved by deducing a comprehensive cholesterol catabolic pathway, developing a software tool for extracting homologous protein data and using protein structure and functional data. Based on the presence of cholesterol catabolic homologous proteins proven or predicted to be either essential or specifically required for the growth of M. tuberculosis H37Rv on cholesterol, we predict that among 93 mycobacterial species, 51 species will be able to utilize cholesterol as a carbon source. This study’s predictions need further experimental validation and the results should be taken as a source of information on cholesterol catabolism and genes/proteins involved in this process among mycobacterial species.

LR_Gapcloser: a tiling path-based gap closer that uses long reads to complete genome assembly

Background

Completing a genome is an important goal of genome assembly. However, many assemblies, including reference assemblies, are unfinished and have a number of gaps. Long reads obtained from third-generation sequencing (TGS) platforms can help close these gaps and improve assembly contiguity. However, current gap-closure approaches using long reads require extensive runtime and high memory usage. Thus, a fast and memory-efficient approach using long reads is needed to obtain complete genomes.

Findings

We developed LR_Gapcloser to rapidly and efficiently close the gaps in genome assembly. This tool utilizes long reads generated from TGS sequencing platforms. Tested on de novo assembled gaps, repeat-derived gaps, and real gaps, LR_Gapcloser closed a higher number of gaps faster and with a lower error rate and a much lower memory usage than two existing, state-of-the art tools. This tool utilized raw reads to fill more gaps than when using error-corrected reads. It is applicable to gaps in the assemblies by different approaches and from large and complex genomes. After performing gap-closure using this tool, the contig N50 size of the human CHM1 genome was improved from 143 kb to 19 Mb, a 132-fold increase. We also closed the gaps in the Triticum urartu genome, a large genome rich in repeats; the contig N50 size was increased by 40%. Further, we evaluated the contiguity and correctness of six hybrid assembly strategies by combining the optimal TGS-based and next-generation sequencing-based assemblers with LR_Gapcloser. A proposed and optimal hybrid strategy generated a new human CHM1 genome assembly with marked contiguity. The contig N50 value was greater than 28 Mb, which is larger than previous non-reference assemblies of the diploid human genome.

Conclusions

LR_Gapcloser is a fast and efficient tool that can be used to close gaps and improve the contiguity of genome assemblies. A proposed hybrid assembly including this tool promises reference-grade assemblies. The software is available at http://www.fishbrowser.org/software/LR_Gapcloser/.

LR_Gapcloser: a tiling path-based gap closer that uses long reads to complete genome assembly

BACKGROUND

Completing a genome is an important goal of genome assembly. However, many assemblies, including reference assemblies, are unfinished and have a number of gaps. Long reads obtained from third-generation sequencing (TGS) platforms can help close these gaps and improve assembly contiguity. However, current gap-closure approaches using long reads require extensive runtime and high memory usage. Thus, a fast and memory-efficient approach using long reads is needed to obtain complete genomes.

FINDINGS

We developed LR_Gapcloser to rapidly and efficiently close the gaps in genome assembly. This tool utilizes long reads generated from TGS sequencing platforms. Tested on de novo assembled gaps, repeat-derived gaps, and real gaps, LR_Gapcloser closed a higher number of gaps faster and with a lower error rate and a much lower memory usage than two existing, state-of-the art tools. This tool utilized raw reads to fill more gaps than when using error-corrected reads. It is applicable to gaps in the assemblies by different approaches and from large and complex genomes. After performing gap-closure using this tool, the contig N50 size of the human CHM1 genome was improved from 143 kb to 19 Mb, a 132-fold increase. We also closed the gaps in the Triticum urartu genome, a large genome rich in repeats; the contig N50 size was increased by 40%. Further, we evaluated the contiguity and correctness of six hybrid assembly strategies by combining the optimal TGS-based and next-generation sequencing-based assemblers with LR_Gapcloser. A proposed and optimal hybrid strategy generated a new human CHM1 genome assembly with marked contiguity. The contig N50 value was greater than 28 Mb, which is larger than previous non-reference assemblies of the diploid human genome.

CONCLUSIONS

LR_Gapcloser is a fast and efficient tool that can be used to close gaps and improve the contiguity of genome assemblies. A proposed hybrid assembly including this tool promises reference-grade assemblies. The software is available at http://www.fishbrowser.org/software/LR_Gapcloser/.

LR_Gapcloser: a tiling path-based gap closer that uses long reads to complete genome assembly

BACKGROUND

Completing a genome is an important goal of genome assembly. However, many assemblies, including reference assemblies, are unfinished and have a number of gaps. Long reads obtained from third-generation sequencing (TGS) platforms can help close these gaps and improve assembly contiguity. However, current gap-closure approaches using long reads require extensive runtime and high memory usage. Thus, a fast and memory-efficient approach using long reads is needed to obtain complete genomes.

FINDINGS

We developed LR_Gapcloser to rapidly and efficiently close the gaps in genome assembly. This tool utilizes long reads generated from TGS sequencing platforms. Tested on de novo assembled gaps, repeat-derived gaps, and real gaps, LR_Gapcloser closed a higher number of gaps faster and with a lower error rate and a much lower memory usage than two existing, state-of-the art tools. This tool utilized raw reads to fill more gaps than when using error-corrected reads. It is applicable to gaps in the assemblies by different approaches and from large and complex genomes. After performing gap-closure using this tool, the contig N50 size of the human CHM1 genome was improved from 143 kb to 19 Mb, a 132-fold increase. We also closed the gaps in the Triticum urartu genome, a large genome rich in repeats; the contig N50 size was increased by 40%. Further, we evaluated the contiguity and correctness of six hybrid assembly strategies by combining the optimal TGS-based and next-generation sequencing-based assemblers with LR_Gapcloser. A proposed and optimal hybrid strategy generated a new human CHM1 genome assembly with marked contiguity. The contig N50 value was greater than 28 Mb, which is larger than previous non-reference assemblies of the diploid human genome.

CONCLUSIONS

LR_Gapcloser is a fast and efficient tool that can be used to close gaps and improve the contiguity of genome assemblies. A proposed hybrid assembly including this tool promises reference-grade assemblies. The software is available at http://www.fishbrowser.org/software/LR_Gapcloser/.

Reconstruction and analysis of carbon metabolic pathway of Ketogulonicigenium vulgare SPU B805 by genome and transcriptome

Ketogulonicigenium vulgare has been widely used in vitamin C two-step fermentation. Four K. vulgare strains (WSH-001, Y25, Hbe602 and SKV) have been completely genome-sequenced, however, less attention was paid to elucidate the reason for the differences in 2-KGA yield on genetic level. Here, a novel K. vulgare SPU B805 with higher 2-keto-L-gulonic acid (2-KGA) yield, was genome-sequenced to confirm harboring one circular chromosome with plasmid free. Comparative genome analyses showed that the absence of plasmid 2 was an important factor for its high 2-KGA productivity. The amino acid biosynthetic pathways in strain SPU B805 are much more complete than those in other K. vulgare strains. Meanwhile, strain SPU B805 harbored a complete PPP and TCA route, as well as a disabled EMP and ED pathway, same as to strain SKV, whereas strain WSH-001, Y25 and Hbe602 harbored complete PPP, ED, TCA pathway and a nonfunctional EMP pathway. The transcriptome of strain SPU B805 validated the carbon metabolism in cytoplasm mainly through the PPP pathway due to its higher transcriptional levels. This is the first time to elucidate the underlying mechanism for the difference in 2-KGA yield, and it is of great significance for strain improvement in the industrial fermentation.

Insights into the ecology, evolution, and metabolism of the widespread Woesearchaeotal lineages

BACKGROUND

As a recently discovered member of the DPANN superphylum, Woesearchaeota account for a wide diversity of 16S rRNA gene sequences, but their ecology, evolution, and metabolism remain largely unknown.

RESULTS

Here, we assembled 133 global clone libraries/studies and 19 publicly available genomes to profile these patterns for Woesearchaeota. Phylogenetic analysis shows a high diversity with 26 proposed subgroups for this recently discovered archaeal phylum, which are widely distributed in different biotopes but primarily in inland anoxic environments. Ecological patterns analysis and ancestor state reconstruction for specific subgroups reveal that oxic status of the environments is the key factor driving the distribution and evolutionary diversity of Woesearchaeota. A selective distribution to different biotopes and an adaptive colonization from anoxic to oxic environments can be proposed and supported by evidence of the presence of ferredoxin-dependent pathways in the genomes only from anoxic biotopes but not from oxic biotopes. Metabolic reconstructions support an anaerobic heterotrophic lifestyle with conspicuous metabolic deficiencies, suggesting the requirement for metabolic complementarity with other microbes. Both lineage abundance distribution and co-occurrence network analyses across diverse biotopes confirmed metabolic complementation and revealed a potential syntrophic relationship between Woesearchaeota and methanogens, which is supported by metabolic modeling. If correct, Woesearchaeota may impact methanogenesis in inland ecosystems.

CONCLUSIONS

The findings provide an ecological and evolutionary framework for Woesearchaeota at a global scale and indicate their potential ecological roles, especially in methanogenesis.

Establishing an innovative carbohydrate metabolic pathway for efficient production of 2-keto-L-gulonic acid in Ketogulonicigenium robustum initiated by intronic promoters

Background

2-Keto-l-gulonic acid (2-KGA), the precursor of vitamin C, is currently produced by two-step fermentation. In the second step, l-sorbose is transformed into 2-KGA by the symbiosis system composed of Ketogulonicigenium vulgare and Bacillus megaterium. Due to the different nutrient requirements and the uncertain ratio of the two strains, the symbiosis system significantly limits strain improvement and fermentation optimization.

Results

In this study, Ketogulonicigenium robustum SPU_B003 was reported for its capability to grow well independently and to produce more 2-KGA than that of K. vulgare in a mono-culture system. The complete genome of K. robustum SPU_B003 was sequenced, and the metabolic characteristics were analyzed. Compared to the four reported K. vulgare genomes, K. robustum SPU_B003 contained more tRNAs, rRNAs, NAD and NADP biosynthetic genes, as well as regulation- and cell signaling-related genes. Moreover, the amino acid biosynthesis pathways were more complete. Two species-specific internal promoters, P1 (orf_01408 promoter) and P2 (orf_02221 promoter), were predicted and validated by detecting their initiation activity. To efficiently produce 2-KGA with decreased CO₂ release, an innovative acetyl-CoA biosynthetic pathway (XFP-PTA pathway) was introduced into K. robustum SPU_B003 by expressing heterologous phosphoketolase (xfp) and phosphotransacetylase (pta) initiated by internal promoters. After gene optimization, the recombinant strain K. robustum/pBBR-P1_xfp2502-P2_pta2145 enhanced acetyl-CoA approximately 2.4-fold and increased 2-KGA production by 22.27% compared to the control strain K. robustum/pBBR1MCS-2. Accordingly, the transcriptional level of the 6-phosphogluconate dehydrogenase (pgd) and pyruvate dehydrogenase genes (pdh) decreased by 24.33 ± 6.67 and 8.67 ± 5.51%, respectively. The key genes responsible for 2-KGA biosynthesis, sorbose dehydrogenase gene (sdh) and sorbosone dehydrogenase gene (sndh), were up-regulated to different degrees in the recombinant strain.

Conclusions

The genome-based functional analysis of K. robustum SPU_B003 provided a new understanding of the specific metabolic characteristics. The new XFP-PTA pathway was an efficient route to enhance acetyl-CoA levels and to therefore promote 2-KGA production.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-0932-9) contains supplementary material, which is available to authorized users.

Natural products discovery from micro-organisms in the post-genome era†

With the decision to award the Nobel Prize in Physiology or Medicine to Drs. S. Ōmura, W.C. Campbell, and Y. Tu, the importance and usefulness of natural drug discovery and development have been revalidated. Since the end of the twentieth century, many genome analyses of organisms have been conducted, and accordingly, numerous microbial genomes have been decoded. In particular, genomic studies of actinomycetes, micro-organisms that readily produce natural products, led to the discovery of biosynthetic gene clusters responsible for producing natural products. New explorations for natural products through a comprehensive approach combining genomic information with conventional methods show great promise for the discovery of new natural products and even systematic generation of unnaturally occurring compounds.

Complete genome sequence of Menghai rhabdovirus, a novel mosquito-borne rhabdovirus from China

Menghai rhabdovirus (MRV) was isolated from Aedes albopictus in Menghai county of Yunnan Province, China, in August 2010. Whole-genome sequencing of MRV was performed using an Ion PGM™ Sequencer. We found that MRV is a single-stranded, negative-sense RNA virus. The complete genome of MRV has 10,744 nt, with short inverted repeat termini, encoding five typical rhabdovirus proteins (N, P, M, G, and L) and an additional small hypothetical protein. Nucleotide BLAST analysis using the BLASTn method showed that the genome sequence most similar to that of MRV is that of Arboretum virus (NC_025393.1), with a Max score of 322, query coverage of 14%, and 66% identity. Genomic and phylogenetic analyses both demonstrated that MRV should be considered a member of a novel species of the family Rhabdoviridae.

SXT/R391 integrative and conjugative elements in Proteus species reveal abundant genetic diversity and multidrug resistance

SXT/R391 integrative and conjugative elements (ICEs) are self-transmissible mobile genetic elements that are found in most members of Enterobacteriaceae. Here, we determined fifteen SXT/R391 ICEs carried by Proteus isolates from food (4.2%) and diarrhoea patients (17.3%). BLASTn searches against GenBank showed that the fifteen SXT/R391 ICEs were closely related to that from different Enterobacteriaceae species, including Proteus mirabilis. Using core gene phylogenetic analysis, the fifteen SXT/R391 ICEs were grouped into six distinct clusters, including a dominant cluster and three clusters that have not been previously reported in Proteus isolates. The SXT/R391 ICEs shared a common structure with a set of conserved genes, five hotspots and two variable regions, which contained more foreign genes, including drug-resistance genes. Notably, a class A β-lactamase gene was identified in nine SXT/R391 ICEs. Collectively, the ICE-carrying isolates carried resistance genes for 20 tested drugs. Six isolates were resistant to chloramphenicol, kanamycin, streptomycin, trimethoprim-sulfamethoxazole, sulfisoxazole and tetracycline, which are drug resistances commonly encoded by ICEs. Our results demonstrate abundant genetic diversity and multidrug resistance of the SXT/R391 ICEs carried by Proteus isolates, which may have significance for public health. It is therefore necessary to continuously monitor the antimicrobial resistance and related mobile elements among Proteus isolates.

A systematic study of the whole genome sequence of Amycolatopsis methanolica strain 239T provides an insight into its physiological and taxonomic properties which correlate with its position in the genus

The complete genome of methanol-utilizing Amycolatopsis methanolica strain 239^T was generated, revealing a single 7,237,391 nucleotide circular chromosome with 7074 annotated protein-coding sequences (CDSs). Comparative analyses against the complete genome sequences of Amycolatopsis japonica strain MG417-CF17^T, Amycolatopsis mediterranei strain U32 and Amycolatopsis orientalis strain HCCB10007 revealed a broad spectrum of genomic structures, including various genome sizes, core/quasi-core/non-core configurations and different kinds of episomes. Although polyketide synthase gene clusters were absent from the A. methanolica genome, 12 gene clusters related to the biosynthesis of other specialized (secondary) metabolites were identified. Complete pathways attributable to the facultative methylotrophic physiology of A. methanolica strain 239^T, including both the mdo/mscR encoded methanol oxidation and the hps/hpi encoded formaldehyde assimilation via the ribulose monophosphate cycle, were identified together with evidence that the latter might be the result of horizontal gene transfer. Phylogenetic analyses based on 16S rDNA or orthologues of AMETH_3452, a novel actinobacterial class-specific conserved gene against 62 or 18 Amycolatopsis type strains, respectively, revealed three major phyletic lineages, namely the mesophilic or moderately thermophilic A. orientalis subclade (AOS), the mesophilic Amycolatopsis taiwanensis subclade (ATS) and the thermophilic A. methanolica subclade (AMS). The distinct growth temperatures of members of the subclades correlated with corresponding genetic variations in their encoded compatible solutes. This study shows the value of integrating conventional taxonomic with whole genome sequence data.

Genome analysis and in vivo virulence of porcine extraintestinal pathogenic Escherichia coli strain PCN033

Background

Strains of extraintestinal pathogenic Escherichia coli (ExPEC) can invade and colonize extraintestinal sites and cause a wide range of infections. Genomic analysis of ExPEC has mainly focused on isolates of human and avian origins, with porcine ExPEC isolates yet to be sequenced. To better understand the genomic attributes underlying the pathogenicity of porcine ExPEC, we isolated two E. coli strains PCN033 and PCN061 from pigs, assessed their in vivo virulence, and completed and compared their genomes.

Results

Animal experiments demonstrated that strain PCN033, but not PCN061, was pathogenic in a pig model. The chromosome of PCN033 was 384 kb larger than that of PCN061. Among the PCN033-specific sequences, genes encoding adhesins, unique lipopolysaccharide, unique capsular polysaccharide, iron acquisition and transport systems, and metabolism were identified. Additionally, a large plasmid PCN033p3 harboring many typical ExPEC virulence factors was identified in PCN033. Based on the genetic variation between PCN033 and PCN061, corresponding phenotypic differences in flagellum-dependent swarming motility and metabolism were verified. Furthermore, the comparative genomic analyses showed that the PCN033 genome shared many similarities with genomic sequences of human ExPEC strains. Additionally, comparison of PCN033 genome with other nine characteristic E. coli genomes revealed 425 PCN033-special coding sequences. Genes of this subset included those encoding type I restriction-modification (R-M) system, type VI secretion system (T6SS) and membrane-associated proteins.

Conclusions

The genetic and phenotypic differences between PCN033 and PCN061 could partially explain their differences in virulence, and also provide insight towards the molecular mechanisms of porcine ExPEC infections. Additionally, the similarities between the genomes of PCN033 and human ExPEC strains suggest that some connections between porcine and human ExPEC strains exist. The first completed genomic sequence for porcine ExPEC and the genomic differences identified by comparative analyses provide a baseline understanding of porcine ExPEC genetics and lay the foundation for their further study.

Electronic supplementary material

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

Metabolic pathway involved in 2-methyl-6-ethylaniline degradation by Sphingobium sp. strain MEA3-1 and cloning of the novel flavin-dependent monooxygenase system meaBA

2-Methyl-6-ethylaniline (MEA) is the main microbial degradation intermediate of the chloroacetanilide herbicides acetochlor and metolachlor. Sphingobium sp. strain MEA3-1 can utilize MEA and various alkyl-substituted aniline and phenol compounds as sole carbon and energy sources for growth. We isolated the mutant strain MEA3-1Mut, which converts MEA only to 2-methyl-6-ethyl-hydroquinone (MEHQ) and 2-methyl-6-ethyl-benzoquinone (MEBQ). MEA may be oxidized by the P450 monooxygenase system to 4-hydroxy-2-methyl-6-ethylaniline (4-OH-MEA), which can be hydrolytically spontaneously deaminated to MEBQ or MEHQ. The MEA microbial metabolic pathway was reconstituted based on the substrate spectra and identification of the intermediate metabolites in both the wild-type and mutant strains. Plasmidome sequencing indicated that both strains harbored 7 plasmids with sizes ranging from 6,108 bp to 287,745 bp. Among the 7 plasmids, 6 were identical, and pMEA02' in strain MEA3-1Mut lost a 37,000-bp fragment compared to pMEA02 in strain MEA3-1. Two-dimensional electrophoresis (2-DE) and protein mass fingerprinting (PMF) showed that MEA3-1Mut lost the two-component flavin-dependent monooxygenase (TC-FDM) MeaBA, which was encoded by a gene in the lost fragment of pMEA02. MeaA shared 22% to 25% amino acid sequence identity with oxygenase components of some TC-FDMs, whereas MeaB showed no sequence identity with the reductase components of those TC-FDMs. Complementation with meaBA in MEA3-1Mut and heterologous expression in Pseudomonas putida strain KT2440 resulted in the production of an active MEHQ monooxygenase.

Visualizing genome and systems biology: technologies, tools, implementation techniques and trends, past, present and future

"Α picture is worth a thousand words." This widely used adage sums up in a few words the notion that a successful visual representation of a concept should enable easy and rapid absorption of large amounts of information. Although, in general, the notion of capturing complex ideas using images is very appealing, would 1000 words be enough to describe the unknown in a research field such as the life sciences? Life sciences is one of the biggest generators of enormous datasets, mainly as a result of recent and rapid technological advances; their complexity can make these datasets incomprehensible without effective visualization methods. Here we discuss the past, present and future of genomic and systems biology visualization. We briefly comment on many visualization and analysis tools and the purposes that they serve. We focus on the latest libraries and programming languages that enable more effective, efficient and faster approaches for visualizing biological concepts, and also comment on the future human-computer interaction trends that would enable for enhancing visualization further.

Complete genome sequence of the glidobactin producing strain [Polyangium] brachysporum DSM 7029

[Polyangium] brachysporum DSM7029 can produce glidobactin, which is a peptide-based proteasome inhibitor that shows great potential as an anticancer drug. To better understand the glidobactin biosynthesis mechanism and to aid further studies of this strain, we report the annotated complete genome sequence of DSM7029, which is 6,476,147 bp in length.

Genomic analysis of a Mycobacterium bovis bacillus [corrected] Calmette-Guérin strain isolated from an adult patient with pulmonary tuberculosis

For years, bacillus Calmette-Guérin (BCG) has served as the unique vaccine against tuberculosis and has generally been regarded as safe. However, a clinical strain labeled 3281 that was isolated from a TB patient was identified to be BCG. Via the combination of next-generation sequencing (NGS) and comparative genomic analysis, unique 3281 genetic characteristics were revealed. A region containing the dnaA and dnaN genes that is closely related to the initial chromosome replication was found to repeat three times on the BCG Pasteur-specific tandem duplication region DU1. Due to the minimum number of epitopes in BCG strains, 3281 was inferred to have a high possibility for immune evasion. Additionally, variations in the virulence genes and predictions for potential virulence factors were analyzed. Overall, we report a pathogen that has never previously been thought to be pathogenic and initial insights that are focused on the genetic characteristics of virulent BCG.

Identification of three extra-chromosomal replicons in Leptospira pathogenic strain and development of new shuttle vectors

Background

The genome of pathogenic Leptospira interrogans contains two chromosomes. Plasmids and prophages are known to play specific roles in gene transfer in bacteria and can potentially serve as efficient genetic tools in these organisms. Although plasmids and prophage remnants have recently been reported in Leptospira species, their characteristics and potential applications in leptospiral genetic transformation systems have not been fully evaluated.

Results

Three extrachromosomal replicons designated lcp1 (65,732 bp), lcp2 (56,757 bp), and lcp3 (54,986 bp) in the L. interrogans serovar Linhai strain 56609 were identified through whole genome sequencing. All three replicons were stable outside of the bacterial chromosomes. Phage particles were observed in the culture supernatant of 56609 after mitomycin C induction, and lcp3, which contained phage-related genes, was considered to be an inducible prophage. L. interrogans–Escherichia coli shuttle vectors, constructed with the predicted replication elements of single rep or rep combined with parAB loci from the three plasmids were shown to successfully transform into both saprophytic and pathogenic Leptospira species, suggesting an essential function for rep genes in supporting auto-replication of the plasmids. Additionally, a wide distribution of homologs of the three rep genes was identified in L. interrogans isolates, and correlation tests showed that the transformability of the shuttle vectors in L. interrogans isolates depended, to certain extent, on genetic compatibility between the rep sequences of both plasmid and host.

Conclusions

Three extrachromosomal replicons co-exist in L. interrogans, one of which we consider to be an inducible prophage. The vectors constructed with the rep genes of the three replicons successfully transformed into saprophytic and pathogenic Leptospira species alike, but this was partly dependent on genetic compatibility between the rep sequences of both plasmid and host.

Electronic supplementary material

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

Genotet: An Interactive Web-based Visual Exploration Framework to Support Validation of Gene Regulatory Networks

Elucidation of transcriptional regulatory networks (TRNs) is a fundamental goal in biology, and one of the most important components of TRNs are transcription factors (TFs), proteins that specifically bind to gene promoter and enhancer regions to alter target gene expression patterns. Advances in genomic technologies as well as advances in computational biology have led to multiple large regulatory network models (directed networks) each with a large corpus of supporting data and gene-annotation. There are multiple possible biological motivations for exploring large regulatory network models, including: validating TF-target gene relationships, figuring out co-regulation patterns, and exploring the coordination of cell processes in response to changes in cell state or environment. Here we focus on queries aimed at validating regulatory network models, and on coordinating visualization of primary data and directed weighted gene regulatory networks. The large size of both the network models and the primary data can make such coordinated queries cumbersome with existing tools and, in particular, inhibits the sharing of results between collaborators. In this work, we develop and demonstrate a web-based framework for coordinating visualization and exploration of expression data (RNA-seq, microarray), network models and gene-binding data (ChIP-seq). Using specialized data structures and multiple coordinated views, we design an efficient querying model to support interactive analysis of the data. Finally, we show the effectiveness of our framework through case studies for the mouse immune system (a dataset focused on a subset of key cellular functions) and a model bacteria (a small genome with high data-completeness).

Complete Genome Sequence of the p-Nitrophenol-Degrading Bacterium Pseudomonas putida DLL-E4

The first complete genome sequence of a p-nitrophenol (PNP)-degrading bacterium is reported here. Pseudomonas putida DLL-E4, a Gram-negative bacterium isolated from methyl-parathion-polluted soil, can utilize PNP as the sole carbon and nitrogen source. P. putida DLL-E4 has a 6,484,062 bp circular chromosome that contains 5,894 genes, with a G+C content of 62.46%.

Complete genome sequence and comparative genomic analyses of the vancomycin-producing Amycolatopsis orientalis

Background

Amycolatopsis orientalis is the type species of the genus and its industrial strain HCCB10007, derived from ATCC 43491, has been used for large-scale production of the vital antibiotic vancomycin. However, to date, neither the complete genomic sequence of this species nor a systemic characterization of the vancomycin biosynthesis cluster (vcm) has been reported. With only the whole genome sequence of Amycolatopsis mediterranei available, additional complete genomes of other species may facilitate intra-generic comparative analysis of the genus.

Results

The complete genome of A. orientalis HCCB10007 comprises an 8,948,591-bp circular chromosome and a 33,499-bp dissociated plasmid. In total, 8,121 protein-coding sequences were predicted, and the species-specific genomic features of A. orientalis were analyzed in comparison with that of A. mediterranei. The common characteristics of Amycolatopsis genomes were revealed via intra- and inter-generic comparative genomic analyses within the domain of actinomycetes, and led directly to the development of sequence-based Amycolatopsis molecular chemotaxonomic characteristics (MCCs). The chromosomal core/quasi-core and non-core configurations of the A. orientalis and the A. mediterranei genome were analyzed reciprocally, with respect to further understanding both the discriminable criteria and the evolutionary implementation. In addition, 26 gene clusters related to secondary metabolism, including the 64-kb vcm cluster, were identified in the genome. Employing a customized PCR-targeting-based mutagenesis system along with the biochemical identification of vancomycin variants produced by the mutants, we were able to experimentally characterize a halogenase, a methyltransferase and two glycosyltransferases encoded in the vcm cluster. The broad substrate spectra characteristics of these modification enzymes were inferred.

Conclusions

This study not only extended the genetic knowledge of the genus Amycolatopsis and the biochemical knowledge of vcm-related post-assembly tailoring enzymes, but also developed methodology useful for in vivo studies in A. orientalis, which has been widely considered as a barrier in this field.

Electronic supplementary material

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

3DScapeCS: application of three dimensional, parallel, dynamic network visualization in Cytoscape

Background

The exponential growth of gigantic biological data from various sources, such as protein-protein interaction (PPI), genome sequences scaffolding, Mass spectrometry (MS) molecular networking and metabolic flux, demands an efficient way for better visualization and interpretation beyond the conventional, two-dimensional visualization tools.

Results

We developed a 3D Cytoscape Client/Server (3DScapeCS) plugin, which adopted Cytoscape in interpreting different types of data, and UbiGraph for three-dimensional visualization. The extra dimension is useful in accommodating, visualizing, and distinguishing large-scale networks with multiple crossed connections in five case studies.

Conclusions

Evaluation on several experimental data using 3DScapeCS and its special features, including multilevel graph layout, time-course data animation, and parallel visualization has proven its usefulness in visualizing complex data and help to make insightful conclusions.