Systems Biology

Functional enrichment analysis is essential for understanding the biological significance of differentially expressed genes. Commonly used tools such as g:Profiler, DAVID, and GOrilla are effective when applied to well-annotated model organisms. However, for non-model organisms, particularly for bacteria and other microorganisms, curated functional annotations are often scarce. In such cases, researchers often rely on homology-based approaches, using tools like BLAST to transfer annotations from closely related species. Although this strategy can yield some insights, it often introduces annotation errors and overlooks unique species-specific functions. To address this limitation, we present a user-friendly and adaptable method for creating custom annotation R packages using genomic data retrieved from NCBI. These packages can be directly imported as libraries into the R environment and are compatible with the clusterProfiler package, enabling effective gene ontology and pathway enrichment analysis. We demonstrate this approach by constructing an R annotation package for Mycobacterium tuberculosis H37Rv, as an example. The annotation package is then utilized to analyze differentially expressed genes from a subset of RNA-seq dataset (GSE292409), which investigates the transcriptional response of M. tuberculosis H37Rv to rifampicin treatment. The chosen dataset includes six samples, with three serving as untreated controls and three exposed to rifampicin for 1 h. Further, enrichment analysis was performed on genes to demonstrate changes in response to the treatment. This workflow provides a reliable and scalable solution for functional enrichment analysis in organisms with limited annotation resources. It also enhances the accuracy and biological relevance of gene expression interpretation in microbial genomics research.

Pathogen gene expression from host samples is often challenging to study due to low signal and high host RNA background. PCR probes have been recently used to hybridize and extract bacterial sequences from next-generation sequencing (NGS) libraries generated from in vitro and animal models of infection; however, these strategies require purchasing commercially synthesized probes that often do not capture the entire transcriptome. Transcript-capture sequencing is a novel capture approach for extracting RNA of a target bacterial species from samples in which there is substantial contamination by the host or other microbes. Biotinylated 150-base-pair DNA probes are generated in-house from bacterial DNA spanning the entire bacterial genome. Probes are hybridized to the cDNA of NGS sequencing libraries prepared from host samples to capture and enrich for bacterial-specific RNA reads before sequencing. This method results in a >200-fold increase in bacterial RNA reads from infected host samples (including in vitro, animal, and human samples) and generates complete bacterial transcriptomes with high gene coverage (>80%). Use of this protocol on infected host samples reveals a snapshot of bacterial activity during disease that may improve understanding of the physiological state of pathogens within their hosts.

RNA alternative splicing (AS) is an essential process that expands transcriptomic and proteomic diversity in eukaryotic cells and contributes to cellular heterogeneity across physiological and pathological conditions in humans. With the advent of single-cell RNA sequencing (scRNA-seq), it has become possible to study AS at cellular resolution, although robust and standardized analytical workflows remain to be developed. Here, we present a stepwise protocol for analyzing AS in single cells from pediatric high-grade gliomas (pHGGs) harboring the histone H3.3 lysine 27-to-methionine (H3.3K27M) mutation using SMART-Seq2 scRNA-seq data. Starting from raw sequencing reads, the workflow includes read alignment, gene-level quantification, splice junction and intron quantification, and single-nucleotide variant-based mutation detection. Gene expression–based clustering and cell-type annotation are performed by using the Seurat R package. AS analysis in tumor cells is then conducted using the MARVEL R package in combination with customized scripts to calculate percent spliced-in (PSI) values, identify variable AS events, perform dimensionality reduction, cluster cells, conduct differential AS analysis, and visualize splicing patterns. This protocol provides a reproducible and comprehensive framework for dissecting AS dynamics at single-cell resolution. It is readily adaptable to other SMART-Seq2 datasets and facilitates systematic investigation of splicing heterogeneity in diverse biological contexts.

