CG Dinucleotides Research Articles

C G composition in transposon-derived genes is increased in FXD with perturbed immune system.

Increasing incidence of Fragile X disorders (FXD) and of immune-mediated disorders in FXD suggests that additional factors besides FMR1 mutations contribute to the pathogenesis. Here, we discovered that the expression levels or splicing of specific transposon element (TE)-derived genes, regulating purine metabolism and immune responses against viral infections are altered in FXD. These genes include HLA genes clustered in chr6p21.3 and viral responsive genes in chr5q15. Remarkably, these TE-derived genes contain a low A T/C G suggesting base substitutions of A T to C G. The TE-derived genes with changed expression levels contained a higher content of 5'-CG-3' dinucleotides in FXD compared to healthy donors. This resembles the genomes of some RNA viruses, which maintain high contents of CG dinucleotides to sustain their latent infection exploiting antiviral responses. Thus, past viral infections may have persisted as TEs, provoking immune-mediated disorders in FXD.

DNA Methylation Participates in Drought Stress Memory and Response to Drought in Medicago ruthenica.

Background: Drought is currently a global environmental problem, which inhibits plant growth and development and seriously restricts crop yields. Many plants exposed to drought stress can generate stress memory, which provides some advantages for resisting recurrent drought. DNA methylation is a mechanism involved in stress memory formation, and many plants can alter methylation levels to form stress memories; however, it remains unclear whether Medicago ruthenica exhibits drought stress memory, as the epigenetic molecular mechanisms underlying this process have not been described in this species. Methods: We conducted methylome and transcriptome sequencing to identify gene methylation and expression changes in plants with a history of two drought stress exposures. Results: Methylation analysis showed that drought stress resulted in an approximately 4.41% decrease in M. ruthenica genome methylation levels. The highest methylation levels were in CG dinucleotide contexts, followed by CHG contexts, with CHH contexts having the lowest levels. Analysis of associations between methylation and transcript levels showed that most DNA methylation was negatively correlated with gene expression except methylation within CHH motifs in gene promoter regions. Genes were divided into four categories according to the relationship between methylation and gene expression; the up-regulation of hypo-methylated gene expression accounted for the vast majority (692 genes) and included genes encoding factors key for abscisic acid (ABA) and proline synthesis. The hypo-methylation of the promoter and body regions of these two gene groups induced increased gene transcription levels. Conclusions: In conclusion, DNA methylation may contribute to drought stress memory formation and maintenance in M. ruthenica by increasing the transcription levels of genes key for ABA and proline biosynthesis.

Distribution rules of 8-mer spectra and characterization of evolution state in animal genome sequences

BackgroundStudying the composition rules and evolution mechanisms of genome sequences are core issues in the post-genomic era, and k-mer spectrum analysis of genome sequences is an effective means to solve this problem.ResultWe divided total 8-mers of genome sequences into 16 kinds of XY-type due to XY dinucleotides number in 8-mers. Previous works explored that the independent unimodal distributions observed only in three CG-type 8-mer spectra, while non-CG type 8-mer spectra have not the universal phenomenon from prokaryotes to eukaryotes. On this basis, we analyzed the distribution variation of non-CG type 8-mer spectra across 889 animal genome sequences. Following the evolutionary order of animals from primitive to more complex, we found that the spectrum distributions gradually transition from unimodal to tri-modal. The relative distance from the average frequency of each non-CG type 8-mers to the center frequency is different within a species and among different species. For the 8-mers contain CG dinucleotides, we further divided these into 16 subsets, where each 8-mer contains both CG and XY dinucleotides, called XY1_CG1 subsets. We found that the separability values of XY1_CG1 spectra are closely related to the evolution and specificity of animals. Considering the constraint of Chargaff’s second parity rule, we finally obtained 10 separability values as the feature set to characterize the evolution state of genome sequences. In order to verify the rationality of the feature set, we used 14 common classification algorithms to perform binary classification tests. The results showed that the accuracy (Acc) ranged between 98.70% and 83.88% among birds, other vertebrates and mammals.ConclusionWe proposed a credible feature set to characterizes the evolution state of genomes and obtained satisfied results by the feature set on large scale classification of animals.

