Long read subcellular fractionation and sequencing reveals the translational fate of full-length mRNA isoforms during neuronal differentiation.

Alternative splicing (AS) alters the cis-regulatory landscape of mRNA isoforms leading to transcripts with distinct localization, stability and translational efficiency. To rigorously investigate mRNA isoform-specific ribosome association, we generated subcellular fractionation and sequencing (Frac-seq) libraries using both conventional short reads and long reads from human embryonic stem cells (ESC) and neural progenitor cells (NPC) derived from the same ESC. We performed de novo transcriptome assembly from high-confidence long reads from cytosolic, monosomal, light and heavy polyribosomal fractions and quantified their abundance using short reads from their respective subcellular fractions. Thousands of transcripts in each cell type exhibited association with particular subcellular fractions relative to the cytosol. Of the multi-isoform genes, 27% and 19% exhibited significant differential isoform sedimentation in ESC and NPC respectively. Alternative promoter usage and internal exon skipping accounted for the majority of differences between isoforms from the same gene. Random forest classifiers implicated coding sequence (CDS) and UTR lengths as important determinants of isoform-specific sedimentation profiles, and motif analyses reveal potential cell type-specific and subcellular fraction-associated RNA-binding protein signatures. Taken together our data demonstrate that alternative mRNA processing within the CDS and UTRs impacts the translational control of mRNA isoforms during stem cell differentiation, and highlights the utility of using a novel long-read sequencing-based method to study translational control.

Sequencing Illumina libraries at high accuracy on the ONT MinION using R2C2.

High-throughput short-read sequencing has taken on a central role in research and diagnostics. Hundreds of different assays take advantage of Illumina short-read sequencers, the predominant short-read sequencing technology available today. Although other short-read sequencing technologies exist, the ubiquity of Illumina sequencers in sequencing core facilities and the high capital costs of these technologies have limited their adoption. Among a new generation of sequencing technologies, Oxford Nanopore Technologies (ONT) holds a unique position because the ONT MinION, an error-prone long-read sequencer, is associated with little to no capital cost. Here we show that we can make short-read Illumina libraries compatible with the ONT MinION by using the rolling circle to concatemeric consensus (R2C2) method to circularize and amplify the short library molecules. This results in longer DNA molecules containing tandem repeats of the original short library molecules. This longer DNA is ideally suited for the ONT MinION, and after sequencing, the tandem repeats in the resulting raw reads can be converted into high-accuracy consensus reads with similar error rates to that of the Illumina MiSeq. We highlight this capability by producing and benchmarking RNA-seq, ChIP-seq, and regular and target-enriched Tn5 libraries. We also explore the use of this approach for rapid evaluation of sequencing library metrics by implementing a real-time analysis workflow.

A conserved long noncoding RNA, GAPLINC, modulates the immune response during endotoxic shock.

Recent studies have identified thousands of long noncoding RNAs (lncRNAs) in mammalian genomes that regulate gene expression in different biological processes. Although lncRNAs have been identified in a variety of immune cells and implicated in immune response, the biological function and mechanism of the majority remain unexplored, especially in sepsis. Here, we identify a role for a lncRNA-gastric adenocarcinoma predictive long intergenic noncoding RNA (GAPLINC)-previously characterized for its role in cancer, now in the context of innate immunity, macrophages, and LPS-induced endotoxic shock. Transcriptome analysis of macrophages from humans and mice reveals that GAPLINC is a conserved lncRNA that is highly expressed following macrophage differentiation. Upon inflammatory activation, GAPLINC is rapidly down-regulated. Macrophages depleted of GAPLINC display enhanced expression of inflammatory genes at baseline, while overexpression of GAPLINC suppresses this response. Consistent with GAPLINC-depleted cells, Gaplinc knockout mice display enhanced basal levels of inflammatory genes and show resistance to LPS-induced endotoxic shock. Mechanistically, survival is linked to increased levels of nuclear NF-κB in Gaplinc knockout mice that drives basal expression of target genes typically only activated following inflammatory stimulation. We show that this activation of immune response genes prior to LPS challenge leads to decreased blood clot formation, which protects Gaplinc knockout mice from multiorgan failure and death. Together, our results identify a previously unknown function for GAPLINC as a negative regulator of inflammation and uncover a key role for this lncRNA in modulating endotoxic shock.

