Molecular Forecasting of Domoic Acid during a Pervasive Toxic Diatom Bloom.

In 2015, the largest recorded harmful algal bloom (HAB) occurred in the Northeast Pacific, causing nearly 100 million dollars in damages to fisheries and killing many protected marine mammals. Dominated by the toxic diatom Pseudo-nitzschia australis , this bloom produced high levels of the neurotoxin domoic acid (DA). Through molecular and transcriptional characterization of 52 near-weekly phytoplankton net-tow samples collected at a bloom hotspot in Monterey Bay, California, we identified active transcription of known DA biosynthesis ( dab ) genes from the three identified toxigenic species, including P. australis as the primary origin of toxicity. Elevated expression of silicon transporters ( sit1 ) during the bloom supports the previously hypothesized role of dissolved silica (Si) exhaustion in contributing to bloom physiology and toxicity. We find that co-expression of the dabA and sit1 genes serves as a robust predictor of DA one week in advance, potentially enabling the forecasting of DA-producing HABs. We additionally present evidence that low levels of iron could have co-limited the diatom population along with low Si. Iron limitation represents a previously unrecognized driver of both toxin production and ecological success of the low iron adapted Pseudo-nitzschia genus during the 2015 bloom, and increasing pervasiveness of iron limitation may fuel the escalating magnitude and frequency of toxic Pseudo-nitzschia blooms globally. Our results advance understanding of bloom physiology underlying toxin production, bloom prediction, and the impact of global change on toxic blooms. Significance: Pseudo-nitzschia diatoms form oceanic harmful algal blooms that threaten human health through production of the neurotoxin domoic acid (DA). DA biosynthetic gene expression is hypothesized to control DA production in the environment, yet what regulates expression of these genes is yet to be discovered. In this study, we uncovered expression of DA biosynthesis genes by multiple toxigenic Pseudo-nitzschia species during an economically impactful bloom along the North American West Coast, and identified genes that predict DA in advance of its production. We discovered that iron and silica co-limitation restrained the bloom and likely promoted toxin production. This work suggests that increasing iron limitation due to global change may play a previously unrecognized role in driving bloom frequency and toxicity.

A comprehensive and contemporary "snapshot" of ß-lactamases in carbapenem resistant Acinetobacter baumannii.

Successful treatment of Acinetobacter baumannii infections require early and appropriate antimicrobial therapy. One of the first steps in this process is understanding which ß-lactamase (bla) alleles are present and in what combinations. Thus, we performed WGS on 98 carbapenem-resistant A. baumannii (CR Ab). In most isolates, an acquired blaOXA carbapenemase was found in addition to the intrinsic blaOXA allele. The most commonly found allele was blaOXA-23 (n = 78/98). In some isolates, blaOXA-23 was found in addition to other carbapenemase alleles: blaOXA-82 (n = 12/78), blaOXA-72 (n = 2/78) and blaOXA-24/40 (n = 1/78). Surprisingly, 20% of isolates carried carbapenemases not routinely assayed for by rapid molecular diagnostic platforms, i.e., blaOXA-82 and blaOXA-172; all had ISAba1 elements. In 8 CR Ab, blaOXA-82 or blaOXA-172 was the only carbapenemase. Both blaOXA-24/40 and its variant blaOXA-72 were each found in 6/98 isolates. The most prevalent ADC variants were blaADC-30 (21%), blaADC-162 (21%), and blaADC-212 (26%). Complete combinations are reported.

Transcriptomic evidence that von Economo neurons are regionally specialized extratelencephalic-projecting excitatory neurons.

von Economo neurons (VENs) are bipolar, spindle-shaped neurons restricted to layer 5 of human frontoinsula and anterior cingulate cortex that appear to be selectively vulnerable to neuropsychiatric and neurodegenerative diseases, although little is known about other VEN cellular phenotypes. Single nucleus RNA-sequencing of frontoinsula layer 5 identifies a transcriptomically-defined cell cluster that contained VENs, but also fork cells and a subset of pyramidal neurons. Cross-species alignment of this cell cluster with a well-annotated mouse classification shows strong homology to extratelencephalic (ET) excitatory neurons that project to subcerebral targets. This cluster also shows strong homology to a putative ET cluster in human temporal cortex, but with a strikingly specific regional signature. Together these results suggest that VENs are a regionally distinctive type of ET neuron. Additionally, we describe the first patch clamp recordings of VENs from neurosurgically-resected tissue that show distinctive intrinsic membrane properties relative to neighboring pyramidal neurons.

