The pan-ErbB negative regulator Lrig1 is an intestinal stem cell marker that functions as a tumor suppressor.

Lineage mapping has identified both proliferative and quiescent intestinal stem cells, but the molecular circuitry controlling stem cell quiescence is incompletely understood. By lineage mapping, we show Lrig1, a pan-ErbB inhibitor, marks predominately noncycling, long-lived stem cells that are located at the crypt base and that, upon injury, proliferate and divide to replenish damaged crypts. Transcriptome profiling of Lrig1(+) colonic stem cells differs markedly from the profiling of highly proliferative, Lgr5(+) colonic stem cells; genes upregulated in the Lrig1(+) population include those involved in cell cycle repression and response to oxidative damage. Loss of Apc in Lrig1(+) cells leads to intestinal adenomas, and genetic ablation of Lrig1 results in heightened ErbB1-3 expression and duodenal adenomas. These results shed light on the relationship between proliferative and quiescent intestinal stem cells and support a model in which intestinal stem cell quiescence is maintained by calibrated ErbB signaling with loss of a negative regulator predisposing to neoplasia.

ArCH: improving the performance of clonal hematopoiesis variant calling and interpretation.

MOTIVATION: The acquisition of somatic mutations in hematopoietic stem and progenitor stem cells with resultant clonal expansion, termed clonal hematopoiesis (CH), is associated with increased risk of hematologic malignancies and other adverse outcomes. CH is generally present at low allelic fractions, but clonal expansion and acquisition of additional mutations leads to hematologic cancers in a small proportion of individuals. With high depth and high sensitivity sequencing, CH can be detected in most adults and its clonal trajectory mapped over time. However, accurate CH variant calling is challenging due to the difficulty in distinguishing low frequency CH mutations from sequencing artifacts. The lack of well-validated bioinformatic pipelines for CH calling may contribute to lack of reproducibility in studies of CH. RESULTS: Here, we developed ArCH, an Artifact filtering Clonal Hematopoiesis variant calling pipeline for detecting single nucleotide variants and short insertions/deletions by combining the output of four variant calling tools and filtering based on variant characteristics and sequencing error rate estimation. ArCH is an end-to-end cloud-based pipeline optimized to accept a variety of inputs with customizable parameters adaptable to multiple sequencing technologies, research questions, and datasets. Using deep targeted sequencing data generated from six acute myeloid leukemia patient tumor: normal dilutions, 31 blood samples with orthogonal validation, and 26 blood samples with technical replicates, we show that ArCH improves the sensitivity and positive predictive value of CH variant detection at low allele frequencies compared to standard application of commonly used variant calling approaches. AVAILABILITY AND IMPLEMENTATION: The code for this workflow is available at: https://github.com/kbolton-lab/ArCH.

Betacoronavirus-specific alternate splicing.

Viruses can subvert a number of cellular processes including splicing in order to block innate antiviral responses, and many viruses interact with cellular splicing machinery. SARS-CoV-2 infection was shown to suppress global mRNA splicing, and at least 10 SARS-CoV-2 proteins bind specifically to one or more human RNAs. Here, we investigate 17 published experimental and clinical datasets related to SARS-CoV-2 infection, datasets from the betacoronaviruses SARS-CoV and MERS, as well as Streptococcus pneumonia, HCV, Zika virus, Dengue virus, influenza H3N2, and RSV. We show that genes showing differential alternative splicing in SARS-CoV-2 have a similar functional profile to those of SARS-CoV and MERS and affect a diverse set of genes and biological functions, including many closely related to virus biology. Additionally, the differentially spliced transcripts of cells infected by coronaviruses were more likely to undergo intron-retention, contain a pseudouridine modification, and have a smaller number of exons as compared with differentially spliced transcripts in the control groups. Viral load in clinical COVID-19 samples was correlated with isoform distribution of differentially spliced genes. A significantly higher number of ribosomal genes are affected by differential alternative splicing and gene expression in betacoronavirus samples, and the betacoronavirus differentially spliced genes are depleted for binding sites of RNA-binding proteins. Our results demonstrate characteristic patterns of differential splicing in cells infected by SARS-CoV-2, SARS-CoV, and MERS. The alternative splicing changes observed in betacoronaviruses infection potentially modify a broad range of cellular functions, via changes in the functions of the products of a diverse set of genes involved in different biological processes.

