High-resolution transcriptomic and epigenetic profiling identifies novel regulators of COPD.

Patients with chronic obstructive pulmonary disease (COPD) are still waiting for curative treatments. Considering its environmental cause, we hypothesized that COPD will be associated with altered epigenetic signaling in lung cells. We generated genome-wide DNA methylation maps at single CpG resolution of primary human lung fibroblasts (HLFs) across COPD stages. We show that the epigenetic landscape is changed early in COPD, with DNA methylation changes occurring predominantly in regulatory regions. RNA sequencing of matched fibroblasts demonstrated dysregulation of genes involved in proliferation, DNA repair, and extracellular matrix organization. Data integration identified 110 candidate regulators of disease phenotypes that were linked to fibroblast repair processes using phenotypic screens. Our study provides high-resolution multi-omic maps of HLFs across COPD stages. We reveal novel transcriptomic and epigenetic signatures associated with COPD onset and progression and identify new candidate regulators involved in the pathogenesis of chronic lung diseases. The presence of various epigenetic factors among the candidates demonstrates that epigenetic regulation in COPD is an exciting research field that holds promise for novel therapeutic avenues for patients.

Whole-genome bisulfite sequencing data analysis learning module on Google Cloud Platform.

This study describes the development of a resource module that is part of a learning platform named 'NIGMS Sandbox for Cloud-based Learning' https://github.com/NIGMS/NIGMS-Sandbox. The overall genesis of the Sandbox is described in the editorial NIGMS Sandbox at the beginning of this Supplement. This module is designed to facilitate interactive learning of whole-genome bisulfite sequencing (WGBS) data analysis utilizing cloud-based tools in Google Cloud Platform, such as Cloud Storage, Vertex AI notebooks and Google Batch. WGBS is a powerful technique that can provide comprehensive insights into DNA methylation patterns at single cytosine resolution, essential for understanding epigenetic regulation across the genome. The designed learning module first provides step-by-step tutorials that guide learners through two main stages of WGBS data analysis, preprocessing and the identification of differentially methylated regions. And then, it provides a streamlined workflow and demonstrates how to effectively use it for large datasets given the power of cloud infrastructure. The integration of these interconnected submodules progressively deepens the user's understanding of the WGBS analysis process along with the use of cloud resources. Through this module, we can enhance the accessibility and adoption of cloud computing in epigenomic research, speeding up the advancements in the related field and beyond. This manuscript describes the development of a resource module that is part of a learning platform named ``NIGMS Sandbox for Cloud-based Learning'' https://github.com/NIGMS/NIGMS-Sandbox. The overall genesis of the Sandbox is described in the editorial NIGMS Sandbox [1] at the beginning of this Supplement. This module delivers learning materials on the analysis of bulk and single-cell ATAC-seq data in an interactive format that uses appropriate cloud resources for data access and analyses.

CRL4DCAF13 E3 ubiquitin ligase targets MeCP2 for degradation to prevent DNA hypermethylation and ensure normal transcription in growing oocytes.

The DNA methylation is gradually acquired during oogenesis, a process sustained by successful follicle development. However, the functional roles of methyl-CpG-binding protein 2 (MeCP2), an epigenetic regulator displaying specifical binding with methylated DNA, remains unknown in oogenesis. In this study, we found MeCP2 protein was highly expressed in primordial and primary follicle, but was almost undetectable in secondary follicles. However, in aged ovary, MeCP2 protein is significantly increased in both oocyte and granulosa cells. Overexpression of MeCP2 in growing oocyte caused transcription dysregulation, DNA hypermethylation, and genome instability, ultimately leading to follicle growth arrest and apoptosis. MeCP2 is targeted by DCAF13, a substrate recognition adaptor of the Cullin 4-RING (CRL4) E3 ligase, and polyubiquitinated for degradation in both cells and oocytes. Dcaf13-null oocyte exhibited an accumulation of MeCP2 protein, and the partial rescue of follicle growth arrest induced by Dcaf13 deletion was observed following MeCP2 knockdown. The RNA-seq results revealed that large amounts of genes were regulated by the DCAF13-MeCP2 axis in growing oocytes. Our study demonstrated that CRL4DCAF13 E3 ubiquitin ligase targets MeCP2 for degradation to ensure normal DNA methylome and transcription in growing oocytes. Moreover, in aged ovarian follicles, deceased DCAF13 and DDB1 protein were observed, indicating a potential novel mechanism that regulates ovary aging.

