Clinical Implementation of MetaFusion for Accurate Cancer-Driving Fusion Detection from RNA Sequencing.

Oncogenic fusion genes may be identified from next-generation sequencing data, typically RNA-sequencing. However, in a clinical setting, identifying these alterations is challenging against a background of nonrelevant fusion calls that reduce workflow precision and specificity. Furthermore, although numerous algorithms have been developed to detect fusions in RNA-sequencing, there are variations in their individual sensitivities. Here this problem was addressed by introducing MetaFusion into clinical use. Its utility was illustrated when applied to both whole-transcriptome and targeted sequencing data sets. MetaFusion combines ensemble fusion calls from eight individual fusion-calling algorithms with practice-informed identification of gene fusions that are known to be clinically relevant. In doing so, it allows oncogenic fusions to be identified with near-perfect sensitivity and high precision and specificity, significantly outperforming the individual fusion callers it uses as well as existing clinical-grade software. MetaFusion enhances clinical yield over existing methods and is able to identify fusions that have patient relevance for the purposes of diagnosis, prognosis, and treatment.

Evaluation of single-cell RNAseq labelling algorithms using cancer datasets.

Single-cell RNA sequencing (scRNA-seq) clustering and labelling methods are used to determine precise cellular composition of tissue samples. Automated labelling methods rely on either unsupervised, cluster-based approaches or supervised, cell-based approaches to identify cell types. The high complexity of cancer poses a unique challenge, as tumor microenvironments are often composed of diverse cell subpopulations with unique functional effects that may lead to disease progression, metastasis and treatment resistance. Here, we assess 17 cell-based and 9 cluster-based scRNA-seq labelling algorithms using 8 cancer datasets, providing a comprehensive large-scale assessment of such methods in a cancer-specific context. Using several performance metrics, we show that cell-based methods generally achieved higher performance and were faster compared to cluster-based methods. Cluster-based methods more successfully labelled non-malignant cell types, likely because of a lack of gene signatures for relevant malignant cell subpopulations. Larger cell numbers present in some cell types in training data positively impacted prediction scores for cell-based methods. Finally, we examined which methods performed favorably when trained and tested on separate patient cohorts in scenarios similar to clinical applications, and which were able to accurately label particularly small or under-represented cell populations in the given datasets. We conclude that scPred and SVM show the best overall performances with cancer-specific data and provide further suggestions for algorithm selection. Our analysis pipeline for assessing the performance of cell type labelling algorithms is available in https://github.com/shooshtarilab/scRNAseq-Automated-Cell-Type-Labelling.

Evaluation of single-cell RNA-seq clustering algorithms on cancer tumor datasets.

Tumors are complex biological entities that comprise cell types of different origins, with different mutational profiles and different patterns of transcriptional dysregulation. The exploration of data related to cancer biology requires careful analytical methods to reflect the heterogeneity of cell populations in cancer samples. Single-cell techniques are now able to capture the transcriptional profiles of individual cells. However, the complexity of RNA-seq data, especially in cancer samples, makes it challenging to cluster single-cell profiles into groups that reflect the underlying cell types. We have developed a framework for a systematic examination of single-cell RNA-seq clustering algorithms for cancer data, which uses a range of well-established metrics to generate a unified quality score and algorithm ranking. To demonstrate this framework, we examined clustering performance of 15 different single-cell RNA-seq clustering algorithms on eight different cancer datasets. Our results suggest that the single-cell RNA-seq clustering algorithms fall into distinct groups by performance, with the highest clustering quality on non-malignant cells achieved by three algorithms: Seurat, bigSCale and Cell Ranger. However, for malignant cells, two additional algorithms often reach a better performance, namely Monocle and SC3. Their ability to detect known rare cell types was also among the best, along with Seurat. Our approach and results can be used by a broad audience of practitioners who analyze single-cell transcriptomic data in cancer research.

Generation of DNA Methylation Signatures and Classification of Variants in Rare Neurodevelopmental Disorders Using EpigenCentral.

