Antagonistic and synergistic epigenetic modulation using orthologous CRISPR/dCas9-based modular system.

Establishing causal relationship between epigenetic marks and gene transcription requires molecular tools, which can precisely modify specific genomic regions. Here, we present a modular and extensible CRISPR/dCas9-based toolbox for epigenetic editing and direct gene regulation. It features a system for expression of orthogonal dCas9 proteins fused to various effector domains and includes a multi-gRNA system for simultaneous targeting dCas9 orthologs to up to six loci. The C- and N-terminal dCas9 fusions with DNMT3A and TET1 catalytic domains were thoroughly characterized. We demonstrated simultaneous use of the DNMT3A-dSpCas9 and TET1-dSaCas9 fusions within the same cells and showed that imposed cytosine hyper- and hypo-methylation altered level of gene transcription if targeted CpG sites were functionally relevant. Dual epigenetic manipulation of the HNF1A and MGAT3 genes, involved in protein N-glycosylation, resulted in change of the glycan phenotype in BG1 cells. Furthermore, simultaneous targeting of the TET1-dSaCas9 and VPR-dSpCas9 fusions to the HNF1A regulatory region revealed strong and persistent synergistic effect on gene transcription, up to 30 days following cell transfection, suggesting involvement of epigenetic mechanisms in maintenance of the reactivated state. Also, modulation of dCas9 expression effectively reduced off-target effects while maintaining the desired effects on target regions.

CRISPR/Cas9-based epigenome editing: An overview of dCas9-based tools with special emphasis on off-target activity.

Molecular tools for gene regulation and epigenome editing consist of two main parts: the targeting moiety binding a specific genomic locus and the effector domain performing the editing or regulatory function. The advent of CRISPR-Cas9 technology enabled easy and flexible targeting of almost any locus by co-expression of a small sgRNA molecule, which is complementary to the target sequence and forms a complex with Cas9, directing it to that particular target. Here, we review strategies for recruitment of effector domains, used in gene regulation and epigenome editing, to the dCas9 DNA-targeting protein. To date, the most important CRISPR-Cas9 applications in gene regulation are CRISPR activation or interference, while epigenome editing focuses on targeted changes in DNA methylation and histone modifications. Several strategies for signal amplification by recruitment of multiple effector domains deserve special focus. While some approaches rely on altering the sgRNA molecule and extending it with aptamers for effector domain recruitment, others use modifications to the Cas9 protein by direct fusions with effector domains or by addition of an epitope tag, which also has the ability to bind multiple effector domains. A major barrier to the widespread use of CRISPR-Cas9 technology for therapeutic purposes is its off-target effect. We review efforts to enhance CRISPR-Cas9 specificity by selection of Cas9 orthologs from various bacterial species and their further refinement by introduction of beneficial mutations. The molecular tools available today enable a researcher to choose the best balance of targeting flexibility, activity amplification, delivery method and specificity.

Unlocking Cancer Glycomes from Histopathological Formalin-fixed and Paraffin-embedded (FFPE) Tissue Microdissections.

N- and O-glycans are attractive clinical biomarkers as glycosylation changes in response to diseases. The limited availability of defined clinical specimens impedes glyco-biomarker identification and validation in large patient cohorts. Formalin-fixed paraffin-embedded (FFPE) clinical specimens are the common form of sample preservation in clinical pathology, but qualitative and quantitative N- and O-glycomics of such samples has not been feasible to date. Here, we report a highly sensitive and glycan isomer selective method for simultaneous N- and O-glycomics from histopathological slides. As few as 2000 cells isolated from FFPE tissue sections by laser capture microdissection were sufficient for in-depth histopathology-glycomics using porous graphitized carbon nanoLC ESI-MS/MS. N- and O-glycan profiles were similar between unstained and hematoxylin and eosin stained FFPE samples but differed slightly compared with fresh tissue. This method provides the key to unlock glyco-biomarker information from FFPE histopathological tissues archived in pathology laboratories worldwide.

Repurposing the CRISPR-Cas9 system for targeted DNA methylation.

