Raf Winand

In vitro screening of embryos by whole-genome sequencing: now, in the future or never?

STUDY QUESTION What are the analytical and clinical validity and the clinical utility of in vitro screening of embryos by whole-genome sequencing? SUMMARY ANSWER At present there are still many limitations in terms of analytical and clinical validity and utility and many ethical questions remain. WHAT IS KNOWN ALREADY Whole-genome sequencing of IVF/ICSI embryos is technically possible. Many loss-of-function mutations exist in the general population without serious effects on the phenotype of the individual. Moreover, annotations of genes and the reference genome are still not 100% correct. STUDY DESIGN, SIZE, DURATION We used publicly available samples from the 1000 Genomes project and Complete Genomics, together with 42 samples from in-house research samples of parents from trios to investigate the presence of loss-of-function mutations in healthy individuals. PARTICIPANTS/MATERIALS, SETTING, METHODS In the samples, we looked for mutations in genes that are associated with a selection of severe Mendelian disorders with a known molecular basis. We looked for mutations predicted to be damaging by PolyPhen and SIFT and for mutations annotated as disease causing in Human Genome Mutation Database (HGMD). MAIN RESULTS AND THE ROLE OF CHANCE More than 40% of individuals who can be considered healthy have mutations that are predicted to be damaging in genes associated with severe Mendelian disorders or are annotated as disease causing. LIMITATIONS, REASONS FOR CAUTION The analysis relies on current knowledge and databases are continuously updated to reflect our increasing knowledge about the genome. In the process of our analysis several updates were already made. WIDER IMPLICATIONS OF THE FINDINGS At this moment it is not advisable to use whole-genome sequencing as a tool to set up health profiles to select embryos for transfer. We also raise some ethical questions that have to be addressed before this technology can be used for embryo selection. TRIAL REGISTRATION NUMBER N/A.

DNA methylation analysis of Homeobox genes implicates HOXB7 hypomethylation as risk factor for neural tube defects.

Neural tube defects (NTDs) are common birth defects of complex etiology. Though family- and population-based studies have confirmed a genetic component, the responsible genes for NTDs are still largely unknown. Based on the hypothesis that folic acid prevents NTDs by stimulating methylation reactions, epigenetic factors, such as DNA methylation, are predicted to be involved in NTDs. Homeobox (HOX) genes play a role in spinal cord development and are tightly regulated in a spatiotemporal and collinear manner, partly by epigenetic modifications. We have quantified DNA methylation for the different HOX genes by subtracting values from a genome-wide methylation analysis using leukocyte DNA from 10 myelomeningocele (MMC) patients and 6 healthy controls. From the 1575 CpGs profiled for the 4 HOX clusters, 26 CpGs were differentially methylated (P-value < 0.05; β-difference > 0.05) between MMC patients and controls. Seventy-seven percent of these CpGs were located in the HOXA and HOXB clusters, with the most profound difference for 3 CpGs within the HOXB7 gene body. A validation case-control study including 83 MMC patients and 30 unrelated healthy controls confirmed a significant association between MMC and HOXB7 hypomethylation (-14.4%; 95% CI: 11.9–16.9%; P-value < 0.0001) independent of the MTHFR 667C>T genotype. Significant HOXB7 hypomethylation was also present in 12 unaffected siblings, each related to a MMC patient, suggestive of an epigenetic change induced by the mother. The inclusion of a neural tube formation model using zebrafish showed that Hoxb7a overexpression but not depletion resulted in deformed body axes with dysmorphic neural tube formation. Our results implicate HOXB7 hypomethylation as risk factor for NTDs and highlight the importance for future genome-wide DNA methylation analyses without preselecting candidate pathways.

Assessing transmissibility of HIV-1 drug resistance mutations from treated and from drug-naive individuals.

Objectives:Surveillance drug resistance mutations (SDRMs) in drug-naive patients are typically used to survey HIV-1-transmitted drug resistance (TDR). We test here how SDRMs in patients failing treatment, the original source of TDR, contribute to assessing TDR, transmissibility and transmission source of SDRMs. Design:This is a retrospective observational study analyzing a Portuguese cohort of HIV-1-infected patients. Methods:The prevalence of SDRMs to protease inhibitors, nucleoside reverse transcriptase inhibitors (NRTIs) and nonnucleoside reverse transcriptase inhibitors (NNRTIs) in drug-naive and treatment-failing patients was measured for 3554 HIV-1 subtype B patients. Transmission ratio (prevalence in drug-naive/prevalence in treatment-failing patients), average viral load and robust linear regression with outlier detection (prevalence in drug-naive versus in treatment-failing patients) were analyzed and used to interpret transmissibility. Results:Prevalence of SDRMs in drug-naive and treatment-failing patients were linearly correlated, but some SDRMs were classified as outliers – above (PRO: D30N, N88D/S, L90 M, RT: G190A/S/E) or below (RT: M184I/V) expectations. The normalized regression slope was 0.073 for protease inhibitors, 0.084 for NRTIs and 0.116 for NNRTIs. Differences between SDRMs transmission ratios were not associated with differences in viral loads. Conclusion:The significant linear correlation between prevalence of SDRMs in drug-naive and in treatment-failing patients indicates that the prevalence in treatment-failing patients can be useful to predict levels of TDR. The slope is a cohort-dependent estimate of rate of TDR per drug class and outlier detection reveals comparative persistence of SDRMs. Outlier SDRMs with higher transmissibility are more persistent and more likely to have been acquired from drug-naive patients. Those with lower transmissibility have faster reversion dynamics after transmission and are associated with acquisition from treatment-failing patients.

