系统生物学

DNA methylation is a fundamental epigenetic mark with critical roles in epigenetic regulation, development, and genome stability across diverse organisms. Whole genome bisulfite sequencing (WGBS) enables single-base resolution mapping of cytosine methylation patterns and has become a standard method in epigenomics. This protocol provides a detailed, step-by-step workflow for WGBS library construction starting from genomic DNA. It includes steps of RNaseA treatment, DNA shearing, end-repair and A-tailing, adapter ligation, bisulfite conversion, library amplification, and quantification. Notably, the method uses self-prepared reagents and customizable index systems, avoiding the constraints of commercial library preparation kits. This flexibility supports cost-effective, scalable methylome profiling, suitable for diverse experimental designs, including high-throughput multiplexed sequencing.

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.

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.

Chromatin accessibility is a key determinant of gene expression that can be altered under different physiological and disease conditions. Skeletal muscle is made up of myofibers that are highly plastic and adaptive. Therefore, assessing the genome-wide chromatin state of myofibers under various conditions is very important to gain insight into the epigenetic state of myonuclei. The rigid nature of myofibers, as well as the low number of myonuclei that they contain, have rendered genome-wide studies with myofibers challenging. In recent years, ATAC-Seq from whole muscle and single nucleus ATAC-Seq have been performed. However, these techniques cannot distinguish between different fiber and cell types present in the muscle. In addition, due to the limited depth capacity obtained from single nucleus ATAC-Seq, an extensive comparative analysis cannot be performed. Here, we introduce a protocol where we combine the isolation of a single myofiber with OMNI ATAC-Seq. This protocol allows for genome-wide analysis of accessible chromatin regions of a selected single myofiber at a sufficient depth for comparative analysis under various physiological and disease conditions. This protocol can also allow for a specific myofiber to be selected, such as a regenerating myofiber. In the future, this protocol can help identify global changes in chromatin state under various conditions, as well as between different types of myofibers.

Graphical abstract:

DNA methylation is a conserved chemical modification, by which methyl groups are added to the cytosine of DNA molecules. Methylation can influence gene expression without changing the sequence of a particular gene. This epigenetic effect is an intriguing phenomenon that has puzzled biologists for years. By probing the temporal and spatial patterns of DNA methylation in genomes, it is possible to learn about the biological role of cytosine methylation, as well as its involvement in gene regulation and transposon silencing. Advances in whole-genome sequencing have led to the widespread adoption of methods that examine genome-wide patterns of DNA methylation. Achieving sufficient sequencing depth in these types of experiments is costly, particularly for pilot studies in organisms with large genome sizes, or incomplete reference genomes. To overcome this issue, assays to determine site-specific DNA methylation can be used. Although often used, these assays are rarely described in detail. Here, we describe a pipeline that applies traditional TA cloning, Sanger sequencing, and online tools to examine DNA methylation. We provide an example of how to use this protocol to examine the pattern of DNA methylation at a specific transposable element in maize.

In recent years, DNA methylation research has been accelerated by the advent of nanopore sequencers. However, read length has been limited by the constraints of base conversion using the bisulfite method, making analysis of chromatin content difficult. The read length of the previous method combining bisulfite conversion and long-read sequencing was ~1.5 kb, even using targeted PCR. In this study, we have improved read length (~5 kb), by converting unmethylated cytosines to uracils with APOBEC enzymes, to reduce DNA fragmentation. The converted DNA was then sequenced using a PromethION nanopore sequencer. We have also developed a new analysis pipeline that accounts for base conversions, which are not present in conventional nanopore sequencing, as well as errors produced by nanopore sequencing.

Characterizing the molecular mechanisms regulating gene expression is crucial for understanding the regulatory processes underlying physiological responses to environmental and developmental signals in eukaryotes. The covalent modification of histones contributes to the compaction levels of chromatin, as well as the recruitment of the transcriptional machinery to specific loci, facilitating metastable changes in gene activity. ChIP-seq (Chromatin Immunoprecipitation followed by sequencing) has become the gold standard method for determining histone modification profiles among different organisms, tissues, and genotypes. In the current protocol, we describe a highly robust method for performing ChIP-seq of histone modifications in Arabidopsis thaliana plantlets. Besides its robustness, this method uses in-house-prepared buffers for chromatin extraction, immunoprecipitation, washing, and elusion, making it cost-effective in contrast to commercial kits.

In the field of chromatin biology, a major goal of understanding the roles of histone post-translational modifications is to identify the proteins and domains that recognize these modifications. Synthetic histone peptides containing one or more modifications are a key tool to probe these interactions in pull-down assays with recombinant proteins or cell lysates. Building on these approaches, the binding specificity of a protein of interest can be screened against many histone peptides in parallel using a peptide array. In this protocol, we describe the expression and purification of a recombinant protein of interest in bacteria, followed by an assay for binding to histone post-translational modifications using a commercially available histone peptide array. The purification uses a versatile dual-tagging and cleavage strategy and equipment commonly available in a molecular biology laboratory.

Graphic abstract:

Overview of protocol for purifying recombinant protein and hybridizing to a histone peptide array.

Primary somatosensory neurons, whose cell bodies reside in the dorsal root ganglion (DRG) and trigeminal ganglion, are specialized to transmit sensory information from the periphery to the central nervous system. Our molecular understanding of peripheral sensory neurons has been limited by both their heterogeneity and low abundance compared with non-neuronal cell types in sensory ganglia. We describe a protocol to isolate nuclei from mouse DRGs using iodixanol density gradient centrifugation, which enriches for neuronal nuclei while still sampling non-neuronal cells such as satellite glia and Schwann cells. This protocol is compatible with a range of downstream applications such as single-nucleus transcriptional and epigenomic assays.

DNA methylation in gene promoters plays a major role in gene expression regulation, and alterations in methylation patterns have been associated with several diseases. In this context, different software suites and statistical methods have been proposed to analyze differentially methylated positions and regions. Among them, the novel statistical method implemented in the mCSEA R package proposed a new framework to detect subtle, but consistent, methylation differences. Here, we provide an easy-to-use pipeline covering all the necessary steps to detect differentially methylated promoters with mCSEA from Illumina 450K and EPIC methylation BeadChips data. This protocol covers the download of data from public repositories, quality control, data filtering and normalization, estimation of cell type proportions, and statistical analysis. In addition, we show the procedure to compare disease vs. normal phenotypes, obtaining differentially methylated regions including promoters or CpG Islands. The entire protocol is based on R programming language, which can be used in any operating system and does not require advanced programming skills.