Single-molecule nanopore sequencing reveals extreme target copy number heterogeneity in arylomycin-resistant mutants.

Tandem gene amplification is a frequent and dynamic source of antibiotic resistance in bacteria. Ongoing expansions and contractions of repeat arrays during population growth are expected to manifest as cell-to-cell differences in copy number (CN). As a result, a clonal bacterial culture could comprise subpopulations of cells with different levels of antibiotic sensitivity that result from variable gene dosage. Despite the high potential for misclassification of heterogenous cell populations as either antibiotic-susceptible or fully resistant in clinical settings, and the concomitant risk of inappropriate treatment, CN distribution among cells has defied analysis. Here, we use the MinION single-molecule nanopore sequencer to uncover CN heterogeneity in clonal populations of Escherichia coli and Acinetobacter baumannii grown from single cells isolated while selecting for resistance to an optimized arylomycin, a member of a recently discovered class of Gram-negative antibiotic. We found that gene amplification of the arylomycin target, bacterial type I signal peptidase LepB, is a mechanism of unstable arylomycin resistance and demonstrate in E. coli that amplification instability is independent of RecA. This instability drives the emergence of a nonuniform distribution of lepB CN among cells with a range of 1 to at least 50 copies of lepB identified in a single clonal population. In sum, this remarkable heterogeneity, and the evolutionary plasticity it fuels, illustrates how gene amplification can enable bacterial populations to respond rapidly to novel antibiotics. This study establishes a rationale for further nanopore-sequencing studies of heterogeneous cell populations to uncover CN variability at single-molecule resolution.

Recovery 101: A Community Psychoeducation Program on Mental Health Recovery and Peer Support Services Created by an Academic Psychiatrist and Delivered by Peers.

OBJECTIVE: Despite the ongoing growing interest and established benefits of mental health recovery and peer support services, little is known about the effects of mental health recovery and peer support services specifically in community education programs. Seeking to further expand this scope, this article details the development, implementation, and evaluation of a recovery and peer services community psychoeducation program created by an academic psychiatrist and delivered by trained mental health peers. METHODS: A community psychoeducation intervention on recovery and peer services was provided in 2 different series: 5-session and 1-session series. Participants in the intervention were recruited from 2 local nonprofit recovery-oriented organizations and a public mental health service agency. Surveys were administered before and after each presentation session and analyzed both quantitatively and qualitatively. RESULTS: Quantitative survey items covering knowledge, awareness, interest, and likelihood of involvement in peer services for recovery all showed statistically significant gains from before to after the sessions. The qualitative data provided information about the participants' interests in the program in terms of knowledge of recovery, self-improvement, advocacy, resources, support, and nonspecific positive comments. CONCLUSIONS: An essential contribution of the program was to introduce participants to recovery and peer services. It also provided them with foundational knowledge that placed participants in a position to begin to practice the concepts of recovery. According to the participants' comments in the surveys, the program imbued a sense of hope and motivation; support and fellowship; and, importantly, role modeling and wisdom from the presenters.

Comparison of hospital course and outcome in hospital-acquired vs pre-existing delirium.

BACKGROUND: Delirium is a major source of morbidity in the inpatient hospital setting. This study examined differences between patients with delirium present prior to hospital admission and those with hospitalacquired delirium in several health outcomes. METHODS: A total of 12,529 patients on 2 inpatient units were included in this retrospective cohort study. Outcomes were assessed using chart review. Other variables were compared across groups and included in multivariate models predicting discharge location within the hospitalacquired delirium group. RESULTS: Of 709 patients with delirium, 83% had pre-admission prevalent and 17% had post-admission incident delirium. Compared with patients with preexisting delirium, patients with hospital-acquired delirium had greater hospital durations and mortality and were more likely to receive ICU care, more likely to receive multiple classes of medications, and less likely to be discharged home without home health services. Multivariate analysis in the hospital-acquired delirium group found that several variables independently predicted discharge location. CONCLUSIONS: Patients with hospital-acquired delirium had worse hospital outcomes and a more complicated hospital course than those with preexisting delirium. Administration of various medications, several demographic variables, and some hospital-related variables were independently associated with worse outcomes within the hospital-acquired delirium group. These results demonstrate that patients with hospitalacquired delirium are a vulnerable subgroup deserving special attention.

