Search

Background: The molecular basis of familial nonmedullary thyroid cancer (FNMTC) is still poorly understood, representing a limitation for molecular diagnosis and clinical management. In this study, we aimed to identify new susceptibility genes for FNMTC through whole-exome sequencing (WES) analysis of leukocyte DNA of patients from a highly informative FNMTC family. Methods: We selected six affected family members to conduct WES analysis. Bioinformatic analyses were undertaken to filter and select the genetic variants shared by the affected members, which were subsequently validated by Sanger sequencing. To select the most likely pathogenic variants, several studies were performed, including family segregation analysis, in silico impact characterization, and gene expression (messenger RNA and protein) depiction in databases. For the most promising variant identified, we performed in vitro studies to validate its pathogenicity. Results: Several potentially pathogenic variants were identified in different candidate genes. After filtering with appropriate criteria, the variant c.701C>T, p.Thr234Met in the SPRY4 gene was prioritized for in vitro functional characterization. This SPRY4 variant led to an increase in cell viability and colony formation, indicating that it confers a proliferative advantage and potentiates clonogenic capacity. Phosphokinase array and Western blot analyses suggested that the effects of the SPRY4 variant were mediated through the mitogen-activated protein kinase/extracellular signal-regulated kinase pathway, which was further supported by a higher responsiveness of thyroid cancer cells with the SPRY4 variant to a MEK inhibitor. Conclusions: WES analysis in one family identified SPRY4 as a likely novel candidate susceptibility gene for FNMTC, allowing a better understanding of the cellular and molecular mechanisms underlying thyroid cancer development.

CpG islands (CGIs) represent a widespread feature of vertebrate genomes, being associated with ~70% of all gene promoters. CGIs control transcription initiation by conferring nearby promoters with unique chromatin properties. In addition, there are thousands of distal or orphan CGIs (oCGIs) whose functional relevance is barely known. Here we show that oCGIs are an essential component of poised enhancers that augment their long-range regulatory activity and control the responsiveness of their target genes. Using a knock-in strategy in mouse embryonic stem cells, we introduced poised enhancers with or without oCGIs within topologically associating domains harboring genes with different types of promoters. Analysis of the resulting cell lines revealed that oCGIs act as tethering elements that promote the physical and functional communication between poised enhancers and distally located genes, particularly those with large CGI clusters in their promoters. Therefore, by acting as genetic determinants of gene–enhancer compatibility, CGIs can contribute to gene expression control under both physiological and potentially pathological conditions., T.P. is supported by a doctoral fellowship from the DAAD (Germany). V.S.-G. is supported by a doctoral fellowship from the University of Cantabria (Spain). Work in the Rada-Iglesias laboratory was supported by the EMBO Young Investigator Programme; CMMC intramural funding (Germany); the German Research Foundation (DFG) (Research Grant no. RA 2547/2-1); ‘Programa STAR-Santander Universidades, Campus Cantabria Internacional de la convocatoria CEI 2015 de Campus de Excelencia Internacional’ (Spain); the Spanish Ministry of Science, Innovation and Universities (Research Grant nos. PGC2018-095301-B-I00 and RED2018-102553-T REDEVNEURAL 3.0); and the European Research Council (ERC CoG ‘PoisedLogic’; grant no. 862022). The Landeira laboratory is funded by grants from the Spanish Ministry of Science and Innovation (grant nos. BFU2016-75233-P and PID2019-108108GB-I00) and the Andalusian Regional Government (grant no. PC-0246-2017).

Background Survival in patients with IDH1/2-mutant (mt) anaplastic astrocytomas is highly variable. We have used the prospective phase 3 CATNON trial to identify molecular factors related to outcome in IDH1/2mt anaplastic astrocytoma patients. Methods The CATNON trial randomized 751 adult patients with newly diagnosed 1p/19q non-codeleted anaplastic glioma to 59.4 Gy radiotherapy +/− concurrent and/or adjuvant temozolomide. The presence of necrosis and/or microvascular proliferation was scored at central pathology review. Infinium MethylationEPIC BeadChip arrays were used for genome-wide DNA methylation analysis and the determination of copy number variations (CNV). Two DNA methylation-based tumor classifiers were used for risk stratification. Next-generation sequencing (NGS) was performed using 1 of the 2 glioma-tailored NGS panels. The primary endpoint was overall survival measured from the date of randomization. Results Full analysis (genome-wide DNA methylation and NGS) was successfully performed on 654 tumors. Of these, 432 tumors were IDH1/2mt anaplastic astrocytomas. Both epigenetic classifiers identified poor prognosis patients that partially overlapped. A predictive prognostic Cox proportional hazard model identified that independent prognostic factors for IDH1/2mt anaplastic astrocytoma patients included; age, mini-mental state examination score, treatment with concurrent and/or adjuvant temozolomide, the epigenetic classifiers, PDGFRA amplification, CDKN2A/B homozygous deletion, PI3K mutations, and total CNV load. Independent recursive partitioning analysis highlights the importance of these factors for patient prognostication. Conclusion Both clinical and molecular factors identify IDH1/2mt anaplastic astrocytoma patients with worse outcome. These results will further refine the current WHO criteria for glioma classification.

Following implantation, the naive pluripotent epiblast of the mouse blastocyst generates a rosette, undergoes lumenogenesis and forms the primed pluripotent egg cylinder, which is able to generate the embryonic tissues. How pluripotency progression and morphogenesis are linked and whether intermediate pluripotent states exist remain controversial. We identify here a rosette pluripotent state defined by the co-expression of naive factors with the transcription factor OTX2. Downregulation of blastocyst WNT signals drives the transition into rosette pluripotency by inducing OTX2. The rosette then activates MEK signals that induce lumenogenesis and drive progression to primed pluripotency. Consequently, combined WNT and MEK inhibition supports rosette-like stem cells, a self-renewing naive-primed intermediate. Rosette-like stem cells erase constitutive heterochromatin marks and display a primed chromatin landscape, with bivalently marked primed pluripotency genes. Nonetheless, WNT induces reversion to naive pluripotency. The rosette is therefore a reversible pluripotent intermediate whereby control over both pluripotency progression and morphogenesis pivots from WNT to MEK signals. Neagu, van Genderen and Escudero et al. show that simultaneous inhibition of WNT and MEK signalling maintains a naive-primed intermediate pluripotency state in vitro, which displays features of the mouse embryonic rosette.

Background Neonatal jaundice is a common finding in newborns in Asia, including Indonesia. In some cases, the serum total bilirubin levels exceeds the 95th percentile for hours of life (neonatal hyperbilirubinemia). Severe neonatal hyperbilirubinemia (NH) could lead to kernicterus and neonatal death. Glucose-6-Phosphage Dehydrogenase (G6PD) genetic variations and deficiency have been reported in several studies to be associated with NH. This study aimed to analyze the G6PD genetic variations and its activity in neonates with and without hyperbilirubinemia in the Deutromalay Indonesian population. Methods Deoxyribose Nucleic Acid (DNA) was isolated from peripheral blood of 116 and 115 healthy term neonates with and without hyperbilirubinemia. All infants underwent the following laboratory examinations: routine hematologic evaluation, Coombs test, G6PD activity measurement using the Randox kit method, and serum total bilirubin level. All exons of the G6PD gene were targeted for deep sequencing using MiSeq (Illumina). An association study of G6PD polymorphisms with NH was performed using PLINK. Results The prevalence of G6PD deficiency in neonates with and without hyperbilirubinemia in Indonesian Deutromalay population were 1.72% (95% Confidence Interval (CI): 0.6–4.1%) and 1.74% (95% CI: 0.7–4.1%), respectively. The most common G6PD polymorphisms, i.e. rs1050757/c.* + 357A > G, rs2230037/c.1311C > T, and rs2071429/c.1365-13 T/IVS11, were identified. However, none of those polymorphisms and their haplotype were associated with NH (p > 0.05, Odds Ratio (OR) ~1.00). The prevalence of G6PD mutations in neonates with and without hyperbilirubinemia were 6.8% (95% CI: 2.3–11.5%) and 6.9% (95% CI: 2.3–11.6%), respectively. The most frequently identified G6PD mutation was the Viangchan variant (p.V291 M), which was followed by the Canton (p.R459L) and Vanua Lava (p.L128P) variants. Two novel mutations were identified both in case (p.V369A, p.I167F) and control (p.L474=, p.I36T) groups. Conclusion The prevalence of G6PD deficiency is low in neonates with or without hyperbilirubinemia in Deutromalay Indonesian population. The majority of G6PD mutations identified among Indonesian Deutromalay population in this study are Viangchan, Canton and Vanua Lava variants.

