Menu
July 7, 2019  |  

Comparative genomics of Burkholderia multivorans, a ubiquitous pathogen with a highly conserved genomic structure.

The natural environment serves as a reservoir of opportunistic pathogens. A well-established method for studying the epidemiology of such opportunists is multilocus sequence typing, which in many cases has defined strains predisposed to causing infection. Burkholderia multivorans is an important pathogen in people with cystic fibrosis (CF) and its epidemiology suggests that strains are acquired from non-human sources such as the natural environment. This raises the central question of whether the isolation source (CF or environment) or the multilocus sequence type (ST) of B. multivorans better predicts their genomic content and functionality. We identified four pairs of B. multivorans isolates, representing distinct STs and consisting of one CF and one environmental isolate each. All genomes were sequenced using the PacBio SMRT sequencing technology, which resulted in eight high-quality B. multivorans genome assemblies. The present study demonstrated that the genomic structure of the examined B. multivorans STs is highly conserved and that the B. multivorans genomic lineages are defined by their ST. Orthologous protein families were not uniformly distributed among chromosomes, with core orthologs being enriched on the primary chromosome and ST-specific orthologs being enriched on the second and third chromosome. The ST-specific orthologs were enriched in genes involved in defense mechanisms and secondary metabolism, corroborating the strain-specificity of these virulence characteristics. Finally, the same B. multivorans genomic lineages occur in both CF and environmental samples and on different continents, demonstrating their ubiquity and evolutionary persistence.


July 7, 2019  |  

Complete genome sequencing and targeted mutagenesis reveal virulence contributions of Tal2 and Tal4b of Xanthomonas translucens pv. undulosa ICMP11055 in bacterial leaf streak of wheat

Bacterial leaf streak caused by Xanthomonas translucens pv. undulosa (Xtu) is an important disease of wheat (Triticum aestivum) and barley (Hordeum vulgare) worldwide. Transcription activator-like effectors (TALEs) play determinative roles in many of the plant diseases caused by the different species and pathovars of Xanthomonas, but their role in this disease has not been characterized. ICMP11055 is a highly virulent Xtu strain from Iran. The aim of this study was to better understand genetic diversity of Xtu and to assess the role of TALEs in bacterial leaf streak of wheat by comparing the genome of this strain to the recently completely sequenced genome of a U.S. Xtu strain, and to several other draft X. translucens genomes, and by carrying out mutational analyses of the TALE (tal) genes the Iranian strain might harbor. The ICMP11055 genome, including its repeat-rich tal genes, was completely sequenced using single molecule, real-time technology (Pacific Biosciences). It consists of a single circular chromosome of 4,561,583 bp, containing 3,953 genes. Whole genome alignment with the genome of the United States Xtu strain XT4699 showed two major re-arrangements, nine genomic regions unique to ICMP11055, and one region unique to XT4699. ICMP110055 harbors 26 non-TALE type III effector genes and seven tal genes, compared to 25 and eight for XT4699. The tal genes occur singly or in pairs across five scattered loci. Four are identical to tal genes in XT4699. In addition to common repeat-variable diresidues (RVDs), the tal genes of ICMP11055, like those of XT4699, encode several RVDs rarely observed in Xanthomonas, including KG, NF, Y*, YD, and YK. Insertion and deletion mutagenesis of ICMP11055 tal genes followed by genetic complementation analysis in wheat cv. Chinese Spring revealed that Tal2 and Tal4b of ICMP11055 each contribute individually to the extent of disease caused by this strain. A largely conserved ortholog of tal2 is present in XT4699, but for tal4b, only a gene with partial, fragmented RVD sequence similarity can be found. Our results lay the foundation for identification of important host genes activated by Xtu TALEs as targets for the development of disease resistant varieties.


July 7, 2019  |  

Emergence and genomic diversification of a virulent serogroup W: ST-2881 (CC175) Neisseria meningitidis clone in the African meningitis belt

