Bioinformatics and Computational Biology

Real-time quantitative PCR (qPCR) is a pivotal technique for analyzing gene expression and DNA copy number variations. However, the limited availability of user-friendly software tools for qPCR data analysis presents a significant challenge for experimental biologists with limited computational skills. To address this issue, we developed Click-qPCR, a user-friendly and web-based Shiny application for qPCR data analysis. Click-qPCR streamlines ΔCq and ΔΔCq calculations using user-uploaded CSV data files. The interactive interface of the application allows users to select genes and experimental groups and perform Welch’s t tests and one-way analysis of variance with Dunnett’s post-hoc test for pairwise and multi-group comparisons, respectively. Results are visualized via interactive bar plots (mean ± standard deviation with individual data points) and can be downloaded as publication-quality images, along with summary statistics. Click-qPCR empowers researchers to efficiently process, interpret, and visualize qPCR data regardless of their programming experience, thereby facilitating routine analysis tasks. Click-qPCR Shiny application is available at https://kubo-azu.shinyapps.io/Click-qPCR/, while its source code and user guide are available at https://github.com/kubo-azu/Click-qPCR.

In neuropharmacology and drug development, in silico methods have become increasingly vital, particularly for studying receptor–ligand interactions at the molecular level. Membrane proteins such as GABA (A) receptors play a central role in neuronal signaling and are key targets for therapeutic intervention. While experimental techniques like electrophysiology and radioligand binding provide valuable functional data, they often fall short in resolving the structural complexity of membrane proteins and can be time-consuming, costly, and inaccessible in many research settings. This study presents a comprehensive computational workflow for investigating membrane protein–ligand interactions, demonstrated using the GABA (A) receptor α5β2γ2 subtype and mitragynine, an alkaloid from Mitragyna speciosa (Kratom), as a case study. The protocol includes homology modeling of the receptor based on a high-resolution template, followed by structure optimization and validation. Ligand docking is then used to predict binding sites and affinities at known modulatory interfaces. Finally, molecular dynamics (MD) simulations assess the stability and conformational dynamics of receptor–ligand complexes over time. Overall, this workflow offers a robust, reproducible approach for structural analysis of membrane protein–ligand interactions, supporting early-stage drug discovery and mechanistic studies across diverse membrane protein targets.

Insects rely on chemosensory proteins, including gustatory receptors, to detect chemical cues that regulate feeding, mating, and oviposition behaviours. Conventional approaches for studying these proteins are limited by the scarcity of experimentally resolved structures, especially in non-model pest species. Here, we present a reproducible computational protocol for the identification, functional annotation, and structural modelling of insect chemosensory proteins, demonstrated using gustatory receptors from the red palm weevil (Rhynchophorus ferrugineus) as an example. The protocol integrates publicly available sequence data with OmicsBox for functional annotation and ColabFold for high-confidence structure prediction, providing a step-by-step framework that can be applied to genome-derived or transcriptomic datasets. The workflow is designed for broad applicability across insect species and generates structurally reliable protein models suitable for downstream applications such as ligand docking or molecular dynamics simulations. By bridging functional annotation with structural characterisation, this protocol enables reproducible studies of chemosensory proteins in agricultural and ecological contexts and supports the development of novel pest management strategies.

The rhizosphere, a 2–10 mm region surrounding the root surface, is colonized by numerous microorganisms, known as the rhizosphere microbiome. These microorganisms interact with each other, leading to emergent properties that affect plant fitness. Mapping these interactions is crucial to understanding microbial ecology in the rhizosphere and predicting and manipulating plant health. However, current methods do not capture the chemistry of the rhizosphere environment, and common plant–microbe interaction study setups do not map bacterial interactions in this niche. Additionally, studying bacterial interactions may require the creation of transgenic bacterial lines with markers for antibiotic resistance/fluorescent probes and even isotope labeling. Here, we describe a protocol for both in silico prediction and in vitro validation of bacterial interactions that closely recapitulate the major chemical constituents of the rhizosphere environment using a widely used Murashige & Skoog (MS)-based gnotobiotic plant growth system. We use the auto-fluorescent Pseudomonas, abundantly found in the rhizosphere, to estimate their interactions with other strains, thereby avoiding the need for the creation of transgenic bacterial strains. By combining artificial root exudate medium, plant cultivation medium, and a synthetic bacterial community (SynCom), we first simulate their interactions using genome-scale metabolic models (GSMMs) and then validate these interactions in vitro, using growth assays. We show that the GSMM-predicted interaction scores correlate moderately, yet significantly, with their in vitro validation. Given the complexity of interactions among rhizosphere microbiome members, this reproducible and efficient protocol will allow confident mapping of interactions of fluorescent Pseudomonas with other bacterial strains within the rhizosphere microbiome.