Codon Usage Analysis Reveals Distinct Evolutionary Patterns and Host Adaptation Strategies in Duck Hepatitis Virus 1 (DHV-1) Phylogroups.

Duck hepatitis virus 1 (DHV-1) is a major threat to the global poultry industry, causing significant economic losses due to high mortality rates in young ducklings. To better understand the evolution and host adaptation strategies of DHV-1, we conducted a comprehensive codon usage analysis of DHV-1 genomes. Our phylogenetic analysis revealed three well-supported DHV-1 phylogroups (Ia, Ib, and II) with distinct genetic diversity patterns. Comparative analyses of the codon usage bias and dinucleotide abundance uncovered a strong preference for A/U-ended codons and a biased pattern of dinucleotide usage in the DHV-1 genome, with CG dinucleotides being extremely underrepresented. Effective number of codons (ENC) analysis indicated a low codon usage bias in the DHV-1 ORF sequences, suggesting adaptation to host codon usage preferences. PR2 bias, ENC plot, and neutrality analyses revealed that both mutation pressure and natural selection influence the codon usage patterns of DHV-1. Notably, the three DHV-1 phylogroups exhibited distinct evolutionary trends, with phylogroups Ia and Ib showing evidence of neutral evolution accompanied by selective pressure, while the phylogroup II evolution was primarily driven by random genetic drift. Comparative analysis of the codon usage indices (CAI, RCDI, and SiD) among the phylogroups highlighted significant differences between subgroups Ia and Ib, suggesting distinct evolutionary pressures or adaptations influencing their codon usage. These findings contribute to our understanding of DHV-1 evolution and host adaptation, with potential implications for the development of effective control measures and vaccines.

Investigation of CaMV-host co-evolution through synonymous codon pattern.

Cauliflower mosaic virus (CaMV) has a double-stranded DNA genome and is globally distributed. The phylogeny tree of 121 CaMV isolates was categorized into two primary groups, with Iranian isolates showing the greatest genetic variations. Nucleotide A demonstrated the highest percentage (36.95%) in the CaMV genome and the dinucleotide odds ratio analysis revealed that TC dinucleotide (1.34 ≥ 1.23) and CG dinucleotide (0.63 ≤ 0.78) are overrepresented and underrepresented, respectively. Relative synonymous codon usage (RSCU) analysis confirmed codon usage bias in CaMV and its hosts. Brassica oleracea and Brassica rapa, among the susceptible hosts of CaMV, showed a codon adaptation index (CAI) value above 0.8. Additionally, relative codon deoptimization index (RCDI) results exhibited the highest degree of deoptimization in Raphanus sativus. These findings suggest that the genes of CaMV underwent codon adaptation with its hosts. Among the CaMV open reading frames (ORFs), genes that produce reverse transcriptase and virus coat proteins showed the highest CAI value of 0.83. These genes are crucial for the creation of new virion particles. The results confirm that CaMV co-evolved with its host to ensure the optimal expression of its genes in the hosts, allowing for easy infection and effective spread. To detect the force behind codon usage bias, an effective number of codons (ENC)-plot and neutrality plot were conducted. The results indicated that natural selection is the primary factor influencing CaMV codon usage bias.

CpG Islands, Gene Expression and Pseudogenization: A Case for a Potential Trilogy.