Countries of the African ‘meningitis belt’ are susceptible to meningococcal meningitis outbreaks. While in the past major epidemics have been primarily caused by serogroup A meningococci, W strains are currently responsible for most of the cases. After an epidemic in Mecca in 2000, W:ST-11 strains have caused many outbreaks worldwide. An unrelated W:ST-2881 clone was described for the first time in 2002, with the first meningitis cases caused by these bacteria reported in 2003. Here we describe results of a comparative whole-genome analysis of 74 W:ST-2881 strains isolated within the framework of two longitudinal colonization and disease studies conducted in Ghana and Burkina Faso. Genomic data indicate that the W:ST-2881 clone has emerged from Y:ST-175(CC175) bacteria by capsule switching. The circulating W:ST-2881 populations were composed of a variety of closely related but distinct genomic variants with no systematic differences between colonization and disease isolates. Two distinct and geographically clustered phylogenetic clonal variants were identified in Burkina Faso and a third in Ghana. On the basis of the presence or absence of 17 recombination fragments, the Ghanaian variant could be differentiated into five clusters. All 25 Ghanaian disease isolates clustered together with 23 out of 40 Ghanaian isolates associated with carriage within one cluster, indicating that W:ST-2881 clusters differ in virulence. More than half of the genes affected by horizontal gene transfer encoded proteins of the ‘cell envelope’ and the ‘transport/binding protein’ categories, which indicates that exchange of non-capsular antigens plays an important role in immune evasion.


July 7, 2019  |  

Insights into Cedecea neteri strain M006 through complete genome sequence, a rare bacterium from aquatic environment.

Cedecea neteri M006 is a rare bacterium typically found as an environmental isolate from the tropical rainforest Sungai Tua waterfall (Gombak, Selangor, Malaysia). It is a Gram-reaction-negative, facultative anaerobic, bacillus. Here, we explore the features of Cedecea neteri M006, together with its genome sequence and annotation. The genome comprised 4,965,436 bp with 4447 protein-coding genes and 103 RNA genes.


July 7, 2019  |  

Paenibacillus ihbetae sp. nov., a cold-adapted antimicrobial producing bacterium isolated from high altitude Suraj Tal Lake in the Indian trans-Himalayas.

The assessment of bacterial diversity and bioprospection of the high-altitude lake Suraj Tal microorganisms for potent antimicrobial activities revealed the presence of two Gram-stain-variable, endospore-forming, rod-shaped, aerobic bacteria, namely IHBB 9852(T) and IHBB 9951. Phylogenetic analysis based on 16S rRNA gene sequence showed the affiliation of strains IHBB 9852(T) and IHBB 9951 within the genus Paenibacillus, exhibiting the highest sequence similarity to Paenibacillus lactis DSM 15596(T) (97.8% and 97.7%) and less than 95.9% similarity to other species of the genus Paenibacillus. DNA-DNA relatedness among strains IHBB 9852(T) and IHBB 9951 was 90.2%, and with P. lactis DSM 15596(T), was 52.7% and 52.4%, respectively. The novel strains contain anteiso-C15:0, iso-C15:0, C16:0 and iso-C16:0 as major fatty acids, and phosphatidylglycerol, phosphatidylethanolamine and diphosphatidylglycerol were predominant polar lipids. The DNA G+C content for IHBB 9852T and IHBB 9951 was 52.1 and 52.2mol%. Based on the results of phenotypic and genomic characterisations, we concluded that strains IHBB 9852(T) and IHBB 9951 belong to a novel Paenibacillus species, for which the name Paenibacillus ihbetae sp. nov. is proposed. The type strain is IHBB 9852(T) (=MTCC 12459(T)=MCC 2795(T)=JCM 31131(T)=KACC 19072(T); DPD TaxonNumber TA00046) and IHBB 9951 (=MTCC 12458=MCC 2794=JCM 31132=KACC 19073) is a reference strain. Copyright © 2017. Published by Elsevier GmbH.


July 7, 2019  |  