Recurrent R-spondin fusions in colon cancer.

Identifying and understanding changes in cancer genomes is essential for the development of targeted therapeutics. Here we analyse systematically more than 70 pairs of primary human colon tumours by applying next-generation sequencing to characterize their exomes, transcriptomes and copy-number alterations. We have identified 36,303 protein-altering somatic changes that include several new recurrent mutations in the Wnt pathway gene TCF7L2, chromatin-remodelling genes such as TET2 and TET3 and receptor tyrosine kinases including ERBB3. Our analysis for significantly mutated cancer genes identified 23 candidates, including the cell cycle checkpoint kinase ATM. Copy-number and RNA-seq data analysis identified amplifications and corresponding overexpression of IGF2 in a subset of colon tumours. Furthermore, using RNA-seq data we identified multiple fusion transcripts including recurrent gene fusions involving R-spondin family members RSPO2 and RSPO3 that together occur in 10% of colon tumours. The RSPO fusions were mutually exclusive with APC mutations, indicating that they probably have a role in the activation of Wnt signalling and tumorigenesis. Consistent with this we show that the RSPO fusion proteins were capable of potentiating Wnt signalling. The R-spondin gene fusions and several other gene mutations identified in this study provide new potential opportunities for therapeutic intervention in colon cancer.

Massively parallel nanowell-based single-cell gene expression profiling.

BACKGROUND: Technological advances have enabled transcriptome characterization of cell types at the single-cell level providing new biological insights. New methods that enable simple yet high-throughput single-cell expression profiling are highly desirable. RESULTS: Here we report a novel nanowell-based single-cell RNA sequencing system, ICELL8, which enables processing of thousands of cells per sample. The system employs a 5,184-nanowell-containing microchip to capture ~1,300 single cells and process them. Each nanowell contains preprinted oligonucleotides encoding poly-d(T), a unique well barcode, and a unique molecular identifier. The ICELL8 system uses imaging software to identify nanowells containing viable single cells and only wells with single cells are processed into sequencing libraries. Here, we report the performance and utility of ICELL8 using samples of increasing complexity from cultured cells to mouse solid tissue samples. Our assessment of the system to discriminate between mixed human and mouse cells showed that ICELL8 has a low cell multiplet rate (< 3%) and low cross-cell contamination. We characterized single-cell transcriptomes of more than a thousand cultured human and mouse cells as well as 468 mouse pancreatic islets cells. We were able to identify distinct cell types in pancreatic islets, including alpha, beta, delta and gamma cells. CONCLUSIONS: Overall, ICELL8 provides efficient and cost-effective single-cell expression profiling of thousands of cells, allowing researchers to decipher single-cell transcriptomes within complex biological samples.

Genome and transcriptome sequencing of lung cancers reveal diverse mutational and splicing events.

Lung cancer is a highly heterogeneous disease in terms of both underlying genetic lesions and response to therapeutic treatments. We performed deep whole-genome sequencing and transcriptome sequencing on 19 lung cancer cell lines and three lung tumor/normal pairs. Overall, our data show that cell line models exhibit similar mutation spectra to human tumor samples. Smoker and never-smoker cancer samples exhibit distinguishable patterns of mutations. A number of epigenetic regulators, including KDM6A, ASH1L, SMARCA4, and ATAD2, are frequently altered by mutations or copy number changes. A systematic survey of splice-site mutations identified 106 splice site mutations associated with cancer specific aberrant splicing, including mutations in several known cancer-related genes. RAC1b, an isoform of the RAC1 GTPase that includes one additional exon, was found to be preferentially up-regulated in lung cancer. We further show that its expression is significantly associated with sensitivity to a MAP2K (MEK) inhibitor PD-0325901. Taken together, these data present a comprehensive genomic landscape of a large number of lung cancer samples and further demonstrate that cancer-specific alternative splicing is a widespread phenomenon that has potential utility as therapeutic biomarkers. The detailed characterizations of the lung cancer cell lines also provide genomic context to the vast amount of experimental data gathered for these lines over the decades, and represent highly valuable resources for cancer biology.

