Similar evolutionary trajectories in an environmental Cryptococcus neoformans isolate after human and murine infection.

A pet cockatoo was the suspected source of Cryptococcus neoformans recovered from an immunocompromised patient with cryptococcosis based on molecular analyses available in 2000. Here, we report whole genome sequence analysis of the clinical and cockatoo strains. Both are closely related MATα strains belonging to the VNII lineage, confirming that the human infection likely originated from pet bird exposure. The two strains differ by 61 single nucleotide polymorphisms, including eight nonsynonymous changes involving seven genes. To ascertain whether changes in these genes are selected for during mammalian infection, we passaged the cockatoo strain in mice. Remarkably, isolates obtained from mouse tissue possess a frameshift mutation in one of the seven genes altered in the human sample (LQVO5_000317), a gene predicted to encode an SWI-SNF chromatin-remodeling complex protein. In addition, both cockatoo and patient strains as well as mouse-passaged isolates obtained from brain tissue had a premature stop codon in a homologue of ZFC3 (LQVO5_004463), a predicted single-zinc finger containing protein, which is associated with larger capsules when deleted and reverted to a full-length protein in the mouse-passaged isolates obtained from lung tissue. The patient strain and mouse-passaged isolates show variability in virulence factors, with differences in capsule size, melanization, rates of nonlytic expulsion from macrophages, and amoeba predation resistance. Our results establish that environmental strains undergo genomic and phenotypic changes during mammalian passage, suggesting that animal virulence can be a mechanism for genetic change and that the genomes of clinical isolates may provide a readout of mutations acquired during infection.

A plasmid locus associated with Klebsiella clinical infections encodes a microbiome-dependent gut fitness factor.

Klebsiella pneumoniae (Kp) is an important cause of healthcare-associated infections, which increases patient morbidity, mortality, and hospitalization costs. Gut colonization by Kp is consistently associated with subsequent Kp disease, and patients are predominantly infected with their colonizing strain. Our previous comparative genomics study, between disease-causing and asymptomatically colonizing Kp isolates, identified a plasmid-encoded tellurite (TeO3-2)-resistance (ter) operon as strongly associated with infection. However, TeO3-2 is extremely rare and toxic to humans. Thus, we used a multidisciplinary approach to determine the biological link between ter and Kp infection. First, we used a genomic and bioinformatic approach to extensively characterize Kp plasmids encoding the ter locus. These plasmids displayed substantial variation in plasmid incompatibility type and gene content. Moreover, the ter operon was genetically independent of other plasmid-encoded virulence and antibiotic resistance loci, both in our original patient cohort and in a large set (n = 88) of publicly available ter operon-encoding Kp plasmids, indicating that the ter operon is likely playing a direct, but yet undescribed role in Kp disease. Next, we employed multiple mouse models of infection and colonization to show that 1) the ter operon is dispensable during bacteremia, 2) the ter operon enhances fitness in the gut, 3) this phenotype is dependent on the colony of origin of mice, and 4) antibiotic disruption of the gut microbiota eliminates the requirement for ter. Furthermore, using 16S rRNA gene sequencing, we show that the ter operon enhances Kp fitness in the gut in the presence of specific indigenous microbiota, including those predicted to produce short chain fatty acids. Finally, administration of exogenous short-chain fatty acids in our mouse model of colonization was sufficient to reduce fitness of a ter mutant. These findings indicate that the ter operon, strongly associated with human infection, encodes factors that resist stress induced by the indigenous gut microbiota during colonization. This work represents a substantial advancement in our molecular understanding of Kp pathogenesis and gut colonization, directly relevant to Kp disease in healthcare settings.

Network Analysis of Protein Adaptation: Modeling the Functional Impact of Multiple Mutations.