Annotated Genome Sequence of Aspergillus tanneri NIH1004.

The annotated genome of Aspergillus tanneri, a recently discovered drug-resistant pathogen, was determined by employing the Oxford Nanopore MinION platform and the Funannotate pipeline. The genome size and the number of protein-coding genes are notably larger than those of the most common etiological agent of aspergillosis, Aspergillus fumigatus.

A Molecular Signature in Blood Reveals a Role for p53 in Regulating Malaria-Induced Inflammation.

Immunity that controls parasitemia and inflammation during Plasmodium falciparum (Pf) malaria can be acquired with repeated infections. A limited understanding of this complex immune response impedes the development of vaccines and adjunctive therapies. We conducted a prospective systems biology study of children who differed in their ability to control parasitemia and fever following Pf infection. By integrating whole-blood transcriptomics, flow-cytometric analysis, and plasma cytokine and antibody profiles, we demonstrate that a pre-infection signature of B cell enrichment, upregulation of T helper type 1 (Th1) and Th2 cell-associated pathways, including interferon responses, and p53 activation associated with control of malarial fever and coordinated with Pf-specific immunoglobulin G (IgG) and Fc receptor activation to control parasitemia. Our hypothesis-generating approach identified host molecules that may contribute to differential clinical outcomes during Pf infection. As a proof of concept, we have shown that enhanced p53 expression in monocytes attenuated Plasmodium-induced inflammation and predicted protection from fever.

Whole-blood transcriptomic signatures induced during immunization by chloroquine prophylaxis and Plasmodium falciparum sporozoites.

A highly effective vaccine that confers sterile protection to malaria is urgently needed. Immunization under chemoprophylaxis with sporozoites (CPS) consistently confers high levels of protection in the Controlled Human Malaria infection (CHMI) model. To provide a broad, unbiased assessment of the composition and kinetics of direct ex vivo human immune responses to CPS, we profiled whole-blood transcriptomes by RNA-seq before and during CPS immunization and following CHMI challenge. Differential expression of genes enriched in modules related to T cells, NK cells, protein synthesis, and mitochondrial processes were detected in fully protected individuals four weeks after the first immunization. Non-protected individuals demonstrated transcriptomic changes after the third immunization and the day of treatment, with upregulation of interferon and innate inflammatory genes and downregulation of B-cell signatures. Protected individuals demonstrated more significant interactions between blood transcription modules compared to non-protected individuals several weeks after the second and third immunizations. These data provide insight into the molecular and cellular basis of CPS-induced immune protection from P. falciparum infection.

Rapid Replacement of Acinetobacter baumannii Strains Accompanied by Changes in Lipooligosaccharide Loci and Resistance Gene Repertoire.

The population structure of health care-associated pathogens reflects patterns of diversification, selection, and dispersal over time. Empirical data detailing the long-term population dynamics of nosocomial pathogens provide information about how pathogens adapt in the face of exposure to diverse antimicrobial agents and other host and environmental pressures and can inform infection control priorities. Extensive sequencing of clinical isolates from one hospital spanning a decade and a second hospital in the Cleveland, OH, metropolitan area over a 3-year time period provided high-resolution genomic analysis of the Acinetobacter baumannii metapopulation. Genomic analysis demonstrated an almost complete replacement of the predominant strain groups with a new, genetically distinct strain group during the study period. The new group, termed clade F, differs from other global clone 2 (GC2) strains of A. baumannii in several ways, including its antibiotic resistance and lipooligosaccharide biosynthesis genes. Clade F strains are part of a large phylogenetic group with broad geographic representation. Phylogenetic analysis of single-nucleotide variants in core genome regions showed that although the Cleveland strains are phylogenetically distinct from those isolated from other locations, extensive intermixing of strains from the two hospital systems was apparent, suggesting either substantial exchange of strains or a shared, but geographically restricted, external pool from which infectious isolates were drawn. These findings document the rapid evolution of A. baumannii strains in two hospitals, with replacement of the predominant clade by a new clade with altered lipooligosaccharide loci and resistance gene repertoires.IMPORTANCE Multidrug-resistant (MDR) A. baumannii is a difficult-to-treat health care-associated pathogen. Knowing the resistance genes present in isolates causing infection aids in empirical treatment selection. Furthermore, knowledge of the genetic background can assist in tracking patterns of transmission to limit the spread of infections in hospitals. The appearance of a new genetic background in A. baumannii strains with a different set of resistance genes and cell surface structures suggests that strong selective pressures exist, even in highly MDR pathogens. Because the new strains have levels of antimicrobial resistance similar to those of the strains that were displaced, we hypothesize that other features, including host colonization and infection, may confer additional selective advantages and contribute to their increased prevalence.

