Simultaneous multiplexed amplicon sequencing and transcriptome profiling in single cells

Mridusmita Saikia, Philip Burnham, Sara H Keshavjee, Michael F Z Wang, Michael Heyang, Pablo Moral-Lopez, Meleana M Hinchman, Charles G Danko, John S L Parker, Iwijn De Vlaminck

Preprint posted on November 13, 2018


High-throughput targeted long-read single cell sequencing reveals the clonal and transcriptional landscape of lymphocytes

Mandeep Singh, Ghamdan Al-Eryani, Shaun Carswell, James M. Ferguson, James Blackburn, Kirston Barton, Daniel Roden, Fabio Luciani, Tri Phan, Simon Junankar, Katherine Jackson, Christopher C Goodnow, Martin A Smith, Alexander Swarbrick

Preprint posted on September 24, 2018

What’s all the RAGE: DARTing into uncharted droplet-based scRNA-seq territory

Selected by Samantha Seah

Categories: genomics


Droplet-based single-cell RNA sequencing (scRNA-seq) is, and continues to be an important tool for estimating changes in gene expression at a single-cell level in a high-throughput manner. However, current techniques are largely limited to studying the 3’ ends of A-tailed mRNA transcripts at a transcriptome-wide level.

These two preprints by Saikia and Burnham et al. (on DART-seq) and Singh and Al-Eryani et al. (on RAGE-seq) tackle various limitations of current droplet-based scRNA-seq. For example, the diversity of B- and T- cell receptors (BCRs and TCRs respectively) generated by V(D)J recombination, and the further random addition or removal of nucleotides cannot be elucidated by traditional droplet-based scRNA-seq. The V(D)J regions of the B- and T-cell receptors are located at the 5’ end of their respective transcripts, while the combination of 3’ end capture, fragmentation and short-read Illumina sequencing in droplet-based scRNA-seq results in one getting sequences largely restricted to the 3’ end. The two preprints both tackle this problem, albeit slightly differently – DART-seq uses specific primers and modified Drop-seq beads to capture the heavy and light chain transcripts just downstream of the variable region, while RAGE-seq uses Oxford Nanopore technology to sequence the captured full-length BCR and TCR transcripts.

Additionally, Saikia and and Burnham et al. used DART-seq to characterise non-A-tailed transcripts of dsRNA viruses, while Singh and Al-Eryani et al.  used RAGE-seq to characterise alternative splicing and gene rearrangements in BCR transcripts, both of which would be extremely difficult with traditional droplet-based scRNA-seq.


DART-seq (Saikia and Burnham et al.)

In DART-seq, customised primers are ligated to a fraction of poly-dTs on Drop-seq beads, with an efficiency of 25-40%, which allows the specific capture and amplification of transcripts of interest. The ligation reaction can be titrated to leave many poly-dTs available, enabling the simultaneous capture of transcripts-of-interest and A-tailed transcripts.

The authors use their technology to specifically capture and sequence heavy and light chain antibody variable regions, while simultaneously obtaining whole-transcriptomic data from those same cells. Antibody variable regions are located at the 5’ end of their transcripts and traditionally cannot be sequenced with droplet-based scRNA-seq, which focuses on the 3’ end of transcripts. By designing primers just downstream of the V(D)J segments, the Ig recovery rate could be increased significantly compared to Drop-seq.

DART-seq was then used to study the B-cell antibody repertoire within human peripheral blood mononuclear cells (PBMCs). Of the 818 B-cells identified, 564 (67%) had immunoglobulin transcripts and the complete heavy and light chain CDR3 regions were obtained for 120 cells (15%). Clone-specific pairing was measured, and the highest frequency was observed between the most highly expressed heavy and light chain transcripts, which supports previous data.

The authors also applied their technology to study the infection of cells with T3D reovirus, an RNA virus with non-A-tailed transcripts. Infected cells were subjected to DART-seq, with beads modified to capture segments of the viral genome. By using primers designed to sequence the entire S2 transcript, the authors were able to analyse point mutations in the viral transcripts. Upon study of the cell transcriptomes, they found four distinct cell subpopulations, one of which was not found in the non-infected control.



RAGE-seq (Singh and Al-Eryani et al.)

In contrast to DART-seq, RAGE-seq uses targeted capture coupled with long-read sequencing with Oxford Nanopore. To obtain full-length BCR and TCR sequences, a capture bait library targeting all V, J and constant (C) region exons was used. The cDNA library is simultaneously sequenced via traditional short-read sequencing (via 10X Genomics), enabling 3’ expression profiling.

Figure 1 of Singh and Al-Eryani et al., made available under a CC-BY-NC-ND 4.0 International License


The authors validated the technology with a mix of Jurkat (T cells), Ramos (B cells) and monocytes (as a negative control). They noted a 13-fold enrichment of nanopore reads aligned to the BCR and TCR, when compared to non-targeted capture, and were able to recover both full length transcripts from a subset of cells – TCRα and β from 18.9% of Jurkat cells and the Ig heavy and light chains from 31% of Ramos cells. By tracking the amino acid mutations present, the authors could plot and follow the evolution of Ramos cells undergoing somatic hypermutation.