Characterisation and reproducibility of the HumanMethylationEPIC v2.0 BeadChip for DNA methylation profiling.

BACKGROUND: The Illumina family of Infinium Methylation BeadChip microarrays has been widely used over the last 15 years for genome-wide DNA methylation profiling, including large-scale and population-based studies, due to their ease of use and cost effectiveness. Succeeding the popular HumanMethylationEPIC BeadChip (EPICv1), the recently released Infinium MethylationEPIC v2.0 BeadChip (EPICv2) claims to extend genomic coverage to more than 935,000 CpG sites. Here, we comprehensively characterise the reproducibility, reliability and annotation of the EPICv2 array, based on bioinformatic analysis of both manifest data and new EPICv2 data from diverse biological samples. RESULTS: We find a high degree of reproducibility with EPICv1, evidenced by comparable sensitivity and precision from empirical cross-platform comparison incorporating whole genome bisulphite sequencing (WGBS), and high correlation between technical sample replicates, including between samples with DNA input levels below the manufacturer's recommendation. We provide a full assessment of probe content, evaluating genomic distribution and changes from previous array versions. We characterise EPICv2's new feature of replicated probes and provide recommendations as to the superior probes. In silico analysis of probe sequences demonstrates that probe cross-hybridisation remains a significant problem in EPICv2. By mapping the off-target sites at single nucleotide resolution and comparing with WGBS we show empirical evidence for preferential off-target binding. CONCLUSIONS: Overall, we find EPICv2 a worthy successor to the previous Infinium methylation microarrays, however some technical issues remain. To support optimal EPICv2 data analysis we provide an expanded version of the EPICv2 manifest to aid researchers in understanding probe design, data processing, choosing appropriate probes for analysis and for integration with methylation datasets from previous versions of the Infinium Methylation BeadChip.

Integrated analysis of the DNA methylome and RNA transcriptome during the development of skeletal muscle in Duroc pigs.

BACKGROUND: Skeletal muscle development plays a crucial role in yield and quality of pork; however, this process is influenced by various factors. In this study, we employed whole-genome bisulfite sequencing (WGBS) and transcriptome sequencing to comprehensively investigate the longissimus dorsi muscle (LDM), aiming to identify key genes that impact the growth and development of Duroc pigs with different average daily gains (ADGs). RESULTS: Eight pigs were selected and divided into two groups based on ADGs: H (774.89 g) group and L (658.77 g) group. Each pair of the H and L groups were half-siblings. The results of methylation sequencing revealed 2631 differentially methylated genes (DMGs) involved in metabolic processes, signalling, insulin secretion, and other biological activities. Furthermore, a joint analysis was conducted on these DMGs and the differentially expressed genes (DEGs) obtained from transcriptome sequencing of the same individual. This analysis identified 316 differentially methylated and differentially expressed genes (DMEGs), including 18 DMEGs in promoter regions and 294 DMEGs in gene body regions. Finally, LPAR1 and MEF2C were selected as candidate genes associated with muscle development. Bisulfite sequencing PCR (BSP) and quantitative real-time PCR (qRT-PCR) revealed that the promoter region of LPAR1 exhibited significantly lower methylation levels (P < 0.05) and greater expression levels (P < 0.05) in the H group than in the L group. Additionally, hypermethylation was observed in the gene body region of MEF2C, as was a low expression level, in the H group (P < 0.05). CONCLUSIONS: These results suggest that the differences in the ADGs of Duroc pigs fed the same diet may be influenced by the methylation levels and expression levels of genes related to skeletal muscle development.

Comparing methylation levels assayed in GC-rich regions with current and emerging methods.