Investigating rare pathogenic/likely pathogenic exonic variation in bipolar disorder.

Bipolar disorder (BD) is a serious mental illness with substantial common variant heritability. However, the role of rare coding variation in BD is not well established. We examined the protein-coding (exonic) sequences of 3,987 unrelated individuals with BD and 5,322 controls of predominantly European ancestry across four cohorts from the Bipolar Sequencing Consortium (BSC). We assessed the burden of rare, protein-altering, single nucleotide variants classified as pathogenic or likely pathogenic (P-LP) both exome-wide and within several groups of genes with phenotypic or biologic plausibility in BD. While we observed an increased burden of rare coding P-LP variants within 165 genes identified as BD GWAS regions in 3,987 BD cases (meta-analysis OR = 1.9, 95% CI = 1.3-2.8, one-sided p = 6.0 × 10-4), this enrichment did not replicate in an additional 9,929 BD cases and 14,018 controls (OR = 0.9, one-side p = 0.70). Although BD shares common variant heritability with schizophrenia, in the BSC sample we did not observe a significant enrichment of P-LP variants in SCZ GWAS genes, in two classes of neuronal synaptic genes (RBFOX2 and FMRP) associated with SCZ or in loss-of-function intolerant genes. In this study, the largest analysis of exonic variation in BD, individuals with BD do not carry a replicable enrichment of rare P-LP variants across the exome or in any of several groups of genes with biologic plausibility. Moreover, despite a strong shared susceptibility between BD and SCZ through common genetic variation, we do not observe an association between BD risk and rare P-LP coding variants in genes known to modulate risk for SCZ.

The whole-genome landscape of Burkitt lymphoma subtypes.

Burkitt lymphoma (BL) is an aggressive, MYC-driven lymphoma comprising 3 distinct clinical subtypes: sporadic BLs that occur worldwide, endemic BLs that occur predominantly in sub-Saharan Africa, and immunodeficiency-associated BLs that occur primarily in the setting of HIV. In this study, we comprehensively delineated the genomic basis of BL through whole-genome sequencing (WGS) of 101 tumors representing all 3 subtypes of BL to identify 72 driver genes. These data were additionally informed by CRISPR screens in BL cell lines to functionally annotate the role of oncogenic drivers. Nearly every driver gene was found to have both coding and non-coding mutations, highlighting the importance of WGS for identifying driver events. Our data implicate coding and non-coding mutations in IGLL5, BACH2, SIN3A, and DNMT1. Epstein-Barr virus (EBV) infection was associated with higher mutation load, with type 1 EBV showing a higher mutational burden than type 2 EBV. Although sporadic and immunodeficiency-associated BLs had similar genetic profiles, endemic BLs manifested more frequent mutations in BCL7A and BCL6 and fewer genetic alterations in DNMT1, SNTB2, and CTCF. Silencing mutations in ID3 were a common feature of all 3 subtypes of BL. In vitro, mass spectrometry-based proteomics demonstrated that the ID3 protein binds primarily to TCF3 and TCF4. In vivo knockout of ID3 potentiated the effects of MYC, leading to rapid tumorigenesis and tumor phenotypes consistent with those observed in the human disease.

Reversible deficits in apical transporter trafficking associated with deficiency in diacylglycerol acyltransferase.

