Molecular Biology

Many organisms alternate the expression of genes from large gene sets or gene families to adapt to environmental cues or immune pressure. The single-celled protozoan pathogen Trypanosoma brucei spp. periodically changes its homogeneous surface coat of variant surface glycoproteins (VSGs) to evade host antibodies during infection. This pathogen expresses one out of ~2,500 VSG genes at a time from telomeric expression sites (ESs) and periodically changes their expression by transcriptional switching or recombination. Attempts to track VSG switching have previously relied on genetic modifications of ES sequences with drug-selectable markers or genes encoding fluorescent proteins. However, genetic modifications of the ESs can interfere with the binding of proteins that control VSG transcription and/or recombination, thus affecting VSG expression and switching. Other approaches include Illumina sequencing of the VSG repertoire, which shows VSGs expressed in the population rather than cell switching; the Illumina short reads often limit the distinction of the large set of VSG genes. Here, we describe a methodology to study antigenic switching without modifications of the ES sequences. Our protocol enables the detection of VSG switching at nucleotide resolution using multiplexed clonal cell barcoding to track cells and nanopore sequencing to identify cell-specific VSG expression. We also developed a computational pipeline that takes DNA sequences and outputs VSGs expressed by cell clones. This protocol can be adapted to study clonal cell expression of large gene families in prokaryotes or eukaryotes.

Key features

• This protocol enables the analysis of variant surface glycoproteins (VSG) switching in T. brucei without modifying the expression site sequences.

• It uses a streamlined computational pipeline that takes fastq DNA sequences and outputs expressed VSG genes by each parasite clone.

• The protocol leverages the long reads sequencing capacity of the Oxford nanopore sequencing technology, which enables accurate identification of the expressed VSGs.

• The protocol requires approximately eight to nine days to complete.

Graphical overview

Rapid development in single-cell chromosome conformation capture technologies has provided valuable insights into the importance of spatial genome architecture for gene regulation. However, a long-standing technical gap remains in the simultaneous characterization of three-dimensional genomes and transcriptomes in the same cell. We have described an assay named Hi-C and RNA-seq employed simultaneously (HiRES), which integrates in situ reverse transcription and chromosome conformation capture (3C) for the parallel analysis of chromatin organization and gene expression. Here, we provide a detailed implementation of the assay, using mouse embryos and cerebral cortices as examples. The versatility of this method extends beyond these two samples, with the potential to be used in various other cell types.

Key features

• A multi-omics sequencing approach to profile 3D genome structure and gene expression simultaneously in single cells.

• Compatible with animal tissues.

• One-tube amplification of both DNA and RNA components.

• Requires three days to complete.

Graphical overview

Schematic illustration for the Hi-C and RNA-seq employed simultaneously (HiRES) workflow

Malaria molecular surveillance has great potential to support national malaria control programs (NMCPs), informing policy for its control and elimination. Here, we present a new three-day workflow for targeted resequencing of markers in 13 resistance-associated genes, histidine rich protein 2 and 3 (hrp2&3), a country (Peru)-specific 28 SNP-barcode for population genetic analysis, and apical membrane antigen 1 (ama1), using Illumina short-read sequencing technology. The assay applies a multiplex PCR approach to amplify all genomic regions of interest in a rapid and easily standardizable procedure and allows simultaneous amplification of a high number of targets at once, therefore having great potential for implementation into routine surveillance practice by NMCPs. The assay can be performed on routinely collected filter paper blood spots and can be easily adapted to different regions to investigate either regional trends or in-country epidemiological changes.

Spontaneous DNA damage frequently occurs on the human genome, and it could alter gene expression by inducing mutagenesis or epigenetic changes. Therefore, it is highly desired to profile DNA damage distribution on the human genome and identify the genes that are prone to DNA damage. Here, we present a novel single-cell whole-genome amplification method which employs linear-copying followed by a split-amplification scheme, to efficiently remove amplification errors and achieve accurate detection of DNA damage in individual cells. In comparison to previous methods that measure DNA damage, our method uses a next-generation sequencing platform to detect misincorporated bases derived from spontaneous DNA damage with single-cell resolution.

