Charting human development using a multi-endodermal organ atlas and organoid models.

Organs are composed of diverse cell types that traverse transient states during organogenesis. To interrogate this diversity during human development, we generate a single-cell transcriptome atlas from multiple developing endodermal organs of the respiratory and gastrointestinal tract. We illuminate cell states, transcription factors, and organ-specific epithelial stem cell and mesenchyme interactions across lineages. We implement the atlas as a high-dimensional search space to benchmark human pluripotent stem cell (hPSC)-derived intestinal organoids (HIOs) under multiple culture conditions. We show that HIOs recapitulate reference cell states and use HIOs to reconstruct the molecular dynamics of intestinal epithelium and mesenchyme emergence. We show that the mesenchyme-derived niche cue NRG1 enhances intestinal stem cell maturation in vitro and that the homeobox transcription factor CDX2 is required for regionalization of intestinal epithelium and mesenchyme in humans. This work combines cell atlases and organoid technologies to understand how human organ development is orchestrated.

Human microglia maturation is underpinned by specific gene regulatory networks.

Microglia phenotypes are highly regulated by the brain environment, but the transcriptional networks that specify the maturation of human microglia are poorly understood. Here, we characterized stage-specific transcriptomes and epigenetic landscapes of fetal and postnatal human microglia and acquired corresponding data in induced pluripotent stem cell (iPSC)-derived microglia, in cerebral organoids, and following engraftment into humanized mice. Parallel development of computational approaches that considered transcription factor (TF) co-occurrence and enhancer activity allowed prediction of shared and state-specific gene regulatory networks associated with fetal and postnatal microglia. Additionally, many features of the human fetal-to-postnatal transition were recapitulated in a time-dependent manner following the engraftment of iPSC cells into humanized mice. These data and accompanying computational approaches will facilitate further efforts to elucidate mechanisms by which human microglia acquire stage- and disease-specific phenotypes.

Unified rhombic lip origins of group 3 and group 4 medulloblastoma.

Medulloblastoma, a malignant childhood cerebellar tumour, segregates molecularly into biologically distinct subgroups, suggesting that a personalized approach to therapy would be beneficial1. Mouse modelling and cross-species genomics have provided increasing evidence of discrete, subgroup-specific developmental origins2. However, the anatomical and cellular complexity of developing human tissues3-particularly within the rhombic lip germinal zone, which produces all glutamatergic neuronal lineages before internalization into the cerebellar nodulus-makes it difficult to validate previous inferences that were derived from studies in mice. Here we use multi-omics to resolve the origins of medulloblastoma subgroups in the developing human cerebellum. Molecular signatures encoded within a human rhombic-lip-derived lineage trajectory aligned with photoreceptor and unipolar brush cell expression profiles that are maintained in group 3 and group 4 medulloblastoma, suggesting a convergent basis. A systematic diagnostic-imaging review of a prospective institutional cohort localized the putative anatomical origins of group 3 and group 4 tumours to the nodulus. Our results connect the molecular and phenotypic features of clinically challenging medulloblastoma subgroups to their unified beginnings in the rhombic lip in the early stages of human development.

De novo missense variants in exon 9 of SEPHS1 cause a neurodevelopmental condition with developmental delay, poor growth, hypotonia, and dysmorphic features.

Selenophosphate synthetase (SEPHS) plays an essential role in selenium metabolism. Two mammalian SEPHS paralogues, SEPHS1 and SEPHS2, share high sequence identity and structural homology with SEPHS. Here, we report nine individuals from eight families with developmental delay, growth and feeding problems, hypotonia, and dysmorphic features, all with heterozygous missense variants in SEPHS1. Eight of these individuals had a recurrent variant at amino acid position 371 of SEPHS1 (p.Arg371Trp, p.Arg371Gln, and p.Arg371Gly); seven of these variants were known to be de novo. Structural modeling and biochemical assays were used to understand the effect of these variants on SEPHS1 function. We found that a variant at residue Trp352 results in local structural changes of the C-terminal region of SEPHS1 that decrease the overall thermal stability of the enzyme. In contrast, variants of a solvent-exposed residue Arg371 do not impact enzyme stability and folding but could modulate direct protein-protein interactions of SEPSH1 with cellular factors in promoting cell proliferation and development. In neuronal SH-SY5Y cells, we assessed the impact of SEPHS1 variants on cell proliferation and ROS production and investigated the mRNA expression levels of genes encoding stress-related selenoproteins. Our findings provided evidence that the identified SEPHS1 variants enhance cell proliferation by modulating ROS homeostasis. Our study supports the hypothesis that SEPHS1 plays a critical role during human development and provides a basis for further investigation into the molecular mechanisms employed by SEPHS1. Furthermore, our data suggest that variants in SEPHS1 are associated with a neurodevelopmental disorder.