Parallel evolution of two clades of a major Atlantic endemic Vibrio parahaemolyticus pathogen lineage by independent acquisition of related pathogenicity islands.

Shellfish-transmitted Vibrio parahaemolyticus infections have recently increased from locations with historically low disease incidence, such as the Northeast United States (US). This change coincided with a bacterial population shift towards human pathogenic variants occurring in part through the introduction of several Pacific native lineages (ST36, ST43 and ST636) to near-shore areas off the Atlantic coast of the Northeast US. Concomitantly, ST631 emerged as a major endemic pathogen. Phylogenetic trees of clinical and environmental isolates indicated that two clades diverged from a common ST631 ancestor, and in each of these clades, a human pathogenic variant evolved independently through acquisition of distinct Vibrio pathogenicity islands (VPaI). These VPaI differ from each other and bear little resemblance to hemolysin-containing VPaI from isolates of the pandemic clonal complex. Clade I ST631 isolates either harbored no hemolysins, or contained a chromosome I-inserted island we call VPaIß that encodes a type three secretion system (T3SS2ß) typical of Trh hemolysin-producers. The more clinically prevalent and clonal ST631 clade II had an island we call VPaI? that encodes both tdh and trh and that was inserted in chromosome II. VPaI? was derived from VPaIß but with some additional acquired elements in common with VPaI carried by pandemic isolates, exemplifying the mosaic nature of pathogenicity islands. Genomics comparisons and amplicon assays identified VPaI?-type islands containing tdh inserted adjacent to the ure cluster in the three introduced Pacific and most other emergent lineages. that collectively cause 67% of Northeast US infections as of 2016.IMPORTANCE The availability of three different hemolysin genotypes in the ST631 lineage provided a unique opportunity to employ genome comparisons to further our understanding of the processes underlying pathogen evolution. The fact that two different pathogenic clades arose in parallel from the same potentially benign lineage by independent VPaI acquisition is surprising considering the historically low prevalence of community members harboring VPaI in waters along the Northeast US coast that could serve as the source of this material. This illustrates a possible predisposition of some lineages to not only acquire foreign DNA but also to become human pathogens. Whereas the underlying cause for the expansion of V. parahaemolyticus lineages harboring VPaI? along the US Atlantic coast and spread of this element to multiple lineages that underlies disease emergence is not known, this work underscores the need to define the environment factors that favor bacteria harboring VPaI in locations of emergent disease. Copyright © 2017 American Society for Microbiology.


July 7, 2019  |  

The Tartary buckwheat genome provides insights into rutin biosynthesis and abiotic stress tolerance.

Tartary buckwheat (Fagopyrum tataricum) is an important pseudocereal crop that is strongly adapted to growth in adverse environments. Its gluten-free grain contains complete proteins with a well-balanced composition of essential amino acids and is a rich source of beneficial phytochemicals that provide significant health benefits. Here, we report a high-quality, chromosome-scale Tartary buckwheat genome sequence of 489.3 Mb that is assembled by combining whole-genome shotgun sequencing of both Illumina short reads and single-molecule real-time long reads, sequence tags of a large DNA insert fosmid library, Hi-C sequencing data, and BioNano genome maps. We annotated 33 366 high-confidence protein-coding genes based on expression evidence. Comparisons of the intra-genome with the sugar beet genome revealed an independent whole-genome duplication that occurred in the buckwheat lineage after they diverged from the common ancestor, which was not shared with rosids or asterids. The reference genome facilitated the identification of many new genes predicted to be involved in rutin biosynthesis and regulation, aluminum stress resistance, and in drought and cold stress responses. Our data suggest that Tartary buckwheat’s ability to tolerate high levels of abiotic stress is attributed to the expansion of several gene families involved in signal transduction, gene regulation, and membrane transport. The availability of these genomic resources will facilitate the discovery of agronomically and nutritionally important genes and genetic improvement of Tartary buckwheat. Copyright © 2017 The Author. Published by Elsevier Inc. All rights reserved.


