ABSTRACT
AIMS/HYPOTHESIS: Disruption of pancreatic islet function and glucose homeostasis can lead to the development of sustained hyperglycaemia, beta cell glucotoxicity and subsequently type 2 diabetes. In this study, we explored the effects of in vitro hyperglycaemic conditions on human pancreatic islet gene expression across 24 h in six pancreatic cell types: alpha; beta; gamma; delta; ductal; and acinar. We hypothesised that genes associated with hyperglycaemic conditions may be relevant to the onset and progression of diabetes. METHODS: We exposed human pancreatic islets from two donors to low (2.8 mmol/l) and high (15.0 mmol/l) glucose concentrations over 24 h in vitro. To assess the transcriptome, we performed single-cell RNA-seq (scRNA-seq) at seven time points. We modelled time as both a discrete and continuous variable to determine momentary and longitudinal changes in transcription associated with islet time in culture or glucose exposure. Additionally, we integrated genomic features and genetic summary statistics to nominate candidate effector genes. For three of these genes, we functionally characterised the effect on insulin production and secretion using CRISPR interference to knock down gene expression in EndoC-ßH1 cells, followed by a glucose-stimulated insulin secretion assay. RESULTS: In the discrete time models, we identified 1344 genes associated with time and 668 genes associated with glucose exposure across all cell types and time points. In the continuous time models, we identified 1311 genes associated with time, 345 genes associated with glucose exposure and 418 genes associated with interaction effects between time and glucose across all cell types. By integrating these expression profiles with summary statistics from genetic association studies, we identified 2449 candidate effector genes for type 2 diabetes, HbA1c, random blood glucose and fasting blood glucose. Of these candidate effector genes, we showed that three (ERO1B, HNRNPA2B1 and RHOBTB3) exhibited an effect on glucose-stimulated insulin production and secretion in EndoC-ßH1 cells. CONCLUSIONS/INTERPRETATION: The findings of our study provide an in-depth characterisation of the 24 h transcriptomic response of human pancreatic islets to glucose exposure at a single-cell resolution. By integrating differentially expressed genes with genetic signals for type 2 diabetes and glucose-related traits, we provide insights into the molecular mechanisms underlying glucose homeostasis. Finally, we provide functional evidence to support the role of three candidate effector genes in insulin secretion and production. DATA AVAILABILITY: The scRNA-seq data from the 24 h glucose exposure experiment performed in this study are available in the database of Genotypes and Phenotypes (dbGap; https://www.ncbi.nlm.nih.gov/gap/ ) with accession no. phs001188.v3.p1. Study metadata and summary statistics for the differential expression, gene set enrichment and candidate effector gene prediction analyses are available in the Zenodo data repository ( https://zenodo.org/ ) under accession number 11123248. The code used in this study is publicly available at https://github.com/CollinsLabBioComp/publication-islet_glucose_timecourse .
ABSTRACT
The hypothalamus, composed of several nuclei, is essential for maintaining our body's homeostasis. The arcuate nucleus (ARC), located in the mediobasal hypothalamus, contains neuronal populations with eminent roles in energy and glucose homeostasis as well as reproduction. These neuronal populations are of great interest for translational research. To fulfill this promise, we used a robotic cell culture platform to provide a scalable and chemically defined approach for differentiating human pluripotent stem cells (hPSCs) into pro-opiomelanocortin (POMC), somatostatin (SST), tyrosine hydroxylase (TH) and gonadotropin-releasing hormone (GnRH) neuronal subpopulations with an ARC-like signature. This robust approach is reproducible across several distinct hPSC lines and exhibits a stepwise induction of key ventral diencephalon and ARC markers in transcriptomic profiling experiments. This is further corroborated by direct comparison to human fetal hypothalamus, and the enriched expression of genes implicated in obesity and type 2 diabetes (T2D). Genome-wide chromatin accessibility profiling by ATAC-seq identified accessible regulatory regions that can be utilized to predict candidate enhancers related to metabolic disorders and hypothalamic development. In depth molecular, cellular, and functional experiments unveiled the responsiveness of the hPSC-derived hypothalamic neurons to hormonal stimuli, such as insulin, neuropeptides including kisspeptin, and incretin mimetic drugs such as Exendin-4, highlighting their potential utility as physiologically relevant cellular models for disease studies. In addition, differential glucose and insulin treatments uncovered adaptability within the generated ARC neurons in the dynamic regulation of POMC and insulin receptors. In summary, the establishment of this model represents a novel, chemically defined, and scalable platform for manufacturing large numbers of hypothalamic arcuate neurons and serves as a valuable resource for modeling metabolic and reproductive disorders.
ABSTRACT
We developed an efficient CRISPR prime editing protocol and generated isogenic-induced pluripotent stem cell (iPSC) lines carrying heterozygous or homozygous alleles for putatively causal single nucleotide variants at six type 2 diabetes loci (ABCC8, MTNR1B, TCF7L2, HNF4A, CAMK1D, and GCK). Our two-step sequence-based approach to first identify transfected cell pools with the highest fraction of edited cells significantly reduced the downstream efforts to isolate single clones of edited cells. We found that prime editing can make targeted genetic changes in iPSC and optimization of system components and guide RNA designs that were critical to achieve acceptable efficiency. Systems utilizing PEmax, epegRNA modifications, and MLH1dn provided significant benefit, producing editing efficiencies of 36-73%. Editing success and pegRNA design optimization required for each variant differed depending on the sequence at the target site. With attention to design, prime editing is a promising approach to generate isogenic iPSC lines, enabling the study of specific genetic changes in a common genetic background.
Subject(s)
Diabetes Mellitus, Type 2 , Induced Pluripotent Stem Cells , Humans , Clustered Regularly Interspaced Short Palindromic Repeats/genetics , CRISPR-Cas Systems/genetics , Gene Editing , RNA, Guide, CRISPR-Cas SystemsABSTRACT
Complete characterization of the genetic effects on gene expression is needed to elucidate tissue biology and the etiology of complex traits. Here, we analyzed 2,344 subcutaneous adipose tissue samples and identified 34K conditionally distinct expression quantitative trait locus (eQTL) signals in 18K genes. Over half of eQTL genes exhibited at least two eQTL signals. Compared to primary signals, non-primary signals had lower effect sizes, lower minor allele frequencies, and less promoter enrichment; they corresponded to genes with higher heritability and higher tolerance for loss of function. Colocalization of eQTL with conditionally distinct genome-wide association study signals for 28 cardiometabolic traits identified 3,605 eQTL signals for 1,861 genes. Inclusion of non-primary eQTL signals increased colocalized signals by 46%. Among 30 genes with ≥2 pairs of colocalized signals, 21 showed a mediating gene dosage effect on the trait. Thus, expanded eQTL identification reveals more mechanisms underlying complex traits and improves understanding of the complexity of gene expression regulation.
ABSTRACT
Genome wide association studies (GWAS) have identified over 100 signals associated with type 1 diabetes (T1D). However, translating any given T1D GWAS signal into mechanistic insights, including putative causal variants and the context (cell type and cell state) in which they function, has been limited. Here, we present a comprehensive multi-omic integrative analysis of single-cell/nucleus resolution profiles of gene expression and chromatin accessibility in healthy and autoantibody+ (AAB+) human islets, as well as islets under multiple T1D stimulatory conditions. We broadly nominate effector cell types for all T1D GWAS signals. We further nominated higher-resolution contexts, including effector cell types, regulatory elements, and genes for three independent T1D risk variants acting through islet cells within the pancreas at the DLK1/MEG3, RASGRP1, and TOX loci. Subsequently, we created isogenic gene knockouts DLK1-/-, RASGRP1-/-, and TOX-/-, and the corresponding regulatory region knockout, RASGRP1Δ, and DLK1Δ hESCs. Loss of RASGRP1 or DLK1, as well as knockout of the regulatory region of RASGRP1 or DLK1, increased ß cell apoptosis. Additionally, pancreatic ß cells derived from isogenic hESCs carrying the risk allele of rs3783355A/A exhibited increased ß cell death. Finally, RNA-seq and ATAC-seq identified five genes upregulated in both RASGRP1-/- and DLK1-/- ß-like cells, four of which are associated with T1D. Together, this work reports an integrative approach for combining single cell multi-omics, GWAS, and isogenic hESC-derived ß-like cells to prioritize the T1D associated signals and their underlying context-specific cell types, genes, SNPs, and regulatory elements, to illuminate biological functions and molecular mechanisms.
ABSTRACT
Genetic studies have identified numerous loci associated with type 2 diabetes (T2D), but the functional roles of many loci remain unexplored. Here, we engineered isogenic knockout human embryonic stem cell lines for 20 genes associated with T2D risk. We examined the impacts of each knockout on ß cell differentiation, functions, and survival. We generated gene expression and chromatin accessibility profiles on ß cells derived from each knockout line. Analyses of T2D-association signals overlapping HNF4A-dependent ATAC peaks identified a likely causal variant at the FAIM2 T2D-association signal. Additionally, the integrative association analyses identified four genes (CP, RNASE1, PCSK1N, and GSTA2) associated with insulin production, and two genes (TAGLN3 and DHRS2) associated with ß cell sensitivity to lipotoxicity. Finally, we leveraged deep ATAC-seq read coverage to assess allele-specific imbalance at variants heterozygous in the parental line and identified a single likely functional variant at each of 23 T2D-association signals.
Subject(s)
Diabetes Mellitus, Type 2 , Human Embryonic Stem Cells , Insulin-Secreting Cells , Humans , Diabetes Mellitus, Type 2/genetics , Diabetes Mellitus, Type 2/metabolism , Human Embryonic Stem Cells/metabolism , Genetic Predisposition to Disease , Genome-Wide Association Study , Insulin-Secreting Cells/metabolism , Polymorphism, Single Nucleotide , Carbonyl Reductase (NADPH)/genetics , Carbonyl Reductase (NADPH)/metabolismABSTRACT
Genetic association studies have identified hundreds of independent signals associated with type 2 diabetes (T2D) and related traits. Despite these successes, the identification of specific causal variants underlying a genetic association signal remains challenging. In this study, we describe a deep learning (DL) method to analyze the impact of sequence variants on enhancers. Focusing on pancreatic islets, a T2D relevant tissue, we show that our model learns islet-specific transcription factor (TF) regulatory patterns and can be used to prioritize candidate causal variants. At 101 genetic signals associated with T2D and related glycemic traits where multiple variants occur in linkage disequilibrium, our method nominates a single causal variant for each association signal, including three variants previously shown to alter reporter activity in islet-relevant cell types. For another signal associated with blood glucose levels, we biochemically test all candidate causal variants from statistical fine-mapping using a pancreatic islet beta cell line and show biochemical evidence of allelic effects on TF binding for the model-prioritized variant. To aid in future research, we publicly distribute our model and islet enhancer perturbation scores across ~67 million genetic variants. We anticipate that DL methods like the one presented in this study will enhance the prioritization of candidate causal variants for functional studies.
Subject(s)
Deep Learning , Diabetes Mellitus, Type 2 , Enhancer Elements, Genetic , Islets of Langerhans , Diabetes Mellitus, Type 2/genetics , Diabetes Mellitus, Type 2/metabolism , Diabetes Mellitus, Type 2/pathology , Islets of Langerhans/metabolism , Islets of Langerhans/pathology , Genetic Variation , Humans , Computer SimulationABSTRACT
Disruption of pancreatic islet function and glucose homeostasis can lead to the development of sustained hyperglycemia, beta cell glucotoxicity, and ultimately type 2 diabetes (T2D). In this study, we sought to explore the effects of hyperglycemia on human pancreatic islet (HPI) gene expression by exposing HPIs from two donors to low (2.8mM) and high (15.0mM) glucose concentrations over 24 hours, assaying the transcriptome at seven time points using single-cell RNA sequencing (scRNA-seq). We modeled time as both a discrete and continuous variable to determine momentary and longitudinal changes in transcription associated with islet time in culture or glucose exposure. Across all cell types, we identified 1,528 genes associated with time, 1,185 genes associated with glucose exposure, and 845 genes associated with interaction effects between time and glucose. We clustered differentially expressed genes across cell types and found 347 modules of genes with similar expression patterns across time and glucose conditions, including two beta cell modules enriched in genes associated with T2D. Finally, by integrating genomic features from this study and genetic summary statistics for T2D and related traits, we nominate 363 candidate effector genes that may underlie genetic associations for T2D and related traits.
ABSTRACT
Hutchinson-Gilford progeria syndrome (HGPS) is a premature aging disorder affecting tissues of mesenchymal origin. Most individuals with HGPS harbor a de novo c.1824C > T (p.G608G) mutation in the gene encoding lamin A (LMNA), which activates a cryptic splice donor site resulting in production of the toxic "progerin" protein. Clinical manifestations include growth deficiency, lipodystrophy, sclerotic dermis, cardiovascular defects, and bone dysplasia. Here we utilized the LmnaG609G knock-in (KI) mouse model of HGPS to further define mechanisms of bone loss associated with normal and premature aging disorders. Newborn skeletal staining of KI mice revealed altered rib cage shape and spinal curvature, and delayed calvarial mineralization with increased craniofacial and mandibular cartilage content. MicroCT analysis and mechanical testing of adult femurs indicated increased fragility associated with reduced bone mass, recapitulating the progressive bone deterioration that occurs in HGPS patients. We investigated mechanisms of bone loss in KI mice at the cellular level in bone cell populations. Formation of wild-type and KI osteoclasts from marrow-derived precursors was inhibited by KI osteoblast-conditioned media in vitro, suggesting a secreted factor(s) responsible for decreased osteoclasts on KI trabecular surfaces in vivo. Cultured KI osteoblasts exhibited abnormal differentiation characterized by reduced deposition and mineralization of extracellular matrix with increased lipid accumulation compared to wild-type, providing a mechanism for altered bone formation. Furthermore, quantitative analyses of KI transcripts confirmed upregulation of adipogenic genes both in vitro and in vivo. Thus, osteoblast phenotypic plasticity, inflammation and altered cellular cross-talk contribute to abnormal bone formation in HGPS mice.
Subject(s)
Aging, Premature , Bone Diseases, Developmental , Progeria , Mice , Animals , Progeria/genetics , Progeria/metabolism , Mutation , Lamin Type A/genetics , Lamin Type A/metabolism , Cell DifferentiationABSTRACT
Hereditary congenital facial paresis type 1 (HCFP1) is an autosomal dominant disorder of absent or limited facial movement that maps to chromosome 3q21-q22 and is hypothesized to result from facial branchial motor neuron (FBMN) maldevelopment. In the present study, we report that HCFP1 results from heterozygous duplications within a neuron-specific GATA2 regulatory region that includes two enhancers and one silencer, and from noncoding single-nucleotide variants (SNVs) within the silencer. Some SNVs impair binding of NR2F1 to the silencer in vitro and in vivo and attenuate in vivo enhancer reporter expression in FBMNs. Gata2 and its effector Gata3 are essential for inner-ear efferent neuron (IEE) but not FBMN development. A humanized HCFP1 mouse model extends Gata2 expression, favors the formation of IEEs over FBMNs and is rescued by conditional loss of Gata3. These findings highlight the importance of temporal gene regulation in development and of noncoding variation in rare mendelian disease.
Subject(s)
Facial Paralysis , Animals , Mice , Facial Paralysis/genetics , Facial Paralysis/congenital , Facial Paralysis/metabolism , GATA2 Transcription Factor/genetics , GATA2 Transcription Factor/metabolism , Motor Neurons/metabolism , Neurogenesis , Neurons, EfferentABSTRACT
Genetic studies have identified numerous loci associated with type 2 diabetes (T2D), but the functional role of many loci has remained unexplored. In this study, we engineered isogenic knockout human embryonic stem cell (hESC) lines for 20 genes associated with T2D risk. We systematically examined ß-cell differentiation, insulin production and secretion, and survival. We performed RNA-seq and ATAC-seq on hESC-ß cells from each knockout line. Analyses of T2D GWAS signals overlapping with HNF4A-dependent ATAC peaks identified a specific SNP as a likely causal variant. In addition, we performed integrative association analyses and identified four genes ( CP, RNASE1, PCSK1N and GSTA2 ) associated with insulin production, and two genes ( TAGLN3 and DHRS2 ) associated with sensitivity to lipotoxicity. Finally, we leveraged deep ATAC-seq read coverage to assess allele-specific imbalance at variants heterozygous in the parental hESC line, to identify a single likely functional variant at each of 23 T2D GWAS signals.
ABSTRACT
Genetic studies have identified ≥240 loci associated with the risk of type 2 diabetes (T2D), yet most of these loci lie in non-coding regions, masking the underlying molecular mechanisms. Recent studies investigating mRNA expression in human pancreatic islets have yielded important insights into the molecular drivers of normal islet function and T2D pathophysiology. However, similar studies investigating microRNA (miRNA) expression remain limited. Here, we present data from 63 individuals, the largest sequencing-based analysis of miRNA expression in human islets to date. We characterized the genetic regulation of miRNA expression by decomposing the expression of highly heritable miRNAs into cis- and trans-acting genetic components and mapping cis-acting loci associated with miRNA expression [miRNA-expression quantitative trait loci (eQTLs)]. We found i) 84 heritable miRNAs, primarily regulated by trans-acting genetic effects, and ii) 5 miRNA-eQTLs. We also used several different strategies to identify T2D-associated miRNAs. First, we colocalized miRNA-eQTLs with genetic loci associated with T2D and multiple glycemic traits, identifying one miRNA, miR-1908, that shares genetic signals for blood glucose and glycated hemoglobin (HbA1c). Next, we intersected miRNA seed regions and predicted target sites with credible set SNPs associated with T2D and glycemic traits and found 32 miRNAs that may have altered binding and function due to disrupted seed regions. Finally, we performed differential expression analysis and identified 14 miRNAs associated with T2D status-including miR-187-3p, miR-21-5p, miR-668, and miR-199b-5p-and 4 miRNAs associated with a polygenic score for HbA1c levels-miR-216a, miR-25, miR-30a-3p, and miR-30a-5p.
Subject(s)
Diabetes Mellitus, Type 2 , Islets of Langerhans , MicroRNAs , Humans , MicroRNAs/metabolism , Diabetes Mellitus, Type 2/genetics , Diabetes Mellitus, Type 2/metabolism , Glycated Hemoglobin , Islets of Langerhans/metabolism , Quantitative Trait Loci/geneticsABSTRACT
Skeletal muscle, the largest human organ by weight, is relevant to several polygenic metabolic traits and diseases including type 2 diabetes (T2D). Identifying genetic mechanisms underlying these traits requires pinpointing the relevant cell types, regulatory elements, target genes, and causal variants. Here, we used genetic multiplexing to generate population-scale single nucleus (sn) chromatin accessibility (snATAC-seq) and transcriptome (snRNA-seq) maps across 287 frozen human skeletal muscle biopsies representing 456,880 nuclei. We identified 13 cell types that collectively represented 983,155 ATAC summits. We integrated genetic variation to discover 6,866 expression quantitative trait loci (eQTL) and 100,928 chromatin accessibility QTL (caQTL) (5% FDR) across the five most abundant cell types, cataloging caQTL peaks that atlas-level snATAC maps often miss. We identified 1,973 eGenes colocalized with caQTL and used mediation analyses to construct causal directional maps for chromatin accessibility and gene expression. 3,378 genome-wide association study (GWAS) signals across 43 relevant traits colocalized with sn-e/caQTL, 52% in a cell-specific manner. 77% of GWAS signals colocalized with caQTL and not eQTL, highlighting the critical importance of population-scale chromatin profiling for GWAS functional studies. GWAS-caQTL colocalization showed distinct cell-specific regulatory paradigms. For example, a C2CD4A/B T2D GWAS signal colocalized with caQTL in muscle fibers and multiple chromatin loop models nominated VPS13C, a glucose uptake gene. Sequence of the caQTL peak overlapping caSNP rs7163757 showed allelic regulatory activity differences in a human myocyte cell line massively parallel reporter assay. These results illuminate the genetic regulatory architecture of human skeletal muscle at high-resolution epigenomic, transcriptomic, and cell state scales and serve as a template for population-scale multi-omic mapping in complex tissues and traits.
ABSTRACT
BACKGROUND: COVID-19 severity varies widely. Although some demographic and cardio-metabolic factors, including age and obesity, are associated with increasing risk of severe illness, the underlying mechanism(s) are uncertain. SUBJECTS/METHODS: In a meta-analysis of three independent studies of 1471 participants in total, we investigated phenotypic and genetic factors associated with subcutaneous adipose tissue expression of Angiotensin I Converting Enzyme 2 (ACE2), measured by RNA-Seq, which acts as a receptor for SARS-CoV-2 cellular entry. RESULTS: Lower adipose tissue ACE2 expression was associated with multiple adverse cardio-metabolic health indices, including type 2 diabetes (T2D) (P = 9.14 × 10-6), obesity status (P = 4.81 × 10-5), higher serum fasting insulin (P = 5.32 × 10-4), BMI (P = 3.94 × 10-4), and lower serum HDL levels (P = 1.92 × 10-7). ACE2 expression was also associated with estimated proportions of cell types in adipose tissue: lower expression was associated with a lower proportion of microvascular endothelial cells (P = 4.25 × 10-4) and higher proportion of macrophages (P = 2.74 × 10-5). Despite an estimated heritability of 32%, we did not identify any proximal or distal expression quantitative trait loci (eQTLs) associated with adipose tissue ACE2 expression. CONCLUSIONS: Our results demonstrate that individuals with cardio-metabolic features known to increase risk of severe COVID-19 have lower background ACE2 levels in this highly relevant tissue. Reduced adipose tissue ACE2 expression may contribute to the pathophysiology of cardio-metabolic diseases, as well as the associated increased risk of severe COVID-19.
Subject(s)
Adipose Tissue , Angiotensin-Converting Enzyme 2 , COVID-19 , Adipose Tissue/metabolism , Angiotensin-Converting Enzyme 2/genetics , Angiotensin-Converting Enzyme 2/metabolism , COVID-19/complications , COVID-19/genetics , Cardiometabolic Risk Factors , Diabetes Mellitus, Type 2/genetics , Endothelial Cells/metabolism , Humans , Obesity , SARS-CoV-2ABSTRACT
Recent clinical data have suggested a correlation between coronavirus disease 2019 (COVID-19) and diabetes. Here, we describe the detection of SARS-CoV-2 viral antigen in pancreatic beta cells in autopsy samples from individuals with COVID-19. Single-cell RNA sequencing and immunostaining from ex vivo infections confirmed that multiple types of pancreatic islet cells were susceptible to SARS-CoV-2, eliciting a cellular stress response and the induction of chemokines. Upon SARS-CoV-2 infection, beta cells showed a lower expression of insulin and a higher expression of alpha and acinar cell markers, including glucagon and trypsin1, respectively, suggesting cellular transdifferentiation. Trajectory analysis indicated that SARS-CoV-2 induced eIF2-pathway-mediated beta cell transdifferentiation, a phenotype that could be reversed with trans-integrated stress response inhibitor (trans-ISRIB). Altogether, this study demonstrates an example of SARS-CoV-2 infection causing cell fate change, which provides further insight into the pathomechanisms of COVID-19.
Subject(s)
COVID-19/virology , Cell Transdifferentiation , Insulin-Secreting Cells/virology , SARS-CoV-2/pathogenicity , Acetamides/pharmacology , Adolescent , Adult , Aged , Aged, 80 and over , Animals , COVID-19/mortality , Cell Transdifferentiation/drug effects , Chlorocebus aethiops , Cyclohexylamines/pharmacology , Cytokines/metabolism , Eukaryotic Initiation Factor-2/metabolism , Female , Glucagon , Host-Pathogen Interactions , Humans , Insulin/metabolism , Insulin-Secreting Cells/drug effects , Insulin-Secreting Cells/metabolism , Insulin-Secreting Cells/pathology , Male , Middle Aged , Phenotype , Signal Transduction , Tissue Culture Techniques , Trypsin/metabolism , Vero Cells , Young AdultABSTRACT
Identifying the molecular mechanisms by which genome-wide association study (GWAS) loci influence traits remains challenging. Chromatin accessibility quantitative trait loci (caQTLs) help identify GWAS loci that may alter GWAS traits by modulating chromatin structure, but caQTLs have been identified in a limited set of human tissues. Here we mapped caQTLs in human liver tissue in 20 liver samples and identified 3,123 caQTLs. The caQTL variants are enriched in liver tissue promoter and enhancer states and frequently disrupt binding motifs of transcription factors expressed in liver. We predicted target genes for 861 caQTL peaks using proximity, chromatin interactions, correlation with promoter accessibility or gene expression, and colocalization with expression QTLs. Using GWAS signals for 19 liver function and/or cardiometabolic traits, we identified 110 colocalized caQTLs and GWAS signals, 56 of which contained a predicted caPeak target gene. At the LITAF LDL-cholesterol GWAS locus, we validated that a caQTL variant showed allelic differences in protein binding and transcriptional activity. These caQTLs contribute to the epigenomic characterization of human liver and help identify molecular mechanisms and genes at GWAS loci.
Subject(s)
Chromatin/metabolism , Liver/metabolism , Quantitative Trait Loci , Amino Acid Motifs , Binding Sites , Chromatin Assembly and Disassembly , Enhancer Elements, Genetic , Genetic Variation , Genome-Wide Association Study , Humans , Promoter Regions, Genetic , Protein Binding , Transcription Factors/chemistry , Transcription Factors/metabolism , TranscriptomeABSTRACT
Identifying the tissue-specific molecular signatures of active regulatory elements is critical to understand gene regulatory mechanisms. Here, we identify transcription start sites (TSS) using cap analysis of gene expression (CAGE) across 57 human pancreatic islet samples. We identify 9,954 reproducible CAGE tag clusters (TCs), â¼20% of which are islet specific and occur mostly distal to known gene TSS. We integrated islet CAGE data with histone modification and chromatin accessibility profiles to identify epigenomic signatures of transcription initiation. Using a massively parallel reporter assay, we validated the transcriptional enhancer activity for 2,279 of 3,378 (â¼68%) tested islet CAGE elements (5% false discovery rate). TCs within accessible enhancers show higher enrichment to overlap type 2 diabetes genome-wide association study (GWAS) signals than existing islet annotations, which emphasizes the utility of mapping CAGE profiles in disease-relevant tissue. This work provides a high-resolution map of transcriptional initiation in human pancreatic islets with utility for dissecting active enhancers at GWAS loci.
Subject(s)
Islets of Langerhans/physiology , Transcription Initiation Site , Enhancer Elements, Genetic , Genome-Wide Association Study , Humans , Polymorphism, Single Nucleotide , Quantitative Trait LociABSTRACT
Hutchinson-Gilford progeria syndrome (HGPS) is a rare accelerated aging disorder characterized by premature death from myocardial infarction or stroke. It is caused by de novo single-nucleotide mutations in the LMNA gene that activate a cryptic splice donor site, resulting in the production of a toxic form of lamin A, which is termed progerin. Here we present a potential genetic therapeutic strategy that utilizes antisense peptide-conjugated phosphorodiamidate morpholino oligomers (PPMOs) to block pathogenic splicing of mutant transcripts. Of several candidates, PPMO SRP-2001 provided the most significant decrease in progerin transcripts in patient fibroblasts. Intravenous delivery of SRP-2001 to a transgenic mouse model of HGPS produced significant reduction of progerin transcripts in the aorta, a particularly critical target tissue in HGPS. Long-term continuous treatment with SRP-2001 yielded a 61.6% increase in lifespan and rescue of vascular smooth muscle cell loss in large arteries. These results provide a rationale for proceeding to human trials.
Subject(s)
Oligonucleotides, Antisense/therapeutic use , Progeria/drug therapy , Animals , Disease Models, Animal , Humans , Mice , Mice, Transgenic , Morpholinos/chemistryABSTRACT
Hutchinson-Gilford progeria syndrome (HGPS or progeria) is typically caused by a dominant-negative Câ¢G-to-Tâ¢A mutation (c.1824 C>T; p.G608G) in LMNA, the gene that encodes nuclear lamin A. This mutation causes RNA mis-splicing that produces progerin, a toxic protein that induces rapid ageing and shortens the lifespan of children with progeria to approximately 14 years1-4. Adenine base editors (ABEs) convert targeted Aâ¢T base pairs to Gâ¢C base pairs with minimal by-products and without requiring double-strand DNA breaks or donor DNA templates5,6. Here we describe the use of an ABE to directly correct the pathogenic HGPS mutation in cultured fibroblasts derived from children with progeria and in a mouse model of HGPS. Lentiviral delivery of the ABE to fibroblasts from children with HGPS resulted in 87-91% correction of the pathogenic allele, mitigation of RNA mis-splicing, reduced levels of progerin and correction of nuclear abnormalities. Unbiased off-target DNA and RNA editing analysis did not detect off-target editing in treated patient-derived fibroblasts. In transgenic mice that are homozygous for the human LMNA c.1824 C>T allele, a single retro-orbital injection of adeno-associated virus 9 (AAV9) encoding the ABE resulted in substantial, durable correction of the pathogenic mutation (around 20-60% across various organs six months after injection), restoration of normal RNA splicing and reduction of progerin protein levels. In vivo base editing rescued the vascular pathology of the mice, preserving vascular smooth muscle cell counts and preventing adventitial fibrosis. A single injection of ABE-expressing AAV9 at postnatal day 14 improved vitality and greatly extended the median lifespan of the mice from 215 to 510 days. These findings demonstrate the potential of in vivo base editing as a possible treatment for HGPS and other genetic diseases by directly correcting their root cause.
Subject(s)
Adenine/metabolism , Gene Editing/methods , Mutation , Progeria/genetics , Progeria/therapy , Alleles , Alternative Splicing , Animals , Aorta/pathology , Base Pairing , Child , DNA/genetics , Disease Models, Animal , Female , Fibroblasts/metabolism , Humans , Lamin Type A/chemistry , Lamin Type A/genetics , Lamin Type A/metabolism , Longevity , Male , Mice , Mice, Transgenic , Mutant Proteins/chemistry , Mutant Proteins/genetics , Mutant Proteins/metabolism , Progeria/pathology , RNA/geneticsABSTRACT
Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, key tissues and cell-types required for functional inference are absent from large-scale resources. Here we explore the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using data from 420 donors. We find: (a) 7741 cis-eQTLs in islets with a replication rate across 44 GTEx tissues between 40% and 73%; (b) marked overlap between islet cis-eQTL signals and active regulatory sequences in islets, with reduced eQTL effect size observed in the stretch enhancers most strongly implicated in GWAS signal location; (c) enrichment of islet cis-eQTL signals with T2D risk variants identified in genome-wide association studies; and (d) colocalization between 47 islet cis-eQTLs and variants influencing T2D or glycemic traits, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in disease relevant tissues.