Pancreatic progenitor epigenome maps prioritize type 2 diabetes risk genes with roles in development

  1. Ryan J Geusz
  2. Allen Wang
  3. Joshua Chiou
  4. Joseph J Lancman
  5. Nichole Wetton
  6. Samy Kefalopoulou
  7. Jinzhao Wang
  8. Yunjiang Qui
  9. Jian Yan
  10. Anthony Aylward
  11. Bing Ren
  12. P Duc Si Dong
  13. Kyle J Gaulton  Is a corresponding author
  14. Maike Sander  Is a corresponding author
  1. Department of Pediatrics, Pediatric Diabetes Research Center, University of California, San Diego, United States
  2. Department of Cellular & Molecular Medicine, University of California, San Diego, United States
  3. Sanford Consortium for Regenerative Medicine, United States
  4. Biomedical Graduate Studies Program, University of California, San Diego, United States
  5. Human Genetics Program, Sanford Burnham Prebys Medical Discovery Institute, United States
  6. Graduate School of Biomedical Sciences, Sanford Burnham Prebys Medical Discovery Institute, United States
  7. Ludwig Institute for Cancer Research, United States
6 figures, 1 table and 1 additional file

Figures

Figure 1 with 1 supplement
T2D-associated risk variants are enriched in stretch enhancers of pancreatic progenitors independent of islet stretch enhancers.

(A) Schematic illustrating the stepwise differentiation of human embryonic stem cells (hES) into pancreatic progenitors (solid arrows) and lineage relationship to islets (dotted arrow). …

Figure 1—figure supplement 1
Characterization of typical and stretch enhancers in pancreatic developmental intermediates and islets.

(A) Diagram illustrating incorporation of histone modification and CTCF ChIP-seq data to generate chromatin state calls via ChromHMM. (B) Percentage of the genome covered by defined chromatin states …

Figure 2 with 1 supplement
Candidate target genes of pancreatic progenitor-specific stretch enhancers regulate developmental processes.

(A) Schematic illustrating identification of pancreatic progenitor-specific stretch enhancers (PSSE). (B) Heatmap showing density of H3K27ac ChIP-seq and ATAC-seq reads at PSSE, centered on …

Figure 2—source data 1

Chromosomal coordinates of pancreatic progenitor-specific stretch enhancers (PSSE).

https://cdn.elifesciences.org/articles/59067/elife-59067-fig2-data1-v1.xlsx
Figure 2—source data 2

Enriched gene ontology terms for PSSE-associated genes.

https://cdn.elifesciences.org/articles/59067/elife-59067-fig2-data2-v1.xlsx
Figure 2—source data 3

Proportion of variants nominally associated with beta cell functional traits.

https://cdn.elifesciences.org/articles/59067/elife-59067-fig2-data3-v1.xlsx
Figure 2—source data 4

Tissue identity of downloaded data from ROADMAP consortium.

https://cdn.elifesciences.org/articles/59067/elife-59067-fig2-data4-v1.xlsx
Figure 2—figure supplement 1
Characterization of pancreatic progenitor-specific stretch enhancers.

(A) Percentage of PSSE overlapping with at least one ATAC-seq peak at PP2. Enrichment analysis comparing observed and expected overlap based on random genomic regions of the same size and located on …

Figure 3 with 1 supplement
Identification of T2D risk variants associated with pancreatic progenitor-specific stretch enhancers.

(A) Manhattan plot showing T2D association p-values (from Mahajan et al., 2018) for 10,738 variants mapping within PSSE. The dotted line shows the threshold for Bonferroni correction (p=4.66 × 10−6).…

Figure 3—figure supplement 1
Activity of T2D risk-associated pancreatic progenitor-specific stretch enhancers across human tissues.

(A) mRNA levels (measured in fragments per kilobase per million fragments mapped [FPKM]) at PP2 and in human embryonic pancreas (54 and 58 days gestation, Emb Pan) of all genes expressed (FPKM ≥1) …

Figure 4 with 2 supplements
A T2D risk-associated LAMA1 pancreatic progenitor-specific stretch enhancer regulates LAMA1 expression specifically in pancreatic progenitors.