The authors then analysed lymphocytes from a human lymph node, where they similarly identified T cells and B cells, and sequenced their respective TCR and immunoglobulin chains. They quantified differences in Ig heavy chain transcripts in different cell types, noting that naive and memory B cells have both membrane and secretory IGH isoforms, while plasmablasts and plasma cells only have the secretory form. In addition, many plasmablasts and plasma cells were assigned to IGHA1, which is consistent with differentiation to high-rate antibody secreting cells.

The lymph node analysis was combined with analysis of lymphocytes from the primary breast tumour of the same patient, to enable the tracking of clonally-related T and B cells across different tissues. 7 clones were found to be shared between the tumour and lymph node, of which 6 were found within the CD8 T cell cluster. These clonal expanded cells have a discrete gene signature associated to active tissue resident cytotoxic lymphocytes, but each clone also seemed to express unique sets of genes.


What I like about these works

High throughput droplet-based transcriptomic sequencing is an extremely powerful technique, but currently suffers from certain limitations, such as the inability to sequence nucleotides located at the 5’ end of transcripts and the lack of sequencing strategies to pick up non-A-tailed transcripts.

The developments proposed by Saikia and Burnham et al. and Singh and Al-Eryani et al. go a long way to solve these problems. The two groups of authors have come up with innovative and contrasting technologies that expand the single-cell transcriptomic toolkit. In particular, both preprints have demonstrated an ability to recover natively-paired heavy and light chains of antibody variable regions, which is difficult to obtain with traditional droplet-based single-cell transcriptomic sequencing.

I’m also enthralled by the recent increase in utilisation of Oxford Nanopore sequencing, and I find the combination of Oxford Nanopore and Illumina sequencing in RAGE-seq (Singh and Al-Eryani et al.) exciting and potentially useful for many applications.

Tags: 10x genomics, antibody sequencing, drop-seq, droplet microfluidics, rna sequencing, single-cell sequencing

Posted on: 20th November 2018

(No Ratings Yet)

  • Have your say

    Your email address will not be published. Required fields are marked *

    This site uses Akismet to reduce spam. Learn how your comment data is processed.

    Sign up to customise the site to your preferences and to receive alerts

    Register here

    Also in the genomics category:

    The microbial basis of impaired wound healing: differential roles for pathogens, "bystanders", and strain-level diversification in clinical outcomes

    Lindsay Kalan, Jacquelyn S Meisel, Michael A Loesche, et al.

    Selected by Snehal Kadam

    Comparative analysis of droplet-based ultra-high-throughput single-cell RNA-seq systems

    Xiannian Zhang, Tianqi Li, Feng Liu, et al.

    Selected by Samantha Seah

    PUMILIO hyperactivity drives premature aging of Norad-deficient mice

    Florian Kopp, Mehmet Yalvac, Beibei Chen, et al.

    Selected by Carmen Adriaens

    LCM-seq reveals unique transcriptional adaption mechanisms of resistant neurons in spinal muscular atrophy

    Susanne Nichterwitz, Helena Storvall, Jik Nijssen, et al.


    Axon-seq decodes the motor axon transcriptome and its modulation in response to ALS

    Jik Nijssen, Julio Cesar Aguila Benitez, Rein Hoogstraaten, et al.

    Selected by Yen-Chung Chen

    LADL: Light-activated dynamic looping for endogenous gene expression control

    Mayuri Rege, Ji Hun Kim, Jacqueline Valeri, et al.

    Selected by Ivan Candido-Ferreira

    Precise tuning of gene expression output levels in mammalian cells

    Yale S. Michaels, Mike B Barnkob, Hector Barbosa, et al.

    Selected by Tim Fessenden


    Template switching causes artificial junction formation and false identification of circular RNAs

    Chong Tang, Tian Yu, Yeming Xie, et al.

    Selected by Fabio Liberante

    The genomic basis of colour pattern polymorphism in the harlequin ladybird

    Mathieu Gautier, Junichi Yamaguchi, Julien Foucaud, et al.

    Selected by Fillip Port

    Widespread inter-individual gene expression variability in Arabidopsis thaliana

    Sandra Cortijo, Zeynep Aydin, Sebastian Ahnert, et al.

    Selected by Martin Balcerowicz

    Single-cell Map of Diverse Immune Phenotypes Driven by the Tumor Microenvironment

    Elham Azizi, Ambrose J. Carr, George Plitas, et al.

    Selected by Tim Fessenden

    Cell type-specific interchromosomal interactions as a mechanism for transcriptional diversity

    Adan Horta, Kevin Monahan, Lisa Bashkirova, et al.

    Selected by Boyan Bonev

    PDX Finder: A Portal for Patient-Derived tumor Xenograft Model Discovery

    Nathalie Conte, Jeremy Mason, Csaba Halmagyi, et al.

    Selected by Carmen Adriaens

    An atlas of silencer elements for the human and mouse genomes

    Naresh Doni Jayavelu, Ajay Jajodia, Arpit Mishra, et al.

    Selected by Rafael Galupa


    Capturing the onset of PRC2-mediated repressive domain formation

    Ozgur Oksuz, Varun Narendra, Chul-Hwan Lee, et al.

    Selected by Boyan Bonev

    Heterochromatin drives organization of conventional and inverted nuclei

    Martin Falk, Yana Feodorova, Natasha Naumova, et al.

    Selected by Boyan Bonev

    The ancestral animal genetic toolkit revealed by diverse choanoflagellate transcriptomes

    Daniel Richter, Parinaz Fozouni, Michael Eisen, et al.

    Selected by Rafael Galupa