系统生物学

Glomerular diseases characterized by injury to post-mitotic epithelial cells called podocytes are a leading cause of chronic kidney disease. Yet, isolating podocytes from the kidney for transcriptomic, proteomic, and metabolomic studies has been a major technical challenge. Protocols utilizing glomerular sieving and laser capture methods are of limited use because they are not podocyte-specific but instead capture all four glomerular cell types. Here, we present a magnetic-activated cell sorting (MACS) method where podocytes are isolated from digested whole kidneys using antibodies specific to extracellular antigens on podocytes. Using microbeaded secondary antibodies binding to the podocyte-specific primary antibodies allows sorting of the podocytes using a magnet. This podocyte-only cell fraction is a unique source of in vivo–derived cells for molecular and cellular experiments.

The complexity of the human transcriptome poses significant challenges for complete annotation. Traditional RNA-seq, often limited by sensitivity and short read lengths, is frequently inadequate for identifying low-abundant transcripts and resolving complex populations of transcript isoforms. Direct long-read sequencing, while offering full-length information, suffers from throughput limitations, hindering the capture of low-abundance transcripts. To address these challenges, we introduce a targeted RNA enrichment strategy, rapid amplification of cDNA ends coupled with Nanopore sequencing (RACE-Nano-Seq). This method unravels the deep complexity of transcripts containing anchor sequences—specific regions of interest that might be exons of annotated genes, in silico predicted exons, or other sequences. RACE-Nano-Seq is based on inverse PCR with primers targeting these anchor regions to enrich the corresponding transcripts in both 5' and 3' directions. This method can be scaled for high-throughput transcriptome profiling by using multiplexing strategies. Through targeted RNA enrichment and full-length sequencing, RACE-Nano-Seq enables accurate and comprehensive profiling of low-abundance transcripts, often revealing complex transcript profiles at the targeted loci, both annotated and unannotated.

Single-cell RNA sequencing has revolutionized molecular cell biology by enabling the identification of unique transcription profiles and cell transcription states within the same tissue. However, tissue dissociation presents a challenge for non-model organisms, as commercial kits are often incompatible, and current protocols rely on tissue enzymatic digestion for extended periods. Tissue digestion can alter cell transcription in response to temperature and the stress caused by enzymatic treatment. Here, we propose a protocol to stabilize RNA using a deep eutectic solvent (Vivophix, Rapid Labs) prior to tissue dissociation, thereby avoiding transcription changes induced by the process and preventing RNase activity during incubation. We validated this methodology for three medically important insect vectors: Anopheles gambiae, Aedes aegypti, and Lutzomyia longipalpis. Single-cell RNA sequencing using our insect midgut dissociation protocol yielded high-quality sequencing results, with a high number of cells recovered, a low percentage of mitochondrial reads, and a low percentage of ambient RNA—two hallmark standards of cell quality.

Capturing produced, consumed, or exchanged metabolites (metabolomics) and the result of gene expression (transcriptomics) require the extraction of metabolites and RNA. Multi-omics approaches and, notably, the combination of metabolomics and transcriptomic analyses are required for understanding the functional changes and adaptation of microorganisms to different physico-chemical and environmental conditions. A protocol was developed to extract total RNA and metabolites from less than 6 mg of a kind of phototrophic biofilm: oxygenic photogranules. These granules are aggregates of several hundred micrometers up to several millimeters. They harbor heterotrophic bacteria and phototrophs. After a common step for cell disruption by bead-beating, a part of the volume was recovered for RNA extraction, and the other half was used for the methanol- and dichloromethane-based extraction of metabolites. The solvents enabled the separation of two phases (aqueous and lipid) containing hydrophilic and lipophilic metabolites, respectively. The ¹H nuclear magnetic resonance (NMR) analysis of these extracts produced spectra that contained over a hundred signals with a signal-to-noise ratio higher than 10. The quality of the spectra enabled the identification of dozens of metabolites per sample. Total RNA was purified using a commercially available kit, yielding sufficient concentration and quality for metatranscriptomic analysis. This novel method enables the co-extraction of RNA and metabolites from the same sample, as opposed to the parallel extraction from two samples. Using the same sample for both extractions is particularly advantageous when working with inherently heterogeneous complex biofilm. In heterogeneous systems, differences between samples may be substantial. The co-extraction will enable a holistic analysis of the metabolomics and metatranscriptomics data generated, minimizing experimental biases, including technical variations and, notably, biological variability. As a result, it will ensure more robust multi-omics analyses, particularly by improving the correlation between metabolic changes and transcript modifications.