DNA methylation is an epigenetic mechanism that regulates gene expression, and for mammals typically occurs on cytosines within CpG dinucleotides. A significant challenge for methylation detection methods is accurately measuring methylation levels within GC-rich regions such as gene promoters, as inaccuracies compromise downstream biological interpretation of the data. To address this challenge, we compared methylation levels assayed using four different Methods Enzymatic Methyl-seq (EM-seq), whole genome bisulphite sequencing (WGBS), Infinium arrays (Illumina MethylationEPIC, "EPIC"), and Oxford Nanopore Technologies nanopore sequencing (ONT) applied to human DNA. Overall, all methods produced comparable and consistent methylation readouts across the human genome. The flexibility offered by current gold standard WGBS in interrogating genome-wide cytosines is surpassed technically by both EM-seq and ONT, as their coverages and methylation readouts are less prone to GC bias. These advantages are tempered by increased laboratory time (EM-seq) and higher complexity (ONT). We further assess the strengths and weaknesses of each method, and provide recommendations in choosing the most appropriate methylation method for specific scientific questions or translational needs.

Single generation epigenetic change in captivity and reinforcement in subsequent generations in a delta smelt (Hypomesus transpacificus) conservation hatchery.

A refugial population of the endangered delta smelt (Hypomesus transpacificus) has been maintained at the Fish Conservation and Culture Laboratory (FCCL) at UC Davis since 2008. Despite intense genetic management, fitness differences between wild and cultured fish have been observed at the FCCL. To investigate the molecular underpinnings of hatchery domestication, we used whole-genome bisulfite sequencing to quantify epigenetic differences between wild and hatchery-origin delta smelt. Differentially methylated regions (DMRs) were identified from 104 individuals by comparing the methylation patterns in different generations of hatchery fish (G1, G2, G3) with their wild parents (G0). We discovered a total of 132 significant DMRs (p < .05) between G0 and G1, 132 significant DMRs between G0 and G2, and 201 significant DMRs between G0 and G3. Our results demonstrate substantial differences in methylation patterns emerged between the wild and hatchery-reared fish in the early generations in the hatchery, with a higher proportion of hypermethylated DMRs in hatchery-reared fish. The rearing environment was found to be a stronger predictor of individual clustering based on methylation patterns than family, sex or generation. Our study indicates a reinforcement of the epigenetic status with successive generations in the hatchery environment, as evidenced by an increase in methylation in hypermethylated DMRs and a decrease in methylation in hypomethylated DMRs over time. Lastly, our results demonstrated heterogeneity in inherited methylation pattern in families across generations. These insights highlight the long-term consequences of hatchery practices on the epigenetic landscape, potentially impacting wild fish populations.

A Strategy for Studying Epigenetic Diversity in Natural Populations: Proof of Concept in Poplar and Oak.

These last 20 years, several techniques have been developed for quantifying DNA methylation, the most studied epigenetic marks in eukaryotes, including the gold standard method, whole-genome bisulphite sequencing (WGBS). WGBS quantifies genome-wide DNA methylation but has several inconveniences rendering it less suitable for population-scale epigenetic studies. The high cost of deep sequencing and the large amounts of data generated prompted us to seek an alternative approach. Restricting studies to parts of the genome would be a satisfactory alternative had there not been a major limitation: the need to select upstream targets corresponding to differentially methylated regions (DMRs) as targets. Given the need to study large numbers of samples, we propose a strategy for investigating DNA methylation variation in natural populations, considering the structural complexity of the genomes with their size and their content in unique as coding regions versus repeated regions as transposable elements. We first identified regions of highly variable DNA methylation in a representative subset of genotypes representative of the biological diversity in the population by WGBS. We then analysed the variations of DNA methylation in these targeted regions at the population level by Sequencing Capture Bisulphite (SeqCapBis). The entire strategy was then validated by applying it to another species. Our strategy was developed as a proof of concept on natural populations of two forest species: Populus nigra and Quercus petraea.

High-resolution methylome analysis uncovers stress-responsive genomic hotspots and drought-sensitive TE superfamilies in the clonal Lombardy poplar.

