Insights into imprinting from parent-of-origin phased methylomes and transcriptomes

Non-coding 886 (nc886/vtRNA2-1), the epigenetic odd duck - implications for future studies

Non-coding 886 (nc886, vtRNA2-1) is the only human polymorphically imprinted gene, in which the methylation status is not determined by genetics. Existing literature regarding the establishment, stability and consequences of the methylation pattern, as well as the nature and function of the nc886 RNAs transcribed from the locus, are contradictory. For example, the methylation status of the locus has been reported to be stable through life and across somatic tissues, but also susceptible to environmental effects. The nature of the produced nc886 RNA(s) has been redefined multiple times, and in carcinogenesis, these RNAs have been reported to have conflicting roles. In addition, due to the bimodal methylation pattern of the nc886 locus, traditional genome-wide methylation analyses can lead to false-positive results, especially in smaller datasets. Herein, we aim to summarize the existing literature regarding nc886, discuss how the characteristics of nc886 give rise to contradictory results, as well as to reinterpret, reanalyse and, where possible, replicate the results presented in the current literature. We also introduce novel findings on how the distribution of the nc886 methylation pattern is associated with the geographical origins of the population and describe the methylation changes in a large variety of human tumours. Through the example of this one peculiar genetic locus and RNA, we aim to highlight issues in the analysis of DNA methylation and non-coding RNAs in general and offer our suggestions for what should be taken into consideration in future analyses.

AIMER: A SNP-independent software for identifying imprinting-like allelic methylated regions from DNA methylome