Disruption of PH-kinase domain interactions leads to oncogenic activation of AKT in human cancers.

The protein kinase v-akt murine thymoma viral oncogene homolog (AKT), a key regulator of cell survival and proliferation, is frequently hyperactivated in human cancers. Intramolecular pleckstrin homology (PH) domain-kinase domain (KD) interactions are important in maintaining AKT in an inactive state. AKT activation proceeds after a conformational change that dislodges the PH from the KD. To understand these autoinhibitory interactions, we generated mutations at the PH-KD interface and found that most of them lead to constitutive activation of AKT. Such mutations are likely another mechanism by which activation may occur in human cancers and other diseases. In support of this likelihood, we found somatic mutations in AKT1 at the PH-KD interface that have not been previously described in human cancers. Furthermore, we show that the AKT1 somatic mutants are constitutively active, leading to oncogenic signaling. Additionally, our studies show that the AKT1 mutants are not effectively inhibited by allosteric AKT inhibitors, consistent with the requirement for an intact PH-KD interface for allosteric inhibition. These results have important implications for therapeutic intervention in patients with AKT mutations at the PH-KD interface.

Argonaute3-SF3B3 complex controls pre-mRNA splicing to restrain type 2 immunity.

Argonaute (AGO) proteins execute microRNA (miRNA)-mediated gene silencing. However, it is unclear whether all 4 mammalian AGO proteins (AGO1, AGO2, AGO3, and AGO4) are required for miRNA activity. We generate Ago1, Ago3, and Ago4-deficient mice (Ago134Δ) and find AGO1/3/4 to be redundant for miRNA biogenesis, homeostasis, or function, a role that is carried out by AGO2. Instead, AGO1/3/4 regulate the expansion of type 2 immunity via precursor mRNA splicing in CD4+ T helper (Th) lymphocytes. Gain- and loss-of-function experiments demonstrate that nuclear AGO3 interacts directly with SF3B3, a component of the U2 spliceosome complex, to aid global mRNA splicing, and in particular the isoforms of the gene Nisch, resulting in a dysregulated Nisch isoform ratio. This work uncouples AGO1, AGO3, and AGO4 from miRNA-mediated RNA interference, identifies an AGO3:SF3B3 complex in the nucleus, and reveals a mechanism by which AGO proteins regulate inflammatory diseases.

Accelerated identification of disease-causing variants with ultra-rapid nanopore genome sequencing.

Whole-genome sequencing (WGS) can identify variants that cause genetic disease, but the time required for sequencing and analysis has been a barrier to its use in acutely ill patients. In the present study, we develop an approach for ultra-rapid nanopore WGS that combines an optimized sample preparation protocol, distributing sequencing over 48 flow cells, near real-time base calling and alignment, accelerated variant calling and fast variant filtration for efficient manual review. Application to two example clinical cases identified a candidate variant in <8 h from sample preparation to variant identification. We show that this framework provides accurate variant calls and efficient prioritization, and accelerates diagnostic clinical genome sequencing twofold compared with previous approaches.

Accurate assembly of the olive baboon (Papio anubis) genome using long-read and Hi-C data.

BACKGROUND: Baboons are a widely used nonhuman primate model for biomedical, evolutionary, and basic genetics research. Despite this importance, the genomic resources for baboons are limited. In particular, the current baboon reference genome Panu_3.0 is a highly fragmented, reference-guided (i.e., not fully de novo) assembly, and its poor quality inhibits our ability to conduct downstream genomic analyses. FINDINGS: Here we present a de novo genome assembly of the olive baboon (Papio anubis) that uses data from several recently developed single-molecule technologies. Our assembly, Panubis1.0, has an N50 contig size of ∼1.46 Mb (as opposed to 139 kb for Panu_3.0) and has single scaffolds that span each of the 20 autosomes and the X chromosome. CONCLUSIONS: We highlight multiple lines of evidence (including Bionano Genomics data, pedigree linkage information, and linkage disequilibrium data) suggesting that there are several large assembly errors in Panu_3.0, which have been corrected in Panubis1.0.

The Indian cobra reference genome and transcriptome enables comprehensive identification of venom toxins.

