Menu
September 22, 2019

Occurrence, evolution, and functions of DNA phosphorothioate epigenetics in bacteria.

The chemical diversity of physiological DNA modifications has expanded with the identification of phosphorothioate (PT) modification in which the nonbridging oxygen in the sugar-phosphate backbone of DNA is replaced by sulfur. Together with DndFGH as cognate restriction enzymes, DNA PT modification, which is catalyzed by the DndABCDE proteins, functions as a bacterial restriction-modification (R-M) system that protects cells against invading foreign DNA. However, the occurrence of dnd systems across a large number of bacterial genomes and their functions other than R-M are poorly understood. Here, a genomic survey revealed the prevalence of bacterial dnd systems: 1,349 bacterial dnd systems were observed to occur sporadically across diverse phylogenetic groups, and nearly half of these occur in the form of a solitary dndBCDE gene cluster that lacks the dndFGH restriction counterparts. A phylogenetic analysis of 734 complete PT R-M pairs revealed the coevolution of M and R components, despite the observation that several PT R-M pairs appeared to be assembled from M and R parts acquired from distantly related organisms. Concurrent epigenomic analysis, transcriptome analysis, and metabolome characterization showed that a solitary PT modification contributed to the overall cellular redox state, the loss of which perturbed the cellular redox balance and induced Pseudomonas fluorescens to reconfigure its metabolism to fend off oxidative stress. An in vitro transcriptional assay revealed altered transcriptional efficiency in the presence of PT DNA modification, implicating its function in epigenetic regulation. These data suggest the versatility of PT in addition to its involvement in R-M protection.


September 22, 2019

Long-read genome sequence and assembly of Leptopilina boulardi: a specialist Drosophila parasitoid

Background: Leptopilina boulardi is a specialist parasitoid belonging to the order Hymenoptera, which attacks the larval stages of Drosophila. The Leptopilina genus has enormous value in the biological control of pests as well as in understanding several aspects of host-parasitoid biology. However, none of the members of Figitidae family has their genomes sequenced. In order to improve the understanding of the parasitoid wasps by generating genomic resources, we sequenced the whole genome of L. boulardi. Findings: Here, we report a high quality genome of L. boulardi, assembled from 70Gb of Illumina reads and 10.5Gb of PacBio reads, forming a total coverage of 230X. The 375Mb draft genome has an N50 of 275Kb with 6315 scaffolds >500bp, and encompasses >95% complete BUSCOs. The GC% of the genome is 28.26%, and RepeatMasker identified 868105 repeat elements covering 43.9% of the assembly. A total of 25259 protein-coding genes were predicted using a combination of ab-initio and RNA-Seq based methods, with an average gene size of 3.9Kb. 78.11% of the predicted genes could be annotated with at least one function. Conclusion: Our study provides a highly reliable assembly of this parasitoid wasp, which will be a valuable resource to researchers studying parasitoids. In particular, it can help delineate the host-parasitoid mechanisms that are part of the Drosophila-Leptopilina model system.


September 22, 2019

Stress-adaptive responses associated with high-level carbapenem resistance in KPC-producing Klebsiella pneumoniae.

Carbapenem-resistant Enterobacteriaceae (CRE) organisms have emerged to become a major global public health threat among antimicrobial resistant bacterial human pathogens. Little is known about how CREs emerge. One characteristic phenotype of CREs is heteroresistance, which is clinically associated with treatment failure in patients given a carbapenem. Through in vitro whole-transcriptome analysis we tracked gene expression over time in two different strains (BR7, BR21) of heteroresistant KPC-producing Klebsiella pneumoniae, first exposed to a bactericidal concentration of imipenem followed by growth in drug-free medium. In both strains, the immediate response was dominated by a shift in expression of genes involved in glycolysis toward those involved in catabolic pathways. This response was followed by global dampening of transcriptional changes involving protein translation, folding and transport, and decreased expression of genes encoding critical junctures of lipopolysaccharide biosynthesis. The emerged high-level carbapenem-resistant BR21 subpopulation had a prophage (IS1) disrupting ompK36 associated with irreversible OmpK36 porin loss. On the other hand, OmpK36 loss in BR7 was reversible. The acquisition of high-level carbapenem resistance by the two heteroresistant strains was associated with distinct and shared stepwise transcriptional programs. Carbapenem heteroresistance may emerge from the most adaptive subpopulation among a population of cells undergoing a complex set of stress-adaptive responses.