There are more than 700 genes that encode proteins that function in epigenetic regulation and chromatin modification. Germline variants in these genes (typically heterozygous) are associated with rare neurodevelopmental disorders (NDDs) characterized by growth abnormalities and intellectual and developmental delay. Advancements in next-generation sequencing have dramatically increased the detection of pathogenic sequence variants in genes encoding epigenetic machinery associated with NDDs and, concurrently, the number of clinically uninterpretable variants classified as variants of uncertain significance (VUS). Recently, DNA methylation (DNAm) signatures, disorder-specific patterns of DNAm change, have emerged as a functional tool that provides insights into disorder pathophysiology and can classify pathogenicity of variants in NDDs. To date, our group and others have identified DNAm signatures for more than 60 Mendelian neurodevelopmental disorders caused by variants in genes encoding epigenetic machinery. There is broad interest in both the research and clinical communities to develop and catalog DNAm signatures in rare NDDs, but there are challenges in optimizing study design considerations and availability of platforms that integrate bioinformatics tools with the appropriate statistical framework required to analyze genome-wide DNAm data. We previously published EpigenCentral, a platform for analysis of DNAm data in rare NDDs. In this article, we utilize the published Weaver syndrome dataset to provide step-by-step protocols for using EpigenCentral for exploratory analysis to identify DNAm signatures and for classification of NDD variants. We also provide important considerations for experimental design and interpretation of DNAm results. © 2022 Wiley Periodicals LLC. Basic Protocol 1: Exploratory analysis to identify disorder-specific DNAm signatures Basic Protocol 2: Classification of variants associated with neurodevelopmental disorders.

Genomics4RD: An integrated platform to share Canadian deep-phenotype and multiomic data for international rare disease gene discovery.

Despite recent progress in the understanding of the genetic etiologies of rare diseases (RDs), a significant number remain intractable to diagnostic and discovery efforts. Broad data collection and sharing of information among RD researchers is therefore critical. In 2018, the Care4Rare Canada Consortium launched the project C4R-SOLVE, a subaim of which was to collect, harmonize, and share both retrospective and prospective Canadian clinical and multiomic data. Here, we introduce Genomics4RD, an integrated web-accessible platform to share Canadian phenotypic and multiomic data between researchers, both within Canada and internationally, for the purpose of discovering the mechanisms that cause RDs. Genomics4RD has been designed to standardize data collection and processing, and to help users systematically collect, prioritize, and visualize participant information. Data storage, authorization, and access procedures have been developed in collaboration with policy experts and stakeholders to ensure the trusted and secure access of data by external researchers. The breadth and standardization of data offered by Genomics4RD allows researchers to compare candidate disease genes and variants between participants (i.e., matchmaking) for discovery purposes, while facilitating the development of computational approaches for multiomic data analyses and enabling clinical translation efforts for new genetic technologies in the future.

Anatomy of DNA methylation signatures: Emerging insights and applications.

DNA methylation (DNAm) signatures are unique patterns of DNAm alterations defined for rare disorders caused by pathogenic variants in epigenetic regulatory genes. The potential of DNAm signatures (also known as "episignatures") is just beginning to emerge as there are >300 known epigenetic regulatory genes, ∼100 of which are linked to neurodevelopmental disorders. To date, approximately 50 signatures have been identified, which have proven unexpectedly successful as predictive tools for classifying variants of uncertain significance as pathogenic or benign. The molecular basis of these signatures is poorly understood. Furthermore, their relationships to primary disease pathophysiology have yet to be adequately investigated, despite clear demonstrations of potential connections. There are currently no published guidelines for signature development. As signatures are highly dependent on the samples and methods used to derive them, we propose a framework for consideration in signature development including sample size, statistical parameters, cell type of origin, and the value of detailed clinical and molecular information. We illustrate the relationship between signature output/efficacy and sample size by generating and testing 837 DNAm signatures of Kleefstra syndrome using downsampling analysis. Our findings highlight that no single DNAm signature encompasses all DNAm alterations present in a rare disorder, and that a substandard study design can generate a DNAm signature that misclassifies variants. Finally, we discuss the importance of further investigating DNAm signatures to inform disease pathophysiology and broaden their scope as a functional assay.