DNA methylation is environment-sensitive and can mediate stress responses. In long-lived trees, changing environments might cumulatively shape the methylome landscape over their lifetime. However, because high-resolution methylome studies usually focus on single environmental cues, it remains unclear to what extent the methylation responses are generic or stress-specific, and how this relates to their long-term stability. Here, we studied the methylome plasticity of a Populus nigra cv. 'Italica' clone that is widespread across Europe. Adult trees from a variety of geographic locations were clonally propagated in a common garden experiment, and the ramets were exposed to cold, heat, drought, herbivory, rust infection, and salicylic acid treatments. Through comprehensive whole-genome bisulfite sequencing, we analyzed stress-induced and naturally occurring DNA methylation variants. Stress-induced methylation changes predominantly targeted transposable elements. When occurring in CG/CHG contexts, the same regions were often affected by multiple stresses, suggesting a generic response of the methylome. Drought stress caused a distinct CHH hypermethylation response in transposable elements, affecting entire TE superfamilies near drought-responsive genes. Methylation differences in CG/CHG contexts that were induced by stress treatments showed striking overlap with methylation differences observed between trees from distinct geographical locations. Thus, we revealed genomic hotspots of methylation change that are not stress-specific and that contribute to natural DNA methylation variation, and we identified specific transposable element superfamilies that respond to a specific stress with possible functional consequences. Our results underscore the importance of studying the effects of multiple stressors in a single experiment for recognizing general versus stress-specific methylome responses.

Vpr driving DNA methylation variation of CD4 + T cells in HIV-1 infection.

BACKGROUND: Despite the existence of available therapeutic interventions for HIV-1, this virus remains a significant global threat, leading to substantial morbidity and mortality. Within HIV-1-infected cells, the accessory viral protein r (Vpr) exerts control over diverse biological processes, including cell cycle progression, DNA repair, and apoptosis. The regulation of gene expression through DNA methylation plays a crucial role in physiological processes, exerting its influence without altering the underlying DNA sequence. However, a thorough examination of the impact of Vpr on DNA methylation in human CD4 + T cells has not been conducted. METHODS: In this study, we employed base-resolution whole-genome bisulfite sequencing (WGBS), real-time quantitative RCR and western blot to explore the effect of Vpr on DNA methylation of host cells under HIV-1 infection. RESULTS: We observed that HIV-1 infection leads to elevated levels of global DNA methylation in primary CD4 + T cells. Specifically, Vpr induces significant modifications in DNA methylation patterns, particularly affecting regions within promoters and gene bodies. These alterations notably influence genes related to immune-related pathways and olfactory receptor activity. Moreover, Vpr demonstrates a distinct ability to diminish the levels of methylation in histone genes. CONCLUSIONS: These findings emphasize the significant involvement of Vpr in regulating transcription through the modulation of DNA methylation patterns. Together, the results of this investigation will considerably enhance our understanding of the influence of HIV-1 Vpr on the DNA methylation of host cells, offer potential avenues for the development of more effective treatments.

Genome-wide DNA methylation changes in Oryzias melastigma embryos exposed to the water accommodated fraction of crude oil.

The water accommodated fraction (WAF) of crude oil exerts considerable impacts on marine fish during embryonic stage. Clarifying changes in epigenetic modifications is helpful for understanding the molecular mechanism underlying the toxicity of embryonic WAF exposure. The aim of this study was to explore genome-wide DNA methylation changes in Oryzias melastigma embryos after exposure to the nominal total petroleum hydrocarbon concentration of 500 µg/L in WAF for 7 days. Whole-genome bisulfite sequencing revealed that 8.47 % and 8.46 % of all the genomic C sites were methylated in the control and WAF-exposed groups, respectively. Among the three sequence contexts, methylated CG site had the largest number in both the two groups. The sequence preferences of nearby methylated cytosines were consistent between the two groups. A total of 4798 differentially methylated regions (DMRs) were identified in the promoter region. Furthermore, Gene Ontology analysis revealed that DMR-related genes were enriched mainly for functions related to development and nervous system. Additionally, the Kyoto Encyclopedia of Genes and Genomes pathways enriched in DMR-related genes were related to nervous system and endocrine system. These novel findings provide comprehensive insights into the genome-wide DNA methylation landscape of O. melastigma following embryonic WAF exposure, shedding light on the epigenetic regulatory mechanisms underlying WAF-induced toxicity.

Epigenome-wide analysis of aging effects on liver regeneration.