Snakebite envenoming is a serious and neglected tropical disease that kills ~100,000 people annually. High-quality, genome-enabled comprehensive characterization of toxin genes will facilitate development of effective humanized recombinant antivenom. We report a de novo near-chromosomal genome assembly of Naja naja, the Indian cobra, a highly venomous, medically important snake. Our assembly has a scaffold N50 of 223.35 Mb, with 19 scaffolds containing 95% of the genome. Of the 23,248 predicted protein-coding genes, 12,346 venom-gland-expressed genes constitute the 'venom-ome' and this included 139 genes from 33 toxin families. Among the 139 toxin genes were 19 'venom-ome-specific toxins' (VSTs) that showed venom-gland-specific expression, and these probably encode the minimal core venom effector proteins. Synthetic venom reconstituted through recombinant VST expression will aid in the rapid development of safe and effective synthetic antivenom. Additionally, our genome could serve as a reference for snake genomes, support evolutionary studies and enable venom-driven drug discovery.

Barcoded sequencing workflow for high throughput digitization of hybridoma antibody variable domain sequences.

Since the invention of Hybridoma technology by Milstein and Köhler in 1975, its application has greatly advanced the antibody discovery process. The technology enables both functional screening and long-term archival of the immortalized monoclonal antibody producing B cells. Despite the dependable cryopreservation technology for hybridoma cells, practicality of long-term storage has been outpaced by recent progress in robotics and automations, which enables routine identification of thousands of antigen specific hybridoma clones. Such throughput increase imposes two nascent challenges in the antibody discovery process, namely limited cryopreservation storage space and limited throughput in conventional antibody sequencing. We herein provide a barcoded sequencing workflow that utilizes next generation sequencing to expand the conventional sequencing capacity. Accompanied with the bioinformatics tools we describe, the barcoded sequencing workflow robustly reports unambiguous antibody sequences as confirmed with Sanger sequencing controls. In complement with the commonly accessible recombinant DNA technology, the barcoded sequencing workflow allows for high throughput digitization of the antibody sequences and provides an effective solution to the limitations imposed by physical storage and sequencing capacity.

Author Correction: Comprehensive analysis of single molecule sequencing-derived complete genome and whole transcriptome of Hyposidra talaca nuclear polyhedrosis virus.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

Comprehensive analysis of single molecule sequencing-derived complete genome and whole transcriptome of Hyposidra talaca nuclear polyhedrosis virus.

We sequenced the Hyposidra talaca NPV (HytaNPV) double stranded circular DNA genome using PacBio single molecule sequencing technology. We found that the HytaNPV genome is 139,089 bp long with a GC content of 39.6%. It encodes 141 open reading frames (ORFs) including the 37 baculovirus core genes, 25 genes conserved among lepidopteran baculoviruses, 72 genes known in baculovirus, and 7 genes unique to the HytaNPV genome. It is a group II alphabaculovirus that codes for the F protein and lacks the gp64 gene found in group I alphabaculovirus viruses. Using RNA-seq, we confirmed the expression of the ORFs identified in the HytaNPV genome. Phylogenetic analysis showed HytaNPV to be closest to BusuNPV, SujuNPV and EcobNPV that infect other tea pests, Buzura suppressaria, Sucra jujuba, and Ectropis oblique, respectively. We identified repeat elements and a conserved non-coding baculovirus element in the genome. Analysis of the putative promoter sequences identified motif consistent with the temporal expression of the genes observed in the RNA-seq data.

The kinase TPL2 activates ERK and p38 signaling to promote neutrophilic inflammation.

Tumor progression locus 2 (TPL2; also known as MAP3K8) is a mitogen-activated protein kinase (MAPK) kinase kinase (MAP3K) that phosphorylates the MAPK kinases MEK1 and MEK2 (MEK1/2), which, in turn, activate the MAPKs extracellular signal-regulated kinase 1 (ERK1) and ERK2 (ERK1/2) in macrophages stimulated through the interleukin-1 receptor (IL-1R), Toll-like receptors (TLRs), or the tumor necrosis factor receptor (TNFR). We describe a conserved and critical role for TPL2 in mediating the effector functions of neutrophils through the activation of the p38 MAPK signaling pathway. Gene expression profiling and functional studies of neutrophils and monocytes revealed a MEK1/2-independent branch point downstream of TPL2 in neutrophils. Biochemical analyses identified the MAPK kinases MEK3 and MEK6 and the MAPKs p38α and p38δ as downstream effectors of TPL2 in these cells. Genetic ablation of the catalytic activity of TPL2 or therapeutic intervention with a TPL2-specific inhibitor reduced the production of inflammatory mediators by neutrophils in response to stimulation with the TLR4 agonist lipopolysaccharide (LPS) in vitro, as well as in rodent models of inflammatory disease. Together, these data suggest that TPL2 is a drug target that activates not only MEK1/2-dependent but also MEK3/6-dependent signaling to promote inflammatory responses.

