Bulk and single-cell transcriptomics identify gene signatures of stem cell-derived NK cell donors with superior cytolytic activity.

Allogeneic natural killer (NK) cell therapies are a valuable treatment option for cancer, given their remarkable safety and favorable efficacy profile. Although the use of allogeneic donors allows for off-the-shelf and timely patient treatment, intrinsic interindividual differences put clinical efficacy at risk. The identification of donors with superior anti-tumor activity is essential to ensure the success of adoptive NK cell therapies. Here, we investigated the heterogeneity of 10 umbilical cord blood stem cell-derived NK cell batches. First, we evaluated the donors' cytotoxic potential against tumor cell lines from solid and hematological cancer indications, to distinguish a group of superior, "excellent" killers (4/10), compared with "good" killers (6/10). Next, bulk and single-cell RNA sequencing, performed at different stages of NK differentiation, revealed distinct transcriptomic features of the two groups. Excellent donors showed an enrichment in cytotoxicity pathways and a depletion of myeloid traits, linked to the presence of a larger population of effector-like NK cells early on during differentiation. Consequently, we defined a multi-factorial gene expression signature able to predict the donors' cytotoxic potential. Our study contributes to the identification of key traits of superior NK cell batches, supporting the development of efficacious NK therapeutics and the achievement of durable anti-tumor responses.

BMP4 and Temozolomide Synergize in the Majority of Patient-Derived Glioblastoma Cultures.

One of the main causes of poor prognoses in patient with glioblastoma (GBM) is drug resistance to current standard treatment, which includes chemoradiation and adjuvant temozolomide (TMZ). In addition, the concept of cancer stem cells provides new insights into therapy resistance and management also in GBM and glioblastoma stem cell-like cells (GSCs), which might contribute to therapy resistance. Bone morphogenetic protein-4 (BMP4) stimulates astroglial differentiation of GSCs and thereby reduces their self-renewal capacity. Exposure of GSCs to BMP4 may also sensitize these cells to TMZ. A recent phase I trial has shown that local delivery of BMP4 is safe, but a large variation in survival is seen in these treated patients and in features of their cultured tumors. We wanted to combine TMZ and BMP4 (TMZ + BMP4) therapy and assess the inter-tumoral variability in response to TMZ + BMP4 in patient-derived GBM cultures. A phase II trial could then benefit a larger group of patients than those treated with BMP4 only. We first show that simultaneous treatment with TMZ + BMP4 is more effective than sequential treatment. Second, when applying our optimized treatment protocol, 70% of a total of 20 GBM cultures displayed TMZ + BMP4 synergy. This combination induces cellular apoptosis and does not inhibit cell proliferation. Comparative bulk RNA-sequencing indicates that treatment with TMZ + BMP4 eventually results in decreased MAPK signaling, in line with previous evidence that increased MAPK signaling is associated with resistance to TMZ. Based on these results, we advocate further clinical trial research to test patient benefit and validate pathophysiological hypothesis.

A three-dimensional vessel-on-chip model to study Puumala orthohantavirus pathogenesis.

Puumala orthohantavirus (PUUV) infection in humans can result in hemorrhagic fever with renal syndrome. Endothelial cells (ECs) are primarily infected with increased vascular permeability as a central aspect of pathogenesis. Historically, most studies included ECs cultured under static two-dimensional (2D) conditions, thereby not recapitulating the physiological environment due to their lack of flow and inherent pro-inflammatory state. Here, we present a high-throughput model for culturing primary human umbilical vein ECs in 3D vessels-on-chip in which we compared host responses of these ECs to those of static 2D-cultured ECs on a transcriptional level. The phenotype of ECs in vessels-on-chip more closely resembled the in vivo situation due to higher similarity in expression of genes encoding described markers for disease severity and coagulopathy, including IDO1, LGALS3BP, IL6 and PLAT, and more diverse endothelial-leukocyte interactions in the context of PUUV infection. In these vessels-on-chip, PUUV infection did not directly increase vascular permeability, but increased monocyte adhesion. This platform can be used for studying pathogenesis and assessment of possible therapeutics for other endotheliotropic viruses even in high biocontainment facilities.

Human Pluripotent Stem Cell-Derived Astrocyte Functionality Compares Favorably with Primary Rat Astrocytes.

