The transition between the proliferation and differentiation of progenitor cells is a key step in organogenesis, and alterations in this process can lead to developmental disorders. The extracellular signal-regulated kinase 1/2 (ERK) signaling pathway is one of the most intensively studied signaling mechanisms that regulates both proliferation and differentiation. How a single molecule (e.g. ERK) can regulate two opposing cellular outcomes is still a mystery. Using both chick and mouse models, we shed light on the mechanism responsible for the switch from proliferation to differentiation of head muscle progenitors and implicate ERK subcellular localization. Manipulation of the fibroblast growth factor (FGF)-ERK signaling pathway in chick embryos in vitro and in vivo demonstrated that blockage of this pathway accelerated myogenic differentiation, whereas its activation diminished it. We next examined whether the spatial subcellular localization of ERK could act as a switch between proliferation (nuclear ERK) and differentiation (cytoplasmic ERK) of muscle progenitors. A myristoylated peptide that blocks importin 7-mediated ERK nuclear translocation induced robust myogenic differentiation of muscle progenitor/stem cells in both head and trunk. In the mouse, analysis of Sprouty mutant embryos revealed that increased ERK signaling suppressed both head and trunk myogenesis. Our findings, corroborated by mathematical modeling, suggest that ERK shuttling between the nucleus and the cytoplasm provides a switch-like transition between proliferation and differentiation of muscle progenitors.

Vitiligo is an autoimmune disease in which depigmented skin results from the destruction of melanocytes, with epidemiological association with other autoimmune diseases. In previous linkage and genome-wide association studies (GWAS1 and GWAS2), we identified 27 vitiligo susceptibility loci in patients of European ancestry. We carried out a third GWAS (GWAS3) in European-ancestry subjects, with augmented GWAS1 and GWAS2 controls, genome-wide imputation, and meta-analysis of all three GWAS, followed by an independent replication. The combined analyses, with 4,680 cases and 39,586 controls, identified 23 new significantly associated loci and 7 suggestive loci. Most encode immune and apoptotic regulators, with some also associated with other autoimmune diseases, as well as several melanocyte regulators. Bioinformatic analyses indicate a predominance of causal regulatory variation, some of which corresponds to expression quantitative trait loci (eQTLs) at these loci. Together, the identified genes provide a framework for the genetic architecture and pathobiology of vitiligo, highlight relationships with other autoimmune diseases and melanoma, and offer potential targets for treatment.

The FaceBase Consortium, funded by the National Institute of Dental and Craniofacial Research, National Institutes of Health, is designed to accelerate understanding of craniofacial developmental biology by generating comprehensive data resources to empower the research community, exploring high-throughput technology, fostering new scientific collaborations among researchers and human/computer interactions, facilitating hypothesis-driven research and translating science into improved health care to benefit patients. The resources generated by the FaceBase projects include a number of dynamic imaging modalities, genome-wide association studies, software tools for analyzing human facial abnormalities, detailed phenotyping, anatomical and molecular atlases, global and specific gene expression patterns, and transcriptional profiling over the course of embryonic and postnatal development in animal models and humans. The integrated data visualization tools, faceted search infrastructure, and curation provided by the FaceBase Hub offer flexible and intuitive ways to interact with these multidisciplinary data. In parallel, the datasets also offer unique opportunities for new collaborations and training for researchers coming into the field of craniofacial studies. Here, we highlight the focus of each spoke project and the integration of datasets contributed by the spokes to facilitate craniofacial research.

Numerous lines of evidence point to a genetic basis for facial morphology in humans, yet little is known about how specific genetic variants relate to the phenotypic expression of many common facial features. We conducted genome-wide association meta-analyses of 20 quantitative facial measurements derived from the 3D surface images of 3118 healthy individuals of European ancestry belonging to two US cohorts. Analyses were performed on just under one million genotyped SNPs (Illumina OmniExpress+Exome v1.2 array) imputed to the 1000 Genomes reference panel (Phase 3). We observed genome-wide significant associations (p < 5 x 10-8) for cranial base width at 14q21.1 and 20q12, intercanthal width at 1p13.3 and Xq13.2, nasal width at 20p11.22, nasal ala length at 14q11.2, and upper facial depth at 11q22.1. Several genes in the associated regions are known to play roles in craniofacial development or in syndromes affecting the face: MAFB, PAX9, MIPOL1, ALX3, HDAC8, and PAX1. We also tested genotype-phenotype associations reported in two previous genome-wide studies and found evidence of replication for nasal ala length and SNPs in CACNA2D3 and PRDM16. These results provide further evidence that common variants in regions harboring genes of known craniofacial function contribute to normal variation in human facial features. Improved understanding of the genes associated with facial morphology in healthy individuals can provide insights into the pathways and mechanisms controlling normal and abnormal facial morphogenesis.