MetaFusion: a high-confidence metacaller for filtering and prioritizing RNA-seq gene fusion candidates.

MOTIVATION: Current fusion detection tools use diverse calling approaches and provide varying results, making selection of the appropriate tool challenging. Ensemble fusion calling techniques appear promising; however, current options have limited accessibility and function. RESULTS: MetaFusion is a flexible metacalling tool that amalgamates outputs from any number of fusion callers. Individual caller results are standardized by conversion into the new file type Common Fusion Format. Calls are annotated, merged using graph clustering, filtered and ranked to provide a final output of high-confidence candidates. MetaFusion consistently achieves higher precision and recall than individual callers on real and simulated datasets, and reaches up to 100% precision, indicating that ensemble calling is imperative for high-confidence results. MetaFusion uses FusionAnnotator to annotate calls with information from cancer fusion databases and is provided with a Benchmarking Toolkit to calibrate new callers. AVAILABILITY AND IMPLEMENTATION: MetaFusion is freely available at https://github.com/ccmbioinfo/MetaFusion. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Navigating the Global Protein-Protein Interaction Landscape Using iRefWeb.

iRefWeb is a resource that provides web interface to a large collection of protein-protein interactions aggregated from major primary databases. The underlying data-consolidation process, called iRefIndex, implements a rigorous methodology of identifying redundant protein sequences and integrating disparate data records that reference the same peptide sequences, despite many potential differences in data identifiers across various source databases. iRefWeb offers a unified user interface to all interaction records and associated information collected by iRefIndex, in addition to a number of data filters and visual features that present the supporting evidence. Users of iRefWeb can explore the consolidated landscape of protein-protein interactions, establish the provenance and reliability of each data record, and compare annotations performed by different data curator teams. The iRefWeb portal is freely available at http://wodaklab.org/iRefWeb .

Obsessive-compulsive disorder and attention-deficit/hyperactivity disorder: distinct associations with DNA methylation and genetic variation.

BACKGROUND: A growing body of research has demonstrated associations between specific neurodevelopmental disorders and variation in DNA methylation (DNAm), implicating this molecular mark as a possible contributor to the molecular etiology of these disorders and/or as a novel disease biomarker. Furthermore, genetic risk variants of neurodevelopmental disorders have been found to be enriched at loci associated with DNAm patterns, referred to as methylation quantitative trait loci (mQTLs). METHODS: We conducted two epigenome-wide association studies in individuals with attention-deficit/hyperactivity disorder (ADHD) or obsessive-compulsive disorder (OCD) (aged 4-18 years) using DNA extracted from saliva. DNAm data generated on the Illumina Human Methylation 450 K array were used to examine the interaction between genetic variation and DNAm patterns associated with these disorders. RESULTS: Using linear regression followed by principal component analysis, individuals with the most endorsed symptoms of ADHD or OCD were found to have significantly more distinct DNAm patterns from controls, as compared to all cases. This suggested that the phenotypic heterogeneity of these disorders is reflected in altered DNAm at specific sites. Further investigations of the DNAm sites associated with each disorder revealed that despite little overlap of these DNAm sites across the two disorders, both disorders were significantly enriched for mQTLs within our sample. CONCLUSIONS: Our DNAm data provide insights into the regulatory changes associated with genetic variation, highlighting their potential utility both in directing GWAS and in elucidating the pathophysiology of neurodevelopmental disorders.

EpigenCentral: Portal for DNA methylation data analysis and classification in rare diseases.