Transcriptomic and morphophysiological evidence for a specialized human cortical GABAergic cell type.

We describe convergent evidence from transcriptomics, morphology, and physiology for a specialized GABAergic neuron subtype in human cortex. Using unbiased single-nucleus RNA sequencing, we identify ten GABAergic interneuron subtypes with combinatorial gene signatures in human cortical layer 1 and characterize a group of human interneurons with anatomical features never described in rodents, having large 'rosehip'-like axonal boutons and compact arborization. These rosehip cells show an immunohistochemical profile (GAD1+CCK+, CNR1-SST-CALB2-PVALB-) matching a single transcriptomically defined cell type whose specific molecular marker signature is not seen in mouse cortex. Rosehip cells in layer 1 make homotypic gap junctions, predominantly target apical dendritic shafts of layer 3 pyramidal neurons, and inhibit backpropagating pyramidal action potentials in microdomains of the dendritic tuft. These cells are therefore positioned for potent local control of distal dendritic computation in cortical pyramidal neurons.

Draft Genome Sequence and Annotation of the Apicomplexan Parasite Besnoitia besnoiti.

The apicomplexan parasite Besnoitia besnoiti is the causative agent of bovine besnoitiosis that affects livestock, particularly cattle. The definitive host of B. besnoiti is unknown and its transmission only partially understood. Here, we report the first draft genome sequence, assembly, and annotation of this parasite.

Corrected Genome Sequence of Acinetobacter baumannii Strain AB0057, an Antibiotic-Resistant Isolate from Lineage 1 of Global Clone 1.

Extensively antibiotic-resistant Acinetobacter baumannii isolate AB0057 recovered in the United States in 2004 was one of the first global clone 1 isolates to be completely sequenced. Here, the complete 4.05-Mb genome sequence (chromosome and one plasmid) has been revised using Illumina HiSeq data and targeted sequencing of PCR products.

LOCUST: a custom sequence locus typer for classifying microbial isolates.

SUMMARY: LOCUST is a custom sequence locus typer tool for classifying microbial genomes. It provides a fully automated opportunity to customize the classification of genome-wide nucleotide variant data most relevant to biological research. AVAILABILITY AND IMPLEMENTATION: Source code, demo data, and detailed documentation are freely available at http://sourceforge.net/projects/locustyper . CONTACT: lbrinkac@jcvi.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

PRODUCTION OF A PRELIMINARY QUALITY CONTROL PIPELINE FOR SINGLE NUCLEI RNA-SEQ AND ITS APPLICATION IN THE ANALYSIS OF CELL TYPE DIVERSITY OF POST-MORTEM HUMAN BRAIN NEOCORTEX.

