ABSTRACT
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
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
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
Transcriptomics data have been integrated with genome-wide association studies (GWASs) to help understand disease/trait molecular mechanisms. The utility of metabolomics, integrated with transcriptomics and disease GWASs, to understand molecular mechanisms for metabolite levels or diseases has not been thoroughly evaluated. We performed probabilistic transcriptome-wide association and locus-level colocalization analyses to integrate transcriptomics results for 49 tissues in 706 individuals from the GTEx project, metabolomics results for 1,391 plasma metabolites in 6,136 Finnish men from the METSIM study, and GWAS results for 2,861 disease traits in 260,405 Finnish individuals from the FinnGen study. We found that genetic variants that regulate metabolite levels were more likely to influence gene expression and disease risk compared to the ones that do not. Integrating transcriptomics with metabolomics results prioritized 397 genes for 521 metabolites, including 496 previously identified gene-metabolite pairs with strong functional connections and suggested 33.3% of such gene-metabolite pairs shared the same causal variants with genetic associations of gene expression. Integrating transcriptomics and metabolomics individually with FinnGen GWAS results identified 1,597 genes for 790 disease traits. Integrating transcriptomics and metabolomics jointly with FinnGen GWAS results helped pinpoint metabolic pathways from genes to diseases. We identified putative causal effects of UGT1A1/UGT1A4 expression on gallbladder disorders through regulating plasma (E,E)-bilirubin levels, of SLC22A5 expression on nasal polyps and plasma carnitine levels through distinct pathways, and of LIPC expression on age-related macular degeneration through glycerophospholipid metabolic pathways. Our study highlights the power of integrating multiple sets of molecular traits and GWAS results to deepen understanding of disease pathophysiology.
Subject(s)
Genome-Wide Association Study , Transcriptome , Bilirubin , Carnitine , Glycerophospholipids , Humans , Male , Metabolomics , Quantitative Trait Loci/genetics , Solute Carrier Family 22 Member 5/genetics , Transcriptome/geneticsABSTRACT
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
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
Hutchinson-Gilford progeria syndrome (HGPS) is a uniformly fatal condition that is especially prevalent in skin, cardiovascular, and musculoskeletal systems. A wide gap exists between our knowledge of the disease and a promising treatment or cure. The aim of this study was to first characterize the musculoskeletal phenotype of the homozygous G608G BAC-transgenic progeria mouse model, and to determine the phenotype changes of HGPS mice after a five-arm preclinical trial of different treatment combinations with lonafarnib, pravastatin, and zoledronic acid. Microcomputed tomography and CT-based rigidity analyses were performed to assess cortical and trabecular bone structure, density, and rigidity. Bones were loaded to failure with three-point bending to assess strength. Contrast-enhanced µCT imaging of mouse femurs was performed to measure glycosaminoglycan content, thickness, and volume of the femoral head articular cartilage. Advanced glycation end products were assessed with a fluorometric assay. The changes demonstrated in the cortical bone structure, rigidity, stiffness, and modulus of the HGPS G608G mouse model may increase the risk for bending and deformation, which could result in the skeletal dysplasia characteristic of HGPS. Cartilage abnormalities seen in this HGPS model resemble changes observed in the age-matched WT controls, including early loss of glycosaminoglycans, and decreased cartilage thickness and volume. Such changes might mimic prevalent degenerative joint diseases in the elderly. Lonafarnib monotherapy did not improve bone or cartilage parameters, but treatment combinations with pravastatin and zoledronic acid significantly improved bone structure and mechanical properties and cartilage structural parameters, which ameliorate the musculoskeletal phenotype of the disease.
Subject(s)
Bone Density Conservation Agents/therapeutic use , Disease Models, Animal , Lamin Type A/genetics , Progeria , Aging/drug effects , Aging/pathology , Animals , Bone and Bones/drug effects , Bone and Bones/pathology , Cartilage/drug effects , Cartilage/pathology , Femur/drug effects , Femur/pathology , Glycosaminoglycans/analysis , Joints/drug effects , Joints/pathology , Lamin Type A/metabolism , Mice , Mice, Transgenic , Mutation , Osteoarthritis/drug therapy , Osteoarthritis/pathology , Phenotype , Piperidines/therapeutic use , Pravastatin/therapeutic use , Progeria/drug therapy , Progeria/genetics , Protein Processing, Post-Translational/drug effects , Pyridines/therapeutic use , X-Ray Microtomography , Zoledronic Acid/therapeutic useABSTRACT
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
We integrate comeasured gene expression and DNA methylation (DNAme) in 265 human skeletal muscle biopsies from the FUSION study with >7 million genetic variants and eight physiological traits: height, waist, weight, waist-hip ratio, body mass index, fasting serum insulin, fasting plasma glucose, and type 2 diabetes. We find hundreds of genes and DNAme sites associated with fasting insulin, waist, and body mass index, as well as thousands of DNAme sites associated with gene expression (eQTM). We find that controlling for heterogeneity in tissue/muscle fiber type reduces the number of physiological trait associations, and that long-range eQTMs (>1 Mb) are reduced when controlling for tissue/muscle fiber type or latent factors. We map genetic regulators (quantitative trait loci; QTLs) of expression (eQTLs) and DNAme (mQTLs). Using Mendelian randomization (MR) and mediation techniques, we leverage these genetic maps to predict 213 causal relationships between expression and DNAme, approximately two-thirds of which predict methylation to causally influence expression. We use MR to integrate FUSION mQTLs, FUSION eQTLs, and GTEx eQTLs for 48 tissues with genetic associations for 534 diseases and quantitative traits. We identify hundreds of genes and thousands of DNAme sites that may drive the reported disease/quantitative trait genetic associations. We identify 300 gene expression MR associations that are present in both FUSION and GTEx skeletal muscle and that show stronger evidence of MR association in skeletal muscle than other tissues, which may partially reflect differences in power across tissues. As one example, we find that increased RXRA muscle expression may decrease lean tissue mass.
Subject(s)
DNA Methylation/genetics , Gene Expression/genetics , Muscle, Skeletal , Blood Glucose/analysis , Body Weights and Measures , Diabetes Mellitus, Type 2 , Genome-Wide Association Study/methods , Genomics/methods , Humans , Insulin/analysis , Muscle, Skeletal/chemistry , Muscle, Skeletal/physiology , Quantitative Trait Loci/geneticsABSTRACT
The high-throughput detection of organelle composition and proteomic mapping of protein environment directly from primary tissue as well as the identification of interactors of insoluble proteins that form higher-order structures have remained challenges in biological research. We report a proximity-based labeling approach that uses an antibody to a target antigen to guide biotin deposition onto adjacent proteins in fixed cells and primary tissues, which allows proteins in close proximity to the target antigen to be captured and identified by mass spectrometry. We demonstrated the specificity and sensitivity of our method by examining the well-studied mitochondrial matrix. We then used the method to profile the dynamic interactome of lamin A/C in multiple cell and tissue types under various treatment conditions. The ability to detect proximal proteins and putative interactors in intact tissues, and to quantify changes caused by different conditions or in the presence of disease mutations, can provide a window into cell biology and disease pathogenesis.
Subject(s)
Antibodies/metabolism , Biotin/metabolism , Biotinylation/methods , Isotope Labeling/methods , Protein Interaction Mapping/methods , Proteins/metabolism , Proteomics/methods , Animals , HeLa Cells , Humans , Mice , Mice, Inbred C57BLABSTRACT
Enhancers regulate spatiotemporal gene expression and impart cell-specific transcriptional outputs that drive cell identity. Super-enhancers (SEs), also known as stretch-enhancers, are a subset of enhancers especially important for genes associated with cell identity and genetic risk of disease. CD4(+) T cells are critical for host defence and autoimmunity. Here we analysed maps of mouse T-cell SEs as a non-biased means of identifying key regulatory nodes involved in cell specification. We found that cytokines and cytokine receptors were the dominant class of genes exhibiting SE architecture in T cells. Nonetheless, the locus encoding Bach2, a key negative regulator of effector differentiation, emerged as the most prominent T-cell SE, revealing a network in which SE-associated genes critical for T-cell biology are repressed by BACH2. Disease-associated single-nucleotide polymorphisms for immune-mediated disorders, including rheumatoid arthritis, were highly enriched for T-cell SEs versus typical enhancers or SEs in other cell lineages. Intriguingly, treatment of T cells with the Janus kinase (JAK) inhibitor tofacitinib disproportionately altered the expression of rheumatoid arthritis risk genes with SE structures. Together, these results indicate that genes with SE architecture in T cells encompass a variety of cytokines and cytokine receptors but are controlled by a 'guardian' transcription factor, itself endowed with an SE. Thus, enumeration of SEs allows the unbiased determination of key regulatory nodes in T cells, which are preferentially modulated by pharmacological intervention.
Subject(s)
Arthritis, Rheumatoid/genetics , Enhancer Elements, Genetic/genetics , T-Lymphocytes, Helper-Inducer/metabolism , T-Lymphocytes, Helper-Inducer/pathology , Animals , Arthritis, Rheumatoid/immunology , Arthritis, Rheumatoid/pathology , Basic-Leucine Zipper Transcription Factors/metabolism , Cell Differentiation/genetics , Cell Lineage/genetics , Gene Expression Regulation/genetics , Genetic Predisposition to Disease/genetics , Janus Kinase 3/antagonists & inhibitors , Mice , Mice, Inbred C57BL , Piperidines/pharmacology , Pyrimidines/pharmacology , Pyrroles/pharmacology , RNA, Untranslated/genetics , T-Lymphocytes, Helper-Inducer/immunology , Transcription, Genetic/genetics , p300-CBP Transcription Factors/metabolismABSTRACT
LMNA encodes the A-type lamins that are part of the nuclear scaffold. Mutations in LMNA can cause a variety of disorders called laminopathies, including Hutchinson-Gilford progeria syndrome (HGPS), atypical Werner syndrome, and Emery-Dreifuss muscular dystrophy. Previous work has shown that treatment of HGPS cells with the mTOR inhibitor rapamycin or with the rapamycin analog everolimus corrects several of the phenotypes seen at the cellular level-at least in part by increasing autophagy and reducing the amount of progerin, the toxic form of lamin A that is overproduced in HGPS patients. Since other laminopathies also result in production of abnormal and potentially toxic lamin proteins, we hypothesized that everolimus would also be beneficial in those disorders. To test this, we applied everolimus to fibroblast cell lines from six laminopathy patients, each with a different mutation in LMNA Everolimus treatment increased proliferative ability and delayed senescence in all cell lines. In several cell lines, we observed that with treatment, there is a significant improvement in nuclear blebbing, which is a cellular hallmark of HGPS and other lamin disorders. These preclinical results suggest that everolimus might have clinical benefit for multiple laminopathy syndromes.
Subject(s)
Everolimus/pharmacology , Fibroblasts/drug effects , Lamin Type A/deficiency , Muscular Dystrophy, Emery-Dreifuss/genetics , Progeria/genetics , TOR Serine-Threonine Kinases/antagonists & inhibitors , Werner Syndrome/genetics , Biomarkers , Cell Division/drug effects , Cell Line , Cell Nucleus/drug effects , Cell Nucleus/ultrastructure , Cellular Senescence/drug effects , Humans , Lamin Type A/genetics , Muscular Dystrophy, Emery-Dreifuss/pathology , Mutation , Phosphorylation/drug effects , Progeria/pathology , Protein Processing, Post-Translational/drug effects , Ribosomal Protein S6/metabolism , Werner Syndrome/pathologyABSTRACT
Genome-wide association studies (GWAS) have identified >100 independent SNPs that modulate the risk of type 2 diabetes (T2D) and related traits. However, the pathogenic mechanisms of most of these SNPs remain elusive. Here, we examined genomic, epigenomic, and transcriptomic profiles in human pancreatic islets to understand the links between genetic variation, chromatin landscape, and gene expression in the context of T2D. We first integrated genome and transcriptome variation across 112 islet samples to produce dense cis-expression quantitative trait loci (cis-eQTL) maps. Additional integration with chromatin-state maps for islets and other diverse tissue types revealed that cis-eQTLs for islet-specific genes are specifically and significantly enriched in islet stretch enhancers. High-resolution chromatin accessibility profiling using assay for transposase-accessible chromatin sequencing (ATAC-seq) in two islet samples enabled us to identify specific transcription factor (TF) footprints embedded in active regulatory elements, which are highly enriched for islet cis-eQTL. Aggregate allelic bias signatures in TF footprints enabled us de novo to reconstruct TF binding affinities genetically, which support the high-quality nature of the TF footprint predictions. Interestingly, we found that T2D GWAS loci were strikingly and specifically enriched in islet Regulatory Factor X (RFX) footprints. Remarkably, within and across independent loci, T2D risk alleles that overlap with RFX footprints uniformly disrupt the RFX motifs at high-information content positions. Together, these results suggest that common regulatory variations have shaped islet TF footprints and the transcriptome and that a confluent RFX regulatory grammar plays a significant role in the genetic component of T2D predisposition.
Subject(s)
Diabetes Mellitus, Type 2/genetics , Genetic Predisposition to Disease , Genome, Human , Islets of Langerhans/metabolism , Quantitative Trait Loci , Transcriptome , Alleles , Base Sequence , Binding Sites , Chromatin/chemistry , Chromatin/metabolism , Diabetes Mellitus, Type 2/metabolism , Diabetes Mellitus, Type 2/pathology , Epigenesis, Genetic , Gene Expression Profiling , Genetic Variation , Genome-Wide Association Study , Genomic Imprinting , Humans , Islets of Langerhans/pathology , Polymorphism, Single Nucleotide , Protein Binding , Protein Isoforms/genetics , Protein Isoforms/metabolism , Regulatory Factor X Transcription Factors/genetics , Regulatory Factor X Transcription Factors/metabolismABSTRACT
BACKGROUND: Bisulfite sequencing is widely employed to study the role of DNA methylation in disease; however, the data suffer from biases due to coverage depth variability. Imputation of methylation values at low-coverage sites may mitigate these biases while also identifying important genomic features associated with predictive power. RESULTS: Here we describe BoostMe, a method for imputing low-quality DNA methylation estimates within whole-genome bisulfite sequencing (WGBS) data. BoostMe uses a gradient boosting algorithm, XGBoost, and leverages information from multiple samples for prediction. We find that BoostMe outperforms existing algorithms in speed and accuracy when applied to WGBS of human tissues. Furthermore, we show that imputation improves concordance between WGBS and the MethylationEPIC array at low WGBS depth, suggesting improved WGBS accuracy after imputation. CONCLUSIONS: Our findings support the use of BoostMe as a preprocessing step for WGBS analysis.
Subject(s)
Computational Biology/methods , DNA Methylation/drug effects , Sulfites/pharmacology , Whole Genome Sequencing , Algorithms , High-Throughput Nucleotide Sequencing , HumansABSTRACT
Hutchinson-Gilford progeria syndrome (HGPS) is a premature aging disease that is frequently caused by a de novo point mutation at position 1824 in LMNA. This mutation activates a cryptic splice donor site in exon 11, and leads to an in-frame deletion within the prelamin A mRNA and the production of a dominant-negative lamin A protein, known as progerin. Here we show that primary HGPS skin fibroblasts experience genome-wide correlated alterations in patterns of H3K27me3 deposition, DNA-lamin A/C associations, and, at late passages, genome-wide loss of spatial compartmentalization of active and inactive chromatin domains. We further demonstrate that the H3K27me3 changes associate with gene expression alterations in HGPS cells. Our results support a model that the accumulation of progerin in the nuclear lamina leads to altered H3K27me3 marks in heterochromatin, possibly through the down-regulation of EZH2, and disrupts heterochromatin-lamina interactions. These changes may result in transcriptional misregulation and eventually trigger the global loss of spatial chromatin compartmentalization in late passage HGPS fibroblasts.
Subject(s)
Genome, Human , Histones/metabolism , Lamins/metabolism , Progeria/genetics , Progeria/metabolism , Cell Line , Chromatin Immunoprecipitation , Fibroblasts/metabolism , Gene Expression Regulation , Heterochromatin/metabolism , Humans , Methylation , Protein Binding , Sequence Analysis, DNAABSTRACT
Chromatin-based functional genomic analyses and genomewide association studies (GWASs) together implicate enhancers as critical elements influencing gene expression and risk for common diseases. Here, we performed systematic chromatin and transcriptome profiling in human pancreatic islets. Integrated analysis of islet data with those from nine cell types identified specific and significant enrichment of type 2 diabetes and related quantitative trait GWAS variants in islet enhancers. Our integrated chromatin maps reveal that most enhancers are short (median = 0.8 kb). Each cell type also contains a substantial number of more extended (≥ 3 kb) enhancers. Interestingly, these stretch enhancers are often tissue-specific and overlap locus control regions, suggesting that they are important chromatin regulatory beacons. Indeed, we show that (i) tissue specificity of enhancers and nearby gene expression increase with enhancer length; (ii) neighborhoods containing stretch enhancers are enriched for important cell type-specific genes; and (iii) GWAS variants associated with traits relevant to a particular cell type are more enriched in stretch enhancers compared with short enhancers. Reporter constructs containing stretch enhancer sequences exhibited tissue-specific activity in cell culture experiments and in transgenic mice. These results suggest that stretch enhancers are critical chromatin elements for coordinating cell type-specific regulatory programs and that sequence variation in stretch enhancers affects risk of major common human diseases.
Subject(s)
Cell Differentiation/physiology , Chromatin/physiology , Diabetes Mellitus, Type 2/physiopathology , Enhancer Elements, Genetic/genetics , Epigenomics/methods , Gene Expression Regulation/physiology , Insulin-Secreting Cells/metabolism , Animals , Chromatin Immunoprecipitation , Diabetes Mellitus, Type 2/genetics , Enhancer Elements, Genetic/physiology , Gene Expression Profiling , Gene Expression Regulation/genetics , Genome-Wide Association Study , High-Throughput Nucleotide Sequencing , Humans , Insulin-Secreting Cells/physiology , Luciferases , Mice , Mice, TransgenicABSTRACT
Nuclear lamin B1 (LMNB1) constitutes one of the major structural proteins in the lamina mesh. We silenced the expression of LMNB1 by RNA interference in the colon cancer cell line DLD-1 and showed a dramatic redistribution of H3K27me3 from the periphery to a more homogeneous nuclear dispersion. In addition, we observed telomere attrition and an increased frequency of micronuclei and nuclear blebs. By 3D-FISH analyses, we demonstrated that the volume and surface of chromosome territories were significantly larger in LMNB1-depleted cells, suggesting that LMNB1 is required to maintain chromatin condensation in interphase nuclei. These changes led to a prolonged S phase due to activation of Chk1. Finally, silencing of LMNB1 resulted in extensive changes in alternative splicing of multiple genes and in a higher number of enlarged nuclear speckles. Taken together, our results suggest a mechanistic role of the nuclear lamina in the organization of chromosome territories, maintenance of genome integrity and proper gene splicing.
Subject(s)
Lamin Type B/physiology , S Phase/physiology , Adenocarcinoma/pathology , Cell Line, Tumor , Cell Nucleus/ultrastructure , Chromosomes, Human, Pair 18/ultrastructure , Chromosomes, Human, Pair 19/ultrastructure , Colonic Neoplasms/pathology , Heterochromatin/physiology , Histones/metabolism , Humans , In Situ Hybridization, Fluorescence , Interphase , Lamin Type B/deficiency , Methylation , Protein Processing, Post-Translational , RNA Interference , RNA, Small Interfering/pharmacology , S Phase/drug effects , Telomere Shortening/drug effectsABSTRACT
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.