The promoters of mammalian genes contain clusters of CG dinucleotides known as CpG islands. Most mammalian housekeeping genes predominantly contain CpG islands (CGIs), facilitating gene transcription. Numerous studies have explored the physiological implications of the relationship between CGIs and gene expression. However, the evolutionary implications of this relationship remain largely unexplored. Pseudogenes, in contrast, are genomic remnants that have lost their function over evolutionary time. In our current research, we employed comparative genomic techniques to demonstrate a correlation between the absence of gene expression due to a lack of CGIs in the gene promoters and pseudogenization. We showed that there is a significant enrichment of tissue-specific genes in the functional orthologs of pseudogenes. We also found a significant correlation between the lack of CGIs and enriched tissue specificity in these functional orthologs of pseudogenes. We inferred that perhaps tissue-specific genes are more prone to the process of pseudogenization. In this way, because of their impact on gene expression, CGIs may affect the fate of a gene. To our knowledge, this is the first study to propose a connection between CGIs, gene expression, and the pseudogenization process and discuss the evolutionary implications of this potential trilogy.

DNA conformational dynamics in the context-dependent non-CG CHH methylation by plant methyltransferase DRM2

DNA methylation provides an important epigenetic mechanism that critically regulates gene expression, genome imprinting, and retrotransposon silencing. In plants, DNA methylation is prevalent not only in a CG dinucleotide context, but also in non-CG contexts, namely CHG and CHH (H = C, T, or A) methylation. It has been established that plant non-CG DNA methylation is highly context dependent, with the +1- and +2-flanking sequences enriched with A/T nucleotides. How DNA sequence, conformation, and dynamics influence non-CG methylation remains elusive. Here, we report structural and biochemical characterizations of the intrinsic substrate preference of DOMAINS REARRANGED METHYLTRANSFERASE 2 (DRM2), a plant DNA methyltransferase responsible for establishing all cytosine methylation and maintaining CHH methylation. Among nine CHH motifs, the DRM2 methyltransferase (MTase) domain shows marked substrate preference toward CWW (W = A or T) motifs, correlating well with their relative abundance in planta. Furthermore, we report the crystal structure of DRM2 MTase in complex with a DNA duplex containing a flexible TpA base step at the +1/+2-flanking sites of the target nucleotide. Comparative structural analysis of the DRM2-DNA complexes provides a mechanism by which flanking nucleotide composition impacts DRM2-mediated DNA methylation. Furthermore, the flexibility of the TpA step gives rise to two alternative DNA conformations, resulting in different interactions with DRM2 and consequently temperature-dependent shift of the substrate preference of DRM2. Together, this study provides insights into how the interplay between the conformational dynamics of DNA and temperature as an environmental factor contributes to the context-dependent CHH methylation by DRM2.

Millennia-long epigenetic fluctuations generate intragenic DNA methylation variance in Arabidopsis populations

Methylation of CG dinucleotides (mCGs), which regulates eukaryotic genome functions, is epigenetically propagated by Dnmt1/MET1 methyltransferases. How mCG is established and transmitted across generations despite imperfect enzyme fidelity is unclear. Whether mCG variation in natural populations is governed by genetic or epigenetic inheritance also remains mysterious. Here, we show that MET1 de novo activity, which is enhanced by existing proximate methylation, seeds and stabilizes mCG in Arabidopsis thaliana genes. MET1 activity is restricted by active demethylation and suppressed by histone variant H2A.Z, producing localized mCG patterns. Based on these observations, we develop a stochastic mathematical model that precisely recapitulates mCG inheritance dynamics and predicts intragenic mCG patterns and their population-scale variation given only CG site spacing. Our results demonstrate that intragenic mCG establishment, inheritance, and variance constitute a unified epigenetic process, revealing that intragenic mCG undergoes large, millennia-long epigenetic fluctuations and can therefore mediate evolution on this timescale.

Exploring objective feature sets in constructing the evolution relationship of animal genome sequences