Evasion from drug-induced apoptosis is a crucial mechanism of cancer treatment resistance. The pro-apoptotic protein NOXA marks an aggressive pancreatic ductal adenocarcinoma (PDAC) subtype. To identify drugs that unleash the death-inducing potential of NOXA, we performed an unbiased drug screening experiment. In NOXA-deficient isogenic cellular models we identified an inhibitor of the transcription factor heterodimer CBFβ/RUNX1. By genetic gain and loss of function experiments we validated that the mode of action depends on RUNX1 and NOXA. Of note, RUNX1 expression is significantly higher in PDACs compared to normal pancreas. We show that pharmacological RUNX1 inhibition significantly blocks tumor growth in vivo and in primary patient-derived PDAC organoids. Through genome wide analysis, we detected that RUNX1-loss reshapes the epigenetic landscape, which gains H3K27ac enrichment at the NOXA promoter. Our study demonstrates a previously unknown mechanism of NOXA-dependent cell death, which can be triggered pharmaceutically. Therefore, our data show a novel way to target a therapy resistant PDAC, an unmet clinical need.SignificanceRecent evidence demonstrated the existence of molecular subtypes in pancreatic ductal adenocarcinoma (PDAC), which resist all current therapies. The paucity of therapeutic options, including a complete lack of targeted therapies, underscore the urgent and unmet medical need for the identification of targets and novel treatment strategies for PDAC. Our study unravels a function of the transcription factor RUNX1 in apoptosis regulation in PDAC. We show that pharmacological RUNX1 inhibition in PDAC is feasible and leads to NOXA-dependent apoptosis. The development of targeted therapies that influence the transcriptional landscape of PDAC might have great benefits for patients who are resistant to conventional therapies. RUNX1 Inhibition as a new therapeutic intervention offers an attractive strategy for future therapies.

Background Alveolar capillary dysplasia with or without misalignment of the pulmonary veins (ACD/MPV) is a lethal congenital lung disorder associated with a variety of heterozygous genomic alterations in the FOXF1 gene or its 60 kb enhancer. Cases without a genomic alteration in the FOXF1 locus have been described as well. The mechanisms responsible for FOXF1 haploinsufficiency and the cause of ACD/MPV in patients without a genomic FOXF1 variant are poorly understood, complicating the search for potential therapeutic targets for ACD/MPV. To investigate the contribution of aberrant DNA methylation, genome wide methylation patterns of ACD/MPV lung tissues were compared with methylation patterns of control lung tissues using the recently developed technique Methylated DNA sequencing (MeD-seq). Results Eight ACD/MPV lung tissue samples and three control samples were sequenced and their mutual comparison resulted in identification of 319 differentially methylated regions (DMRs) genome wide, involving 115 protein coding genes. The potentially upregulated genes were significantly enriched in developmental signalling pathways, whereas potentially downregulated genes were mainly enriched in O-linked glycosylation. In patients with a large maternal deletion encompassing the 60 kb FOXF1 enhancer, DNA methylation patterns in this FOXF1 enhancer were not significantly different compared to controls. However, two hypermethylated regions were detected in the 60 kb FOXF1 enhancer of patients harbouring a FOXF1 point mutation. Lastly, a large hypermethylated region overlapping the first FOXF1 exon was found in one of the ACD/MPV patients without a known pathogenic FOXF1 variation. Conclusion This is the first study providing genome wide methylation data on lung tissue of ACD/MPV patients. DNA methylation analyses in the FOXF1 locus excludes maternal imprinting of the 60 kb FOXF1 enhancer. Hypermethylation at the 60 kb FOXF1 enhancer might contribute to FOXF1 haploinsufficiency caused by heterozygous mutations in the FOXF1 coding region. Interestingly, DNA methylation analyses of patients without a genomic FOXF1 variant suggest that abnormal hypermethylation of exon 1 might play a role in some ACD/MPV in patients.

An innovative approach to eliminate HIV-1-infected cells emerging out of latency, the major hurdle to HIV-1 cure, is to pharmacologically reactivate viral expression and concomitantly trigger intracellular pro-apoptotic pathways in order to selectively induce cell death (ICD) of infected cells, without reliance on the extracellular immune system. In this work, we demonstrate the effect of DDX3 inhibitors on selectively inducing cell death in latent HIV-1-infected cell lines, primary CD4+ T cells and in CD4+ T cells from cART-suppressed people living with HIV-1 (PLWHIV). We used single-cell FISH-Flow technology to characterise the contribution of viral RNA to inducing cell death. The pharmacological targeting of DDX3 induced HIV-1 RNA expression, resulting in phosphorylation of IRF3 and upregulation of IFNβ. DDX3 inhibition also resulted in the downregulation of BIRC5, critical to cell survival during HIV-1 infection, and selectively induced apoptosis in viral RNA-expressing CD4+ T cells but not bystander cells. DDX3 inhibitor treatment of CD4+ T cells from PLWHIV resulted in an approximately 50% reduction of the inducible latent HIV-1 reservoir by quantitation of HIV-1 RNA, by FISH-Flow, RT-qPCR and TILDA. This study provides proof of concept for pharmacological reversal of latency coupled to induction of apoptosis towards the elimination of the inducible reservoir., DEAD-box polypeptide 3 (DDX3) is a host protein belonging to the family of ATP-dependent RNA helicases. Here, the authors demonstrate that DDX3 inhibitors reverse HIV-1 latency and selectively induce cell death in HIV-1-infected cell lines, primary CD4+ T cells and in CD4+ T cells from cART-suppressed people living with HIV-1.

Aims Hepatic capillaries are lined with specialized liver sinusoidal endothelial cells (LSECs) which support macromolecule passage to hepatocytes and prevent fibrosis by keeping hepatic stellate cells (HSCs) quiescent. LSEC specialization is co-determined by transcription factors. The zinc-finger E-box-binding homeobox (Zeb)2 transcription factor is enriched in LSECs. Here, we aimed to elucidate the endothelium-specific role of Zeb2 during maintenance of the liver and in liver fibrosis. Methods and results To study the role of Zeb2 in liver endothelium we generated EC-specific Zeb2 knock-out (ECKO) mice. Sequencing of liver EC RNA revealed that deficiency of Zeb2 results in prominent expression changes in angiogenesis-related genes. Accordingly, the vascular area was expanded and the presence of pillars inside ECKO liver vessels indicated that this was likely due to increased intussusceptive angiogenesis. LSEC marker expression was not profoundly affected and fenestrations were preserved upon Zeb2 deficiency. However, an increase in continuous EC markers suggested that Zeb2-deficient LSECs are more prone to dedifferentiation, a process called ‘capillarization’. Changes in the endothelial expression of ligands that may be involved in HSC quiescence together with significant changes in the expression profile of HSCs showed that Zeb2 regulates LSEC–HSC communication and HSC activation. Accordingly, upon exposure to the hepatotoxin carbon tetrachloride (CCl4), livers of ECKO mice showed increased capillarization, HSC activation, and fibrosis compared to livers from wild-type littermates. The vascular maintenance and anti-fibrotic role of endothelial Zeb2 was confirmed in mice with EC-specific overexpression of Zeb2, as the latter resulted in reduced vascularity and attenuated CCl4-induced liver fibrosis. Conclusion Endothelial Zeb2 preserves liver angioarchitecture and protects against liver fibrosis. Zeb2 and Zeb2-dependent genes in liver ECs may be exploited to design novel therapeutic strategies to attenuate hepatic fibrosis.

Group 2 innate lymphoid cells (ILC2s) orchestrate protective type 2 immunity and have been implicated in various immune disorders. In the mouse, circulatory inflammatory ILC2s (iILC2s) were identified as a major source of type 2 cytokines. The human equivalent of the iILC2 subset remains unknown. Here, we identify a human inflammatory ILC2 population that resides in inflamed mucosal tissue and is specifically marked by surface CD45RO expression. CD45RO+ ILC2s are derived from resting CD45RA+ ILC2s upon activation by epithelial alarmins such as IL-33 and TSLP, which is tightly linked to STAT5 activation and up-regulation of the IRF4/BATF transcription factors. Transcriptome analysis reveals marked similarities between human CD45RO+ ILC2s and mouse iILC2s. Frequencies of CD45RO+ inflammatory ILC2 are increased in inflamed mucosal tissue and in the circulation of patients with chronic rhinosinusitis or asthma, correlating with disease severity and resistance to corticosteroid therapy. CD45RA-to-CD45RO ILC2 conversion is suppressed by corticosteroids via induction of differentiation toward an immunomodulatory ILC2 phenotype characterized by low type 2 cytokine and high amphiregulin expression. Once converted, however, CD45RO+ ILC2s are resistant to corticosteroids, which is associated with metabolic reprogramming resulting in the activation of detoxification pathways. Our combined data identify CD45RO+ inflammatory ILC2s as a human analog of mouse iILC2s linked to severe type 2 inflammatory disease and therapy resistance.