Next generation sequencing of the RNA content of single cells or single nuclei (sc/nRNA-seq) has become a powerful approach to understand the cellular complexity and diversity of multicellular organisms and environmental ecosystems. However, the fact that the procedure begins with a relatively small amount of starting material, thereby pushing the limits of the laboratory procedures required, dictates that careful approaches for sample quality control (QC) are essential to reduce the impact of technical noise and sample bias in downstream analysis applications. Here we present a preliminary framework for sample level quality control that is based on the collection of a series of quantitative laboratory and data metrics that are used as features for the construction of QC classification models using random forest machine learning approaches. We've applied this initial framework to a dataset comprised of 2272 single nuclei RNA-seq results and determined that ~79% of samples were of high quality. Removal of the poor quality samples from downstream analysis was found to improve the cell type clustering results. In addition, this approach identified quantitative features related to the proportion of unique or duplicate reads and the proportion of reads remaining after quality trimming as useful features for pass/fail classification. The construction and use of classification models for the identification of poor quality samples provides for an objective and scalable approach to sc/nRNA-seq quality control.

Transcriptomic evidence for modulation of host inflammatory responses during febrile Plasmodium falciparum malaria.

Identifying molecular predictors and mechanisms of malaria disease is important for understanding how Plasmodium falciparum malaria is controlled. Transcriptomic studies in humans have so far been limited to retrospective analysis of blood samples from clinical cases. In this prospective, proof-of-principle study, we compared whole-blood RNA-seq profiles at pre-and post-infection time points from Malian adults who were either asymptomatic (n = 5) or febrile (n = 3) during their first seasonal PCR-positive P. falciparum infection with those from malaria-naïve Dutch adults after a single controlled human malaria infection (n = 5). Our data show a graded activation of pathways downstream of pro-inflammatory cytokines, with the highest activation in malaria-naïve Dutch individuals and significantly reduced activation in malaria-experienced Malians. Newly febrile and asymptomatic infections in Malians were statistically indistinguishable except for genes activated by pro-inflammatory cytokines. The combined data provide a molecular basis for the development of a pyrogenic threshold as individuals acquire immunity to clinical malaria.

Nuclear RNA-seq of single neurons reveals molecular signatures of activation.

Single-cell sequencing methods have emerged as powerful tools for identification of heterogeneous cell types within defined brain regions. Application of single-cell techniques to study the transcriptome of activated neurons can offer insight into molecular dynamics associated with differential neuronal responses to a given experience. Through evaluation of common whole-cell and single-nuclei RNA-sequencing (snRNA-seq) methods, here we show that snRNA-seq faithfully recapitulates transcriptional patterns associated with experience-driven induction of activity, including immediate early genes (IEGs) such as Fos, Arc and Egr1. SnRNA-seq of mouse dentate granule cells reveals large-scale changes in the activated neuronal transcriptome after brief novel environment exposure, including induction of MAPK pathway genes. In addition, we observe a continuum of activation states, revealing a pseudotemporal pattern of activation from gene expression alone. In summary, snRNA-seq of activated neurons enables the examination of gene expression beyond IEGs, allowing for novel insights into neuronal activation patterns in vivo.

Using single nuclei for RNA-seq to capture the transcriptome of postmortem neurons.

A protocol is described for sequencing the transcriptome of a cell nucleus. Nuclei are isolated from specimens and sorted by FACS, cDNA libraries are constructed and RNA-seq is performed, followed by data analysis. Some steps follow published methods (Smart-seq2 for cDNA synthesis and Nextera XT barcoded library preparation) and are not described in detail here. Previous single-cell approaches for RNA-seq from tissues include cell dissociation using protease treatment at 30 °C, which is known to alter the transcriptome. We isolate nuclei at 4 °C from tissue homogenates, which cause minimal damage. Nuclear transcriptomes can be obtained from postmortem human brain tissue stored at -80 °C, making brain archives accessible for RNA-seq from individual neurons. The method also allows investigation of biological features unique to nuclei, such as enrichment of certain transcripts and precursors of some noncoding RNAs. By following this procedure, it takes about 4 d to construct cDNA libraries that are ready for sequencing.

Genetic Analysis Using an Isogenic Mating Pair of Aspergillus fumigatus Identifies Azole Resistance Genes and Lack of MAT Locus's Role in Virulence.

