Mirjam C. G. N. van den Hout

Dynamic long-range chromatin interactions control Myb proto-oncogene transcription during erythroid development

The key haematopoietic regulator Myb is essential for coordinating proliferation and differentiation. ChIP‐Sequencing and Chromosome Conformation Capture (3C)‐Sequencing were used to characterize the structural and protein‐binding dynamics of the Myb locus during erythroid differentiation. In proliferating cells expressing Myb, enhancers within the Myb‐Hbs1l intergenic region were shown to form an active chromatin hub (ACH) containing the Myb promoter and first intron. This first intron was found to harbour the transition site from transcription initiation to elongation, which takes place around a conserved CTCF site. Upon erythroid differentiation, Myb expression is downregulated and the ACH destabilized. We propose a model for Myb activation by distal enhancers dynamically bound by KLF1 and the GATA1/TAL1/LDB1 complex, which primarily function as a transcription elongation element through chromatin looping.

Whole-transcriptome analysis of endothelial to hematopoietic stem cell transition reveals a requirement for Gpr56 in HSC generation.

Using highly sensitive RNAseq to examine the whole transcriptome of enriched aortic hematopoietic stem cells and endothelial cells, the authors find G-protein–coupled receptor, Gpr56, is required to generate the first HSCs during endothelial to hematopoietic cell transition.

Next-generation sequencing-based genome diagnostics across clinical genetics centers: implementation choices and their effects

Implementation of next-generation DNA sequencing (NGS) technology into routine diagnostic genome care requires strategic choices. Instead of theoretical discussions on the consequences of such choices, we compared NGS-based diagnostic practices in eight clinical genetic centers in the Netherlands, based on genetic testing of nine pre-selected patients with cardiomyopathy. We highlight critical implementation choices, including the specific contributions of laboratory and medical specialists, bioinformaticians and researchers to diagnostic genome care, and how these affect interpretation and reporting of variants. Reported pathogenic mutations were consistent for all but one patient. Of the two centers that were inconsistent in their diagnosis, one reported to have found ‘no causal variant’, thereby underdiagnosing this patient. The other provided an alternative diagnosis, identifying another variant as causal than the other centers. Ethical and legal analysis showed that informed consent procedures in all centers were generally adequate for diagnostic NGS applications that target a limited set of genes, but not for exome- and genome-based diagnosis. We propose changes to further improve and align these procedures, taking into account the blurring boundary between diagnostics and research, and specific counseling options for exome- and genome-based diagnostics. We conclude that alternative diagnoses may infer a certain level of ‘greediness’ to come to a positive diagnosis in interpreting sequencing results. Moreover, there is an increasing interdependence of clinic, diagnostics and research departments for comprehensive diagnostic genome care. Therefore, we invite clinical geneticists, physicians, researchers, bioinformatics experts and patients to reconsider their role and position in future diagnostic genome care.

Exome sequencing and functional analyses suggest that SIX6 is a gene involved in an altered proliferation-differentiation balance early in life and optic nerve degeneration at old age.

Primary open-angle glaucoma (POAG) is a hereditary neurodegenerative disease, characterized by optic nerve changes including increased excavation, notching and optic disc hemorrhages. The excavation can be described by the vertical cup-disc ratio (VCDR). Previously, genome-wide significant evidence for the association of rs10483727 in SIX1-SIX6 locus with VCDR and subsequent POAG was found. Using 1000 genomes-based imputation of four independent population-based cohorts in the Netherlands, we identified a missense variant rs33912345 (His141Asn) in SIX6 associated with VCDR (Pmeta = 7.74 × 10(-7), n = 11 473) and POAG (Pmeta = 6.09 × 10(-3), n = 292). Exome sequencing analysis revealed another missense variant rs146737847 (Glu129Lys) also in SIX6 associated with VCDR (P = 5.09 × 10(-3), n = 1208). These two findings point to SIX6 as the responsible gene for the previously reported association signal. Functional characterization of SIX6 in zebrafish revealed that knockdown of six6b led to a small eye phenotype. Histological analysis showed retinal lamination, implying an apparent normal development of the eye, but an underdeveloped lens, and reduced optic nerve diameter. Expression analysis of morphants at 3 dpf showed a 5.5-fold up-regulation of cdkn2b, a cyclin-dependent kinase inhibitor, involved in cell cycle regulation and previously associated with VCDR and POAG in genome-wide association studies (GWASs). Since both six6b and cdkn2b play a key role in cell proliferation, we assessed the proliferative activity in the eye of morphants and found an alteration in the proliferative pattern of retinal cells. Our findings in humans and zebrafish suggest a functional involvement of six6b in early eye development, and open new insights into the genetic architecture of POAG.