Retinoic acid signaling regulates spatiotemporal specification of human green and red cones.

Trichromacy is unique to primates among placental mammals, enabled by blue (short/S), green (medium/M), and red (long/L) cones. In humans, great apes, and Old World monkeys, cones make a poorly understood choice between M and L cone subtype fates. To determine mechanisms specifying M and L cones, we developed an approach to visualize expression of the highly similar M- and L-opsin mRNAs. M-opsin was observed before L-opsin expression during early human eye development, suggesting that M cones are generated before L cones. In adult human tissue, the early-developing central retina contained a mix of M and L cones compared to the late-developing peripheral region, which contained a high proportion of L cones. Retinoic acid (RA)-synthesizing enzymes are highly expressed early in retinal development. High RA signaling early was sufficient to promote M cone fate and suppress L cone fate in retinal organoids. Across a human population sample, natural variation in the ratios of M and L cone subtypes was associated with a noncoding polymorphism in the NR2F2 gene, a mediator of RA signaling. Our data suggest that RA promotes M cone fate early in development to generate the pattern of M and L cones across the human retina.

Regionally distinct progenitor cells in the lower airway give rise to neuroendocrine and multiciliated cells in the developing human lung.

Using scRNA-seq and microscopy, we describe a cell that is enriched in the lower airways of the developing human lung and identified by the unique coexpression of SCGB3A2/SFTPB/CFTR. To functionally interrogate these cells, we apply a single-cell barcode-based lineage tracing method, called CellTagging, to track the fate of SCGB3A2/SFTPB/CFTR cells during airway organoid differentiation in vitro. Lineage tracing reveals that these cells have a distinct differentiation potential from basal cells, giving rise predominantly to pulmonary neuroendocrine cells and a subset of multiciliated cells distinguished by high C6 and low MUC16 expression. Lineage tracing results are supported by studies using organoids and isolated cells from the lower noncartilaginous airway. We conclude that SCGB3A2/SFTPB/CFTR cells are enriched in the lower airways of the developing human lung and contribute to the epithelial diversity and heterogeneity in this region.

MerTK is a mediator of alpha-synuclein fibril uptake by human microglia.

Mer tyrosine kinase (MerTK) is a receptor tyrosine kinase that mediates non-inflammatory, homeostatic phagocytosis of diverse types of cellular debris. Highly expressed on the surface of microglial cells, MerTK is of importance in brain development, homeostasis, plasticity and disease. Yet, involvement of this receptor in the clearance of protein aggregates that accumulate with ageing and in neurodegenerative diseases has yet to be defined. The current study explored the function of MerTK in the microglial uptake of alpha-synuclein fibrils which play a causative role in the pathobiology of synucleinopathies. Using human primary and induced pluripotent stem cell-derived microglia, the MerTK-dependence of alpha-synuclein fibril internalization was investigated in vitro. Relevance of this pathway in synucleinopathies was assessed through burden analysis of MERTK variants and analysis of MerTK expression in patient-derived cells and tissues. Pharmacological inhibition of MerTK and siRNA-mediated MERTK knockdown both caused a decreased rate of alpha-synuclein fibril internalization by human microglia. Consistent with the non-inflammatory nature of MerTK-mediated phagocytosis, alpha-synuclein fibril internalization was not observed to induce secretion of pro-inflammatory cytokines such as IL-6 or TNF, and downmodulated IL-1ß secretion from microglia. Burden analysis in two independent patient cohorts revealed a significant association between rare functionally deleterious MERTK variants and Parkinson's disease in one of the cohorts (P = 0.002). Despite a small upregulation in MERTK mRNA expression in nigral microglia from Parkinson's disease/Lewy body dementia patients compared to those from non-neurological control donors in a single-nuclei RNA-sequencing dataset (P = 5.08 × 10-21), no significant upregulation in MerTK protein expression was observed in human cortex and substantia nigra lysates from Lewy body dementia patients compared to controls. Taken together, our findings define a novel role for MerTK in mediating the uptake of alpha-synuclein fibrils by human microglia, with possible involvement in limiting alpha-synuclein spread in synucleinopathies such as Parkinson's disease. Upregulation of this pathway in synucleinopathies could have therapeutic values in enhancing alpha-synuclein fibril clearance in the brain.