Deficiency in diacylglycerol acyltransferase (DGAT1) is a rare cause of neonatal diarrhea, without a known mechanism or in vitro model. A patient presenting at our institution at 7 weeks of life with failure to thrive and diarrhea was found by whole-exome sequencing to have a homozygous DGAT1 truncation mutation. Duodenal biopsies showed loss of DGAT1 and deficits in apical membrane transporters and junctional proteins in enterocytes. When placed on a very low-fat diet, the patient's diarrhea resolved with normalization of brush border transporter localization in endoscopic biopsies. DGAT1 knockdown in Caco2-BBe cells modeled the deficits in apical trafficking, with loss of apical DPPIV and junctional occludin. Elevation in cellular lipid levels, including diacylglycerol (DAG) and phospholipid metabolites of DAG, was documented by lipid analysis in DGAT1 knockdown cells. Culture of the DGAT1 knockdown cells in lipid-depleted media led to re-establishment of occludin and return of apical DPPIV. DGAT1 loss appears to elicit global changes in enterocyte polarized trafficking that could account for deficits in absorption seen in the patient. The in vitro modeling of this disease should allow for investigation of possible therapeutic targets.

Functional variants in TBX2 are associated with a syndromic cardiovascular and skeletal developmental disorder.

The 17 genes of the T-box family are transcriptional regulators that are involved in all stages of embryonic development, including craniofacial, brain, heart, skeleton and immune system. Malformation syndromes have been linked to many of the T-box genes. For example, haploinsufficiency of TBX1 is responsible for many structural malformations in DiGeorge syndrome caused by a chromosome 22q11.2 deletion. We report four individuals with an overlapping spectrum of craniofacial dysmorphisms, cardiac anomalies, skeletal malformations, immune deficiency, endocrine abnormalities and developmental impairments, reminiscent of DiGeorge syndrome, who are heterozygotes for TBX2 variants. The p.R20Q variant is shared by three affected family members in an autosomal dominant manner; the fourth unrelated individual has a de novo p.R305H mutation. Bioinformatics analyses indicate that these variants are rare and predict them to be damaging. In vitro transcriptional assays in cultured cells show that both variants result in reduced transcriptional repressor activity of TBX2. We also show that the variants result in reduced protein levels of TBX2. Heterologous over-expression studies in Drosophila demonstrate that both p.R20Q and p.R305H function as partial loss-of-function alleles. Hence, these and other data suggest that TBX2 is a novel candidate gene for a new multisystem malformation disorder.

Estimated number of adult survivors of childhood cancer in United States with cancer-predisposing germline variants.

PURPOSE: To estimate the absolute number of adult survivors of childhood cancer in the U.S. population who carry a pathogenic or likely pathogenic variant in a cancer predisposition gene. METHODS: Using the Surveillance, Epidemiology, and End Results (SEER) Program, we estimated the number of childhood cancer survivors on December 31, 2016 for each childhood cancer diagnosis, multiplied this by the proportion of carriers of pathogenic/likely pathogenic variants in the St. Jude Lifetime Cohort (SJLIFE) study, and projected the resulting number onto the U.S. RESULTS: Based on genome sequence data, 11.8% of 2450 SJLIFE participants carry a pathogenic/likely pathogenic variant in one of 156 cancer predisposition genes. Given this information, we estimate that 21 800 adult survivors of childhood cancer in the United States carry a pathogenic/likely pathogenic variant in one of these genes. The highest estimated absolute number of variant carriers are among survivors of central nervous system tumors (n = 4300), particularly astrocytoma (n = 1800) and other gliomas (n = 1700), acute lymphoblastic leukemia (n = 4300), and retinoblastoma (n = 3500). The most frequently mutated genes are RB1 (n = 3000), NF1 (n = 2300), and BRCA2 (n = 800). CONCLUSION: Given the increasing number of childhood cancer survivors in the United States, clinicians should counsel survivors regarding their potential genetic risk, consider referral for genetic counseling and testing, and, as appropriate, implement syndrome-specific cancer surveillance or risk-reducing measures.