Epigenetic studies relied so far on correlations between epigenetic marks and gene expression pattern. Technologies developed for epigenome editing now enable direct study of functional relevance of precise epigenetic modifications and gene regulation. The reversible nature of epigenetic modifications, including DNA methylation, has been already exploited in cancer therapy for remodeling the aberrant epigenetic landscape. However, this was achieved non-selectively using epigenetic inhibitors. Epigenetic editing at specific loci represents a novel approach that might selectively and heritably alter gene expression. Here, we developed a CRISPR-Cas9-based tool for specific DNA methylation consisting of deactivated Cas9 (dCas9) nuclease and catalytic domain of the DNA methyltransferase DNMT3A targeted by co-expression of a guide RNA to any 20 bp DNA sequence followed by the NGG trinucleotide. We demonstrated targeted CpG methylation in a ∼35 bp wide region by the fusion protein. We also showed that multiple guide RNAs could target the dCas9-DNMT3A construct to multiple adjacent sites, which enabled methylation of a larger part of the promoter. DNA methylation activity was specific for the targeted region and heritable across mitotic divisions. Finally, we demonstrated that directed DNA methylation of a wider promoter region of the target loci IL6ST and BACH2 decreased their expression.

Glyco-genes change expression in cancer through aberrant methylation.

BACKGROUND: Most eukaryotic proteins are modified by covalent addition of glycan molecules that considerably influence their function. Aberrant glycosylation is profoundly involved in malignant transformation, tumor progression and metastasis. Some glycan structures are tumor-specific and reflect disturbed glycan biosynthesis pathways. METHODS: We analyzed DNA methylation and expression of 86 glyco-genes in melanoma, hepatocellular, breast and cervical cancers using data from publicly available databases. We also analyzed methylation datasets without the available matching expression data for glyco-genes in lung cancer, and progression of melanoma into lymph node and brain metastases. RESULTS: Ten glyco-genes (GALNT3, GALNT6, GALNT7, GALNT14, MGAT3, MAN1A1, MAN1C1, ST3GAL2, ST6GAL1, ST8SIA3) showing changes in both methylation and expression in the same type of cancer belong to GalNAc transferases, GlcNAc transferases, mannosidases and sialyltransferases, which is in line with changes in glycan structures already reported in the same type of tumors. Some of those genes were additionally identified as potentially valuable for disease prognosis. The MGAT5B gene, so far identified as specifically expressed in brain, emerged as a novel candidate gene that is epigenetically dysregulated in different cancers other than brain cancer. We also report for the first time aberrant expression of the GALNT and MAN genes in cancer by aberrant promoter methylation. CONCLUSIONS: Aberrant expression of glyco-genes due to aberrant promoter methylation could be a way leading to characteristic glycosylation profiles commonly described in cancer. GENERAL SIGNIFICANCE: Methylation status in promoters of candidate glyco-genes might serve as prognostic markers for specific tumors and point to potential novel targets for epigenetic drugs. This article is part of a Special Issue entitled "Glycans in personalised medicine" Guest Editor: Professor Gordan Lauc.

Epigenetic regulation of glycosylation is the quantum mechanics of biology.

BACKGROUND: Most proteins are glycosylated, with glycans being integral structural and functional components of a glycoprotein. In contrast to polypeptides, which are fully encoded by the corresponding gene, glycans result from a dynamic interaction between the environment and a network of hundreds of genes. SCOPE OF REVIEW: Recent developments in glycomics, genomics and epigenomics are discussed in the context of an evolutionary advantage for higher eukaryotes over microorganisms, conferred by the complexity and adaptability which glycosylation adds to their proteome. MAJOR CONCLUSIONS: Inter-individual variation of glycome composition in human population is large; glycome composition is affected by both genes and environment; epigenetic regulation of "glyco-genes" has been demonstrated; and several mechanisms for transgenerational inheritance of epigenetic marks have been documented. GENERAL SIGNIFICANCE: Epigenetic recording of acquired characteristics and their transgenerational inheritance could be important mechanisms used by higher organisms to compete or collaborate with microorganisms.

Epigenetic modulation of the HeLa cell membrane N-glycome.

