系统生物学

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.

Trypanosoma cruzi, the causative agent of Chagas disease, faces significant metabolic challenges due to fluctuating nutrient availability and oxidative stress within its insect vector. Metabolomic techniques, such as gas chromatography–mass spectrometry (GC–MS), have been widely used to study the adaptive mechanisms of the parasite. This article describes a standardized method for the untargeted metabolomics analysis of T. cruzi epimastigote, covering parasite cultivation, sample deproteinization with methanol, metabolite extraction, derivatization with BSTFA, and GC–MS analysis. To ensure robustness and reproducibility, statistical analysis uses univariate tests, as well as multivariate approaches such as principal component analysis (PCA) and partial least squares (PLS) regression. The protocol offers a reliable and sensitive method to study metabolic responses in T. cruzi under environmental stress, with low biological variability and high reproducibility.

This protocol provides a step-by-step approach for generating single-gene knockout in hard-to-transfect suspension immune cell lines like THP1, specifically demonstrated by knocking out the GSDMD gene. By employing CRISPR-Cas9 system delivered via lentivirus, this protocol enables precise gene disruption through targeted single-guide RNAs (sgRNAs). Key steps include designing specific sgRNAs, cloning them into a CRISPR vector, viral packaging, and transducing the target cells, followed by selection and validation. This optimized protocol is particularly useful for functional studies in immune cells, allowing researchers to reliably explore gene function in complex cellular pathways.

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.

N⁶-methyladenosine (m6A) is an abundant internal mRNA modification with roles in regulating cellular and organismal physiology, including development, differentiation, and disease. The deposition of m6A is highly regulated, with various m6A levels across different environmental conditions, cellular states, and cell types. Available methods for measuring bulk m6A levels are often time-consuming, have low throughput, and/or require specialized instrumentation or data analyses. Here, we present a detailed protocol for measuring bulk m6A levels in purified poly(A) RNA samples with m6A-ELISA using a standard-based approach. Critical steps of the protocol are highlighted and optimized, including poly(A) RNA quality controls and antibody specificity testing. The protocol is fast, scalable, adaptable, and cost-effective. It does not require specialized instrumentation, training, or skills in data analysis. We have successfully tested this protocol on mRNAs isolated from budding yeast and mouse cell lines.

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.

With reduced genotyping costs, genome-wide association studies (GWAS) face more challenges in diverse populations with complex structures to map genes of interest. The complex structure demands sophisticated statistical models, and increased marker density and population size require efficient computing tools. Many statistical models and computing tools have been developed with varied properties in statistical power, computing efficiency, and user-friendly accessibility. Some statistical models were developed with dedicated computing tools, such as efficient mixed model analysis (EMMA), multiple loci mixed model (MLMM), fixed and random model circulating probability unification (FarmCPU), and Bayesian-information and linkage-disequilibrium iteratively nested keyway (BLINK). However, there are computing tools (e.g., GAPIT) that implement multiple statistical models, retain a constant user interface, and maintain enhancement on input data and result interpretation. In this study, we developed a protocol utilizing a minimal set of software tools (BEAGLE, BLINK, and GAPIT) to perform a variety of analyses including file format conversion, missing genotype imputation, GWAS, and interpretation of input data and outcome results. We demonstrated the protocol by reanalyzing data from the Rice 3000 Genomes Project and highlighting advancements in GWAS model development.

Bayesian phylogenetic analysis is essential for elucidating evolutionary relationships among organisms. Traditional methods often rely on fixed models and manual parameter settings, which can limit accuracy and efficiency. This protocol presents an integrated workflow that leverages GUIDANCE2 for rigorous sequence alignment, ProtTest and MrModeltest for robust model selection, and MrBayes for phylogenetic tree estimation through Bayesian inference. By automating key steps and providing detailed command-line instructions, this protocol enhances the reliability and reproducibility of phylogenetic studies.

Zebrafish genetic mutants have emerged as a valuable model system for studying various aspects of disease and developmental biology. Mutant zebrafish embryos are generally identified based on phenotypic defects at later developmental stages, making it difficult to investigate underlying molecular mechanisms at earlier stages. This protocol presents a PCR-based genotyping method that enables the identification of wild-type, heterozygous, and homozygous zebrafish genetic mutants at any developmental stage, even when they are phenotypically indistinguishable. The approach involves the amplification of specific genomic regions using carefully designed primers, followed by gel electrophoresis. This genotyping method facilitates the investigation of the molecular mechanisms driving phenotypic defects that are observed at later timepoints. This protocol allows researchers to perform analyses such as immunofluorescence, RT-PCR, RNA sequencing, and other molecular experiments on early developmental stages of mutants. The availability of this protocol expands the utility of zebrafish genetic mutants for elucidating the molecular underpinnings of various biological processes throughout development.

This manuscript details two modified protocols for the isolation of long-stranded or high molecular weight (HMW) DNA from Magnaporthaceae (Ascomycota) fungal mycelium intended for whole genome sequencing. The Cytiva Nucleon PhytoPure and the Macherey-Nagel NucleoBond HMW DNA kits were selected because the former requires lower amounts of starting material and the latter utilizes gentler methods to maximize DNA length, albeit at a higher requirement for input material. The Cytiva Nucleon PhytoPure kit successfully recovered HMW DNA for half of our fungal species by increasing the amount of RNase A treatment and adding in a proteinase K treatment. To reduce the impact of pigmentation development, which occurs toward later stages of culturing, extractions were run in quadruplicate to increase overall DNA concentration. We also adapted the Macherey-Nagel NucleoBond HMW DNA kit for high-quality HMW DNA by grinding the sample to a fine powder, overnight lysis, and splitting the sample before washing the precipitated DNA. For both kits, precipitated DNA was spooled out pre-washing, ensuring a higher percentage of high-integrity long strands. The Macherey-Nagel protocol offers advantages over the first through the utilization of gravity columns that provide gentler treatment, yielding >50% of high-purity DNA strands exceeding 40 kbp. The limitation of this method is the requirement for a large quantity of starting material (1 g). By triaging samples based on the rate of growth relative to the accumulation of secondary metabolites, our methodologies hold promise for yielding reliable and high-quality HMW DNA from a variety of fungal samples, improving sequencing outcomes.