Qiyu Bao

MagicViewer: integrated solution for next-generation sequencing data visualization and genetic variation detection and annotation

New sequencing technologies, such as Roche 454, ABI SOLiD and Illumina, have been increasingly developed at an astounding pace with the advantages of high throughput, reduced time and cost. To satisfy the impending need for deciphering the large-scale data generated from next-generation sequencing, an integrated software MagicViewer is developed to easily visualize short read mapping, identify and annotate genetic variation based on the reference genome. MagicViewer provides a user-friendly environment in which large-scale short reads can be displayed in a zoomable interface under user-defined color scheme through an operating system-independent manner. Meanwhile, it also holds a versatile computational pipeline for genetic variation detection, filtration, annotation and visualization, providing details of search option, functional classification, subset selection, sequence association and primer design. In conclusion, MagicViewer is a sophisticated assembly visualization and genetic variation annotation tool for next-generation sequencing data, which can be widely used in a variety of sequencing-based researches, including genome re-sequencing and transcriptome studies. MagicViewer is freely available at http://bioinformatics.zj.cn/magicviewer/.

Molecular Characterization and Antimicrobial Susceptibility of Nasal Staphylococcus aureus Isolates from a Chinese Medical College Campus

Staphylococcus aureus colonization and infection occur more commonly among persons living or working in crowded conditions, but characterization of S. aureus colonization within medical communities in China is lacking. A total of 144 (15.4%, 144/935) S. aureus isolates, including 28 (3.0%, 28/935) MRSA isolates, were recovered from the nares of 935 healthy human volunteers residing on a Chinese medical college campus. All S. aureus isolates were susceptible to vancomycin, quinupristin/dalfopristin and linezolid but the majority were resistant to penicillin (96.5%), ampicillin/sulbactam (83.3%) and trimethoprim/sulfamethoxazole (93.1%). 82%, (23/28) of the MRSA isolates and 66% (77/116) of the MSSA isolates were resistant to multiple antibiotics, and 3 MRSA isolates were resistant to mupirocin—an agent commonly used for nasal decolonization. 16 different sequence types (STs), as well as SCCmec genes II, III, IVd, and V, were represented among MRSA isolates. We also identified, for the first time, two novel STs (ST1778 and ST1779) and 5 novel spa types for MRSA. MRSA isolates were distributed in different sporadic clones, and ST59-MRSA-VId- t437 was found within 3 MRSA isolates. Moreover, one isolate with multidrug resistance belonging to ST398-MRSA-V- t571 associated with animal infections was identified, and 3 isolates distributed in three different clones harbored PVL genes. Collectively, these data indicate a high prevalence of nasal MRSA carriage and molecular heterogeneity of S. aureus isolates among persons residing on a Chinese medical college campus. Identification of epidemic MRSA clones associated with community infection supports the need for more effective infection control measures to reduce nasal carriage and prevent dissemination of MRSA to hospitalized patients and health care workers in this community.

Sequencing and Genetic Variation of Multidrug Resistance Plasmids in Klebsiella pneumoniae

Background The development of multidrug resistance is a major problem in the treatment of pathogenic microorganisms by distinct antimicrobial agents. Characterizing the genetic variation among plasmids from different bacterial species or strains is a key step towards understanding the mechanism of virulence and their evolution. Results We applied a deep sequencing approach to 206 clinical strains of Klebsiella pneumoniae collected from 2002 to 2008 to understand the genetic variation of multidrug resistance plasmids, and to reveal the dynamic change of drug resistance over time. First, we sequenced three plasmids (70 Kb, 94 Kb, and 147 Kb) from a clonal strain of K. pneumoniae using Sanger sequencing. Using the Illumina sequencing technology, we obtained more than 17 million of short reads from two pooled plasmid samples. We mapped these short reads to the three reference plasmid sequences, and identified a large number of single nucleotide polymorphisms (SNPs) in these pooled plasmids. Many of these SNPs are present in drug-resistance genes. We also found that a significant fraction of short reads could not be mapped to the reference sequences, indicating a high degree of genetic variation among the collection of K. pneumoniae isolates. Moreover, we identified that plasmid conjugative transfer genes and antibiotic resistance genes are more likely to suffer from positive selection, as indicated by the elevated rates of nonsynonymous substitution. Conclusion These data represent the first large-scale study of genetic variation in multidrug resistance plasmids and provide insight into the mechanisms of plasmid diversification and the genetic basis of antibiotic resistance.

ampG Gene of Pseudomonas aeruginosa and Its Role in β-Lactamase Expression