INTRODUCTION: The majority of desmoid-type fibromatosis (DTF) tumors harbor a β-catenin mutation, affecting specific codons in CTNNB1 exon 3. S45F tumors are reported to have a higher chance of recurrence after surgery and more resistance to systemic treatments compared to wild-type (WT) and T41A tumors. The aim of this pilot study was to examine the genome-wide DNA methylation profiles of S45F and T41A mutated DTF, to explain the observed differences in clinical behavior between these DTF subtypes.MATERIAL AND METHODS: Genome-wide analysis of DNA methylation was performed using MeD-seq on formalin-fixed, paraffin-embedded primary DTF samples harboring a S45F (n = 14) or a T41A (n = 15) mutation. Differentially methylated regions (DMRs) between S45F and T41A DTF were identified and used for a supervised hierarchical cluster analysis. DMRs with a fold-change ≥1.5 were considered to be differentially methylated and differences between S45F and T41A tumors were quantitatively assessed. The effect of DMRs on the expression of associated genes was assessed using an independent mRNA expression dataset. Protein-protein interactions between WT β-catenin and mutant variants and DNA methyltransferase 1 (DNMT1) were examined by immunoprecipitation experiments.RESULTS: MeD-seq analyses indicated 354 regions that displayed differential methylation. Cluster analysis yielded no distinct clusters based on mutation, sex, tumor site or tumor size. A supervised clustering based on DMRs between small (≤34 mm) and large (>87 mm) DTF distinguished the two groups. Only ten DMRs displayed a fold change of ≥1.5 and six of them were found associated with the following genes: NLRP4, FOXK2, PERM1, CCDC6, NOC4L, and DUX4L6. The effects of DMRs on gene expression yielded a significant difference (p < 0.05) in the expression between S45F and T41A for CCDC6 and FOXK2 but not for all Affymetrix probe-sets used to detect these genes. Immunoprecipitations did not reveal an association of WT β-catenin or mutant variants with DNMT1.CONCLUSION: This study demonstrated that S45F and T41A DTF tumors did not exhibit gross differences in DNA methylation patterns. This implies that distinct DNA methylation profiles are not the sole determinant for the divergent clinical behavior of these different DTF mutant subtypes.

Transcriptomics can be combined with TRA and TRB clonotype analysis at the single cell level. The aim of this study was to validate this approach on the ICELL8 Single-Cell system and to evaluate its usefulness to analyse clinical paucicellular samples. For this purpose, we carefully selected T cell lines with defined TRA/TRB clonotypes as well as clinical samples enriched for CD3+ T cells that possess a complex TCR repertoire. Low cell numbers of the different samples were dispensed in a chip on the ICELL8 Single-Cell System. Two sequencing libraries were generated from each single cell cDNA preparation, one for the TRA/TRB repertoire and one for the 5' ends of transcripts, and subsequently sequenced. Transcriptome analysis revealed that the cell lines on average express 2,268 unique genes/cell and T cells of clinical samples 770 unique genes/cell. The expected combined TRA/TRB clonotype was determined for on average 71% of the cells of the cell lines. In the clinical samples the TRA/TRB repertoire was more complex than those of the cell lines. Furthermore, the TRB clonotype distribution of the clinical samples was positively correlated to frequencies of TCRVβ families in CD3+ T cells obtained by a flow cytometry-based approach (Spearman's Rho correlation coefficient 0.81, P = 6.49 * 10-7). Combined analyses showed that transcriptome-based cell type-specific clusters in clinical samples corresponded to clinical features such as CMV status. In conclusion, we showed that the ICELL8 Single-Cell System enabled combined interrogation of both TRA/TRB repertoire and transcriptome of paucicellular clinical samples. This opens the way to study the response of single T cells within heterogeneous samples for both their transcriptome and TRA/TRB clonotypes in disease or upon treatment.

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. ispartof: HUMAN MOLECULAR GENETICS vol:29 issue:15 pages:2535-2550 ispartof: location:England status: published

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.

Autosomal‐dominant polycystic kidney disease (ADPKD) is the most common inherited kidney disease, leading to kidney failure in most patients. In approximately 85% of cases, the disease is caused by mutations in PKD1. How dysregulation of PKD1 leads to cyst formation on a molecular level is unknown. Induced pluripotent stem cells (iPSCs) are a powerful tool for in vitro modeling of genetic disorders. Here, we established ADPKD patient‐specific iPSCs to study the function of PKD1 in kidney development and cyst formation in vitro. Somatic mutations are proposed to be the initiating event of cyst formation, and therefore, iPSCs were derived from cystic renal epithelial cells rather than fibroblasts. Mutation analysis of the ADPKD iPSCs revealed germline mutations in PKD1 but no additional somatic mutations in PKD1/PKD2. Although several somatic mutations in other genes implicated in ADPKD were identified in cystic renal epithelial cells, only few of these mutations were present in iPSCs, indicating a heterogeneous mutational landscape, and possibly in vitro cell selection before and during the reprogramming process. Whole‐genome DNA methylation analysis indicated that iPSCs derived from renal epithelial cells maintain a kidney‐specific DNA methylation memory. In addition, comparison of PKD1+/− and control iPSCs revealed differences in DNA methylation associated with the disease history. In conclusion, we generated and characterized iPSCs derived from cystic and healthy control renal epithelial cells, which can be used for in vitro modeling of kidney development in general and cystogenesis in particular., Autosomal dominant polycystic kidney disease (ADPKD) is the most common inherited kidney disease. Here, we generated induced pluripotent stem cells (iPSCs) of ADPKD patients to study kidney development and cyst formation in vitro. These iPSCs revealed germline and autosomal mutations implicated in ADPKD and displayed an epigenetic memory of kidney epithelial cells, providing powerful models to study ADPKD in vitro.

The transcription factor Zeb2 controls fate specification and subsequent differentiation and maturation of multiple cell types in various embryonic tissues. It binds many protein partners, including activated Smad proteins and the NuRD co-repressor complex. How Zeb2 subdomains support cell differentiation in various contexts has remained elusive. Here, we have studied the role of Zeb2 and its domains in neurogenesis and neural differentiation in the young postnatal ventricular-subventricular zone (V-SVZ), where neural stem cells generate olfactory bulb-destined interneurons. Conditional Zeb2 knockouts and separate acute loss- and gain-of-function approaches indicated that Zeb2 is essential to control apoptosis and neuronal differentiation of V-SVZ progenitors before and after birth, and identified Sox6 as Zeb2-dependent and potential downstream target gene. Zeb2 genetic inactivation impaired the differentiation potential of the V-SVZ niche in a cell-autonomous fashion. We also provide evidence that its normal function in the V-SVZ involves non-autonomous mechanisms as well. Additionally, we could demonstrate distinct roles for Zeb2 protein-binding domains, suggesting that Zeb2 partners co-determine neuronal output from the mouse V-SVZ in both quantitative and qualitative manners in early postnatal life. ispartof: Development vol:147 issue:10 ispartof: location:England status: published

Developmental and/or epileptic encephalopathies (DEEs) are a group of devastating genetic disorders, resulting in early-onset, therapy-resistant seizures and developmental delay. Here we report on 22 individuals from 15 families presenting with a severe form of intractable epilepsy, severe developmental delay, progressive microcephaly, visual disturbance and similar minor dysmorphisms. Whole exome sequencing identified a recurrent, homozygous variant (chr2:64083454A > G) in the essential UDP-glucose pyrophosphorylase (UGP2) gene in all probands. This rare variant results in a tolerable Met12Val missense change of the longer UGP2 protein isoform but causes a disruption of the start codon of the shorter isoform, which is predominant in brain. We show that the absence of the shorter isoform leads to a reduction of functional UGP2 enzyme in neural stem cells, leading to altered glycogen metabolism, upregulated unfolded protein response and premature neuronal differentiation, as modeled during pluripotent stem cell differentiation in vitro. In contrast, the complete lack of all UGP2 isoforms leads to differentiation defects in multiple lineages in human cells. Reduced expression of Ugp2a/Ugp2b in vivo in zebrafish mimics visual disturbance and mutant animals show a behavioral phenotype. Our study identifies a recurrent start codon mutation in UGP2 as a cause of a novel autosomal recessive DEE syndrome. Importantly, it also shows that isoform-specific start-loss mutations causing expression loss of a tissue-relevant isoform of an essential protein can cause a genetic disease, even when an organism-wide protein absence is incompatible with life. We provide additional examples where a similar disease mechanism applies. Electronic supplementary material The online version of this article (10.1007/s00401-019-02109-6) contains supplementary material, which is available to authorized users.

Dendritic cells (DCs) are key immune modulators and are able to mount immune responses or tolerance. DC differentiation and activation imply a plethora of molecular and cellular responses, including transcriptional changes. PU.1 is a highly expressed transcription factor in DCs and coordinates relevant aspects of DC biology. Due to their role as immune regulators, DCs pose as a promising immunotherapy tool. However, some of their functional features, such as survival, activation, or migration, are compromised due to the limitations to simulate in vitro the physiologic DC differentiation process. A better knowledge of transcriptional programs would allow the identification of potential targets for manipulation with the aim of obtaining “qualified” DCs for immunotherapy purposes. Most of the current knowledge regarding DC biology derives from studies using mouse models, which not always find a parallel in human. In the present study, we dissect the PU.1 transcriptional regulome and interactome in mouse and human DCs, in the steady state or LPS activated. The PU.1 transcriptional regulome was identified by performing PU.1 chromatin immunoprecipitation followed by high-throughput sequencing and pairing these data with RNAsequencing data. The PU.1 interactome was identified by performing PU.1 immunoprecipitation followed by mass spectrometry analysis. Our results portray PU.1 as a pivotal factor that plays an important role in the regulation of genes required for proper DC activation and function, and assures the repression of nonlineage genes. The interspecies differences between human and mouse DCs are surprisingly substantial, highlighting the need to study the biology of human DCs.