AS03-Adjuvanted H5N1 Avian Influenza Vaccine Modulates Early Innate Immune Signatures in Human Peripheral Blood Mononuclear Cells.

BACKGROUND: Adjuvant System 03 (AS03) markedly enhances responses to influenza A/H5N1 vaccines, but the mechanisms of this enhancement are incompletely understood. METHODS: Using ribonucleic acid sequencing on peripheral blood mononuclear cells (PBMCs) from AS03-adjuvanted and unadjuvanted inactivated H5N1 vaccine recipients, we identified differentially expressed genes, enriched pathways, and genes that correlated with serologic responses. We compared bulk PBMC findings with our previously published assessments of flow-sorted immune cell types. RESULTS: AS03-adjuvanted vaccine induced the strongest differential signals on day 1 postvaccination, activating multiple innate immune pathways including interferon and JAK-STAT signaling, Fcγ receptor (FcγR)-mediated phagocytosis, and antigen processing and presentation. Changes in signal transduction and immunoglobulin genes predicted peak hemagglutinin inhibition (HAI) titers. Compared with individual immune cell types, activated PBMC genes and pathways were most similar to innate immune cells. However, several pathways were unique to PBMCs, and several pathways identified in individual cell types were absent in PBMCs. CONCLUSIONS: Transcriptomic analysis of PBMCs after AS03-adjuvanted H5N1 vaccination revealed early activation of innate immune signaling, including a 5- to 8-fold upregulation of FcγR1A/1B/1C genes. Several early gene responses were correlated with HAI titer, indicating links with the adaptive immune response. Although PBMCs and cell-specific results shared key innate immune signals, unique signals were identified by both approaches.

O-GlcNAc transferase missense mutations linked to X-linked intellectual disability deregulate genes involved in cell fate determination and signaling.

It is estimated that ∼1% of the world's population has intellectual disability, with males affected more often than females. OGT is an X-linked gene encoding for the enzyme O-GlcNAc transferase (OGT), which carries out the reversible addition of N-acetylglucosamine (GlcNAc) to Ser/Thr residues of its intracellular substrates. Three missense mutations in the tetratricopeptide (TPR) repeats of OGT have recently been reported to cause X-linked intellectual disability (XLID). Here, we report the discovery of two additional novel missense mutations (c.775 G>A, p.A259T, and c.1016 A>G, p.E339G) in the TPR domain of OGT that segregate with XLID in affected families. Characterization of all five of these XLID missense variants of OGT demonstrates modest declines in thermodynamic stability and/or activities of the variants. We engineered each of the mutations into a male human embryonic stem cell line using CRISPR/Cas9. Investigation of the global O-GlcNAc profile as well as OGT and O-GlcNAc hydrolase levels by Western blotting showed no gross changes in steady-state levels in the engineered lines. However, analyses of the differential transcriptomes of the OGT variant-expressing stem cells revealed shared deregulation of genes involved in cell fate determination and liver X receptor/retinoid X receptor signaling, which has been implicated in neuronal development. Thus, here we reveal two additional mutations encoding residues in the TPR regions of OGT that appear causal for XLID and provide evidence that the relatively stable and active TPR variants may share a common, unelucidated mechanism of altering gene expression profiles in human embryonic stem cells.

Advancements in Next-Generation Sequencing.

The term next-generation sequencing is almost a decade old, but it remains the colloquial way to describe highly parallel or high-output sequencing methods that produce data at or beyond the genome scale. Since the introduction of these technologies, the number of applications and methods that leverage the power of genome-scale sequencing has increased at an exponential pace. This review highlights recent concepts, technologies, and methods from next-generation sequencing to illustrate the breadth and depth of the applications and research areas that are driving progress in genomics.

The Airplane Cabin Microbiome.

