Emergence and Evolution of Multidrug-Resistant Klebsiella pneumoniae with both blaKPC and blaCTX-M Integrated in the Chromosome

The extended-spectrum-β-lactamase (ESBL)- and Klebsiella pneumoniaecarbapenemase (KPC)-producing Enterobacteriaceae represent serious and urgent threats to public health. In a retrospective study of multidrug-resistant K. pneumoniae, we identified three clinical isolates, CN1, CR14, and NY9, carrying both blaCTX-M and blaKPC genes. The complete genomes of these three K. pneumoniae isolates were de novo assembled by using both short- and long-read whole-genome sequencing. In CR14 and NY9, blaCTX-M and blaKPC were carried on two different plasmids. In contrast, CN1 had one copy of blaKPC-2 and three copies of blaCTX-M-15 integrated in the chromosome, for which the blaCTX-M-15 genes were linked to an insertion sequence, ISEcp1, whereas the blaKPC-2 gene was in the context of a Tn4401a transposition unit conjugated with a PsP3-like prophage. Intriguingly, downstream of the Tn4401a-blaKPC-2-prophage genomic island, CN1 also carried a clustered regularly interspaced short palindromic repeat (CRISPR)-cas array with four spacers targeting a variety of K. pneumoniae plasmids harboring antimicrobial resistance genes. Comparative genomic analysis revealed that there were two subtypes of type I-E CRISPR-cas in K. pneumoniae strains and suggested that the evolving CRISPR-cas, with its acquired novel spacer, induced the mobilization of antimicrobial resistance genes from plasmids into the chromosome. The integration and dissemination of multiple copies of blaCTX-Mand blaKPC from plasmids to chromosome depicts the complex pandemic scenario of multidrug-resistant K. pneumoniae. Additionally, the implications from this study also raise concerns for the application of a CRISPR-cas strategy against antimicrobial resistance

Single-Molecule Real-Time (SMRT) Full-Length RNA-Sequencing Reveals Novel and Distinct mRNA Isoforms in Human Bone Marrow Cell Subpopulations.

Hematopoietic cells are continuously replenished from progenitor cells that reside in the bone marrow. To evaluate molecular changes during this process, we analyzed the transcriptomes of freshly harvested human bone marrow progenitor (lineage-negative) and differentiated (lineage-positive) cells by single-molecule real-time (SMRT) full-length RNA-sequencing. This analysis revealed a ~5-fold higher number of transcript isoforms than previously detected and showed a distinct composition of individual transcript isoforms characteristic for bone marrow subpopulations. A detailed analysis of messenger RNA (mRNA) isoforms transcribed from the ANXA1 and EEF1A1 loci confirmed their distinct composition. The expression of proteins predicted from the transcriptome analysis was evaluated by mass spectrometry and validated previously unknown protein isoforms predicted e.g., for EEF1A1. These protein isoforms distinguished the lineage negative cell population from the lineage positive cell population. Finally, transcript isoforms expressed from paralogous gene loci (e.g., CFD, GATA2, HLA-A, B, and C) also distinguished cell subpopulations but were only detectable by full-length RNA sequencing. Thus, qualitatively distinct transcript isoforms from individual genomic loci separate bone marrow cell subpopulations indicating complex transcriptional regulation and protein isoform generation during hematopoiesis

Comparative Analysis of Tandem Repeats from Hundreds of Species Reveals Unique Insights into Centromere Evolution

Centromeres are essential for chromosome segregation, yet their DNA sequences evolve rapidly. In most animals and plants that have been studied, centromeres contain megabase-scale arrays of tandem repeats. Despite their importance, very little is known about the degree to which centromere tandem repeats share common properties between different species across different phyla. We used bioinformatic methods to identify high-copy tandem repeats from 282 species using publicly available genomic sequence and our own data. The assumption that the most abundant tandem repeat is the centromere DNA was true for most species whose centromeres have been previously characterized, suggesting this is a general property of genomes. Our methods are compatible with all current sequencing technologies. Long Pacific Biosciences sequence reads allowed us to find tandem repeat monomers up to 1,419 bp. High-copy centromere tandem repeats were found in almost all animal and plant genomes, but repeat monomers were highly variable in sequence composition and in length. Furthermore, phylogenetic analysis of sequence homology showed little evidence of sequence conservation beyond ~50 million years of divergence. We find that despite an overall lack of sequence conservation, centromere tandem repeats from diverse species showed similar modes of evolution, including the appearance of higher order repeat structures in which several polymorphic monomers make up a larger repeating unit. While centromere position in most eukaryotes is epigenetically determined, our results indicate that tandem repeats are highly prevalent at centromeres of both animals and plants. This suggests a functional role for such repeats, perhaps in promoting concerted evolution of centromere DNA across chromosomes

