A prenatal skin atlas reveals immune regulation of human skin morphogenesis.

Human prenatal skin is populated by innate immune cells, including macrophages, but whether they act solely in immunity or have additional functions in morphogenesis is unclear. Here we assembled a comprehensive multi-omics reference atlas of prenatal human skin (7-17 post-conception weeks), combining single-cell and spatial transcriptomics data, to characterize the microanatomical tissue niches of the skin. This atlas revealed that crosstalk between non-immune and immune cells underpins the formation of hair follicles, is implicated in scarless wound healing and is crucial for skin angiogenesis. We systematically compared a hair-bearing skin organoid (SkO) model derived from human embryonic stem cells and induced pluripotent stem cells to prenatal and adult skin1. The SkO model closely recapitulated in vivo skin epidermal and dermal cell types during hair follicle development and expression of genes implicated in the pathogenesis of genetic hair and skin disorders. However, the SkO model lacked immune cells and had markedly reduced endothelial cell heterogeneity and quantity. Our in vivo prenatal skin cell atlas indicated that macrophages and macrophage-derived growth factors have a role in driving endothelial development. Indeed, vascular network remodelling was enhanced following transfer of autologous macrophages derived from induced pluripotent stem cells into SkO cultures. Innate immune cells are therefore key players in skin morphogenesis beyond their conventional role in immunity, a function they achieve through crosstalk with non-immune cells.

Dual diagnosis of UQCRFS1-related mitochondrial complex III deficiency and recessive GJA8-related cataracts.

Biallelic pathogenic variants in UQCRFS1 underlie a rare form of isolated mitochondrial complex III deficiency associated with lactic acidosis and a distinctive scalp alopecia previously described in two unrelated probands. Here, we describe a participant in the Undiagnosed Diseases Network (UDN) with a dual diagnosis of two autosomal recessive disorders revealed by genome sequencing: UQCRFS1-related mitochondrial complex III deficiency and GJA8-related cataracts. Both pathogenic variants have been reported before: UQCRFS1 (NM_006003.3:c.215-1 G>C, p.Val72_Thr81del10) in a case with mitochondrial complex III deficiency and GJA8 (NM 005267.5:c.736 G>T, p.Glu246*) as a somatic change in aged cornea leading to decreased junctional coupling. A multi-modal approach combining enzyme assays and cellular proteomics analysis provided clear evidence of complex III respiratory chain dysfunction and low abundance of the Rieske iron-sulfur protein, validating the pathogenic effect of the UQCRFS1 variant. This report extends the genotypic and phenotypic spectrum for these two rare disorders and highlights the utility of deep phenotyping and genomics data to achieve diagnosis and insights into rare disease.

The Human Neural Cell Atlas of Zika Infection in developing human brain tissue: viral pathogenesis, innate immunity, and lineage reprogramming.

Zika virus (ZIKV) infection during pregnancy can lead to fetal brain infection and developmental anomalies collectively known as congenital Zika syndrome (CZS). To define the molecular features underlying CZS in a relevant human cell model, we evaluated ZIKV infection and neurodevelopment in primary fetal brain explants and induced pluripotent stem cell-derived mixed neural cultures at single cell resolution. We identified astrocytes as key innate immune sentinel cells detecting ZIKV and producing IFN-ß. In contrast, neural progenitor cells displayed impaired innate immunity and supported high levels of viral replication. ZIKV infection of neurons suppressed differentiation and synaptic signaling networks and programmed a molecular switch from neurogenesis to astrogliogenesis. We identified a universal ZIKV-driven cellular stress response linked to intrinsic apoptosis and regulated by IFN-ß. These findings reveal how innate immune signaling intersects with ZIKV-driven perturbations in cellular function to influence CZS outcomes including neuron developmental dysfunction and apoptotic cell death.

Multi-omics identification of extracellular components of the fetal monkey and human neocortex.

