The Chinese pine genome and methylome unveil key features of conifer evolution.

Conifers dominate the world's forest ecosystems and are the most widely planted tree species. Their giant and complex genomes present great challenges for assembling a complete reference genome for evolutionary and genomic studies. We present a 25.4-Gb chromosome-level assembly of Chinese pine (Pinus tabuliformis) and revealed that its genome size is mostly attributable to huge intergenic regions and long introns with high transposable element (TE) content. Large genes with long introns exhibited higher expressions levels. Despite a lack of recent whole-genome duplication, 91.2% of genes were duplicated through dispersed duplication, and expanded gene families are mainly related to stress responses, which may underpin conifers' adaptation, particularly in cold and/or arid conditions. The reproductive regulation network is distinct compared with angiosperms. Slow removal of TEs with high-level methylation may have contributed to genomic expansion. This study provides insights into conifer evolution and resources for advancing research on conifer adaptation and development.

Reconstructing Denisovan Anatomy Using DNA Methylation Maps.

Denisovans are an extinct group of humans whose morphology remains unknown. Here, we present a method for reconstructing skeletal morphology using DNA methylation patterns. Our method is based on linking unidirectional methylation changes to loss-of-function phenotypes. We tested performance by reconstructing Neanderthal and chimpanzee skeletal morphologies and obtained >85% precision in identifying divergent traits. We then applied this method to the Denisovan and offer a putative morphological profile. We suggest that Denisovans likely shared with Neanderthals traits such as an elongated face and a wide pelvis. We also identify Denisovan-derived changes, such as an increased dental arch and lateral cranial expansion. Our predictions match the only morphologically informative Denisovan bone to date, as well as the Xuchang skull, which was suggested by some to be a Denisovan. We conclude that DNA methylation can be used to reconstruct anatomical features, including some that do not survive in the fossil record.

Cell-of-Origin Patterns Dominate the Molecular Classification of 10,000 Tumors from 33 Types of Cancer.

We conducted comprehensive integrative molecular analyses of the complete set of tumors in The Cancer Genome Atlas (TCGA), consisting of approximately 10,000 specimens and representing 33 types of cancer. We performed molecular clustering using data on chromosome-arm-level aneuploidy, DNA hypermethylation, mRNA, and miRNA expression levels and reverse-phase protein arrays, of which all, except for aneuploidy, revealed clustering primarily organized by histology, tissue type, or anatomic origin. The influence of cell type was evident in DNA-methylation-based clustering, even after excluding sites with known preexisting tissue-type-specific methylation. Integrative clustering further emphasized the dominant role of cell-of-origin patterns. Molecular similarities among histologically or anatomically related cancer types provide a basis for focused pan-cancer analyses, such as pan-gastrointestinal, pan-gynecological, pan-kidney, and pan-squamous cancers, and those related by stemness features, which in turn may inform strategies for future therapeutic development.

Genome-wide DNA methylation changes in human spermatogenesis.

Sperm production and function require the correct establishment of DNA methylation patterns in the germline. Here, we examined the genome-wide DNA methylation changes during human spermatogenesis and its alterations in disturbed spermatogenesis. We found that spermatogenesis is associated with remodeling of the methylome, comprising a global decline in DNA methylation in primary spermatocytes followed by selective remethylation, resulting in a spermatids/sperm-specific methylome. Hypomethylated regions in spermatids/sperm were enriched in specific transcription factor binding sites for DMRT and SOX family members and spermatid-specific genes. Intriguingly, while SINEs displayed differential methylation throughout spermatogenesis, LINEs appeared to be protected from changes in DNA methylation. In disturbed spermatogenesis, germ cells exhibited considerable DNA methylation changes, which were significantly enriched at transposable elements and genes involved in spermatogenesis. We detected hypomethylation in SVA and L1HS in disturbed spermatogenesis, suggesting an association between the abnormal programming of these regions and failure of germ cells progressing beyond meiosis.