A dyadic approach to the delineation of diagnostic entities in clinical genomics.

The delineation of disease entities is complex, yet recent advances in the molecular characterization of diseases provide opportunities to designate diseases in a biologically valid manner. Here, we have formalized an approach to the delineation of Mendelian genetic disorders that encompasses two distinct but inter-related concepts: (1) the gene that is mutated and (2) the phenotypic descriptor, preferably a recognizably distinct phenotype. We assert that only by a combinatorial or dyadic approach taking both of these attributes into account can a unitary, distinct genetic disorder be designated. We propose that all Mendelian disorders should be designated as "GENE-related phenotype descriptor" (e.g., "CFTR-related cystic fibrosis"). This approach to delineating and naming disorders reconciles the complexity of gene-to-phenotype relationships in a simple and clear manner yet communicates the complexity and nuance of these relationships.

Targeted long-read sequencing identifies missing disease-causing variation.

Despite widespread clinical genetic testing, many individuals with suspected genetic conditions lack a precise diagnosis, limiting their opportunity to take advantage of state-of-the-art treatments. In some cases, testing reveals difficult-to-evaluate structural differences, candidate variants that do not fully explain the phenotype, single pathogenic variants in recessive disorders, or no variants in genes of interest. Thus, there is a need for better tools to identify a precise genetic diagnosis in individuals when conventional testing approaches have been exhausted. We performed targeted long-read sequencing (T-LRS) using adaptive sampling on the Oxford Nanopore platform on 40 individuals, 10 of whom lacked a complete molecular diagnosis. We computationally targeted up to 151 Mbp of sequence per individual and searched for pathogenic substitutions, structural variants, and methylation differences using a single data source. We detected all genomic aberrations-including single-nucleotide variants, copy number changes, repeat expansions, and methylation differences-identified by prior clinical testing. In 8/8 individuals with complex structural rearrangements, T-LRS enabled more precise resolution of the mutation, leading to changes in clinical management in one case. In ten individuals with suspected Mendelian conditions lacking a precise genetic diagnosis, T-LRS identified pathogenic or likely pathogenic variants in six and variants of uncertain significance in two others. T-LRS accurately identifies pathogenic structural variants, resolves complex rearrangements, and identifies Mendelian variants not detected by other technologies. T-LRS represents an efficient and cost-effective strategy to evaluate high-priority genes and regions or complex clinical testing results.

Differential Regulation of Immune-Related Genes in the Developing Heart.

In many congenital heart defects, it can be difficult to ascertain primary pathology from secondary consequences from altered flow through the developing heart. The molecular differences between the growing right and left ventricles (RV and LV, respectively) following the completion of septation and the impact of sex on these mechanisms have not been investigated. We analyzed RNA-seq data derived from twelve RV and LVs, one with Hypoplastic Left Heart Syndrome (HLHS), to compare the transcriptomic landscape between the ventricles during development. Differential gene expression analysis revealed a large proportion of genes unique to either the RV or LV as well as sex bias. Our GO enrichment and network analysis strategy highlighted the differential role of immune functions between the RV and LV in the developing heart. Comparatively, RNA-seq analysis of data from C57Bl6/J mice hearts collected at E14 resulted in the enrichment of similar processes related to T cells and leukocyte migration and activation. Differential gene expression analysis of an HLHS case highlighted significant downregulation of chromatin organization pathways and upregulation of genes involved in muscle organ development. This analysis also identified previously unreported upregulation of genes involved in IL-17 production pathways. In conclusion, differences exist between the gene expression profiles of RV versus LV with the expression of immune-related genes being significantly different between these two chambers. The pathogenesis of HLHS may involve alterations in the expression of chromatin and muscle gene organization as well as upregulation of the IL-17 response pathway.