September 22, 2019

The global distribution and spread of the mobilized colistin resistance gene mcr-1.

Colistin represents one of the few available drugs for treating infections caused by carbapenem-resistant Enterobacteriaceae. As such, the recent plasmid-mediated spread of the colistin resistance gene mcr-1 poses a significant public health threat, requiring global monitoring and surveillance. Here, we characterize the global distribution of mcr-1 using a data set of 457 mcr-1-positive sequenced isolates. We find mcr-1 in various plasmid types but identify an immediate background common to all mcr-1 sequences. Our analyses establish that all mcr-1 elements in circulation descend from the same initial mobilization of mcr-1 by an ISApl1 transposon in the mid 2000s (2002-2008; 95% highest posterior density), followed by a marked demographic expansion, which led to its current global distribution. Our results provide the first systematic phylogenetic analysis of the origin and spread of mcr-1, and emphasize the importance of understanding the movement of antibiotic resistance genes across multiple levels of genomic organization.


September 22, 2019

CliqueSNV: Scalable reconstruction of intra-host viral populations from NGS reads

Highly mutable RNA viruses such as influenza A virus, human immunodeficiency virus and hepatitis C virus exist in infected hosts as highly heterogeneous populations of closely related genomic variants. The presence of low-frequency variants with few mutations with respect to major strains may result in an immune escape, emergence of drug resistance, and an increase of virulence and infectivity. Next-generation sequencing technologies permit detection of sample intra-host viral population at extremely great depth, thus providing an opportunity to access low-frequency variants. Long read lengths offered by single-molecule sequencing technologies allow all viral variants to be sequenced in a single pass. However, high sequencing error rates limit the ability to study heterogeneous viral populations composed of rare, closely related variants. In this article, we present CliqueSNV, a novel reference-based method for reconstruction of viral variants from NGS data. It efficiently constructs an allele graph based on linkage between single nucleotide variations and identifies true viral variants by merging cliques of that graph using combinatorial optimization techniques. The new method outperforms existing methods in both accuracy and running time on experimental and simulated NGS data for titrated levels of known viral variants. For PacBio reads, it accurately reconstructs variants with frequency as low as 0.1%. For Illumina reads, it fully reconstructs main variants. The open source implementation of CliqueSNV is freely available for download at https://github.com/vyacheslav-tsivina/CliqueSNV


September 22, 2019

Dynamic evolution of a-gliadin prolamin gene family in homeologous genomes of hexaploid wheat.

Wheat Gli-2 loci encode complex groups of a-gliadin prolamins that are important for breadmaking, but also major triggers of celiac disease (CD). Elucidation of a-gliadin evolution provides knowledge to produce wheat with better end-use properties and reduced immunogenic potential. The Gli-2 loci contain a large number of tandemly duplicated genes and highly repetitive DNA, making sequence assembly of their genomic regions challenging. Here, we constructed high-quality sequences spanning the three wheat homeologous a-gliadin loci by aligning PacBio-based sequence contigs with BioNano genome maps. A total of 47 a-gliadin genes were identified with only 26 encoding intact full-length protein products. Analyses of a-gliadin loci and phylogenetic tree reconstruction indicate significant duplications of a-gliadin genes in the last ~2.5 million years after the divergence of the A, B and D genomes, supporting its rapid lineage-independent expansion in different Triticeae genomes. We showed that dramatic divergence in expression of a-gliadin genes could not be attributed to sequence variations in the promoter regions. The study also provided insights into the evolution of CD epitopes and identified a single indel event in the hexaploid wheat D genome that likely resulted in the generation of the highly toxic 33-mer CD epitope.


September 22, 2019

Genomic diversity of Taylorella equigenitalis introduced into the United States from 1978 to 2012.

