The lion's share of bacteria in various environments cannot be cloned in the laboratory and thus cannot be sequenced using existing technologies. A major goal of single-cell genomics is to complement gene-centric metagenomic data with whole-genome assemblies of uncultivated organisms. Assembly of single-cell data is challenging because of highly non-uniform read coverage as well as elevated levels of sequencing errors and chimeric reads. We describe SPAdes, a new assembler for both single-cell and standard (multicell) assembly, and demonstrate that it improves on the recently released E+V-SC assembler (specialized for single-cell data) and on popular assemblers Velvet and SoapDeNovo (for multicell data). SPAdes generates single-cell assemblies, providing information about genomes of uncultivatable bacteria that vastly exceeds what may be obtained via traditional metagenomics studies. SPAdes is available online ( http://bioinf.spbau.ru/spades ). It is distributed as open source software.

SPAdes, a new genome assembly algorithm and its applications to single-cell sequencing ( 7th Annual SFAF Meeting, 2012)

Motivation Recent advances in sequencing technologies promise ultra-long reads of ∼100 kb in average, full-length mRNA or cDNA reads in high throughput and genomic contigs over 100 Mb in length. Existing alignment programs are unable or inefficient to process such data at scale, which presses for the development of new alignment algorithms. Results Minimap2 is a general-purpose alignment program to map DNA or long mRNA sequences against a large reference database. It works with accurate short reads of ≥100 bp in length, ≥1 kb genomic reads at error rate ∼15%, full-length noisy Direct RNA or cDNA reads and assembly contigs or closely related full chromosomes of hundreds of megabases in length. Minimap2 does split-read alignment, employs concave gap cost for long insertions and deletions and introduces new heuristics to reduce spurious alignments. It is 3-4 times as fast as mainstream short-read mappers at comparable accuracy, and is ≥30 times faster than long-read genomic or cDNA mappers at higher accuracy, surpassing most aligners specialized in one type of alignment. Availability and implementation https://github.com/lh3/minimap2. Supplementary information Supplementary data are available at Bioinformatics online.

Minimap2: pairwise alignment for nucleotide sequences

RNA abundance is a powerful indicator of the state of individual cells. Single-cell RNA sequencing can reveal RNA abundance with high quantitative accuracy, sensitivity and throughput1. However, this approach captures only a static snapshot at a point in time, posing a challenge for the analysis of time-resolved phenomena such as embryogenesis or tissue regeneration. Here we show that RNA velocity-the time derivative of the gene expression state-can be directly estimated by distinguishing between unspliced and spliced mRNAs in common single-cell RNA sequencing protocols. RNA velocity is a high-dimensional vector that predicts the future state of individual cells on a timescale of hours. We validate its accuracy in the neural crest lineage, demonstrate its use on multiple published datasets and technical platforms, reveal the branching lineage tree of the developing mouse hippocampus, and examine the kinetics of transcription in human embryonic brain. We expect RNA velocity to greatly aid the analysis of developmental lineages and cellular dynamics, particularly in humans.

/pdf/rna-velocity-of-single-cells-10pdgilahj.pdf

RNA velocity of single cells

Most physiology and behavior of mammalian organisms follow daily oscillations. These rhythmic processes are governed by environmental cues (e.g., fluctuations in light intensity and temperature), an internal circadian timing system, and the interaction between this timekeeping system and environmental signals. In mammals, the circadian timekeeping system has a complex architecture, composed of a central pacemaker in the brain's suprachiasmatic nuclei (SCN) and subsidiary clocks in nearly every body cell. The central clock is synchronized to geophysical time mainly via photic cues perceived by the retina and transmitted by electrical signals to SCN neurons. In turn, the SCN influences circadian physiology and behavior via neuronal and humoral cues and via the synchronization of local oscillators that are operative in the cells of most organs and tissues. Thus, some of the SCN output pathways serve as input pathways for peripheral tissues. Here we discuss knowledge acquired during the past few years on the complex structure and function of the mammalian circadian timing system.

/pdf/the-mammalian-circadian-timing-system-organization-and-54vfplckht.pdf

The mammalian circadian timing system: organization and coordination of central and peripheral clocks.

The circadian system of mammals is composed of a hierarchy of oscillators that function at the cellular, tissue, and systems levels. A common molecular mechanism underlies the cell-autonomous circadian oscillator throughout the body, yet this clock system is adapted to different functional contexts. In the central suprachiasmatic nucleus (SCN) of the hypothalamus, a coupled population of neuronal circadian oscillators acts as a master pacemaker for the organism to drive rhythms in activity and rest, feeding, body temperature, and hormones. Coupling within the SCN network confers robustness to the SCN pacemaker, which in turn provides stability to the overall temporal architecture of the organism. Throughout the majority of the cells in the body, cell-autonomous circadian clocks are intimately enmeshed within metabolic pathways. Thus, an emerging view for the adaptive significance of circadian clocks is their fundamental role in orchestrating metabolism.

/pdf/central-and-peripheral-circadian-clocks-in-mammals-2qtn9mbv5i.pdf