ABSTRACT In enterobacteria, the ampG gene encodes a transmembrane protein (permease) that transports 1,6-GlcNAc-anhydro-MurNAc and the 1,6-GlcNAc-anhydro-MurNAc peptide from the periplasm to the cytoplasm, which serve as signal molecules for the induction of ampC β-lactamase. The role of AmpG as a transporter is also essential for cell wall recycling. Pseudomonas aeruginosa carries two AmpG homologues, AmpG (PA4393) and AmpGh1 (PA4218), with 45 and 41% amino acid sequence identity, respectively, to Escherichia coli AmpG, while the two homologues share only 19% amino acid identity. In P. aeruginosa strains PAO1 and PAK, inactivation of ampG drastically repressed the intrinsic β-lactam resistance while ampGh1 deletion had little effect on the resistance. Further, deletion of ampG in an ampD-null mutant abolished the high-level β-lactam resistance that is associated with the loss of AmpD activity. The cloned ampG gene is able to complement both the P. aeruginosa and the E. coli ampG mutants, while that of ampGh1 failed to do so, suggesting that PA4393 encodes the only functional AmpG protein in P. aeruginosa. We also demonstrate that the function of AmpG in laboratory strains of P. aeruginosa can effectively be inhibited by carbonyl cyanide m-chlorophenylhydrazone (CCCP), causing an increased sensitivity to β-lactams among laboratory as well as clinical isolates of P. aeruginosa. Our results suggest that inhibition of the AmpG activity is a potential strategy for enhancing the efficacy of β-lactams against P. aeruginosa, which carries inducible chromosomal ampC, especially in AmpC-hyperproducing clinical isolates.

Effects of curcumin on levels of nitric oxide synthase and AQP-4 in a rat model of hypoxia-ischemic brain damage.

This study examines the preventive and therapeutic effects of curcumin on brain edema after hypoxic-ischemic brain damage (HIBD) in a rat model. Male Sprague-Dawley rats were divided into four groups: a sham group (SH), a hypoxic-ischemic group (HI) without drug treatment, a hypoxic-ischemic group (CU) with curcumin injection, and a hypoxic-ischemic group with DMSO injection (solvent control, SC). HIBD treatment led to edema and ultrastructural changes in the hippocampus, increased the activity levels of nitric oxide synthase (NOS) in the brain (P<0.05), and raised the expression of water channel protein 4 (Aquaporin-4, AQP-4) in the blood-brain barrier (BBB) (P<0.05). Curcumin injection, but not the control DMSO injection, partially reversed HIBD-induced brain edema and morphological changes, as well as HIBD-induced increase in NOS activities and AQP-4 expression (P<0.05). In conclusion, our results showed that BBB ultrastructural changes may play an important role in the formation and development of brain edema after HIBD. Curcumin may protect the BBB ultrastructure and thus decrease brain edema following HIBD by down-regulating HIBD-induced increase in NOS activities and AQP-4 protein expression.

cTFbase: a database for comparative genomics of transcription factors in cyanobacteria

BackgroundComprehensive identification and classification of the transcription factors (TFs) in a given genome is an important aspect in understanding transcriptional regulatory networks of a specific organism. Cyanobacteria are an ancient group of gram-negative bacteria with strong variation in genome size ranging from about 1.6 to 9.1 Mb and little is known about their TF repertoires. Therefore, we constructed the cTFbase database to classify and analyze all the putative TFs in cyanobacterial genomes, followed by genome-wide comparative analysis.DescriptionIn the current release, cTFbase contains 1288 putative TFs identified from 21 fully sequenced cyanobacterial genomes. Through its user-friendly interactive interface, users can employ various criteria to retrieve all TF sequences and their detailed annotation information, including sequence features, domain architecture and sequence similarity against the linked databases. Furthermore, cTFbase provides phylogenetic trees of individual TF family, multiple sequence alignments of the DNA-binding domains and ortholog identification from any selected genomes. Comparative analysis revealed great variability of the TF sequences in cyanobacterial genomes. The high variance on the gene number and domain organization would be related to their diverse biological functions and their adaptation to various environmental conditions.ConclusioncTFbase provides a centralized warehouse for comparative analysis of putative TFs in cyanobacterial genomes. The availability of such an extensive database would be of great interest for the community of researchers working on TFs or transcriptional regulatory networks in cyanobacteria. cTFbase can be freely accessible at http://cegwz.com/ and will be continuously updated when the newly sequenced cyanobacterial genomes are available.

Sequence Analysis of pKF3-70 in Klebsiella pneumoniae: Probable Origin from R100-Like Plasmid of Escherichia coli

Background Klebsiella pneumoniae is a clinically significant species of bacterium which causes a variety of diseases. Clinical treatment of this bacterial infection is greatly hindered by the emergence of multidrug-resistant strains. The resistance is largely due to the acquisition of plasmids carrying drug-resistant as well as pathogenic genes, and its conjugal transfer facilitates the spread of resistant phenotypes. Methodology/Principal Findings The 70,057 bp plasmid pKF3-70, commonly found in Klebsiella pneumoniae, is composed of five main functional modules, including regions involved in replication, partition, conjugation, transfer leading, and variable regions. This plasmid is more similar to several Escherichia coli plasmids than any previously reported K. pneumoniae plasmids and pKF3-70 like plasmids share a common and conserved backbone sequence. The replication system of the pKF3-70 is 100% identical to that of RepFII plasmid R100 from E. coli. A beta-lactamase gene ctx-m-14 with its surrounding insertion elements (ISEcp1, truncated IS903 and a 20 bp inverted repeat sequence) may compose an active transposon which is directly bordered by two putative target repeats “ATTAC.” Conclusions/Significance The K. pneumoniae plasmid pKF3-70 carries an extended-spectrum beta-lactamase gene, ctx-m-14. The conjugative characteristic makes it a widespread plasmid among genetically relevant genera which poses significant threat to public health.