BACKGROUND: Epigenetic changes play a role in all major events during tumorigenesis and changes in glycan structures are hallmarks of virtually every cancer. Also, proper N-glycosylation of membrane receptors is important in cell to cell and cell-environment communication. To study how modulation of epigenetic information can affect N-glycan expression we analyzed effects of epigenetic inhibitors on HeLa cell membrane N-glycome. METHODS: HeLa cells were treated with DNA methylation (zebularin and 5-aza-2-deoxycytidine) and histone deacetylation (trichostatin A and Na-butyrate) inhibitors. The effects on HeLa cell membrane N-glycome were analyzed by hydrophilic interaction high performance liquid chromatography (HILIC). RESULTS: Each of the four epigenetic inhibitors induced changes in the expression of HeLa cell membrane N-glycans that were seen either as an increase or a decrease of individual glycans in the total N-glycome. Compared to DNA methylation inhibitors, histone deacetylation inhibitors showed more moderate changes, probably due to their higher gene target selectivity. CONCLUSIONS: The results clearly show that composition of HeLa cell membrane N-glycome can be specifically altered by epigenetic inhibitors. GENERAL SIGNIFICANCE: Glycans on the cell membrane are essential elements of tumor cell's metastatic potential and are also an entry point for nearly all pathogenic microorganisms. Since epigenetic inhibitors used in this work are registered drugs, our results provide a new line of research in the application of these drugs as anticancer and antimicrobial agents. This article is part of a Special Issue entitled Glycoproteomics.

Genomics and epigenomics of the human glycome.

The majority of all proteins are glycosylated and glycans have numerous important structural, functional and regulatory roles in various physiological processes. While structure of the polypeptide part of a glycoprotein is defined by the sequence of nucleotides in the corresponding gene, structure of a glycan part results from dynamic interactions between hundreds of genes, their protein products and environmental factors. The composition of the glycome attached to an individual protein, or to a complex mixture of proteins, like human plasma, is stable within an individual, but very variable between individuals. This variability stems from numerous common genetic polymorphisms reflecting in changes in the complex biosynthetic pathway of glycans, but also from the interaction with the environment. Environment can affect glycan biosynthesis at the level of substrate availability, regulation of enzyme activity and/or hormonal signals, but also through gene-environment interactions. Epigenetics provides a molecular basis how the environment can modify phenotype of an individual. The epigenetic information (DNA methylation pattern and histone code) is especially vulnerable to environmental effects in the early intrauterine and neo-natal development and many common late-onset diseases take root already at that time. The evidences showing the link between epigenetics and glycosylation are accumulating. Recent progress in high-throughput glycomics, genomics and epigenomics enabled first epidemiological and genome-wide association studies of the glycome, which are presented in this mini-review.

Long-Term Culturing of FreeStyle 293-F Cells Affects Immunoglobulin G Glycome Composition.

Glycosylation of IgG regulates the effector function of this antibody in the immune response. Glycosylated IgG is a potent therapeutic used for both research and clinical purposes. While there is ample research on how different cell culture conditions affect IgG glycosylation, the data are missing on the stability of IgG glycome during long cell passaging, i.e., cell "aging". To test this, we performed three independent time course experiments in FreeStyle 293-F cells, which secrete IgG with a human-like glycosylation pattern and are frequently used to generate defined IgG glycoforms. During long-term cell culturing, IgG glycome stayed fairly stable except for galactosylation, which appeared extremely variable. Cell transcriptome analysis revealed no correlation in galactosyltransferase B4GALT1 expression with galactosylation change, but with expression of EEF1A1 and SLC38A10, genes previously associated with IgG galactosylation through GWAS. The FreeStyle 293-F cell-based system for IgG production is a good model for studies of mechanisms underlying IgG glycosylation, but results from the present study point to the utmost importance of the need to control IgG galactosylation in both in vitro and in vivo systems. This is especially important for improving the production of precisely glycosylated IgG for therapeutic purposes, since IgG galactosylation affects the inflammatory potential of IgG.

Mapping of the gene network that regulates glycan clock of ageing.

Glycans are an essential structural component of immunoglobulin G (IgG) that modulate its structure and function. However, regulatory mechanisms behind this complex posttranslational modification are not well known. Previous genome-wide association studies (GWAS) identified 29 genomic regions involved in regulation of IgG glycosylation, but only a few were functionally validated. One of the key functional features of IgG glycosylation is the addition of galactose (galactosylation), a trait which was shown to be associated with ageing. We performed GWAS of IgG galactosylation (N=13,705) and identified 16 significantly associated loci, indicating that IgG galactosylation is regulated by a complex network of genes that extends beyond the galactosyltransferase enzyme that adds galactose to IgG glycans. Gene prioritization identified 37 candidate genes. Using a recently developed CRISPR/dCas9 system we manipulated gene expression of candidate genes in the in vitro IgG expression system. Upregulation of three genes, EEF1A1, MANBA and TNFRSF13B, changed the IgG glycome composition, which confirmed that these three genes are involved in IgG galactosylation in this in vitro expression system.

Lentivirus- or AAV-mediated gene therapy interventions in ischemic stroke: A systematic review of preclinical in vivo studies.