We previously reported a genome-wide association study (GWAS) identifying 14 susceptibility loci for generalized vitiligo. We report here a second GWAS (450 individuals with vitiligo (cases) and 3,182 controls), an independent replication study (1,440 cases and 1,316 controls) and a meta-analysis (3,187 cases and 6,723 controls) identifying 13 additional vitiligo-associated loci. These include OCA2-HERC2 (combined P = 3.80 × 10(-8)), MC1R (P = 1.82 × 10(-13)), a region near TYR (P = 1.57 × 10(-13)), IFIH1 (P = 4.91 × 10(-15)), CD80 (P = 3.78 × 10(-10)), CLNK (P = 1.56 × 10(-8)), BACH2 (P = 2.53 × 10(-8)), SLA (P = 1.58 × 10(-8)), CASP7 (P = 3.56 × 10(-8)), CD44 (P = 1.78 × 10(-9)), IKZF4 (P = 2.75 × 10(-14)), SH2B3 (P = 3.54 × 10(-18)) and TOB2 (P = 6.81 × 10(-10)). Most vitiligo susceptibility loci encode immunoregulatory proteins or melanocyte components that likely mediate immune targeting and the relationships among vitiligo, melanoma, and eye, skin and hair coloration.

Vitiligo is an autoimmune disease in which melanocyte destruction causes skin depigmentation, with 49 loci known from previous GWAS. Aiming to define vitiligo subtypes, we discovered that age-of-onset is bimodal; one-third of cases have early onset (mean 10.3 years) and two-thirds later onset (mean 34.0 years). In the early-onset subgroup we found novel association with MHC class II region indel rs145954018, and independent association with the principal MHC class II locus from previous GWAS, represented by rs9271597; greatest association was with rs145954018del-rs9271597A haplotype (P = 2.40 × 10-86, OR = 8.10). Both rs145954018 and rs9271597 are located within lymphoid-specific enhancers, and the rs145954018del-rs9271597A haplotype is specifically associated with increased expression of HLA-DQB1 mRNA and HLA-DQ protein by monocytes and dendritic cells. Thus, for vitiligo, MHC regulatory variation confers extreme risk, more important than HLA coding variation. MHC regulatory variation may represent a significant component of genetic risk for other autoimmune diseases.

Tissue homeostasis requires the production of newly differentiated cells from resident adult stem cells. Central to this process is the expansion of undifferentiated intermediates known as transit-amplifying (TA) cells, but how stem cells are triggered to enter this proliferative TA state remains an important open question. Using the continuously growing mouse incisor as a model of stem cell-based tissue renewal, we found that the transcriptional cofactors YAP and TAZ are required both to maintain TA cell proliferation and to inhibit differentiation. Specifically, we identified a pathway involving activation of integrin α3 in TA cells that signals through an LATS-independent FAK/CDC42/PP1A cascade to control YAP-S397 phosphorylation and nuclear localization. This leads to Rheb expression and potentiates mTOR signaling to drive the proliferation of TA cells. These findings thus reveal a YAP/TAZ signaling mechanism that coordinates stem cell expansion and differentiation during organ renewal.

Investigations into stem cell-fueled renewal of an organ benefit from an inventory of cell type-specific markers and a deep understanding of the cellular diversity within stem cell niches. Using the adult mouse incisor as a model for a continuously renewing organ, we performed an unbiased analysis of gene co-expression relationships to identify modules of co-expressed genes that represent differentiated cells, transit-amplifying cells, and residents of stem cell niches. Through in vivo lineage tracing, we demonstrated the power of this approach by showing that co-expression module members Lrig1 and Igfbp5 define populations of incisor epithelial and mesenchymal stem cells. We further discovered that two adjacent mesenchymal tissues, the periodontium and dental pulp, are maintained by distinct pools of stem cells. These findings reveal novel mechanisms of incisor renewal and illustrate how gene co-expression analysis of intact biological systems can provide insights into the transcriptional basis of cellular identity.