Genomic imprinting is essential for mammalian growth and embryogenesis. High-throughput bisulfite sequencing accompanied with parental haplotype-specific information allows analysis of imprinted genes and imprinting control regions (ICRs) on a large scale. Currently, although several allelic methylated regions (AMRs) detection software were developed, methods for detecting imprinted AMRs is still limited. Here, we developed a SNP-independent statistical approach, AIMER, to detect imprinting-like AMRs. By using the mouse frontal cortex methylome as input, we demonstrated that AIMER performs very well in detecting known germline ICRs compared with other methods. Furthermore, we found the putative parental AMRs AIMER detected could be distinguished from sequence-dependent AMRs. Finally, we found a novel germline imprinting-like AMR using WGBS data from 17 distinct mouse tissue samples. The results indicate that AIMER is a good choice for detecting imprinting-like (parent-of-origin-dependent) AMRs. We hope this method will be helpful for future genomic imprinting studies. The Python source code for our project is now publicly available on both GitHub (https://github.com/ZhaoLab-TMU/AIMER) and Gitee (https://gitee.com/zhaolab_tmu/AIMER).

Expression of the non-coding RNA nc886 facilitates the development of tyrosine kinase inhibitor resistance in EGFR-mutated non-small-cell lung cancer cells

Treatment of non-small-cell lung cancer (NSCLC) patients possessing EGFR-activating mutations with tyrosine kinase inhibitors (TKIs) can confer an initial promising response. However, TKI resistance inevitably arises. Numerous TKI resistance mechanisms are identified including EGFR secondary mutations, bypass receptor tyrosine kinase (RTK) signaling, and cellular transition e.g. epithelial-mesenchymal transition (EMT). To increase the knowledge of TKI resistance we performed an epigenetic screen to identify small non-coding (nc) genes with DNA methylation alterations in HCC827 NSCLC EGFR-mutated cells with acquired TKI resistance. We analyzed Infinium Methylation EPIC 850K Array data for DNA methylation changes present in both TKI-resistant HCC827 cells with EMT and MET-amplification. Hereby, we identified that the polymorphic maternal imprinted gene nc886 (vtRNA2-1) has a decrease in promoter DNA methylation in TKI-resistant cells. This epigenetic change was associated with an increase in the expression of nc886. The induction of EMT did not affect nc886 expression. CRISPR/Cas9-mediated distortion of the nc886 sequence increased the sensitivity of HCC827 cells towards TKI. Finally, nc886 sequence distortion hindered MET RTK activation and instead was EMT the endpoint TKI resistance mechanism. In conclusion, the expression of nc886 contributes to TKI resistance in the HCC827 NSCLC cell line by supporting cell survival and selection of the endpoint TKI resistance mechanism. We propose DNA methylation and expression changes for nc886 to constitute a novel TKI resistance contributing mechanism in NSCLC.

The correlation between CpG methylation and gene expression is driven by sequence variants

Gene promoter and enhancer sequences are bound by transcription factors and are depleted of methylated CpG sites (cytosines preceding guanines in DNA). The absence of methylated CpGs in these sequences typically correlates with increased gene expression, indicating a regulatory role for methylation. We used nanopore sequencing to determine haplotype-specific methylation rates of 15.3 million CpG units in 7,179 whole-blood genomes. We identified 189,178 methylation depleted sequences where three or more proximal CpGs were unmethylated on at least one haplotype. A total of 77,789 methylation depleted sequences (~41%) associated with 80,503 cis-acting sequence variants, which we termed allele-specific methylation quantitative trait loci (ASM-QTLs). RNA sequencing of 896 samples from the same blood draws used to perform nanopore sequencing showed that the ASM-QTL, that is, DNA sequence variability, drives most of the correlation found between gene expression and CpG methylation. ASM-QTLs were enriched 40.2-fold (95% confidence interval 32.2, 49.9) among sequence variants associating with hematological traits, demonstrating that ASM-QTLs are important functional units in the noncoding genome.

MEG3 in hematologic malignancies: from the role of disease biomarker to therapeutic target

Maternally expressed gene 3 ( MEG3 ) is a noncoding RNA that is known as a tumor suppressor in solid cancers. Recently, a line of studies has emphasized its potential role in hematological malignancies in terms of tumorigenesis, metastasis, and drug resistance. Similar to solid cancers, MEG3 can regulate various cancer hallmarks via sponging miRNA, transcriptional, or posttranslational regulation mechanisms, but may regulate different key elements. In contrast with solid cancers, in some subtypes of leukemia, MEG3 has been found to be upregulated and oncogenic. In this review, we systematically describe the role and underlying mechanisms of MEG3 in multiple types of hematological malignancies. Particularly, we highlight the role of MEG3 in drug resistance and as a novel therapeutic target.

Characterization of the genetic determinants of context-specific DNA methylation in primary monocytes

To better understand inter-individual variation in sensitivity of DNA methylation (DNAm) to immune activity, we characterized effects of inflammatory stimuli on primary monocyte DNAm (n = 190). We find that monocyte DNAm is site-dependently sensitive to lipopolysaccharide (LPS), with LPS-induced demethylation occurring following hydroxymethylation. We identify 7,359 high-confidence immune-modulated CpGs (imCpGs) that differ in genomic localization and transcription factor usage according to whether they represent a gain or loss in DNAm. Demethylated imCpGs are profoundly enriched for enhancers and colocalize to genes enriched for disease associations, especially cancer. DNAm is age associated, and we find that 24-h LPS exposure triggers approximately 6 months of gain in epigenetic age, directly linking epigenetic aging with innate immune activity. By integrating LPS-induced changes in DNAm with genetic variation, we identify 234 imCpGs under local genetic control. Exploring shared causal loci between LPS-induced DNAm responses and human disease traits highlights examples of disease-associated loci that modulate imCpG formation.

A novel imprinted locus on bovine chromosome 18 homologous with human chromosome 16q24.1

Genomic imprinting is an epigenetic regulation mechanism in mammals resulting in the parentally dependent monoallelic expression of genes. Imprinting disorders in humans are associated with several congenital syndromes and cancers and remain the focus of many medical studies. Cattle is a better model organism for investigating human embryo development than mice. Imprinted genes usually cluster on chromosomes and are regulated by different methylation regions (DMRs) located in imprinting control regions that control gene expression in cis. There is an imprinted locus on human chromosome 16q24.1 associated with congenital lethal developmental lung disease in newborns. However, genomic imprinting on bovine chromosome 18, which is homologous with human chromosome 16 has not been systematically studied. The aim of this study was to analyze the allelic expressions of eight genes (CDH13, ATP2C2, TLDC1, COTL1, CRISPLD2, ZDHHC7, KIAA0513, and GSE1) on bovine chromosome 18 and to search the DMRs associated gene allelic expression. Three transcript variants of the ZDHHC7 gene (X1, X2, and X5) showed maternal imprinting in bovine placentas. In addition, the monoallelic expression of X2 and X5 was tissue-specific. Five transcripts of the KIAA0513 gene showed tissue- and isoform-specific monoallelic expression. The CDH13, ATP2C2, and TLDC1 genes exhibited tissue-specific imprinting, however, COTL1, CRISLPLD2, and GSE1 escaped imprinting. Four DMRs, established after fertilization, were found in this region. Two DMRs were located between the ZDHHC7 and KIAA0513 genes, and two were in exon 1 of the CDH13 and ATP2C2 genes, respectively. The results from this study support future studies on the molecular mechanism to regulate the imprinting of candidate genes on bovine chromosome 18.

A comparison of methods for detecting DNA methylation from long-read sequencing of human genomes

BACKGROUND

Long-read sequencing can enable the detection of base modifications, such as CpG methylation, in single molecules of DNA. The most commonly used methods for long-read sequencing are nanopore developed by Oxford Nanopore Technologies (ONT) and single molecule real-time (SMRT) sequencing developed by Pacific Bioscience (PacBio). In this study, we systematically compare the performance of CpG methylation detection from long-read sequencing.

RESULTS

We demonstrate that CpG methylation detection from 7179 nanopore-sequenced DNA samples is highly accurate and consistent with 132 oxidative bisulfite-sequenced (oxBS) samples, isolated from the same blood draws. We introduce quality filters for CpGs that further enhance the accuracy of CpG methylation detection from nanopore-sequenced DNA, while removing at most 30% of CpGs. We evaluate the per-site performance of CpG methylation detection across different genomic features and CpG methylation rates and demonstrate how the latest R10.4 flowcell chemistry and base-calling algorithms improve methylation detection from nanopore sequencing. Additionally, we show how the methylation detection of 50 SMRT-sequenced genomes compares to nanopore sequencing and oxBS.

CONCLUSIONS

This study provides the first systematic comparison of CpG methylation detection tools for long-read sequencing methods. We compare two commonly used computational methods for the detection of CpG methylation in a large number of nanopore genomes, including samples sequenced using the latest R10.4 nanopore flowcell chemistry and 50 SMRT sequenced samples. We provide insights into the strengths and limitations of each sequencing method as well as recommendations for standardization and evaluation of tools designed for genome-scale modified base detection using long-read sequencing.

Roles of endogenous retroviral elements in the establishment and maintenance of imprinted gene expression

DNA methylation (DNAme) has long been recognized as a host defense mechanism, both in the restriction modification systems of prokaryotes as well as in the transcriptional silencing of repetitive elements in mammals. When DNAme was shown to be implicated as a key epigenetic mechanism in the regulation of imprinted genes in mammals, a parallel with host defense mechanisms was drawn, suggesting perhaps a common evolutionary origin. Here we review recent work related to this hypothesis on two different aspects of the developmental imprinting cycle in mammals that has revealed unexpected roles for long terminal repeat (LTR) retroelements in imprinting, both canonical and noncanonical. These two different forms of genomic imprinting depend on different epigenetic marks inherited from the mature gametes, DNAme and histone H3 lysine 27 trimethylation (H3K27me3), respectively. DNAme establishment in the maternal germline is guided by transcription during oocyte growth. Specific families of LTRs, evading silencing mechanisms, have been implicated in this process for specific imprinted genes. In noncanonical imprinting, maternally inherited histone marks play transient roles in transcriptional silencing during preimplantation development. These marks are ultimately translated into DNAme, notably over LTR elements, for the maintenance of silencing of the maternal alleles in the extraembryonic trophoblast lineage. Therefore, LTR retroelements play important roles in both establishment and maintenance of different epigenetic pathways leading to imprinted expression during development. Because such elements are mobile and highly polymorphic among different species, they can be coopted for the evolution of new species-specific imprinted genes.

Allelic reprogramming of chromatin states in human early embryos

The reprogramming of parental epigenomes in human early embryos remains elusive. To what extent the characteristics of parental epigenomes are conserved between humans and mice is currently unknown. Here, we mapped parental haploid epigenomes using human parthenogenetic and androgenetic embryos. Human embryos have a larger portion of genome with parentally specific epigenetic states than mouse embryos. The allelic patterns of epigenetic states for orthologous regions are not conserved between humans and mice. Nevertheless, it is conserved that maternal DNA methylation and paternal H3K27me3 are associated with the repression of two alleles in humans and mice. In addition, for DNA-methylation-dependent imprinting, we report 19 novel imprinted genes and their associated germline differentially methylated regions. Unlike in mice, H3K27me3-dependent imprinting is not observed in human early embryos. Collectively, allele-specific epigenomic reprogramming is different in humans and mice.

The role of long non-coding RNAs and circular RNAs in cervical cancer: modulating miRNA function

Cervical cancer (CC) is a primary global health concern, ranking as the fourth leading cause of cancer-related death in women. Despite advancements in prognosis, long-term outcomes remained poor. Beyond HPV, cofactors like dietary deficiencies, immunosuppression, hormonal contraceptives, co-infections, and genetic variations are involved in CC progression. The pathogenesis of various diseases, including cancer, has brought to light the critical regulatory roles of microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs). The aberrant expression of these miRNAs, lncRNAs, and circRNAs plays a pivotal role in the initiation and progression of CC. This review provides a comprehensive summary of the recent literature regarding the involvement of lncRNAs and circRNAs in modulating miRNA functions in cervical neoplasia and metastasis. Studies have shown that lncRNAs and circRNAs hold great potential as therapeutic agents and innovative biomarkers in CC. However, more clinical research is needed to advance our understanding of the therapeutic benefits of circRNAs and lncRNAs in CC.

Epigenome editing for targeted DNA (de)methylation: a new perspective in modulating gene expression

Traditionally, it has been believed that inheritance is driven as phenotypic variations resulting from changes in DNA sequence. However, this paradigm has been challenged and redefined in the contemporary era of epigenetics. The changes in DNA methylation, histone modification, non-coding RNA biogenesis, and chromatin remodeling play crucial roles in genomic functions and regulation of gene expression. More importantly, some of these changes are inherited to the next generations as a part of epigenetic memory and play significant roles in gene expression. The sum total of all changes in DNA bases, histone proteins, and ncRNA biogenesis constitutes the epigenome. Continuous progress in deciphering epigenetic regulations and the existence of heritable epigenetic/epiallelic variations associated with trait of interest enables to deploy epigenome editing tools to modulate gene expression. DNA methylation marks can be utilized in epigenome editing for the manipulation of gene expression. Initially, genome/epigenome editing technologies relied on zinc-finger protein or transcriptional activator-like effector protein. However, the discovery of clustered regulatory interspaced short palindromic repeats CRISPR)/deadCRISPR-associated protein 9 (dCas9) enabled epigenome editing to be more specific/efficient for targeted DNA (de)methylation. One of the major concerns has been the off-target effects, wherein epigenome editing may unintentionally modify gene/regulatory element which may cause unintended change/harmful effects. Moreover, epigenome editing of germline cell raises several ethical/safety issues. This review focuses on the recent developments in epigenome editing tools/techniques, technological limitations, and future perspectives of this emerging technology in therapeutics for human diseases as well as plant improvement to achieve sustainable developmental goals.