BioM2: biologically informed multi-stage machine learning for phenotype prediction using omics data.

Navigating the complex landscape of high-dimensional omics data with machine learning models presents a significant challenge. The integration of biological domain knowledge into these models has shown promise in creating more meaningful stratifications of predictor variables, leading to algorithms that are both more accurate and generalizable. However, the wider availability of machine learning tools capable of incorporating such biological knowledge remains limited. Addressing this gap, we introduce BioM2, a novel R package designed for biologically informed multistage machine learning. BioM2 uniquely leverages biological information to effectively stratify and aggregate high-dimensional biological data in the context of machine learning. Demonstrating its utility with genome-wide DNA methylation and transcriptome-wide gene expression data, BioM2 has shown to enhance predictive performance, surpassing traditional machine learning models that operate without the integration of biological knowledge. A key feature of BioM2 is its ability to rank predictor variables within biological categories, specifically Gene Ontology pathways. This functionality not only aids in the interpretability of the results but also enables a subsequent modular network analysis of these variables, shedding light on the intricate systems-level biology underpinning the predictive outcome. We have proposed a biologically informed multistage machine learning framework termed BioM2 for phenotype prediction based on omics data. BioM2 has been incorporated into the BioM2 CRAN package (https://cran.r-project.org/web/packages/BioM2/index.html).

Paternal aging impacts expression and epigenetic markers as early as the first embryonic tissue lineage differentiation.

BACKGROUND: Advanced paternal age (APA) is associated with adverse outcomes to offspring health, including increased risk for neurodevelopmental disorders. The aim of this study was to investigate the methylome and transcriptome of the first two early embryonic tissue lineages, the inner cell mass (ICM) and the trophectoderm (TE), from human blastocysts in association with paternal age and disease risk. High quality human blastocysts were donated with patient consent from donor oocyte IVF cycles from either APA (≥ 50 years) or young fathers. Blastocysts were mechanically separated into ICM and TE lineage samples for both methylome and transcriptome analyses. RESULTS: Significant differential methylation and transcription was observed concurrently in ICM and TE lineages of APA-derived blastocysts compared to those from young fathers. The methylome revealed significant enrichment for neuronal signaling pathways, as well as an association with neurodevelopmental disorders and imprinted genes, largely overlapping within both the ICM and TE lineages. Significant enrichment of neurodevelopmental signaling pathways was also observed for differentially expressed genes, but only in the ICM. In stark contrast, no significant signaling pathways or gene ontology terms were identified in the trophectoderm. Despite normal semen parameters in aged fathers, these significant molecular alterations can adversely contribute to downstream impacts on offspring health, in particular neurodevelopmental disorders like autism spectrum disorder and schizophrenia. CONCLUSIONS: An increased risk for neurodevelopmental disorders is well described in children conceived by aged fathers. Using blastocysts derived from donor oocyte IVF cycles to strategically control for maternal age, our data reveals evidence of methylation dysregulation in both tissue lineages, as well as transcription dysregulation in neurodevelopmental signaling pathways associated with APA fathers. This data also reveals that embryos derived from APA fathers do not appear to be compromised for initial implantation potential with no significant pathway signaling disruption in trophectoderm transcription. Collectively, our work provides insights into the complex molecular mechanisms that occur upon paternal aging during the first lineage differentiation in the preimplantation embryo. Early expression and epigenetic markers of APA-derived preimplantation embryos highlight the susceptibility of the future fetus to adverse health outcomes.

Genetic and epigenetic alterations in precursor lesions of endometrial endometrioid carcinoma.

The hyperplasia-carcinoma sequence is a stepwise tumourigenic programme towards endometrial cancer in which normal endometrial epithelium becomes neoplastic through non-atypical endometrial hyperplasia (NAEH) and atypical endometrial hyperplasia (AEH), under the influence of unopposed oestrogen. NAEH and AEH are known to exhibit polyclonal and monoclonal cell growth, respectively; yet, aside from focal PTEN protein loss, the genetic and epigenetic alterations that occur during the cellular transition remain largely unknown. We sought to explore the potential molecular mechanisms that promote the NAEH-AEH transition and identify molecular markers that could help to differentiate between these two states. We conducted target-panel sequencing on the coding exons of 596 genes, including 96 endometrial cancer driver genes, and DNA methylome microarrays for 48 NAEH and 44 AEH lesions that were separately collected via macro- or micro-dissection from the endometrial tissues of 30 cases. Sequencing analyses revealed acquisition of the PTEN mutation and the clonal expansion of tumour cells in AEH samples. Further, across the transition, alterations to the DNA methylome were characterised by hypermethylation of promoter/enhancer regions and CpG islands, as well as hypo- and hyper-methylation of DNA-binding regions for transcription factors relevant to endometrial cell differentiation and/or tumourigenesis, including FOXA2, SOX17, and HAND2. The identified DNA methylation signature distinguishing NAEH and AEH lesions was reproducible in a validation cohort with modest discriminative capability. These findings not only support the concept that the transition from NAEH to AEH is an essential step within neoplastic cell transformation of endometrial epithelium but also provide deep insight into the molecular mechanism of the tumourigenic programme. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Single-base resolution quantitative genome methylation analysis in the model bacterium Helicobacter pylori by enzymatic methyl sequencing (EM-Seq) reveals influence of strain, growth phase, and methyl homeostasis.

BACKGROUND: Bacterial epigenetics is a rapidly expanding research field. DNA methylation by diverse bacterial methyltransferases (MTases) contributes to genomic integrity and replication, and many recent studies extended MTase function also to global transcript regulation and phenotypic variation. Helicobacter pylori is currently one of those bacterial species which possess the highest number and the most variably expressed set of DNA MTases. Next-generation sequencing technologies can directly detect DNA base methylation. However, they still have limitations in their quantitative and qualitative performance, in particular for cytosine methylation. RESULTS: As a complementing approach, we used enzymatic methyl sequencing (EM-Seq), a technology recently established that has not yet been fully evaluated for bacteria. Thereby, we assessed quantitatively, at single-base resolution, whole genome cytosine methylation for all methylated cytosine motifs in two different H. pylori strains and isogenic MTase mutants. EM-Seq reliably detected both m5C and m4C methylation. We demonstrated that three different active cytosine MTases in H. pylori provide considerably different levels of average genome-wide single-base methylation, in contrast to isogenic mutants which completely lost specific motif methylation. We found that strain identity and changed environmental conditions, such as growth phase and interference with methyl donor homeostasis, significantly influenced quantitative global and local genome-wide methylation in H. pylori at specific motifs. We also identified significantly hyper- or hypo-methylated cytosines, partially linked to overlapping MTase target motifs. Notably, we revealed differentially methylated cytosines in genome-wide coding regions under conditions of methionine depletion, which can be linked to transcript regulation. CONCLUSIONS: This study offers new knowledge on H. pylori global and local genome-wide methylation and establishes EM-Seq for quantitative single-site resolution analyses of bacterial cytosine methylation.

The Helicobacter pylori methylome is acid-responsive due to regulation by the two-component system ArsRS and the type I DNA methyltransferase HsdM1 (HP0463).

In addition to its role in genome protection, DNA methylation can regulate gene expression. In this study, we characterized the impact of acidity, phase variation, and the ArsRS TCS on the expression of the Type I m6A DNA methyltransferase HsdM1 (HP0463) of Helicobacter pylori 26695 and their subsequent effects on the methylome. Transcription of hsdM1 increases at least fourfold in the absence of the sensory histidine kinase ArsS, the major acid-sensing protein of H. pylori. hsdM1 exists in the phase-variable operon hsdR1-hsdM1. Phase-locking hsdR1 (HP0464), the restriction endonuclease gene, has significant impacts on the transcription of hsdM1. To determine the impacts of methyltransferase transcription patterns on the methylome, we conducted methylome sequencing on samples cultured at pH 7 or pH 5. We found differentially methylated motifs between these growth conditions and that deletions of arsS and/or hsdM1 interfere with the epigenetic acid response. Deletion of arsS leads to altered activity of HsdM1 and multiple other methyltransferases under both pH conditions indicating that the ArsRS TCS, in addition to direct effects on regulon transcription during acid acclimation, may also indirectly impact gene expression via regulation of the methylome. We determined the target motif of HsdM1 (HP0463) to be the complementary bipartite sequence pair 5'-TCAm6AVN6TGY-3' and 3'-AGTN6GAm6ACA-5'. This complex regulation of DNA methyltransferases, and thus differential methylation patterns, may have implications for the decades-long persistent infection by H. pylori. IMPORTANCE This study expands the possibilities for complex, epigenomic regulation in Helicobacter pylori. We demonstrate that the H. pylori methylome is plastic and acid sensitive via the two-component system ArsRS and the DNA methyltransferase HsdM1. The control of a methyltransferase by ArsRS may allow for a layered response to changing acidity. Likely, an early response whereby ArsR~P affects regulon expression, including the methyltransferase hsdM1. Then, a somewhat later effect as the altered methylome, due to altered HsdM1 expression, subsequently alters the expression of other genes involved in acclimation. The intermediate methylation of certain motifs supports the hypothesis that methyltransferases play a regulatory role. Untangling this additional web of regulation could play a key role in understanding H. pylori colonization and persistence.

Transcriptomic and methylomic analyses provide insights into the molecular mechanism and prediction of heterosis in rice.

Heterosis has been widely used in multiple crops. However, the molecular mechanism and prediction of heterosis remains elusive. We generated five F1 hybrids [four showing better-parent heterosis (BPH) and one showing mid-parent heterosis], and performed the transcriptomic and methylomic analyses to identify the candidate genes for BPH and explore the molecular mechanism of heterosis and the potential predictors for heterosis. Transcriptomic results showed that most of the differentially expressed genes shared in the four better-parent hybrids were significantly enriched into the terms of molecular function, and the additive and dominant effects played crucial roles for BPH. DNA methylation level, especially in CG context, significantly and positively correlated with grain yield per plant. The ratios of differentially methylated regions in CG context in exons to transcription start sites between the parents exhibited significantly negative correlation with the heterosis levels of their hybrids, as was further confirmed in 24 pairwise comparisons of other rice lines, implying that this ratio could be a feasible predictor for heterosis level, and this ratio of less than 5 between parents in early growth stages might be a critical index for judging that their F1 hybrids would show BPH. Additionally, we identified some important genes showing differential expression and methylation, such as OsDCL2, Pi5, DTH2, DTH8, Hd1 and GLW7 in the four better-parent hybrids as the candidate genes for BPH. Our findings helped shed more light on the molecular mechanism and heterosis prediction.

Sensory plastid-associated PsbP DOMAIN-CONTAINING PROTEIN 3 triggers plant growth- and defense-related epigenetic responses.

Sensory plastids are important in plant responses to environmental changes. Previous studies show that MutS HOMOLOG 1 (MSH1) perturbation in sensory plastids induces heritable epigenetic phenotype adjustment. Previously, the PsbP homolog DOMAIN-CONTAINING PROTEIN 3 (PPD3), a protein of unknown function, was postulated to be an interactor with MSH1. This study investigates the relationship of PPD3 with MSH1 and with plant environmental sensing. The ppd3 mutant displays a whole-plant phenotype variably altered in growth rate, flowering time, reactive oxygen species (ROS) modulation and response to salt, with effects on meristem growth. Present in both chloroplasts and sensory plastids, PPD3 colocalized with MSH1 in root tips but not in leaf tissues. The suppression or overexpression of PPD3 affected the plant growth rate and stress tolerance, and led to a heritable, heterogenous 'memory' state with both dwarfed and vigorous growth phenotypes. Gene expression and DNA methylome data sets from PPD3-OX and derived memory states showed enrichment in growth versus defense networks and meristem effects. Our results support a model of sensory plastid influence on nuclear epigenetic behavior and ppd3 as a second trigger, functioning within meristem plastids to recalibrate growth plasticity.

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.

Cell-type-specific meQTLs extend melanoma GWAS annotation beyond eQTLs and inform melanocyte gene-regulatory mechanisms.

Although expression quantitative trait loci (eQTLs) have been powerful in identifying susceptibility genes from genome-wide association study (GWAS) findings, most trait-associated loci are not explained by eQTLs alone. Alternative QTLs, including DNA methylation QTLs (meQTLs), are emerging, but cell-type-specific meQTLs using cells of disease origin have been lacking. Here, we established an meQTL dataset by using primary melanocytes from 106 individuals and identified 1,497,502 significant cis-meQTLs. Multi-QTL colocalization with meQTLs, eQTLs, and mRNA splice-junction QTLs from the same individuals together with imputed methylome-wide and transcriptome-wide association studies identified candidate susceptibility genes at 63% of melanoma GWAS loci. Among the three molecular QTLs, meQTLs were the single largest contributor. To compare melanocyte meQTLs with those from malignant melanomas, we performed meQTL analysis on skin cutaneous melanomas from The Cancer Genome Atlas (n = 444). A substantial proportion of meQTL probes (45.9%) in primary melanocytes is preserved in melanomas, while a smaller fraction of eQTL genes is preserved (12.7%). Integration of melanocyte multi-QTLs and melanoma meQTLs identified candidate susceptibility genes at 72% of melanoma GWAS loci. Beyond GWAS annotation, meQTL-eQTL colocalization in melanocytes suggested that 841 unique genes potentially share a causal variant with a nearby methylation probe in melanocytes. Finally, melanocyte trans-meQTLs identified a hotspot for rs12203592, a cis-eQTL of a transcription factor, IRF4, with 131 candidate target CpGs. Motif enrichment and IRF4 ChIP-seq analysis demonstrated that these target CpGs are enriched in IRF4 binding sites, suggesting an IRF4-mediated regulatory network. Our study highlights the utility of cell-type-specific meQTLs.

SPEN haploinsufficiency causes a neurodevelopmental disorder overlapping proximal 1p36 deletion syndrome with an episignature of X chromosomes in females.

Deletion 1p36 (del1p36) syndrome is the most common human disorder resulting from a terminal autosomal deletion. This condition is molecularly and clinically heterogeneous. Deletions involving two non-overlapping regions, known as the distal (telomeric) and proximal (centromeric) critical regions, are sufficient to cause the majority of the recurrent clinical features, although with different facial features and dysmorphisms. SPEN encodes a transcriptional repressor commonly deleted in proximal del1p36 syndrome and is located centromeric to the proximal 1p36 critical region. Here, we used clinical data from 34 individuals with truncating variants in SPEN to define a neurodevelopmental disorder presenting with features that overlap considerably with those of proximal del1p36 syndrome. The clinical profile of this disease includes developmental delay/intellectual disability, autism spectrum disorder, anxiety, aggressive behavior, attention deficit disorder, hypotonia, brain and spine anomalies, congenital heart defects, high/narrow palate, facial dysmorphisms, and obesity/increased BMI, especially in females. SPEN also emerges as a relevant gene for del1p36 syndrome by co-expression analyses. Finally, we show that haploinsufficiency of SPEN is associated with a distinctive DNA methylation episignature of the X chromosome in affected females, providing further evidence of a specific contribution of the protein to the epigenetic control of this chromosome, and a paradigm of an X chromosome-specific episignature that classifies syndromic traits. We conclude that SPEN is required for multiple developmental processes and SPEN haploinsufficiency is a major contributor to a disorder associated with deletions centromeric to the previously established 1p36 critical regions.

Defining Interactions Between the Genome, Epigenome, and the Environment in Inflammatory Bowel Disease: Progress and Prospects.

Recent advances in our understanding of the pathogenesis of inflammatory bowel disease (IBD) have highlighted the complex interplay between the genome, the epigenome, and the environment. Despite the exciting advances in genomics that have enabled the identification of over 200 susceptibility loci, these only account for a small proportion of the disease variance and the estimated heritability in IBD. It is likely that gene-environment (GxE) interactions contribute to "missing heritability" and these may act through epigenetic mechanisms. Several environmental factors, such as the microbiome, nutrition, and tobacco smoking, induce alterations in the epigenome of children and adults, which may impact disease susceptibility. Other mechanisms for GxE interactions are also directly pertinent in early life. We discuss a model in which environmental factors imprint disease risk in a window of susceptibility during infancy that may contribute to later disease onset, whereas other elements of the exposome act later in life and contribute directly to the pathogenesis and course of the disease. Understanding the mechanisms underlying GxE interactions may provide the basis for new therapeutic targets or preventative strategies for IBD.

Mdivi-1: Effective but complex mitochondrial fission inhibitor.

Mdivi-1, Mitochondrial DIVIsion inhibitor 1, has been widely employed in research under the assumption that it exclusively influences mitochondrial fusion, but effects other than mitochondrial dynamics have been underinvestigated. This paper provides transcriptome and DNA methylome-wide analysis for Mdivi-1 treated SH-SY5Y human neuroblastoma cells using RNA sequencing (RNA-seq) and methyl capture sequencing (MC-seq) methods. Gene ontology analysis of RNA sequences revealed that p53 transcriptional gene network and DNA replication initiation-related genes were significantly up and down-regulated, respectively, showing the correlation with the arrest cell cycle in the G1 phase. MC-seq, a powerful sequencing method for capturing DNA methylation status in CpG sites, revealed that although Mdivi-1 does not induce dramatic DNA methylation change, the subtle alterations were concentrated within the CpG island. Integrative analysis of both sequencing data disclosed that the p53 transcriptional network was activated while the Parkinson's disease pathway was halted. Next, we investigated several changes in mitochondria in response to Mdivi-1. Copy number and transcription of mitochondrial DNA were suppressed. ROS levels increased, and elevated ROS triggered mitochondrial retrograde signaling rather than inducing direct DNA damage. In this study, we could better understand the molecular network of Mdivi-1 by analyzing DNA methylation and mRNA transcription in the nucleus and further investigating various changes in mitochondria, providing inspiration for studying nuclear-mitochondrial communications.

In search of epigenetic hallmarks of different tissues: an integrative omics study of horse liver, lung, and heart.

DNA methylation and microRNA (miRNA) expression are epigenetic mechanisms essential for regulating tissue-specific gene expression and metabolic processes. However, high-resolution transcriptome, methylome, or miRNAome data is only available for a few model organisms and selected tissues. Up to date, only a few studies have reported on gene expression, DNA methylation, or miRNA expression in adult equine tissues at the genome-wide level. In the present study, we used RNA-Seq, miRNA-seq, and reduced representation bisulfite sequencing (RRBS) data from the heart, lung, and liver tissues of healthy cold-blooded horses to identify differentially expressed genes (DEGs), differentially expressed miRNA (DE miRNA) and differentially methylated sites (DMSs) between three types of horse tissues. Additionally, based on integrative omics analysis, we described the observed interactions of epigenetic mechanisms with tissue-specific gene expression alterations. The obtained data allowed identification from 4067 to 6143 DMSs, 9733 to 11,263 mRNAs, and 155 to 185 microRNAs, differentially expressed between various tissues. We pointed out specific genes whose expression level displayed a negative correlation with the level of CpG methylation and miRNA expression and revealed biological processes that they enrich. Furthermore, we confirmed and validated the accuracy of the Next-Generation Sequencing (NGS) results with bisulfite sequencing PCR (BSP) and quantitative PCR (qPCR). This comprehensive analysis forms a strong foundation for exploring the epigenetic mechanisms involved in tissue differentiation, especially the growth and development of the equine heart, lungs, and liver.

A strategy for studying epigenetic diversity in natural populations: proof of concept in poplar and oak.

In the 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 bisulfite 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 as targets. Given the need to study large numbers of samples, we propose a strategy for investigating DNA methylation variation in natural populations, taking into account the structural complexity of genomes, their size, and their content in unique coding regions versus repeated regions as transposable elements. We first identified regions of highly variable DNA methylation in a 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 bisulfite (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.

A role for DNA methylation in bumblebee morphogenesis hints at female-specific developmental erasure.

Epigenetic mechanisms, such as DNA methylation, are crucial factors in animal development. In some mammals, almost all DNA methylation is erased during embryo development and re-established in a sex- and cell-specific manner. This erasure and re-establishment is thought to primarily be a vertebrate-specific trait. Insects are particularly interesting in terms of development as many species often undergo remarkable morphological changes en route to maturity, that is, morphogenesis. However, little is known about the role of epigenetic mechanisms in this process across species. We have used whole-genome bisulfite sequencing to track genome-wide DNA methylation changes through the development of an economically and environmentally important pollinator species, the bumblebee Bombus terrestris (Hymenoptera:Apidae Linnaeus). We find overall levels of DNA methylation vary throughout development, and we find developmentally relevant differentially methylated genes throughout. Intriguingly, we have identified a depletion of DNA methylation in ovaries/eggs and an enrichment of highly methylated genes in sperm. We suggest this could represent a sex-specific DNA methylation erasure event. To our knowledge, this is the first suggestion of possible developmental DNA methylation erasure in an insect species. This study lays the required groundwork for functional experimental work to determine if there is a causal nature to the DNA methylation differences identified. Additionally, the application of single-cell methylation sequencing to this system will enable more accurate identification of if or when DNA methylation is erased during development.

The utility of DNA methylation profiling in the diagnosis of un-, de- and trans-differentiated melanoma: a series of 11 cases.

AIMS: Melanomas are recognised for their remarkable morphological plasticity. Some tumours may lose conventional features and/or acquire non-melanocytic characteristics, referred to as undifferentiated, dedifferentiated and transdifferentiated melanoma. Despite this phenotypical variability, melanomas typically maintain their cancer driver aberrations, affecting genes such as BRAF, NRAS and NF1. Currently, little is known about whether the DNA methylation profile follows the loss or change of differentiation or is retained despite extensive morphological transformation. METHODS AND RESULTS: In this study we analysed 11 melanoma cases, comprising six males and five females, with a median age of 67 years, including five undifferentiated, four trans-differentiated and two de-differentiated melanomas. Undifferentiated and trans-differentiated tumours either arose in a patient with known melanoma and/or presented in the groin/axilla with molecular alterations consistent with melanoma. Cases with heterologous differentiation resembled chondrosarcoma, osteosarcoma, angiosarcoma and rhabdomyosarcoma both morphologically and immunohistochemically, while undifferentiated tumours resembled undifferentiated pleomorphic sarcoma. Methylome profiling was performed, and unsupervised clustering analysis revealed nine cases (five undifferentiated, three trans-differentiated and one de-differentiated) to cluster closely together with conventional melanomas from a reference set. Two cases clustered separately with a distinct group of conventional melanomas exhibiting H3K27me3 loss. CONCLUSIONS: Despite loss of differentiation and phenotypical plasticity, methylation patterns seem to be retained in undifferentiated, de-differentiated and trans-differentiated melanomas and represent useful diagnostic tools to enhance diagnostic precision in these diagnostically challenging cases.