The impact of sequencing depth on the inferred taxonomic composition and AMR gene content of metagenomic samples

Shotgun metagenomics is increasingly used to characterise microbial communities, particularly for the investigation of antimicrobial resistance (AMR) in different animal and environmental contexts. There are many different approaches for inferring the taxonomic composition and AMR gene content of complex community samples from shotgun metagenomic data, but there has been little work establishing the optimum sequencing depth, data processing and analysis methods for these samples. In this study we used shotgun metagenomics and sequencing of cultured isolates from the same samples to address these issues. We sampled three potential environmental AMR gene reservoirs (pig caeca, river sediment, effluent) and sequenced samples with shotgun metagenomics at high depth (~ 200 million reads per sample). Alongside this, we cultured single-colony isolates of Enterobacteriaceae from the same samples and used hybrid sequencing (short- and long-reads) to create high- quality assemblies for comparison to the metagenomic data. To automate data processing, we developed an open- source software pipeline, ‘ResPipe’

Capture of complete ciliate chromosomes in single sequencing reads reveals widespread chromosome isoforms

Background Whole-genome shotgun sequencing, which stitches together millions of short sequencing reads into a single genome, ushered in the era of modern genomics and led to a rapid expansion of the number of genome sequences available. Nevertheless, assembly of short reads remains difficult, resulting in fragmented genome sequences. Ultimately, only a sequencing technology capable of capturing complete chromosomes in a single run could resolve all ambiguities. Even “third generation” sequencing technologies produce reads far shorter than most eukaryotic chromosomes. However, the ciliate Oxytricha trifallax has a somatic genome with thousands of chromosomes averaging only 3.2 kbp, making it an ideal candidate for exploring the benefits of sequencing whole chromosomes without assembly. Results We used single-molecule real-time sequencing to capture thousands of complete chromosomes in single reads and to update the published Oxytricha trifallax JRB310 genome assembly. In this version, over 50% of the completed chromosomes with two telomeres derive from single reads. The improved assembly includes over 12,000 new chromosome isoforms, and demonstrates that somatic chromosomes derive from variable rearrangements between somatic segments encoded up to 191,000 base pairs away. However, while long reads reduce the need for assembly, a hybrid approach that supplements long-read sequencing with short reads for error correction produced the most complete and accurate assembly, overall. Conclusions This assembly provides the first example of complete eukaryotic chromosomes captured by single sequencing reads and demonstrates that traditional approaches to genome assembly can mask considerable structural variation

The Complete Genome Sequence of the Emerging Pathogen Mycobacterium haemophilum Explains Its Unique Culture Requirements

Mycobacterium haemophilum is an emerging pathogen associated with a variety of clinical syndromes, most commonly skin infections in immunocompromised individuals. M. haemophilum exhibits a unique requirement for iron supplementation to support its growth in culture, but the basis for this property and how it may shape pathogenesis is unclear. Using a combination of Illumina, PacBio, and Sanger sequencing, the complete genome sequence of M. haemophilum was determined. Guided by this sequence, experiments were performed to define the basis for the unique growth requirements of M. haemophilum. We found that M. haemophilum, unlike many other mycobacteria, is unable to synthesize iron-binding siderophores known as mycobactins or to utilize ferri-mycobactins to support growth. These differences correlate with the absence of genes associated with mycobactin synthesis, secretion, and uptake. In agreement with the ability of heme to promote growth, we identified genes encoding heme uptake machinery. Consistent with its propensity to infect the skin, we show at the whole-genome level the genetic closeness of M. haemophilumwith Mycobacterium leprae, an organism which cannot be cultivated in vitro, and we identify genes uniquely shared by these organisms. Finally, we identify means to express foreign genes in M. haemophilum. These data explain the unique culture requirements for this important pathogen, provide a foundation upon which the genome sequence can be exploited to improve diagnostics and therapeutics, and suggest use of M. haemophilum as a tool to elucidate functions of genes shared with M. leprae. IMPORTANCE Mycobacterium haemophilum is an emerging pathogen with an unknown natural reservoir that exhibits unique requirements for iron supplementation to grow in vitro. Understanding the basis for this iron requirement is important because it is fundamental to isolation of the organism from clinical samples and environmental sources. Defining the molecular basis for M. haemophilium\u27s growth requirements will also shed new light on mycobacterial strategies to acquire iron and can be exploited to define how differences in such strategies influence pathogenesis. Here, through a combination of sequencing and experimental approaches, we explain the basis for the iron requirement. We further demonstrate the genetic closeness of M. haemophilum and Mycobacterium leprae, the causative agent of leprosy which cannot be cultured in vitro, and we demonstrate methods to genetically manipulate M. haemophilum. These findings pave the way for the use of M. haemophilum as a model to elucidate functions of genes shared with M. leprae