Next generation sequencing of SNPs for non-invasive prenatal diagnosis: challenges and feasibility as illustrated by an application to β-thalassaemia.

β-Thalassaemia is one of the most common autosomal recessive single-gene disorder worldwide, with a carrier frequency of 12% in Cyprus. Prenatal tests for at risk pregnancies use invasive methods and development of a non-invasive prenatal diagnostic (NIPD) method is of paramount importance to prevent unnecessary risks inherent to invasive methods. Here, we describe such a method by assessing a modified version of next generation sequencing (NGS) using the Illumina platform, called ‘targeted sequencing’, based on the detection of paternally inherited fetal alleles in maternal plasma. We selected four single-nucleotide polymorphisms (SNPs) located in the β-globin locus with a high degree of heterozygosity in the Cypriot population. Spiked genomic samples were used to determine the specificity of the platform. We could detect the minor alleles in the expected ratio, showing the specificity of the platform. We then developed a multiplexed format for the selected SNPs and analysed ten maternal plasma samples from pregnancies at risk. The presence or absence of the paternal mutant allele was correctly determined in 27 out of 34 samples analysed. With haplotype analysis, NIPD was possible on eight out of ten families. This is the first study carried out for the NIPD of β-thalassaemia using targeted NGS and haplotype analysis. Preliminary results show that NGS is effective in detecting paternally inherited alleles in the maternal plasma.

Whole exome sequencing coupled with unbiased functional analysis reveals new Hirschsprung disease genes.

BackgroundHirschsprung disease (HSCR), which is congenital obstruction of the bowel, results from a failure of enteric nervous system (ENS) progenitors to migrate, proliferate, differentiate, or survive within the distal intestine. Previous studies that have searched for genes underlying HSCR have focused on ENS-related pathways and genes not fitting the current knowledge have thus often been ignored. We identify and validate novel HSCR genes using whole exome sequencing (WES), burden tests, in silico prediction, unbiased in vivo analyses of the mutated genes in zebrafish, and expression analyses in zebrafish, mouse, and human.ResultsWe performed de novo mutation (DNM) screening on 24 HSCR trios. We identify 28 DNMs in 21 different genes. Eight of the DNMs we identified occur in RET, the main HSCR gene, and the remaining 20 DNMs reside in genes not reported in the ENS. Knockdown of all 12 genes with missense or loss-of-function DNMs showed that the orthologs of four genes (DENND3, NCLN, NUP98, and TBATA) are indispensable for ENS development in zebrafish, and these results were confirmed by CRISPR knockout. These genes are also expressed in human and mouse gut and/or ENS progenitors. Importantly, the encoded proteins are linked to neuronal processes shared by the central nervous system and the ENS.ConclusionsOur data open new fields of investigation into HSCR pathology and provide novel insights into the development of the ENS. Moreover, the study demonstrates that functional analyses of genes carrying DNMs are warranted to delineate the full genetic architecture of rare complex diseases.

Genetic variants in RBFOX3 are associated with sleep latency

