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.

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.

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.

RNA-sequencing from single nuclei.

It has recently been established that synthesis of double-stranded cDNA can be done from a single cell for use in DNA sequencing. Global gene expression can be quantified from the number of reads mapping to each gene, and mutations and mRNA splicing variants determined from the sequence reads. Here we demonstrate that this method of transcriptomic analysis can be done using the extremely low levels of mRNA in a single nucleus, isolated from a mouse neural progenitor cell line and from dissected hippocampal tissue. This method is characterized by excellent coverage and technical reproducibility. On average, more than 16,000 of the 24,057 mouse protein-coding genes were detected from single nuclei, and the amount of gene-expression variation was similar when measured between single nuclei and single cells. Several major advantages of the method exist: first, nuclei, compared with whole cells, have the advantage of being easily isolated from complex tissues and organs, such as those in the CNS. Second, the method can be widely applied to eukaryotic species, including those of different kingdoms. The method also provides insight into regulatory mechanisms specific to the nucleus. Finally, the method enables dissection of regulatory events at the single-cell level; pooling of 10 nuclei or 10 cells obscures some of the variability measured in transcript levels, implying that single nuclei and cells will be extremely useful in revealing the physiological state and interconnectedness of gene regulation in a manner that avoids the masking inherent to conventional transcriptomics using bulk cells or tissues.

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.

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 .

Sequence analysis of a complete 1.66 Mb Prochlorococcus marinus MED4 genome cloned in yeast.

Marine cyanobacteria of the genus Prochlorococcus represent numerically dominant photoautotrophs residing throughout the euphotic zones in the open oceans and are major contributors to the global carbon cycle. Prochlorococcus has remained a genetically intractable bacterium due to slow growth rates and low transformation efficiencies using standard techniques. Our recent successes in cloning and genetically engineering the AT-rich, 1.1 Mb Mycoplasma mycoides genome in yeast encouraged us to explore similar methods with Prochlorococcus. Prochlorococcus MED4 has an AT-rich genome, with a GC content of 30.8%, similar to that of Saccharomyces cerevisiae (38%), and contains abundant yeast replication origin consensus sites (ACS) evenly distributed around its 1.66 Mb genome. Unlike Mycoplasma cells, which use the UGA codon for tryptophane, Prochlorococcus uses the standard genetic code. Despite this, we observed no toxic effects of several partial and 15 whole Prochlorococcus MED4 genome clones in S. cerevisiae. Sequencing of a Prochlorococcus genome purified from yeast identified 14 single base pair missense mutations, one frameshift, one single base substitution to a stop codon and one dinucleotide transversion compared to the donor genomic DNA. We thus provide evidence of transformation, replication and maintenance of this 1.66 Mb intact bacterial genome in S. cerevisiae.

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.

Experimental annotation of post-translational features and translated coding regions in the pathogen Salmonella Typhimurium.

BACKGROUND: Complete and accurate genome annotation is crucial for comprehensive and systematic studies of biological systems. However, determining protein-coding genes for most new genomes is almost completely performed by inference using computational predictions with significant documented error rates (> 15%). Furthermore, gene prediction programs provide no information on biologically important post-translational processing events critical for protein function. RESULTS: We experimentally annotated the bacterial pathogen Salmonella Typhimurium 14028, using "shotgun" proteomics to accurately uncover the translational landscape and post-translational features. The data provide protein-level experimental validation for approximately half of the predicted protein-coding genes in Salmonella and suggest revisions to several genes that appear to have incorrectly assigned translational start sites, including a potential novel alternate start codon. Additionally, we uncovered 12 non-annotated genes missed by gene prediction programs, as well as evidence suggesting a role for one of these novel ORFs in Salmonella pathogenesis. We also characterized post-translational features in the Salmonella genome, including chemical modifications and proteolytic cleavages. We find that bacteria have a much larger and more complex repertoire of chemical modifications than previously thought including several novel modifications. Our in vivo proteolysis data identified more than 130 signal peptide and N-terminal methionine cleavage events critical for protein function. CONCLUSION: This work highlights several ways in which application of proteomics data can improve the quality of genome annotations to facilitate novel biological insights and provides a comprehensive proteome map of Salmonella as a resource for systems analysis.

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.

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.

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.

A bioinformatic filter for improved base-call accuracy and polymorphism detection using the Affymetrix GeneChip whole-genome resequencing platform.

DNA resequencing arrays enable rapid acquisition of high-quality sequence data. This technology represents a promising platform for rapid high-resolution genotyping of microorganisms. Traditional array-based resequencing methods have relied on the use of specific PCR-amplified fragments from the query samples as hybridization targets. While this specificity in the target DNA population reduces the potential for artifacts caused by cross-hybridization, the subsampling of the query genome limits the sequence coverage that can be obtained and therefore reduces the technique's resolution as a genotyping method. We have developed and validated an Affymetrix Inc. GeneChip(R) array-based, whole-genome resequencing platform for Francisella tularensis, the causative agent of tularemia. A set of bioinformatic filters that targeted systematic base-calling errors caused by cross-hybridization between the whole-genome sample and the array probes and by deletions in the sample DNA relative to the chip reference sequence were developed. Our approach eliminated 91% of the false-positive single-nucleotide polymorphism calls identified in the SCHU S4 query sample, at the cost of 10.7% of the true positives, yielding a total base-calling accuracy of 99.992%.

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.

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.

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.