BackgroundExploring evolution regularities of genome sequences and constructing more objective species evolution relationships at the genomic level are high-profile topics. Based on the evolution mechanism of genome sequences proposed in our previous research, we found that only the 8-mers containing CG or TA dinucleotides correlate directly with the evolution of genome sequences, and the relative frequency rather than the actual frequency of these 8-mers is more suitable to characterize the evolution of genome sequences.ResultTherefore, two types of feature sets were obtained, they are the relative frequency sets of CG1 + CG2 8-mers and TA1 + TA2 8-mers. The evolution relationships of mammals and reptiles were constructed by the relative frequency set of CG1 + CG2 8-mers, and two types of evolution relationships of insects were constructed by the relative frequency sets of CG1 + CG2 8-mers and TA1 + TA2 8-mers respectively. Through comparison and analysis, we found that evolution relationships are consistent with the known conclusions. According to the evolution mechanism, we considered that the evolution relationship constructed by CG1 + CG2 8-mers reflects the evolution state of genome sequences in current time, and the evolution relationship constructed by TA1 + TA2 8-mers reflects the evolution state in the early stage.ConclusionOur study provides objective feature sets in constructing evolution relationships at the genomic level.

Targeted DNA Methylation Analysis Facilitates Leukocyte Counts in Dried Blood Samples.

Cell-type specific DNA methylation (DNAm) can be employed to determine the numbers of leukocyte subsets in blood. In contrast to conventional methods for leukocyte counts, which are based on cellular morphology or surface marker protein expression, the cellular deconvolution based on DNAm levels is applicable for frozen or dried blood. Here, we further enhanced targeted DNAm assays for leukocyte counts in clinical application. DNAm profiles of 40 different studies were compiled to identify CG dinucleotides (CpGs) with cell-type specific DNAm using a computational framework, CimpleG. DNAm levels at these CpGs were then measured with digital droplet PCR in venous blood from 160 healthy donors and 150 patients with various hematological disorders. Deconvolution was further validated with venous blood (n = 75) and capillary blood (n = 31) that was dried on Whatman paper or on Mitra microsampling devices. In venous blood, automated cell counting or flow cytometry correlated well with epigenetic estimates of relative leukocyte counts for granulocytes (r = 0.95), lymphocytes (r = 0.97), monocytes (r = 0.82), CD4 T cells (r = 0.84), CD8 T cells (r = 0.94), B cells (r = 0.96), and NK cells (r = 0.72). Similar correlations and precisions were achieved for dried blood samples. Spike-in with a reference plasmid enabled accurate epigenetic estimation of absolute leukocyte counts from dried blood samples, correlating with conventional venous (r = 0.86) and capillary (r = 0.80) blood measurements. The advanced selection of cell-type specific CpGs and utilization of digital droplet PCR analysis provided accurate epigenetic blood counts. Analysis of dried blood facilitates self-sampling with a finger prick, thereby enabling easier accessibility to testing.

Genome-wide impact of cytosine methylation and DNA sequence context on UV-induced CPD formation.

Exposure to ultraviolet (UV) light is the primary etiological agent for skin cancers because UV damages cellular DNA. The most frequent form of UV damage is the cyclobutane pyrimidine dimer (CPD), which consists of covalent linkages between neighboring pyrimidine bases in DNA. In human cells, the 5' position of cytosine bases in CG dinucleotides is frequently methylated, and methylated cytosines in the TP53 tumor suppressor are often sites of mutation hotspots in skin cancers. It has been argued that this is because cytosine methylation promotes UV-induced CPD formation; however, the effects of cytosine methylation on CPD formation are controversial, with conflicting results from previous studies. Here, we use a genome-wide method known as CPD-seq to map UVB- and UVC-induced CPDs across the yeast genome in the presence or absence in vitro methylation by the CpG methyltransferase M.SssI. Our data indicate that cytosine methylation increases UVB-induced CPD formation nearly 2-fold relative to unmethylated DNA, but the magnitude of induction depends on the flanking sequence context. Sequence contexts with a 5' guanine base (e.g., GCCG and GTCG) show the strongest induction due to cytosine methylation, potentially because these sequence contexts are less efficient at forming CPD lesions in the absence of methylation. We show that cytosine methylation also modulates UVC-induced CPD formation, albeit to a lesser extent than UVB. These findings can potentially reconcile previous studies, and define the impact of cytosine methylation on UV damage across a eukaryotic genome.