Complete characterization of the human immune cell transcriptome using accurate full-length cDNA sequencing.

The human immune system relies on highly complex and diverse transcripts and the proteins they encode. These include transcripts encoding human leukocyte antigen (HLA) receptors as well as B cell and T cell receptors (BCR and TCR). Determining which alleles an individual possesses for each HLA gene (high-resolution HLA typing) is essential to establish donor-recipient compatibility in organ and bone marrow transplantations. In turn, the repertoires of millions of unique BCR and TCR transcripts in each individual carry a vast amount of health-relevant information. Both short-read RNA-seq-based HLA typing and BCR/TCR repertoire sequencing (AIRR-seq) currently rely on our incomplete knowledge of the genetic diversity at HLA and BCR/TCR loci. Here, we generated over 10,000,000 full-length cDNA sequences at a median accuracy of 97.9% using our nanopore sequencing-based Rolling Circle Amplification to Concatemeric Consensus (R2C2) protocol. We used this data set to (1) show that deep and accurate full-length cDNA sequencing can be used to provide isoform-level transcriptome analysis for more than 9000 loci, (2) generate accurate sequences of HLA alleles, and (3) extract detailed AIRR data for the analysis of the adaptive immune system. The HLA and AIRR analysis approaches we introduce here are untargeted and therefore do not require prior knowledge of the composition or genetic diversity of HLA and BCR/TCR loci.

Whole-genome sequence and assembly of the Javan gibbon (Hylobates moloch).

The Javan gibbon, Hylobates moloch, is an endangered gibbon species restricted to the forest remnants of western and central Java, Indonesia, and one of the rarest of the Hylobatidae family. Hylobatids consist of 4 genera (Holoock, Hylobates, Symphalangus, and Nomascus) that are characterized by different numbers of chromosomes, ranging from 38 to 52. The underlying cause of this karyotype plasticity is not entirely understood, at least in part, due to the limited availability of genomic data. Here we present the first scaffold-level assembly for H. moloch using a combination of whole-genome Illumina short reads, 10X Chromium linked reads, PacBio, and Oxford Nanopore long reads and proximity-ligation data. This Hylobates genome represents a valuable new resource for comparative genomics studies in primates.

Horizontal transmission and recombination maintain forever young bacterial symbiont genomes.

Bacterial symbionts bring a wealth of functions to the associations they participate in, but by doing so, they endanger the genes and genomes underlying these abilities. When bacterial symbionts become obligately associated with their hosts, their genomes are thought to decay towards an organelle-like fate due to decreased homologous recombination and inefficient selection. However, numerous associations exist that counter these expectations, especially in marine environments, possibly due to ongoing horizontal gene flow. Despite extensive theoretical treatment, no empirical study thus far has connected these underlying population genetic processes with long-term evolutionary outcomes. By sampling marine chemosynthetic bacterial-bivalve endosymbioses that range from primarily vertical to strictly horizontal transmission, we tested this canonical theory. We found that transmission mode strongly predicts homologous recombination rates, and that exceedingly low recombination rates are associated with moderate genome degradation in the marine symbionts with nearly strict vertical transmission. Nonetheless, even the most degraded marine endosymbiont genomes are occasionally horizontally transmitted and are much larger than their terrestrial insect symbiont counterparts. Therefore, horizontal transmission and recombination enable efficient natural selection to maintain intermediate symbiont genome sizes and substantial functional genetic variation.

Generation of an isoform-level transcriptome atlas of macrophage activation.

