Comprehensive multi-omics integration identifies differentially active enhancers during human brain development with clinical relevance

Abstract Background Non-coding regulatory elements (NCREs), such as enhancers, play a crucial role in gene regulation, and genetic aberrations in NCREs can lead to human disease, including brain disorders. The human brain is a complex organ that is susceptible to numerous disorders; many of these are caused by genetic changes, but a multitude remain currently unexplained. Understanding NCREs acting during brain development has the potential to shed light on previously unrecognized genetic causes of human brain disease. Despite immense community-wide efforts to understand the role of the non-coding genome and NCREs, annotating functional NCREs remains challenging. Methods Here we performed an integrative computational analysis of virtually all currently available epigenome data sets related to human fetal brain. Results Our in-depth analysis unravels 39,709 differentially active enhancers (DAEs) that show dynamic epigenomic rearrangement during early stages of human brain development, indicating likely biological function. Many of these DAEs are linked to clinically relevant genes, and functional validation of selected DAEs in cell models and zebrafish confirms their role in gene regulation. Compared to enhancers without dynamic epigenomic rearrangement, DAEs are subjected to higher sequence constraints in humans, have distinct sequence characteristics and are bound by a distinct transcription factor landscape. DAEs are enriched for GWAS loci for brain-related traits and for genetic variation found in individuals with neurodevelopmental disorders, including autism. Conclusion This compendium of high-confidence enhancers will assist in deciphering the mechanism behind developmental genetics of human brain and will be relevant to uncover missing heritability in human genetic brain disorders

Hexb enzyme deficiency leads to lysosomal abnormalities in radial glia and microglia in zebrafish brain development

Sphingolipidoses are severe, mostly infantile lysosomal storage disorders (LSDs) caused by defective glycosphingolipid degradation. Two of these sphingolipidoses, Tay Sachs and Sandhoff diseases, are caused by β-Hexosaminidase (HEXB) enzyme deficiency, resulting in ganglioside (GM2) accumulation and neuronal loss. The precise sequence of cellular events preceding, and leading to, neuropathology remains unclear, but likely involves inflammation and lysosomal accumulation of GM2 in multiple cell types. We aimed to determine the consequences of Hexb activity loss for different brain cell types using zebrafish. Hexb deficient zebrafish (hexb−/− ) showed lysosomal abnormalities already early in development both in radial glia, which are the neuronal and glial progenitors, and in microglia. Additionally, at 5 days postfertilization, hexb−/− zebrafish showed reduced locomotor activity. Although specific oligosaccharides accumulate in the adult brain, hexb−/− ) zebrafish are viable and apparently resistant to Hexb deficiency. In all, we identified cellular consequences of loss of Hexb enzyme activity during embryonic brain development, showing early effects on glia, which possibly underlie the behavioral aberrations. Hereby, we identified clues into the contribution of non-neuronal lysosomal abnormalities in LSDs affecting the brain and provide a tool to further study what underlies the relative resistance to Hexb deficiency in vivo

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

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

Background: Hirschsprung 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. Results: We 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. Conclusions: Our 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

New insights into the genetics of primary open-angle glaucoma based on meta-analyses of intraocular pressure and optic disc characteristics

Primary open-angle glaucoma (POAG), the most common optic neuropathy, is a highly heritable disease (h2 = 0.42 ± 0.09). Siblings of POAG cases have a ten-fold increase risk of developing the disease. Intraocular pressure (IOP) and optic nerve head characteristics are used clinically to predict POAG risk. We conducted a genome-wide association meta-analysis of IOP and optic disc parameters and validated our findings in multiple sets of POAG cases and controls. Using imputation to the 1000 genomes (1000G) reference set, we identified 9 new genomic regions associated with vertical cup disc ratio (VCDR) and 1 new region associated with IOP. Additionally, we found 5 novel loci for optic nerve cup area and 6 for disc area. Previously it was assumed that genetic variation influenced POAG either through IOP or via changes to the optic nerve head; here we present evidence that some genomic regions affect both IOP and the disc parameters. We characterized the effect of the novel loci through pathway analysis and found that pathways involved are not entirely distinct as assumed. Further, we identified a novel association between CDKN1A and POAG. Using a zebrafish model we show that six6b (associated with POAG and optic nerve head variation) alters the expression of cdkn1a. In summary, we have identified several novel genes influencing the major clinical risk predictors of POAG and showed that genetic variation in CDKN1A is important in POAG ris

Reverse genetic screen reveals that Il34 facilitates yolk sac macrophage distribution and seeding of the brain

Microglia are brain resident macrophages, which have specialized functions important in brain development and in disease. They colonize the brain in early embryonic stages, but few factors that drive the migration of yolk sac macrophages (YSMs) into the embryonic brain , or regulate their acquisition of specialized properties are currently known.Here, we present a CRISPR/Cas9-based in vivo reverse genetic screening pipeline to identify new microglia regulators using zebrafish. Zebrafish larvae are particularly suitable due to their external development, transparency and conserved microglia features. We targeted putative microglia regulators, by Cas9/gRNA-complex injections, followed by neutral red-based visualization of microglia. Microglia were quantified automatically in 3-day-old larvae using a software tool we called SpotNGlia. We identified that loss of the zebrafish colony stimulating factor 1 receptor (CSF1R) ligand IL34, caused reduced microglia numbers. Previous studies on the role of the IL34 on microglia development in vivo were ambiguous. Our data, and a concurrent paper, show that in zebrafish, il34 is required during the earliest seeding of the brain by microglia. Our data also indicate that Il34 is required for YSM distribution to other organs. Disruption of the other CSF1R ligand, Csf1, did not reduce microglia numbers in mutants, whereas overexpression increased the number of microglia. This shows Csf1 can influence microglia numbers, but might not be essential for the early seeding of the brain. In all, we identified il34 as a modifier of microglia colonization, by affecting distribution of YSMs to target organs, validating our reverse genetic screening pipeline in zebrafish