July 7, 2019  |  

Genomic comparison between Staphylococcus aureus GN strains clinically isolated from a familial infection case: IS1272 transposition through a novel inverted repeat-replacing mechanism.

A bacterial insertion sequence (IS) is a mobile DNA sequence carrying only the transposase gene (tnp) that acts as a mutator to disrupt genes, alter gene expressions, and cause genomic rearrangements. “Canonical” ISs have historically been characterized by their terminal inverted repeats (IRs), which may form a stem-loop structure, and duplications of a short (non-IR) target sequence at both ends, called target site duplications (TSDs). The IS distributions and virulence potentials of Staphylococcus aureus genomes in familial infection cases are unclear. Here, we determined the complete circular genome sequences of familial strains from a Panton-Valentine leukocidin (PVL)-positive ST50/agr4 S. aureus (GN) infection of a 4-year old boy with skin abscesses. The genomes of the patient strain (GN1) and parent strain (GN3) were rich for “canonical” IS1272 with terminal IRs, both having 13 commonly-existing copies (ce-IS1272). Moreover, GN1 had a newly-inserted IS1272 (ni-IS1272) on the PVL-converting prophage, while GN3 had two copies of ni-IS1272 within the DNA helicase gene and near rot. The GN3 genome also had a small deletion. The targets of ni-IS1272 transposition were IR structures, in contrast with previous “canonical” ISs. There were no TSDs. Based on a database search, the targets for ce-IS1272 were IRs or “non-IRs”. IS1272 included a larger structure with tandem duplications of the left (IRL) side sequence; tnp included minor cases of a long fusion form and truncated form. One ce-IS1272 was associated with the segments responsible for immune evasion and drug resistance. Regarding virulence, GN1 expressed cytolytic peptides (phenol-soluble modulin a and d-hemolysin) and PVL more strongly than some other familial strains. These results suggest that IS1272 transposes through an IR-replacing mechanism, with an irreversible process unlike that of “canonical” transpositions, resulting in genomic variations, and that, among the familial strains, the patient strain has strong virulence potential based on community-associated virulence factors.


July 7, 2019  |  

DNA methylation profiling using long-read Single Molecule Real-Time bisulfite sequencing (SMRT-BS).

For the past two decades, bisulfite sequencing has been a widely used method for quantitative CpG methylation detection of genomic DNA. Coupled with PCR amplicon cloning, bisulfite Sanger sequencing allows for allele-specific CpG methylation assessment; however, its time-consuming protocol and inability to multiplex has recently been overcome by next-generation bisulfite sequencing techniques. Although high-throughput sequencing platforms have enabled greater accuracy in CpG methylation quantitation as a result of increased bisulfite sequencing depth, most common sequencing platforms generate reads that are similar in length to the typical bisulfite PCR size range (~300-500 bp). Using the Pacific Biosciences (PacBio) sequencing platform, we developed single molecule real-time bisulfite sequencing (SMRT-BS), which is an accurate targeted CpG methylation analysis method capable of a high degree of multiplexing and long read lengths. SMRT-BS is reproducible and was found to be concordant with other lower throughput quantitative CpG methylation methods. Moreover, the ability to sequence up to ~1.5-2.0 kb amplicons, when coupled with an optimized bisulfite-conversion protocol, allows for more thorough assessment of CpG islands and increases the capacity for studying the relationship between single nucleotide variants and allele-specific CpG methylation.


July 7, 2019  |  

HISEA: HIerarchical SEed Aligner for PacBio data.

