Spatial Sequencing: A Perspective.

Fluorescent spatial sequencing brings next-generation sequencing into a new realm capable of identifying nucleic acids in the cell's natural environment. For the first time, scientists are able to multiplex the assignment of specific locations to hundreds of transcriptional targets and lay the foundation for understanding how genetic changes control the fate of each cell within the tissue microenvironment. In this perspective, we discuss the capabilities of fluorescent spatial sequencing in the context of other spatial imaging technologies and describe how these new technologies offer a data-rich, multiomic solution to many research applications. Fluorescent spatial sequencing has opened options for exploring many fundamental questions in biology, helping us gain a better understanding of cell and tissue development and disease progression.

H3K27me3-mediated silencing of structural genes is required for zebrafish heart regeneration.

Deciphering the genetic and epigenetic regulation of cardiomyocyte proliferation in organisms that are capable of robust cardiac renewal, such as zebrafish, represents an attractive inroad towards regenerating the human heart. Using integrated high-throughput transcriptional and chromatin analyses, we have identified a strong association between H3K27me3 deposition and reduced sarcomere and cytoskeletal gene expression in proliferative cardiomyocytes following cardiac injury in zebrafish. To move beyond an association, we generated an inducible transgenic strain expressing a mutant version of histone 3, H3.3K27M, that inhibits H3K27me3 catalysis in cardiomyocytes during the regenerative window. Hearts comprising H3.3K27M-expressing cardiomyocytes fail to regenerate, with wound edge cells showing heightened expression of structural genes and prominent sarcomeres. Although cell cycle re-entry was unperturbed, cytokinesis and wound invasion were significantly compromised. Collectively, our study identifies H3K27me3-mediated silencing of structural genes as requisite for zebrafish heart regeneration and suggests that repression of similar structural components in the border zone of an infarcted human heart might improve its regenerative capacity.

Deep learning using tumor HLA peptide mass spectrometry datasets improves neoantigen identification.

Neoantigens, which are expressed on tumor cells, are one of the main targets of an effective antitumor T-cell response. Cancer immunotherapies to target neoantigens are of growing interest and are in early human trials, but methods to identify neoantigens either require invasive or difficult-to-obtain clinical specimens, require the screening of hundreds to thousands of synthetic peptides or tandem minigenes, or are only relevant to specific human leukocyte antigen (HLA) alleles. We apply deep learning to a large (N = 74 patients) HLA peptide and genomic dataset from various human tumors to create a computational model of antigen presentation for neoantigen prediction. We show that our model, named EDGE, increases the positive predictive value of HLA antigen prediction by up to ninefold. We apply EDGE to enable identification of neoantigens and neoantigen-reactive T cells using routine clinical specimens and small numbers of synthetic peptides for most common HLA alleles. EDGE could enable an improved ability to develop neoantigen-targeted immunotherapies for cancer patients.

Author Correction: Comprehensive comparative analysis of 5'-end RNA-sequencing methods.

The original version of this paper contained an incorrect primer sequence. In the Methods subsection "Rampage libraries," the text for modification 3 stated that the reverse primer used for library indexing was 5'-CAAGCAGAAGACGGCATACGAGATXXXXXXXXGTGACTGGAGT-3'. The correct sequence of the oligonucleotide used is 5'-CAAGCAGAAGACGGCATACGAGATXXXXXXXXGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT-3'. This error has been corrected in the PDF and HTML versions of the paper.

Comprehensive comparative analysis of 5'-end RNA-sequencing methods.

Specialized RNA-seq methods are required to identify the 5' ends of transcripts, which are critical for studies of gene regulation, but these methods have not been systematically benchmarked. We directly compared six such methods, including the performance of five methods on a single human cellular RNA sample and a new spike-in RNA assay that helps circumvent challenges resulting from uncertainties in annotation and RNA processing. We found that the 'cap analysis of gene expression' (CAGE) method performed best for mRNA and that most of its unannotated peaks were supported by evidence from other genomic methods. We applied CAGE to eight brain-related samples and determined sample-specific transcription start site (TSS) usage, as well as a transcriptome-wide shift in TSS usage between fetal and adult brain.

Evidence of Ebola Virus Replication and High Concentration in Semen of a Patient During Recovery.