Due to the limited therapeutic options after ischemic stroke, gene therapy has emerged as a promising choice, especially with recent advances in viral vector delivery systems. Therefore, we aimed to provide the current state of the art of lentivirus (LV) and adeno-associated virus (AAV) mediated gene interventions in preclinical ischemic stroke models. A systematic analysis including qualitative and quantitative syntheses of studies published until December 2020 was performed. Most of the 87 selected publications used adult male rodents and the preferred stroke model was transient middle cerebral artery occlusion. LV and AAV vectors were equally used for transgene delivery, however loads of AAVs were higher than LVs. Serotypes having broad cell tropism, the use of constitutive promoters, and virus delivery before the stroke induction via stereotaxic injection in the cortex and striatum were preferred in the analyzed studies. The meta-analysis based on infarct volume as the primary outcome confirmed the efficacy of the preclinical interventions. The quality assessment exposed publication bias and setbacks in regard to risks of bias and study relevance. The translational potential could increase by using specific cell targeting, post-stroke interventions, non-invasive systematic delivery, and use of large animals.

Genetic and Epigenetic Association of Hepatocyte Nuclear Factor-1α with Glycosylation in Post-Traumatic Stress Disorder.

Post-traumatic stress disorder (PTSD) is a complex trauma-related disorder, the etiology and underlying molecular mechanisms of which are still unclear and probably involve different (epi)genetic and environmental factors. Protein N-glycosylation is a common post-translational modification that has been associated with several pathophysiological states, including inflammation and PTSD. Hepatocyte nuclear factor-1α (HNF1A) is a transcriptional regulator of many genes involved in the inflammatory processes, and it has been identified as master regulator of plasma protein glycosylation. The aim of this study was to determine the association between N-glycan levels in plasma and immunoglobulin G, methylation at four CpG positions in the HNF1A gene, HNF1A antisense RNA 1 (HNF1A-AS1), rs7953249 and HNF1A rs735396 polymorphisms in a total of 555 PTSD and control subjects. We found significant association of rs7953249 and rs735396 polymorphisms, as well as HNF1A gene methylation at the CpG3 site, with highly branched, galactosylated and sialyated plasma N-glycans, mostly in patients with PTSD. HNF1A-AS1 rs7953249 polymorphism was also associated with PTSD; however, none of the polymorphisms were associated with HNF1A gene methylation. These results indicate a possible regulatory role of the investigated HNF1A polymorphisms with respect to the abundance of complex plasma N-glycans previously associated with proinflammatory response, which could contribute to the clinical manifestation of PTSD and its comorbidities.

A Transient Expression System with Stably Integrated CRISPR-dCas9 Fusions for Regulation of Genes Involved in Immunoglobulin G Glycosylation.

Alternative glycosylation of immunoglobulin G (IgG) is functionally important in multiple human physiological and pathological states. Our understanding of molecular mechanisms that regulate IgG glycosylation is vague because of the complexity of this process, which involves hundreds of genes. Several genome-wide association (GWA) studies have revealed a network of genes associated with IgG glycosylation that are pleiotropic for a number of diseases. Here, we report a design of a versatile system for IgG production and gene manipulations that can be used for in vitro functional follow-up of GWA hits or any gene of interest. The system is based on CRISPR-dCas9, extended by a piggyBac integrase compatible vector, and drives IgG production in HEK-293F cells. We validated our systems that stably express VPR-dCas9 and KRAB-dCas9 by manipulation of four glyco-genes with a known role in IgG glycosylation, and then functionally validated three GWAS hits for IgG glycosylation with an as-yet-unknown role in this process.

Heritability of the glycan clock of biological age.