The evolution of new biochemical activities frequently involves complex dependencies between mutations and rapid evolutionary radiation. Mutation co-occurrence and covariation have previously been used to identify compensating mutations that are the result of physical contacts and preserve protein function and fold. Here, we model pairwise functional dependencies and higher order interactions that enable evolution of new protein functions. We use a network model to find complex dependencies between mutations resulting from evolutionary trade-offs and pleiotropic effects. We present a method to construct these networks and to identify functionally interacting mutations in both extant and reconstructed ancestral sequences (Network Analysis of Protein Adaptation). The time ordering of mutations can be incorporated into the networks through phylogenetic reconstruction. We apply NAPA to three distantly homologous ß-lactamase protein clusters (TEM, CTX-M-3, and OXA-51), each of which has experienced recent evolutionary radiation under substantially different selective pressures. By analyzing the network properties of each protein cluster, we identify key adaptive mutations, positive pairwise interactions, different adaptive solutions to the same selective pressure, and complex evolutionary trajectories likely to increase protein fitness. We also present evidence that incorporating information from phylogenetic reconstruction and ancestral sequence inference can reduce the number of spurious links in the network, whereas preserving overall network community structure. The analysis does not require structural or biochemical data. In contrast to function-preserving mutation dependencies, which are frequently from structural contacts, gain-of-function mutation dependencies are most commonly between residues distal in protein structure.

Applying Rapid Whole-Genome Sequencing To Predict Phenotypic Antimicrobial Susceptibility Testing Results among Carbapenem-Resistant Klebsiella pneumoniae Clinical Isolates.

Standard antimicrobial susceptibility testing (AST) approaches lead to delays in the selection of optimal antimicrobial therapy. Here, we sought to determine the accuracy of antimicrobial resistance (AMR) determinants identified by Nanopore whole-genome sequencing in predicting AST results. Using a cohort of 40 clinical isolates (21 carbapenemase-producing carbapenem-resistant Klebsiella pneumoniae, 10 non-carbapenemase-producing carbapenem-resistant K. pneumoniae, and 9 carbapenem-susceptible K. pneumoniae isolates), three separate sequencing and analysis pipelines were performed, as follows: (i) a real-time Nanopore analysis approach identifying acquired AMR genes, (ii) an assembly-based Nanopore approach identifying acquired AMR genes and chromosomal mutations, and (iii) an approach using short-read correction of Nanopore assemblies. The short-read correction of Nanopore assemblies served as the reference standard to determine the accuracy of Nanopore sequencing results. With the real-time analysis approach, full annotation of acquired AMR genes occurred within 8 h from subcultured isolates. Assemblies sufficient for full resistance gene and single-nucleotide polymorphism annotation were available within 14 h from subcultured isolates. The overall agreement of genotypic results and anticipated AST results for the 40 K. pneumoniae isolates was 77% (range, 30% to 100%) and 92% (range, 80% to 100%) for the real-time approach and the assembly approach, respectively. Evaluating the patients contributing the 40 isolates, the real-time approach and assembly approach could shorten the median time to effective antibiotic therapy by 20 h and 26 h, respectively, compared to standard AST. Nanopore sequencing offers a rapid approach to both accurately identify resistance mechanisms and to predict AST results for K. pneumoniae isolates. Bioinformatics improvements enabling real-time alignment, coupled with rapid extraction and library preparation, will further enhance the accuracy and workflow of the Nanopore real-time approach.

Multiplexed analysis of fixed tissue RNA using Ligation in situ Hybridization.

Clinical tissues are prepared for histological analysis and long-term storage via formalin fixation and paraffin embedding (FFPE). The FFPE process results in fragmentation and chemical modification of RNA, rendering it less suitable for analysis by techniques that rely on reverse transcription (RT) such as RT-qPCR and RNA-Seq. Here we describe a broadly applicable technique called 'Ligation in situ Hybridization' ('LISH'), which is an alternative methodology for the analysis of FFPE RNA. LISH utilizes the T4 RNA Ligase 2 to efficiently join adjacent chimeric RNA-DNA probe pairs hybridized in situ on fixed RNA target sequences. Subsequent treatment with RNase H releases RNA-templated ligation products into solution for downstream analysis. We demonstrate several unique advantages of LISH-based assays using patient-derived FFPE tissue. These include >100-plex capability, compatibility with common histochemical stains and suitability for analysis of decade-old materials and exceedingly small microdissected tissue fragments. High-throughput DNA sequencing modalities, including single molecule sequencing, can be used to analyze ligation products from complex panels of LISH probes ('LISH-seq'), which can be amplified efficiently and with negligible bias. LISH analysis of FFPE RNA is a novel methodology with broad applications that range from multiplexed gene expression analysis to the sensitive detection of infectious organisms.