In one patient over time, we found that concentration of Ebola virus RNA in semen during recovery is remarkably higher than blood at peak illness. Virus in semen is replication-competent with no change in viral genome over time. Presence of sense RNA suggests replication in cells present in semen.

Systematic comparison of monoclonal versus polyclonal antibodies for mapping histone modifications by ChIP-seq.

BACKGROUND: The robustness of ChIP-seq datasets is highly dependent upon the antibodies used. Currently, polyclonal antibodies are the standard despite several limitations: They are non-renewable, vary in performance between lots and need to be validated with each new lot. In contrast, monoclonal antibody lots are renewable and provide consistent performance. To increase ChIP-seq standardization, we investigated whether monoclonal antibodies could replace polyclonal antibodies. We compared monoclonal antibodies that target five key histone modifications (H3K4me1, H3K4me3, H3K9me3, H3K27ac and H3K27me3) to their polyclonal counterparts in both human and mouse cells. RESULTS: Overall performance was highly similar for four monoclonal/polyclonal pairs, including when we used two distinct lots of the same monoclonal antibody. In contrast, the binding patterns for H3K27ac differed substantially between polyclonal and monoclonal antibodies. However, this was most likely due to the distinct immunogen used rather than the clonality of the antibody. CONCLUSIONS: Altogether, we found that monoclonal antibodies as a class perform equivalently to polyclonal antibodies for the detection of histone post-translational modifications in both human and mouse. Accordingly, we recommend the use of monoclonal antibodies in ChIP-seq experiments.

Simultaneous generation of many RNA-seq libraries in a single reaction.

Although RNA-seq is a powerful tool, the considerable time and cost associated with library construction has limited its utilization for various applications. RNAtag-Seq, an approach to generate multiple RNA-seq libraries in a single reaction, lowers time and cost per sample, and it produces data on prokaryotic and eukaryotic samples that are comparable to those generated by traditional strand-specific RNA-seq approaches.

Enhanced methods for unbiased deep sequencing of Lassa and Ebola RNA viruses from clinical and biological samples.

We have developed a robust RNA sequencing method for generating complete de novo assemblies with intra-host variant calls of Lassa and Ebola virus genomes in clinical and biological samples. Our method uses targeted RNase H-based digestion to remove contaminating poly(rA) carrier and ribosomal RNA. This depletion step improves both the quality of data and quantity of informative reads in unbiased total RNA sequencing libraries. We have also developed a hybrid-selection protocol to further enrich the viral content of sequencing libraries. These protocols have enabled rapid deep sequencing of both Lassa and Ebola virus and are broadly applicable to other viral genomics studies.

Functional optimization of gene clusters by combinatorial design and assembly.

Large microbial gene clusters encode useful functions, including energy utilization and natural product biosynthesis, but genetic manipulation of such systems is slow, difficult and complicated by complex regulation. We exploit the modularity of a refactored Klebsiella oxytoca nitrogen fixation (nif) gene cluster (16 genes, 103 parts) to build genetic permutations that could not be achieved by starting from the wild-type cluster. Constraint-based combinatorial design and DNA assembly are used to build libraries of radically different cluster architectures by varying part choice, gene order, gene orientation and operon occupancy. We construct 84 variants of the nifUSVWZM operon, 145 variants of the nifHDKY operon, 155 variants of the nifHDKYENJ operon and 122 variants of the complete 16-gene pathway. The performance and behavior of these variants are characterized by nitrogenase assay and strand-specific RNA sequencing (RNA-seq), and the results are incorporated into subsequent design cycles. We have produced a fully synthetic cluster that recovers 57% of wild-type activity. Our approach allows the performance of genetic parts to be quantified simultaneously in hundreds of genetic contexts. This parallelized design-build-test-learn cycle, which can access previously unattainable regions of genetic space, should provide a useful, fast tool for genetic optimization and hypothesis testing.

Comparative analysis of RNA sequencing methods for degraded or low-input samples.