BACKGROUND: Aging is known to exert an effect on liver regeneration, with the ability of liver to regenerate displaying a significant decline over time. Liver physiological parameters such as liver volume, blood flow, and metabolism, as well as the ability to regenerate after injury have all been shown to decrease at old age in humans and model systems, with a number of molecular mechanisms proposed to be involved, including DNA methylation-dependent genome remodeling. To address how changes in DNA methylation mediate the adverse aging effect on liver regeneration, we searched for differentially methylated genomic regions (DMRs) in mouse livers co-regulated by aging and regeneration and determined their associated genes and enriched pathways. RESULTS: DMRs were identified using whole-genome bisulfite sequencing (WGBS). Pathway analysis of aging DMR-mapped genes revealed two distinct phases of aging, 2-to-8 and 8-to-16 months old (m/o). Regenerative DMR-mapped differentially expressed genes (DEGs) were enriched in pathways controlling cell proliferation and differentiation. Most DMRs shared by both aging and regeneration changed in the same methylation direction between 2 and 8 m/o but in the opposite direction between 8 and 16 m/o. Regenerative DMRs inversely affected by aging during 8-to-16 m/o were found in the promoter/gene regions of 12 genes. Four regenerative DEGs were synchronously regulated by early aging and inversely regulated by mid-to-late aging DMRs. Lead DMR-mapped genes were validated by their expression profiles in liver aging and regeneration. CONCLUSIONS: Our study has uncovered new DMRs and gene targets inversely affected by liver aging and regeneration to explain the adverse aging effect on liver regeneration. These findings will be of fundamental importance to understand the epigenomic changes underlying the biology of aging on liver regeneration.

Intestinal epigenotype of Atlantic salmon (Salmo salar) associates with tenacibaculosis and gut microbiota composition.

It remains a challenge to obtain the desired phenotypic traits in aquacultural production of Atlantic salmon, and part of the challenge might come from the effect that host-associated microorganisms have on the fish phenotype. To manipulate the microbiota towards the desired host traits, it is critical to understand the factors that shape it. The bacterial gut microbiota composition can vary greatly among fish, even when reared in the same closed system. While such microbiota differences can be linked to diseases, the molecular effect of disease on host-microbiota interactions and the potential involvement of epigenetic factors remain largely unknown. The aim of this study was to investigate the DNA methylation differences associated with a tenacibaculosis outbreak and microbiota displacement in the gut of Atlantic salmon. Using Whole Genome Bisulfite Sequencing (WGBS) of distal gut tissue from 20 salmon, we compared the genome-wide DNA methylation levels between uninfected individuals and sick fish suffering from tenacibaculosis and microbiota displacement. We discovered >19,000 differentially methylated cytosine sites, often located in differentially methylated regions, and aggregated around genes. The 68 genes connected to the most significant regions had functions related to the ulcerous disease such as epor and slc48a1a but also included prkcda and LOC106590732 whose orthologs are linked to microbiota changes in other species. Although the expression level was not analysed, our epigenetic analysis suggests specific genes potentially involved in host-microbiota interactions and more broadly it highlights the value of considering epigenetic factors in efforts to manipulate the microbiota of farmed fish.

Transcriptomic and Epigenomic Responses to Cortisol-Mediated Stress in Rainbow Trout (Oncorhynchus mykiss) Skeletal Muscle.

The production and release of cortisol during stress responses are key regulators of growth in teleosts. Understanding the molecular responses to cortisol is crucial for the sustainable farming of rainbow trout (Oncorhynchus mykiss) and other salmonid species. While several studies have explored the genomic and non-genomic impacts of cortisol on fish growth and skeletal muscle development, the long-term effects driven by epigenetic mechanisms, such as cortisol-induced DNA methylation, remain unexplored. In this study, we analyzed the transcriptome and genome-wide DNA methylation in the skeletal muscle of rainbow trout seven days after cortisol administration. We identified 550 differentially expressed genes (DEGs) by RNA-seq and 9059 differentially methylated genes (DMGs) via whole-genome bisulfite sequencing (WGBS) analysis. KEGG enrichment analysis showed that cortisol modulates the differential expression of genes associated with nucleotide metabolism, ECM-receptor interaction, and the regulation of actin cytoskeleton pathways. Similarly, cortisol induced the differential methylation of genes associated with focal adhesion, adrenergic signaling in cardiomyocytes, and Wnt signaling. Through integrative analyses, we determined that 126 genes showed a negative correlation between up-regulated expression and down-regulated methylation. KEGG enrichment analysis of these genes indicated participation in ECM-receptor interaction, regulation of actin cytoskeleton, and focal adhesion. Using RT-qPCR, we confirmed the differential expression of lamb3, itga6, limk2, itgb4, capn2, and thbs1. This study revealed for the first time the molecular responses of skeletal muscle to cortisol at the transcriptomic and whole-genome DNA methylation levels in rainbow trout.

