Search

Background: Tracking the spread of the Neisseria gonorrhoeae strains with decreased susceptibility or resistance to cephalosporins is a major priority for global surveillance programmes. Whole-genome sequencing (WGS) has been widely used by increasing countries in North America, Europe, and Pacific to determine the decreased susceptible or resistance determinants of Neisseria gonorrhoeae, track the spread of these determinants throughout the gonococcal population at national or regional level. However, no studies to date have examined the genomic epidemiology of gonorrhea in Asia where the antimicrobial resistant strains of Neisseria gonorrhoeae appears to have emerged before disseminating the strains globally. Methods: We obtained clinical isolates and data from the China Gonococcal Resistance Surveillance Programme (China-GRSP) from 2012 to 2013. We sequenced the genomes of 435 clinical isolates of Neisseria gonorrhoeae, including 112 (25.6%) isolates with decreased susceptibility to ceftriaxone (Cfx-DS). We assessed the association between antimicrobial resistance genotype and phenotype. We also compared our data with the whole genome data of the isolates from the USA and the UK in the GenBank. Findings: The most prevalent MLST STs in our gonococcal population were MLST ST7827 (n = 74), followed by ST7365 (n = 58), ST1600 (n = 38), ST7367 (n = 35), and ST7363 (n = 29). MLST ST1901 which was reported as the predominant ST in the US was not found in our population. A total of 2512 strains, including additional 2077 published NG strains, were further included for phylogenetic analysis. It generated two distinct lineages - lineage 1 and lineage 2. Analysis of MLST ST1901 in the database indicate that most of MLST ST1901 isolates in the lineage2.6 were Cfx-DS isolates while all isolates in the lineage 2.1 were sensitive to ceftriaxone (77/110 vs. 0/13; p

Background Circular RNAs (circRNAs) are a group of noncoding RNAs that participate in gene expression regulation in various pathways. The essential roles of circRNAs have been revealed in many species. However, knowledge of circRNAs in fungi is still not comprehensive. Results Trichophyton rubrum (T. rubrum) is considered a model organism of human pathogenic filamentous fungi and dermatophytes. In this study, we performed a genome-wide investigation of circRNAs in T. rubrum based on high-throughput sequencing and ultimately identified 4254 circRNAs. Most of these circRNAs were specific to the conidial or mycelial stage, revealing a developmental stage-specific expression pattern. In addition, 940 circRNAs were significantly differentially expressed between the conidial and mycelial stages. PCR experiments conducted on seven randomly selected differentially expressed (DE-) circRNAs confirmed the circularized structures and relative expression levels of these circRNAs. Based on their genome locations, most circRNAs originated from intergenic regions, unlike those in plants and animals. Furthermore, we constructed circRNA-miRNA-mRNA regulatory networks that included 661 DE-circRNAs targeting 140 miRNAs and further regulating 2753 mRNAs. The relative expression levels of two randomly selected circRNA-miRNA-mRNA axes were investigated by qRT-PCR, and the competing endogenous RNA (ceRNA) network theory was validated. Functional enrichment analysis of the target genes suggested that they were significantly involved in posttranscriptional processes and protein synthesis as well as some small-molecule metabolism processes. CircRNAs are relatively more conserved in closely related dermatophytes but rarely conserved in distantly related species. Tru_circ07138_001 is a highly conserved circRNA that was conserved in all ten dermatophytes analyzed in our study and three distantly related species. Its host gene TERG_07138 was also highly conserved in two of these distantly related species Gallus gallus and Caenorhabditis elegans. The specific role of this circRNA deserves further exploration. Conclusions Our study is the first to provide a global profile of circRNAs in T. rubrum as well as dermatophytes. These results could serve as valuable resources for research on circRNA regulatory mechanisms in fungi and reveal new insights for further investigation of the physical characteristics of these significant human fungal pathogens.

Whole-genome sequencing was used to analyze the profiles of isoniazid (INH) resistance-related mutations among 188 multidrug-resistant strains of Mycobacterium tuberculosis (MDR-TB) and mono-INH-resistant isolates collected in a recent Chinese national survey. Mutations were detected in 18 structural genes and two promoter regions in 96.8% of 188 resistant isolates. There were high mutation frequencies in katG, the inhA promoter, and ahpC-oxyR regulator regions in INH-resistant isolates with frequencies of 86.2%, 19.6%, and 18.6%, respectively. Moreover, a high diversity of mutations was identified as 102 mutants contained various types of single or combined gene mutations in the INH-resistant group of isolates. The cumulative frequencies of katG 315 or inhA-P/inhA mutations was 68.1% (128/188) for the INH-resistant isolates. Of these isolates, 46 isolates (24.5% of 188) exhibited a high level of resistance. A high level of resistance was also observed in 21 isolates (11.2% of 188) with single ahpC-oxyR mutations or a combination of ahpC-oxyR and katG non-315 mutations. The remaining 17 mutations occurred sporadically and emerged in isolates with combined katG mutations. Such development of INH resistance is likely due to an accumulation of mutations under the pressure of drug selection. Thus, these findings provided insights on the levels of INH resistance and its correlation with the combinatorial mutation effect resulting from less frequent genes (inhA and/or ahpC). Such knowledge of other genes (apart from katG) in high-level resistance will aid in developing better strategies for the diagnosis and management of TB.

The global epidemic of drug-resistant tuberculosis is a catastrophic example of how antimicrobial resistance is undermining the public health gains made possible by combination drug therapy. Recent evidence points to unappreciated bacterial factors that accelerate the emergence of drug resistance. In a genome-wide association study of Mycobacterium tuberculosis isolates from China, we find mutations in the gene encoding the transcription factor prpR enriched in drug-resistant strains. prpR mutations confer conditional drug tolerance to three of the most effective classes of antibiotics by altering propionyl-CoA metabolism. prpR-mediated drug tolerance is carbon-source dependent, and while readily detectable during infection of human macrophages, is not captured by standard susceptibility testing. These data define a previously unrecognized and clinically prevalent class of M. tuberculosis variants that undermine antibiotic efficacy and drive drug resistance.

Background Recent outbreaks of severe pneumonia or acute respiratory distress syndrome have attracted much public interest. Rapid and accurate diagnosis of the causative agent is key for an adequate response to suspected outbreaks. Objectives We report a case that highlights the potential of semiconductor sequencing to rapidly determine the novel virus genome sequences. Study design We have developed a method for rapid de novo assembly of the novel influenza A H7N9 virus genome directly from the tracheal aspirate of a patient using semiconductor sequencer without culture and prior sequence information. Further, characteristic amino acids were analyzed and phylogenetic analysis were done for key genes of the influenza A virus. Results Deep sequencing yielded 435,239 reads assigned to H7N9 viruses, with an average length of 172 bp, accounting for 18.6% of total reads (2,339,680). Complete genome of the virus was obtained by de novo assembly method within 2 days. Genomic average depth of coverage of the Ion Torrent PGM was up to 5679 fold. Selected characteristic amino acids were observed, and phylogenetic analyses showed that the novel H7 virus was genetically close to 2011 duck H7N3 viruses in Zhejiang. The novel N9 sequences were most closely related to gene sequences of N9 derived from ducks H11N9 in 2011 in Jiangxi and H2N9 sequences from Hong Kong in 2010, in China, and therefore they may share a common ancestor. Conclusions The sequence-independent semiconductor sequencing is a powerful tool to investigate outbreak of a novel pathogen.

The Shigella bacteria cause bacillary dysentery, which remains a significant threat to public health. The genus status and species classification appear no longer valid, as compelling evidence indicates that Shigella, as well as enteroinvasive Escherichia coli, are derived from multiple origins of E.coli and form a single pathovar. Nevertheless, Shigella dysenteriae serotype 1 causes deadly epidemics but Shigella boydii is restricted to the Indian subcontinent, while Shigella flexneri and Shigella sonnei are prevalent in developing and developed countries respectively. To begin to explain these distinctive epidemiological and pathological features at the genome level, we have carried out comparative genomics on four representative strains. Each of the Shigella genomes includes a virulence plasmid that encodes conserved primary virulence determinants. The Shigella chromosomes share most of their genes with that of E.coli K12 strain MG1655, but each has over 200 pseudogenes, 300 approximately 700 copies of insertion sequence (IS) elements, and numerous deletions, insertions, translocations and inversions. There is extensive diversity of putative virulence genes, mostly acquired via bacteriophage-mediated lateral gene transfer. Hence, via convergent evolution involving gain and loss of functions, through bacteriophage-mediated gene acquisition, IS-mediated DNA rearrangements and formation of pseudogenes, the Shigella spp. became highly specific human pathogens with variable epidemiological and pathological features.

Rapid determination of full genomes of influenza A viruses from clinical samples is vital to outbreak preparedness. We benchmarked Ion Torrent PGM and Illumina MiSeq to the standard Sanger sequencing on two real clinical samples and demonstrated their performance for full-genome determination and polymorphism identification.

Serogroup A meningococcal strains have been involved in several pandemics and a series of epidemics worldwide in the past. Determination of the genome sequence of the prevalent genotype strain will help us understand the genetic background of the evolutionary and epidemiological properties of these bacteria. We sequenced the complete genome of Neisseria meningitidis NMA510612, a clinical isolate from a patient with meningococcal meningitis.

To the Editor: Middle East respiratory syndrome coronavirus (MERS-CoV), a novel lineage C betacoronavirus, was first described in September 2012, and by April 16, 2014, the virus had caused 238 infections and 92 deaths in humans worldwide (1). Antibodies against MERS-CoV in dromedary camels were recently reported (2), as was the full genome of MERS-CoV from dromedary camels (3). Finding the natural reservoir of MERS-CoV is fundamental to our ability to control transmission of this virus to humans (4). We report a novel lineage C betacoronavirus identified from Vespertilio superans bats in China. The full-length genome of this betacoronavirus showed close genetic relationship with MERS-CoV. Together with other evidence of MERS-CoV–related viruses in bats (5–8), our findings suggest that bats might be the natural reservoirs of MERS-related CoVs. In June 2013, we collected anal swab samples from 32 V. superans bats from southwestern China. A small proportion of each sample was pooled (without barcoding) and processed by using virus particle–protected nucleic acid purification and sequence-independent PCR for next-generation sequencing analysis with the Illumina (Solexa) Genome Analyzer II (Illumina, San Diego, CA, USA). Redundant reads were filtered, as described (9), from the raw sequencing reads generated by the genome analyzer and then aligned with the nonredundant protein database of the National Center for Biotechnology Information (ftp://ftp.ncbi.nlm.nih.gov/blast/db/) by using BLAST (http://blast.ncbi.nlm.nih.gov). The taxonomy of these aligned reads was parsed by using MEGAN 4 (http://ab.inf.uni-tuebingen.de/software/megan/). On the basis of the BLAST results, 8,751,354 sequence reads 81 nt in length were aligned with the protein sequences of the nonredundant protein database: 72,084 of the reads were uniquely matched with virus proteins. Of these 72,084 reads, 32,365 were assigned to the family Coronaviridae, primarily to lineage C of the genus Betacoronavirus, and found to share 60%–97% aa identity with MERS-CoV. The MERS-CoV–related reads were extracted and assembled by using SeqMan software from the Lasergene 7.1.0 program (DNASTAR, Madison, WI, USA), resulting in a draft CoV genome. Reverse transcription PCR selective for the partial RNA-dependent RNA polymerase (RdRp) gene of this novel lineage C betacoronavirus suggested that 5 of the 32 samples (≈16%) were positive for the novel betacoronavirus, and the PCR amplicons shared >98% nt identity with each other. Using a set of overlapped nested PCRs and the rapid amplification of cDNA ends method, we determined the full-length genome of 1 strain of this V. superans bat–derived betacoronavirus (referred to as BtVs-BetaCoV/SC2013, GenBank accession no. {"type":"entrez-nucleotide","attrs":{"text":"KJ473821","term_id":"627792518","term_text":"KJ473821"}}KJ473821). The betacoronavirus strain had a genome length of 30,413 nt, excluding the 3′ poly (A) tails, and a G+C content of 43.1%. Pairwise genome sequence alignment, conducted by the EMBOSS Needle software (http://www.ebi.ac.uk/Tools/psa/emboss_needle/) with default parameters, suggested that the genome sequence of BtVs-BetaCoV/SC2013 showed 75.7% nt identity with that of human MERS-CoV (hCoV-MERS); this shared identity is higher than that for other lineage C betacoronaviruses (from bats and hedgehogs) with full genomes available. hCoV-MERS showed 69.9% nt identity with bat CoV (BtCoV) HKU4-1, 70.1% nt identity with BtCoV-HKU5-1, and 69.6% nt identity with hedgehog CoV EriCoV-2012–174. Compared with those lineage C betacoronaviruses, which had an 816-bp partial RdRp sequence fragment available, BtVs-BetaCoV/SC2013 shared 96.7 % aa identity with hCoV-MERS. Pipistrellus BtCoVs found in Europe (BtCoV-8-724, BtCoV-8-691, BtCoV-UKR-G17) shared 98.2 % aa identity with hCoV-MERS, and Eptesicus BtCoV found in Italy (BtCoV-ITA26/384/2012) and other lineage C betacoronaviruses shared 96.3 % aa and

An epidemic caused by influenza A (H7N9) virus was recently reported in China. Deep sequencing revealed the full genome of the virus obtained directly from a patient’s sputum without virus culture. The full genome showed substantial sequence heterogeneity and large differences compared with that from embryonated chicken eggs.

Characterization of expression profile of organisms in response to antimicrobials provides important information on the potential mechanism of action of the drugs. The special expression signature can be used to predict whether other drugs act on the same target. Here, the common response of Shigella flexneri to two inhibitors of RNA polymerase was examined using gene expression profiling. Consistent with similar effects of the two drugs, the gene expression profiles indicated that responses of the bacteria to these drugs were roughly the same, with 225 genes affected commonly. Of them, 88 were induced and 137 were repressed. Real-time PCR was performed for selected genes to verify the microarray results. Analysis of the expression data revealed that more than 30% of the plasmid-encoded genes on the array were up-regulated by the antibiotics including virF regulon, other virulence-related genes, and genes responsible for plasmid replication, maintenance, and transfer. In addition, some chromosome-encoded genes involved in virulence and genes acquired from horizontal transfer were also significantly up-regulated. However, the expression of genes encoding the beta-subunit of RNA polymerase was increased moderately. The repressed genes include those that code for products associated with the ribosome, citrate cycle, glycolysis, thiamine biosynthesis, purine metabolism, fructose metabolism, mannose metabolism, and cold shock proteins. This study demonstrates that the two antibiotics induce rapid cessation of RNA synthesis resulting in inhibition of translation components. It also indicates that the production of virulence factors involved in intercellular dissemination, tissue invasion and inflammatory destruction may be enhanced through derepressing horizontal transfer genes by the drugs.

Bacillus cereus strain Q1 was isolated from a deep-subsurface oil reservoir in the Daqing oil field in northeastern China. This strain is able to produce biosurfactants and to survive in extreme environments. Here we report the finished and annotated genome sequence of this organism.

Summary Objective Information regarding the different types of FetA and PorB meningococci that circulate in various regions of the world is still scarce. The present study investigated the distribution of FetA and PorB variable region (VR) types among meningococci belonging to hyperinvasive lineages circulating in China. Methods The approach consisted of genotypic analysis of 201 Neisseria meningitidis strains belonging to hyperinvasive lineages isolated in China during the period 1956–2006. Results Sixteen different PorB types were found, 8 of which were newly identified. Of the 24 different FetA VR types, 3 were determined to be novel. Particular combinations of FetA and PorB types associated with distinct clonal complexes were also observed. Most cases of invasive disease were caused by five individual clones: A: P1.7-1,10: F5-5: ST-3 (cc1) with P3.6,11,10,7 (class 3 PorB protein; VR1-6, VR2-11, VR3-10, and VR4-7); A: P1.20,9: F3-1: ST-5 (cc5) with P3.4,11,10,7; A: P1.20,9: F3-1: ST-5 (cc5) with P3.9,11,10,7; A: P1.20,9: F3-1: ST-7 (cc5) with P3.4,11,10,7; and C: P1.7-2,14: F3-3: ST-4821 (cc4821) with P3.9,15,6,7. Conclusion A number of antigen–gene variants and combinations exhibited broad temporal and geographic distributions, although several invasive clones were mainly associated with a specified timeframe. The changes that are increasingly emerging in circulating strains and the prevalent clone replacement describe the molecular epidemiology of meningococcal disease in China. Our findings have implications for both public-health monitoring and further study of this organism.

Background Dermatophytes are the primary causative agent of dermatophytoses, a disease that affects billions of individuals worldwide. Trichophyton rubrum is the most common of the superficial fungi. Although T. rubrum is a recognized pathogen for humans, little is known about how its transcriptional pattern is related to development of the fungus and establishment of disease. It is therefore necessary to identify genes whose expression is relevant to growth, metabolism and virulence of T. rubrum. Results We generated 10 cDNA libraries covering nearly the entire growth phase and used them to isolate 11,085 unique expressed sequence tags (ESTs), including 3,816 contigs and 7,269 singletons. Comparisons with the GenBank non-redundant (NR) protein database revealed putative functions or matched homologs from other organisms for 7,764 (70%) of the ESTs. The remaining 3,321 (30%) of ESTs were only weakly similar or not similar to known sequences, suggesting that these ESTs represent novel genes. Conclusion The present data provide a comprehensive view of fungal physiological processes including metabolism, sexual and asexual growth cycles, signal transduction and pathogenic mechanisms.

We determined and analyzed the Shigella flexneri serotype 5 (pSF5) and S. dysenteriae serotype 1 (pSD1) virulence plasmid genomes. The total length of pSF5 is 136513 bp, including 165 open reading frames (ORFs). Of these ORFs, 133 were identified and 32 of those had no significant homology to proteins with known functions. The length of pSD1 is 182545 bp, including 224 ORFs, of which we identified 181. The remaining 43 ORFs were not significantly homologous to proteins with known functions. The insertion sequence (IS) elements are 53787 bp in pSF5, and 49616 bp in pSD1, which represents 39.4% and 27.1% of the genome, respectively. There are 22 IS element types in pSF5 and pSD1, among which we report ISEc8 and ISSbo6 for the first time in the Shigella virulence plasmid. Compared to pCP301, there are a large number of deleted genes and gene inversions in both pSF5 and pSD1. The ipa-mxi-spa locus in pSF5 is completely absent, and the genes related to the O-antigen biosynthesis are partially missing. In contrast, the above genes in pSD1 are integral, with the exception of virF. The whole genome analysis of the two plasmids shows that the loss of genes related to gene invasion or regulation also obliterates the ability of pPF5 and pSD1 to bind Congo red (Crb). Whether these genes determine the Crb function requires continued investigation.

The theory that Shigella is derived from multiple independent origins of Escherichia coli (Pupo et al. 2000) has been challenged by recent findings that the virulence plasmids (VPs) and the chromosomes share a similar evolutionary history (Escobar-Paramo et al. 2003), which suggests that an ancestral VP entered an E. coli strain only once, which gave rise to Shigella spp. In an attempt to resolve these conflicting theories, we constructed three phylogenetic trees in this study: a robust chromosomal tree using 23 housekeeping genes from 46 strains of Shigella and enteroinvasive E. coli (EIEC), a chromosomal tree using 4 housekeeping genes from 19 EcoR strains and 46 Shigella/EIEC strains, and a VP tree using 5 genes outside of the VP cell-entry region from 38 Shigella/EIEC strains. Both chromosomal trees group Shigella into three main clusters and five outliers, and strongly suggest that Shigella has multiple origins within E. coli. Most strikingly, the VP tree shows that the VPs from two main Shigella clusters, C1 and C2, are more closely related, which contradicts the chromosomal trees that place C2 and C3 next to each other but C1 at a distance. Additionally, we have identified a complete tra operon of the F-plasmid in the genome sequence of an EIEC strain and found that two other EIEC strains are also likely to possess a complete tra operon. All lines of evidence support an alternative multiorigin theory that transferable diverse ancestral VPs entered diverse origins of E. coli multiple times during a prolonged period of time, resulting in Shigella species with diverse genomes but similar pathogenic properties.

The genetic and phylogenetic characteristics of human enterovirus 71 (EV71) and coxsackievirus A16 (CA16) sampled from children with hand, foot, and mouth disease in Shenzhen, People's Republic of China, over a 6-year period (1999 to 2004) were examined with reverse transcription-PCR and DNA sequencing. Out of 147 stool specimens, 60 showed positive signals when screened with EV71- and CA16-specific primers. EV71 was identified in 19 specimens, and CA16 was identified in 41 specimens; coinfection by EV71 and CA16 was not observed. Phylogenetic analysis of all EV71 strains isolated from the mainland Chinese samples established C4 as the predominant genotype. Only one other known strain (3254-TAI-98; AF286531), isolated in Taiwan in 1998, was identified as belonging to genotype C4. Phylogenetic analysis of CA16 strains allowed us to identify three new genetic lineages (A, B, and C), with lineage C recently predominating in Asian countries, such as the People's Republic of China, Malaysia, and Japan. These new observations indicate that CA16 circulating in the People's Republic of China is genetically diverse, and additional surveillance is warranted.

The complete sequence of pSS, which is the large virulence plasmid of Shigella sonnei, was determined. The 214-kb plasmid is composed of segments of virulence-associated genes, the O-antigen gene clusters, a range of replication and maintenance genes, and large numbers of insertion sequence (IS) elements. Two hundred and forty-one open reading frames (ORFs) were identified, of which 117 are highly homologous to IS elements or transposases, 57 are homologous to known pathogenesis-associated proteins, and 30 are related to replication, plasmid maintenance, or other metabolic functions. Thirty-seven ORFs have no similarity to proteins with a known function, including two with no significant similarity to any hypothetical proteins. Interestingly, 10 ORFs encoding O-antigen gene clusters were identified on the plasmid and this is markedly different from most other Shigella spp. virulent plasmids. A novel toxin-antitoxin system, a series of stbDE homologs, was found on the plasmid immediately downstream of the replication region; the sole segregation stability system may be responsible for the instability of pSS. The pSS plasmid is a mixture of genes with different origins and functions. The sequence suggests a remarkable history of IS-mediated recombination and acquisition of DNA across a range of bacterial species.

Background Conidia are considered to be the primary cause of infections by Trichophyton rubrum. Results We have developed a cDNA microarray containing 10250 ESTs to monitor the transcriptional strategy of conidial germination. A total of 1561 genes that had their expression levels specially altered in the process were obtained and hierarchically clustered with respect to their expression profiles. By functional analysis, we provided a global view of an important biological system related to conidial germination, including characterization of the pattern of gene expression at sequential developmental phases, and changes of gene expression profiles corresponding to morphological transitions. We matched the EST sequences to GO terms in the Saccharomyces Genome Database (SGD). A number of homologues of Saccharomyces cerevisiae genes related to signalling pathways and some important cellular processes were found to be involved in T. rubrum germination. These genes and signalling pathways may play roles in distinct steps, such as activating conidial germination, maintenance of isotropic growth, establishment of cell polarity and morphological transitions. Conclusion Our results may provide insights into molecular mechanisms of conidial germination at the cell level, and may enhance our understanding of regulation of gene expression related to the morphological construction of T. rubrum.