Bone repair is a complex regenerative process relying on skeletal stem/progenitor cells (SSPCs) recruited predominantly from the periosteum. Activation and differentiation of periosteal SSPCs occur in a heterogeneous environment, raising the need for single cell/nucleus transcriptomics to decipher the response of the periosteum to injury. Enzymatic cell dissociation can induce a stress response affecting the transcriptome and lead to overrepresentation of certain cell types (i.e., immune and endothelial cells) and low coverage of other cell types of interest. To counteract these limitations, we optimized a protocol to isolate nuclei directly from the intact periosteum and from the fracture callus to perform single-nucleus RNA sequencing. This protocol is adapted for fresh murine periosteum, fracture callus, and frozen human periosteum. Nuclei are isolated using mechanical extraction combined with fluorescence-based nuclei sorting to obtain high-quality nucleus suspensions. This protocol allows the capture of the full diversity of cell types in the periosteum and fracture environment to better reflect the in vivo tissue composition.

Dual RNA-Seq technology has significantly advanced the study of biological interactions between two organisms by allowing parallel transcriptomic analysis. Existing analysis methods employ various combinations of open-source bioinformatics tools to process dual RNA-Seq data. Upon reviewing these methods, we intend to explore crucial criteria for selecting standard tools and methods, especially focusing on critical steps such as trimming and mapping reads to the reference genome. In order to validate the different combinatorial approaches, we performed benchmarking using top-ranking tools and a publicly available dual RNA-Seq Sequence Read Archive (SRA) dataset. An important observation while evaluating the mapping approach is that when the adapter trimmed reads are first mapped to the pathogen genome, more reads align to the pathogen genome than the unmapped reads derived from the traditional host-first mapping approach. This mapping method prevents the misalignment of pathogen reads to the host genome due to their shorter length. In this way, the pathogenic read information found at lesser proportions in a complex eukaryotic dataset is precisely obtained. This protocol presents a comprehensive comparison of these possible approaches, resulting in a robust unified standard methodology.

Information on RNA localisation is essential for understanding physiological and pathological processes, such as gene expression, cell reprogramming, host–pathogen interactions, and signalling pathways involving RNA transactions at the level of membrane-less or membrane-bounded organelles and extracellular vesicles. In many cases, it is important to assess the topology of RNA localisation, i.e., to distinguish the transcripts encapsulated within an organelle of interest from those merely attached to its surface. This allows establishing which RNAs can, in principle, engage in local molecular interactions and which are prevented from interacting by membranes or other physical barriers. The most widely used techniques interrogating RNA localisation topology are based on the treatment of isolated organelles with RNases with subsequent identification of the surviving transcripts by northern blotting, qRT-PCR, or RNA-seq. However, this approach produces incoherent results and many false positives. Here, we describe Controlled Level of Contamination coupled to deep sequencing (CoLoC-seq), a more refined subcellular transcriptomics approach that overcomes these pitfalls. CoLoC-seq starts by the purification of organelles of interest. They are then either left intact or lysed and subjected to a gradient of RNase concentrations to produce unique RNA degradation dynamics profiles, which can be monitored by northern blotting or RNA-seq. Through straightforward mathematical modelling, CoLoC-seq distinguishes true membrane-enveloped transcripts from degradable and non-degradable contaminants of any abundance. The method has been implemented in the mitochondria of HEK293 cells, where it outperformed alternative subcellular transcriptomics approaches. It is applicable to other membrane-bounded organelles, e.g., plastids, single-membrane organelles of the vesicular system, extracellular vesicles, or viral particles.