Serving over three billion passengers annually, air travel serves as a conduit for infectious disease spread, including emerging infections and pandemics. Over two dozen cases of in-flight transmissions have been documented. To understand these risks, a characterization of the airplane cabin microbiome is necessary. Our study team collected 229 environmental samples on ten transcontinental US flights with subsequent 16S rRNA sequencing. We found that bacterial communities were largely derived from human skin and oral commensals, as well as environmental generalist bacteria. We identified clear signatures for air versus touch surface microbiome, but not for individual types of touch surfaces. We also found large flight-to-flight beta diversity variations with no distinguishing signatures of individual flights, rather a high between-flight diversity for all touch surfaces and particularly for air samples. There was no systematic pattern of microbial community change from pre- to post-flight. Our findings are similar to those of other recent studies of the microbiome of built environments. In summary, the airplane cabin microbiome has immense airplane to airplane variability. The vast majority of airplane-associated microbes are human commensals or non-pathogenic, and the results provide a baseline for non-crisis-level airplane microbiome conditions.

Comparative analysis of chimeric ZFP-, TALE- and Cas9-piggyBac transposases for integration into a single locus in human cells.

Integrating DNA delivery systems hold promise for many applications including treatment of diseases; however, targeted integration is needed for improved safety. The piggyBac (PB) transposon system is a highly active non-viral gene delivery system capable of integrating defined DNA segments into host chromosomes without requiring homologous recombination. We systematically compared four different engineered zinc finger proteins (ZFP), four transcription activator-like effector proteins (TALE), CRISPR associated protein 9 (SpCas9) and the catalytically inactive dSpCas9 protein fused to the amino-terminus of the transposase enzyme designed to target the hypoxanthine phosphoribosyltransferase (HPRT) gene located on human chromosome X. Chimeric transposases were evaluated for expression, transposition activity, chromatin immunoprecipitation at the target loci, and targeted knockout of the HPRT gene in human cells. One ZFP-PB and one TALE-PB chimera demonstrated notable HPRT gene targeting. In contrast, Cas9/dCas9-PB chimeras did not result in gene targeting. Instead, the HPRT locus appeared to be protected from transposon integration. Supplied separately, PB permitted highly efficient isolation of Cas9-mediated knockout of HPRT, with zero transposon integrations in HPRT by deep sequencing. In summary, these tools may allow isolation of 'targeted-only' cells, be utilized to protect a genomic locus from transposon integration, and enrich for Cas9-mutated cells.

A signalling cascade of IL-33 to IL-13 regulates metaplasia in the mouse stomach.

OBJECTIVE: Alternatively activated macrophages (M2) are associated with the progression of spasmolytic polypeptide-expressing metaplasia (SPEM) in the stomach. However, the precise mechanism(s) and critical mediators that induce SPEM are unknown. DESIGN: To determine candidate genes important in these processes, macrophages from the stomach corpus of mice with SPEM (DMP-777-treated) or advanced SPEM (L635-treated) were isolated and RNA sequenced. Effects on metaplasia development after acute parietal cell loss induced by L635 were evaluated in interleukin (IL)-33, IL-33 receptor (ST2) and IL-13 knockout (KO) mice. RESULTS: Profiling of metaplasia-associated macrophages in the stomach identified an M2a-polarised macrophage population. Expression of IL-33 was significantly upregulated in macrophages associated with advanced SPEM. L635 induced metaplasia in the stomachs of wild-type mice, but not in the stomachs of IL-33 and ST2 KO mice. While IL-5 and IL-9 were not required for metaplasia induction, IL-13 KO mice did not develop metaplasia in response to L635. Administration of IL-13 to ST2 KO mice re-established the induction of metaplasia following acute parietal cell loss. CONCLUSIONS: Metaplasia induction and macrophage polarisation after parietal cell loss is coordinated through a cytokine signalling network of IL-33 and IL-13, linking a combined response to injury by both intrinsic mucosal mechanisms and infiltrating M2 macrophages.

Defects of CRB2 cause steroid-resistant nephrotic syndrome.