Cell-type specific transcriptional programs underlie the development and maintenance of organs. Not only distinct cell types within a tissue, even cells with supposedly identical cell fates show a high degree of transcriptional heterogeneity. Inevitable, low cell numbers are a major hurdle to study transcriptomes of pure cell populations. Here we describe DigiTAG, a high-throughput method that combines transposase fragmentation and molecular barcoding to retrieve high quality transcriptome data of rare cell types in Drosophila melanogaster. The protocol showcases how DigiTAG can be used to analyse the transcriptome of rare neural stem cells (type II neuroblasts) of Drosophila larval brains, but can also be utilized for other cell types or model systems.

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive hematological malignancy that arises from transformation of T-cell primed hematopoietic progenitors. Although T-ALL is a heterogenous and molecularly complex disease, more than 65% of T-ALL patients carry activating mutations in the NOTCH1 gene. The majority of T-ALL–associated NOTCH1 mutations either disrupt the negative regulatory region, allowing signal activation in the absence of ligand binding, or result in truncation of the C-terminal PEST domain involved in the termination of NOTCH1 signaling by proteasomal degradation. To date, retroviral transduction models have relied heavily on the overexpression of aggressively truncated variants of NOTCH1 (such as ICN1 or ΔE-NOTCH1), which result in supraphysiological levels of signaling activity and are rarely found in human T-ALL. The current protocol describes the method for mouse bone marrow isolation, hematopoietic stem and progenitor cell (HSC) enrichment, followed by retroviral transduction with an oncogenic mutant form of the NOTCH1 receptor (NOTCH1-L1601P-ΔP) that closely resembles the gain-of-function mutations most commonly found in patient samples. A hallmark of this forced expression of constitutively active NOTCH1 is a transient wave of extrathymic immature T-cell development, which precedes oncogenic transformation to T-ALL. Furthermore, this approach models leukemic transformation and progression in vivo by allowing for crosstalk between leukemia cells and the microenvironment, an aspect unaccounted for in cell-line based in vitro studies. Thus, the HSC transduction and transplantation model more faithfully recapitulates development of the human disease, providing a highly comprehensive and versatile tool for further in vivo and ex vivo functional studies.

Chromatin immunoprecipitation is extensively used to investigate the epigenetic profile and transcription factor binding sites in the genome. However, when the starting material is limited, the conventional ChIP-Seq approach cannot be implemented. This protocol describes a method that can be used to generate the chromatin profiles from as low as 100 human or 1,000 Drosophila cells. The method employs tagmentation to fragment the chromatin with concomitant addition of sequencing adaptors. The method generates datasets with high signal to noise ratio and can be subjected to standard tools for ChIP-Seq analysis

Most bacterial genomes have biased nucleotide composition, and the asymmetry is considered to be caused by a single-stranded DNA (ssDNA) deamination arising from the bacterial replication machinery. In order to evaluate the relationship experimentally, the position and frequency of ssDNA formed during replication must be verified clearly. Although many ssDNA detection technologies exist, almost all methods have been developed for eukaryotic genomes. To apply these to bacterial genomes, which harbor a smaller amount of DNA than those of eukaryotes, more efficient, new methods are required. Therefore, we developed a novel strand-specific ssDNA sequencing method, called 4S-seq, for the bacterial genome. The 4S-seq method enriches ssDNA in the extracted genomic DNA by a dsDNA-specific nuclease and implements a strand-specific library using a biotin label with a customized tag. As a result, the 4S-seq is able to calculate the ssDNA content in each strand (Watson/Crick) at each position of the genome efficiently.

Detecting heteroplasmies in the mitochondrial DNA (mtDNA) has been a challenge for many years. In the past, Sanger sequencing was the main option to perform this analysis, however, this method could not detect low frequency heteroplasmies. Massive Parallel Sequencing (MPS) provides the opportunity to study the mtDNA in depth, but a controlled pipeline is necessary to reliably retrieve and quantify the low frequency variants. It has been shown that differences in methods can significantly affect the number and frequency of the retrieved variants.

In this protocol, we present a method involving both wet lab and bioinformatics that allows identifying and quantifying single nucleotide variants in the full mtDNA sequence, down to a heteroplasmic load of 1.5%. For this, we set up a PCR-based amplification of the mtDNA, followed by MPS using Illumina chemistry, and variant calling with two different algorithms, mtDNA server and Mutect.

The PCR amplification is used to enrich the mitochondrial fraction, while the bioinformatic processing with two algorithms is used to discriminate the true heteroplasmies from background noise. The protocol described here allows for deep sequencing of the mitochondrial DNA in bulk DNA samples as well as single cells (both large cells such as human oocytes, and small-sized single cells such as human embryonic stem cells) with minor modifications to the protocol.

Precise genome editing is essential for scientific research and clinical application. At present, linear amplification-mediated high-throughput genome-wide translocation sequencing (LAM-HTGTS) is one of most effective methods to evaluate the off-target activity of CRISPR-Cas9, which is based on chromosomal translocation and employs a “bait” DNA double-stranded break (DSB) to capture genome-wide “prey” DNA DSBs. Here, we described an improved HTGTS (iHTGTS) method, in which size-selection beads were used to enhance reaction efficiency and a new primer system was designed to be compatible with Illumina Hiseq sequencing. Compared with LAM-HTGTS, iHTGTS is lower cost and has much higher sensitivity for off-target detection in HEK293T, K562, U2OS and HCT116 cell lines. So we believe that iHTGTS is a powerful method for comprehensively assessing Cas9 off-target effect.