RNA-seq is routinely used to measure gene expression changes in response to cell perturbation. Genes upregulated or downregulated following some perturbation are designated as genes of interest, and their most expressed isoform(s) would then be selected for follow-up experimentation. However, because of its need to fragment RNA molecules, RNA-seq is limited in its ability to capture gene isoforms and their expression patterns. This lack of isoform-specific data means that isoforms would be selected based on annotation databases that are incomplete, not tissue specific, or do not provide key information on expression levels. As a result, minority or nonexistent isoforms might be selected for follow-up, leading to loss in valuable resources and time. There is therefore a great need to comprehensively identify gene isoforms along with their corresponding levels of expression. Using the long-read nanopore-based R2C2 method, which does not fragment RNA molecules, we generated an Isoform-level transcriptome Atlas of Macrophage Activation that identifies full-length isoforms in primary human monocyte-derived macrophages. Macrophages are critical innate immune cells important for recognizing pathogens through binding of pathogen-associated molecular patterns to toll-like receptors, culminating in the initiation of host defense pathways. We characterized isoforms for most moderately-to-highly expressed genes in resting and toll-like receptor-activated monocyte-derived macrophages, identified isoforms differentially expressed between conditions, and validated these isoforms by RT-qPCR. We compiled these data into a user-friendly data portal within the UCSC Genome Browser (https://genome.ucsc.edu/s/vollmers/IAMA). Our atlas represents a valuable resource for innate immune research, providing unprecedented isoform information for primary human macrophages.

One fly-one genome: chromosome-scale genome assembly of a single outbred Drosophila melanogaster.

A high quality genome assembly is a vital first step for the study of an organism. Recent advances in technology have made the creation of high quality chromosome scale assemblies feasible and low cost. However, the amount of input DNA needed for an assembly project can be a limiting factor for small organisms or precious samples. Here we demonstrate the feasibility of creating a chromosome scale assembly using a hybrid method for a low input sample, a single outbred Drosophila melanogaster. Our approach combines an Illumina shotgun library, Oxford nanopore long reads, and chromosome conformation capture for long range scaffolding. This single fly genome assembly has a N50 of 26 Mb, a length that encompasses entire chromosome arms, contains 95% of expected single copy orthologs, and a nearly complete assembly of this individual's Wolbachia endosymbiont. The methods described here enable the accurate and complete assembly of genomes from small, field collected organisms as well as precious clinical samples.

Signatures of selection in the human antibody repertoire: Selective sweeps, competing subclones, and neutral drift.

Antibodies are created and refined by somatic evolution in B cell populations, which endows the human immune system with the ability to recognize and eliminate diverse pathogens. However, the evolutionary processes that sculpt antibody repertoires remain poorly understood. Here, using an unbiased repertoire-scale approach, we show that the population genetic signatures of evolution are evident in human B cell lineages and reveal how antibodies evolve somatically. We measured the dynamics and genetic diversity of B cell responses in five adults longitudinally before and after influenza vaccination using high-throughput antibody repertoire sequencing. We identified vaccine-responsive B cell lineages that carry signatures of selective sweeps driven by positive selection, and discovered that they often display evidence for selective sweeps favoring multiple subclones. We also found persistent B cell lineages that exhibit stable population dynamics and carry signatures of neutral drift. By exploiting the relationship between B cell fitness and antibody binding affinity, we demonstrate the potential for using phylogenetic approaches to identify antibodies with high binding affinity. This quantitative characterization reveals that antibody repertoires are shaped by an unexpectedly broad spectrum of evolutionary processes and shows how signatures of evolutionary history can be harnessed for antibody discovery and engineering.

An Annotated Draft Genome for the Andean Bear, Tremarctos ornatus.

The Andean bear is the only extant member of the Tremarctine subfamily and the only extant ursid species to inhabit South America. Here, we present an annotated de novo assembly of a nuclear genome from a captive-born female Andean bear, Mischief, generated using a combination of short and long DNA and RNA reads. Our final assembly has a length of 2.23 Gb, and a scaffold N50 of 21.12 Mb, contig N50 of 23.5 kb, and BUSCO score of 88%. The Andean bear genome will be a useful resource for exploring the complex phylogenetic history of extinct and extant bear species and for future population genetics studies of Andean bears.

Improving nanopore read accuracy with the R2C2 method enables the sequencing of highly multiplexed full-length single-cell cDNA.

High-throughput short-read sequencing has revolutionized how transcriptomes are quantified and annotated. However, while Illumina short-read sequencers can be used to analyze entire transcriptomes down to the level of individual splicing events with great accuracy, they fall short of analyzing how these individual events are combined into complete RNA transcript isoforms. Because of this shortfall, long-distance information is required to complement short-read sequencing to analyze transcriptomes on the level of full-length RNA transcript isoforms. While long-read sequencing technology can provide this long-distance information, there are issues with both Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) long-read sequencing technologies that prevent their widespread adoption. Briefly, PacBio sequencers produce low numbers of reads with high accuracy, while ONT sequencers produce higher numbers of reads with lower accuracy. Here, we introduce and validate a long-read ONT-based sequencing method. At the same cost, our Rolling Circle Amplification to Concatemeric Consensus (R2C2) method generates more accurate reads of full-length RNA transcript isoforms than any other available long-read sequencing method. These reads can then be used to generate isoform-level transcriptomes for both genome annotation and differential expression analysis in bulk or single-cell samples.