Astrocytes are essential for the formation and maintenance of neural networks. However, a major technical challenge for investigating astrocyte function and disease-related pathophysiology has been the limited ability to obtain functional human astrocytes. Despite recent advances in human pluripotent stem cell (hPSC) techniques, primary rodent astrocytes remain the gold standard in coculture with human neurons. We demonstrate that a combination of leukemia inhibitory factor (LIF) and bone morphogenetic protein-4 (BMP4) directs hPSC-derived neural precursor cells to a highly pure population of astroglia in 28 d. Using single-cell RNA sequencing, we confirm the astroglial identity of these cells and highlight profound transcriptional adaptations in cocultured hPSC-derived astrocytes and neurons, consistent with their further maturation. In coculture with human neurons, multielectrode array recordings revealed robust network activity of human neurons in a coculture with hPSC-derived or rat astrocytes [3.63 ± 0.44 min-1 (hPSC-derived), 2.86 ± 0.64 min-1 (rat); p = 0.19]. In comparison, we found increased spike frequency within network bursts of human neurons cocultured with hPSC-derived astrocytes [56.31 ± 8.56 Hz (hPSC-derived), 24.77 ± 4.04 Hz (rat); p < 0.01], and whole-cell patch-clamp recordings revealed an increase of postsynaptic currents [2.76 ± 0.39 Hz (hPSC-derived), 1.07 ± 0.14 Hz (rat); p < 0.001], consistent with a corresponding increase in synapse density [14.90 ± 1.27/100 µm2 (hPSC-derived), 8.39 ± 0.63/100 µm2 (rat); p < 0.001]. Taken together, we show that hPSC-derived astrocytes compare favorably with rat astrocytes in supporting human neural network activity and maturation, providing a fully human platform for investigating astrocyte function and neuronal-glial interactions.

Unraveling the impact of AXIN1 mutations on HCC development: Insights from CRISPR/Cas9 repaired AXIN1-mutant liver cancer cell lines.

BACKGROUND: Hepatocellular carcinoma (HCC) is a highly aggressive liver cancer with significant morbidity and mortality rates. AXIN1 is one of the top-mutated genes in HCC, but the mechanism by which AXIN1 mutations contribute to HCC development remains unclear. METHODS: In this study, we utilized CRISPR/Cas9 genome editing to repair AXIN1-truncated mutations in five HCC cell lines. RESULTS: For each cell line we successfully obtained 2-4 correctly repaired clones, which all show reduced ß-catenin signaling accompanied with reduced cell viability and colony formation. Although exposure of repaired clones to Wnt3A-conditioned medium restored ß-catenin signaling, it did not or only partially recover their growth characteristics, indicating the involvement of additional mechanisms. Through RNA-sequencing analysis, we explored the gene expression patterns associated with repaired AXIN1 clones. Except for some highly-responsive ß-catenin target genes, no consistent alteration in gene/pathway expression was observed. This observation also applies to the Notch and YAP/TAZ-Hippo signaling pathways, which have been associated with AXIN1-mutant HCCs previously. The AXIN1-repaired clones also cannot confirm a recent observation that AXIN1 is directly linked to YAP/TAZ protein stability and signaling. CONCLUSIONS: Our study provides insights into the effects of repairing AXIN1 mutations on ß-catenin signaling, cell viability, and colony formation in HCC cell lines. However, further investigations are necessary to understand the complex mechanisms underlying HCC development associated with AXIN1 mutations.

Species-specific responses during Seoul orthohantavirus infection in human and rat lung microvascular endothelial cells.

Seoul orthohantavirus (SEOV) is a rat-borne zoonotic virus that is transmitted via inhalation of aerosolized infectious excreta, and can cause hemorrhagic fever with renal syndrome (HFRS) in humans worldwide. In rats, SEOV predominantly exists as a persistent infection in the absence of overt clinical signs. Lack of disease in rats is attributed to downregulation of pro-inflammatory and upregulation of regulatory host responses. As lung microvascular endothelial cells (LMECs) represent a primary target of infection in both human and rats, infections in these cells provide a unique opportunity to study the central role of LMECs in the dichotomy between pathogenicity in both species. In this study, host responses to SEOV infection in primary human and rat LMECs were directly compared on a transcriptional level. As infection of rat LMECs was more efficient than human LMECs, the majority of anti-viral defense responses were observed earlier in rat LMECs. Most prominently, SEOV-induced processes in both species included responses to cytokine stimulus, negative regulation of innate immune responses, responses to type I and II interferons, regulation of pattern recognition receptor signaling and MHC-I signaling. However, over time, in the rat LMECs, responses shifted from an anti-viral state towards a more immunotolerant state displayed by a PD-L1, B2M-, JAK2-focused interaction network aiding in negative regulation of cytotoxic CD8-positive T cell activation. This suggests a novel mechanism by which species-specific orthohantavirus-induced endothelium and T cell crosstalk may play a crucial role in the development of acute disease in humans and persistence in rodents.