The next generation sequencing (NGS) techniques have been around for over a decade. Many of their fundamental applications rely on the ability to compute good genome assemblies. As the technology evolves, the assembly algorithms and tools have to continuously adjust and improve. The currently dominant technology of Illumina produces reads that are too short to bridge many repeats, setting limits on what can be successfully assembled. The emerging SMRT (Single Molecule, Real-Time) sequencing technique from Pacific Biosciences produces uniform coverage and long reads of length up to sixty thousand base pairs, enabling significantly better genome assemblies. However, SMRT reads are much more expensive and have a much higher error rate than Illumina’s – around 10-15% – mostly due to indels. New algorithms are very much needed to take advantage of the long reads while mitigating the effect of high error rate and lowering the required coverage.An essential step in assembling SMRT data is the detection of alignments, or overlaps, between reads. High error rate and very long reads make this a much more challenging problem than for Illumina data. We present a new pairwise read aligner, or overlapper, HISEA (Hierarchical SEed Aligner) for SMRT sequencing data. HISEA uses a novel two-step k-mer search, employing consistent clustering, k-mer filtering, and read alignment extension.We compare HISEA against several state-of-the-art programs – BLASR, DALIGNER, GraphMap, MHAP, and Minimap – on real datasets from five organisms. We compare their sensitivity, precision, specificity, F1-score, as well as time and memory usage. We also introduce a new, more precise, evaluation method. Finally, we compare the two leading programs, MHAP and HISEA, for their genome assembly performance in the Canu pipeline.Our algorithm has the best alignment detection sensitivity among all programs for SMRT data, significantly higher than the current best. The currently best assembler for SMRT data is the Canu program which uses the MHAP aligner in its pipeline. We have incorporated our new HISEA aligner in the Canu pipeline and benchmarked it against the best pipeline for multiple datasets at two relevant coverage levels: 30x and 50x. Our assemblies are better than those using MHAP for both coverage levels. Moreover, Canu+HISEA assemblies for 30x coverage are comparable with Canu+MHAP assemblies for 50x coverage, while being faster and cheaper.The HISEA algorithm produces alignments with highest sensitivity compared with the current state-of-the-art algorithms. Integrated in the Canu pipeline, currently the best for assembling PacBio data, it produces better assemblies than Canu+MHAP.


July 7, 2019  |  

Characterization of Fusobacterium varium Fv113-g1 isolated from a patient with ulcerative colitis based on complete genome sequence and transcriptome analysis.

Fusobacterium spp. present in the oral and gut flora is carcinogenic and is associated with the risk of pancreatic and colorectal cancers. Fusobacterium spp. is also implicated in a broad spectrum of human pathologies, including Crohn’s disease and ulcerative colitis (UC). Here we report the complete genome sequence of Fusobacterium varium Fv113-g1 (genome size, 3.96 Mb) isolated from a patient with UC. Comparative genome analyses totally suggested that Fv113-g1 is basically assigned as F. varium, in particular, it could be reclassified as notable F. varium subsp. similar to F. ulcerans because of partial shared orthologs. Compared with the genome sequences of F. varium ATCC 27725 (genome size, 3.30 Mb) and other strains of Fusobacterium spp., Fv113-g1 possesses many accessary pan-genome sequences with noteworthy multiple virulence factors, including 44 autotransporters (type V secretion system, T5SS) and 13 Fusobacterium adhesion (FadA) paralogs involved in potential mucosal inflammation. Indeed, transcriptome analysis demonstrated that Fv113-g1-specific accessary genes, such as multiple T5SS and fadA paralogs, showed notably increased expression with D-MEM cultivation than with brain heart infusion broth. This implied that growth condition may enhance the expression of such potential virulence factors, leading to remarkable survival against other gut microorganisms and to the pathogenicity to human intestinal epithelium.


July 7, 2019  |  

A 3-way hybrid approach to generate a new high-quality chimpanzee reference genome (Pan_tro_3.0).