Immunoglobulin G is posttranslationally modified by the addition of complex N-glycans affecting its function and mediating inflammation at multiple levels. IgG glycome composition changes with age and health in a predictive pattern, presumably due to inflammaging. As a result, a novel biological aging biomarker, glycan clock of age, was developed. Glycan clock of age is the first of biological aging clocks for which multiple studies showed a possibility of clock reversal even with simple lifestyle interventions. However, none of the previous studies determined to which extent the glycan clock can be turned, and how much is fixed by genetic predisposition. To determine the contribution of genetic and environmental factors to phenotypic variation of the glycan clock, we performed heritability analysis on two TwinsUK female cohorts. IgG glycans from monozygotic and dizygotic twin pairs were analyzed by UHPLC and glycan age was calculated using the glycan clock. In order to determine additive genetic, shared, and unique environmental contributions, a classical twin design was applied. Heritability of the glycan clock was calculated for participants of one cross-sectional and one longitudinal cohort with three time points to assess the reliability of measurements. Heritability estimate for the glycan clock was 39% on average, suggesting a moderate contribution of additive genetic factors (A) to glycan clock variation. Remarkably, heritability estimates remained approximately the same in all time points of the longitudinal study, even though IgG glycome composition changed substantially. Most environmental contributions came from shared environmental factors (C), with unique environmental factors (E) having a minor role. Interestingly, heritability estimates nearly doubled, to an average of 71%, when we included age as a covariant. This intervention also inflated the estimates of unique environmental factors contributing to glycan clock variation. A complex interplay between genetic and environmental factors defines alternative IgG glycosylation during aging and, consequently, dictates the glycan clock's ticking. Apparently, environmental factors (including lifestyle choices) have a strong impact on the biological age measured with the glycan clock, which additionally clarifies why this aging clock is one of the most potent biomarkers of biological aging.

Epigenetics of Immunoglobulin G Glycosylation.

Alternative glycosylation of immunoglobulin G (IgG) affects its effector functions during the immune response. IgG glycosylation is altered in many diseases, but also during a healthy life of an individual. Currently, there is limited knowledge of factors that alter IgG glycosylation in the healthy state and factors involved in specific IgG glycosylation patterns associated with pathophysiology. Genetic background plays an important role, but epigenetic mechanisms also contribute to the alteration of IgG glycosylation patterns in healthy life and in disease. It is known that the expression of many glycosyltransferases is regulated by DNA methylation and by microRNA (miRNA) molecules, but the involvement of other epigenetic mechanisms, such as histone modifications, in the regulation of glycosylation-related genes (glycogenes) is still poorly understood. Recent studies have identified several differentially methylated loci associated with IgG glycosylation, but the mechanisms involved in the formation of specific IgG glycosylation patterns remain poorly understood.

N-glycans as functional effectors of genetic and epigenetic disease risk.

N-glycosylation is a frequent modification of proteins, essential for all domains of life. N-glycan biosynthesis is a dynamic, complex, non-templated process, wherein specific glycoforms are modulated by various microenvironmental cues, cellular signals and local availability of dedicated enzymes and sugar precursors. This intricate regulatory network comprises hundreds of proteins, whose activity is dependent on both sequence of implicated genes and the regulation of their expression. In this regard, variation in N-glycosylation patterns stems from either gene polymorphisms or from stable epigenetic regulation of gene expression in different individuals. Moreover, epigenome alters in response to various environmental factors, representing a direct link between environmental exposure and changes in gene expression, that are subsequently reflected through altered N-glycosylation. N-glycosylation itself has a fundamental role in numerous biological processes, ranging from protein folding, cellular homeostasis, adhesion and immune regulation, to the effector functions in multiple diseases. Moreover, specific modification of the glycan structure can modulate glycoprotein's biological function or direct the faith of the entire cell, as seen on the examples of antibodies and T cells, respectively. Since immunoglobulin G is one of the most profoundly studied glycoproteins in general, the focus of this review will be on its N-glycosylation changes and their functional implications. By deepening the knowledge on the mechanistic roles that certain glycoforms exert in differential pathological processes, valuable insight into molecular perturbations occurring during disease development could be obtained. The prospect of resolving the exact biological pathways involved offers a potential for the development of new therapeutic interventions and molecular tools that would aid in prognosis, early referral and timely treatment of multiple disease conditions.

Effects of Estradiol on Immunoglobulin G Glycosylation: Mapping of the Downstream Signaling Mechanism.

Glycans attached to immunoglobulin G (IgG) directly affect this antibody effector functions and regulate inflammation at several levels. The composition of IgG glycome changes significantly with age. In women, the most notable change coincides with the perimenopausal period. Aiming to investigate the effect of estrogen on IgG glycosylation, we analysed IgG and total serum glycomes in 36 healthy premenopausal women enrolled in a randomized controlled trial of the gonadotropin-releasing hormone analogue (GnRHAG) leuprolide acetate to lower gonadal steroids to postmenopausal levels and then randomized to transdermal placebo or estradiol (E2) patch. The suppression of gonadal hormones induced significant changes in the IgG glycome, while E2 supplementation was sufficient to prevent changes. The observed glycan changes suggest that depletion of E2 primarily affects B cell glycosylation, while liver glycosylation stays mostly unchanged. To determine whether previously identified IgG GWAS hits RUNX1, RUNX3, SPINK4, and ELL2 are involved in downstream signaling mechanisms, linking E2 with IgG glycosylation, we used the FreeStyle 293-F transient system expressing IgG antibodies with stably integrated CRISPR/dCas9 expression cassettes for gene up- and downregulation. RUNX3 and SPINK4 upregulation using dCas9-VPR resulted in a decreased IgG galactosylation and, in the case of RUNX3, a concomitant increase in IgG agalactosylation.