Linked genetic variation and not genome structure causes widespread differential expression associated with chromosomal inversions.

Chromosomal inversions are widely thought to be favored by natural selection because they suppress recombination between alleles that have higher fitness on the same genetic background or in similar environments. Nonetheless, few selected alleles have been characterized at the molecular level. Gene expression profiling provides a powerful way to identify functionally important variation associated with inversions and suggests candidate phenotypes. However, altered genome structure itself might also impact gene expression by influencing expression profiles of the genes proximal to inversion breakpoint regions or by modifying expression patterns genome-wide due to rearranging large regulatory domains. In natural inversions, genetic differentiation and genome structure are inextricably linked. Here, we characterize differential expression patterns associated with two chromosomal inversions found in natural Drosophila melanogaster populations. To isolate the impacts of genome structure, we engineered synthetic chromosomal inversions on controlled genetic backgrounds with breakpoints that closely match each natural inversion. We find that synthetic inversions have negligible effects on gene expression. Nonetheless, natural inversions have broad-reaching regulatory impacts in cis and trans Furthermore, we find that differentially expressed genes associated with both natural inversions are enriched for loci associated with immune response to bacterial pathogens. Our results support the idea that inversions in D. melanogaster experience natural selection to maintain associations between functionally related alleles to produce complex phenotypic outcomes.

Tn5Prime, a Tn5 based 5' capture method for single cell RNA-seq.

RNA-sequencing (RNA-seq) is a powerful technique to investigate and quantify entire transcriptomes. Recent advances in the field have made it possible to explore the transcriptomes of single cells. However, most widely used RNA-seq protocols fail to provide crucial information regarding transcription start sites. Here we present a protocol, Tn5Prime, that takes advantage of the Tn5 transposase-based Smart-seq2 protocol to create RNA-seq libraries that capture the 5' end of transcripts. The Tn5Prime method dramatically streamlines the 5' capture process and is both cost effective and reliable. By applying Tn5Prime to bulk RNA and single cell samples, we were able to define transcription start sites as well as quantify transcriptomes at high accuracy and reproducibility. Additionally, similar to 3' end-based high-throughput methods like Drop-seq and 10× Genomics Chromium, the 5' capture Tn5Prime method allows the introduction of cellular identifiers during reverse transcription, simplifying the analysis of large numbers of single cells. In contrast to 3' end-based methods, Tn5Prime also enables the assembly of the variable 5' ends of the antibody sequences present in single B-cell data. Therefore, Tn5Prime presents a robust tool for both basic and applied research into the adaptive immune system and beyond.

Phylogenetic analysis of the human antibody repertoire reveals quantitative signatures of immune senescence and aging.