Comparison of Single Cell Transcriptome Sequencing Methods: Of Mice and Men.

Single cell RNAseq has been a big leap in many areas of biology. Rather than investigating gene expression on a whole organism level, this technology enables scientists to get a detailed look at rare single cells or within their cell population of interest. The field is growing, and many new methods appear each year. We compared methods utilized in our core facility: Smart-seq3, PlexWell, FLASH-seq, VASA-seq, SORT-seq, 10X, Evercode, and HIVE. We characterized the equipment requirements for each method. We evaluated the performances of these methods based on detected features, transcriptome diversity, mitochondrial RNA abundance and multiplets, among others and benchmarked them against bulk RNA sequencing. Here, we show that bulk transcriptome detects more unique transcripts than any single cell method. While most methods are comparable in many regards, FLASH-seq and VASA-seq yielded the best metrics, e.g., in number of features. If no equipment for automation is available or many cells are desired, then HIVE or 10X yield good results. In general, more recently developed methods perform better. This also leads to the conclusion that older methods should be phased out, and that the development of single cell RNAseq methods is still progressing considerably.

Zeb2 DNA-Binding Sites in Neuroprogenitor Cells Reveal Autoregulation and Affirm Neurodevelopmental Defects, Including in Mowat-Wilson Syndrome.

Functional perturbation and action mechanism studies have shown that the transcription factor Zeb2 controls cell fate decisions, differentiation, and/or maturation in multiple cell lineages in embryos and after birth. In cultured embryonic stem cells (ESCs), Zeb2's mRNA/protein upregulation is necessary for the exit from primed pluripotency and for entering general and neural differentiation. We edited mouse ESCs to produce Flag-V5 epitope-tagged Zeb2 protein from one endogenous allele. Using chromatin immunoprecipitation coupled with sequencing (ChIP-seq), we mapped 2432 DNA-binding sites for this tagged Zeb2 in ESC-derived neuroprogenitor cells (NPCs). A new, major binding site maps promoter-proximal to Zeb2 itself. The homozygous deletion of this site demonstrates that autoregulation of Zeb2 is necessary to elicit the appropriate Zeb2-dependent effects in ESC-to-NPC differentiation. We have also cross-referenced all the mapped Zeb2 binding sites with previously obtained transcriptome data from Zeb2 perturbations in ESC-derived NPCs, GABAergic interneurons from the ventral forebrain of mouse embryos, and stem/progenitor cells from the post-natal ventricular-subventricular zone (V-SVZ) in mouse forebrain, respectively. Despite the different characteristics of each of these neurogenic systems, we found interesting target gene overlaps. In addition, our study also contributes to explaining developmental disorders, including Mowat-Wilson syndrome caused by ZEB2 deficiency, and also other monogenic syndromes.

CTCF chromatin residence time controls three-dimensional genome organization, gene expression and DNA methylation in pluripotent cells.

The 11 zinc finger (ZF) protein CTCF regulates topologically associating domain formation and transcription through selective binding to thousands of genomic sites. Here, we replaced endogenous CTCF in mouse embryonic stem cells with green-fluorescent-protein-tagged wild-type or mutant proteins lacking individual ZFs to identify additional determinants of CTCF positioning and function. While ZF1 and ZF8-ZF11 are not essential for cell survival, ZF8 deletion strikingly increases the DNA binding off-rate of mutant CTCF, resulting in reduced CTCF chromatin residence time. Loss of ZF8 results in widespread weakening of topologically associating domains, aberrant gene expression and increased genome-wide DNA methylation. Thus, important chromatin-templated processes rely on accurate CTCF chromatin residence time, which we propose depends on local sequence and chromatin context as well as global CTCF protein concentration.

Molecular analysis of the erythroid phenotype of a patient with BCL11A haploinsufficiency.