During development, precursor cells are continuously and intimately interacting with their extracellular environment, which guides their ability to generate functional tissues and organs. Much is known about the development of the neocortex in mammals. This information has largely been derived from histological analyses, heterochronic cell transplants, and genetic manipulations in mice, and to a lesser extent from transcriptomic and histological analyses in humans. However, these approaches have not led to a characterization of the extracellular composition of the developing neocortex in any species. Here, using a combination of single-cell transcriptomic analyses from published datasets, and our proteomics and immunohistofluorescence analyses, we provide a more comprehensive and unbiased picture of the early developing fetal neocortex in humans and non-human primates. Our findings provide a starting point for further hypothesis-driven studies on structural and signaling components in the developing cortex that had previously not been identified.

APOE-NOTCH axis governs elastogenesis during human cardiac valve remodeling.

Valve remodeling is a process involving extracellular matrix organization and elongation of valve leaflets. Here, through single-cell RNA sequencing of human fetal valves, we identified an elastin-producing valve interstitial cell (VIC) subtype (apolipoprotein E (APOE)+, elastin-VICs) spatially located underneath valve endothelial cells (VECs) sensing unidirectional flow. APOE knockdown in fetal VICs resulted in profound elastogenesis defects. In valves with pulmonary stenosis (PS), we observed elastin fragmentation and decreased expression of APOE along with other genes regulating elastogenesis. Cell-cell interaction analysis revealed that jagged 1 (JAG1) from unidirectional VECs activates elastogenesis in elastin-VICs through NOTCH2. Similar observations were made in VICs cocultured with VECs under unidirectional flow. Notably, a drastic reduction of JAG1-NOTCH2 was also observed in PS valves. Lastly, we found that APOE controls JAG1-induced NOTCH activation and elastogenesis in VICs through the extracellular signal-regulated kinase pathway. Our study suggests important roles of both APOE and NOTCH in regulating elastogenesis during human valve remodeling.

Developmental control of rod number via a light-dependent retrograde pathway from intrinsically photosensitive retinal ganglion cells.

Photoreception is essential for the development of the visual system, shaping vision's first synapse to cortical development. Here, we find that the lighting environment controls developmental rod apoptosis via Opn4-expressing intrinsically photosensitive retinal ganglion cells (ipRGCs). Using genetics, sensory environment manipulations, and computational approaches, we establish a pathway where light-dependent glutamate released from ipRGCs is detected via a transiently expressed glutamate receptor (Grik3) on rod precursors within the inner retina. Communication between these cells is mediated by hybrid neurites on ipRGCs that sense light before eye opening. These structures span the ipRGC-rod precursor distance over development and contain the machinery for photoreception (Opn4) and neurotransmitter release (Vglut2 & Syp). Assessment of the human gestational retina identifies conserved hallmarks of an ipRGC-to-rod axis, including displaced rod precursors, transient GRIK3 expression, and ipRGCs with deep-projecting neurites. This analysis defines an adaptive retrograde pathway linking the sensory environment to rod precursors via ipRGCs prior to eye opening.

Lipid Nanoparticle-Mediated Delivery of mRNA Into the Mouse and Human Retina and Other Ocular Tissues.