Epigenetic processes play a key role in regulating gene expression. Genetic variants that disrupt chromatin-modifying proteins are associated with a broad range of diseases, some of which have specific epigenetic patterns, such as aberrant DNA methylation (DNAm), which may be used as disease biomarkers. While much of the epigenetic research has focused on cancer, there is a paucity of resources devoted to neurodevelopmental disorders (NDDs), which include autism spectrum disorder and many rare, clinically overlapping syndromes. To address this challenge, we created EpigenCentral, a free web resource for biomedical researchers, molecular diagnostic laboratories, and clinical practitioners to perform the interactive classification and analysis of DNAm data related to NDDs. It allows users to search for known disease-associated patterns in their DNAm data, classify genetic variants as pathogenic or benign to assist in molecular diagnostics, or analyze patterns of differential methylation in their data through a simple web form. EpigenCentral is freely available at http://epigen.ccm.sickkids.ca/.

DNA Methylation Signature for EZH2 Functionally Classifies Sequence Variants in Three PRC2 Complex Genes.

Weaver syndrome (WS), an overgrowth/intellectual disability syndrome (OGID), is caused by pathogenic variants in the histone methyltransferase EZH2, which encodes a core component of the Polycomb repressive complex-2 (PRC2). Using genome-wide DNA methylation (DNAm) data for 187 individuals with OGID and 969 control subjects, we show that pathogenic variants in EZH2 generate a highly specific and sensitive DNAm signature reflecting the phenotype of WS. This signature can be used to distinguish loss-of-function from gain-of-function missense variants and to detect somatic mosaicism. We also show that the signature can accurately classify sequence variants in EED and SUZ12, which encode two other core components of PRC2, and predict the presence of pathogenic variants in undiagnosed individuals with OGID. The discovery of a functionally relevant signature with utility for diagnostic classification of sequence variants in EZH2, EED, and SUZ12 supports the emerging paradigm shift for implementation of DNAm signatures into diagnostics and translational research.

Prevalence and Clinical Features of Inflammatory Bowel Diseases Associated With Monogenic Variants, Identified by Whole-Exome Sequencing in 1000 Children at a Single Center.

BACKGROUND & AIMS: A proportion of infants and young children with inflammatory bowel diseases (IBDs) have subtypes associated with a single gene variant (monogenic IBD). We aimed to determine the prevalence of monogenic disease in a cohort of pediatric patients with IBD. METHODS: We performed whole-exome sequencing analyses of blood samples from an unselected cohort of 1005 children with IBD, aged 0-18 years (median age at diagnosis, 11.96 years) at a single center in Canada and their family members (2305 samples total). Variants believed to cause IBD were validated using Sanger sequencing. Biopsies from patients were analyzed by immunofluorescence and histochemical analyses. RESULTS: We identified 40 rare variants associated with 21 monogenic genes among 31 of the 1005 children with IBD (including 5 variants in XIAP, 3 in DOCK8, and 2 each in FOXP3, GUCY2C, and LRBA). These variants occurred in 7.8% of children younger than 6 years and 2.3% of children aged 6-18 years. Of the 17 patients with monogenic Crohn's disease, 35% had abdominal pain, 24% had nonbloody loose stool, 18% had vomiting, 18% had weight loss, and 5% had intermittent bloody loose stool. The 14 patients with monogenic ulcerative colitis or IBD-unclassified received their diagnosis at a younger age, and their most predominant feature was bloody loose stool (78%). Features associated with monogenic IBD, compared to cases of IBD not associated with a single variant, were age of onset younger than 2 years (odds ratio [OR], 6.30; P = .020), family history of autoimmune disease (OR, 5.12; P = .002), extra-intestinal manifestations (OR, 15.36; P < .0001), and surgery (OR, 3.42; P = .042). Seventeen patients had variants in genes that could be corrected with allogeneic hematopoietic stem cell transplantation. CONCLUSIONS: In whole-exome sequencing analyses of more than 1000 children with IBD at a single center, we found that 3% had rare variants in genes previously associated with pediatric IBD. These were associated with different IBD phenotypes, and 1% of the patients had variants that could be potentially corrected with allogeneic hematopoietic stem cell transplantation. Monogenic IBD is rare, but should be considered in analysis of all patients with pediatric onset of IBD.