Successful Treatment of Persistent Burkholderia cepacia Complex Bacteremia with Ceftazidime-Avibactam.

We report our clinical experience treating a 2-month-old infant with congenital diaphragmatic hernia who experienced prolonged bacteremia with Burkholderia cepacia complex (Bcc) despite conventional antibiotic therapy and appropriate source control measures. The infection resolved after initiation of ceftazidime-avibactam. Whole-genome sequencing revealed that the isolate most closely resembled B. contaminans and identified the mechanism of resistance that likely contributed to clinical cure with this agent. Ceftazidime-avibactam should be considered salvage therapy for Bcc infections if other treatment options have been exhausted.

Comprehensive Screening of a Light-Inducible Split Cre Recombinase with Domain Insertion Profiling.

Splitting proteins with light- or chemically inducible dimers provides a mechanism for post-translational control of protein function. However, current methods for engineering stimulus-responsive split proteins often require significant protein engineering expertise and the laborious screening of individual constructs. To address this challenge, we use a pooled library approach that enables rapid generation and screening of nearly all possible split protein constructs in parallel, where results can be read out by using sequencing. We perform our method on Cre recombinase with optogenetic dimers as a proof of concept, resulting in comprehensive data on the split sites throughout the protein. To improve the accuracy in predicting split protein behavior, we develop a Bayesian computational approach to contextualize errors inherent to experimental procedures. Overall, our method provides a streamlined approach for achieving inducible post-translational control of a protein of interest.

Comprehensive screening of a light-inducible split Cre recombinase with domain insertion profiling.

Splitting proteins with light- or chemically-inducible dimers provides a mechanism for post-translational control of protein function. However, current methods for engineering stimulus-responsive split proteins often require significant protein engineering expertise and laborious screening of individual constructs. To address this challenge, we use a pooled library approach that enables rapid generation and screening of nearly all possible split protein constructs in parallel, where results can be read out using sequencing. We perform our method on Cre recombinase with optogenetic dimers as a proof of concept, resulting in comprehensive data on split sites throughout the protein. To improve accuracy in predicting split protein behavior, we develop a Bayesian computational approach to contextualize errors inherent to experimental procedures. Overall, our method provides a streamlined approach for achieving inducible post-translational control of a protein of interest.

Correlations among the plasma concentrations of first-line anti-tuberculosis drugs and the physiological parameters influencing concentrations.

Background: The plasma concentrations of the four most commonly used first-line anti-tuberculosis (TB) drugs, isoniazid (INH), rifampicin (RMP), ethambutol (EMB), and pyrazinamide (PZA), are often not within the therapeutic range. Insufficient drug exposure could lead to drug resistance and treatment failure, while excessive drug levels may lead to adverse reactions. The purpose of this study was to identify the physiological parameters influencing anti-TB drug concentrations. Methods: A retrospective cohort study was conducted. The 2-h plasma concentrations of the four drugs were measured by using the high-performance liquid chromatography-tandem mass spectrometry method. Results: A total of 317 patients were included in the study. The proportions of patients with INH, RMP, EMB, and PZA concentrations within the therapeutic range were 24.3%, 31.5%, 27.8%, and 18.6%, respectively. There were positive associations between the concentrations of INH and PZA and RMP and EMB, but negative associations were observed between the concentrations of INH and RMP, INH and EMB, RMP and PZA, and EMB and PZA. In the multivariate analysis, the influencing factors of the INH concentration were the PZA concentration, total bile acid (TBA), serum potassium, dose, direct bilirubin, prealbumin (PA), and albumin; those of the RMP concentration were PZA and EMB concentrations, weight, α-l-fucosidase (AFU), drinking, and dose; those of the EMB concentration were the RMP and PZA concentrations, creatinine, TBA and indirect bilirubin; and those of the PZA concentration were INH, RMP and EMB concentrations, sex, weight, uric acid and drinking. Conclusion: The complex correlations between the concentrations of the four first-line anti-TB drugs lead to a major challenge in dose adjustment to maintain all drugs within the therapeutic window. Levels of TBA, PA, AFU, and serum potassium should also be considered when adjusting the dose of the four drugs.