Purpose: Lipid nanoparticles (LNPs) show promise in their ability to introduce mRNA to drive protein expression in specific cell types of the mammalian eye. Here, we examined the ability of mRNA encapsulated in LNPs with two distinct formulations to drive gene expression in mouse and human retina and other ocular tissues. Methods: We introduced mRNA-carrying LNPs into two biological systems. Intravitreal injections were tested to deliver LNPs into the mouse eye. Human retinal pigment epithelium (RPE) and retinal explants were used to assess mRNA expression in human tissue. We analyzed specificity of expression using histology, immunofluorescence, and imaging. Results: In mice, mRNAs encoding GFP and ciliary neurotrophic factor (CNTF) were specifically expressed by Müller glia and RPE. Acute inflammatory changes measured by microglia distribution (Iba-1) or interleukin-6 (IL-6) expression were not observed 6 hours post-injection. Human RPE also expressed high levels of GFP. Human retinal explants expressed GFP in cells with apical and basal processes consistent with Müller glia and in perivascular cells consistent with macrophages. Conclusions: We demonstrated the ability to reliably transfect subpopulations of retinal cells in mouse eye tissues in vivo and in human ocular tissues. Of significance, intravitreal injections were sufficient to transfect the RPE in mice. To our knowledge, we demonstrate delivery of mRNA using LNPs in human ocular tissues for the first time. Translational Relevance: Ocular gene-replacement therapies using non-viral vector methods are a promising alternative to adeno-associated virus (AAV) vectors. Our studies show that mRNA LNP delivery can be used to transfect retinal cells in both mouse and human tissues without inducing significant inflammation. This methodology could be used to transfect retinal cell lines, tissue explants, mice, or potentially as gene-replacement therapy in a clinical setting in the future.

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.

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.

De novo variants in PLCG1 are associated with hearing impairment, ocular pathology, and cardiac defects.

Phospholipase C isozymes (PLCs) hydrolyze phosphatidylinositol 4,5-bisphosphate into inositol 1,4,5-trisphosphate and diacylglycerol, important signaling molecules involved in many cellular processes. PLCG1 encodes the PLCγ1 isozyme that is broadly expressed. Hyperactive somatic mutations of PLCG1 are observed in multiple cancers, but only one germline variant has been reported. Here we describe three unrelated individuals with de novo heterozygous missense variants in PLCG1 (p.Asp1019Gly, p.His380Arg, and p.Asp1165Gly) who exhibit variable phenotypes including hearing loss, ocular pathology and cardiac septal defects. To model these variants in vivo, we generated the analogous variants in the Drosophila ortholog, small wing (sl). We created a null allele slT2A and assessed the expression pattern. sl is broadly expressed, including in wing discs, eye discs, and a subset of neurons and glia. Loss of sl causes wing size reductions, ectopic wing veins and supernumerary photoreceptors. We document that mutant flies exhibit a reduced lifespan and age-dependent locomotor defects. Expressing wild-type sl in slT2A mutant rescues the loss-of-function phenotypes whereas expressing the variants causes lethality. Ubiquitous overexpression of the variants also reduces viability, suggesting that the variants are toxic. Ectopic expression of an established hyperactive PLCG1 variant (p.Asp1165His) in the wing pouch causes severe wing phenotypes, resembling those observed with overexpression of the p.Asp1019Gly or p.Asp1165Gly variants, further arguing that these two are gain-of-function variants. However, the wing phenotypes associated with p.His380Arg overexpression are mild. Our data suggest that the PLCG1 de novo heterozygous missense variants are pathogenic and contribute to the features observed in the probands.

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.

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.

Sweat gland development requires an eccrine dermal niche and couples two epidermal programs.

Eccrine sweat glands are indispensable for human thermoregulation and, similar to other mammalian skin appendages, form from multipotent epidermal progenitors. Limited understanding of how epidermal progenitors specialize to form these vital organs has precluded therapeutic efforts toward their regeneration. Herein, we applied single-nucleus transcriptomics to compare the expression content of wild-type, eccrine-forming mouse skin to that of mice harboring a skin-specific disruption of Engrailed 1 (En1), a transcription factor that promotes eccrine gland formation in humans and mice. We identify two concurrent but disproportionate epidermal transcriptomes in the early eccrine anlagen: one that is shared with hair follicles and one that is En1 dependent and eccrine specific. We demonstrate that eccrine development requires the induction of a dermal niche proximal to each developing gland in humans and mice. Our study defines the signatures of eccrine identity and uncovers the eccrine dermal niche, setting the stage for targeted regeneration and comprehensive skin repair.