Key features

• Tested on human mitochondria; potentially applicable to cell cultures, non-model organisms, extracellular vesicles, enveloped viruses, tissues; does not require genetic manipulations or highly pure organelles.

• In the case of human cells, the required amount of starting material is ~2,500 cm² of 80% confluent cells (or ~3 × 10⁸ HEK293 cells).

• CoLoC-seq implements a special RNA-seq strategy to selectively capture intact transcripts, which requires RNases generating 5′-hydroxyl and 2′/3′-phosphate termini (e.g., RNase A, RNase I).

• Relies on nonlinear regression software with customisable exponential functions.

Graphical overview

Accessible chromatin regions modulate gene expression by acting as cis-regulatory elements. Understanding the epigenetic landscape by mapping accessible regions of DNA is therefore imperative to decipher mechanisms of gene regulation under specific biological contexts of interest. The assay for transposase-accessible chromatin sequencing (ATAC-seq) has been widely used to detect accessible chromatin and the recent introduction of single-cell technology has increased resolution to the single-cell level. In a recent study, we used droplet-based, single-cell ATAC-seq technology (scATAC-seq) to reveal the epigenetic profile of the transit-amplifying subset of thymic epithelial cells (TECs), which was identified previously using single-cell RNA-sequencing technology (scRNA-seq). This protocol allows the preparation of nuclei from TECs in order to perform droplet-based scATAC-seq and its integrative analysis with scRNA-seq data obtained from the same cell population. Integrative analysis has the advantage of identifying cell types in scATAC-seq data based on cell cluster annotations in scRNA-seq analysis.

Due to overlapping sequences with linear cognates, identifying internal sequences of circular RNA (circRNA) remains a challenge. Recently, we have developed a full-length circRNA sequencing method (circFL-seq) and computational pipeline, to profile ordinary and fusion circRNA at the isoform level. Compared to short-read RNA-seq, rolling circular reverse transcription and nanopore long-read sequencing of circFL-seq make circRNA reads more than tenfold enriched, and show advantages for identification of both short (<100 nt) and long (>2,000 nt) circRNA transcripts. circFL-seq allows identification of differential alternative splicing suggested wide application prospects for functional studies of internal sequences in circRNAs. In addition, the experimental protocol and computational pipeline of circFL-seq shows better sensitivity and precision for identification of back-splicing junctions than current long-read sequencing methods. Together, the accurate identification and quantification of full-length circRNAs makes circFL-seq a potential tool for large-scale screening of functional circRNAs.

In neurons, local translation in dendritic and axonal compartments allows for the fast and on-demand modification of the local proteome. As the last few years have witnessed dramatic advancements in our appreciation of the brain’s neuronal diversity, it is increasingly relevant to understand how local translation is regulated according to cell type. To this end, both sequencing-based and imaging-based techniques have recently been reported. Here, we present a subcellular single cell RNA sequencing protocol that allows molecular quantification from the soma and dendrites of single neurons, and which can be scaled up for the characterization of several hundreds to thousands of neurons. Somata and dendrites of cultured neurons are dissected using laser capture microdissection, followed by cell lysis to release mRNA content. Reverse transcription is then conducted using an indexed primer that allows the downstream pooling of samples. The pooled cDNA library is prepared for and sequenced in an Illumina platform. Finally, the data generated are processed and converted into a gene vs. cells digital expression table. This protocol provides detailed instructions for both wet lab and bioinformatic steps, as well as insights into controls, data analysis, interpretations, and ways to achieve robust and reproducible results.

Graphic abstract:

Subcellular Single Cell RNA-seq in Neurons.