Stem cells in stratified epithelia are generally believed to adhere to a non-hierarchical single-progenitor model. Using lineage tracing and genetic label-retention assays, we show that the hard palatal epithelium of the oral cavity is unique in displaying marked proliferative heterogeneity. We identify a previously uncharacterized, infrequently-dividing stem cell population that resides within a candidate niche, the junctional zone (JZ). JZ stem cells tend to self-renew by planar symmetric divisions, respond to masticatory stresses, and promote wound healing, whereas frequently-dividing cells reside outside the JZ, preferentially renew through perpendicular asymmetric divisions, and are less responsive to injury. LRIG1 is enriched in the infrequently-dividing population in homeostasis, dynamically changes expression in response to tissue stresses, and promotes quiescence, whereas Igfbp5 preferentially labels a rapidly-growing, differentiation-prone population. These studies establish the oral mucosa as an important model system to study epithelial stem cell populations and how they respond to tissue stresses.

Many immune diseases occur at different rates among people with schizophrenia compared to the general population. Here, we evaluated whether this phenomenon might be explained by shared genetic risk factors. We used data from large genome-wide association studies to compare the genetic architecture of schizophrenia to 19 immune diseases. First, we evaluated the association with schizophrenia of 581 variants previously reported to be associated with immune diseases at genome-wide significance. We identified five variants with potentially pleiotropic effects. While colocalization analyses were inconclusive, functional characterization of these variants provided the strongest evidence for a model in which genetic variation at rs1734907 modulates risk of schizophrenia and Crohn's disease via altered methylation and expression of EPHB4-a gene whose protein product guides the migration of neuronal axons in the brain and the migration of lymphocytes towards infected cells in the immune system. Next, we investigated genome-wide sharing of common variants between schizophrenia and immune diseases using cross-trait LD score regression. Of the 11 immune diseases with available genome-wide summary statistics, we observed genetic correlation between six immune diseases and schizophrenia: inflammatory bowel disease (rg = 0.12 ± 0.03, P = 2.49 × 10-4), Crohn's disease (rg = 0.097 ± 0.06, P = 3.27 × 10-3), ulcerative colitis (rg = 0.11 ± 0.04, P = 4.05 × 10-3), primary biliary cirrhosis (rg = 0.13 ± 0.05, P = 3.98 × 10-3), psoriasis (rg = 0.18 ± 0.07, P = 7.78 × 10-3) and systemic lupus erythematosus (rg = 0.13 ± 0.05, P = 3.76 × 10-3). With the exception of ulcerative colitis, the degree and direction of these genetic correlations were consistent with the expected phenotypic correlation based on epidemiological data. Our findings suggest shared genetic risk factors contribute to the epidemiological association of certain immune diseases and schizophrenia.

Deep phenotyping is an emerging trend in precision medicine for genetic disease. The shape of the face is affected in 30-40% of known genetic syndromes. Here, we determine whether syndromes can be diagnosed from 3D images of human faces.

Collective migration is essential for development, wound repair, and cancer metastasis. For most collective systems, "leader cells" determine both the direction and the power of the migration. It has remained unclear, however, how the highly polarized vertebrate epithelium migrates directionally during branching morphogenesis. We show here that, unlike in other systems, front-rear polarity of the mammary epithelium is set up by preferential cell proliferation in the front in response to the FGF10 gradient. This leads to frontal stratification, loss of apicobasal polarity, and leader cell formation. Leader cells are a dynamic population and move faster and more directionally toward the FGF10 signal than do follower cells, partly because of their intraepithelial protrusions toward the signal. Together, our data show that directional migration of the mammary epithelium is a unique multistep process and that, despite sharing remarkable cellular and molecular similarities, vertebrate and invertebrate epithelial branching are fundamentally distinct processes.