A Spatiotemporal and Machine-Learning Platform Accelerates the Manufacturing of hPSC-derived Esophageal Mucosa.

Human pluripotent stem cell-derived tissue engineering offers great promise in designer cell-based personalized therapeutics. To harness such potential, a broader approach requires a deeper understanding of tissue-level interactions. We previously developed a manufacturing system for the ectoderm-derived skin epithelium for cell replacement therapy. However, it remains challenging to manufacture the endoderm-derived esophageal epithelium, despite both possessing similar stratified structure. Here we employ single cell and spatial technologies to generate a spatiotemporal multi-omics cell atlas for human esophageal development. We illuminate the cellular diversity, dynamics and signal communications for the developing esophageal epithelium and stroma. Using the machine-learning based Manatee, we prioritize the combinations of candidate human developmental signals for in vitro derivation of esophageal basal cells. Functional validation of the Manatee predictions leads to a clinically-compatible system for manufacturing human esophageal mucosa. Our approach creates a versatile platform to accelerate human tissue manufacturing for future cell replacement therapies to treat human genetic defects and wounds.

Species-specific FMRP regulation of RACK1 is critical for prenatal cortical development.

Fragile X messenger ribonucleoprotein 1 protein (FMRP) deficiency leads to fragile X syndrome (FXS), an autism spectrum disorder. The role of FMRP in prenatal human brain development remains unclear. Here, we show that FMRP is important for human and macaque prenatal brain development. Both FMRP-deficient neurons in human fetal cortical slices and FXS patient stem cell-derived neurons exhibit mitochondrial dysfunctions and hyperexcitability. Using multiomics analyses, we have identified both FMRP-bound mRNAs and FMRP-interacting proteins in human neurons and unveiled a previously unknown role of FMRP in regulating essential genes during human prenatal development. We demonstrate that FMRP interaction with CNOT1 maintains the levels of receptor for activated C kinase 1 (RACK1), a species-specific FMRP target. Genetic reduction of RACK1 leads to both mitochondrial dysfunctions and hyperexcitability, resembling FXS neurons. Finally, enhancing mitochondrial functions rescues deficits of FMRP-deficient cortical neurons during prenatal development, demonstrating targeting mitochondrial dysfunction as a potential treatment.

Synchronized long-read genome, methylome, epigenome, and transcriptome for resolving a Mendelian condition.

Resolving the molecular basis of a Mendelian condition (MC) remains challenging owing to the diverse mechanisms by which genetic variants cause disease. To address this, we developed a synchronized long-read genome, methylome, epigenome, and transcriptome sequencing approach, which enables accurate single-nucleotide, insertion-deletion, and structural variant calling and diploid de novo genome assembly, and permits the simultaneous elucidation of haplotype-resolved CpG methylation, chromatin accessibility, and full-length transcript information in a single long-read sequencing run. Application of this approach to an Undiagnosed Diseases Network (UDN) participant with a chromosome X;13 balanced translocation of uncertain significance revealed that this translocation disrupted the functioning of four separate genes (NBEA, PDK3, MAB21L1, and RB1) previously associated with single-gene MCs. Notably, the function of each gene was disrupted via a distinct mechanism that required integration of the four 'omes' to resolve. These included nonsense-mediated decay, fusion transcript formation, enhancer adoption, transcriptional readthrough silencing, and inappropriate X chromosome inactivation of autosomal genes. Overall, this highlights the utility of synchronized long-read multi-omic profiling for mechanistically resolving complex phenotypes.

Developmental adaptation of rod photoreceptor number via photoreception in melanopsin (OPN4) retinal ganglion cells.