Allergic asthma is mediated by Th2 responses to inhaled allergens. Although previous experiments indicated that Notch signaling activates expression of the key Th2 transcription factor Gata3, it remains controversial how Notch promotes allergic airway inflammation. Here we show that T cell-specific Notch deficiency in mice prevented house dust mite-driven eosinophilic airway inflammation and significantly reduced Th2 cytokine production, serum IgE levels, and airway hyperreactivity. However, transgenic Gata3 overexpression in Notch-deficient T cells only partially rescued this phenotype. We found that Notch signaling was not required for T cell proliferation or Th2 polarization. Instead, Notch-deficient in vitro-polarized Th2 cells showed reduced accumulation in the lungs upon in vivo transfer and allergen challenge, as Notch-deficient Th2 cells were retained in the lung-draining lymph nodes. Transcriptome analyses and sequential adoptive transfer experiments revealed that while Notch-deficient lymph node Th2 cells established competence for lung migration, they failed to upregulate sphingosine-1-phosphate receptor 1 (S1PR1) and its critical upstream transcriptional activator Krüppel-like factor 2 (KLF2). As this KLF2/S1PR1 axis represents the essential cell-intrinsic regulator of T cell lymph node egress, we conclude that the druggable Notch signaling pathway licenses the Th2 response in allergic airway inflammation via promoting lymph node egress.

Background Approximately one-third of the patients with well-differentiated liposarcoma (WDLPS) will develop a local recurrence after resection of the primary tumor. Not much is known about the molecular relationship between the primary tumor and the recurrent tumor, which is important to reveal potential drivers of recurrence. Here we investigated the biology of recurrent WDLPS by comparing paired primary and recurrent WDLPS using microRNA profiling and genome-wide DNA methylation analyses. Methods Paired primary and recurrent WDLPS formalin-fixed and paraffin-embedded tumor samples were collected (N = 27 pairs). MicroRNA expression profiles were determined using TaqMan® Low Density Array (TLDA) cards, capable of detecting 754 microRNAs. Genome-wide DNA methylation and differentially methylated regions were assessed by methylated DNA sequencing (MeD-seq). Results A supervised cluster analysis based on differentially expressed microRNAs between paired primary and recurrent WDLPS did not reveal a clear cluster pattern separating the primary from the recurrent tumors. The clustering was also not based on tumor localization, time to recurrence, age or the status of the resection margins. Changes in DNA methylation between primary and recurrent tumors were extremely variable, and no consistent DNA methylation changes were found. As a result, a supervised clustering analysis based on differentially methylated regions between primary and recurrent tumors did not show a distinct cluster pattern based on any of the features. Subgroup analysis for tumors localized in the extremity or the retroperitoneum also did not yield a clear distinction between primary and recurrent WDLPS samples. Conclusions Primary and recurrent WDLPS could not be distinguished based on microRNA expression and DNA methylation profiles and no common microRNA and DNA methylation alterations for recurrence could be identified. Legal entity responsible for the study The authors. Funding Has not received any funding. Disclosure All authors have declared no conflicts of interest.

Background The genetic bases of PD in sub-Saharan African (SSA) populations remain poorly characterized, and analysis of SSA families with PD might lead to the discovery of novel disease-related genes. Objectives To investigate the clinical features and identify the disease-causing gene in a black South African family with 3 members affected by juvenile-onset parkinsonism and intellectual disability. Methods Clinical evaluation, neuroimaging studies, whole-exome sequencing, homozygosity mapping, two-point linkage analysis, and Sanger sequencing of candidate variants. Result A homozygous 28-nucleotide frameshift deletion in the PTRHD1 coding region was identified in the 3 affected family members and linked to the disease with genome-wide significant evidence. PTRHD1 was recently nominated as the disease-causing gene in two Iranian families, each containing 2 siblings with similar phenotypes and homozygous missense mutations. Conclusion Together with the previous reports, we provide conclusive evidence that loss-of-function mutations in PTRHD1 cause autosomal-recessive juvenile parkinsonism and intellectual disability. © 2018 International Parkinson and Movement Disorder Society.

Melanoma is one of the most aggressive cancers, and patients with distant metastases have dire outcomes. We observed previously that melanoma progression is driven by a high migratory potential of melanoma cells, which survive and proliferate under harsh environmental conditions. In this study, we report that CREPT (cell-cycle related and expression-elevated protein in tumor), an oncoprotein highly expressed in other cancers, is overexpressed in melanoma cells but not melanocytes. Overexpression of CREPT stimulates cell proliferation, migration, and invasion in several melanoma cell lines. Further, we show that CREPT enhances melanoma progression through upregulating and activating Ras homolog family member A (RhoA)-induced actin organization and focal adhesion assembly. Our study reveals a novel role of CREPT in promoting melanoma progression. Targeting CREPT may be a promising strategy for melanoma treatment.

Summary Background Most patients with Parkinson's disease, Parkinson's disease dementia, and dementia with Lewy bodies do not carry mutations in known disease-causing genes. The aim of this study was to identify a novel gene implicated in the development of these disorders. Methods Our study was done in three stages. First, we did genome-wide linkage analysis of an Italian family with dominantly inherited Parkinson's disease to identify the disease locus. Second, we sequenced the candidate gene in an international multicentre series of unrelated probands who were diagnosed either clinically or pathologically with Parkinson's disease, Parkinson's disease dementia, or dementia with Lewy bodies. As a control, we used gene sequencing data from individuals with abdominal aortic aneurysms (who were not examined neurologically). Third, we enrolled an independent series of patients diagnosed clinically with Parkinson's disease and controls with no signs or family history of Parkinson's disease, Parkinson's disease dementia, or dementia with Lewy bodies from centres in Portugal, Sardinia, and Taiwan, and screened them for specific variants. We also did mRNA and brain pathology studies in three patients from the international multicentre series carrying disease-associated variants, and we did functional protein studies in in-vitro models, including neurons from induced pluripotent stem-like cells. Findings Molecular studies were done between Jan 1, 2008, and Dec 31, 2017. In the initial kindred of ten affected Italian individuals (mean age of disease onset 59·8 years [SD 8·7]), we detected significant linkage of Parkinson's disease to chromosome 14 and nominated LRP10 as the disease-causing gene. Among the international series of 660 probands, we identified eight individuals (four with Parkinson's disease, two with Parkinson's disease dementia, and two with dementia with Lewy bodies) who carried different, rare, potentially pathogenic LRP10 variants; one carrier was found among 645 controls with abdominal aortic aneurysms. In the independent series, two of these eight variants were detected in three additional Parkinson's disease probands (two from Sardinia and one from Taiwan) but in none of the controls. Of the 11 probands from the international and independent cohorts with LRP10 variants, ten had a positive family history of disease and DNA was available from ten affected relatives (in seven of these families). The LRP10 variants were present in nine of these ten relatives, providing independent—albeit limited—evidence of co-segregation with disease. Post-mortem studies in three patients carrying distinct LRP10 variants showed severe Lewy body pathology. Of nine variants identified in total (one in the initial family and eight in stage 2), three severely affected LRP10 expression and mRNA stability (1424+5delG, 1424+5G→A, and Ala212Serfs*17, shown by cDNA analysis), four affected protein stability (Tyr307Asn, Gly603Arg, Arg235Cys, and Pro699Ser, shown by cycloheximide-chase experiments), and two affected protein localisation (Asn517del and Arg533Leu; shown by immunocytochemistry), pointing to loss of LRP10 function as a common pathogenic mechanism. Interpretation Our findings implicate LRP10 gene defects in the development of inherited forms of α-synucleinopathies. Future elucidation of the function of the LRP10 protein and pathways could offer novel insights into mechanisms, biomarkers, and therapeutic targets. Funding Stichting ParkinsonFonds, Dorpmans-Wigmans Stichting, Erasmus Medical Center, ZonMw—Memorabel programme, EU Joint Programme Neurodegenerative Disease Research (JPND), Parkinson's UK, Avtal om Lakarutbildning och Forskning (ALF) and Parkinsonfonden (Sweden), Lijf and Leven foundation, and cross-border grant of Alzheimer Netherlands–Ligue Europeene Contre la Maladie d'Alzheimer (LECMA).