Multicentre genetic diversity study of carbapenem-resistant Enterobacterales: predominance of untypeable pUVA-like blaKPC bearing plasmids.

Objectives: Carbapenem-resistant Enterobacterales (CRE) are an urgent public health threat. A better understanding of the molecular epidemiology and transmission dynamics of CRE is necessary to limit their dissemination within healthcare settings. We sought to investigate the mechanisms of resistance and spread of CRE within multiple hospitals in Maryland. Methods: From 2016 to 2018, all CRE were collected from any specimen source from The Johns Hopkins Medical Institutions. The isolates were further characterized using both phenotypic and genotypic approaches, including short- and/or long-read WGS. Results: From 2016 to 2018, 302 of 40 908 (0.7%) unique Enterobacterales isolates were identified as CRE. Of CRE, 142 (47%) were carbapenemase-producing CRE with KPC (80.3%) predominating among various genera. Significant genetic diversity was identified among all CRE with high-risk clones serving as major drivers of clonal clusters. Further, we found the predominance of pUVA-like plasmids, with a subset harbouring resistance genes to environmental cleaning agents, involved in intergenus dissemination of blaKPC genes. Conclusions: Our findings provide valuable data to understand the transmission dynamics of all CRE within the greater Maryland region. These data can help guide targeted interventions to limit CRE transmission in healthcare facilities.

Treatment of Sindbis Virus-Infected Neurons with Antibody to E2 Alters Synthesis of Complete and nsP1-Expressing Defective Viral RNAs.

Alphaviruses are positive-sense RNA viruses that are important causes of viral encephalomyelitis. Sindbis virus (SINV), the prototype alphavirus, preferentially infects neurons in mice and is a model system for studying mechanisms of viral clearance from the nervous system. Antibody specific to the SINV E2 glycoprotein plays an important role in SINV clearance, and this effect is reproduced in cultures of infected mature neurons. To determine how anti-E2 antibody affects SINV RNA synthesis, Oxford Nanopore Technologies direct long-read RNA sequencing was used to sequence viral RNAs following antibody treatment of infected neurons. Differentiated AP-7 rat olfactory neuronal cells, an in vitro model for mature neurons, were infected with SINV and treated with anti-E2 antibody. Whole-cell RNA lysates were collected for sequencing of poly(A)-selected RNA 24, 48, and 72 h after infection. Three primary species of viral RNA were produced: genomic, subgenomic, and defective viral genomes (DVGs) encoding the RNA capping protein nsP1. Antibody treatment resulted in overall lower production of SINV RNA, decreased synthesis of subgenomic RNA relative to genomic RNA, and suppressed production of the nsP1 DVG. The nsP1 DVG was packaged into virus particles and could be translated. Because antibody-treated cells released a higher proportion of virions with noncapped genomes and transient transfection to express the nsP1 DVG improved viral RNA capping in antibody-treated cells, we postulate that one mechanism by which antibody inhibits SINV replication in neurons is to suppress DVG synthesis and thus decrease production of infectious virions containing capped genomes. IMPORTANCE Alphaviruses are important causes of viral encephalomyelitis without approved treatments or vaccines. Antibody to the Sindbis virus (SINV) E2 glycoprotein is required for immune-mediated noncytolytic virus clearance from neurons. We used direct RNA nanopore sequencing to evaluate how anti-E2 antibody affects SINV replication at the RNA level. Antibody altered the viral RNAs produced by decreasing the proportion of subgenomic relative to genomic RNA and suppressing production of a previously unrecognized defective viral genome (DVG) encoding nsP1, the viral RNA capping enzyme. Antibody-treated neurons released a lower proportion of SINV particles with capped genomes necessary for translation and infection. Decreased nsP1 DVG production in antibody-treated neurons led to lower expression of nsP1 protein, decreased genome capping efficiency, and release of fewer infectious virus particles. Capping was increased with exogenous expression of the nsP1 DVG. These studies identify a novel alphavirus DVG function and new mechanism for antibody-mediated control of virus replication.