Nephrotic syndrome (NS), the association of gross proteinuria, hypoalbuminaemia, edema, and hyperlipidemia, can be clinically divided into steroid-sensitive (SSNS) and steroid-resistant (SRNS) forms. SRNS regularly progresses to end-stage renal failure. By homozygosity mapping and whole exome sequencing, we here identify recessive mutations in Crumbs homolog 2 (CRB2) in four different families affected by SRNS. Previously, we established a requirement for zebrafish crb2b, a conserved regulator of epithelial polarity, in podocyte morphogenesis. By characterization of a loss-of-function mutation in zebrafish crb2b, we now show that zebrafish crb2b is required for podocyte foot process arborization, slit diaphragm formation, and proper nephrin trafficking. Furthermore, by complementation experiments in zebrafish, we demonstrate that CRB2 mutations result in loss of function and therefore constitute causative mutations leading to NS in humans. These results implicate defects in podocyte apico-basal polarity in the pathogenesis of NS.

Genomic sequencing identifies secondary findings in a cohort of parent study participants.

PURPOSE: Clinically relevant secondary variants were identified in parents enrolled with a child with developmental delay and intellectual disability. METHODS: Exome/genome sequencing and analysis of 789 "unaffected" parents was performed. RESULTS: Pathogenic/likely pathogenic variants were identified in 21 genes within 25 individuals (3.2%), with 11 (1.4%) participants harboring variation in a gene defined as clinically actionable by the American College of Medical Genetics and Genomics. These 25 individuals self-reported either relevant clinical diagnoses (5); relevant family history or symptoms (13); or no relevant family history, symptoms, or clinical diagnoses (7). A limited carrier screen was performed yielding 15 variants in 48 (6.1%) parents. Parents were also analyzed as mate pairs (n = 365) to identify cases in which both parents were carriers for the same recessive disease, yielding three such cases (0.8%), two of which had children with the relevant recessive disease. Four participants had two findings (one carrier and one noncarrier variant). In total, 71 of the 789 enrolled parents (9.0%) received secondary findings. CONCLUSION: We provide an overview of the rates and types of clinically relevant secondary findings, which may be useful in the design and implementation of research and clinical sequencing efforts to identify such findings.

Patterns and rates of exonic de novo mutations in autism spectrum disorders.

Autism spectrum disorders (ASD) are believed to have genetic and environmental origins, yet in only a modest fraction of individuals can specific causes be identified. To identify further genetic risk factors, here we assess the role of de novo mutations in ASD by sequencing the exomes of ASD cases and their parents (n = 175 trios). Fewer than half of the cases (46.3%) carry a missense or nonsense de novo variant, and the overall rate of mutation is only modestly higher than the expected rate. In contrast, the proteins encoded by genes that harboured de novo missense or nonsense mutations showed a higher degree of connectivity among themselves and to previous ASD genes as indexed by protein-protein interaction screens. The small increase in the rate of de novo events, when taken together with the protein interaction results, are consistent with an important but limited role for de novo point mutations in ASD, similar to that documented for de novo copy number variants. Genetic models incorporating these data indicate that most of the observed de novo events are unconnected to ASD; those that do confer risk are distributed across many genes and are incompletely penetrant (that is, not necessarily sufficient for disease). Our results support polygenic models in which spontaneous coding mutations in any of a large number of genes increases risk by 5- to 20-fold. Despite the challenge posed by such models, results from de novo events and a large parallel case-control study provide strong evidence in favour of CHD8 and KATNAL2 as genuine autism risk factors.

Targeted next-generation sequencing reveals MODY in up to 6.5% of antibody-negative diabetes cases listed in the Norwegian Childhood Diabetes Registry.