Grafting based DNA methylation alteration of snoRNAs in upland cotton (Gossypium L.).

The effects of grafting in response to various biotic and abiotic stressors have been studied, however, the methylation status of small nucleolar RNA (snoRNA) genes in heterograft and homograft cotton needs investigation. This study was undertaken to determine grafting effects on DNA methylation of snoRNA genes in Upland cotton. Rootstocks used were Pima 3-79 (Gossypium barbadense acc. Pima 3-79) and Texas Marker-1 (G. hirsutum acc. TM-1), representing two different species with different fiber properties, adaptations, and morphologies. The methylation ratio and differently methylated cytosines (DMCs) of 10935 snoRNA genes in mature seeds of heterograft and homograft cotton samples were studied using the whole genome bisulfite sequencing method. Seedling vigor and seed weight were studied to investigate phenotype alterations that might be associated with altered methylation levels among grafts. Statistically significant DMC differences among gene elements of snoRNA genes and between homograft and heterograft cotton samples were identified in the absence of DNA sequence alterations. DNA methylation alterations of snoRNA genes associated with seedling vigor and 100 seed weight. The majority of snoRNA genes showed higher numbers of mCG + mCHG-DMCs with increased methylation levels in heterograft, while there were higher numbers of mCG + mCHG-DMCs with decreased methylation levels in homograft. Since snoRNAs regulate essential genes for plant growth and development and plant adaptation to different habitats or extreme environments, their altered methylation levels should be related with plant physiology.

BiocMAP: a Bioconductor-friendly, GPU-accelerated pipeline for bisulfite-sequencing data.

BACKGROUND: Bisulfite sequencing is a powerful tool for profiling genomic methylation, an epigenetic modification critical in the understanding of cancer, psychiatric disorders, and many other conditions. Raw data generated by whole genome bisulfite sequencing (WGBS) requires several computational steps before it is ready for statistical analysis, and particular care is required to process data in a timely and memory-efficient manner. Alignment to a reference genome is one of the most computationally demanding steps in a WGBS workflow, taking several hours or even days with commonly used WGBS-specific alignment software. This naturally motivates the creation of computational workflows that can utilize GPU-based alignment software to greatly speed up the bottleneck step. In addition, WGBS produces raw data that is large and often unwieldy; a lack of memory-efficient representation of data by existing pipelines renders WGBS impractical or impossible to many researchers. RESULTS: We present BiocMAP, a Bioconductor-friendly methylation analysis pipeline consisting of two modules, to address the above concerns. The first module performs computationally-intensive read alignment using Arioc, a GPU-accelerated short-read aligner. Since GPUs are not always available on the same computing environments where traditional CPU-based analyses are convenient, the second module may be run in a GPU-free environment. This module extracts and merges DNA methylation proportions-the fractions of methylated cytosines across all cells in a sample at a given genomic site. Bioconductor-based output objects in R utilize an on-disk data representation to drastically reduce required main memory and make WGBS projects computationally feasible to more researchers. CONCLUSIONS: BiocMAP is implemented using Nextflow and available at http://research.libd.org/BiocMAP/ . To enable reproducible analysis across a variety of typical computing environments, BiocMAP can be containerized with Docker or Singularity, and executed locally or with the SLURM or SGE scheduling engines. By providing Bioconductor objects, BiocMAP's output can be integrated with powerful analytical open source software for analyzing methylation data.

LiBis: an ultrasensitive alignment augmentation for low-input bisulfite sequencing.

The cell-free DNA (cfDNA) methylation profile in liquid biopsy has been utilized to diagnose early-stage disease and estimate therapy response. However, typical clinical procedures are capable of purifying only very small amounts of cfDNA. Whole-genome bisulfite sequencing (WGBS) is the gold standard for measuring DNA methylation; however, WGBS using small amounts of fragmented DNA introduces a critical challenge for data analysis, namely a low-mapping ratio. The resulting low sequencing depth and low coverage of CpG sites genome-wide is a bottleneck for the clinical application of cfDNA-based WGBS assays. We developed LiBis (Low-input Bisulfite Sequencing), a novel method for low-input WGBS data alignment. By dynamically clipping initially unmapped reads and remapping clipped fragments, we judiciously rescued those reads and uniquely aligned them to the genome. By substantially increasing the mapping ratio by up to 88%, LiBis dramatically improved the number of informative CpGs and the precision in quantifying the methylation status of individual CpG sites. LiBis significantly improved the cost efficiency of low-input WGBS experiments by dynamically removing contamination introduced by random priming. The high sensitivity and cost effectiveness afforded by LiBis for low-input samples will allow the discovery of genetic and epigenetic features suitable for downstream analysis and biomarker identification using liquid biopsy.