The BCL11A gene encodes a transcriptional repressor with essential functions in multiple tissues during human development. Haploinsufficiency for BCL11A causes Dias-Logan syndrome (OMIM 617101), an intellectual developmental disorder with hereditary persistence of fetal hemoglobin (HPFH). Due to the severe phenotype, disease-causing variants in BCL11A occur de novo. We describe a patient with a de novo heterozygous variant, c.1453G>T, in the BCL11A gene, resulting in truncation of the BCL11A-XL protein (p.Glu485X). The truncated protein lacks the 3 C-terminal DNA-binding zinc fingers and the nuclear localization signal, rendering it inactive. The patient displayed high fetal hemoglobin (HbF) levels (12.1-18.7% of total hemoglobin), in contrast to the parents who had HbF levels of 0.3%. We used cultures of patient-derived erythroid progenitors to determine changes in gene expression and chromatin accessibility. In addition, we investigated DNA methylation of the promoters of the γ-globin genes HBG1 and HBG2. HUDEP1 and HUDEP2 cells were used as models for fetal and adult human erythropoiesis, respectively. Similar to HUDEP1 cells, the patient's cells displayed Assay for Transposase-Accessible Chromatin (ATAC) peaks at the HBG1/2 promoters and significant expression of HBG1/2 genes. In contrast, HBG1/2 promoter methylation and genome-wide gene expression profiling were consistent with normal adult erythropoiesis. We conclude that HPFH is the major erythroid phenotype of constitutive BCL11A haploinsufficiency. Given the essential functions of BCL11A in other hematopoietic lineages and the neuronal system, erythroid-specific targeting of the BCL11A gene has been proposed for reactivation of γ-globin expression in ß-hemoglobinopathy patients. Our data strongly support this approach.

Targeted chromatin conformation analysis identifies novel distal neural enhancers of ZEB2 in pluripotent stem cell differentiation.

The transcription factor zinc finger E-box binding protein 2 (ZEB2) controls embryonic and adult cell fate decisions and cellular maturation in many stem/progenitor cell types. Defects in these processes in specific cell types underlie several aspects of Mowat-Wilson syndrome (MOWS), which is caused by ZEB2 haplo-insufficiency. Human ZEB2, like mouse Zeb2, is located on chromosome 2 downstream of a ±3.5 Mb-long gene-desert, lacking any protein-coding gene. Using temporal targeted chromatin capture (T2C), we show major chromatin structural changes based on mapping in-cis proximities between the ZEB2 promoter and this gene desert during neural differentiation of human-induced pluripotent stem cells, including at early neuroprogenitor cell (NPC)/rosette state, where ZEB2 mRNA levels increase significantly. Combining T2C with histone-3 acetylation mapping, we identified three novel candidate enhancers about 500 kb upstream of the ZEB2 transcription start site. Functional luciferase-based assays in heterologous cells and NPCs reveal co-operation between these three enhancers. This study is the first to document in-cis Regulatory Elements located in ZEB2's gene desert. The results further show the usability of T2C for future studies of ZEB2 REs in differentiation and maturation of multiple cell types and the molecular characterization of newly identified MOWS patients that lack mutations in ZEB2 protein-coding exons.

Exome Sequencing Analysis Identifies Rare Variants in ATM and RPL8 That Are Associated With Shorter Telomere Length.

Telomeres are important for maintaining genomic stability. Telomere length has been associated with aging, disease, and mortality and is highly heritable (∼82%). In this study, we aimed to identify rare genetic variants associated with telomere length using whole-exome sequence data. We studied 1,303 participants of the Erasmus Rucphen Family (ERF) study, 1,259 of the Rotterdam Study (RS), and 674 of the British Heart Foundation Family Heart Study (BHF-FHS). We conducted two analyses, first we analyzed the family-based ERF study and used the RS and BHF-FHS for replication. Second, we combined the summary data of the three studies in a meta-analysis. Telomere length was measured by quantitative polymerase chain reaction in blood. We identified nine rare variants significantly associated with telomere length (p-value < 1.42 × 10-7, minor allele frequency of 0.2-0.5%) in the ERF study. Eight of these variants (in C11orf65, ACAT1, NPAT, ATM, KDELC2, and EXPH5) were located on chromosome 11q22.3 that contains ATM, a gene involved in telomere maintenance. Although we were unable to replicate the variants in the RS and BHF-FHS (p-value ≥ 0.21), segregation analysis showed that all variants segregate with shorter telomere length in a family. In the meta-analysis of all studies, a nominally significant association with LTL was observed with a rare variant in RPL8 (p-value = 1.48 × 10-6), which has previously been associated with age. Additionally, a novel rare variant in the known RTEL1 locus showed suggestive evidence for association (p-value = 1.18 × 10-4) with LTL. To conclude, we identified novel rare variants associated with telomere length. Larger samples size are needed to confirm these findings and to identify additional variants.