(A) (Top) Locus plots showing T2D association p-values for variants in a 35 kb window (hg19 chr18:7,050,000–7,085,000) at the LAMA1 locus and LAMA1 PSSE (red box). Fine mapped variants within the …

Figure 4—source data 1

Genes downregulated in ∆LAMA1Enh PP2 stage cells compared to control cells (p adj. <0.05).

https://cdn.elifesciences.org/articles/59067/elife-59067-fig4-data1-v1.xlsx
Figure 4—source data 2

Genes upregulated in ∆LAMA1Enh PP2 stage cells compared to control cells (p adj. <0.05).

https://cdn.elifesciences.org/articles/59067/elife-59067-fig4-data2-v1.xlsx
Figure 4—figure supplement 1
Deletion of the LAMA1-associated pancreatic progenitor-specific enhancer does not affect pancreatic lineage specification.

(A) Odds ratio estimates (points) and 95% CIs (lines) for rs10502347 association with T2D and metabolic GWAS traits. Significant associations are colored black, non-significant are colored light …

Figure 4—figure supplement 2
Deletion of LAMA1 does not affect pancreatic lineage specification.

(A) Schematic illustrating CRISPR-Cas9-mediated deletion strategy of LAMA1 to generate ∆LAMA1 hESC clonal line. (B) Immunofluorescent staining for Laminin in control and ∆LAMA1 PP2 cells …

Figure 5 with 3 supplements
A T2D risk-associated CRB2 pancreatic progenitor-specific stretch enhancer regulates CRB2 expression specifically in pancreatic progenitors.

(A) (Top) Locus plots showing T2D association p-values for variants in a 35 kb window (hg19 chr9:126,112,000–126,147,000) at the CRB2 locus and CRB2 PSSE (red box). Fine mapped variants within the …

Figure 5—source data 1

Genes downregulated in ∆CRB2Enh PP2 stage cells compared to control cells (p adj. <0.05).

https://cdn.elifesciences.org/articles/59067/elife-59067-fig5-data1-v1.xlsx
Figure 5—figure supplement 1
Activity of CRB2- and PGM1-associated pancreatic progenitor-specific stretch enhancers across human tissues.

(A) (Top) Locus plots showing T2D association p-values for variants in a 43 kb window (hg19 chr1:64,084,000–64,127,000) at the PGM1 locus and PGM1 PSSE (red box). Fine mapped variants within the 99% …

Figure 5—figure supplement 2
Deletion of the CRB2-associated pancreatic progenitor-specific enhancer does not affect pancreatic lineage specification.

(A) Schematic illustrating CRISPR-Cas9-mediated deletion strategy of CRB2-associated PSSE to generate independent ∆CRB2Enh hESC clones with different DNA cleavage products. (B) Flow cytometry …

Figure 5—figure supplement 3
Deletion of CRB2 does not affect pancreatic lineage specification.