Tuning Methylation-Dependent Silencing Dynamics by Synthetic Modulation of CpG Density.

Methylation of cytosines in CG dinucleotides (CpGs) within promoters has been shown to lead to gene silencing in mammals in natural contexts. Recently, engineered recruitment of methyltransferases (DNMTs) at specific loci was shown to be sufficient to silence synthetic and endogenous gene expression through this mechanism. A critical parameter for DNA methylation-based silencing is the distribution of CpGs within the target promoter. However, how the number or density of CpGs in the target promoter affects the dynamics of silencing by DNMT recruitment has remained unclear. Here, we constructed a library of promoters with systematically varying CpG content, and analyzed the rate of silencing in response to recruitment of DNMT. We observed a tight correlation between silencing rate and CpG content. Further, methylation-specific analysis revealed a constant accumulation rate of methylation at the promoter after DNMT recruitment. We identified a single CpG site between TATA box and transcription start site (TSS) that accounted for a substantial part of the difference in silencing rates between promoters with differing CpG content, indicating that certain residues play disproportionate roles in controlling silencing. Together, these results provide a library of promoters for synthetic epigenetic and gene regulation applications, as well as insights into the regulatory link between CpG content and silencing rate.

Culture expansion of CAR T cells results in aberrant DNA methylation that is associated with adverse clinical outcome.

Chimeric antigen receptor (CAR) T cells provide new perspectives for treatment of hematological malignancies. Manufacturing of these cellular products includes culture expansion procedures, which may affect cellular integrity and therapeutic outcome. In this study, we investigated culture-associated epigenetic changes in CAR T cells and found continuous gain of DNAm, particularly within genes that are relevant for T cell function. Hypermethylation in many genes, such as TCF7, RUNX1, and TOX, was reflected by transcriptional downregulation. 332 CG dinucleotides (CpGs) showed an almost linear gain in methylation with cell culture time, albeit neighboring CpGs were not coherently regulated on the same DNA strands. An epigenetic signature based on 14 of these culture-associated CpGs predicted cell culture time across various culture conditions. Notably, even in CAR T cell products of similar culture time higher DNAm levels at these CpGs were associated with significantly reduced long-term survival post transfusion. Our data demonstrate that cell culture expansion of CAR T cells evokes DNA hypermethylation at specific sites in the genome and the signature may also reflect loss of potential in CAR T cell products. Hence, reduced cultivation periods are beneficial to avoid dysfunctional methylation programs that seem to be associated with worse therapeutic outcome.

Quantification of hematopoietic stem and progenitor cells by targeted DNA methylation analysis

Hematopoietic stem and progenitor cells (HSPCs) are quantified in daily clinical practice by flow cytometry. In this study, we provide proof of concept that HSPCs can also be estimated by targeted DNA methylation (DNAm) analysis. The DNAm levels at three individual CG dinucleotides (CpG sites) in the genes MYO1D, STK17A, and SP140 correlated with CD34+ cell numbers in mobilized peripheral blood and with blast counts in leukemia. In the future, such epigenetic biomarkers can support the evaluation of stem cell mobilization, HSPC harvesting, or blast count in leukemia.

Extensive de novo activity stabilizes epigenetic inheritance of CG methylation in Arabidopsis transposons

Cytosine methylation within CG dinucleotides (mCG) can be epigenetically inherited over many generations. Such inheritance is thought to be mediated by a semiconservative mechanism that produces binary present/absent methylation patterns. However, we show here that, in Arabidopsis thaliana h1ddm1 mutants, intermediate heterochromatic mCG is stably inherited across many generations and is quantitatively associated with transposon expression. We develop a mathematical model that estimates the rates of semiconservative maintenance failure and de novo methylation at each transposon, demonstrating that mCG can be stably inherited at any level via a dynamic balance of these activities. We find that DRM2-the core methyltransferase of the RNA-directed DNA methylation pathway-catalyzes most of the heterochromatic de novo mCG, with de novo rates orders of magnitude higher than previously thought, whereas chromomethylases make smaller contributions. Our results demonstrate that stable epigenetic inheritance of mCG in plant heterochromatin is enabled by extensive de novo methylation.