Epithelial organoids derived from intestinal tissue, called enteroids, recapitulate many aspects of the organ in vitro and can be used for biological discovery, personalized medicine, and drug development. Here, we interrogated the cell signaling environment within the developing human intestine to identify niche cues that may be important for epithelial development and homeostasis. We identified an EGF family member, EPIREGULIN (EREG), which is robustly expressed in the developing human crypt. Enteroids generated from the developing human intestine grown in standard culture conditions, which contain EGF, are dominated by stem and progenitor cells and feature little differentiation and no spatial organization. Our results demonstrate that EREG can replace EGF in vitro, and EREG leads to spatially resolved enteroids that feature budded and proliferative crypt domains and a differentiated villus-like central lumen. Multiomic (transcriptome plus epigenome) profiling of native crypts, EGF-grown enteroids, and EREG-grown enteroids showed that EGF enteroids have an altered chromatin landscape that is dependent on EGF concentration, downregulate the master intestinal transcription factor CDX2, and ectopically express stomach genes, a phenomenon that is reversible. This is in contrast to EREG-grown enteroids, which remain intestine like in culture. Thus, EREG creates a homeostatic intestinal niche in vitro, enabling interrogation of stem cell function, cellular differentiation, and disease modeling.

Costello syndrome (CS) is a congenital disorder caused by heterozygous activating germline HRAS mutations in the canonical Ras/mitogen-activated protein kinase (Ras/MAPK) pathway. CS is one of the RASopathies, a large group of syndromes caused by mutations within various components of the Ras/MAPK pathway. An important part of the phenotype that greatly impacts quality of life is hypotonia. To gain a better understanding of the mechanisms underlying hypotonia in CS, a mouse model with an activating HrasG12V allele was utilized. We identified a skeletal myopathy that was due, in part, to inhibition of embryonic myogenesis and myofiber formation, resulting in a reduction in myofiber size and number that led to reduced muscle mass and strength. In addition to hyperactivation of the Ras/MAPK and PI3K/AKT pathways, there was a significant reduction in p38 signaling, as well as global transcriptional alterations consistent with the myopathic phenotype. Inhibition of Ras/MAPK pathway signaling using a MEK inhibitor rescued the HrasG12V myopathy phenotype both in vitro and in vivo, demonstrating that increased MAPK signaling is the main cause of the muscle phenotype in CS.

Niche-derived factors regulate tissue stem cells, but apart from the mechanosensory pathways, the effect of niche geometry is not well understood. We used organoids and bioengineered tissue culture platforms to demonstrate that the conical shape of Lgr5+ small intestinal stem cells (ISCs) facilitate their self-renewal and function. Inhibition of non-muscle myosin II (NM II)-driven apical constriction altered ISC shape and reduced niche curvature and stem cell capacity. Niche curvature is decreased in aged mice, suggesting that suboptimal interactions between old ISCs and their niche develop with age. We show that activation of NM IIC or physical restriction to young topology improves in vitro regeneration by old epithelium. We propose that the increase in lateral surface area of ISCs induced by apical constriction promotes interactions between neighboring cells, and the curved topology of the intestinal niche has evolved to maximize signaling between ISCs and neighboring cells.

When evolution leads to differences in body size, organs generally scale along. A well-known example of the tight relationship between organ and body size is the scaling of mammalian molar teeth. To investigate how teeth scale during development and evolution, we compared molar development from initiation through final size in the mouse and the rat. Whereas the linear dimensions of the rat molars are twice that of the mouse molars, their shapes are largely the same. Here, we focus on the first lower molars that are considered the most reliable dental proxy for size-related patterns due to their low within-species variability. We found that scaling of the molars starts early, and that the rat molar is patterned equally as fast but in a larger size than the mouse molar. Using transcriptomics, we discovered that a known regulator of body size, insulin-like growth factor 1 (Igf1), is more highly expressed in the rat molars compared to the mouse molars. Ex vivo and in vivo mouse models demonstrated that modulation of the IGF pathway reproduces several aspects of the observed scaling process. Furthermore, analysis of IGF1-treated mouse molars and computational modeling indicate that IGF signaling scales teeth by simultaneously enhancing growth and by inhibiting the cusp-patterning program, thereby providing a relatively simple mechanism for scaling teeth during development and evolution. Finally, comparative data from shrews to elephants suggest that this scaling mechanism regulates the minimum tooth size possible, as well as the patterning potential of large teeth.