Metastable epialleles in humans

First identified in isogenic mice, metastable epialleles (MEs) are loci where the extent of DNA methylation (DNAm) is variable between individuals but correlates across tissues derived from different germ layers within a given individual. This property, termed systemic interindividual variation (SIV), is attributed to stochastic methylation establishment before germ layer differentiation. Evidence suggests that some putative human MEs are sensitive to environmental exposures in early development. In this review we introduce key concepts pertaining to human MEs, describe methods used to identify MEs in humans, and review their genomic features. We also highlight studies linking DNAm at putative human MEs to early environmental exposures and postnatal (including disease) phenotypes.

Case report: Complete paternal isodisomy on chromosome 18 induces methylation changes in PARD6G-AS1 promotor in a case with arthrogryposis

Uniparental disomy (UPD) is the inheritance of both alleles of a chromosome from only one parent. So far, the detection of UPDs in sequencing data is not well established and a known gap in next-generation sequencing (NGS) diagnostics. By developing a new tool for UPD detection, we re-evaluated an eight-year-old individual presenting with scoliosis, muscle weakness and global developmental delay. Previous panel analysis identified a homozygous likely pathogenic loss-of-function variant in the PIEZO2-gene associated with arthrogryposis (OMIM # 617146). Interestingly, during a re-evaluation process, we identified a region of homozygosity (ROH) covering over 95% of chromosome 18. Segregation and microsatellite analysis within the family revealed that only the father is a heterozygous carrier of the variant in PIEZO2 and confirmed paternal uniparental isodisomy (iUPD) on chromosome 18 in the individual. Further methylation analysis indicated demethylation of the promotor region of PARD6G-AS1, which is described to be maternally imprinted and could possibly influence the individuals' phenotype. Our report describes the first complete iUPD on chromosome 18 and highlights that UPDs can be a cause for homozygous pathogenic variants, which reduces the risk of reoccurrence in case of a new pregnancy in comparison to an autosomal recessive inheritance trait significantly.

Genetics and epidemiology of mutational barcode-defined clonal hematopoiesis

Clonal hematopoiesis (CH) arises when a substantial proportion of mature blood cells is derived from a single hematopoietic stem cell lineage. Using whole-genome sequencing of 45,510 Icelandic and 130,709 UK Biobank participants combined with a mutational barcode method, we identified 16,306 people with CH. Prevalence approaches 50% in elderly participants. Smoking demonstrates a dosage-dependent impact on risk of CH. CH associates with several smoking-related diseases. Contrary to published claims, we find no evidence that CH is associated with cardiovascular disease. We provide evidence that CH is driven by genes that are commonly mutated in myeloid neoplasia and implicate several new driver genes. The presence and nature of a driver mutation alters the risk profile for hematological disorders. Nevertheless, most CH cases have no known driver mutations. A CH genome-wide association study identified 25 loci, including 19 not implicated previously in CH. Splicing, protein and expression quantitative trait loci were identified for CD164 and TCL1A.

DNA 5-methylcytosine detection and methylation phasing using PacBio circular consensus sequencing

Long single-molecular sequencing technologies, such as PacBio circular consensus sequencing (CCS) and nanopore sequencing, are advantageous in detecting DNA 5-methylcytosine in CpGs (5mCpGs), especially in repetitive genomic regions. However, existing methods for detecting 5mCpGs using PacBio CCS are less accurate and robust. Here, we present ccsmeth, a deep-learning method to detect DNA 5mCpGs using CCS reads. We sequence polymerase-chain-reaction treated and M.SssI-methyltransferase treated DNA of one human sample using PacBio CCS for training ccsmeth. Using long (≥10 Kb) CCS reads, ccsmeth achieves 0.90 accuracy and 0.97 Area Under the Curve on 5mCpG detection at single-molecule resolution. At the genome-wide site level, ccsmeth achieves >0.90 correlations with bisulfite sequencing and nanopore sequencing using only 10× reads. Furthermore, we develop a Nextflow pipeline, ccsmethphase, to detect haplotype-aware methylation using CCS reads, and then sequence a Chinese family trio to validate it. ccsmeth and ccsmethphase can be robust and accurate tools for detecting DNA 5-methylcytosines.

The influence of early environment and micronutrient availability on developmental epigenetic programming: lessons from the placenta

DNA methylation is the most commonly studied epigenetic mark in humans, as it is well recognised as a stable, heritable mark that can affect genome function and influence gene expression. Somatic DNA methylation patterns that can persist throughout life are established shortly after fertilisation when the majority of epigenetic marks, including DNA methylation, are erased from the pre-implantation embryo. Therefore, the period around conception is potentially critical for influencing DNA methylation, including methylation at imprinted alleles and metastable epialleles (MEs), loci where methylation varies between individuals but is correlated across tissues. Exposures before and during conception can affect pregnancy outcomes and health throughout life. Retrospective studies of the survivors of famines, such as those exposed to the Dutch Hunger Winter of 1944-45, have linked exposures around conception to later disease outcomes, some of which correlate with DNA methylation changes at certain genes. Animal models have shown more directly that DNA methylation can be affected by dietary supplements that act as cofactors in one-carbon metabolism, and in humans, methylation at birth has been associated with peri-conceptional micronutrient supplementation. However, directly showing a role of micronutrients in shaping the epigenome has proven difficult. Recently, the placenta, a tissue with a unique hypomethylated methylome, has been shown to possess great inter-individual variability, which we highlight as a promising target tissue for studying MEs and mixed environmental exposures. The placenta has a critical role shaping the health of the fetus. Placenta-associated pregnancy complications, such as preeclampsia and intrauterine growth restriction, are all associated with aberrant patterns of DNA methylation and expression which are only now being linked to disease risk later in life.

Comparative analysis reveals epigenomic evolution related to species traits and genomic imprinting in mammals

DNA methylation is an epigenetic modification that plays a crucial role in various regulatory processes, including gene expression regulation, transposable element repression, and genomic imprinting. However, most studies on DNA methylation have been conducted in humans and other model species, whereas the dynamics of DNA methylation across mammals remain poorly explored, limiting our understanding of epigenomic evolution in mammals and the evolutionary impacts of conserved and lineage-specific DNA methylation. Here, we generated and gathered comparative epigenomic data from 13 mammalian species, including two marsupial species, to demonstrate that DNA methylation plays critical roles in several aspects of gene evolution and species trait evolution. We found that the species-specific DNA methylation of promoters and noncoding elements correlates with species-specific traits such as body patterning, indicating that DNA methylation might help establish or maintain interspecies differences in gene regulation that shape phenotypes. For a broader view, we investigated the evolutionary histories of 88 known imprinting control regions across mammals to identify their evolutionary origins. By analyzing the features of known and newly identified potential imprints in all studied mammals, we found that genomic imprinting may function in embryonic development through the binding of specific transcription factors. Our findings show that DNA methylation and the complex interaction between the genome and epigenome have a significant impact on mammalian evolution, suggesting that evolutionary epigenomics should be incorporated to develop a unified evolutionary theory.

Phenome-wide analyses identify an association between the parent-of-origin effects dependent methylome and the rate of aging in humans

BACKGROUND

The variation in the rate at which humans age may be rooted in early events acting through the genomic regions that are influenced by such events and subsequently are related to health phenotypes in later life. The parent-of-origin-effect (POE)-regulated methylome includes regions enriched for genetically controlled imprinting effects (the typical type of POE) and regions influenced by environmental effects associated with parents (the atypical POE). This part of the methylome is heavily influenced by early events, making it a potential route connecting early exposures, the epigenome, and aging. We aim to test the association of POE-CpGs with early and later exposures and subsequently with health-related phenotypes and adult aging.

RESULTS

We perform a phenome-wide association analysis for the POE-influenced methylome using GS:SFHS (Ndiscovery = 5087, Nreplication = 4450). We identify and replicate 92 POE-CpG-phenotype associations. Most of the associations are contributed by the POE-CpGs belonging to the atypical class where the most strongly enriched associations are with aging (DNAmTL acceleration), intelligence, and parental (maternal) smoking exposure phenotypes. A proportion of the atypical POE-CpGs form co-methylation networks (modules) which are associated with these phenotypes, with one of the aging-associated modules displaying increased within-module methylation connectivity with age. The atypical POE-CpGs also display high levels of methylation heterogeneity, fast information loss with age, and a strong correlation with CpGs contained within epigenetic clocks.

CONCLUSIONS

These results identify the association between the atypical POE-influenced methylome and aging and provide new evidence for the "early development of origin" hypothesis for aging in humans.

Haplotype-aware pantranscriptome analyses using spliced pangenome graphs

Pangenomics is emerging as a powerful computational paradigm in bioinformatics. This field uses population-level genome reference structures, typically consisting of a sequence graph, to mitigate reference bias and facilitate analyses that were challenging with previous reference-based methods. In this work, we extend these methods into transcriptomics to analyze sequencing data using the pantranscriptome: a population-level transcriptomic reference. Our toolchain, which consists of additions to the VG toolkit and a standalone tool, RPVG, can construct spliced pangenome graphs, map RNA sequencing data to these graphs, and perform haplotype-aware expression quantification of transcripts in a pantranscriptome. We show that this workflow improves accuracy over state-of-the-art RNA sequencing mapping methods, and that it can efficiently quantify haplotype-specific transcript expression without needing to characterize the haplotypes of a sample beforehand.

Haplotype-aware pantranscriptome analyses using spliced pangenome graphs

Pangenomics is emerging as a powerful computational paradigm in bioinformatics. This field uses population-level genome reference structures, typically consisting of a sequence graph, to mitigate reference bias and facilitate analyses that were challenging with previous reference-based methods. In this work, we extend these methods into transcriptomics to analyze sequencing data using the pantranscriptome: a population-level transcriptomic reference. Our toolchain, which consists of additions to the VG toolkit and a standalone tool, RPVG, can construct spliced pangenome graphs, map RNA sequencing data to these graphs, and perform haplotype-aware expression quantification of transcripts in a pantranscriptome. We show that this workflow improves accuracy over state-of-the-art RNA sequencing mapping methods, and that it can efficiently quantify haplotype-specific transcript expression without needing to characterize the haplotypes of a sample beforehand.

Phenome-wide analysis identifies parent-of-origin effects on the human methylome associated with changes in the rate of aging

Variation in the rate at which humans age may be rooted in early life events acting through genomic regions that are influenced by such events and subsequently are related to health phenotypes in later life. The parent-of-origin-effect (POE)-regulated methylome includes regions either enriched for genetically controlled imprinting effects (the typical type of POE) or atypical POE introduced by environmental effects associated with parents. This part of the methylome is heavily influenced by early life events, making it a potential route connecting early environmental exposures, the epigenome and the rate of aging. Here, we aim to test the association of POE-influenced methylation of CpG dinucleotides (POE-CpG sites) with early and later environmental exposures and subsequently with health-related phenotypes and adult aging phenotypes. We do this by performing phenome-wide association analyses of the POE-influenced methylome using a large family-based population cohort (GS:SFHS, Ndiscovery=5,087, Nreplication=4,450). At the single CpG level, 92 associations of POE-CpGs with phenotypic variation were identified and replicated. Most of the associations were contributed by POE-CpGs belonging to the atypical class and the most strongly enriched associations were with aging (DNAmTL acceleration), intelligence and parental (maternal) smoking exposure phenotypes. We further found that a proportion of the atypical-POE-CpGs formed co-methylation networks (modules) which are associated with these phenotypes, with one of the aging-associated modules displaying increased internal module connectivity (strength of methylation correlation across constituent CpGs) with age. Atypical POE-CpGs also displayed high levels of methylation heterogeneity and epigenetic drift (i.e. information loss with age) and a strong correlation with CpGs contained within epigenetic clocks. These results identified associations between the atypical-POE-influenced methylome and aging and provided new evidence for the "early development of origin" hypothesis for aging in humans.

Inference of putative cell-type-specific imprinted regulatory elements and genes during human neuronal differentiation

Genomic imprinting results in gene expression bias caused by parental chromosome of origin and occurs in genes with important roles during human brain development. However, the cell-type and temporal specificity of imprinting during human neurogenesis is generally unknown. By detecting within-donor allelic biases in chromatin accessibility and gene expression that are unrelated to cross-donor genotype, we inferred imprinting in both primary human neural progenitor cells and their differentiated neuronal progeny from up to 85 donors. We identified 43/20 putatively imprinted regulatory elements (IREs) in neurons/progenitors, and 133/79 putatively imprinted genes in neurons/progenitors. Although 10 IREs and 42 genes were shared between neurons and progenitors, most putative imprinting was only detected within specific cell types. In addition to well-known imprinted genes and their promoters, we inferred novel putative IREs and imprinted genes. Consistent with both DNA methylation-based and H3K27me3-based regulation of imprinted expression, some putative IREs also overlapped with differentially methylated or histone-marked regions. Finally, we identified a progenitor-specific putatively imprinted gene overlapping with copy number variation that is associated with uniparental disomy-like phenotypes. Our results can therefore be useful in interpreting the function of variants identified in future parent-of-origin association studies.

Parent-of-origin detection and chromosome-scale haplotyping using long-read DNA methylation sequencing and Strand-seq

Hundreds of loci in human genomes have alleles that are methylated differentially according to their parent of origin. These imprinted loci generally show little variation across tissues, individuals, and populations. We show that such loci can be used to distinguish the maternal and paternal homologs for all human autosomes without the need for the parental DNA. We integrate methylation-detecting nanopore sequencing with the long-range phase information in Strand-seq data to determine the parent of origin of chromosome-length haplotypes for both DNA sequence and DNA methylation in five trios with diverse genetic backgrounds. The parent of origin was correctly inferred for all autosomes with an average mismatch error rate of 0.31% for SNVs and 1.89% for insertions or deletions (indels). Because our method can determine whether an inherited disease allele originated from the mother or the father, we predict that it will improve the diagnosis and management of many genetic diseases.

Phasing DNA Methylation

Haplotyping enables the study of allele-specific events. Heterozygous variants, primarily single nucleotide variants (SNVs), enable the assignment of the paternal and maternal origin of the chromosomes and are widely employed to phase sequencing reads to their haplotype of origin. Certain long-read technologies enable the detection of both the DNA sequence and DNA modifications. These long reads and their inherent methylation information are suitable for genome-wide haplotyping and allele-specific DNA methylation analysis. Here, we describe the workflow to phase reads and DNA methylation using nanopore sequencing.

Parent-of-Origin inference for biobanks

Identical genetic variations can have different phenotypic effects depending on their parent of origin. Yet, studies focusing on parent-of-origin effects have been limited in terms of sample size due to the lack of parental genomes or known genealogies. We propose a probabilistic approach to infer the parent-of-origin of individual alleles that does not require parental genomes nor prior knowledge of genealogy. Our model uses Identity-By-Descent sharing with second- and third-degree relatives to assign alleles to parental groups and leverages chromosome X data in males to distinguish maternal from paternal groups. We combine this with robust haplotype inference and haploid imputation to infer the parent-of-origin for 26,393 UK Biobank individuals. We screen 99 phenotypes for parent-of-origin effects and replicate the discoveries of 6 GWAS studies, confirming signals on body mass index, type 2 diabetes, standing height and multiple blood biomarkers, including the known maternal effect at the MEG3/DLK1 locus on platelet phenotypes. We also report a novel maternal effect at the TERT gene on telomere length, thereby providing new insights on the heritability of this phenotype. All our summary statistics are publicly available to help the community to better characterize the molecular mechanisms leading to parent-of-origin effects and their implications for human health.

Tissue- and ethnicity-independent hypervariable DNA methylation states show evidence of establishment in the early human embryo

We analysed DNA methylation data from 30 datasets comprising 3474 individuals, 19 tissues and 8 ethnicities at CpGs covered by the Illumina450K array. We identified 4143 hypervariable CpGs ('hvCpGs') with methylation in the top 5% most variable sites across multiple tissues and ethnicities. hvCpG methylation was influenced but not determined by genetic variation, and was not linked to probe reliability, epigenetic drift, age, sex or cell heterogeneity effects. hvCpG methylation tended to covary across tissues derived from different germ-layers and hvCpGs were enriched for proximity to ERV1 and ERVK retrovirus elements. hvCpGs were also enriched for loci previously associated with periconceptional environment, parent-of-origin-specific methylation, and distinctive methylation signatures in monozygotic twins. Together, these properties position hvCpGs as strong candidates for studying how stochastic and/or environmentally influenced DNA methylation states which are established in the early embryo and maintained stably thereafter can influence life-long health and disease.

Genome-wide detection of imprinted differentially methylated regions using nanopore sequencing

Imprinting is a critical part of normal embryonic development in mammals, controlled by defined parent-of-origin (PofO) differentially methylated regions (DMRs) known as imprinting control regions. Direct nanopore sequencing of DNA provides a means to detect allelic methylation and to overcome the drawbacks of methylation array and short-read technologies. Here, we used publicly available nanopore sequencing data for 12 standard B-lymphocyte cell lines to acquire the genome-wide mapping of imprinted intervals in humans. Using the sequencing data, we were able to phase 95% of the human methylome and detect 94% of the previously well-characterized, imprinted DMRs. In addition, we found 42 novel imprinted DMRs (16 germline and 26 somatic), which were confirmed using whole-genome bisulfite sequencing (WGBS) data. Analysis of WGBS data in mouse (Mus musculus), rhesus monkey (Macaca mulatta), and chimpanzee (Pan troglodytes) suggested that 17 of these imprinted DMRs are conserved. Some of the novel imprinted intervals are within or close to imprinted genes without a known DMR. We also detected subtle parental methylation bias, spanning several kilobases at seven known imprinted clusters. At these blocks, hypermethylation occurs at the gene body of expressed allele(s) with mutually exclusive H3K36me3 and H3K27me3 allelic histone marks. These results expand upon our current knowledge of imprinting and the potential of nanopore sequencing to identify imprinting regions using only parent-offspring trios, as opposed to the large multi-generational pedigrees that have previously been required.

TARSII and CARSII: Two approaches for SNP-independent identification of germline differentially methylated regions in mammals

Identifying germline differentially methylated regions (DMRs) in outbred mammals remains a challenge because of difficulty in obtaining single-nucleotide polymorphisms (SNPs). To overcome this difficulty, we developed two computational approaches, TARSII and CARSII, which allow accurate prediction of germline DMRs from DNA methylomes independent of SNPs. Furthermore, we introduce an easy and quick way to validate the predicted germline DMRs with allelic DNA methylation using CGmapTools. Collectively, our strategy can greatly facilitate de novo identification of germline DMRs in outbred mammals. For complete details on the use and execution of this protocol, please refer to Chu et al. (2021).

Methylation pattern of polymorphically imprinted nc886 is not conserved across mammalia

Background

In humans, the nc886 locus is a polymorphically imprinted metastable epiallele. Periconceptional conditions have an effect on the methylation status of nc886, and further, this methylation status is associated with health outcomes in later life, in line with the Developmental Origins of Health and Disease (DOHaD) hypothesis. Animal models would offer opportunities to study the associations between periconceptional conditions, nc886 methylation status and metabolic phenotypes further. Thus, we set out to investigate the methylation pattern of the nc886 locus in non-human mammals.

Data

We obtained DNA methylation data from the data repository GEO for mammals, whose nc886 gene included all three major parts of nc886 and had sequency similarity of over 80% with the human nc886. Our final sample set consisted of DNA methylation data from humans, chimpanzees, bonobos, gorillas, orangutangs, baboons, macaques, vervets, marmosets and guinea pigs.

Results

In human data sets the methylation pattern of nc886 locus followed the expected bimodal distribution, indicative of polymorphic imprinting. In great apes, we identified a unimodal DNA methylation pattern with 50% methylation level in all individuals and in all subspecies. In Old World monkeys, the between individual variation was greater and methylation on average was close to 60%. In guinea pigs the region around the nc886 homologue was non-methylated. Results obtained from the sequence comparison of the CTCF binding sites flanking the nc886 gene support the results on the DNA methylation data.

Conclusions

Our results indicate that unlike in humans, nc886 is not a polymorphically imprinted metastable epiallele in non-human primates or in guinea pigs, thus implying that animal models are not applicable for nc886 research. The obtained data suggests that the nc886 region may be classically imprinted in great apes, and potentially also in Old World monkeys, but not in guinea pigs.

New insights regarding origin of monosomy occurrence in early developing embryos as demonstrated in preimplantation genetic testing

Introduction

Analyses of miscarriage products indicate that the majority of aneuploidies in early developing embryos derive from errors occurring during maternal meiosis and the paternal contribution is less than 10%. Our aim was to assess the aneuploidy (mainly monosmies) frequencies at the earliest stages of embryo development, 3 days following fertilization during In vitro fertilization (IVF) treatments and to elucidate their parental origin. Later, we compared monosomies rates of day 3 to those of day 5 as demonstrated from Preimplantation Genetic Testing for Structural chromosomal Rearrangement (PGT-SR) results.

Methods

For a retrospective study, we collected data of 210 Preimplantation Genetic Testing for Monogenic Disorder (PGT-M) cycles performed between years 2008 and 2019.This study includes 2083 embryos, of 113 couples. It also included 432 embryos from 90 PGT-SR cycles of other 45 patients, carriers of balanced translocations. Defining the parental origin of aneuploidy in cleavage stage embryos was based on haplotypes analysis of at least six informative markers flanking the analyzed gene. For comprehensive chromosomal screening (CCS), chromosomal microarray (CMA) and next generation sequencing (NGS) was used.

Results

We inspected haplotype data of 40 genomic regions, flanking analyzed genes located on 9 different chromosomes.151 (7.2%) embryos presented numerical alterations in the tested chromosomes. We found similar paternal and maternal contribution to monosomy at cleavage stage. We demonstrated paternal origin in 51.5% of the monosomy, and maternal origin in 48.5% of the monosomies cases.

Conclusion

In our study, we found equal parental contribution to monosomies in cleavage-stage embryos. Comparison to CCS analyses of PGT-SR patients revealed a lower rate of monosomy per chromosome in embryos at day 5 of development. This is in contrast to the maternal dominancy described in studies of early miscarriage. Mitotic errors and paternal involvement in chemical pregnancies and IVF failure should be re-evaluated. Our results show monosomies are relatively common and may play a role in early development of ART embryos.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13039-022-00582-5.

Genomic Imprinting in the New Omics Era: A Model for Systems-Level Approaches

Genomic imprinting represents a noteworthy inheritance mechanism leading to allele-specific regulations dependent of the parental origin. Imprinted loci are especially involved in essential mammalian functions related to growth, development and behavior. In this mini-review, we first offer a summary of current representations associated with genomic imprinting through key results of the three last decades. We then outline new perspectives allowed by the spread of new omics technologies tackling various interacting levels of imprinting regulations, including genomics, transcriptomics and epigenomics. We finally discuss the expected contribution of new omics data to unresolved big questions in the field.

Environmentally sensitive hotspots in the methylome of the early human embryo

In humans, DNA methylation marks inherited from gametes are largely erased following fertilisation, prior to construction of the embryonic methylome. Exploiting a natural experiment of seasonal variation including changes in diet and nutritional status in rural Gambia, we analysed three datasets covering two independent child cohorts and identified 259 CpGs showing consistent associations between season of conception (SoC) and DNA methylation. SoC effects were most apparent in early infancy, with evidence of attenuation by mid-childhood. SoC-associated CpGs were enriched for metastable epialleles, parent-of-origin-specific methylation and germline differentially methylated regions, supporting a periconceptional environmental influence. Many SoC-associated CpGs overlapped enhancers or sites of active transcription in H1 embryonic stem cells and fetal tissues. Half were influenced but not determined by measured genetic variants that were independent of SoC. Environmental 'hotspots' providing a record of environmental influence at periconception constitute a valuable resource for investigating epigenetic mechanisms linking early exposures to lifelong health and disease.

Lifestyle and Genetic Factors Modify Parent-of-Origin Effects on the Human Methylome

BACKGROUND

parent-of-origin effects (POE) play important roles in complex disease and thus understanding their regulation and associated molecular and phenotypic variation are warranted. Previous studies mainly focused on the detection of genomic regions or phenotypes regulated by POE. Understanding whether POE may be modified by environmental or genetic exposures is important for understanding of the source of POE-associated variation, but only a few case studies addressing modifiable POE exist.

METHODS

in order to understand this high order of POE regulation, we screened 101 genetic and environmental factors such as 'predicted mRNA expression levels' of DNA methylation/imprinting machinery genes and environmental exposures. POE-mQTL-modifier interaction models were proposed to test the potential of these factors to modify POE at DNA methylation using data from Generation Scotland: The Scottish Family Health Study(N=2315).

FINDINGS

a set of vulnerable/modifiable POE-CpGs were identified (modifiable-POE-regulated CpGs, N=3). Four factors, 'lifetime smoking status' and 'predicted mRNA expression levels' of TET2, SIRT1 and KDM1A, were found to significantly modify the POE on the three CpGs in both discovery and replication datasets. We further identified plasma protein and health-related phenotypes associated with the methylation level of one of the identified CpGs.

INTERPRETATION

the modifiable POE identified here revealed an important yet indirect path through which genetic background and environmental exposures introduce their effect on DNA methylation, motivating future comprehensive evaluation of the role of these modifiers in complex diseases.

FUNDING

NSFC (81971270),H2020-MSCA-ITN(721815), Wellcome (204979/Z/16/Z,104036/Z/14/Z), MRC (MC_UU_00007/10, MC_PC_U127592696), CSO (CZD/16/6,CZB/4/276, CZB/4/710), SFC (HR03006), EUROSPAN (LSHG-CT-2006-018947), BBSRC (BBS/E/D/30002276), SYSU, Arthritis Research UK, NHLBI, NIH.

The Role of Long Non-coding RNAs in Human Imprinting Disorders: Prospective Therapeutic Targets

Genomic imprinting is a term used for an intergenerational epigenetic inheritance and involves a subset of genes expressed in a parent-of-origin-dependent way. Imprinted genes are expressed preferentially from either the paternally or maternally inherited allele. Long non-coding RNAs play essential roles in regulating this allele-specific expression. In several well-studied imprinting clusters, long non-coding RNAs have been found to be essential in regulating temporal- and spatial-specific establishment and maintenance of imprinting patterns. Furthermore, recent insights into the epigenetic pathological mechanisms underlying human genomic imprinting disorders suggest that allele-specific expressed imprinted long non-coding RNAs serve as an upstream regulator of the expression of other protein-coding or non-coding imprinted genes in the same cluster. Aberrantly expressed long non-coding RNAs result in bi-allelic expression or silencing of neighboring imprinted genes. Here, we review the emerging roles of long non-coding RNAs in regulating the expression of imprinted genes, especially in human imprinting disorders, and discuss three strategies targeting the central long non-coding RNA UBE3A-ATS for the purpose of developing therapies for the imprinting disorders Prader-Willi syndrome and Angelman syndrome. In summary, a better understanding of long non-coding RNA-related mechanisms is key to the development of potential therapeutic targets for human imprinting disorders.

Analysis of developmental imprinting dynamics in primates using SNP-free methods to identify imprinting defects in cloned placenta

Our knowledge of genomic imprinting in primates is lagging behind that of mice largely because of the difficulties of allelic analyses in outbred animals. To understand imprinting dynamics in primates, we profiled transcriptomes, DNA methylomes, and H3K27me3 in uniparental monkey embryos. We further developed single-nucleotide-polymorphism (SNP)-free methods, TARSII and CARSII, to identify germline differentially methylated regions (DMRs) in somatic tissues. Our comprehensive analyses showed that allelic DNA methylation, but not H3K27me3, is a major mark that correlates with paternal-biasedly expressed genes (PEGs) in uniparental monkey embryos. Interestingly, primate germline DMRs are different from PEG-associated DMRs in early embryos and are enriched in placenta. Strikingly, most placenta-specific germline DMRs are lost in placenta of cloned monkeys. Collectively, our study establishes SNP-free germline DMR identification methods, defines developmental imprinting dynamics in primates, and demonstrates imprinting defects in cloned monkey placenta, which provides important clues for improving primate cloning.

Parent-of-origin effects on nuclear chromatin organization and behavior in a Drosophila model for Williams-Beuren Syndrome

Prognosis of neuropsychiatric disorders in progeny requires consideration of individual (1) parent-of-origin effects (POEs) relying on (2) the nerve cell nuclear 3D chromatin architecture and (3) impact of parent-specific miRNAs. Additionally, the shaping of cognitive phenotypes in parents depends on both learning acquisition and forgetting, or memory erasure. These processes are independent and controlled by different signal cascades: the first is cAMPdependent, the second relies on actin remodeling by small GTPase Rac1 - LIMK1 (LIM-kinase 1). Simple experimental model systems such as Drosophila help probe the causes and consequences leading to human neurocognitive pathologies. Recently, we have developed a Drosophila model for Williams-Beuren Syndrome (WBS): a mutant agn^ts3 of the agnostic locus (X:11AB) harboring the dlimk1 gene. The agn^ts3 mutation drastically increases the frequency of ectopic contacts (FEC) in specific regions of intercalary heterochromatin, suppresses learning/memory and affects locomotion. As is shown in this study, the polytene X chromosome bands in reciprocal hybrids between agn^ts3 and the wild type strain Berlin are heterogeneous in modes of FEC regulation depending either on maternal or paternal gene origin. Bioinformatic analysis reveals that FEC between X:11AB and the other X chromosome bands correlates with the occurrence of short (~30 bp) identical DNA fragments partly homologous to Drosophila 372-bp satellite DNA repeat. Although learning acquisition in a conditioned courtship suppression paradigm is similar in hybrids, the middle-term memory formation shows patroclinic inheritance. Seemingly, this depends on changes in miR-974 expression. Several parameters of locomotion demonstrate heterosis. Our data indicate that the agn^ts3 locus is capable of trans-regulating gene activity via POEs on the chromatin nuclear organization, thereby affecting behavior.

Targeted RNA-seq improves efficiency, resolution, and accuracy of allele specific expression for human term placentas

Genomic imprinting is an epigenetic mechanism that results in allele-specific expression (ASE) based on the parent of origin. It is known to play a role in the prenatal and postnatal allocation of maternal resources in mammals. ASE detected by whole transcriptome RNA-seq (wht-RNAseq) has been widely used to analyze imprinted genes using reciprocal crosses in mice to generate large numbers of informative SNPs. Studies in humans are more challenging due to the paucity of SNPs and the poor preservation of RNA in term placentas and other tissues. Targeted RNA-seq (tar-RNAseq) can potentially mitigate these challenges by focusing sequencing resources on the regions of interest in the transcriptome. Here, we compared tar-RNAseq and wht-RNAseq in a study of ASE in known imprinted genes in placental tissue collected from a healthy human cohort in Mali, West Africa. As expected, tar-RNAseq substantially improved the coverage of SNPs. Compared to wht-RNAseq, tar-RNAseq produced on average four times more SNPs in twice as many genes per sample and read depth at the SNPs increased fourfold. In previous research on humans, discordant ASE values for SNPs of the same gene have limited the ability to accurately quantify ASE. We show that tar-RNAseq reduces this limitation as it unexpectedly increased the concordance of ASE between SNPs of the same gene, even in cases of degraded RNA. Studies aimed at discovering associations between individual variation in ASE and phenotypes in mammals and flowering plants will benefit from the improved power and accuracy of tar-RNAseq.

Methylation status of nc886 epiallele reflects periconceptional conditions and is associated with glucose metabolism through nc886 RNAs

BACKGROUND

Non-coding RNA 886 (nc886) is coded from a maternally inherited metastable epiallele. We set out to investigate the determinants and dynamics of the methylation pattern at the nc886 epiallele and how this methylation status associates with nc886 RNA expression. Furthermore, we investigated the associations between the nc886 methylation status or the levels of nc886 RNAs and metabolic traits in the YFS and KORA cohorts. The association between nc886 epiallele methylation and RNA expression was also validated in induced pluripotent stem cell (iPSC) lines.

RESULTS

We confirm that the methylation status of the nc886 epiallele is mostly binomial, with individuals displaying either a non- or hemi-methylated status, but we also describe intermediately and close to fully methylated individuals. We show that an individual's methylation status is associated with the mother's age and socioeconomic status, but not with the individual's own genetics. Once established, the methylation status of the nc886 epiallele remains stable for at least 25 years. This methylation status is strongly associated with the levels of nc886 non-coding RNAs in serum, blood, and iPSC lines. In addition, nc886 methylation status associates with glucose and insulin levels during adolescence but not with the indicators of glucose metabolism or the incidence of type 2 diabetes in adulthood. However, the nc886-3p RNA levels also associate with glucose metabolism in adulthood.

CONCLUSIONS

These results indicate that nc886 metastable epiallele methylation is tuned by the periconceptional conditions and it associates with glucose metabolism through the expression of the ncRNAs coded in the epiallele region.

Methylation status of nc886 epiallele reflects periconceptional conditions and is associated with glucose metabolism through nc886 RNAs

BACKGROUND

Non-coding RNA 886 (nc886) is coded from a maternally inherited metastable epiallele. We set out to investigate the determinants and dynamics of the methylation pattern at the nc886 epiallele and how this methylation status associates with nc886 RNA expression. Furthermore, we investigated the associations between the nc886 methylation status or the levels of nc886 RNAs and metabolic traits in the YFS and KORA cohorts. The association between nc886 epiallele methylation and RNA expression was also validated in induced pluripotent stem cell (iPSC) lines.

RESULTS

We confirm that the methylation status of the nc886 epiallele is mostly binomial, with individuals displaying either a non- or hemi-methylated status, but we also describe intermediately and close to fully methylated individuals. We show that an individual's methylation status is associated with the mother's age and socioeconomic status, but not with the individual's own genetics. Once established, the methylation status of the nc886 epiallele remains stable for at least 25 years. This methylation status is strongly associated with the levels of nc886 non-coding RNAs in serum, blood, and iPSC lines. In addition, nc886 methylation status associates with glucose and insulin levels during adolescence but not with the indicators of glucose metabolism or the incidence of type 2 diabetes in adulthood. However, the nc886-3p RNA levels also associate with glucose metabolism in adulthood.

CONCLUSIONS

These results indicate that nc886 metastable epiallele methylation is tuned by the periconceptional conditions and it associates with glucose metabolism through the expression of the ncRNAs coded in the epiallele region.

Epigenome editing by CRISPR/Cas9 in clinical settings: possibilities and challenges

Epigenome editing is a promising approach for both basic research and clinical application. With the convergence of techniques from different fields, regulating gene expression artificially becomes possible. From a clinical point of view, targeted epigenome editing by CRISPR/Cas9 of disease-related genes offers novel therapeutic avenues for many diseases. In this review, we summarize the EpiEffectors used in epigenome editing by CRISPR/Cas9, current applications of epigenome editing and progress made in this field. Moreover, application challenges such as off-target effects, inefficient delivery, stability and immunogenicity are discussed. In conclusion, epigenome editing by CRISPR/Cas9 has broad prospects in the clinic, and future work will promote the application of this technology.

Distinction between the effects of parental and fetal genomes on fetal growth

Birth weight is a common measure of fetal growth that is associated with a range of health outcomes. It is directly affected by the fetal genome and indirectly by the maternal genome. We performed genome-wide association studies on birth weight in the genomes of the child and parents and further analyzed birth length and ponderal index, yielding a total of 243 fetal growth variants. We clustered those variants based on the effects of transmitted and nontransmitted alleles on birth weight. Out of 141 clustered variants, 22 were consistent with parent-of-origin-specific effects. We further used haplotype-specific polygenic risk scores to directly test the relationship between adult traits and birth weight. Our results indicate that the maternal genome contributes to increased birth weight through blood-glucose-raising alleles while blood-pressure-raising alleles reduce birth weight largely through the fetal genome.

Genome-wide survey of parent-of-origin-specific associations across clinical traits derived from electronic health records

Parent-of-origin (PoO) effects refer to the differential phenotypic impacts of genetic variants dependent on their parental inheritance due to imprinting. While PoO effects can influence complex traits, they may be poorly captured by models that do not differentiate the parental origin of the variant. The aim of this study was to conduct a genome-wide screen for PoO effects on a broad range of clinical traits derived from electronic health records (EHR) in the DiscovEHR study enriched with familial relationships. Using pairwise kinship estimates from genetic data and demographic data, we identified 22,051 offspring among 134,049 individuals in the DiscovEHR study. PoO of ~9 million variants was assigned in the offspring by comparing offspring and parental genotypes and haplotypes. We then performed genome-wide PoO association analyses across 154 quantitative and 611 binary traits extracted from EHR. Of the 732 significant PoO associations identified (p < 5 × 10^-8), we attempted to replicate 274 PoO associations in the UK Biobank study with 5,015 offspring and replicated 9 PoO associations (p < 0.05). In summary, our study implements a bioinformatic and statistical approach to examine PoO effects genome-wide in a large population study enriched with familial relationships and systematically characterizes PoO effects on hundreds of clinical traits derived from EHR. Our results suggest that, while the statistical power to detect PoO effects remains modest yet, accurately modeling PoO effects has the potential to find new associations that may have been missed by the standard additive model, further enhancing the mechanistic understanding of genetic influence on complex traits.

Imprinted Gene Expression and Function of the Dopa Decarboxylase Gene in the Developing Heart

Dopa decarboxylase (DDC) synthesizes serotonin in the developing mouse heart where it is encoded by Ddc_exon1a, a tissue-specific paternally expressed imprinted gene. Ddc_exon1a shares an imprinting control region (ICR) with the imprinted, maternally expressed (outside of the central nervous system) Grb10 gene on mouse chromosome 11, but little else is known about the tissue-specific imprinted expression of Ddc_exon1a. Fluorescent immunostaining localizes DDC to the developing myocardium in the pre-natal mouse heart, in a region susceptible to abnormal development and implicated in congenital heart defects in human. Ddc_exon1a and Grb10 are not co-expressed in heart nor in brain where Grb10 is also paternally expressed, despite sharing an ICR, indicating they are mechanistically linked by their shared ICR but not by Grb10 gene expression. Evidence from a Ddc_exon1a gene knockout mouse model suggests that it mediates the growth of the developing myocardium and a thinning of the myocardium is observed in a small number of mutant mice examined, with changes in gene expression detected by microarray analysis. Comparative studies in the human developing heart reveal a paternal expression bias with polymorphic imprinting patterns between individual human hearts at DDC_EXON1a, a finding consistent with other imprinted genes in human.

Genome-wide detection of cytosine methylation by single molecule real-time sequencing

Single molecule real-time (SMRT) sequencing theoretically offers the opportunity to directly assess certain base modifications of native DNA molecules without any prior chemical/enzymatic conversions and PCR amplification, using kinetic signals of a DNA polymerase. However, the kinetic signal changes caused by 5mC modification are extremely subtle. Hence, the robust genome-wide measurement of 5mC modification has not been achieved. We enhanced 5mC detection using SMRT sequencing by holistically analyzing kinetic signals of a DNA polymerase and sequence context for every base within a measurement window. We employed a convolutional neural network to train a methylation classification model, leading to genome-wide 5mC detection. The sensitivity and specificity reached 90% and 94%, with a 99% correlation of overall methylation level with bisulfite sequencing.

5-Methylcytosine (5mC) is an important type of epigenetic modification. Bisulfite sequencing (BS-seq) has limitations, such as severe DNA degradation. Using single molecule real-time sequencing, we developed a methodology to directly examine 5mC. This approach holistically examined kinetic signals of a DNA polymerase (including interpulse duration and pulse width) and sequence context for every nucleotide within a measurement window, termed the holistic kinetic (HK) model. The measurement window of each analyzed double-stranded DNA molecule comprised 21 nucleotides with a cytosine in a CpG site in the center. We used amplified DNA (unmethylated) and M.SssI-treated DNA (methylated) (M.SssI being a CpG methyltransferase) to train a convolutional neural network. The area under the curve for differentiating methylation states using such samples was up to 0.97. The sensitivity and specificity for genome-wide 5mC detection at single-base resolution reached 90% and 94%, respectively. The HK model was then tested on human–mouse hybrid fragments in which each member of the hybrid had a different methylation status. The model was also tested on human genomic DNA molecules extracted from various biological samples, such as buffy coat, placental, and tumoral tissues. The overall methylation levels deduced by the HK model were well correlated with those by BS-seq (r = 0.99; P < 0.0001) and allowed the measurement of allele-specific methylation patterns in imprinted genes. Taken together, this methodology has provided a system for simultaneous genome-wide genetic and epigenetic analyses.

Features and mechanisms of canonical and noncanonical genomic imprinting

Genomic imprinting is the monoallelic expression of a gene based on parent of origin and is a consequence of differential epigenetic marking between the male and female germlines. Canonically, genomic imprinting is mediated by allelic DNA methylation. However, recently it has been shown that maternal H3K27me3 can result in DNA methylation-independent imprinting, termed "noncanonical imprinting." In this review, we compare and contrast what is currently known about the underlying mechanisms, the role of endogenous retroviral elements, and the conservation of canonical and noncanonical genomic imprinting.

Megabase-scale methylation phasing using nanopore long reads and NanoMethPhase

The ability of nanopore sequencing to simultaneously detect modified nucleotides while producing long reads makes it ideal for detecting and phasing allele-specific methylation. However, there is currently no complete software for detecting SNPs, phasing haplotypes, and mapping methylation to these from nanopore sequence data. Here, we present NanoMethPhase, a software tool to phase 5-methylcytosine from nanopore sequencing. We also present SNVoter, which can post-process nanopore SNV calls to improve accuracy in low coverage regions. Together, these tools can accurately detect allele-specific methylation genome-wide using nanopore sequence data with low coverage of about ten-fold redundancy.

A meta-analysis uncovers the first sequence variant conferring risk of Bell's palsy

Bell's palsy is the most common cause of unilateral facial paralysis and is defined as an idiopathic and acute inability to control movements of the facial muscles on the affected side. While the pathogenesis remains unknown, previous studies have implicated post-viral inflammation and resulting compression of the facial nerve. Reported heritability estimates of 4-14% suggest a genetic component in the etiology and an autosomal dominant inheritance has been proposed. Here, we report findings from a meta-analysis of genome-wide association studies uncovering the first unequivocal association with Bell's palsy (rs9357446-A; P = 6.79 × 10^-23, OR = 1.23; N_cases = 4714, N_controls = 1,011,520). The variant also confers risk of intervertebral disc disorders (P = 2.99 × 10^-11, OR = 1.04) suggesting a common pathogenesis in part or a true pleiotropy.

Inducible uniparental chromosome disomy to probe genomic imprinting at single-cell level in brain and beyond

Genomic imprinting is an epigenetic mechanism that results in parental allele-specific expression of ~1% of all genes in mouse and human. Imprinted genes are key developmental regulators and play pivotal roles in many biological processes such as nutrient transfer from the mother to offspring and neuronal development. Imprinted genes are also involved in human disease, including neurodevelopmental disorders, and often occur in clusters that are regulated by a common imprint control region (ICR). In extra-embryonic tissues ICRs can act over large distances, with the largest surrounding Igf2r spanning over 10 million base-pairs. Besides classical imprinted expression that shows near exclusive maternal or paternal expression, widespread biased imprinted expression has been identified mainly in brain. In this review we discuss recent developments mapping cell type specific imprinted expression in extra-embryonic tissues and neocortex in the mouse. We highlight the advantages of using an inducible uniparental chromosome disomy (UPD) system to generate cells carrying either two maternal or two paternal copies of a specific chromosome to analyze the functional consequences of genomic imprinting. Mosaic Analysis with Double Markers (MADM) allows fluorescent labeling and concomitant induction of UPD sparsely in specific cell types, and thus to over-express or suppress all imprinted genes on that chromosome. To illustrate the utility of this technique, we explain how MADM-induced UPD revealed new insights about the function of the well-studied Cdkn1c imprinted gene, and how MADM-induced UPDs led to identification of highly cell type specific phenotypes related to perturbed imprinted expression in the mouse neocortex. Finally, we give an outlook on how MADM could be used to probe cell type specific imprinted expression in other tissues in mouse, particularly in extra-embryonic tissues.