Differential methylation of microRNA encoding genes may contribute to high myopia.

Introduction: High myopia (HM), an eye disorder with a refractive error ≤-6.0 diopters, has multifactorial etiology with environmental and genetic factors involved. Recent studies confirm the impact of alterations in DNA methylation and microRNAs (miRNAs) on myopia. Here, we studied the combined aspects evaluating to the role of methylation of miRNA encoding genes in HM. Materials and Methods: From the genome-wide DNA methylation data of 18 Polish children with HM and 18 matched controls, we retrieved differentially methylated CG dinucleotides localized in miRNA encoding genes. Putative target genes of the highest-ranked miRNAs were obtained from the miRDB and included in overrepresentation analyses in the ConsensusPathDB. Expression of target genes was assessed using the RNA sequencing data of retinal ARPE-19 cell line. Results: We identified differential methylation of CG dinucleotides in promoter regions of MIR3621, MIR34C, MIR423 (increased methylation level), and MIR1178, MIRLET7A2, MIR885, MIR548I3, MIR6854, MIR675, MIRLET7C, MIR99A (decreased methylation level) genes. Several targets of these miRNAs, e.g. GNAS, TRAM1, CTNNB1, EIF4B, TENM3 and RUNX were previously associated with myopia/HM/refractive error in Europeans in genome-wide association studies. Overrepresentation analyses of miRNAs' targets revealed enrichment in pathways/processes related to eye structure/function, such as axon guidance, transcription, focal adhesion, and signaling pathways of TGF-β, insulin, MAPK and EGF-EGFR. Conclusion: Differential methylation of indicated miRNA encoding genes might influence their expression and contribute to HM pathogenesis via disrupted regulation of transcription of miRNAs' target genes. Methylation of genes encoding miRNAs may be a new direction in research on both the mechanisms determining HM and non-invasive indicators in diagnostics.

Epigenetic biomarkers to track differentiation of pluripotent stem cells.

Quality control of induced pluripotent stem cells remains a challenge. For validation of the pluripotent state, it is crucial to determine trilineage differentiation potential toward endoderm, mesoderm, and ectoderm. Here, we report GermLayerTracker, a combination of site-specific DNA methylation (DNAm) assays that serve as biomarker for early germ layer specification. CG dinucleotides (CpGs) were identified with characteristic DNAm at pluripotent state and after differentiation into endoderm, mesoderm, and ectoderm. Based on this, a pluripotency score was derived that tracks reprogramming and may indicate differentiation capacity, as well as lineage-specific scores to monitor either directed differentiation or self-organized multilineage differentiation in embryoid bodies. Furthermore, we established pyrosequencing assays for fast and cost-effective analysis. In the future, the GermLayerTracker could be used for quality control of pluripotent cells and to estimate lineage-specific commitment during initial differentiation events.

Identification, expression, and purification of DNA cytosine 5-methyltransferases with short recognition sequences

BackgroundDNA methyltransferases (MTases) are enzymes that induce methylation, one of the representative epigenetic modifications of DNA, and are also useful tools for analyzing epigenomes. However, regarding DNA cytosine 5-methylation, MTases identified so far have drawbacks in that their recognition sequences overlap with those for intrinsic DNA methylation in mammalian cells and/or that the recognition sequence is too long for fine epigenetic mapping. To identify MTases with short recognition sequences that never overlap with the CG dinucleotide, we systematically investigated the 25 candidate enzymes identified using a database search, which showed high similarity to known cytosine 5-MTases recognizing short sequences.ResultsWe identified MTases with six new recognition sequences, including TCTG, CC, CNG, TCG, GCY, and GGCA. Because the recognition sequence never overlapped with the CG dinucleotide, MTases recognizing the CC dinucleotide were promising.ConclusionsIn the current study, we established a procedure for producing active CC-methylating MTases and applied it to nucleosome occupancy and methylome sequencing to prove the usefulness of the enzyme for fine epigenetic mapping. MTases that never overlap with CG dinucleotides would allow us to profile multiple epigenomes simultaneously.