Glycosylation of immunoglobulin G is regulated by a large network of genes pleiotropic with inflammatory diseases.

Effector functions of immunoglobulin G (IgG) are regulated by the composition of a glycan moiety, thus affecting activity of the immune system. Aberrant glycosylation of IgG has been observed in many diseases, but little is understood about the underlying mechanisms. We performed a genome-wide association study of IgG N-glycosylation (N = 8090) and, using a data-driven network approach, suggested how associated loci form a functional network. We confirmed in vitro that knockdown of IKZF1 decreases the expression of fucosyltransferase FUT8, resulting in increased levels of fucosylated glycans, and suggest that RUNX1 and RUNX3, together with SMARCB1, regulate expression of glycosyltransferase MGAT3. We also show that variants affecting the expression of genes involved in the regulation of glycoenzymes colocalize with variants affecting risk for inflammatory diseases. This study provides new evidence that variation in key transcription factors coupled with regulatory variation in glycogenes modifies IgG glycosylation and has influence on inflammatory diseases.

Active fusions of Cas9 orthologs.

Many recent epigenetic studies utilize the advantages of CRISPR/dCas9 based tools in linking certain epigenetic modification with gene expression regulation. Various multifactorial diseases often contain changed epigenetic signatures at many loci, so tools for simultaneously targeting different loci would significantly facilitate the understanding of disease pathogenesis. We tested different dCas9 orthologs (dCjCas9, dNmCas9, dSt1Cas9, dFnCas9, dSaCas9 and dSpCas9) in C-terminal fusion with DNMT3A effector domain to find candidates that potentiate effector domain to perform its function at the target site. We demonstrated that nuclear localization signals (NLS) at both termini of fusion constructs is crucial for both proper nuclear import of such large constructs as well as for maximization of targeted DNA methylation activity. We identified SpCas9, SaCas9 and CjCas9 as potential candidates for the fusion constructs. With further optimization of the SaCas9 ortholog, due to less complex PAM requirements in contrast to CjCas9, we showed that N-terminal fusion with DNMT3A (dSaCas9-DNMT3A) is optimal to exert targeted DNA methylation activity comparable to the dSpCas9-DNMT3A construct. N-terminal fusions showed better results for both Cas9 orthologs, SaCas9 and SpCas9, so it can be used as universal approach for linking different effector domains in order to obtain highly active fusions.

Epigenetic Differentiation of Natural Populations of Lilium bosniacum Associated with Contrasting Habitat Conditions.

Epigenetic variation in natural populations with contrasting habitats might be an important element, in addition to the genetic variation, in plant adaptation to environmental stress. Here, we assessed genetic, epigenetic, and cytogenetic structure of the three Lilium bosniacum populations growing on distinct habitats. One population was growing under habitual ecological conditions for this species and the other two were growing under stress associated with high altitude and serpentine soil. Amplified fragment length polymorphism and methylation-sensitive amplification polymorphism analyses revealed that the three populations did not differentiate genetically, but were clearly separated in three distinct clusters according to DNA methylation profiles. Principal coordinate analysis showed that overall epigenetic variation was closely related to habitat conditions. A new methylation-sensitive amplification polymorphism scoring approach allowed identification of mainly unmethylated (φST = 0.190) and fully CpG methylated (φST = 0.118) subepiloci playing a role in overall population differentiation, in comparison with hemimethylated sites (φST = 0.073). In addition, unusual rDNA repatterning and the presence of B chromosomes bearing 5S rDNA loci were recorded in the population growing on serpentine soil, suggesting dynamic chromosome rearrangements probably linked to global genome demethylation, which might have reactivated some mobile elements. We discuss our results considering our earlier data on morphology and leaf anatomy of several L. bosniacum populations, and suggest a possible role of epigenetics as a key element in population differentiation associated with environmental stress in these particular lily populations.