Time to fall asleep (sleep latency) is a major determinant of sleep quality. Chronic, long sleep latency is a major characteristic of sleep-onset insomnia and/or delayed sleep phase syndrome. In this study we aimed to discover common polymorphisms that contribute to the genetics of sleep latency. We performed a meta-analysis of genome-wide association studies (GWAS) including 2 572 737 single nucleotide polymorphisms (SNPs) established in seven European cohorts including 4242 individuals. We found a cluster of three highly correlated variants (rs9900428, rs9907432 and rs7211029) in the RNA-binding protein fox-1 homolog 3 gene (RBFOX3) associated with sleep latency (P-values=5.77 × 10−08, 6.59 × 10−08 and 9.17 × 10−08). These SNPs were replicated in up to 12 independent populations including 30 377 individuals (P-values=1.5 × 10−02, 7.0 × 10−03 and 2.5 × 10−03; combined meta-analysis P-values=5.5 × 10−07, 5.4 × 10−07 and 1.0 × 10−07). A functional prediction of RBFOX3 based on co-expression with other genes shows that this gene is predominantly expressed in brain (P-value=1.4 × 10−316) and the central nervous system (P-value=7.5 × 10−321). The predicted function of RBFOX3 based on co-expression analysis with other genes shows that this gene is significantly involved in the release cycle of neurotransmitters including gamma-aminobutyric acid and various monoamines (P-values<2.9 × 10−11) that are crucial in triggering the onset of sleep. To conclude, in this first large-scale GWAS of sleep latency we report a novel association of variants in RBFOX3 gene. Further, a functional prediction of RBFOX3 supports the involvement of RBFOX3 with sleep latency.

Paroxysmal exercise-induced dystonia within the phenotypic spectrum of ECHS1 deficiency

ECHS1 encodes a mitochondrial enzyme involved in the degradation of essential amino acids and fatty acids. Recently, ECHS1 mutations were shown to cause a new severe metabolic disorder presenting as Leigh or Leigh‐like syndromes. The objective of this study was to describe a family with 2 siblings affected by different dystonic disorders as a resulting phenotype of ECHS1 mutations.

Deciphering the RNA landscape by RNAome sequencing

Current RNA expression profiling methods rely on enrichment steps for specific RNA classes, thereby not detecting all RNA species in an unperturbed manner. We report strand-specific RNAome sequencing that determines expression of small and large RNAs from rRNA-depleted total RNA in a single sequence run. Since current analysis pipelines cannot reliably analyze small and large RNAs simultaneously, we developed TRAP, Total Rna Analysis Pipeline, a robust interface that is also compatible with existing RNA sequencing protocols. RNAome sequencing quantitatively preserved all RNA classes, allowing cross-class comparisons that facilitates the identification of relationships between different RNA classes. We demonstrate the strength of RNAome sequencing in mouse embryonic stem cells treated with cisplatin. MicroRNA and mRNA expression in RNAome sequencing significantly correlated between replicates and was in concordance with both existing RNA sequencing methods and gene expression arrays generated from the same samples. Moreover, RNAome sequencing also detected additional RNA classes such as enhancer RNAs, anti-sense RNAs, novel RNA species and numerous differentially expressed RNAs undetectable by other methods. At the level of complete RNA classes, RNAome sequencing also identified a specific global repression of the microRNA and microRNA isoform classes after cisplatin treatment whereas all other classes such as mRNAs were unchanged. These characteristics of RNAome sequencing will significantly improve expression analysis as well as studies on RNA biology not covered by existing methods.

Erratum : Next-generation sequencing-based genome diagnostics across clinical genetics centers: Implementation choices and their effects (European Journal of Human Genetics (2015) 23 (1142-1150) DOI:10.1038/ejhg.2014.279)

Terry Vrijenhoek, Ken Kraaijeveld, Martin Elferink, Joep de Ligt, Elcke Kranendonk, Gijs Santen, Isaac J Nijman, Derek Butler, Godelieve Claes, Adalberto Costessi, Wim Dorlijn, Winfried van Eyndhoven, Dicky JJ Halley, Mirjam CGN van den Hout, Steven van Hove, Lennart F Johansson, Jan DH Jongbloed, Rick Kamps, Christel EM Kockx, Bart de Koning, Marjolein Kriek, Ronald Lekanne dit Deprez, Hans Lunstroo, Marcel Mannens, Olaf R Mook, Marcel Nelen, Corrette Ploem, Marco Rijnen, Jasper J Saris, Richard Sinke, Erik Sistermans, Marjon van Slegtenhorst, Frank Sleutels, Nienke van der Stoep, Marianne van Tienhoven, Martijn Vermaat, Maartje Vogel, Quinten Waisfisz, Janneke Marjan Weiss, Arthur van den Wijngaard, Wilbert van Workum, Helger Ijntema, Bert van der Zwaag, Wilfred FJ van IJcken, Johan T den Dunnen, Joris A Veltman, Raoul Hennekam and Edwin Cuppen