Background Group 2 innate lymphoid cells (ILC2s) are major producers of the cytokines driving allergic asthma, and increased ILC2 numbers have been detected in blood and sputum of asthmatic patients. Asthma susceptibility has a strong genetic component, but the underlying mechanisms and whether asthma genetics affect ILC2 biology remain unclear. Objective We sought to study the ILC2 transcriptome and epigenome during airway inflammation (AI) to couple these to genes and genetic variants associated with asthma pathogenesis. Methods Mice harboring a reporter for the key ILC2 transcription factor GATA-3 were subjected to IL-33–driven AI, and ILC2s were isolated from bronchoalveolar lavage fluid and mediastinal lymph nodes. Human ILC2s were purified from peripheral blood and activated in vitro . We used RNA sequencing, genome-wide identification of histone-3 lysine-4 dimethylation–marked chromatin, and computational approaches to study the ILC2 transcriptome and epigenome. Results Activated ILC2s in mice displayed a tissue-specific gene expression signature that emerged from remarkably similar epigenomes. We identified superenhancers implicated in controlling ILC2 identity and asthma-associated genes. More than 300 asthma-associated genetic polymorphisms identified in genome-wide association studies localized to H3K4Me2 + gene regulatory elements in ILC2s. A refined set of candidate causal asthma-associated variants was uniquely enriched in ILC2, but not T H 2 cell, regulatory regions. Conclusions ILC2s in AI use a flexible epigenome that couples adaptation to new microenvironments with functional plasticity. Importantly, we reveal strong correlations between gene regulatory mechanisms in ILC2s and the genetic basis of asthma, supporting a pathogenic role for ILC2s in patients with allergic asthma.

Prenatal diagnostics has been impacted by technological changes in the past decade, which have affected the diagnostic yield. The aim of this study was to evaluate the impact of SNP array and noninvasive prenatal testing (NIPT) on the diagnostic yield and the number of invasive tests in our center. The frequency of pathogenic fetal unbalanced chromosome aberrations was studied in 10,005 cases referred for prenatal testing in 2009-2015. Chromosomal SNP microarray analysis replaced karyotyping in all invasively tested pregnancies and since 2014 a choice between NIPT and diagnostic testing with microarray was offered to women with an increased risk for common aneuploidy. The introduction of microarray led to an additional yield of submicroscopic pathogenic chromosome aberrations: 3.6% in fetuses with ultrasound anomalies and 1.9% in fetuses without ultrasound anomalies. The introduction of NIPT led to a decrease of invasive tests and of diagnostic yield. Moreover, a diagnostic delay in about 1:350 cases was observed. Since 20%-33% of pathogenic fetal chromosome aberrations are different from the common aneuploidies and triploidy, whole-genome analysis should be offered after invasive sampling. Because NIPT (as a second screening) has led to a decreased diagnostic yield, it should be accompanied by an appropriate pretest counseling.

Background Despite high heritability, little success was achieved in mapping genetic determinants of depression-related traits by means of genome-wide association studies. Methods To identify genes associated with depressive symptomology, we performed a gene-based association analysis of nonsynonymous variation captured using exome-sequencing and exome-chip genotyping in a genetically isolated population from the Netherlands ( n = 1999). Finally, we reproduced our significant findings in an independent population-based cohort ( n = 1604). Results We detected significant association of depressive symptoms with a gene NKPD1 ( p = 3.7 × 10 −08 ). Nonsynonymous variants in the gene explained 0.9% of sex- and age-adjusted variance of depressive symptoms in the discovery study, which is translated into 3.8% of the total estimated heritability ( h 2 = 0.24). Significant association of depressive symptoms with NKPD1 was also observed ( n = 1604; p = 1.5 × 10 −03 ) in the independent replication sample despite little overlap with the discovery cohort in the set of nonsynonymous genetic variants observed in the NKPD1 gene. Meta-analysis of the discovery and replication studies improved the association signal ( p = 1.0 × 10 −09 ). Conclusions Our study suggests that nonsynonymous variation in the gene NKPD1 affects depressive symptoms in the general population. NKPD1 is predicted to be involved in the de novo synthesis of sphingolipids, which have been implicated in the pathogenesis of depression.

BACKGROUND & AIMS: Hirschsprung disease (HSCR) is an inherited congenital disorder characterized by absence of enteric ganglia in the distal part of the gut. Variants in ret proto-oncogene (RET) have been associated with up to 50% of familial and 35% of sporadic cases. We searched for variants that affect disease risk in a large, multigenerational family with history of HSCR in a linkage region previously associated with the disease (4q31.3-q32.3) and exome wide. METHODS: We performed exome sequencing analyses of a family in the Netherlands with 5 members diagnosed with HSCR and 2 members diagnosed with functional constipation. We initially focused on variants in genes located in 4q31.3-q32.3; however, we also performed an exome-wide analysis in which known HSCR or HSCR-associated gene variants predicted to be deleterious were prioritized for further analysis. Candidate genes were expressed in HEK293, COS-7, and Neuro-2a cells and analyzed by luciferase and immunoblot assays. Morpholinos were designed to target exons of candidate genes and injected into 1-cell stage zebrafish embryos. Embryos were allowed to develop and stained for enteric neurons. RESULTS: Within the linkage region, we identified 1 putative splice variant in the lipopolysaccharide responsive beige-like anchor protein gene (LRBA). Functional assays could not confirm its predicted effect on messenger RNA splicing or on expression of the mab-21 like 2 gene (MAB21L2), which is embedded in LRBA. Zebrafish that developed following injection of the lrba morpholino had a shortened body axis and subtle gut morphological defects, but no significant reduction in number of enteric neurons compared with controls. Outside the linkage region, members of 1 branch of the family carried a previously unidentified RET variant or an in-frame deletion in the glial cell line derived neurotrophic factor gene (GDNF), which encodes a ligand of RET. This deletion was located 6 base pairs before the last codon. We also found variants in the Indian hedgehog gene (IHH) and its mediator, the transcription factor GLI family zinc finger 3 (GLI3). When expressed in cells, the RET-P399L variant disrupted protein glycosylation and had altered phosphorylation following activation by GDNF. The deletion in GDNF prevented secretion of its gene product, reducing RET activation, and the IHH-Q51K variant reduced expression of the transcription factor GLI1. Injection of morpholinos that target ihh reduced the number of enteric neurons to 13% +/- 1.4% of control zebrafish. CONCLUSIONS: In a study of a large family with history of HSCR, we identified variants in LRBA, RET, the gene encoding the RET ligand (GDNF), IHH, and a gene encoding a mediator of IHH signaling (GLI3). These variants altered functions of the gene products when expressed in cells and knockout of ihh reduced the number of enteric neurons in the zebrafish gut.

How cytokine‐driven changes in chromatin topology are converted into gene regulatory circuits during inflammation still remains unclear. Here, we show that interleukin (IL)‐1α induces acute and widespread changes in chromatin accessibility via the TAK1 kinase and NF‐κB at regions that are highly enriched for inflammatory disease‐relevant SNPs. Two enhancers in the extended chemokine locus on human chromosome 4 regulate the IL‐1α‐inducible IL8 and CXCL1‐3 genes. Both enhancers engage in dynamic spatial interactions with gene promoters in an IL‐1α/TAK1‐inducible manner. Microdeletions of p65‐binding sites in either of the two enhancers impair NF‐κB recruitment, suppress activation and biallelic transcription of the IL8/CXCL2 genes, and reshuffle higher‐order chromatin interactions as judged by i4C interactome profiles. Notably, these findings support a dominant role of the IL8 “master” enhancer in the regulation of sustained IL‐1α signaling, as well as for IL‐8 and IL‐6 secretion. CRISPR‐guided transactivation of the IL8 locus or cross‐TAD regulation by TNFα‐responsive enhancers in a different model locus supports the existence of complex enhancer hierarchies in response to cytokine stimulation that prime and orchestrate proinflammatory chromatin responses downstream of NF‐κB., ATAC‐seq, native 4C‐seq and CRISPR genome editing reveal functional hierarchies and chromatin looping dynamics of cytokine‐inducible “master” enhancers that govern gene regulation during the early proinflammatory response.

Developmental and/or epileptic encephalopathies (DEEs) are a group of devastating genetic disorders, resulting in early onset, therapy resistant seizures and developmental delay. Here we report on 12 individuals from 10 families presenting with a severe form of intractable epilepsy, severe developmental delay, progressive microcephaly and visual disturbance. Whole exome sequencing identified a recurrent, homozygous variant (chr2:64083454A>G) in the essentialUDP-glucose pyrophosphorylase(UGP2) gene in all probands. This rare variant results in a tolerable Met12Val missense change of the longer UGP2 protein isoform but causes a disruption of the start codon of the shorter isoform. We show that the absence of the shorter isoform leads to a reduction of functional UGP2 enzyme in brain cell types, leading to altered glycogen metabolism, upregulated unfolded protein response and premature neuronal differentiation, as modelled during pluripotent stem cell differentiationin vitro. In contrast, the complete lack of all UGP2 isoforms leads to differentiation defects in multiple lineages in human cells. Reduced expression of Ugp2a/Ugp2bin vivoin zebrafish mimics visual disturbance and mutant animals show a behavioral phenotype. Our study identifies a recurrent start codon mutation inUGP2as a cause of a novel autosomal recessive DEE. Importantly, it also shows that isoform specific start-loss mutations causing expression loss of a tissue relevant isoform of an essential protein can cause a genetic disease, even when an organism-wide protein absence is incompatible with life. We provide additional examples where a similar disease mechanism applies.