Characterization of the ferret TRB locus guided by V, D, J, and C gene expression analysis.

The domestic ferret, Mustela putorius furo, is an important mammalian animal model to study human respiratory infection. However, insufficient genomic annotation hampers detailed studies of ferret T cell responses. In this study, we analyzed the published T cell receptor beta (TRB) locus and performed high-throughput sequencing (HTS) of peripheral blood of four healthy adult ferrets to identify expressed V, D, J, and C genes. The HTS data is used as a guide to manually curate the expressed V, D, J, and C genes. The ferret locus appears to be most similar to that of the dog. Like other mammalian TRB loci, the ferret TRB locus contains a library of variable genes located upstream of two D-J-C gene clusters, followed by a (in the ferret non-functional) V gene with an inverted transcriptional orientation. All TRB genes (expressed or not) reported here have been approved by the IMGT/WHO-IUIS nomenclature committee.

A functional variant in the miR-142 promoter modulating its expression and conferring risk of Alzheimer disease.

Noncoding RNAs have been widely recognized as essential mediators of gene regulation. However, in contrast to protein-coding genes, much less is known about the influence of noncoding RNAs on human diseases. Here we examined the association of genetic variants located in primary microRNA sequences and long noncoding RNAs (lncRNAs) with Alzheimer disease (AD) by leveraging data from the largest genome-wide association meta-analysis of late-onset AD. Variants annotated to 5 miRNAs and 10 lncRNAs (in seven distinct loci) exceeded the Bonferroni-corrected significance threshold (p < 1.02 × 10-6 ). Among these, a leading variant (rs2526377:A>G) at the 17q22 locus annotated to two noncoding RNAs (MIR142 and BZRAP1-AS) was significantly associated with a reduced risk of AD and fulfilled predefined criteria for being a functional variant. Our functional genomic analyses revealed that rs2526377 affects the promoter activity and decreases the expression of miR-142. Moreover, differential expression analysis by RNA-Seq in human iPSC-derived neural progenitor cells and the hippocampus of miR-142 knockout mice demonstrated multiple target genes of miR-142 in the brain that are likely to be involved in the inflammatory and neurodegenerative manifestations of AD. These include TGFBR1 and PICALM, of which their derepression in the brain due to reduced expression levels of miR-142-3p may reduce the risk of AD.

Engram-specific transcriptome profiling of contextual memory consolidation.

Sparse populations of neurons in the dentate gyrus (DG) of the hippocampus are causally implicated in the encoding of contextual fear memories. However, engram-specific molecular mechanisms underlying memory consolidation remain largely unknown. Here we perform unbiased RNA sequencing of DG engram neurons 24 h after contextual fear conditioning to identify transcriptome changes specific to memory consolidation. DG engram neurons exhibit a highly distinct pattern of gene expression, in which CREB-dependent transcription features prominently (P = 6.2 × 10-13), including Atf3 (P = 2.4 × 10-41), Penk (P = 1.3 × 10-15), and Kcnq3 (P = 3.1 × 10-12). Moreover, we validate the functional relevance of the RNAseq findings by establishing the causal requirement of intact CREB function specifically within the DG engram during memory consolidation, and identify a novel group of CREB target genes involved in the encoding of long-term memory.

Whole-Genome Linkage Scan Combined With Exome Sequencing Identifies Novel Candidate Genes for Carotid Intima-Media Thickness.

Carotid intima-media thickness (cIMT) is an established heritable marker for subclinical atherosclerosis. In this study, we aim to identify rare variants with large effects driving differences in cIMT by performing genome-wide linkage analysis of individuals in the extremes of cIMT trait distribution (>90th percentile) in a large family-based study from a genetically isolated population in the Netherlands. Linked regions were subsequently explored by fine-mapping using exome sequencing. We observed significant evidence of linkage on chromosomes 2p16.3 [rs1017418, heterogeneity LOD (HLOD) = 3.35], 19q13.43 (rs3499, HLOD = 9.09), 20p13 (rs1434789, HLOD = 4.10), and 21q22.12 (rs2834949, HLOD = 3.59). Fine-mapping using exome sequencing data identified a non-coding variant (rs62165235) in PNPT1 gene under the linkage peak at chromosome 2 that is likely to have a regulatory function. The variant was associated with quantitative cIMT in the family-based study population (effect = 0.27, p-value = 0.013). Furthermore, we identified several genes under the linkage peak at chromosome 21 highly expressed in tissues relevant for atherosclerosis. To conclude, our linkage analysis identified four genomic regions significantly linked to cIMT. Further analyses are needed to demonstrate involvement of identified candidate genes in development of atherosclerosis.