Comprehensive genomic analysis of malignant pleural mesothelioma identifies recurrent mutations, gene fusions and splicing alterations.

We analyzed transcriptomes (n = 211), whole exomes (n = 99) and targeted exomes (n = 103) from 216 malignant pleural mesothelioma (MPM) tumors. Using RNA-seq data, we identified four distinct molecular subtypes: sarcomatoid, epithelioid, biphasic-epithelioid (biphasic-E) and biphasic-sarcomatoid (biphasic-S). Through exome analysis, we found BAP1, NF2, TP53, SETD2, DDX3X, ULK2, RYR2, CFAP45, SETDB1 and DDX51 to be significantly mutated (q-score ≥ 0.8) in MPMs. We identified recurrent mutations in several genes, including SF3B1 (∼2%; 4/216) and TRAF7 (∼2%; 5/216). SF3B1-mutant samples showed a splicing profile distinct from that of wild-type tumors. TRAF7 alterations occurred primarily in the WD40 domain and were, except in one case, mutually exclusive with NF2 alterations. We found recurrent gene fusions and splice alterations to be frequent mechanisms for inactivation of NF2, BAP1 and SETD2. Through integrated analyses, we identified alterations in Hippo, mTOR, histone methylation, RNA helicase and p53 signaling pathways in MPMs.

Spectrum of diverse genomic alterations define non-clear cell renal carcinoma subtypes.

To further understand the molecular distinctions between kidney cancer subtypes, we analyzed exome, transcriptome and copy number alteration data from 167 primary human tumors that included renal oncocytomas and non-clear cell renal cell carcinomas (nccRCCs), consisting of papillary (pRCC), chromophobe (chRCC) and translocation (tRCC) subtypes. We identified ten significantly mutated genes in pRCC, including MET, NF2, SLC5A3, PNKD and CPQ. MET mutations occurred in 15% (10/65) of pRCC samples and included previously unreported recurrent activating mutations. In chRCC, we found TP53, PTEN, FAAH2, PDHB, PDXDC1 and ZNF765 to be significantly mutated. Gene expression analysis identified a five-gene set that enabled the molecular classification of chRCC, renal oncocytoma and pRCC. Using RNA sequencing, we identified previously unreported gene fusions, including ACTG1-MITF fusion. Ectopic expression of the ACTG1-MITF fusion led to cellular transformation and induced the expression of downstream target genes. Finally, we observed upregulation of the anti-apoptotic factor BIRC7 in MiTF-high RCC tumors, suggesting a potential therapeutic role for BIRC7 inhibitors.

Comparison of sequencing platforms for single nucleotide variant calls in a human sample.

Next-generation sequencings platforms coupled with advanced bioinformatic tools enable re-sequencing of the human genome at high-speed and large cost savings. We compare sequencing platforms from Roche/454(GS FLX), Illumina/HiSeq (HiSeq 2000), and Life Technologies/SOLiD (SOLiD 3 ECC) for their ability to identify single nucleotide substitutions in whole genome sequences from the same human sample. We report on significant GC-related bias observed in the data sequenced on Illumina and SOLiD platforms. The differences in the variant calls were investigated with regards to coverage, and sequencing error. Some of the variants called by only one or two of the platforms were experimentally tested using mass spectrometry; a method that is independent of DNA sequencing. We establish several causes why variants remained unreported, specific to each platform. We report the indel called using the three sequencing technologies and from the obtained results we conclude that sequencing human genomes with more than a single platform and multiple libraries is beneficial when high level of accuracy is required.