Correction to: DNA methylation signature is prognostic of choroid plexus tumor aggressiveness.

After publication of the original article [1], authors have requested to add a 'J' as middle name for Richard Gilbertson. Hence, full name should be Richard J Gilbertson.

DNA methylation signature is prognostic of choroid plexus tumor aggressiveness.

BACKGROUND: Histological grading of choroid plexus tumors (CPTs) remains the best prognostic tool to distinguish between aggressive choroid plexus carcinoma (CPC) and the more benign choroid plexus papilloma (CPP) or atypical choroid plexus papilloma (aCPP); however, these distinctions can be challenging. Standard treatment of CPC is very aggressive and often leads to severe damage to the young child's brain. Therefore, it is crucial to distinguish between CPC and less aggressive entities (CPP or aCPP) to avoid unnecessary exposure of the young patient to neurotoxic therapy. To better stratify CPTs, we utilized DNA methylation (DNAm) to identify prognostic epigenetic biomarkers for CPCs. METHODS: We obtained DNA methylation profiles of 34 CPTs using the HumanMethylation450 BeadChip from Illumina, and the data was analyzed using the Illumina Genome Studio analysis software. Validation of differentially methylated CpG sites chosen as biomarkers was performed using pyrosequencing analysis on additional 22 CPTs. Sensitivity testing of the CPC DNAm signature was performed on a replication cohort of 61 CPT tumors obtained from Neuropathology, University Hospital Münster, Germany. RESULTS: Generated genome-wide DNAm profiles of CPTs showed significant differences in DNAm between CPCs and the CPPs or aCPPs. The prediction of clinical outcome could be improved by combining the DNAm profile with the mutational status of TP53. CPCs with homozygous TP53 mutations clustered as a group separate from those carrying a heterozygous TP53 mutation or CPCs with wild type TP53 (TP53-wt) and showed the worst survival outcome. Specific DNAm signatures for CPCs revealed AK1, PER2, and PLSCR4 as potential biomarkers for CPC that can be used to improve molecular stratification for diagnosis and treatment. CONCLUSIONS: We demonstrate that combining specific DNAm signature for CPCs with histological approaches better differentiate aggressive tumors from those that are not life threatening. These findings have important implications for future prognostic risk prediction in clinical disease management.

New insights into DNA methylation signatures: SMARCA2 variants in Nicolaides-Baraitser syndrome.

BACKGROUND: Nicolaides-Baraitser syndrome (NCBRS) is a neurodevelopmental disorder caused by pathogenic sequence variants in SMARCA2 which encodes the catalytic component of the chromatin remodeling BAF complex. Pathogenic variants in genes that encode epigenetic regulators have been associated with genome-wide changes in DNA methylation (DNAm) in affected individuals termed DNAm signatures. METHODS: Genome-wide DNAm was assessed in whole-blood samples from the individuals with pathogenic SMARCA2 variants and NCBRS diagnosis (n = 8) compared to neurotypical controls (n = 23) using the Illumina MethylationEPIC array. Differential methylated CpGs between groups (DNAm signature) were identified and used to generate a model enabling classification variants of uncertain significance (VUS; n = 9) in SMARCA2 as "pathogenic" or "benign". A validation cohort of NCBRS cases (n = 8) and controls (n = 96) demonstrated 100% model sensitivity and specificity. RESULTS: We identified a DNAm signature of 429 differentially methylated CpG sites in individuals with NCBRS. The genes to which these CpG sites map are involved in cell differentiation, calcium signaling, and neuronal function consistent with NCBRS pathophysiology. DNAm model classifications of VUS were concordant with the clinical phenotype; those within the SMARCA2 ATPase/helicase domain classified as "pathogenic". A patient with a mild neurodevelopmental NCBRS phenotype and a VUS distal to the ATPase/helicase domain did not score as pathogenic, clustering away from cases and controls. She demonstrated an intermediate DNAm profile consisting of one subset of signature CpGs with methylation levels characteristic of controls and another characteristic of NCBRS cases; each mapped to genes with ontologies consistent with the patient's unique clinical presentation. CONCLUSIONS: Here we find that a DNAm signature of SMARCA2 pathogenic variants in NCBRS maps to CpGs relevant to disorder pathophysiology, classifies VUS, and is sensitive to the position of the variant in SMARCA2. The patient with an intermediate model score demonstrating a unique genotype-epigenotype-phenotype correlation underscores the potential utility of this signature as a functionally relevant VUS classification system scalable beyond binary "benign" versus "pathogenic" scoring. This is a novel feature of DNAm signatures that could enable phenotypic predictions from genotype data. Our findings also demonstrate that DNAm signatures can be domain-specific, highlighting the precision with which they can reflect genotypic variation.