Bisulfite-Free and Single-Base Resolution Detection of Epigenetic DNA Modification of 5-Methylcytosine by Methyltransferase-Directed Labeling with APOBEC3A Deamination Sequencing.

DNA methylation (5-methylcytosine, 5mC) is the most prevalent epigenetic modification that is predominantly found in CG dinucleotides in mammalian genomes. In-depth investigation of the functions of 5mC heavily relies on the quantitative measurement of 5mC at single-base resolution in genomes. Here, we proposed a methyltransferase-directed labeling with APOBEC3A (A3A) deamination sequencing (MLAD-seq) method for the single-base resolution and quantitative detection of 5mC in DNA. In MLAD-seq, a mutant of DNA methyltransferase, M.MpeI-N374K, is utilized to selectively transfer a carboxymethyl group to the 5 position of cytosine in the CG dinucleotide to form 5-carboxymethylcytosine (5camC) using carboxy-S-adenosyl-l-methionine (caSAM) as the cofactor. After A3A treatment, 5camC is resistant to the deamination and base pairs with guanine. Thus, the cytosines in CG sites are read as C in sequencing. On the contrary, the methyl group in 5mC inhibits its carboxymethylcytosine by M.MpeI-N374K and therefore is readily deaminated by A3A to produce thymine that pairs with adenine and is read as T in sequencing. The differential readouts from C and 5mC in the MLAD-seq enable the single-base resolution mapping of 5mC in CG sites in DNA. With the developed MLAD-seq method, we observed the hypermethylation in the promoter region of retinoic acid receptor β (RARB) gene from human nonsmall cell lung tumor tissue. Compared to harsh reaction conditions in bisulfite sequencing that could lead to significant degradation of DNA, the whole procedure of MLAD-seq is carried out under mild conditions, which will avoid DNA damage. Thus, MLAD-seq is more suitable in the scenario where only limited input DNA is available. Taken together, the MLAD-seq offers a valuable tool for bisulfite-free, single-base resolution and quantitative detection of 5mC in limited DNA.

Decreased Levels of DNA Methylation in the PCDHA Gene Cluster as a Risk Factor for Early-Onset High Myopia in Young Children.

PurposeHigh myopia (HM), an eye disorder with at least –6.0 diopters refractive error, has a complex etiology with environmental, genetic, and likely epigenetic factors involved. To complement the DNA methylation assessment in children with HM, we analyzed genes that had significantly lower DNA methylation levels.MethodsThe DNA methylation pattern was studied based on the genome-wide methylation data of 18 Polish children with HM paired with 18 controls. Genes overlapping CG dinucleotides with decreased methylation level in HM cases were assessed by enrichment analyses. From those, genes with CG dinucleotides in promoter regions were further evaluated based on exome sequencing (ES) data of 16 patients with HM from unrelated Polish families, Sanger sequencing data of the studied children, and the RNA sequencing data of human retinal ARPE-19 cells.ResultsThe CG dinucleotide with the most decreased methylation level in cases was identified in a promoter region of PCDHA10 that overlaps intronic regions of PCDHA1–9 of the PCDHA gene cluster in myopia 5q31 locus. Also, two single nucleotide variants, rs200661444, detected in our ES, and rs246073, previously found as associated with a refractive error in a genome-wide association study, were revealed within this gene cluster. Additionally, genes previously linked to ocular phenotypes, myopia-related traits, or loci, including ADAM20, ZFAND6, ETS1, ABHD13, SBSPON, SORBS2, LMOD3, ATXN1, and FARP2, were found to have decreased methylation.ConclusionsAlterations in the methylation pattern of specific CG dinucleotides may be associated with early-onset HM, so this could be used to develop noninvasive biomarkers of HM in children and adolescents.