Bioinformatic analysis of the Acinetobacter baumannii phage AB1 genome

As one of the pathogens of hospital-acquired infections, Acinetobacter baumannii poses great challenges to the public health. A. baumannii phage could be an effective way to fight multi-resistant A. baumannii. Here, we completed the whole genome sequencing of the complete genome of A. baumannii phage AB1, which consists of 45,159 bp and is a double-stranded DNA molecule with an average GC content of 37.7%. The genome encodes one tRNA gene and 85 open reading frames (ORFs) and the average size of the ORF is 531 bp in length. Among 85 ORFs, only 14 have been identified to share significant sequence similarities to the genes with known functions, while 28 are similar in sequence to the genes with function-unknown genes in the database and 43 ORFs are uniquely present in the phage AB1 genome. Fourteen function-assigned genes with putative functions include five phage structure proteins, an RNA polymerase, a big sub-unit and a small sub-unit of a terminase, a methylase and a recombinase and the proteins involved in DNA replication and so on. Multiple sequence alignment was conducted among those homologous proteins and the phylogenetic trees were reconstructed to analyze the evolutionary courses of these essential genes. From comparative genomics analysis, it turned out clearly that the frame of the phage genome mainly consisted of genes from Xanthomonas phages, Burkholderia ambifaria phages and Enterobacteria phages and while it comprises genes of its host A. baumannii only sporadically. The mosaic feature of the phage genome suggested that the horizontal gene transfer occurred among the phage genomes and between the phages and the host bacterium genomes. Analyzing the genome sequences of the phages should lay sound foundation to investigate how phages adapt to the environment and infect their hosts, and even help to facilitate the development of biological agents to deal with pathogenic bacteria.

Identification of differentially expressed proteins of Arthrospira (Spirulina) plantensis-YZ under salt-stress conditions by proteomics and qRT-PCR analysis.

Arthrospira (Spirulina) platensis as a representative species of cyanobacteria has been recognized and used worldwide as a source of protein in the food, which possesses some unusual and valuable physiological characteristics, such as alkali and salt tolerance. Based on complete genome sequencing of Arthrospira (Spirulina) plantensis-YZ, we compared the protein expression profiles of this organism under different salt-stress conditions (i.e. 0.02 M, 0.5 M and 1.0 M NaCl, respectively), using 2-D electrophoresis and peptide mass fingerprinting, and retrieved 141 proteins showing significantly differential expression in response to salt-stress. Of the 141 proteins, 114 Arthrospira (Spirulina) plantensis-YZ proteins were found with significant homology to those found in Arthrospira (76 proteins in Arthrospira platensis str. Paraca and 38 in Arthrospira maxima CS-328). The remaining 27 proteins belong to other bacteria. Subsequently, we determined the transcriptional level of 29 genes in vivo in response to NaCl treatments and verified them by qRT-PCR. We found that 12 genes keep consistency at both transcription and protein levels, and transcription of all of them but one were up-regulated. We classified the 141 differentially expressed proteins into 18 types of function categories using COG database, and linked them to their respective KEGG metabolism pathways. These proteins are involved in 31 metabolism pathways, such as photosynthesis, glucose metabolism, cysteine and methionine metabolism, lysine synthesis, fatty acid metabolism, glutathione metabolism. Additionally, the SRPs, heat shock protein and ABC transporter proteins were identified, which probably render Arthrospira (Spirulina) plantensis’s resistance against high salt stress.

Codon usage patterns and adaptive evolution of marine unicellular cyanobacteria Synechococcus and Prochlorococcus

Marine unicellular cyanobacteria, represented by Synechococcus and Prochlorococcus, dominate the total phytoplankton biomass and production in oligotrophic ocean. In this study, we employed comparative genomics approaches to extensively investigate synonymous codon usage bias and evolutionary rates in a large number of closely related species of marine unicellular cyanobacteria. Although these two groups of marine cyanobacteria have a close phylogenetic relationship, we find that they are highly divergent not only in codon usage patterns but also in the driving forces behind the diversification. It is revealed that in Prochlorococcus, mutation and genome compositional constraints are the main forces contributing to codon usage bias, whereas in Synechococcus, translational selection. In addition, nucleotide substitution rate analysis indicates that they are not evolving at a constant rate after the divergence and that the average d(N)/d(S) values of core genes in Synechococcus are significantly higher than those in Prochlorococcus. Our evolutionary genomic analysis provides the first insight into codon usage, evolutionary genetic mechanisms and environmental adaptation of Synechococcus and Prochlorococcus after divergence.