DNA methylation is a crucial epigenetic modification that influences gene expression and plays a role in various biological processes. High-throughput sequencing techniques, such as bisulfite sequencing (BS-seq) and enzymatic methyl sequencing (EM-seq), enable genome-wide profiling of DNA methylation patterns with single-base resolution. In this protocol, we present a bioinformatics pipeline for analyzing genome-wide DNA methylation. We outline the step-by-step process of the essential analyses, including quality control using FASTQ for BS- and EM-seqs raw reads, read alignment with commonly used aligners such as Bowtie2 and BS-Seeker2, DNA methylation calling to generate CGmap files, identification of differentially methylated regions (DMRs) using tools including MethylC-analyzer and HOME, data visualization, and post-alignment analyses. Compared to existing workflows, this pipeline integrates multiple steps into a single protocol, lowering the technical barrier, improving reproducibility, and offering flexibility for both plant and animal methylome studies. To illustrate the application of BS-seq and EM-seq, we demonstrate a case study on analyzing a mutant in Arabidopsis thaliana with a mutation in the met1 gene, which encodes a DNA methyltransferase, and results in global CG hypomethylation and altered gene regulation. This example highlights the biological insights that can be gained through systematic methylome analysis. Our workflow is adaptable to any organism with a reference genome and provides a robust framework for uncovering methylation-associated regulatory mechanisms. All scripts and detailed instructions are provided in GitHub repository: https://github.com/PaoyangLab/Methylation_Analysis.

High-content analysis (HCA) is a powerful image-based approach for phenotypic profiling and drug discovery, enabling the extraction of multiparametric data from individual cells. Traditional HCA protocols often rely on fixed-cell imaging, with assays like cell painting widely adopted as standard. While these methods provide rich morphological information, the integration of live-cell imaging expands analytical capabilities by enabling the study of dynamic biological processes and real-time cellular responses. This protocol presents a simple, cost-effective, and scalable method for live-cell HCA using acridine orange (AO), a metachromatic fluorescent dye that highlights cellular organization by staining nucleic acids and acidic compartments. The assay provides visualization of distinct subcellular structures, including nuclei and cytoplasmic organelles, using a two-channel fluorescence readout. Compatible with high-throughput microscopy and computational analysis, the method supports diverse applications such as phenotypic screening, cytotoxicity assessment, and morphological profiling. By preserving cell viability and enabling dynamic, real-time measurements, this live-cell imaging approach complements existing fixed-cell assays and offers a versatile platform for uncovering complex cellular phenotypes.

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.

Chromatin-associated RNAs (caRNAs) have been increasingly recognized as key regulators of gene expression and genome architecture. A few technologies, such as ChRD-PET and RedChIP, have emerged to assess protein-mediated RNA–chromatin interactions, but each has limitations. Here, we describe the TaDRIM-seq (targeted DNA-associated RNA and RNA–RNA interaction mapping by sequencing) technique, which combines Protein G (PG)-Tn5-targeted DNA tagmentation with in situ proximity ligation to simultaneously profile caRNAs across genomic regions and capture global RNA–RNA interactions within intact nuclei. This approach reduces the required cell input, shortens the experimental duration compared to existing protocols, and is applicable to both mammalian and plant systems.

Thousands of RNAs are localized to specific subcellular locations, and these localization patterns are often required for optimal cell function. However, the sequences within RNAs that direct their transport are unknown for almost all localized transcripts. Similarly, the RNA content of most subcellular locations remains unknown. To facilitate the study of subcellular transcriptomes, we developed the RNA proximity labeling method OINC-seq. OINC-seq utilizes photoactivatable, spatially restricted RNA oxidation to specifically label RNA in proximity to a subcellularly localized bait protein. After labeling, these oxidative RNA marks are then read out via high-throughput sequencing due to their ability to induce predictable misincorporation events by reverse transcriptase. These induced mutations are then quantitatively assessed for each gene using our software package PIGPEN. The observed mutation rate for a given RNA species is therefore related to its proximity to the localized bait protein. This protocol describes procedures for assaying RNA localization via OINC-seq experiments as well as computational procedures for analyzing the resulting data using PIGPEN.

Transcriptional pausing dynamically regulates spatiotemporal gene expression during cellular differentiation, development, and environmental adaptation. Precise measurement of pausing duration, a critical parameter in transcriptional control, has been challenging due to limitations in resolution and confounding factors. We introduce Fast TV-PRO-seq, an optimized protocol built on time-variant precision run-on sequencing (TV-PRO-seq), which enables genome-wide, single-base resolution mapping of RNA polymerase II pausing times. Unlike standard PRO-seq, Fast TV-PRO-seq employs sarkosyl-free biotin-NTP run-on with time gradients and integrates on-bead enzymatic reactions to streamline workflows. Key improvements include (1) reducing experimental time from 4 to 2 days, (2) reducing cell input requirements, and (3) improved process efficiency and simplified command-line operations through the use of bash scripts.