Invasive aspergillosis (IA) due to Aspergillus fumigatus is a major cause of mortality in immunocompromised patients. The discovery of highly fertile strains of A. fumigatus opened the possibility to merge classical and contemporary genetics to address key questions about this pathogen. The merger involves sexual recombination, selection of desired traits, and genomics to identify any associated loci. We constructed a highly fertile isogenic pair of A. fumigatus strains with opposite mating types and used them to investigate whether mating type is associated with virulence and to find the genetic loci involved in azole resistance. The pair was made isogenic by 9 successive backcross cycles of the foundational strain AFB62 (MAT1-1) with a highly fertile (MAT1-2) progeny. Genome sequencing showed that the F9 MAT1-2 progeny was essentially identical to the AFB62. The survival curves of animals infected with either strain in three different animal models showed no significant difference, suggesting that virulence in A. fumigatus was not associated with mating type. We then employed a relatively inexpensive, yet highly powerful strategy to identify genomic loci associated with azole resistance. We used traditional in vitro drug selection accompanied by classical sexual crosses of azole-sensitive with resistant isogenic strains. The offspring were plated under varying drug concentrations and pools of resulting colonies were analyzed by whole genome sequencing. We found that variants in 5 genes contributed to azole resistance, including mutations in erg11A (cyp51A), as well as multi-drug transporters, erg25, and in HMG-CoA reductase. The results demonstrated that with minimal investment into the sequencing of three pools from a cross of interest, the variation(s) that contribute any phenotype can be identified with nucleotide resolution. This approach can be applied to multiple areas of interest in A. fumigatus or other heterothallic pathogens, especially for virulence associated traits.

NeatFreq: reference-free data reduction and coverage normalization for De Novo sequence assembly.

BACKGROUND: Deep shotgun sequencing on next generation sequencing (NGS) platforms has contributed significant amounts of data to enrich our understanding of genomes, transcriptomes, amplified single-cell genomes, and metagenomes. However, deep coverage variations in short-read data sets and high sequencing error rates of modern sequencers present new computational challenges in data interpretation, including mapping and de novo assembly. New lab techniques such as multiple displacement amplification (MDA) of single cells and sequence independent single primer amplification (SISPA) allow for sequencing of organisms that cannot be cultured, but generate highly variable coverage due to amplification biases. RESULTS: Here we introduce NeatFreq, a software tool that reduces a data set to more uniform coverage by clustering and selecting from reads binned by their median kmer frequency (RMKF) and uniqueness. Previous algorithms normalize read coverage based on RMKF, but do not include methods for the preferred selection of (1) extremely low coverage regions produced by extremely variable sequencing of random-primed products and (2) 2-sided paired-end sequences. The algorithm increases the incorporation of the most unique, lowest coverage, segments of a genome using an error-corrected data set. NeatFreq was applied to bacterial, viral plaque, and single-cell sequencing data. The algorithm showed an increase in the rate at which the most unique reads in a genome were included in the assembled consensus while also reducing the count of duplicative and erroneous contigs (strings of high confidence overlaps) in the deliverable consensus. The results obtained from conventional Overlap-Layout-Consensus (OLC) were compared to simulated multi-de Bruijn graph assembly alternatives trained for variable coverage input using sequence before and after normalization of coverage. Coverage reduction was shown to increase processing speed and reduce memory requirements when using conventional bacterial assembly algorithms. CONCLUSIONS: The normalization of deep coverage spikes, which would otherwise inhibit consensus resolution, enables High Throughput Sequencing (HTS) assembly projects to consistently run to completion with existing assembly software. The NeatFreq software package is free, open source and available at https://github.com/bioh4x/NeatFreq .

PGP: parallel prokaryotic proteogenomics pipeline for MPI clusters, high-throughput batch clusters and multicore workstations.

SUMMARY: We present the first public release of our proteogenomic annotation pipeline. We have previously used our original unreleased implementation to improve the annotation of 46 diverse prokaryotic genomes by discovering novel genes, post-translational modifications and correcting the erroneous annotations by analyzing proteomic mass-spectrometry data. This public version has been redesigned to run in a wide range of parallel Linux computing environments and provided with the automated configuration, build and testing facilities for easy deployment and portability. AVAILABILITY AND IMPLEMENTATION: Source code is freely available from https://bitbucket.org/andreyto/proteogenomics under GPL license. It is implemented in Python and C++. It bundles the Makeflow engine to execute the workflows. CONTACT: atovtchi@jcvi.org.