The gastrointestinal tract relies on the production, maturation, and transit of mucin to protect against pathogens and to lubricate the epithelial lining. Although the molecular and cellular mechanisms that regulate mucin production and movement are beginning to be understood, the upstream epithelial signals that contribute to mucin regulation remain unclear. Here, we report that the inflammatory cytokine tumor necrosis factor (TNF), generated by the epithelium, contributes to mucin homeostasis by regulating both cell differentiation and cystic fibrosis transmembrane conductance regulator (CFTR) activity. We used genetic mouse models and noninflamed samples from patients with inflammatory bowel disease (IBD) undergoing anti-TNF therapy to assess the effect of in vivo perturbation of TNF. We found that inhibition of epithelial TNF promotes the differentiation of secretory progenitor cells into mucus-producing goblet cells. Furthermore, TNF treatment and CFTR inhibition in intestinal organoids demonstrated that TNF promotes ion transport and luminal flow via CFTR. The absence of TNF led to slower gut transit times, which we propose results from increased mucus accumulation coupled with decreased luminal fluid pumping. These findings point to a TNF/CFTR signaling axis in the adult intestine and identify epithelial cell-derived TNF as an upstream regulator of mucin homeostasis.

Deciphering how signaling pathways interact during development is necessary for understanding the etiopathogenesis of congenital malformations and disease. In several embryonic structures, components of the Hedgehog and retinoic acid pathways, two potent players in development and disease are expressed and operate in the same or adjacent tissues and cells. Yet whether and, if so, how these pathways interact during organogenesis is, to a large extent, unclear. Using genetic and experimental approaches in the mouse, we show that during development of ontogenetically different organs, including the tail, genital tubercle, and secondary palate, Sonic hedgehog (SHH) loss-of-function causes anomalies phenocopying those induced by enhanced retinoic acid signaling and that SHH is required to prevent supraphysiological activation of retinoic signaling through maintenance and reinforcement of expression of the Cyp26 genes. Furthermore, in other tissues and organs, disruptions of the Hedgehog or the retinoic acid pathways during development generate similar phenotypes. These findings reveal that rigidly calibrated Hedgehog and retinoic acid activities are required for normal organogenesis and tissue patterning.

The oral mucosa is one of the most rapidly dividing tissues in the body and serves as a barrier to physical and chemical insults from mastication, food, and microorganisms. Breakdown of this barrier can lead to significant morbidity and potentially life-threatening infections for patients. Determining the identity and organization of oral epithelial progenitor cells (OEPCs) is therefore paramount to understanding their roles in homeostasis and disease. Using lineage tracing and label retention experiments, we show that rapidly dividing OEPCs are located broadly within the basal layer of the mucosa throughout the oral cavity. Quantitative clonal analysis demonstrated that OEPCs undergo population-asymmetrical divisions following neutral drift dynamics and that they respond to chemotherapy-induced damage by altering daughter cell fates. Finally, using single-cell RNA-seq, we establish the basal layer population structure and propose a model that defines the organization of cells within the basal layer.

Epithelial dysfunction and crypt destruction are defining features of inflammatory bowel disease (IBD). However, current IBD therapies targeting epithelial dysfunction are lacking. The nuclear receptor LRH-1 (NR5A2) is expressed in intestinal epithelium and thought to contribute to epithelial renewal. Here we show that LRH-1 maintains intestinal epithelial health and protects against inflammatory damage. Knocking out LRH-1 in murine intestinal organoids reduces Notch signaling, increases crypt cell death, distorts the cellular composition of the epithelium, and weakens the epithelial barrier. Human LRH-1 (hLRH-1) rescues epithelial integrity and when overexpressed, mitigates inflammatory damage in murine and human intestinal organoids, including those derived from IBD patients. Finally, hLRH-1 greatly reduces disease severity in T-cell-mediated murine colitis. Together with the failure of a ligand-incompetent hLRH-1 mutant to protect against TNFα-damage, these findings provide compelling evidence that hLRH-1 mediates epithelial homeostasis and is an attractive target for intestinal disease.