Contagious equine metritis is a disease of worldwide concern in equids. The United States is considered to be free of the disease although sporadic outbreaks have occurred over the last few decades that were thought to be associated with the importation of horses. The objective of this study was to create finished, reference quality genomes that characterize the diversity of Taylorella equigenitalis isolates introduced into the USA, and identify their differences. Five isolates of T. equigenitalis associated with introductions into the USA from unique sources were sequenced using both short and long read chemistries allowing for complete assembly and annotation. These sequences were compared to previously published genomes as well as the short read sequences of the 200 isolates in the National Veterinary Services Laboratories’ diagnostic repository to identify unique regions and genes, potential virulence factors, and characterize diversity. The 5 genomes varied in size by up to 100,000 base pairs, but averaged 1.68 megabases. The majority of that diversity in size can be explained by repeat regions and 4 main regions of difference, which ranged in size from 15,000 to 45,000 base pairs. The first region of difference contained mostly hypothetical proteins, the second contained the CRISPR, the third contained primarily hemagglutinin proteins, and the fourth contained primarily segments of a type IV secretion system. As expected and previously reported, little evidence of recombination was found within these genomes. Several additional areas of interest were also observed including a mechanism for streptomycin resistance and other virulence factors. A SNP distance comparison of the T. equigenitalis isolates and Mycobacterium tuberculosis complex (MTBC) showed that relatively, T. equigenitalis was a more diverse species than the entirety of MTBC.


September 22, 2019

New Delhi metallo-beta-lactamase-producing Enterobacteriaceae in South Korea between 2010 and 2015.

This study was carried out to investigate the epidemiological time-course of New Delhi metallo-beta-lactamase- (NDM-) mediated carbapenem resistance in Enterobacteriaceae in South Korea. A total of 146 non-duplicate NDM-producing Enterobacteriaceae recovered between 2010 and 2015 were voluntarily collected from 33 general hospitals and confirmed by PCR. The species were identified by sequences of the 16S rDNA. Antimicrobial susceptibility was determined either by the disk diffusion method or by broth microdilution, and the carbapenem MICs were determined by agar dilution. Then, multilocus sequence typing and PCR-based replicon typing was carried out. Co-carried genes for drug resistance were identified by PCR and sequencing. The entire genomes of eight random selected NDM producers were sequenced. A total of 69 Klebsiella pneumoniae of 12 sequence types (STs), 34 Escherichia coli of 15 STs, 28 Enterobacter spp. (including one Enterobacter aerogenes), nine Citrobacter freundii, four Raoultella spp., and two Klebsiella oxytoca isolates produced either NDM-1 (n = 126), NDM-5 (n = 18), or NDM-7 (n = 2). The isolates co-produced CTX-M-type ESBL (52.1%), AmpCs (27.4%), additional carbapenemases (7.1%), and/or 16S rRNA methyltransferases (4.8%), resulting in multidrug-resistance (47.9%) or extensively drug-resistance (52.1%). Among plasmids harboring blaNDM, IncX3 was predominant (77.4%), followed by the IncFII type (5.8%). Genome analysis revealed inter-species and inter-strain horizontal gene transfer of the plasmid. Both clonal dissemination and plasmid transfer contributed to the wide dissemination of NDM producers in South Korea.


September 22, 2019

Cross-species comparison of the gut: Differential gene expression sheds light on biological differences in closely related tenebrionids.

The gut is one of the primary interfaces between an insect and its environment. Understanding gene expression profiles in the insect gut can provide insight into interactions with the environment as well as identify potential control methods for pests. We compared the expression profiles of transcripts from the gut of larval stages of two coleopteran insects, Tenebrio molitor and Tribolium castaneum. These tenebrionids have different life cycles, varying in the duration and number of larval instars. T. castaneum has a sequenced genome and has been a model for coleopterans, and we recently obtained a draft genome for T. molitor. We assembled gut transcriptome reads from each insect to their respective genomes and filtered mapped reads to RPKM>1, yielding 11,521 and 17,871 genes in the T. castaneum and T. molitor datasets, respectively. There were identical GO terms in each dataset, and enrichment analyses also identified shared GO terms. From these datasets, we compiled an ortholog list of 6907 genes; 45% of the total assembled reads from T. castaneum were found in the top 25 orthologs, but only 27% of assembled reads were found in the top 25 T. molitor orthologs. There were 2281 genes unique to T. castaneum, and 2088 predicted genes unique to T. molitor, although improvements to the T. molitor genome will likely reduce these numbers as more orthologs are identified. We highlight a few unique genes in T. castaneum or T. molitor that may relate to distinct biological functions. A large number of putative genes expressed in the larval gut with uncharacterized functions (36 and 68% from T. castaneum and T. molitor, respectively) support the need for further research. These data are the first step in building a comprehensive understanding of the physiology of the gut in tenebrionid insects, illustrating commonalities and differences that may be related to speciation and environmental adaptation. Published by Elsevier Ltd.


September 22, 2019

Genetic characterization of a VanG-type vancomycin-resistant Enterococcus faecium clinical isolate.