Targeted nanopore sequencing by real-time mapping of raw electrical signal with UNCALLED.

Conventional targeted sequencing methods eliminate many of the benefits of nanopore sequencing, such as the ability to accurately detect structural variants or epigenetic modifications. The ReadUntil method allows nanopore devices to selectively eject reads from pores in real time, which could enable purely computational targeted sequencing. However, this requires rapid identification of on-target reads while most mapping methods require computationally intensive basecalling. We present UNCALLED ( https://github.com/skovaka/UNCALLED ), an open source mapper that rapidly matches streaming of nanopore current signals to a reference sequence. UNCALLED probabilistically considers k-mers that could be represented by the signal and then prunes the candidates based on the reference encoded within a Ferragina-Manzini index. We used UNCALLED to deplete sequencing of known bacterial genomes within a metagenomics community, enriching the remaining species 4.46-fold. UNCALLED also enriched 148 human genes associated with hereditary cancers to 29.6× coverage using one MinION flowcell, enabling accurate detection of single-nucleotide polymorphisms, insertions and deletions, structural variants and methylation in these genes.

De novo genome assembly of the tobacco hornworm moth (Manduca sexta).

The tobacco hornworm, Manduca sexta, is a lepidopteran insect that is used extensively as a model system for studying insect biology, development, neuroscience, and immunity. However, current studies rely on the highly fragmented reference genome Msex_1.0, which was created using now-outdated technologies and is hindered by a variety of deficiencies and inaccuracies. We present a new reference genome for M. sexta, JHU_Msex_v1.0, applying a combination of modern technologies in a de novo assembly to increase continuity, accuracy, and completeness. The assembly is 470 Mb and is ∼20× more continuous than the original assembly, with scaffold N50 > 14 Mb. We annotated the assembly by lifting over existing annotations and supplementing with additional supporting RNA-based data for a total of 25,256 genes. The new reference assembly is accessible in annotated form for public use. We demonstrate that improved continuity of the M. sexta genome improves resequencing studies and benefits future research on M. sexta as a model organism.

Genome and transcriptome of a pathogenic yeast, Candida nivariensis.

We present a highly contiguous genome and transcriptome of the pathogenic yeast, Candida nivariensis. We sequenced both the DNA and RNA of this species using both the Oxford Nanopore Technologies and Illumina platforms. We assembled the genome into an 11.8 Mb draft composed of 16 contigs with an N50 of 886 Kb, including a circular mitochondrial sequence of 28 Kb. Using direct RNA nanopore sequencing and Illumina cDNA sequencing, we constructed an annotation of our new assembly, supplemented by lifting over genes from Saccharomyces cerevisiae and Candida glabrata.

Screening of 20 Mycobacterium tuberculosis sRNAs in plasma for detection of active pulmonary tuberculosis.

BACKGROUND: Mycobacterium tuberculosis (MTB) sRNAs are abundant. However, the level of MTB sRNA in peripheral blood remains elusive. METHODS: Twenty MTB sRNAs annotated in the reference genome of H37Rv were detected in the plasma of 170 active pulmonary tuberculosis patients and 124 healthy people by qRT-PCR detection system. The differential expression of sRNAs were analyzed in two groups. The value of sRNAs for diagnosis of active tuberculosis were evaluated by ROC curve analysis. RESULTS: Eight of the 20 sRNAs (MTS2823, MTS0997, MTS1338, ASdes, G2, C8, mcr15 and MTS1082) were found in at least 50% of the samples detected. The relative expression levels of MTS2823, MTS0997, MTS1338 and ASdes in plasma of tuberculosis patients were statistically higher than those in healthy controls. ROC curve analysis showed that the AUC of MTS0997, MTS1338, MTS2823 and ASdes were 0.8935 (95% CI 0.8109-0.9760), 0.8722 (95% CI 0.7862-0.9581), 0.8208 (95% CI 0.7246-0.9170) and 0.5792 (95% CI 0.4240-0.7344), respectively. The AUC value of combination of MTS0997, MTS1338 and MTS2823 was 0.914 (95% CI 0.8281-0.9926). CONCLUSIONS: MTB sRNAs MTS2823, MTS0997 and MTS1338 have the potential to be plasma biomarkers for active pulmonary tuberculosis.

Genomic diversity of SARS-CoV-2 during early introduction into the Baltimore-Washington metropolitan area.

The early COVID-19 pandemic was characterized by rapid global spread. In Maryland and Washington, DC, United States, more than 2500 cases were reported within 3 weeks of the first COVID-19 detection in March 2020. We aimed to use genomic sequencing to understand the initial spread of SARS-CoV-2 - the virus that causes COVID-19 - in the region. We analyzed 620 samples collected from the Johns Hopkins Health System during March 11-31, 2020, comprising 28.6% of the total cases in Maryland and Washington, DC. From these samples, we generated 114 complete viral genomes. Analysis of these genomes alongside a subsampling of over 1000 previously published sequences showed that the diversity in this region rivaled global SARS-CoV-2 genetic diversity at that time and that the sequences belong to all of the major globally circulating lineages, suggesting multiple introductions into the region. We also analyzed these regional SARS-CoV-2 genomes alongside detailed clinical metadata and found that clinically severe cases had viral genomes belonging to all major viral lineages. We conclude that efforts to control local spread of the virus were likely confounded by the number of introductions into the region early in the epidemic and the interconnectedness of the region as a whole.

Genomic Diversity of SARS-CoV-2 During Early Introduction into the United States National Capital Region.

BACKGROUND: The early COVID-19 pandemic has been characterized by rapid global spread. In the United States National Capital Region, over 2,000 cases were reported within three weeks of its first detection in March 2020. We aimed to use genomic sequencing to understand the initial spread of SARS-CoV-2, the virus that causes COVID-19, in the region. By correlating genetic information to disease phenotype, we also aimed to gain insight into any correlation between viral genotype and case severity or transmissibility. METHODS: We performed whole genome sequencing of clinical SARS-CoV-2 samples collected in March 2020 by the Johns Hopkins Health System. We analyzed these regional SARS-CoV-2 genomes alongside detailed clinical metadata and the global phylogeny to understand early establishment of the virus within the region. RESULTS: We analyzed 620 samples from the Johns Hopkins Health System collected between March 11-31, 2020, comprising 37.3% of the total cases in Maryland during this period. We selected 143 of these samples for sequencing, generating 114 complete viral genomes. These genomes belong to all five major Nextstrain-defined clades, suggesting multiple introductions into the region and underscoring the diversity of the regional epidemic. We also found that clinically severe cases had genomes belonging to all of these clades. CONCLUSIONS: We established a pipeline for SARS-CoV-2 sequencing within the Johns Hopkins Health system, which enabled us to capture the significant viral diversity present in the region as early as March 2020. Efforts to control local spread of the virus were likely confounded by the number of introductions into the region early in the epidemic and interconnectedness of the region as a whole.

Transmission and clearance of potential procarcinogenic bacteria during fecal microbiota transplantation for recurrent Clostridioides difficile.

BACKGROUNDFecal microbiota transplantation (FMT) is an effective treatment for recurrent Clostridioides difficile infection (rCDI) in adults and children, but donor stool samples are currently screened for only a limited number of potential pathogens. We sought to determine whether putative procarcinogenic bacteria (enterotoxigenic Bacteroides fragilis, Fusobacterium nucleatum, and Escherichia coli harboring the colibactin toxin) could be durably transmitted from donors to patients during FMT.METHODSStool samples were collected from 11 pediatric rCDI patients and their respective FMT donors prior to FMT as well as from the patients at 2-10 weeks, 10-20 weeks, and 6 months after FMT. Bacterial virulence factors in stool DNA extracts and stool cultures were measured by quantitative PCR: Bacteroides fragilis toxin (bft), Fusobacterium adhesin A (fadA), and Escherichia coli colibactin (clbB).RESULTSFour of 11 patients demonstrated sustained acquisition of a procarcinogenic bacteria. Whole genome sequencing was performed on colony isolates from one of these donor/recipient pairs and demonstrated that clbB+ E. coli strains present in the recipient after FMT were identical to a strain present in the donor, confirming strain transmission. Conversely, 2 patients exhibited clearance of procarcinogenic bacteria following FMT from a negative donor.CONCLUSIONBoth durable transmission and clearance of procarcinogenic bacteria occurred following FMT, suggesting that additional studies on appropriate screening measures for FMT donors and the long-term consequences and/or benefits of FMT are warranted.FUNDINGCrohn's & Colitis Foundation, the Bloomberg~Kimmel Institute for Cancer Immunotherapy at Johns Hopkins University School of Medicine, the National Cancer Institute, and the Canadian Institutes of Health Research.

First Draft Genome Sequence of the Pathogenic Fungus Lomentospora prolificans (Formerly Scedosporium prolificans).

Here we describe the sequencing and assembly of the pathogenic fungus Lomentospora prolificans using a combination of short, highly accurate Illumina reads and additional coverage in very long Oxford Nanopore reads. The resulting assembly is highly contiguous, containing a total of 37,627,092 bp with over 98% of the sequence in just 26 scaffolds. Annotation identified 8896 protein-coding genes. Pulsed-field gel analysis suggests that this organism contains at least 7 and possibly 11 chromosomes, the two longest of which have sizes corresponding closely to the sizes of the longest scaffolds, at 6.6 and 5.7 Mb.

Nanopore sequencing detects structural variants in cancer.

Despite advances in sequencing, structural variants (SVs) remain difficult to reliably detect due to the short read length (<300 bp) of 2nd generation sequencing. Not only do the reads (or paired-end reads) need to straddle a breakpoint, but repetitive elements often lead to ambiguities in the alignment of short reads. We propose to use the long-reads (up to 20 kb) possible with 3rd generation sequencing, specifically nanopore sequencing on the MinION. Nanopore sequencing relies on a similar concept to a Coulter counter, reading the DNA sequence from the change in electrical current resulting from a DNA strand being forced through a nanometer-sized pore embedded in a membrane. Though nanopore sequencing currently has a relatively high mismatch rate that precludes base substitution and small frameshift mutation detection, its accuracy is sufficient for SV detection because of its long reads. In fact, long reads in some cases may improve SV detection efficiency. We have tested nanopore sequencing to detect a series of well-characterized SVs, including large deletions, inversions, and translocations that inactivate the CDKN2A/p16 and SMAD4/DPC4 tumor suppressor genes in pancreatic cancer. Using PCR amplicon mixes, we have demonstrated that nanopore sequencing can detect large deletions, translocations and inversions at dilutions as low as 1:100, with as few as 500 reads per sample. Given the speed, small footprint, and low capital cost, nanopore sequencing could become the ideal tool for the low-level detection of cancer-associated SVs needed for molecular relapse, early detection, or therapeutic monitoring.