Don't brush off buccal data heterogeneity.

Buccal epithelial cells are among the most clinically accessible tissues and are increasingly being used to identify epigenetic disease patterns. However, substantial variation in buccal DNA methylation patterns indicates heterogeneity of cell types within and between samples, raising questions of data quality. We systematically estimated cell-type composition for a large collection of buccal and saliva samples from 11 published studies of DNA methylation. In these we identified numerous cases of buccal samples with questionable purity, which may be affected by sampling from individuals with neurodevelopmental disorders, and by the brushes used for sample collection. Further challenges are involved in comparisons with tissues such as saliva, in which buccal component varies widely. We propose a reference-based method of correcting for buccal purity that reduces unwanted variation while preserving cross-tissue differences. Our work demonstrates the wide variation of buccal quality in epigenetic studies and suggests a possible approach to overcome this issue.

Impact of assisted reproduction, infertility, sex and paternal factors on the placental DNA methylome.

Children conceived using Assisted Reproductive Technologies (ART) have a higher incidence of growth and birth defects, attributable in part to epigenetic perturbations. Both ART and germline defects associated with parental infertility could interfere with epigenetic reprogramming events in germ cells or early embryos. Mouse models indicate that the placenta is more susceptible to the induction of epigenetic abnormalities than the embryo, and thus the placental methylome may provide a sensitive indicator of 'at risk' conceptuses. Our goal was to use genome-wide profiling to examine the extent of epigenetic abnormalities in matched placentas from an ART/infertility group and control singleton pregnancies (n = 44/group) from a human prospective longitudinal birth cohort, the Design, Develop, Discover (3D) Study. Principal component analysis revealed a group of ART outliers. The ART outlier group was enriched for females and a subset of placentas showing loss of methylation of several imprinted genes including GNAS, SGCE, KCNQT1OT1 and BLCAP/NNAT. Within the ART group, placentas from pregnancies conceived with in vitro fertilization (IVF)/intracytoplasmic sperm injection (ICSI) showed distinct epigenetic profiles as compared to those conceived with less invasive procedures (ovulation induction, intrauterine insemination). Male factor infertility and paternal age further differentiated the IVF/ICSI group, suggesting an interaction of infertility and techniques in perturbing the placental epigenome. Together, the results suggest that the human placenta is sensitive to the induction of epigenetic defects by ART and/or infertility, and we stress the importance of considering both sex and paternal factors and that some but not all ART conceptuses will be susceptible.

CHARGE and Kabuki Syndromes: Gene-Specific DNA Methylation Signatures Identify Epigenetic Mechanisms Linking These Clinically Overlapping Conditions.