FGF18 promotes human lung branching morphogenesis through regulating mesenchymal progenitor cells.

Fibroblast growth factor (FGF) signaling is known to play an important role in lung organogenesis. However, we recently demonstrated that FGF10 fails to induce branching in human fetal lungs as is observed in mouse. Our previous human fetal lung RNA sequencing data exhibited increased FGF18 during the pseudoglandular stage of development, suggestive of its importance in human lung branching morphogenesis. Whereas it has been previously reported that FGF18 is critical during alveologenesis, few studies have described its implication in lung branching, specifically in human. Therefore, we aimed to determine the role of FGF18 in human lung branching morphogenesis. Human fetal lung explants within the pseudoglandular stage of development were treated with recombinant human FGF18 in air-liquid interface culture. Explants were analyzed grossly to assess differences in branching pattern, as well as at the cellular and molecular levels. FGF18 treatment promoted branching in explant cultures and demonstrated increased epithelial proliferation as well as maintenance of the double positive SOX2/SOX9 distal bud progenitor cells, confirming its role in human lung branching morphogenesis. In addition, FGF18 treated explants displayed increased expression of SOX9, FN1, and COL2A1 within the mesenchyme, all factors that are important to chondrocyte differentiation. In humans, cartilaginous airways extend deep into the lung up to the 12th generation of branching whereas in mouse these are restricted to the trachea and main bronchi. Therefore, our data suggest that FGF18 promotes human lung branching morphogenesis through regulating mesenchymal progenitor cells.

The contribution of mosaicism to genetic diseases and de novo pathogenic variants.

The contribution of mosaicism to diagnosed genetic disease and presumed de novo variants (DNV) is under investigated. We determined the contribution of mosaic genetic disease (MGD) and diagnosed parental mosaicism (PM) in parents of offspring with reported DNV (in the same variant) in the (1) Undiagnosed Diseases Network (UDN) (N = 1946) and (2) in 12,472 individuals electronic health records (EHR) who underwent genetic testing at an academic medical center. In the UDN, we found 4.51% of diagnosed probands had MGD, and 2.86% of parents of those with DNV exhibited PM. In the EHR, we found 6.03% and 2.99% and (of diagnosed probands) had MGD detected on chromosomal microarray and exome/genome sequencing, respectively. We found 2.34% (of those with a presumed pathogenic DNV) had a parent with PM for the variant. We detected mosaicism (regardless of pathogenicity) in 4.49% of genetic tests performed. We found a broad phenotypic spectrum of MGD with previously unknown phenotypic phenomena. MGD is highly heterogeneous and provides a significant contribution to genetic diseases. Further work is required to improve the diagnosis of MGD and investigate how PM contributes to DNV risk.

Redefining the Etiologic Landscape of Cerebellar Malformations.

Cerebellar malformations are diverse congenital anomalies frequently associated with developmental disability. Although genetic and prenatal non-genetic causes have been described, no systematic analysis has been performed. Here, we present a large-exome sequencing study of Dandy-Walker malformation (DWM) and cerebellar hypoplasia (CBLH). We performed exome sequencing in 282 individuals from 100 families with DWM or CBLH, and we established a molecular diagnosis in 36 of 100 families, with a significantly higher yield for CBLH (51%) than for DWM (16%). The 41 variants impact 27 neurodevelopmental-disorder-associated genes, thus demonstrating that CBLH and DWM are often features of monogenic neurodevelopmental disorders. Though only seven monogenic causes (19%) were identified in more than one individual, neuroimaging review of 131 additional individuals confirmed cerebellar abnormalities in 23 of 27 genetic disorders (85%). Prenatal risk factors were frequently found among individuals without a genetic diagnosis (30 of 64 individuals [47%]). Single-cell RNA sequencing of prenatal human cerebellar tissue revealed gene enrichment in neuronal and vascular cell types; this suggests that defective vasculogenesis may disrupt cerebellar development. Further, de novo gain-of-function variants in PDGFRB, a tyrosine kinase receptor essential for vascular progenitor signaling, were associated with CBLH, and this discovery links genetic and non-genetic etiologies. Our results suggest that genetic defects impact specific cerebellar cell types and implicate abnormal vascular development as a mechanism for cerebellar malformations. We also confirmed a major contribution for non-genetic prenatal factors in individuals with cerebellar abnormalities, substantially influencing diagnostic evaluation and counseling regarding recurrence risk and prognosis.