To characterize, phenotypically and genotypically, the first Enterococcus faecium clinical isolate harbouring a vanG operon.The antibiotic resistance profile of E. faecium 16-346 was determined and its whole genome sequenced using PacBio technology. Attempts to transfer vancomycin resistance by filter mating were performed and the inducibility of expression of the vanG operon was studied by reverse-transcription quantitative PCR (RT-qPCR) in the presence or absence of subinhibitory concentrations of vancomycin.E. faecium 16-346 was resistant to rifampicin (MIC >4?mg/L), erythromycin (MIC >4?mg/L), tetracycline (MIC >16?mg/L) and vancomycin (MIC 8?mg/L), but susceptible to teicoplanin (MIC 0.5?mg/L). The strain harboured the vanG operon in its chromosome, integrated in a 45.5?kb putative mobile genetic element, similar to that of Enterococcus faecalis BM4518. We were unable to transfer vancomycin resistance from E. faecium 16-346 to E. faecium BM4107 and E. faecalis JH2-2. Lastly, transcription of the vanG gene was inducible by vancomycin.This is, to the best of our knowledge, the first report of a VanG-type vancomycin-resistant strain of E. faecium. Despite the alarm pulled because of the therapeutic problems caused by VRE, our work shows that new resistant loci can still be found in E. faecium.


September 22, 2019

Characterization of Lactobacillus amylolyticus L6 as potential probiotics based on genome sequence and corresponding phenotypes

The potential of newly isolated Lactobacillus amylolyticus L6 as probiotics was investigated based on the whole genome sequence and corresponding phenotypes. With Lactobacillus acidophilus NCFM as positive control, several established methods of evaluating potential probiotics were performed on L. amylolyticus L6. The results indicated that L. amylolyticus L6 retained higher viability in human gastrointestinal (GI) tract and it also had strong inhibitory effect on pathogenic bacteria. Meanwhile, the candidate probiotics exhibited similar adhesion level as that of L. acidophilus NCFM in vitro test. As for carbohydrate utilization profile, L. amylolyticus L6 had high ability of utilizing raffinose and stachyose which were known as flatulence factors in soybean products. And this strain could also utilize starch. Besides, the mechanisms of probiotic and metabolic properties for L. amylolyticus L6 were further illustrated with the identification of related genes through the analysis of genome sequence. Therefore, we proposed that L. amylolyticus L6 have the potential to be used as probiotics from phenotypes to genotypes. And it is the first time that the complete genome sequence of L. amylolyticus L6 and the potential of this strain to be used as probiotics were reported in this study.


September 22, 2019

Comparative genome analysis reveals a complex population structure of Legionella pneumophila subspecies.

The majority of Legionnaires’ disease (LD) cases are caused by Legionella pneumophila, a genetically heterogeneous species composed of at least 17 serogroups. Previously, it was demonstrated that L. pneumophila consists of three subspecies: pneumophila, fraseri and pascullei. During an LD outbreak investigation in 2012, we detected that representatives of both subspecies fraseri and pascullei colonized the same water system and that the outbreak-causing strain was a new member of the least represented subspecies pascullei. We used partial sequence based typing consensus patterns to mine an international database for additional representatives of fraseri and pascullei subspecies. As a result, we identified 46 sequence types (STs) belonging to subspecies fraseri and two STs belonging to subspecies pascullei. Moreover, a recent retrospective whole genome sequencing analysis of isolates from New York State LD clusters revealed the presence of a fourth L. pneumophila subspecies that we have termed raphaeli. This subspecies consists of 15 STs. Comparative analysis was conducted using the genomes of multiple members of all four L. pneumophila subspecies. Whereas each subspecies forms a distinct phylogenetic clade within the L. pneumophila species, they share more average nucleotide identity with each other than with other Legionella species. Unique genes for each subspecies were identified and could be used for rapid subspecies detection. Improved taxonomic classification of L. pneumophila strains may help identify environmental niches and virulence attributes associated with these genetically distinct subspecies. Published by Elsevier B.V.


September 22, 2019

Genetic basis of chromosomally-encoded mcr-1 gene.