Background Aging is known to induce immunosenescence, resulting in alterations in both the innate and adaptive immune system. Here we evaluated the effects of aging on B cell subsets in peripheral blood of 155 immunologically healthy individuals in four age categories (range 20-95y) via multi-parameter flow cytometry. Furthermore, we studied the naive and antigen-experienced B cell receptor (BCR) repertoire of different age groups and compared it to the clonal BCR repertoire of chronic lymphocytic leukemia (CLL), a disease typically presenting in elderly individuals. Results Total numbers and relative frequencies of B cells were found to decline upon aging, with reductions in transitional B cells, memory cell types, and plasma blasts in the 70 + y group. The BCR repertoire of naive mature B cells and antigen-experienced B cells did not clearly alter until age 70y. Clear changes in IGHV gene usage were observed in naive mature B cells of 70 + y individuals, with a transitional pattern in the 50-70y group. IGHV gene usage of naive mature B cells of the 50-70y, but not the 70 + y, age group resembled that of both younger (50-70y) and older (70 + y) CLL patients. Additionally, CLL-associated stereotypic BCR were found as part of the healthy control BCR repertoire, with an age-associated increase in frequency of several stereotypic BCR (particularly subsets #2 and #5). Conclusion Composition of the peripheral B cell compartment changes with ageing, with clear reductions in non-switched and CD27 + IgG+ switched memory B cells and plasma blasts in especially the 70 + y group. The BCR repertoire is relatively stable until 70y, whereafter differences in IGHV gene usage are seen. Upon ageing, an increasing trend in the occurrence of particular CLL-associated stereotypic BCR is observed. Electronic supplementary material The online version of this article (10.1186/s12979-019-0163-x) contains supplementary material, which is available to authorized users.

Objective Recessive mutations in the Gap Junction Protein Gamma 2 (GJC2) gene cause Pelizaeus-Merzbacher-like disease type 1, a severe infantile-onset hypomyelinating leukodystrophy. Milder, late-onset phenotypes including complicated spastic paraplegia in one family (SPG44), and mild tremor in one case, were reported associated to GJC2 homozygous missense mutations. Here, we report a new family with two siblings carrying a different homozygous GJC2 mutation, presenting with late-onset ataxic and pyramidal disturbances, and parkinsonism in one of them. Methods Two affected siblings were studied by neurological examination and brain MRI. Genetic analyses included genome-wide homozygosity mapping in both siblings, and whole exome sequencing in one sib. The resulting candidate gene variant was validated by Sanger sequencing. Results The affected siblings share a novel homozygous GJC2 missense mutation (c.820G>C, p.Val274Leu), predicted as pathogenic by all used in-silico tools. Brain MRI showed hyperintense signal in T2-weighted images in the internal capsule and subcortical and periventricular white matter, consistent with hypomyelination. Conclusions Our findings confirm and further expand the late-onset phenotypes of GJC2 mutations, to include prominent ataxia, pyramidal disturbances and mild parkinsonism, and confirm the distinctive associated MRI pattern.

Cohesin is a ring-shaped multiprotein complex that is crucial for 3D genome organization and transcriptional regulation during differentiation and development. It also confers sister chromatid cohesion and facilitates DNA damage repair. Besides its core subunits SMC3, SMC1A and RAD21, cohesin contains in somatic cells one of two orthologous STAG subunits, SA1 or SA2. How these variable subunits affect the function of the cohesin complex is still unclear. SA1- and SA2-cohesin were initially proposed to organize cohesion at telomeres and centromeres, respectively. Here, we uncover redundant and specific roles of SA1 and SA2 in gene regulation and chromatin looping using HCT116 cells with an auxin-inducible degron (AID) tag fused to either SA1 or SA2. Following rapid depletion of either subunit, we perform high resolution Hi-C, RNA-sequencing and sequential ChIP studies to show that SA1 and SA2 do not co-occupy individual binding sites and have distinct ways how they affect looping and gene expression. These findings are supported at the single cell level by single-molecule localizations via dSTORM super-resolution imaging. Since somatic and congenital mutations of the SA subunits are associated with cancer (SA2) and intellectual disability syndromes with congenital abnormalities (SA1 and SA2), we verified SA1-/SA2-dependencies using human neural stem cells, hence highlighting their importance for understanding particular disease contexts.

Krüppel-like factor 1 (KLF1) is an essential transcription factor for erythroid development, as demonstrated by Klf1 knockout mice which die around E14 due to severe anemia. In humans, >140 KLF1 variants, causing different erythroid phenotypes, have been described. The KLF1 Nan variant, a single amino acid substitution (p.E339D) in the DNA binding domain, causes hemolytic anemia and is dominant over wildtype KLF1. Here we describe the effects of the KLF1 Nan variant during fetal development. We show that Nan embryos have defects in erythroid maturation. RNA-sequencing of the KLF1 Nan fetal liver cells revealed that Exportin 7 (Xpo7) was among the 782 deregulated genes. This nuclear exportin is implicated in terminal erythroid differentiation; in particular it is involved in nuclear condensation. Indeed, KLF1 Nan fetal liver cells had larger nuclei and reduced chromatin condensation. Knockdown of XPO7 in wildtype erythroid cells caused a similar phenotype. We propose that reduced expression of XPO7 is partially responsible for the erythroid defects observed in KLF1 Nan erythroid cells.

Mammalian cytoplasmic linker associated protein 1 and -2 (CLASP1 and -2) are microtubule (MT) plus-end tracking proteins that selectively stabilize MTs at the edge of cells and that promote MT nucleation and growth at the Golgi, thereby sustaining cell polarity. In vitro analysis has shown that CLASPs are MT growth promoting factors. To date, a single CLASP1 isoform (called CLASP1α) has been described, whereas three CLASP2 isoforms are known (CLASP2α, -β, and -γ). Although CLASP2β/γ are enriched in neurons, suggesting isoform-specific functions, it has been proposed that during neurite outgrowth CLASP1 and -2 act in a redundant fashion by modulating MT dynamics downstream of glycogen synthase kinase 3 (GSK3). Here, we show that in differentiating N1E-115 neuroblastoma cells CLASP1 and CLASP2 differ in their accumulation at MT plus-ends and display different sensitivity to GSK3-mediated phosphorylation, and hence regulation. More specifically, western blot (WB) analysis suggests that pharmacological inhibition of GSK3 affects CLASP2 but not CLASP1 phosphorylation and fluorescence-based microscopy data show that GSK3 inhibition leads to an increase in the number of CLASP2-decorated MT ends, as well as to increased CLASP2 staining of individual MT ends, whereas a reduction in the number of CLASP1-decorated ends is observed. Thus, in N1E-115 cells CLASP2 appears to be a prominent target of GSK3 while CLASP1 is less sensitive. Surprisingly, knockdown of either CLASP causes phosphorylation of GSK3, pointing to the existence of feedback loops between CLASPs and GSK3. In addition, CLASP2 depletion also leads to the activation of protein kinase C (PKC). We found that these differences correlate with opposite functions of CLASP1 and CLASP2 during neuronal differentiation, i.e., CLASP1 stimulates neurite extension, whereas CLASP2 inhibits it. Consistent with knockdown results in N1E-115 cells, primary Clasp2 knockout (KO) neurons exhibit early accelerated neurite and axon outgrowth, showing longer axons than control neurons. We propose a model in which neurite outgrowth is fine-tuned by differentially posttranslationally modified isoforms of CLASPs acting at distinct intracellular locations, thereby targeting MT stabilizing activities of the CLASPs and controlling feedback signaling towards upstream kinases. In summary, our findings provide new insight into the roles of neuronal CLASPs, which emerge as regulators acting in different signaling pathways and locally modulating MT behavior during neurite/axon outgrowth., European Reintegration Grant (ERG) from the European Commission (FP6-2002-MOBILITY-11, grant agreement: ID: 517582)

Parkinson's disease, the most common age-related movement disorder, is a progressive neurodegenerative disease with unclear etiology. Key neuropathological hallmarks are Lewy bodies and Lewy neurites: neuronal inclusions immunopositive for the protein alpha-synuclein. In-depth ultrastructural analysis of Lewy pathology is crucial to understanding pathogenesis of this disease. Using correlative light and electron microscopy and tomography on postmortem human brain tissue from Parkinson's disease brain donors, we identified alpha-synuclein immunopositive Lewy pathology and show a crowded environment of membranes therein, including vesicular structures and dysmorphic organelles. Filaments interspersed between the membranes and organelles were identifiable in many but not all alpha-synuclein inclusions. Crowding of organellar components was confirmed by stimulated emission depletion (STED)-based super-resolution microscopy, and high lipid content within alpha-synuclein immunopositive inclusions was corroborated by confocal imaging, Fourier-transform coherent anti-Stokes Raman scattering infrared imaging and lipidomics. Applying such correlative high-resolution imaging and biophysical approaches, we discovered an aggregated protein-lipid compartmentalization not previously described in the Parkinsons' disease brain.

Background 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. Methods Clinical evaluation, MRI imaging, genome-wide linkage, exome sequencing, urine metabolite profiling, and protein expression studies were performed. Results The first sibling is 17 years old and presents with generalized dystonia and severe bilateral pallidal MRI lesions after 1 episode of infantile subacute metabolic encephalopathy (Leigh-like syndrome). In contrast, the younger sibling (15 years old) only suffers from paroxysmal exercise-induced dystonia and has very mild pallidal MRI abnormalities. Both patients carry compound heterozygous ECHS1 mutations: c.232G>T (predicted protein effect: p.Glu78Ter) and c.518C>T (p.Ala173Val). Linkage analysis, exome sequencing, cosegregation, expression studies, and metabolite profiling support the pathogenicity of these mutations. Expression studies in patients' fibroblasts showed mitochondrial localization and severely reduced levels of ECHS1 protein. Increased urinary S-(2-carboxypropyl)cysteine and N-acetyl-S-(2-carboxypropyl)cysteine levels, proposed metabolic markers of this disorder, were documented in both siblings. Sequencing ECHS1 in 30 unrelated patients with paroxysmal dyskinesias revealed no further mutations. Conclusions The phenotype associated with ECHS1 mutations might be milder than reported earlier, compatible with prolonged survival, and also includes isolated paroxysmal exercise-induced dystonia. ECHS1 screening should be considered in patients with otherwise unexplained paroxysmal exercise-induced dystonia, in addition to those with Leigh and Leigh-like syndromes. Diet regimens and detoxifying agents represent potential therapeutic strategies. © 2016 International Parkinson and Movement Disorder Society

BACKGROUND Cardiomyopathies are usually inherited and predominantly affect adults, but they can also present in childhood. Although our understanding of the molecular basis of pediatric cardiomyopathy has improved, the underlying mechanism remains elusive in a substantial proportion of cases.OBJECTIVES This study aimed to identify new genes involved in pediatric cardiomyopathy.METHODS The authors performed homozygosity mapping and whole-exome sequencing in 2 consanguineous families with idiopathic pediatric cardiomyopathy. Sixty unrelated patients with pediatric cardiomyopathy were subsequently screened for mutations in a candidate gene. First-degree relatives were submitted to cardiac screening and cascade genetic testing. Myocardial samples from 2 patients were processed for histological and immunohistochemical studies.RESULTS We identified 5 patients from 3 unrelated families with pediatric cardiomyopathy caused by homozygous truncating mutations in ALPK3, a gene encoding a nuclear kinase that plays an essential role in early differentiation of cardiomyocytes. All patients with biallelic mutations presented with severe hypertrophic and/or dilated cardiomyopathy in utero, at birth, or in early childhood. Three patients died from heart failure within the first week of life. Moreover, 2 of 10 (20%) heterozygous family members showed hypertrophic cardiomyopathy with an atypical distribution of hypertrophy. Deficiency of alpha-kinase 3 has previously been associated with features of both hypertrophic and dilated cardiomyopathy in mice. Consistent with studies in knockout mice, we provide microscopic evidence for intercalated disc remodeling.CONCLUSIONS Biallelic truncating mutations in the newly identified gene ALPK3 give rise to severe, early-onset cardiomyopathy in humans. Our findings highlight the importance of transcription factor pathways in the molecular mechanisms underlying human cardiomyopathies. (C) 2016 by the American College of Cardiology Foundation.

The X and Y sex chromosomes of placental mammals show hallmarks of a tumultuous evolutionary past. The X Chromosome has a rich and conserved gene content, while the Y Chromosome has lost most of its genes. In the Transcaucasian mole vole Ellobius lutescens, the Y Chromosome including Sry has been lost, and both females and males have a 17,X diploid karyotype. Similarly, the closely related Ellobius talpinus, has a 54,XX karyotype in both females and males. Here, we report the sequencing and assembly of the E. lutescens and E. talpinus genomes. The results indicate that the loss of the Y Chromosome in E. lutescens and E. talpinus occurred in two independent events. Four functional homologs of mouse Y-Chromosomal genes were detected in both female and male E. lutescens, of which three were also detected in the E. talpinus genome. One of these is Eif2s3y, known as the only Y-derived gene that is crucial for successful male meiosis. Female and male E. lutescens can carry one and the same X Chromosome with a largely conserved gene content, including all genes known to function in X Chromosome inactivation. The availability of the genomes of these mole vole species provides unique models to study the dynamics of sex chromosome evolution.

Esophageal atresia (EA) and tracheoesophageal fistula (TEF) are relatively frequently occurring foregut malformations. EA/TEF is thought to have a strong genetic component. Not much is known regarding the biological processes disturbed or which cell type is affected in patients. This hampers the detection of the responsible culprits (genetic or environmental) for the origin of these congenital anatomical malformations. Therefore, we examined gene expression patterns in the TEF and compared them to the patterns in esophageal, tracheal and lung control samples. We studied tissue organization and key proteins using immunohistochemistry. There were clear differences between TEF and control samples. Based on the number of differentially expressed genes as well as histological characteristics, TEFs were most similar to normal esophagus. The BMP-signaling pathway, actin cytoskeleton and extracellular matrix pathways are downregulated in TEF. Genes involved in smooth muscle contraction are overexpressed in TEF compared to esophagus as well as trachea. These enriched pathways indicate myofibroblast activated fibrosis. TEF represents a specific tissue type with large contributions of intestinal smooth muscle cells and neurons. All major cell types present in esophagus are present—albeit often structurally disorganized—in TEF, indicating that its etiology should not be sought in cell fate specification.

Krüppel-like factor 1 (KLF1) is an essential transcription factor for erythroid development, as demonstrated by Klf1 knockout mice which die around E14 due to severe anemia. In humans, >65 KLF1 variants, causing different erythroid phenotypes, have been described. The Klf1 Nan variant, a single amino acid substitution (p.E339D) in the DNA binding domain, causes hemolytic anemia and is dominant over wildtype KLF1. Here we describe the effects of the Nan variant during fetal development. We show that Nan embryos have defects in erythroid maturation. RNA-sequencing of the Nan fetal liver cells revealed that Exportin 7 (Xpo7) was among the ~780 deregulated genes. This nuclear exportin is implicated in terminal erythroid differentiation; in particular it is involved in nuclear condensation. Indeed, KLF1 Nan fetal liver cells had larger nuclei and reduced chromatin condensation. Knockdown of XPO7 in wildtype erythroid cells caused a similar phenotype. We conclude that reduced expression of XPO7 is partially responsible for the erythroid defects observed in Nan erythroid cells.

ObjectiveTo identify the clinical characteristics and genetic etiology of a family affected with hereditary spastic paraplegia (HSP).MethodsClinical, genetic, and functional analyses involving genome-wide linkage coupled to whole-exome sequencing in a consanguineous family with complicated HSP.ResultsA homozygous missense mutation was identified in the ACO2 gene (c.1240T>G p.Phe414Val) that segregated with HSP complicated by intellectual disability and microcephaly. Lymphoblastoid cell lines of homozygous carrier patients revealed significantly decreased activity of the mitochondrial aconitase enzyme and defective mitochondrial respiration. ACO2 encodes mitochondrial aconitase, an essential enzyme in the Krebs cycle. Recessive mutations in this gene have been previously associated with cerebellar ataxia.ConclusionsOur findings nominate ACO2 as a disease-causing gene for autosomal recessive complicated HSP and provide further support for the central role of mitochondrial defects in the pathogenesis of HSP.

Context: Despite the well-recognized clinical features resulting from insufficient or excessive thyroid hormone (TH) levels in humans, it is largely unknown which genes are regulated by TH in human tissues. Objective: To study the effect of TH on human gene expression profiles in whole blood, mainly consisting of T3 receptor (TR) a-expressing cells. Methods: We performed next-generation RNA sequencing on whole blood samples from eight athyroid patients (four females) on and after 4 weeks off levothyroxine replacement. Gene expression changes were analyzed through paired differential expression analysis and confirmed in a validation cohort. Weighted gene coexpression network analysis (WGCNA) was applied to identify thyroid state-related networks. Results: We detected 486 differentially expressed genes (fold-change .1.5; multiple testing corrected P value, 0.05), of which 76% were positively and 24% were negatively regulated. Gene ontology (GO) enrichment analysis revealed that three biological processes were significantly overrepresented, of which the process translational elongation showed the highest fold enrichment (7.3-fold, P = 1.8 3 1026). WGCNA analysis independently identified various gene clusters that correlated with thyroid state. Further GO analysis suggested that thyroid state affects platelet function. Conclusions: Changes in thyroid state regulate numerous genes in human whole blood, predominantly TRa-expressing leukocytes. In addition, TH may regulate gene transcripts in platelets.

This study links a noncoding DNA variant to long range regulation of IL-32 isoform expression, modulating susceptibility to HIV., We integrated data obtained from HIV-1 genome-wide association studies with T cell–derived epigenome data and found that the noncoding intergenic variant rs4349147, which is statistically associated with HIV-1 acquisition, is located in a CD4+ T cell–specific deoxyribonuclease I hypersensitive region, suggesting regulatory potential for this variant. Deletion of the rs4349147 element in Jurkat cells strongly reduced expression of interleukin-32 (IL-32), approximately 10-kb upstream, and chromosome conformation capture assays identified a chromatin loop between rs4349147 and the IL-32 promoter validating its function as a long-distance enhancer. We generated single rs4349147-A or rs4349147-G allele clones and demonstrated that IL-32 enhancer activity and interaction with the IL-32 promoter are strongly allele dependent; rs4349147 −/A cells display reduced IL-32 expression and altered chromatin conformation as compared to rs4349147 G/− cells. Moreover, RNA sequencing demonstrated that rs4349147 G/− cells express a lower relative ratio of IL-32α to non-α isoforms than rs4349147 −/A cells and display increased expression of lymphocyte activation factors rendering them more prone to infection with HIV-1. In agreement, in primary CD4+ T cells, both treatment with recombinant IL-32γ (rIL-32γ) but not rIL-32α, and exogenous lentiviral overexpression of IL-32γ or IL-32β but not IL-32α resulted in a proinflammatory T cell cytokine environment concomitant with increased susceptibility to HIV infection. Our data demonstrate that rs4349147-G promotes transcription of non–IL-32α isoforms, generating a proinflammatory environment more conducive to HIV infection. This study provides a mechanistic link between a HIV-associated noncoding DNA variant and the expression of different IL-32 isoforms that display discrete anti-HIV properties.

Neonatal hyperbilirubinemia (NH) is a common finding in newborn babies in Indonesia. Common and rare variants of UGT1A1 have been known to contribute to NH etiology. This study aims to identify UGT1A1 genetic variation and haplotype associated with NH in Indonesian population. DNA was isolated from 116 cases and 115 controls and a targeted-deep sequencing approach was performed on the promoter, UTRs, and exonic regions of UGT1A1. Determining association of common variants and haplotype analysis were performed using PLINK and Haploview. Ten and 4 rare variants were identified in cases and controls, respectively. The UGT1A1 rare variants frequency in cases (5.17%) was higher than that in controls (1.7%). Four of those rare variants in cases (p.Ala61Thr, p.His300Arg, p.Lys407Asn, and p.Tyr514Asn) and three in controls (p.Tyr79X, p.Ala346Val, and p.Thr412Ser) are novel variants. The frequencies of p.Gly71Arg, p.Pro229Gln, and TA7 common variants were not significantly different between cases and controls. A haplotype, consisting of 3 major alleles of 3′ UTRs common variants (rs8330C>G, rs10929303C>T, and rs1042640C>G), was associated with NH incidence (p=0.025) in this population. Using targeted-deep sequencing and haplotype analysis, we identified novel UGT1A1 rare variants and disease-associated haplotype in NH in Indonesian population.

Chromosome conformation capture (3C) and its derivatives (e.g., 4C, 5C and Hi-C) are used to analyze the 3D organization of genomes. We recently developed targeted chromatin capture (T2C), an inexpensive method for studying the 3D organization of genomes, interactomes and structural changes associated with gene regulation, the cell cycle, and cell survival and development. Here, we present the protocol for T2C based on capture, describing all experimental steps and bio-informatic tools in full detail. T2C offers high resolution, a large dynamic interaction frequency range and a high signal-to-noise ratio. Its resolution is determined by the resulting fragment size of the chosen restriction enzyme, which can lead to sub-kilobase-pair resolution. T2C's high coverage allows the identification of the interactome of each individual DNA fragment, which makes binning of reads (often used in other methods) basically unnecessary. Notably, T2C requires low sequencing efforts. T2C also allows multiplexing of samples for the direct comparison of multiple samples. It can be used to study topologically associating domains (TADs), determining their position, shape, boundaries, and intra- and inter-domain interactions, as well as the composition of aggregated loops, interactions between nucleosomes, individual transcription factor binding sites, and promoters and enhancers. T2C can be performed by any investigator with basic skills in molecular biology techniques in ∼7-8 d. Data analysis requires basic expertise in bioinformatics and in Linux and Python environments.

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. © Kasper WJ Derks, Branislav Misovic, Mirjam CGN van den Hout, Christel EM Kockx, Cesar Payan Gomez, Rutger WW Brouwer, Harry Vrieling, Jan HJ Hoeijmakers, Wilfred FJ van IJcken, and Joris Pothof.

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, smeta = 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.

West Nile virus (WNV) and chikungunya virus (CHIKV) are arboviruses that are constantly (re-)emerging and expanding their territory. Both viruses often cause a mild form of disease, but severe forms of the disease can consist of neurological symptoms, most often observed in the elderly and young children, respectively, for which the mechanisms are poorly understood. To further elucidate the mechanisms responsible for end-stage WNV and CHIKV neuroinvasive disease, we used transcriptomics to compare the induction of effector pathways in the brain during the early and late stage of disease in young mice. In addition to the more commonly described cell death pathways such as apoptosis and autophagy, we also found evidence for the differential expression of pyroptosis and necroptosis cell death markers during both WNV and CHIKV neuroinvasive disease. In contrast, no evidence of cell dysfunction was observed, indicating that cell death may be the most important mechanism of disease. Interestingly, there was overlap when comparing immune markers involved in neuroinvasive disease to those seen in neurodegenerative diseases. Nonetheless, further validation studies are needed to determine the activation and involvement of these effector pathways at the end stage of disease. Furthermore, evidence for a strong inflammatory response was found in mice infected with WNV and CHIKV. The transcriptomics profile measured in mice with WNV and CHIKV neuroinvasive disease in our study showed strong overlap with the mRNA profile described in the literature for other viral neuroinvasive diseases. More studies are warranted to decipher the role of cell inflammation and cell death in viral neuroinvasive disease and whether common mechanisms are active in both neurodegenerative and brain infectious diseases.

BackgroundThe dynamic three-dimensional chromatin architecture of genomes and its co-evolutionary connection to its function—the storage, expression, and replication of genetic information—is still one of the central issues in biology. Here, we describe the much debated 3D architecture of the human and mouse genomes from the nucleosomal to the megabase pair level by a novel approach combining selective high-throughput high-resolution chromosomal interaction capture (T2C), polymer simulations, and scaling analysis of the 3D architecture and the DNA sequence.ResultsThe genome is compacted into a chromatin quasi-fibre with ~5 ± 1 nucleosomes/11 nm, folded into stable ~30–100 kbp loops forming stable loop aggregates/rosettes connected by similar sized linkers. Minor but significant variations in the architecture are seen between cell types and functional states. The architecture and the DNA sequence show very similar fine-structured multi-scaling behaviour confirming their co-evolution and the above.ConclusionsThis architecture, its dynamics, and accessibility, balance stability and flexibility ensuring genome integrity and variation enabling gene expression/regulation by self-organization of (in)active units already in proximity. Our results agree with the heuristics of the field and allow “architectural sequencing” at a genome mechanics level to understand the inseparable systems genomic properties.

The pathological consequences of structural variants disrupting 3D genome organization can be difficult to elucidate in vivo due to differences in gene dosage sensitivity between mice and humans. This is illustrated by branchiooculofacial syndrome (BOFS), a rare congenital disorder caused by heterozygous mutations within TFAP2A, a neural crest regulator for which humans, but not mice, are haploinsufficient. Here, we present a BOFS patient carrying a heterozygous inversion with one breakpoint located within a topologically associating domain (TAD) containing enhancers essential for TFAP2A expression in human neural crest cells (hNCCs). Using patient-specific hiPSCs, we show that, although the inversion shuffles the TFAP2A hNCC enhancers with novel genes within the same TAD, this does not result in enhancer adoption. Instead, the inversion disconnects one TFAP2A allele from its cognate enhancers, leading to monoallelic and haploinsufficient TFAP2A expression in patient hNCCs. Our work illustrates the power of hiPSC differentiation to unveil long-range pathomechanisms., The Tomo Šarić laboratory participates in the EBiSC project, which received support from the Innovative Medicines Initiative Joint Undertaking under grant agreement no. 115582, resources of which are composed of financial contributions from the European Union’s Seventh Framework Programme (FP7/2007-2013) and EFPIA companies’ in-kind contributions. K.K.’s work was supported by the Deutsche Forschungsgemeinschaft (DFG) under grant agreement KO 3588/2-1. Work in the Rada-Iglesias laboratory was supported by CMMC intramural funding and grants from Fritz Thyssen Stiftung and Else Kröner Fresenius Stiftung.