Comprehensive benchmarking of software for mapping whole genome bisulfite data: from read alignment to DNA methylation analysis.

Whole genome bisulfite sequencing is currently at the forefront of epigenetic analysis, facilitating the nucleotide-level resolution of 5-methylcytosine (5mC) on a genome-wide scale. Specialized software have been developed to accommodate the unique difficulties in aligning such sequencing reads to a given reference, building on the knowledge acquired from model organisms such as human, or Arabidopsis thaliana. As the field of epigenetics expands its purview to non-model plant species, new challenges arise which bring into question the suitability of previously established tools. Herein, nine short-read aligners are evaluated: Bismark, BS-Seeker2, BSMAP, BWA-meth, ERNE-BS5, GEM3, GSNAP, Last and segemehl. Precision-recall of simulated alignments, in comparison to real sequencing data obtained from three natural accessions, reveals on-balance that BWA-meth and BSMAP are able to make the best use of the data during mapping. The influence of difficult-to-map regions, characterized by deviations in sequencing depth over repeat annotations, is evaluated in terms of the mean absolute deviation of the resulting methylation calls in comparison to a realistic methylome. Downstream methylation analysis is responsive to the handling of multi-mapping reads relative to mapping quality (MAPQ), and potentially susceptible to bias arising from the increased sequence complexity of densely methylated reads.

Pig fetal skeletal muscle development is associated with genome-wide DNA hypomethylation and corresponding alterations in transcript and microRNA expression.

Fetal myogenesis represents a critical period of porcine skeletal muscle development and requires coordinated expression of thousands of genes. Epigenetic mechanisms, including DNA methylation, drive transcriptional regulation during development; however, these processes are understudied in developing porcine tissues. We performed bisulfite sequencing to assess DNA methylation in pig longissimus dorsi muscle at 41- and 70-days gestation (dg), as well as RNA- and small RNA-sequencing to identify coordinated changes in methylation and expression between myogenic stages. We identified 45 739 differentially methylated regions (DMRs) between stages, and the majority (N = 34 232) were hypomethylated at 70 versus 41 dg. Integration of methylation and transcriptomic data revealed strong associations between differential gene methylation and expression. Differential miRNA methylation was significantly negatively correlated with abundance, and dynamic expression of assayed miRNAs persisted postnatally. Motif analysis revealed significant enrichment of myogenic regulatory factor motifs among hypomethylated regions, suggesting that DNA hypomethylation may function to increase accessibility of muscle-specific transcription factors. We show that developmental DMRs are enriched for GWAS SNPs for muscle- and meat-related traits, demonstrating the potential for epigenetic processes to influence phenotypic diversity. Our results enhance understanding of DNA methylation dynamics of porcine myogenesis and reveal putative cis-regulatory elements governed by epigenetic processes.

Regulation of DNA methylation during the testicular development of Shaziling pigs.

DNA methylation is one of the key epigenetic regulatory mechanisms in development and spermatogenesis. However, the dynamic regulatory mechanisms of genome-wide DNA methylation during testicular development remain largely unknown. Herein, we generated a single-base resolution DNA methylome and transcriptome atlas of precocious porcine testicular tissues across three developmental stages (1, 75, and 150 days old). The results showed that the dynamic methylation patterns were directly related to the expression of the DNMT3A gene. Conjoint analysis revealed a negative regulatory pattern between promoter methylation and the positive regulation of 3'-untranslated region (3'UTR) methylation. Mechanistically, the decrease in promoter methylation affected the upregulation of meiosis-related genes, such as HORMAD1, SPO11, and SYCE1. Demethylation in the 3'UTR induced the downregulation of the INHBA gene and knockdown of INHBA inhibited cell proliferation by reducing the synthesis of activin A. These findings contribute to exploring the regulatory mechanisms of DNA methylation in testicular development.