(A) Schematic illustrating CRISPR-Cas9-mediated deletion strategy of CRB2 to generate ∆CRB2 hESC clonal line. (B) Immunofluorescent staining for CRB2 in control and ∆CRB2 PP2 cells (representative …

Figure 6 with 4 supplements
lama1 and crb2 regulate pancreas morphogenesis and beta cell differentiation.

(A,B) Representative 3D renderings of Tg(ptf1a:eGFP)jh1 control zebrafish embryos (A,A’) and lama1 morphants (B,B’) stained with DAPI (nuclei, blue) and antibody against insulin (red); n ≥ 15 …

Figure 6—figure supplement 1
Laminin and Crb are expressed in zebrafish pancreas progenitors.

(A,B) Confocal images of 36 hr post fertilization (hpf) Tg(ptf1a:eGFP)jh1 zebrafish foregut endoderm labeled with DAPI (nuclei, grey) and antibodies against Nkx6.1 (blue; pancreatic progenitors) and …

Figure 6—figure supplement 2
Validation of morpholinos targeting lama1.

Confocal images (3D rendered) of 32–34 hpf Tg(ptf1a:eGFP)jh1 zebrafish foregut endoderm labeled with DAPI (nuclei, grey) and anti-Laminin antibodies (red). (A–A’’) Z-focal plane image of pancreatic …

Figure 6—figure supplement 3
Validation of morpholinos targeting crb2a and crb2b.

(A,B) 3D renderings of 32–34 hpf Tg(ptf1a:eGFP)jh1 zebrafish foregut endoderm labeled with DAPI (nuclei, grey) and anti-Crb2a antibodies (red). (A–A’’) Z-focal plane image showing Crb2a labeled …

Figure 6—figure supplement 4
crb2b but not crb2a regulates pancreatic beta cell differentiation.

Quantification of beta (insulin+) cell nuclei per embryo in Tg(ptf1a:eGFP)jh1 control zebrafish embryos and crb2a or crb2b morphants at 78 hr post fertilization (hpf). p adj. = 0.91 and 4.4 × 10−2

Tables

Key resources table
Reagent type (species)
or resource
DesignationSource or referenceIdentifiersAdditional information
AntibodyAPC Mouse monoclonal IgG1, κ Isotype ControlBD PharmingenCat# 555751, RRID:AB_398613Flow cytometry (1:100)
AntibodyChicken polyclonal anti-GFPAves LabsCat# GFP-1020, RRID:AB_10000240Immunohistochemistry (1:200)
AntibodyCy3-conjugated donkey polyclonal anti-mouseJackson ImmunoResearch LabsCat# 715-165-150, RRID:AB_2340813Immunofluorescence (1:1000)
AntibodyDyLight 488-conjugated donkey polyconal anti-goatJackson ImmunoResearch LabsCat# 705-545-003, RRID:AB_2340428Immunofluorescence (1:500)
AntibodyGoat polyclonal anti-CTCFSanta Cruz BiotechnologyCat# SC-15914X, RRID:AB_2086899ChIP-seq (4 ug)
AntibodyGoat polyclonal anti-FOXA1AbcamCat# ab5089, RRID:AB_304744ChIP-seq (4 ug)
AntibodyGoat polyclonal anti-FOXA2Santa Cruz BiotechnologyCat# sc-6554, RRID:AB_2262810ChIP-seq (4 ug)
AntibodyGoat polyclonal anti-GATA4Santa Cruz BiotechnologyCat# sc-1237, RRID:AB_2108747ChIP-seq (4 ug)
AntibodyGoat polyclonal anti-PDX1AbcamCat# ab47383, RRID:AB_2162359Immunofluorescence (1:500)
AntibodyGuinea pig polyclonal anti-InsulinBiomedaCat# v2024Immunohistochemistry (1:200)
AntibodyMouse monoclonal anti-Crb2aZIRCCat# Zs-4Immunohistochemistry (1:100)
AntibodyMouse polyclonal anti-GATA6Santa Cruz BiotechnologyCat# sc-9055, RRID:AB_2108768ChIP-seq (4 ug)
AntibodyMouse monoclonal anti-NKX6.1Developmental Studies Hybridoma BankCat# F64A6B4, RRID:AB_532380Immunofluorescence (1:300)
AntibodyMouse monoclonal anti-NKX6.1-Alexa Fluor 647BD BiosciencesCat# 563338, RRID:AB_2738144Flow cytometry (1:5)
AntibodyMouse monoclonal anti-NKX6.1Developmental Studies Hybridoma BankCat# F55A10, RRID:AB_532378Immunohistochemistry (1:10)
AntibodyMouse monoclonal anti-PDX1-PEBD BiosciencesCat# 562161, RRID:AB_10893589Flow cytometry (1:10)
AntibodyPE Mouse monoclonal IgG1, κ Isotype ControlBD PharmingenCat# 555749, RRID:AB_396091Flow cytometry (1:100)
AntibodyRabbit polyclonal anti-CRB2SigmaCat # SAB1301340Immunofluorescence (1:500)
AntibodyRabbit polyclonal anti-H3K27acActive MotifCat# 39133, RRID:AB_2561016ChIP-seq (4 ug)
AntibodyRabbit polyclonal anti-H3K4me1AbcamCat# ab8895, RRID:AB_306847ChIP-seq (4 ug)
AntibodyRabbit polyclonal anti-HNF6Santa Cruz BiotechnologyCat# sc-13050, RRID:AB_2251852ChIP-seq (4 ug)
AntibodyRabbit polyclonal anti-lamininSigmaCat# L9393, RRID:AB_477163Immunohistochemistry (1:100)
Immunofluorescence (1:30)
AntibodyRabbit monoclonal anti-panCrbJensen Laboratory, University of Massachusetts, AmherstN/AImmunohistochemistry (1:100)
AntibodyRabbit polyclonal anti-PDX1Beta Cell Biology ConsortiumAB1068ChIP-seq (4 ug)
AntibodyRabbit polyclonal anti-SOX9ChemiconCat# 5535, RRID:AB_2239761ChIP-seq (4 ug)
Cell line (Homo-sapiens)CyT49ViaCyte, IncNIHhESC-10–0041,
RRID:CVCL_B850
Male
Cell line (Homo-sapiens)H1WiCell Research
Institute
NIHhESC-10–0043,
RRID:CVCL_9771
Male
Chemical compound, drug2-MercaptoethanolThermo Fisher ScientificCat# 21985023
Chemical compound, drugAccutaseThermo Fisher ScientificCat# 00-4555-56
Chemical compound, drugB-27 supplementThermo Fisher ScientificCat# 17504044
Chemical compound, drugBovine Albumin Fraction VLife TechnologiesCat# 15260037
Chemical compound, drugD-(+)-Glucose Solution, 45%Sigma-AldrichCat# G8769
Chemical compound, drugDAPIInvitrogenCat# D1306Immunohistochemistry (1:200)
Chemical compound, drugDMEM High GlucoseVWRCat# 16750–082
Chemical compound, drugDMEM/F12 [-] L-glutamineVWRCat# 15–090-CV
Chemical compound, drugDMEM/F12 with L-Glutamine, HEPESCorningCat# 45000–350
Chemical compound, drugDMFEMD MilliporeCat# DX1730
Chemical compound, drugDPBSThermo Fisher ScientificCat# 21–031-CV
Chemical compound, drugDTTSigmaCat# D9779
Chemical compound, drugFetal Bovine SerumThermo Fisher ScientificCat# MT35011CV
Chemical compound, drugGlutamaxThermo Fisher ScientificCat# 35050–079
Chemical compound, drugGlutaMAXThermo Fisher ScientificCat# 35050061
Chemical compound, drugHoechst 33342Thermo Fisher ScientificCat# H3570
Chemical compound, drugHyClone Dulbecco’s Modified Eagles MediumThermo Fisher
Scientific
Cat# SH30081.FS
Chemical compound, drugIGEPAL-CA630SigmaCat# I8896
Chemical compound, drugIllumina tagmentation enzymeIlluminaCat# FC-121–1030
Chemical compound, drugInsulin-Transferrin-Selenium (ITS)Thermo Fisher ScientificCat# 41400045
Chemical compound, drugInsulin-Transferrin-Selenium-Ethanolamine (ITS-X)Thermo Fisher
Scientific
Cat# 51500–056
Chemical compound, drugKAAD-CyclopamineToronto Research ChemicalsCat# K171000
Chemical compound, drugK-acetateSigmaCat# P5708
Chemical compound, drugKnockOut SR XenoFreeThermo Fisher ScientificCat# A1099202
Chemical compound, drugLDN-193189StemgentCat# 04–0074
Chemical compound, drugMatrigelCorningCat# 356231
Chemical compound, drugMCDB 131Thermo Fisher
Scientific
Cat# 10372–019
Chemical compound, drugMg-acetateSigmaCat# M2545
Chemical compound, drugmTeSR1 Complete Kit - GMPSTEMCELL
Technologies
Cat# 85850
Chemical compound, drugNEBNext High-Fidelity 2X PCR Master MixNEBCat# M0541
Chemical compound, drugNon-Essential Amino AcidsThermo Fisher ScientificCat# 11140050
Chemical compound, drugO.C.T. CompoundSakura Finetek USACat# 25608–930
Chemical compound, drugPenicillin-StreptomycinThermo Fisher ScientificCat# 15140122
Chemical compound, drugPolyethylenimine (PEI)PolysciencesCat# 23966–1
Chemical compound, drugProtease inhibitorRocheCat# 05056489001
Chemical compound, drugRetinoic acidSigma-AldrichCat# R2625
Chemical compound, drugRNA ScreenTape Sample BufferAgilent TechnologiesCat# 5067–5577
Chemical compound, drugROCK Inhibitor Y-27632STEMCELL
Technologies
Cat# 72305
Chemical compound, drugRPMI 1640 [-] L-glutamineVWRCat# 15–040-CV
Chemical compound, drugSANT-1Sigma-AldrichCat# S4572
Chemical compound, drugSodium BicarbonateSigma-AldrichCat# NC0564699
Chemical compound, drugTamoxifenSigmaCat# T5648
Chemical compound, drugTGF-β RI Kinase Inhibitor IVCalbiochemCat# 616454
Chemical compound, drugTPBCalbiochemCat# 565740
Chemical compound, drugTranylcypromineCayman ChemicalCat# 10010494
Chemical compound, drugTris-acetateThermo Fisher ScientificCat# BP-152
Chemical compound, drugTTNPBEnzo Life SciencesCat# BML-GR105
Chemical compound, drugVectashield Antifade Mounting MediumVector LaboratoriesCat# H-1000
Chemical compound, drugXtremeGene 9RocheCat# 6365787001
Commercial assayHigh Sensitivity D1000 ScreenTapeAgilent TechnologiesCat# 5067–5584
Commercial assay, kitRNA ScreenTapeAgilent TechnologiesCat# 5067–5576
Commercial assay, kitRNA ScreenTape LadderAgilent TechnologiesCat# 5067–5578
Commercial assay, kitBD Cytofix/Cytoperm Plus Fixation/Permeabilization Solution KitBD BiosciencesCat# 554715
Commercial assay, kitChIP-IT High Sensitivity KitActive MotifCat# 53040
Commercial assay, kitiQ SYBR Green SupermixBio-RadCat# 1708884
Commercial assay, kitiScript cDNA Synthesis KitBio-RadCat# 1708891
Commercial assay, kitKAPA Library Preparation Kit (Illumina)Kapa BiosystemsCat# KK8234
Commercial assay, kitKAPA Stranded mRNA-Seq KitsKapa BiosystemsCat# KK8401
Commercial assay, kitMinElute PCR purification kitQIAGENCat# 28004
Commercial assay, kitQubit ssDNA assay kitThermo Fisher ScientificCat# Q10212
Commercial assay, kitRNeasy Micro KitQIAGENCat# 74004
Genetic reagent (D. rerio)Tg(ptf1a:eGFP)jh1PMID:16258076N/A
OtherSPRIselect beadBeckman CoulterCat# B23317
Recombinant proteinActivin AR and D SystemsCat# 338-AC/CF
Recombinant proteinHuman AB SerumValley BiomedicalCat# HP1022
Recombinant proteinRecombinant EGFR and D SystemsCat# 236-EG
Recombinant proteinRecombinant Heregulinβ−1PeprotechCat# 100–03
Recombinant proteinRecombinant KGF/FGF7R and D SystemsCat# 251 KG
Recombinant proteinRecombinant Mouse Wnt3AR and D SystemsCat# 1324-WN/CF
Recombinant proteinRecombinant NogginR and D SystemsCat# 3344 NG
Sequence-based reagentPx333 Plasmidhttp://www.addgene.org/64073/RRID:Addgene_64073
Sequence-based reagentLAMA1 ForwardThis paperqPCR primersGTG ATG GCA ACA GCG CAA A
Sequence-based reagentLAMA1 ReverseThis paperqPCR primersGAC CCA GTG ATA TTC TCT CCC A
Sequence-based reagentCRB2 ForwardThis paperqPCR primersACC ACT GTG CTT GTC CTG AG
Sequence-based reagentCRB2 ReverseThis paperqPCR primersTCC AGG GTC GCT AGA TGG AG
Sequence-based reagentTBP ForwardThis paperqPCR primersTGT GCA CAG GAG CCA AGA GT
Sequence-based reagentTBP ReverseThis paperqPCR primersATT TTC TTG CTG CCA GTC TGG
Sequence-based reagentLAMA1Enh Upstream GuideThis paperCRISPR sgRNAGTC AAA TTG CTA TAA CAC GG
Sequence-based reagentLAMA1Enh Downstream GuideThis paperCRISPR sgRNACCA CTT TAA GTA TCT CAG CA
Sequence-based reagentCRB2Enh Upstream GuideThis paperCRISPR sgRNAATA CAA AGC ACG TGA GA
Sequence-based reagentCRB2Enh Downstream GuideThis paperCRISPR sgRNAGAA TGC GGA TGA CGC CTG AG
Sequence-based reagentlama1-ATGPMID:16321372MorpholinoTCA TCC TCA TCT CCA TCA TCG CTC A
Obtained from GeneTools, LLC
Sequence-based reagentcrb2a-SPPMID:16713951MorpholinoACG TTG CCA GTA CCT GTG TAT CCT G
Obtained from GeneTools, LLC
Sequence-based reagentcrb2b-SPPMID:16713951MorpholinoTAA AGA TGT CCT ACC CAG CTT GAA C
Obtained from GeneTools, LLC
Sequence-based reagentstandard control MON/AMorpholinoCCT CTT ACC TCA GTT ACA ATT TAT A
Obtained from GeneTools, LLC
Software, algorithmAdobe Illustrator v 5.1http://www.adobe.com/products/illustrator.htmlRRID:SCR_014198
Software, algorithmAdobe Photoshop v 5.1http://www.adobe.com/products/photoshop.htmlRRID:SCR_014199
Software, algorithmBEDtools v 2.26.0https://github.com/arq5x/bedtools2RRID:SCR_006646
Software, algorithmBioconductorhttps://www.bioconductor.org/RRID: SCR_006442
Software, algorithmBurrows-Wheeler Aligner v 0.7.13http://bio-bwa.sourceforge.net/RRID:SCR_010910
Software, algorithmCENTIPEDE v 1.2http://centipede.uchicago.edu/N/A
Software, algorithmCufflinks v 2.2.1http://cole-trapnell-lab.github.io/cufflinks/RRID:SCR_014597
Software, algorithmdeepTools2 v 3.1.3https://deeptools.readthedocs.io/en/develop/content/installation.htmlN/A
Software, algorithmDESeq2 v 3.10https://bioconductor.org/packages/release/bioc/html/DESeq2.htmlRRID:SCR_015687
Software, algorithmFlowJo v10 softwarehttps://www.flowjo.com/solutions/flowjoRRID: SCR_008520
Software, algorithmGraphPad Prism v 8.1.2https://www.graphpad.com/scientific-software/prism/RRID: SCR_002798
Software, algorithmHOMER v 4.10.4http://homer.ucsd.edu/homer/RRID: SCR_010881
Software, algorithmJuicebox Tools v 1.4https://github.com/aidenlab/Juicebox/wiki/Juicebox-Assembly-ToolsN/A
Software, algorithmMACS2 v 2.1.4http://liulab.dfci.harvard.edu/MACS/RRID:SCR_013291
Software, algorithmMEME suite v 5.1.1http://meme-suite.org/RRID:SCR_001783
Software, algorithmMetascapehttp://metscape.ncibi.orgRRID:SCR_014687
Software, algorithmPicard Tools v 1.131http://broadinstitute.github.io/picard/RRID:SCR_006525
Software, algorithmR Project for Statistical Computing v 3.6.1http://www.r-project.org/RRID:SCR_001905
Software, algorithmSAMtools v 1.5http://samtools.sourceforge.netRRID:SCR_002105
Software, algorithmSTAR v 2.4https://github.com/alexdobin/STARN/A
Software, algorithmUCSC Genome Browserhttp://genome.ucsc.edu/RRID:SCR_005780
Software, algorithmvcf2diploid v 0.2.6ahttps://github.com/abyzovlab/vcf2diploidN/A
Software, algorithmZEISS ZEN Digital Imaging for Light Microscopyhttp://www.zeiss.com/microscopy/en_us/products/microscope-software/zen.html#introductionRRID:SCR_013672

Additional files

Download links