The elderly have reduced humoral immunity, as manifested by increased susceptibility to infections and impaired vaccine responses. To investigate the effects of aging on B-cell receptor (BCR) repertoire evolution during an immunological challenge, we used a phylogenetic distance metric to analyze Ig heavy-chain transcript sequences in both young and elderly individuals before and after influenza vaccination. We determined that BCR repertoires become increasingly specialized over a span of decades, but less plastic. In 50% of the elderly individuals, a large space in the repertoire was occupied by a small number of recall lineages that did not decline during vaccine response and contained hypermutated IgD+ B cells. Relative to their younger counterparts, older subjects demonstrated a contracted naive repertoire and diminished intralineage diversification, signifying a reduced substrate for mounting novel responses and decreased fine-tuning of BCR specificities by somatic hypermutation. Furthermore, a larger proportion of the repertoire exhibited premature stop codons in some elderly subjects, indicating that aging may negatively affect the ability of B cells to discriminate between functional and nonfunctional receptors. Finally, we observed a decreased incidence of radical mutations compared with conservative mutations in elderly subjects' vaccine responses, which suggests that accumulating original antigenic sin may be limiting the accessible space for paratope evolution. Our findings shed light on the complex interplay of environmental and gerontological factors affecting immune senescence, and provide direct molecular characterization of the effects of senescence on the immune repertoire.

Combining ATAC-seq with nuclei sorting for discovery of cis-regulatory regions in plant genomes.

Chromatin structure plays a pivotal role in facilitating proper control of gene expression. Transcription factor (TF) binding of cis-elements is often associated with accessible chromatin regions. Therefore, the ability to identify these accessible regions throughout plant genomes will advance understanding of the relationship between TF binding, chromatin status and the regulation of gene expression. Assay for Transposase Accessible Chromatin sequencing (ATAC-seq) is a recently developed technique used to map open chromatin zones in animal genomes. However, in plants, the existence of cell walls, subcellular organelles and the lack of stable cell lines have prevented routine application of this technique. Here, we describe an assay combining ATAC-seq with fluorescence-activated nuclei sorting (FANS) to identify and map open chromatin and TF-binding sites in plant genomes. FANS-ATAC-seq compares favorably with published DNaseI sequencing (DNase-seq) results and it requires less than 50 000 nuclei for accurate identification of accessible genomic regions. SUMMARY: Application of ATAC-seq to sorted nuclei identifies accessible regions genome-wide.

Highly Accurate Sequencing of Full-Length Immune Repertoire Amplicons Using Tn5-Enabled and Molecular Identifier-Guided Amplicon Assembly.

Ab repertoire sequencing is a powerful tool to analyze the adaptive immune system. To sequence entire Ab repertoires, amplicons are created from Ab H chain (IgH) transcripts and sequenced on a high-throughput sequencer. The field of immune repertoire sequencing is growing rapidly and the protocols used are steadily improving; however, thus far, immune repertoire sequencing protocols have not been able to sequence full-length immune repertoires including the entire IgH V region and enough of the IgH C region to identify isotype subtypes. In this study, we present a method that combines Tn5 transposase and molecular identifiers for the highly accurate sequencing of amplicons >500 bp using Illumina short read paired-end sequencing. We then apply this method to Ab H chain amplicons to sequence the first, to our knowledge, highly accurate full-length immune repertoire.

Local circadian clock gates cell cycle progression of transient amplifying cells during regenerative hair cycling.

Regenerative cycling of hair follicles offers an unique opportunity to explore the role of circadian clock in physiological tissue regeneration. We focused on the role of circadian clock in actively proliferating transient amplifying cells, as opposed to quiescent stem cells. We identified two key sites of peripheral circadian clock activity specific to regenerating anagen hair follicles, namely epithelial matrix and mesenchymal dermal papilla. We showed that peripheral circadian clock in epithelial matrix cells generates prominent daily mitotic rhythm. As a consequence of this mitotic rhythmicity, hairs grow faster in the morning than in the evening. Because cells are the most susceptible to DNA damage during mitosis, this cycle leads to a remarkable time-of-day-dependent sensitivity of growing hair follicles to genotoxic stress. Same doses of γ-radiation caused dramatic hair loss in wild-type mice when administered in the morning, during mitotic peak, compared with the evening, when hair loss is minimal. This diurnal radioprotective effect becomes lost in circadian mutants, consistent with asynchronous mitoses in their hair follicles. Clock coordinates cell cycle progression with genotoxic stress responses by synchronizing Cdc2/Cyclin B-mediated G2/M checkpoint. Our results uncover diurnal mitotic gating as the essential protective mechanism in highly proliferative hair follicles and offer strategies for minimizing or maximizing cytotoxicity of radiation therapies.

Genetic measurement of memory B-cell recall using antibody repertoire sequencing.