Paralog Studies Augment Gene Discovery: DDX and DHX Genes.

Members of a paralogous gene family in which variation in one gene is known to cause disease are eight times more likely to also be associated with human disease. Recent studies have elucidated DHX30 and DDX3X as genes for which pathogenic variant alleles are involved in neurodevelopmental disorders. We hypothesized that variants in paralogous genes encoding members of the DExD/H-box RNA helicase superfamily might also underlie developmental delay and/or intellectual disability (DD and/or ID) disease phenotypes. Here we describe 15 unrelated individuals who have DD and/or ID, central nervous system (CNS) dysfunction, vertebral anomalies, and dysmorphic features and were found to have probably damaging variants in DExD/H-box RNA helicase genes. In addition, these individuals exhibit a variety of other tissue and organ system involvement including ocular, outer ear, hearing, cardiac, and kidney tissues. Five individuals with homozygous (one), compound-heterozygous (two), or de novo (two) missense variants in DHX37 were identified by exome sequencing. We identified ten total individuals with missense variants in three other DDX/DHX paralogs: DHX16 (four individuals), DDX54 (three individuals), and DHX34 (three individuals). Most identified variants are rare, predicted to be damaging, and occur at conserved amino acid residues. Taken together, these 15 individuals implicate the DExD/H-box helicases in both dominantly and recessively inherited neurodevelopmental phenotypes and highlight the potential for more than one disease mechanism underlying these disorders.

De Novo ZMYND8 variants result in an autosomal dominant neurodevelopmental disorder with cardiac malformations.

PURPOSE: ZMYND8 encodes a multidomain protein that serves as a central interactive hub for coordinating critical roles in transcription regulation, chromatin remodeling, regulation of super-enhancers, DNA damage response and tumor suppression. We delineate a novel neurocognitive disorder caused by variants in the ZMYND8 gene. METHODS: An international collaboration, exome sequencing, molecular modeling, yeast two-hybrid assays, analysis of available transcriptomic data and a knockdown Drosophila model were used to characterize the ZMYND8 variants. RESULTS: ZMYND8 variants were identified in 11 unrelated individuals; 10 occurred de novo and one suspected de novo; 2 were truncating, 9 were missense, of which one was recurrent. The disorder is characterized by intellectual disability with variable cardiovascular, ophthalmologic and minor skeletal anomalies. Missense variants in the PWWP domain of ZMYND8 abolish the interaction with Drebrin and missense variants in the MYND domain disrupt the interaction with GATAD2A. ZMYND8 is broadly expressed across cell types in all brain regions and shows highest expression in the early stages of brain development. Neuronal knockdown of the DrosophilaZMYND8 ortholog results in decreased habituation learning, consistent with a role in cognitive function. CONCLUSION: We present genomic and functional evidence for disruption of ZMYND8 as a novel etiology of syndromic intellectual disability.

A unique cardiovascular presentation of Marfan syndrome.

We report a neonate with severe Marfan syndrome (MS), prenatally identified to have persistent atrial tachycardia, biventricular dysfunction, and an unusual structure within the atria. Detailed postnatal echocardiographic evaluation and cross-sectional imaging confirmed congenital pseudoaneurysm of the mitral-aortic intervalvular fibrosa. Emergent testing by next-generation sequencing identified a FBN1 pathological variant, key to establishing goals of care. To our knowledge, this is the first reported case of a congenital pseudoaneurysm of the mitral-aortic intervalvular fibrosa in MS.

Hedgehog Signaling Pathway Orchestrates Human Lung Branching Morphogenesis.