Regulation of the cohesin-loading factor NIPBL: Role of the lncRNA NIPBL-AS1 and identification of a distal enhancer element.

Cohesin is crucial for genome stability, cell division, transcription and chromatin organization. Its functions critically depend on NIPBL, the cohesin-loader protein that is found to be mutated in >60% of the cases of Cornelia de Lange syndrome (CdLS). Other mutations are described in the cohesin subunits SMC1A, RAD21, SMC3 and the HDAC8 protein. In 25-30% of CdLS cases no mutation in the known CdLS genes is detected. Until now, functional elements in the noncoding genome were not characterized in the molecular etiology of CdLS and therefore are excluded from mutation screening, although the impact of such mutations has now been recognized for a wide range of diseases. We have identified different elements of the noncoding genome involved in regulation of the NIPBL gene. NIPBL-AS1 is a long non-coding RNA transcribed upstream and antisense to NIPBL. By knockdown and transcription blocking experiments, we could show that not the NIPBL-AS1 gene product, but its actual transcription is important to regulate NIPBL expression levels. This reveals a possibility to boost the transcriptional activity of the NIPBL gene by interfering with the NIPBL-AS1 lncRNA. Further, we have identified a novel distal enhancer regulating both NIPBL and NIPBL-AS1. Deletion of the enhancer using CRISPR genome editing in HEK293T cells reduces expression of NIPBL, NIPBL-AS1 as well as genes found to be dysregulated in CdLS.

Exome-Wide Meta-Analysis Identifies Rare 3'-UTR Variant in ERCC1/CD3EAP Associated with Symptoms of Sleep Apnea.

Obstructive sleep apnea (OSA) is a common sleep breathing disorder associated with an increased risk of cardiovascular and cerebrovascular diseases and mortality. Although OSA is fairly heritable (~40%), there have been only few studies looking into the genetics of OSA. In the present study, we aimed to identify genetic variants associated with symptoms of sleep apnea by performing a whole-exome sequence meta-analysis of symptoms of sleep apnea in 1,475 individuals of European descent. We identified 17 rare genetic variants with at least suggestive evidence of significance. Replication in an independent dataset confirmed the association of a rare genetic variant (rs2229918; minor allele frequency = 0.3%) with symptoms of sleep apnea (p-valuemeta = 6.98 × 10-9, ßmeta = 0.99). Rs2229918 overlaps with the 3' untranslated regions of ERCC1 and CD3EAP genes on chromosome 19q13. Both genes are expressed in tissues in the neck area, such as the tongue, muscles, cartilage and the trachea. Further, CD3EAP is localized in the nucleus and mitochondria and involved in the tumor necrosis factor-alpha/nuclear factor kappa B signaling pathway. Our results and biological functions of CD3EAP/ERCC1 genes suggest that the 19q13 locus is interesting for further OSA research.

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.

Unbiased Interrogation of 3D Genome Topology Using Chromosome Conformation Capture Coupled to High-Throughput Sequencing (4C-Seq).

The development and widespread implementation of chromosome conformation capture (3C) technology has allowed unprecedented new insight into how chromosomes are folded in three-dimensional (3D) space. 3C and its derivatives have contributed tremendously to the now widely accepted view that genome topology plays an important role in many major cellular processes, at a chromosome-wide scale, but certainly also at the level of individual genetic loci. A particularly popular application of 3C technology is to study transcriptional regulation, allowing researchers to draw maps of gene regulatory connections beyond the linear genome through addition of the third dimension. In this chapter, we provide a highly detailed protocol describing 3C coupled to high-throughput sequencing (referred to as 3C-Seq or more commonly 4C-Seq), allowing the unbiased interrogation of genome-wide chromatin interactions with specific genomic regions of interest. Interactions between spatially clustered DNA fragments are revealed by crosslinking the cells with formaldehyde, digesting the genome with a restriction endonuclease and performing a proximity ligation step to link interacting genomic fragments. Next, interactions with a selected DNA fragment are extracted from the 3C library through a second round of digestion and ligation followed by an inverse PCR. The generated products are immediately compatible with high-throughput sequencing, and amplicons from different PCR reactions can easily be multiplexed to dramatically increase throughput. Finally, we provide suggestions for data analysis and visualization.