AIMS/HYPOTHESIS: MODY can be wrongly diagnosed as type 1 diabetes in children. We aimed to find the prevalence of MODY in a nationwide population-based registry of childhood diabetes. METHODS: Using next-generation sequencing, we screened the HNF1A, HNF4A, HNF1B, GCK and INS genes in all 469 children (12.1%) negative for both GAD and IA-2 autoantibodies and 469 antibody-positive matched controls selected from the Norwegian Childhood Diabetes Registry (3882 children). Variants were classified using clinical diagnostic criteria for pathogenicity ranging from class 1 (neutral) to class 5 (pathogenic). RESULTS: We identified 58 rare exonic and splice variants in cases and controls. Among antibody-negative patients, 6.5% had genetic variants of classes 3-5 (vs 2.4% in controls; p = 0.002). For the stricter classification (classes 4 and 5), the corresponding number was 4.1% (vs 0.2% in controls; p = 1.6 × 10-5). HNF1A showed the strongest enrichment of class 3-5 variants, with 3.9% among antibody-negative patients (vs 0.4% in controls; p = 0.0002). Antibody-negative carriers of variants in class 3 had a similar phenotype to those carrying variants in classes 4 and 5. CONCLUSIONS/INTERPRETATION: This is the first study screening for MODY in all antibody-negative children in a nationwide population-based registry. Our results suggest that the prevalence of MODY in antibody-negative childhood diabetes may reach 6.5%. One-third of these MODY cases had not been recognised by clinicians. Since a precise diagnosis is important for treatment and genetic counselling, molecular screening of all antibody-negative children should be considered in routine diagnostics.

Mutations in EMP2 cause childhood-onset nephrotic syndrome.

Nephrotic syndrome (NS) is a genetically heterogeneous group of diseases that are divided into steroid-sensitive NS (SSNS) and steroid-resistant NS (SRNS). SRNS inevitably leads to end-stage kidney disease, and no curative treatment is available. To date, mutations in more than 24 genes have been described in Mendelian forms of SRNS; however, no Mendelian form of SSNS has been described. To identify a genetic form of SSNS, we performed homozygosity mapping, whole-exome sequencing, and multiplex PCR followed by next-generation sequencing. We thereby detected biallelic mutations in EMP2 (epithelial membrane protein 2) in four individuals from three unrelated families affected by SRNS or SSNS. We showed that EMP2 exclusively localized to glomeruli in the kidney. Knockdown of emp2 in zebrafish resulted in pericardial effusion, supporting the pathogenic role of mutated EMP2 in human NS. At the cellular level, we showed that knockdown of EMP2 in podocytes and endothelial cells resulted in an increased amount of CAVEOLIN-1 and decreased cell proliferation. Our data therefore identify EMP2 mutations as causing a recessive Mendelian form of SSNS.

Vascular endothelial growth factor coordinates islet innervation via vascular scaffolding.

Neurovascular alignment is a common anatomical feature of organs, but the mechanisms leading to this arrangement are incompletely understood. Here, we show that vascular endothelial growth factor (VEGF) signaling profoundly affects both vascularization and innervation of the pancreatic islet. In mature islets, nerves are closely associated with capillaries, but the islet vascularization process during embryonic organogenesis significantly precedes islet innervation. Although a simple neuronal meshwork interconnects the developing islet clusters as they begin to form at E14.5, the substantial ingrowth of nerve fibers into islets occurs postnatally, when islet vascularization is already complete. Using genetic mouse models, we demonstrate that VEGF regulates islet innervation indirectly through its effects on intra-islet endothelial cells. Our data indicate that formation of a VEGF-directed, intra-islet vascular plexus is required for development of islet innervation, and that VEGF-induced islet hypervascularization leads to increased nerve fiber ingrowth. Transcriptome analysis of hypervascularized islets revealed an increased expression of extracellular matrix components and axon guidance molecules, with these transcripts being enriched in the islet-derived endothelial cell population. We propose a mechanism for coordinated neurovascular development within pancreatic islets, in which endocrine cell-derived VEGF directs the patterning of intra-islet capillaries during embryogenesis, forming a scaffold for the postnatal ingrowth of essential autonomic nerve fibers.