The Hedgehog (HH) signaling pathway plays an essential role in mouse lung development. We hypothesize that the HH pathway is necessary for branching during human lung development and is impaired in pulmonary hypoplasia. Single-cell, bulk RNA-sequencing data, and human fetal lung tissues were analyzed to determine the spatiotemporal localization of HH pathway actors. Distal human lung segments were cultured in an air-liquid interface and treated with an SHH inhibitor (5E1) to determine the effect of HH inhibition on human lung branching, epithelial-mesenchymal markers, and associated signaling pathways in vitro. Our results showed an early and regulated expression of HH pathway components during human lung development. Inhibiting HH signaling caused a reduction in branching during development and dysregulated epithelial (SOX2, SOX9) and mesenchymal (ACTA2) progenitor markers. FGF and Wnt pathways were also disrupted upon HH inhibition. Finally, we demonstrated that HH signaling elements were downregulated in lung tissues of patients with a congenital diaphragmatic hernia (CDH). In this study, we show for the first time that HH signaling inhibition alters important genes and proteins required for proper branching of the human developing lung. Understanding the role of the HH pathway on human lung development could lead to the identification of novel therapeutic targets for childhood pulmonary diseases.

Mutations in the Epithelial Cadherin-p120-Catenin Complex Cause Mendelian Non-Syndromic Cleft Lip with or without Cleft Palate.

Non-syndromic cleft lip with or without cleft palate (NS-CL/P) is one of the most common human birth defects and is generally considered a complex trait. Despite numerous loci identified by genome-wide association studies, the effect sizes of common variants are relatively small, with much of the presumed genetic contribution remaining elusive. We report exome-sequencing results in 209 people from 72 multi-affected families with pedigree structures consistent with autosomal-dominant inheritance and variable penetrance. Herein, pathogenic variants are described in four genes encoding components of the p120-catenin complex (CTNND1, PLEKHA7, PLEKHA5) and an epithelial splicing regulator (ESRP2), in addition to the known CL/P-associated gene, CDH1, which encodes E-cadherin. The findings were also validated in a second cohort of 497 people with NS-CL/P, comprising small families and singletons with pathogenic variants in these genes identified in 14% of multi-affected families and 2% of the replication cohort of smaller families. Enriched expression of each gene/protein in human and mouse embryonic oro-palatal epithelia, demonstration of functional impact of CTNND1 and ESRP2 variants, and recapitulation of the CL/P spectrum in Ctnnd1 knockout mice support a causative role in CL/P pathogenesis. These data show that primary defects in regulators of epithelial cell adhesion are the most significant contributors to NS-CL/P identified to date and that inherited and de novo single gene variants explain a substantial proportion of NS-CL/P.

Evidence of disrupted rhombic lip development in the pathogenesis of Dandy-Walker malformation.

Dandy-Walker malformation (DWM) and Cerebellar vermis hypoplasia (CVH) are commonly recognized human cerebellar malformations diagnosed following ultrasound and antenatal or postnatal MRI. Specific radiological criteria are used to distinguish them, yet little is known about their differential developmental disease mechanisms. We acquired prenatal cases diagnosed as DWM and CVH and studied cerebellar morphobiometry followed by histological and immunohistochemical analyses. This was supplemented by laser capture microdissection and RNA-sequencing of the cerebellar rhombic lip, a transient progenitor zone, to assess the altered transcriptome of DWM vs control samples. Our radiological findings confirm that the cases studied fall within the accepted biometric range of DWM. Our histopathological analysis points to reduced foliation and inferior vermian hypoplasia as common features in all examined DWM cases. We also find that the rhombic lip, a dorsal stem cell zone that drives the growth and maintenance of the posterior vermis is specifically disrupted in DWM, with reduced proliferation and self-renewal of the progenitor pool, and altered vasculature, all confirmed by transcriptomics analysis. We propose a unified model for the developmental pathogenesis of DWM. We hypothesize that rhombic lip development is disrupted through either aberrant vascularization and/or direct insult which causes reduced proliferation and failed expansion of the rhombic lip progenitor pool leading to disproportionate hypoplasia and dysplasia of the inferior vermis. Timing of insult to the developing rhombic lip (before or after 14 PCW) dictates the extent of hypoplasia and distinguishes DWM from CVH.