Analysis of the Aedes albopictus C6/36 genome provides insight into cell line utility for viral propagation

BACKGROUND: The 50-year-old Aedes albopictus C6/36 cell line is a resource for the detection, amplification, and analysis of mosquito-borne viruses including Zika, dengue, and chikungunya. The cell line is derived from an unknown number of larvae from an unspecified strain of Aedes albopictus mosquitoes. Toward improved utility of the cell line for research in virus transmission, we present an annotated assembly of the C6/36 genome. RESULTS: The C6/36 genome assembly has the largest contig N50 (3.3 Mbp) of any mosquito assembly, presents the sequences of both haplotypes for most of the diploid genome, reveals independent null mutations in both alleles of the Dicer locus, and indicates a male-specific genome. Gene annotation was computed with publicly available mosquito transcript sequences. Gene expression data from cell line RNA sequence identified enrichment of growth-related pathways and conspicuous deficiency in aquaporins and inward rectifier K+ channels. As a test of utility, RNA sequence data from Zika-infected cells were mapped to the C6/36 genome and transcriptome assemblies. Host subtraction reduced the data set by 89%, enabling faster characterization of nonhost reads. CONCLUSIONS: The C6/36 genome sequence and annotation should enable additional uses of the cell line to study arbovirus vector interactions and interventions aimed at restricting the spread of human disease

Genomic Epidemiology of Complex, Multispecies, Plasmid-Borne bla KPC Carbapenemase in Enterobacterales in the United Kingdom from 2009 to 2014.

Carbapenem resistance in Enterobacterales is a public health threat. Klebsiella pneumoniae carbapenemase (encoded by alleles of the bla KPC family) is one of the most common transmissible carbapenem resistance mechanisms worldwide. The dissemination of bla KPC historically has been associated with distinct K. pneumoniae lineages (clonal group 258 [CG258]), a particular plasmid family (pKpQIL), and a composite transposon (Tn4401). In the United Kingdom, bla KPC has represented a large-scale, persistent management challenge for some hospitals, particularly in North West England. The dissemination of bla KPC has evolved to be polyclonal and polyspecies, but the genetic mechanisms underpinning this evolution have not been elucidated in detail; this study used short-read whole-genome sequencing of 604 bla KPC-positive isolates (Illumina) and long-read assembly (PacBio)/polishing (Illumina) of 21 isolates for characterization. We observed the dissemination of bla KPC (predominantly bla KPC-2; 573/604 [95%] isolates) across eight species and more than 100 known sequence types. Although there was some variation at the transposon level (mostly Tn4401a, 584/604 [97%] isolates; predominantly with ATTGA-ATTGA target site duplications, 465/604 [77%] isolates), bla KPC spread appears to have been supported by highly fluid, modular exchange of larger genetic segments among plasmid populations dominated by IncFIB (580/604 isolates), IncFII (545/604 isolates), and IncR (252/604 isolates) replicons. The subset of reconstructed plasmid sequences (21 isolates, 77 plasmids) also highlighted modular exchange among non-bla KPC and bla KPC plasmids and the common presence of multiple replicons within bla KPC plasmid structures (>60%). The substantial genomic plasticity observed has important implications for our understanding of the epidemiology of transmissible carbapenem resistance in Enterobacterales for the implementation of adequate surveillance approaches and for control

A personalized platform identifies trametinib plus zoledronate for a patient with KRAS-mutant metastatic colorectal cancer

Colorectal cancer remains a leading source of cancer mortality worldwide. Initial response is often followed by emergent resistance that is poorly responsive to targeted therapies, reflecting currently undruggable cancer drivers such as KRAS and overall genomic complexity. Here, we report a novel approach to developing a personalized therapy for a patient with treatment-resistant metastatic KRAS-mutant colorectal cancer. An extensive genomic analysis of the tumor's genomic landscape identified nine key drivers. A transgenic model that altered orthologs of these nine genes in the Drosophila hindgut was developed; a robotics-based screen using this platform identified trametinib plus zoledronate as a candidate treatment combination. Treating the patient led to a significant response: Target and nontarget lesions displayed a strong partial response and remained stable for 11 months. By addressing a disease's genomic complexity, this personalized approach may provide an alternative treatment option for recalcitrant disease such as KRAS-mutant colorectal cancer