Compared with plasmid-borne mcr-1, the occurrence of chromosomally-encoded mcr-1 is rare although it has been reported in several cases. This study aimed to investigate the genetic features of chromosomally-encoded mcr-1 among Escherichia coli strains as well as the potential genetic basis governing mobilisation of mcr-1 in bacterial chromosomes. The genome sequences of 16 E. coli strains containing a chromosomal mcr-1 gene were obtained and analysed. Phylogenetic and whole-genome sequencing (WGS) analysis demonstrated that mcr-1 was associated with four major types of genetic arrangements, namely ISApl1-mcr1-orf, Tn6330, complex Tn6330 and ?Tn6330 in chromosomes of genetically unrelated E. coli strains. The mcr-1-carrying mobile elements were shown to insert into the AT-rich region, which was also the case for ISApl1. Analysis of complete E. coli genome sequences showed that there were multiple copies of ISApl1 present in E. coli chromosomes that also carried mcr-1, whilst all mcr-1-negative chromosomes were absent of any copy of ISApl1, suggesting the strong association of ISApl1 and mcr-1. Insertion of ISApl1 into E. coli chromosomes may be a prerequisite for the insertion of mcr-1-carrying mobile elements. Insertion of mcr-1 into E. coli chromosomes would enable it to become intrinsically resistant, which is expected to become more prevalent. Policy on the prudent use of colistin both in veterinary and clinical settings should be imposed globally to further prevent dissemination of mcr-1 in E. coli and other bacterial pathogens. Copyright © 2017 Elsevier B.V. and International Society of Chemotherapy. All rights reserved.


September 22, 2019

Benefit from decline: the primary transcriptome of Alteromonas macleodii str. Te101 during Trichodesmium demise.

Interactions between co-existing microorganisms deeply affect the physiology of the involved organisms and, ultimately, the function of the ecosystem as a whole. Copiotrophic Alteromonas are marine gammaproteobacteria that thrive during the late stages of phytoplankton blooms in the marine environment and in laboratory co-cultures with cyanobacteria such as Trichodesmium. The response of this heterotroph to the sometimes rapid and transient changes in nutrient supply when the phototroph crashes is not well understood. Here, we isolated and sequenced the strain Alteromonas macleodii str. Te101 from a laboratory culture of Trichodesmium erythraeum IMS101, yielding a chromosome of 4.63?Mb and a single plasmid of 237?kb. Increasing salinities to =43 ppt inhibited the growth of Trichodesmium but stimulated growth of the associated Alteromonas. We characterized the transcriptomic responses of both microorganisms and identified the complement of active transcriptional start sites in Alteromonas at single-nucleotide resolution. In replicate cultures, a similar set of genes became activated in Alteromonas when growth rates of Trichodesmium declined and mortality was high. The parallel activation of fliA, rpoS and of flagellar assembly and growth-related genes indicated that Alteromonas might have increased cell motility, growth, and multiple biosynthetic activities. Genes with the highest expression in the data set were three small RNAs (Aln1a-c) that were identified as analogs of the small RNAs CsrB-C in E. coli or RsmX-Z in pathogenic bacteria. Together with the carbon storage protein A (CsrA) homolog Te101_05290, these RNAs likely control the expression of numerous genes in responding to changes in the environment.


September 22, 2019

Autologous cell therapy approach for Duchenne muscular dystrophy using PiggyBac transposons and mesoangioblasts.

Duchenne muscular dystrophy (DMD) is a lethal muscle-wasting disease currently without cure. We investigated the use of the PiggyBac transposon for full-length dystrophin expression in murine mesoangioblast (MABs) progenitor cells. DMD murine MABs were transfected with transposable expression vectors for full-length dystrophin and transplanted intramuscularly or intra-arterially into mdx/SCID mice. Intra-arterial delivery indicated that the MABs could migrate to regenerating muscles to mediate dystrophin expression. Intramuscular transplantation yielded dystrophin expression in 11%-44% of myofibers in murine muscles, which remained stable for the assessed period of 5 months. The satellite cells isolated from transplanted muscles comprised a fraction of MAB-derived cells, indicating that the transfected MABs may colonize the satellite stem cell niche. Transposon integration site mapping by whole-genome sequencing indicated that 70% of the integrations were intergenic, while none was observed in an exon. Muscle resistance assessment by atomic force microscopy indicated that 80% of fibers showed elasticity properties restored to those of wild-type muscles. As measured in vivo, transplanted muscles became more resistant to fatigue. This study thus provides a proof-of-principle that PiggyBac transposon vectors may mediate full-length dystrophin expression as well as functional amelioration of the dystrophic muscles within a potential autologous cell-based therapeutic approach of DMD. Copyright © 2018 The American Society of Gene and Cell Therapy. Published by Elsevier Inc. All rights reserved.


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.