RNA-seq is an effective method for studying the transcriptome, but it can be difficult to apply to scarce or degraded RNA from fixed clinical samples, rare cell populations or cadavers. Recent studies have proposed several methods for RNA-seq of low-quality and/or low-quantity samples, but the relative merits of these methods have not been systematically analyzed. Here we compare five such methods using metrics relevant to transcriptome annotation, transcript discovery and gene expression. Using a single human RNA sample, we constructed and sequenced ten libraries with these methods and compared them against two control libraries. We found that the RNase H method performed best for chemically fragmented, low-quality RNA, and we confirmed this through analysis of actual degraded samples. RNase H can even effectively replace oligo(dT)-based methods for standard RNA-seq. SMART and NuGEN had distinct strengths for measuring low-quantity RNA. Our analysis allows biologists to select the most suitable methods and provides a benchmark for future method development.

Scotty: a web tool for designing RNA-Seq experiments to measure differential gene expression.

MOTIVATION: A common question arises at the beginning of every experiment where RNA-Seq is used to detect differential gene expression between two conditions: How many reads should we sequence? RESULTS: Scotty is an interactive web-based application that assists biologists to design an experiment with an appropriate sample size and read depth to satisfy the user-defined experimental objectives. This design can be based on data available from either pilot samples or publicly available datasets. AVAILABILITY: Scotty can be freely accessed on the web at http://euler.bc.edu/marthlab/scotty/scotty.php

Targeted proteomic dissection of Toxoplasma cytoskeleton sub-compartments using MORN1.

The basal complex in Toxoplasma functions as the contractile ring in the cell division process. Basal complex contraction tapers the daughter cytoskeleton toward the basal end and is required for daughter segregation. We have previously shown that the protein MORN1 is essential for basal complex assembly and likely acts as a scaffolding protein. To further our understanding of the basal complex, we combined subcellular fractionation with an affinity purification of the MORN1 complex and identified its protein composition. We identified two new components of the basal complex, one of which uniquely associated with the basal complex in mature parasites, the first of its kind. In addition, we identified several other novel cytoskeleton proteins with different spatiotemporal dynamics throughout cell division. Since many of these proteins are unique to Apicomplexa this study significantly contributes to the annotation of their unique cytoskeleton. Furthermore, we show that G-actin binding protein TgCAP is localized at the apical cap region in intracellular parasites, but quickly redistributes to a cytoplasmic localization pattern upon egress. © 2012 Wiley Periodicals, Inc.

Expression divergence measured by transcriptome sequencing of four yeast species.

BACKGROUND: The evolution of gene expression is a challenging problem in evolutionary biology, for which accurate, well-calibrated measurements and methods are crucial. RESULTS: We quantified gene expression with whole-transcriptome sequencing in four diploid, prototrophic strains of Saccharomyces species grown under the same condition to investigate the evolution of gene expression. We found that variation in expression is gene-dependent with large variations in each gene's expression between replicates of the same species. This confounds the identification of genes differentially expressed across species. To address this, we developed a statistical approach to establish significance bounds for inter-species differential expression in RNA-Seq data based on the variance measured across biological replicates. This metric estimates the combined effects of technical and environmental variance, as well as Poisson sampling noise by isolating each component. Despite a paucity of large expression changes, we found a strong correlation between the variance of gene expression change and species divergence (R² = 0.90). CONCLUSION: We provide an improved methodology for measuring gene expression changes in evolutionary diverged species using RNA Seq, where experimental artifacts can mimic evolutionary effects.GEO Accession Number: GSE32679.

A comprehensive map of mobile element insertion polymorphisms in humans.

As a consequence of the accumulation of insertion events over evolutionary time, mobile elements now comprise nearly half of the human genome. The Alu, L1, and SVA mobile element families are still duplicating, generating variation between individual genomes. Mobile element insertions (MEI) have been identified as causes for genetic diseases, including hemophilia, neurofibromatosis, and various cancers. Here we present a comprehensive map of 7,380 MEI polymorphisms from the 1000 Genomes Project whole-genome sequencing data of 185 samples in three major populations detected with two detection methods. This catalog enables us to systematically study mutation rates, population segregation, genomic distribution, and functional properties of MEI polymorphisms and to compare MEI to SNP variation from the same individuals. Population allele frequencies of MEI and SNPs are described, broadly, by the same neutral ancestral processes despite vastly different mutation mechanisms and rates, except in coding regions where MEI are virtually absent, presumably due to strong negative selection. A direct comparison of MEI and SNP diversity levels suggests a differential mobile element insertion rate among populations.