Central and Peripheral Circadian Clocks in Mammals

Methods, machine readable media and kits are provided for analyzing an antibody repertoire of a subject to determine the immunocompetence of the subject. The method may include analyzing the subject's antibody repertoire. Analyzing may include identifying the composition of abundant related IgH sequences of the subject's antibody repertoire. The method may further include providing an assessment of the immunocompetence of the subject, based on the analysis.

Monitoring immunocompetance using immune repertoire sequencing

High-throughput short-read sequencing has taken on a central role in research and diagnostics. Literally hundreds of different assays exist today to 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 inertia associated with the research enterprise as a whole 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 long-read ONT MinION by circularizing and rolling circle amplifying the short library molecules using the R2C2 method. 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 millions of 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, as well as 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.

https://www.biorxiv.org/content/biorxiv/early/2021/10/30/2021.10.30.466545.full.pdf

Illumina But With Nanopore: Sequencing Illumina libraries at high accuracy on the ONT MinION using R2C2

Background RNA-Seq has brought forth significant discoveries regarding aberrations in RNA processing, implicating these RNA variants in a variety of diseases. Aberrant splicing and single nucleotide variants in RNA have been demonstrated to alter transcript stability, localization, and function. In particular, the upregulation of ADAR, an enzyme which mediates adenosine-to-inosine editing, has been previously linked to an increase in the invasiveness of lung ADC cells and associated with splicing regulation. Despite the functional importance of studying splicing and SNVs, short read RNA-Seq has limited the community’s ability to interrogate both forms of RNA variation simultaneously. Results We employed long-read technology to obtain full-length transcript sequences, elucidating cis-effects of variants on splicing changes at a single molecule level. We have developed a computational workflow that augments FLAIR, a tool that calls isoform models expressed in long-read data, to integrate RNA variant calls with the associated isoforms that bear them. We generated nanopore data with high sequence accuracy of H1975 lung adenocarcinoma cells with and without knockdown of ADAR. We applied our workflow to identify key inosine-isoform associations to help clarify the prominence of ADAR in tumorigenesis. Conclusions Ultimately, we find that a long-read approach provides valuable insight toward characterizing the relationship between RNA variants and splicing patterns. Highlights FLAIR2 has improved transcript isoform detection and incorporates sequence variants for haplotype-specific transcript detection. In addition to haplotype-specific variant detection, it identifies transcript-specific RNA editing Able to identify haplotype-specific transcript isoform bias in expression Long-read sequencing identifies hyperedited transcripts that are missed from short-read sequencing methods for a more comprehensive identification of ADAR targets

Detecting haplotype-specific transcript variation in long reads with FLAIR2

High-throughput short-read sequencing has taken on a central role in research and diagnostics. Literally hundreds of different assays exist today to 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 inertia associated with the research enterprise as a whole 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 long-read ONT MinION by circularizing and rolling circle amplifying the short library molecules using the R2C2 method. 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 millions of 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, as well as 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. be used to evaluate RNA-seq, ChIP-seq and enriched Tn5 libraries.

Illumina But With Nanopore

Promoters and the noncoding sequences that drive their function are fundamental aspects of genes that are critical to their regulation. The transcription preinitiation complex binds and assembles on promoters where it facilitates transcription. The transcription start site (TSS) is located downstream of the promoter sequence and is defined as the location in the genome where polymerase begins transcribing DNA into RNA. Knowing the location of TSSs is useful for annotation of genes, identification of non‐coding sequences important to gene regulation, detection of alternative TSSs, and understanding of 5′ UTR content. Several existing techniques make it possible to accurately identify TSSs, but are often difficult to perform experimentally, require large amounts of input RNA, or are unable to identify a large number of TSSs from a single sample. Many of these protocols take advantage of template switching reverse transcriptases (TSRTs), which reliably place an adaptor at the 5′ end of a first strand synthesis of cDNA. Here, we introduce a protocol that exploits TSRT activity combined with rolling circle amplification to identify TSSs with several unique advantages over existing methods. Sequence adaptors are placed on the 5′ and 3′ end of the full‐length cDNA copy of a transcript. A splint compatible with those adaptors is then used to circularize the full‐length cDNA. Linear DNA containing concatemers of the cDNA are generated using rolling circle amplification, and a sequencing library is formed by fragmenting the concatemers. This protocol is straightforward to execute, requiring limited bench time with relatively stable reagents. Using extremely low amounts of RNA input, this protocol produces large numbers of accurate, deduplicated TSSs genome wide. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC.

https://onlinelibrary.wiley.com/doi/pdfdirect/10.1002/cpz1.705

Christopher Vollmers

Papers

Monitoring immunocompetance using immune repertoire sequencing

Illumina But With Nanopore: Sequencing Illumina libraries at high accuracy on the ONT MinION using R2C2

Detecting haplotype-specific transcript variation in long reads with FLAIR2

Illumina But With Nanopore

Smar2C2: A Simple and Efficient Protocol for the Identification of Transcription Start Sites