Annual influenza vaccinations aim to protect against seasonal infections, and vaccine strain compositions are updated every year. This protection is based on antibodies that are produced by either newly activated or memory B cells recalled from previous encounters with influenza vaccination or infection. The extent to which the B-cell repertoire responds to vaccination and recalls antibodies has so far not been analyzed at a genetic level--which is to say, at the level of antibody sequences. Here, we developed a consensus read sequencing approach that incorporates unique barcode labels on each starting RNA molecule. These labels allow one to combine multiple sequencing reads covering the same RNA molecule to reduce the error rate to a desired level, and they also enable accurate quantification of RNA and isotype levels. We validated this approach and analyzed the differential response of the antibody repertoire to live-attenuated or trivalent-inactivated influenza vaccination. Additionally, we analyzed the antibody repertoire in response to repeated yearly vaccinations with trivalent-inactivated influenza vaccination. We found antibody sequences that were present in both years, providing a direct genetic measurement of B-cell recall.

Monitoring pharmacologically induced immunosuppression by immune repertoire sequencing to detect acute allograft rejection in heart transplant patients: a proof-of-concept diagnostic accuracy study.

BACKGROUND: It remains difficult to predict and to measure the efficacy of pharmacological immunosuppression. We hypothesized that measuring the B-cell repertoire would enable assessment of the overall level of immunosuppression after heart transplantation. METHODS AND FINDINGS: In this proof-of-concept study, we implemented a molecular-barcode-based immune repertoire sequencing assay that sensitively and accurately measures the isotype and clonal composition of the circulating B cell repertoire. We used this assay to measure the temporal response of the B cell repertoire to immunosuppression after heart transplantation. We selected a subset of 12 participants from a larger prospective cohort study (ClinicalTrials.gov NCT01985412) that is ongoing at Stanford Medical Center and for which enrollment started in March 2010. This subset of 12 participants was selected to represent post-heart-transplant events, with and without acute rejection (six participants with moderate-to-severe rejection and six without). We analyzed 130 samples from these patients, with an average follow-up period of 15 mo. Immune repertoire sequencing enables the measurement of a patient's net state of immunosuppression (correlation with tacrolimus level, r = -0.867, 95% CI -0.968 to -0.523, p = 0.0014), as well as the diagnosis of acute allograft rejection, which is preceded by increased immune activity with a sensitivity of 71.4% (95% CI 30.3% to 94.9%) and a specificity of 82.0% (95% CI 72.1% to 89.1%) (cell-free donor-derived DNA as noninvasive gold standard). To illustrate the potential of immune repertoire sequencing to monitor atypical post-transplant trajectories, we analyzed two more patients, one with chronic infections and one with amyloidosis. A larger, prospective study will be needed to validate the power of immune repertoire sequencing to predict rejection events, as this proof-of-concept study is limited to a small number of patients who were selected based on several criteria including the availability of a large number of samples and the absence or presence of rejection events. CONCLUSIONS: If confirmed in larger, prospective studies, the method described here has potential applications in the tailored management of post-transplant immunosuppression and, more broadly, as a method for assessing the overall activity of the immune system.

Generation and analysis of a mouse multi-tissue genome annotation atlas.

Generating an accurate and complete genome annotation for an organism is complex because the cells within each tissue can express a unique set of transcript isoforms from a unique set of genes. A comprehensive genome annotation should contain information on what tissues express what transcript isoforms at what level. This tissue-level isoform information can then inform a wide range of research questions as well as experiment designs. Long-read sequencing technology combined with advanced full-length cDNA library preparation methods has now achieved throughput and accuracy where generating these types of annotations is achievable. Here, we show this by generating a genome annotation of the mouse (Mus musculus). We used the nanopore-based R2C2 long-read sequencing method to generate 64 million highly accurate full length cDNA consensus reads - averaging 5.4 million reads per tissue for a dozen tissues. Using the Mandalorion tool we processed these reads to generate the Tissue-level Atlas of Mouse Isoforms (TAMI - available at https://genome.ucsc.edu/s/vollmers/TAMI) which we believe will be a valuable complement to conventional, manually curated reference genome annotations.