The canonical role of transfer RNAs (tRNAs) in protein synthesis has been extensively characterized; however, recent studies have uncovered novel functions for tRNA as a mediator of long-distance signaling in plants. Several studies have identified dicistronic tRNA-mRNA transcripts that contain a tRNA gene and an adjacent protein-coding gene (PCG) that are transcribed as a single unit. These transcripts are associated with RNA systemic mobility through the plant’s vascular tissues, potentially acting as non-cell-autonomous signaling messengers in coordinating development and stress responses. Here, we report a computational pipeline to detect dicistronic tRNA-mRNA transcripts from short-read next-generation RNA-sequencing datasets; to our knowledge, this is the only established pipeline for the systematic identification of such candidates in plants. The dicistronic RNA transcript version 2 (v2) described here improves on the earlier version DiRT v1 by expanding the repertoire of dicistronic transcripts detected to include tRNA-like structures (TLS) as well as functional tRNAs, which were already supported in the pipeline. The updated protocol also includes detection of dicistronic tRNA or TLS sequences within genomic features such as untranslated regions (UTRs). The accurate detection of both tRNAs and UTR-embedded tRNA-like sequences (TLS) is critical, as these RNA structures have been reported to function as mediators of long-distance RNA mobility. Furthermore, as NGS datasets are prone to sequencing artifacts and potential DNA contamination, we improved the pipeline’s statistical robustness by including read coverage of flanking intronic regions as a baseline control. To account for potential DNA contamination during RNA-seq library preparation, detected tRNA-mRNA transcripts are deemed as putatively dicistronic only if the coverage of their intergenic region is significantly higher (Student’s t-test, FDR < 0.05) than flanking intronic regions. Furthermore, the updated pipeline allows this statistical test to be applied to intronless and single-intron genes. Using this updated protocol, we identified novel tRNA and TLS dicistronic transcripts in both grapevine (Vitis spp. Ruggeri 140) and Arabidopsis thaliana datasets and validated in vitro using RT-PCR. We provide a fast and reliable method to detect dicistronic transcripts that can be applied to any short-read RNA-sequencing dataset, fast-tracking the functional characterization of these newly emerging transcripts.

While cell hashing enhances single-cell RNA sequencing (scRNA-seq) efficiency and minimizes batch effects, commercial mouse hashtags often fail in FVB/N and several other strains due to antibody-epitope incompatibility. We describe a robust alternative utilizing biotinylated antibody cocktails and streptavidin-conjugated oligos to enable reliable sample multiplexing. This approach was validated in FVB/N lung tissues, yielding high-quality single-cell libraries. Our protocol offers a practical solution for researchers requiring strain-specific or custom-designed multiplexing strategies for single-cell transcriptomics.

Single-cell RNA sequencing (scRNA-seq) is a powerful technique for exploring cellular heterogeneity and host–pathogen interactions. This protocol details the Zika virus (ZIKV)-targeted scRNA-seq workflow for preparing high-quality single-cell suspensions from the whole brain tissues of neonatal mice, high-quality single-cell sorting, cDNA reverse transcription, amplification, ZIKV enrichment and host transcriptome library preparation, and sequencing dataset integration in downstream analysis to complete the quantification of ZIKV RNA in individual cells.

Weighted gene co-expression network analysis (WGCNA) is widely used in transcriptomic studies to identify groups of highly correlated genes, aiding in the understanding of disease mechanisms. Although numerous protocols exist for constructing WGCNA networks from gene expression data, many focus on single datasets and do not address how to compare module stability across conditions. Here, we present a protocol for constructing and comparing WGCNA modules in paired tumor and normal datasets, enabling the identification of modules involved in both core biological processes and those specifically related to cancer pathogenesis. By incorporating module preservation analysis, this approach allows researchers to gain deeper insights into the molecular underpinnings of oral cancer, as well as other diseases. Overall, this protocol provides a framework for module preservation analysis in paired datasets, enabling researchers to identify which gene co-expression modules are conserved or disrupted between conditions, thereby advancing our understanding of disease-specific vs. universal biological processes.

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.

Formalin-fixed paraffin-embedded (FFPE) samples remain an underutilized resource in single-cell omics due to RNA degradation from formalin fixation. Here, we present snPATHO-seq, a robust and adaptable approach that enables the generation of high-quality single-nucleus (sn) transcriptomic data from FFPE tissues, utilizing advancements in single-cell genomic techniques. The snPATHO-seq workflow integrates optimized nuclei isolation with the 10× Genomics Flex assay, targeting short RNA fragments to mitigate FFPE-related RNA degradation. Benchmarking against standard 10× 3' and Flex assays for fresh/frozen tissues confirmed robust detection of transcriptomic signatures and cell types. snPATHO-seq demonstrated high performance across diverse FFPE samples, including diseased tissues like breast cancer. It seamlessly integrates with FFPE spatial transcriptomics (e.g., FFPE Visium) for multi-modal spatial and single-nucleus profiling. Compared to workflows like 10× Genomics’ snFFPE, snPATHO-seq delivers superior data quality by reducing tissue debris and preserving RNA integrity via nuclei isolation. This cost-effective workflow enables high-resolution transcriptomics of archival FFPE samples, advancing single-cell omics in translational and clinical research.