Epigenetic dysregulation has emerged as a recurring mechanism in the etiology of neurodevelopmental disorders. Two such disorders, CHARGE and Kabuki syndromes, result from loss of function mutations in chromodomain helicase DNA-binding protein 7 (CHD7LOF) and lysine (K) methyltransferase 2D (KMT2DLOF), respectively. Although these two syndromes are clinically distinct, there is significant phenotypic overlap. We therefore expected that epigenetically driven developmental pathways regulated by CHD7 and KMT2D would overlap and that DNA methylation (DNAm) alterations downstream of the mutations in these genes would identify common target genes, elucidating a mechanistic link between these two conditions, as well as specific target genes for each disorder. Genome-wide DNAm profiles in individuals with CHARGE and Kabuki syndromes with CHD7LOF or KMT2DLOF identified distinct sets of DNAm differences in each of the disorders, which were used to generate two unique, highly specific and sensitive DNAm signatures. These DNAm signatures were able to differentiate pathogenic mutations in these two genes from controls and from each other. Analysis of the DNAm targets in each gene-specific signature identified both common gene targets, including homeobox A5 (HOXA5), which could account for some of the clinical overlap in CHARGE and Kabuki syndromes, as well as distinct gene targets. Our findings demonstrate how characterization of the epigenome can contribute to our understanding of disease pathophysiology for epigenetic disorders, paving the way for explorations of novel therapeutics.

Genome-wide placental DNA methylation analysis of severely growth-discordant monochorionic twins reveals novel epigenetic targets for intrauterine growth restriction.

BACKGROUND: Intrauterine growth restriction (IUGR), which refers to reduced fetal growth in the context of placental insufficiency, is etiologically heterogeneous. IUGR is associated not only with perinatal morbidity and mortality but also with adult-onset disorders, such as cardiovascular disease and diabetes, posing a major health burden. Placental epigenetic dysregulation has been proposed as one mechanism that causes IUGR; however, the spectrum of epigenetic pathophysiological mechanisms leading to IUGR remains to be elucidated. Monozygotic monochorionic twins are particularly affected by IUGR, in the setting of severe discordant growth. Because monozygotic twins have the same genotype at conception and a shared maternal environment, they provide an ideal model system for studying epigenetic dysregulation of the placenta. RESULTS: We compared genome-wide placental DNA methylation patterns of severely growth-discordant twins to identify novel candidate genes for IUGR. Snap-frozen placental samples for eight severely growth-discordant monozygotic monochorionic twin pairs were obtained at delivery from each twin. A high-resolution DNA methylation array platform was used to identify methylation differences between IUGR and normal twins. Our analysis revealed differentially methylated regions in the promoters of eight genes: DECR1, ZNF300, DNAJA4, CCL28, LEPR, HSPA1A/L, GSTO1, and GNE. The largest methylation differences between the two groups were in the promoters of DECR1 and ZNF300. The significance of these group differences was independently validated by bisulfite pyrosequencing, implicating aberrations in fatty acid beta oxidation and transcriptional regulation, respectively. Further analysis of the array data identified methylation changes most prominently affecting the Wnt and cadherin pathways in the IUGR cohort. CONCLUSIONS: Our results suggest that IUGR in monozygotic twins is associated with impairments in lipid metabolism and transcriptional regulation as well as cadherin and Wnt signaling. We show that monozygotic monochorionic twins discordant for growth provide a useful model to study one type of the epigenetic placental dysregulation that drives IUGR.

The missing indels: an estimate of indel variation in a human genome and analysis of factors that impede detection.

With the development of High-Throughput Sequencing (HTS) thousands of human genomes have now been sequenced. Whenever different studies analyze the same genome they usually agree on the amount of single-nucleotide polymorphisms, but differ dramatically on the number of insertion and deletion variants (indels). Furthermore, there is evidence that indels are often severely under-reported. In this manuscript we derive the total number of indel variants in a human genome by combining data from different sequencing technologies, while assessing the indel detection accuracy. Our estimate of approximately 1 million indels in a Yoruban genome is much higher than the results reported in several recent HTS studies. We identify two key sources of difficulties in indel detection: the insufficient coverage, read length or alignment quality; and the presence of repeats, including short interspersed elements and homopolymers/dimers. We quantify the effect of these factors on indel detection. The quality of sequencing data plays a major role in improving indel detection by HTS methods. However, many indels exist in long homopolymers and repeats, where their detection is severely impeded. The true number of indel events is likely even higher than our current estimates, and new techniques and technologies will be required to detect them.