Photoreception, a form of sensory experience, is essential for normal development of the mammalian visual system. Detecting photons during development is a prerequisite for visual system function - from vision's first synapse at the cone pedicle and maturation of retinal vascular networks, to transcriptional establishment and maturation of cell types within the visual cortex. Consistent with this theme, we find that the lighting environment regulates developmental rod photoreceptor apoptosis via OPN4-expressing intrinsically photosensitive retinal ganglion cells (ipRGCs). Using a combination of genetics, sensory environment manipulations, and computational approaches, we establish a molecular pathway in which light-dependent glutamate release from ipRGCs is detected via a transiently expressed kainate receptor (GRIK3) in immature rods localized to the inner retina. Communication between ipRGCs and nascent inner retinal rods appears to be mediated by unusual hybrid neurites projecting from ipRGCs that sense light before eye-opening. These structures, previously referred to as outer retinal dendrites (ORDs), span the ipRGC-immature rod distance over the first postnatal week and contain the machinery for sensory detection (melanopsin, OPN4) and axonal/anterograde neurotransmitter release (Synaptophysin, and VGLUT2). Histological and computational assessment of human mid-gestation development reveal conservation of several hallmarks of an ipRGC-to-immature rod pathway, including displaced immature rods, transient GRIK3 expression in the rod lineage, and the presence of ipRGCs with putative neurites projecting deep into the developing retina. Thus, this analysis defines a retinal retrograde signaling pathway that links the sensory environment to immature rods via ipRGC photoreceptors, allowing the visual system to adapt to distinct lighting environments priory to eye-opening.

Opposing roles for TGFß- and BMP-signaling during nascent alveolar differentiation in the developing human lung.

Alveolar type 2 (AT2) cells function as stem cells in the adult lung and aid in repair after injury. The current study aimed to understand the signaling events that control differentiation of this therapeutically relevant cell type during human development. Using lung explant and organoid models, we identified opposing effects of TGFß- and BMP-signaling, where inhibition of TGFß- and activation of BMP-signaling in the context of high WNT- and FGF-signaling efficiently differentiated early lung progenitors into AT2-like cells in vitro. AT2-like cells differentiated in this manner exhibit surfactant processing and secretion capabilities, and long-term commitment to a mature AT2 phenotype when expanded in media optimized for primary AT2 culture. Comparing AT2-like cells differentiated with TGFß-inhibition and BMP-activation to alternative differentiation approaches revealed improved specificity to the AT2 lineage and reduced off-target cell types. These findings reveal opposing roles for TGFß- and BMP-signaling in AT2 differentiation and provide a new strategy to generate a therapeutically relevant cell type in vitro.

Preferential FGF18/FGFR activity in pseudoglandular versus canalicular stage human lung fibroblasts.

Fibroblast growth factor (FGF) signaling is necessary for proper lung branching morphogenesis, alveolarization, and vascular development. Dysregulation of FGF activity has been implicated in various lung diseases. Recently, we showed that FGF18 promotes human lung branching morphogenesis by regulating mesenchymal progenitor cells. However, the underlying mechanisms remain unclear. Thus, we aimed to determine the role of FGF18 and its receptors (FGFR) in regulating mesenchymal cell proliferation, migration, and differentiation from pseudoglandular to canalicular stage. We performed siRNA assays to identify the specific FGFR(s) associated with FGF18-induced biological processes. We found that FGF18 increased proliferation and migration in human fetal lung fibroblasts (HFLF) from both stages. FGFR2/FGFR4 played a significant role in pseudoglandular stage. HFLF proliferation, while FGFR3/FGFR4 were involved in canalicular stage. FGF18 enhanced HFLF migration through FGFR2 and FGFR4 in pseudoglandular and canalicular stage, respectively. Finally, we provide evidence that FGF18 treatment leads to reduced expression of myofibroblast markers (ACTA2 and COL1A1) and increased expression of lipofibroblast markers (ADRP and PPARγ) in both stages HFLF. However, the specific FGF18/FGFR complex involved in this process varies depending on the stage. Our findings suggest that in context of human lung development, FGF18 tends to associate with distinct FGFRs to initiate specific biological processes on mesenchymal cells.