Methylome analysis for spina bifida shows SOX18 hypomethylation as a risk factor with evidence for a complex (epi)genetic interplay to affect neural tube development

BackgroundNeural tube defects (NTDs) are severe congenital malformations that arise from failure of neurulation during early embryonic development. The molecular basis underlying most human NTDs still remains largely unknown. Based on the hypothesis that folic acid prevents NTDs by stimulating methylation reactions, DNA methylation changes could play a role in NTDs. We performed a methylome analysis for patients with myelomeningocele (MMC). Using a candidate CpG analysis for HOX genes, a significant association between HOXB7 hypomethylation and MMC was found.MethodsIn the current study, we analyzed leukocyte methylome data of ten patients with MMC and six controls using Illumina Methylation Analyzer and WateRmelon R-packages and performed validation studies using larger MMC and control cohorts with Sequenom EpiTYPER.ResultsThe methylome analysis showed 75 CpGs in 45 genes that are significantly differentially methylated in MMC patients. CpG-specific methylation differences were next replicated for the top six candidate genes ABAT, CNTNAP1, SLC1A6, SNED1, SOX18, and TEPP but only for the SOX18 locus a significant overall hypomethylation was observed (P value = 0.0003). Chemically induced DNA demethylation in HEK cells resulted in SOX18 hypomethylation and increased expression. Injection of sox18 mRNA in zebrafish resulted in abnormal neural tube formation. Quantification of DNA methylation for the SOX18 locus was also determined for five families where parents had normal methylation values compared to significant lower values for both the MMC as their non-affected child. SOX18 methylation studies were performed for a MMC patient with a paternally inherited chromosomal deletion that includes BMP4. The patient showed extreme SOX18 hypomethylation similar to his healthy mother while his father had normal methylation values.ConclusionsThis is the first genome-wide methylation study in leukocytes for patients with NTDs. We report SOX18 as a novel MMC risk gene but our findings also suggest that SOX18 hypomethylation must interplay with environmental and (epi)genetic factors to cause NTDs. Further studies are needed that combine methylome data with next-generation sequencing approaches to unravel NTD etiology.

dendsort: modular leaf ordering methods for dendrogram representations in R

Dendrograms are graphical representations of binary tree structures resulting from agglomerative hierarchical clustering. In Life Science, a cluster heat map is a widely accepted visualization technique that utilizes the leaf order of a dendrogram to reorder the rows and columns of the data table. The derived linear order is more meaningful than a random order, because it groups similar items together. However, two consecutive items can be quite dissimilar despite proximity in the order. In addition, there are 2 n-1 possible orderings given n input elements as the orientation of clusters at each merge can be flipped without affecting the hierarchical structure. We present two modular leaf ordering methods to encode both the monotonic order in which clusters are merged and the nested cluster relationships more faithfully in the resulting dendrogram structure. We compare dendrogram and cluster heat map visualizations created using our heuristics to the default heuristic in R and seriation-based leaf ordering methods. We find that our methods lead to a dendrogram structure with global patterns that are easier to interpret, more legible given a limited display space, and more insightful for some cases. The implementation of methods is available as an R package, named ”dendsort”, from the CRAN package repository. Further examples, documentations, and the source code are available at [https://bitbucket.org/biovizleuven/dendsort/].

Nanopore sequencing technology: a new route for the fast detection of unauthorized GMO

In order to strengthen the current genetically modified organism (GMO) detection system for unauthorized GMO, we have recently developed a new workflow based on DNA walking to amplify unknown sequences surrounding a known DNA region. This DNA walking is performed on transgenic elements, commonly found in GMO, that were earlier detected by real-time PCR (qPCR) screening. Previously, we have demonstrated the ability of this approach to detect unauthorized GMO via the identification of unique transgene flanking regions and the unnatural associations of elements from the transgenic cassette. In the present study, we investigate the feasibility to integrate the described workflow with the MinION Next-Generation-Sequencing (NGS). The MinION sequencing platform can provide long read-lengths and deal with heterogenic DNA libraries, allowing for rapid and efficient delivery of sequences of interest. In addition, the ability of this NGS platform to characterize unauthorized and unknown GMO without any a priori knowledge has been assessed.