The chimpanzee is arguably the most important species for the study of human origins. A key resource for these studies is a high-quality reference genome assembly; however, as with most mammalian genomes, the current iteration of the chimpanzee reference genome assembly is highly fragmented. In the current iteration of the chimpanzee reference genome assembly (Pan_tro_2.1.4), the sequence is scattered across more then 183 000 contigs, incorporating more than 159 000 gaps, with a genome-wide contig N50 of 51 Kbp. In this work, we produce an extensive and diverse array of sequencing datasets to rapidly assemble a new chimpanzee reference that surpasses previous iterations in bases represented and organized in large scaffolds. To this end, we show substantial improvements over the current release of the chimpanzee genome (Pan_tro_2.1.4) by several metrics, such as increased contiguity by >750% and 300% on contigs and scaffolds, respectively, and closure of 77% of gaps in the Pan_tro_2.1.4 assembly gaps spanning >850 Kbp of the novel coding sequence based on RNASeq data. We further report more than 2700 genes that had putatively erroneous frame-shift predictions to human in Pan_tro_2.1.4 and show a substantial increase in the annotation of repetitive elements. We apply a simple 3-way hybrid approach to considerably improve the reference genome assembly for the chimpanzee, providing a valuable resource for the study of human origins. Furthermore, we produce extensive sequencing datasets that are all derived from the same cell line, generating a broad non-human benchmark dataset.© The Author 2017. Published by Oxford University Press.


July 7, 2019  |  

The plastid genome in Cladophorales green algae is encoded by hairpin chromosomes.

Virtually all plastid (chloroplast) genomes are circular double-stranded DNA molecules, typically between 100 and 200 kb in size and encoding circa 80-250 genes. Exceptions to this universal plastid genome architecture are very few and include the dinoflagellates, where genes are located on DNA minicircles. Here we report on the highly deviant chloroplast genome of Cladophorales green algae, which is entirely fragmented into hairpin chromosomes. Short- and long-read high-throughput sequencing of DNA and RNA demonstrated that the chloroplast genes of Boodlea composita are encoded on 1- to 7-kb DNA contigs with an exceptionally high GC content, each containing a long inverted repeat with one or two protein-coding genes and conserved non-coding regions putatively involved in replication and/or expression. We propose that these contigs correspond to linear single-stranded DNA molecules that fold onto themselves to form hairpin chromosomes. The Boodlea chloroplast genes are highly divergent from their corresponding orthologs, and display an alternative genetic code. The origin of this highly deviant chloroplast genome most likely occurred before the emergence of the Cladophorales, and coincided with an elevated transfer of chloroplast genes to the nucleus. A chloroplast genome that is composed only of linear DNA molecules is unprecedented among eukaryotes, and highlights unexpected variation in plastid genome architecture. Copyright © 2017 Elsevier Ltd. All rights reserved.


July 7, 2019  |  

Complete Sequences and Characterization of Two Novel Plasmids Carrying aac(6′)-Ib-cr and qnrS Gene in Shigella flexneri.

The complete sequences of two previously reported plasmids carrying plasmid-mediated quinolone resistance genes from Shigella flexneri in China have not been available. The present study using the p5-C3 assembly method revealed that (1) the plasmid pSF07201 with aac(6′)-Ib-cr had 75,335?bp with antibiotic resistance genes CTX-M-3, TEM-1, and FosA3; (2) seven fragments of pSF07201 had more than 99% homology with the seven corresponding plasmids; (3) the other plasmid pSF07202 with qnrS had 47,669?bp with antibiotic resistance gene TEM-1 and 99.95% homology with a segment of pKF362122, which has the qnrS gene from location 162,490 to 163,146. A conjugation and electrotransformation experiment suggested that these two plasmids might horizontally transfer between and coexist in Escherichia coli J53 and S. flexneri 2a 301. Either the aac(6′)-Ib-cr or qnrS gene contributed to, but only the coexistence of the two genes conferred to the resistance to ciprofloxacin in these two strains. To the best of our knowledge, this is the first report of the complete sequences of the aac(6′)-Ib-cr- and qnrS-positive plasmids in Shigella isolates. Our findings indicate that two genes probably evolve through horizontal plasmid transfer between the different bacterial types.


Talk with an expert

If you have a question, need to check the status of an order, or are interested in purchasing an instrument, we're here to help.