RESUMO
Understanding the molecular mechanisms of complex traits is essential for developing targeted interventions. We analyzed liver expression quantitative-trait locus (eQTL) meta-analysis data on 1,183 participants to identify conditionally distinct signals. We found 9,013 eQTL signals for 6,564 genes; 23% of eGenes had two signals, and 6% had three or more signals. We then integrated the eQTL results with data from 29 cardiometabolic genome-wide association study (GWAS) traits and identified 1,582 GWAS-eQTL colocalizations for 747 eGenes. Non-primary eQTL signals accounted for 17% of all colocalizations. Isolating signals by conditional analysis prior to coloc resulted in 37% more colocalizations than using marginal eQTL and GWAS data, highlighting the importance of signal isolation. Isolating signals also led to stronger evidence of colocalization: among 343 eQTL-GWAS signal pairs in multi-signal regions, analyses that isolated the signals of interest resulted in higher posterior probability of colocalization for 41% of tests. Leveraging allelic heterogeneity, we predicted causal effects of gene expression on liver traits for four genes. To predict functional variants and regulatory elements, we colocalized eQTL with liver chromatin accessibility QTL (caQTL) and found 391 colocalizations, including 73 with non-primary eQTL signals and 60 eQTL signals that colocalized with both a caQTL and a GWAS signal. Finally, we used publicly available massively parallel reporter assays in HepG2 to highlight 14 eQTL signals that include at least one expression-modulating variant. This multi-faceted approach to unraveling the genetic underpinnings of liver-related traits could lead to therapeutic development.
RESUMO
Osteoarthritis (OA) poses a significant healthcare burden with limited treatment options. While genome-wide association studies (GWAS) have identified over 100 OA-associated loci, translating these findings into therapeutic targets remains challenging. Integrating expression quantitative trait loci (eQTL), 3D chromatin structure, and other genomic approaches with OA GWAS data offers a promising approach to elucidate disease mechanisms; however, comprehensive eQTL maps in OA-relevant tissues and conditions remain scarce. We mapped gene expression, chromatin accessibility, and 3D chromatin structure in primary human articular chondrocytes in both resting and OA-mimicking conditions. We identified thousands of differentially expressed genes, including those associated with differences in sex and age. RNA-seq in chondrocytes from 101 donors across two conditions uncovered 3782 unique eGenes, including 420 that exhibited strong and significant condition-specific effects. Colocalization with OA GWAS signals revealed 13 putative OA risk genes, 10 of which have not been previously identified. Chromatin accessibility and 3D chromatin structure provided insights into the mechanisms and conditional specificity of these variants. Our findings shed light on OA pathogenesis and highlight potential targets for therapeutic development. Highlights: ∘ Comprehensive analysis of sex- and age-related global gene expression in human chondrocytes revealed differences that correlate with osteoarthritis ∘ First response eQTLs in chondrocytes treated with an OA-related stimulus ∘ Deeply sequenced Hi-C in resting and activated chondrocytes helps connect OA risk variants to their putative causal genes ∘ Colocalization analysis reveals 13 (including 10 novel) putative OA risk genes.
RESUMO
Genome-wide association studies (GWAS) have identified numerous body mass index (BMI) loci. However, most underlying mechanisms from risk locus to BMI remain unknown. Leveraging omics data through integrative analyses could provide more comprehensive views of biological pathways on BMI. We analyzed genotype and blood gene expression data in up to 5,619 samples from the Framingham Heart Study (FHS). Using 3,992 single nucleotide polymorphisms (SNPs) at 97 BMI loci and 20,692 transcripts within 1 Mb, we performed separate association analyses of transcript with BMI and SNP with transcript (PBMI and PSNP, respectively) and then a correlated meta-analysis between the full summary data sets (PMETA). We identified transcripts that met Bonferroni-corrected significance for each omic, were more significant in the correlated meta-analysis than each omic, and were at least nominally associated with BMI in FHS data. Among 308 significant SNP-transcript-BMI associations, we identified seven genes (NT5C2, GSTM3, SNAPC3, SPNS1, TMEM245, YPEL3, and ZNF646) in five association regions. Using an independent sample of blood gene expression data, we validated results for SNAPC3 and YPEL3. We tested for generalization of these associations in hypothalamus, nucleus accumbens, and liver and observed significant (PMETA<0.05 & PMETA
RESUMO
Vertebrates transport hydrophobic triglycerides through the circulatory system by packaging them within amphipathic particles called Triglyceride-Rich Lipoproteins. Yet, it remains largely unknown how triglycerides are loaded onto these particles. Mutations in Phospholipase A2 group 12B (PLA2G12B) are known to disrupt lipoprotein homeostasis, but its mechanistic role in this process remains unclear. Here we report that PLA2G12B channels lipids within the lumen of the endoplasmic reticulum into nascent lipoproteins. This activity promotes efficient lipid secretion while preventing excess accumulation of intracellular lipids. We characterize the functional domains, subcellular localization, and interacting partners of PLA2G12B, demonstrating that PLA2G12B is calcium-dependent and tightly associated with the membrane of the endoplasmic reticulum. We also detect profound resistance to atherosclerosis in PLA2G12B mutant mice, suggesting an evolutionary tradeoff between triglyceride transport and cardiovascular disease risk. Here we identify PLA2G12B as a key driver of triglyceride incorporation into vertebrate lipoproteins.
Assuntos
Retículo Endoplasmático , Lipoproteínas , Animais , Camundongos , Transporte Biológico , Retículo Endoplasmático/metabolismo , Lipoproteínas/metabolismo , Triglicerídeos/metabolismoRESUMO
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.
Assuntos
Diabetes Mellitus Tipo 2 , Células-Tronco Pluripotentes Induzidas , Humanos , Repetições Palindrômicas Curtas Agrupadas e Regularmente Espaçadas/genética , Sistemas CRISPR-Cas/genética , Edição de Genes , RNA Guia de Sistemas CRISPR-CasRESUMO
Genome-wide association studies (GWASs) have identified hundreds of risk loci for liver disease and lipid-related metabolic traits, although identifying their target genes and molecular mechanisms remains challenging. We predicted target genes at GWAS signals by integrating them with molecular quantitative trait loci for liver gene expression (eQTL) and liver chromatin accessibility QTL (caQTL). We predicted specific regulatory caQTL variants at four GWAS signals located near EFHD1, LITAF, ZNF329, and GPR180. Using transcriptional reporter assays, we determined that caQTL variants rs13395911, rs11644920, rs34003091, and rs9556404 exhibit allelic differences in regulatory activity. We also performed a protein binding assay for rs13395911 and found that FOXA2 differentially interacts with the alleles of rs13395911. For variants rs13395911 and rs11644920 in putative enhancer regulatory elements, we used CRISPRi to demonstrate that repression of the enhancers altered the expression of the predicted target and/or nearby genes. Repression of the element at rs13395911 reduced the expression of EFHD1, and repression of the element at rs11644920 reduced the expression of LITAF, SNN, and TXNDC11. Finally, we showed that EFHD1 is a metabolically active gene in HepG2 cells. Together, these results provide key steps to connect genetic variants with cellular mechanisms and help elucidate the causes of liver disease.
Assuntos
Estudo de Associação Genômica Ampla , Hepatopatias , Humanos , Sequências Reguladoras de Ácido Nucleico , Lipídeos , Proteínas de TransporteRESUMO
BACKGROUND: Age and obesity are dominant risk factors for several common cardiometabolic disorders, and both are known to impair adipose tissue function. However, the underlying cellular and genetic factors linking aging and obesity on adipose tissue function have remained elusive. Adipose stem and precursor cells (ASPCs) are an understudied, yet crucial adipose cell type due to their deterministic adipocyte differentiation potential, which impacts the capacity to store fat in a metabolically healthy manner. METHODS: We integrated subcutaneous adipose tissue (SAT) bulk (n=435) and large single-nucleus RNA sequencing (n=105) data with the UK Biobank (UKB) (n=391,701) data to study age-obesity interactions originating from ASPCs by performing cell-type decomposition, differential expression testing, cell-cell communication analyses, and construction of polygenic risk scores for body mass index (BMI). RESULTS: We found that the SAT ASPC proportions significantly decrease with age in an obesity-dependent way consistently in two independent cohorts, both showing that the age dependency of ASPC proportions is abolished by obesity. We further identified 76 genes (72 SAT ASPC marker genes and 4 transcription factors regulating ASPC marker genes) that are differentially expressed by age in SAT and functionally enriched for developmental processes and adipocyte differentiation (i.e., adipogenesis). The 76 age-perturbed ASPC genes include multiple negative regulators of adipogenesis, such as RORA, SMAD3, TWIST2, and ZNF521, form tight clusters of longitudinally co-expressed genes during human adipogenesis, and show age-based differences in cellular interactions between ASPCs and adipose cell types. Finally, our genetic data demonstrate that cis-regional variants of these genes interact with age as predictors of BMI in an obesity-dependent way in the large UKB, while no such gene-age interaction on BMI is observed with non-age-dependent ASPC marker genes, thus independently confirming our cellular ASPC results at the biobank level. CONCLUSIONS: Overall, we discover that obesity prematurely induces a decrease in ASPC proportions and identify 76 developmentally important ASPC genes that implicate altered negative regulation of fat cell differentiation as a mechanism for aging and directly link aging to obesity via significant cellular and genetic interactions.
Assuntos
Tecido Adiposo , Obesidade , Humanos , Diferenciação Celular/genética , Obesidade/genética , Obesidade/metabolismo , Tecido Adiposo/metabolismo , Adipócitos/metabolismo , Envelhecimento/genética , Fatores de Transcrição/metabolismoRESUMO
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.
RESUMO
Conventional measurements of fasting and postprandial blood glucose levels investigated in genome-wide association studies (GWAS) cannot capture the effects of DNA variability on 'around the clock' glucoregulatory processes. Here we show that GWAS meta-analysis of glucose measurements under nonstandardized conditions (random glucose (RG)) in 476,326 individuals of diverse ancestries and without diabetes enables locus discovery and innovative pathophysiological observations. We discovered 120 RG loci represented by 150 distinct signals, including 13 with sex-dimorphic effects, two cross-ancestry and seven rare frequency signals. Of these, 44 loci are new for glycemic traits. Regulatory, glycosylation and metagenomic annotations highlight ileum and colon tissues, indicating an underappreciated role of the gastrointestinal tract in controlling blood glucose. Functional follow-up and molecular dynamics simulations of lower frequency coding variants in glucagon-like peptide-1 receptor (GLP1R), a type 2 diabetes treatment target, reveal that optimal selection of GLP-1R agonist therapy will benefit from tailored genetic stratification. We also provide evidence from Mendelian randomization that lung function is modulated by blood glucose and that pulmonary dysfunction is a diabetes complication. Our investigation yields new insights into the biology of glucose regulation, diabetes complications and pathways for treatment stratification.
Assuntos
Diabetes Mellitus Tipo 2 , Glucose , Humanos , Estudo de Associação Genômica Ampla , Glicemia/genética , Diabetes Mellitus Tipo 2/genética , ColoRESUMO
Understanding differences in adipose gene expression between individuals with different levels of clinical traits may reveal the genes and mechanisms leading to cardiometabolic diseases. However, adipose is a heterogeneous tissue. To account for cell-type heterogeneity, we estimated cell-type proportions in 859 subcutaneous adipose tissue samples with bulk RNA sequencing (RNA-seq) using a reference single-nuclear RNA-seq data set. Cell-type proportions were associated with cardiometabolic traits; for example, higher macrophage and adipocyte proportions were associated with higher and lower BMI, respectively. We evaluated cell-type proportions and BMI as covariates in tests of association between >25,000 gene expression levels and 22 cardiometabolic traits. For >95% of genes, the optimal, or best-fit, models included BMI as a covariate, and for 79% of associations, the optimal models also included cell type. After adjusting for the optimal covariates, we identified 2,664 significant associations (P ≤ 2e-6) for 1,252 genes and 14 traits. Among genes proposed to affect cardiometabolic traits based on colocalized genome-wide association study and adipose expression quantitative trait locus signals, 25 showed a corresponding association between trait and gene expression levels. Overall, these results suggest the importance of modeling cell-type proportion when identifying gene expression associations with cardiometabolic traits.
Assuntos
Doenças Cardiovasculares , Estudo de Associação Genômica Ampla , Humanos , Índice de Massa Corporal , Obesidade/genética , Expressão Gênica , Doenças Cardiovasculares/genéticaRESUMO
Importance: Hidradenitis suppurativa (HS) is a common and severely morbid chronic inflammatory skin disease that is reported to be highly heritable. However, the genetic understanding of HS is insufficient, and limited genome-wide association studies (GWASs) have been performed for HS, which have not identified significant risk loci. Objective: To identify genetic variants associated with HS and to shed light on the underlying genes and genetic mechanisms. Design, Setting, and Participants: This genetic association study recruited 753 patients with HS in the HS Program for Research and Care Excellence (HS ProCARE) at the University of North Carolina Department of Dermatology from August 2018 to July 2021. A GWAS was performed for 720 patients (after quality control) with controls from the Add Health study and then meta-analyzed with 2 large biobanks, UK Biobank (247 cases) and FinnGen (673 cases). Variants at 3 loci were tested for replication in the BioVU biobank (290 cases). Data analysis was performed from September 2021 to December 2022. Main Outcomes and Measures: Main outcome measures are loci identified, with association of P < 1 × 10-8 considered significant. Results: A total of 753 patients were recruited, with 720 included in the analysis. Mean (SD) age at symptom onset was 20.3 (10.57) years and at enrollment was 35.3 (13.52) years; 360 (50.0%) patients were Black, and 575 (79.7%) were female. In a meta-analysis of the 4 studies, 2 HS-associated loci were identified and replicated, with lead variants rs10512572 (P = 2.3 × 10-11) and rs17090189 (P = 2.1 × 10-8) near the SOX9 and KLF5 genes, respectively. Variants at these loci are located in enhancer regulatory elements detected in skin tissue. Conclusions and Relevance: In this genetic association study, common variants associated with HS located near the SOX9 and KLF5 genes were associated with risk of HS. These or other nearby genes may be associated with genetic risk of disease and the development of clinical features, such as cysts, comedones, and inflammatory tunnels, that are unique to HS. New insights into disease pathogenesis related to these genes may help predict disease progression and novel treatment approaches in the future.
Assuntos
Acne Vulgar , Hidradenite Supurativa , Humanos , Feminino , Masculino , Hidradenite Supurativa/genética , Hidradenite Supurativa/patologia , Estudo de Associação Genômica Ampla , Pele/patologia , Fatores de RiscoRESUMO
Gene regulatory effects in bulk-post mortem brain tissues are undetected at many non-coding brain trait-associated loci. We hypothesized that context-specific genetic variant function during stimulation of a developmental signaling pathway would explain additional regulatory mechanisms. We measured chromatin accessibility and gene expression following activation of the canonical Wnt pathway in primary human neural progenitors from 82 donors. TCF/LEF motifs, brain structure-, and neuropsychiatric disorder-associated variants were enriched within Wnt-responsive regulatory elements (REs). Genetically influenced REs were enriched in genomic regions under positive selection along the human lineage. Stimulation of the Wnt pathway increased the detection of genetically influenced REs/genes by 66.2%/52.7%, and led to the identification of 397 REs primed for effects on gene expression. Context-specific molecular quantitative trait loci increased brain-trait colocalizations by up to 70%, suggesting that genetic variant effects during early neurodevelopmental patterning lead to differences in adult brain and behavioral traits.
RESUMO
Cleft lip with or without cleft palate (CL/P) is a common birth defect with a complex, heterogeneous etiology. It is well established that common and rare sequence variants contribute to the formation of CL/P, but the contribution of copy-number variants (CNVs) to cleft formation remains relatively understudied. To fill this knowledge gap, we conducted a large-scale comparative analysis of genome-wide CNV profiles of 869 individuals from the Philippines and 233 individuals of European ancestry with CL/P with three primary goals: first, to evaluate whether differences in CNV number, amount of genomic content, or amount of coding genomic content existed within clefting subtypes; second, to assess whether CNVs in our cohort overlapped with known Mendelian clefting loci; and third, to identify unestablished Mendelian clefting genes. Significant differences in CNVs across cleft types or in individuals with non-syndromic versus syndromic clefts were not observed; however, several CNVs in our cohort overlapped with known syndromic and non-syndromic Mendelian clefting loci. Moreover, employing a filtering strategy relying on population genetics data that rare variants are on the whole more deleterious than common variants, we identify several CNV-associated gene losses likely driving non-syndromic clefting phenotypes. By prioritizing genes deleted at a rare frequency across multiple individuals with clefts yet enriched in our cohort of individuals with clefts compared to control subjects, we identify COBLL1, RIC1, and ARHGEF38 as clefting genes. CRISPR-Cas9 mutagenesis of these genes in Xenopus laevis and Danio rerio yielded craniofacial dysmorphologies, including clefts analogous to those seen in human clefting disorders.
Assuntos
Fenda Labial , Fissura Palatina , Variações do Número de Cópias de DNA , Humanos , Fenda Labial/genética , Fissura Palatina/genética , Estudo de Associação Genômica Ampla , Fatores de Troca do Nucleotídeo Guanina/genética , Fenótipo , Fatores de Transcrição/genéticaRESUMO
There is an increasing interest in using multiple types of omics features (e.g., DNA sequences, RNA expressions, methylation, protein expressions, and metabolic profiles) to study how the relationships between phenotypes and genotypes may be mediated by other omics markers. Genotypes and phenotypes are typically available for all subjects in genetic studies, but typically, some omics data will be missing for some subjects, due to limitations such as cost and sample quality. In this article, we propose a powerful approach for mediation analysis that accommodates missing data among multiple mediators and allows for various interaction effects. We formulate the relationships among genetic variants, other omics measurements, and phenotypes through linear regression models. We derive the joint likelihood for models with two mediators, accounting for arbitrary patterns of missing values. Utilizing computationally efficient and stable algorithms, we conduct maximum likelihood estimation. Our methods produce unbiased and statistically efficient estimators. We demonstrate the usefulness of our methods through simulation studies and an application to the Metabolic Syndrome in Men study.
Assuntos
Análise de Mediação , Modelos Genéticos , Humanos , Genótipo , Simulação por Computador , Funções Verossimilhança , AlgoritmosRESUMO
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.
RESUMO
Genome-wide association studies have identified over 100 loci associated with osteoarthritis risk, but the majority of osteoarthritis risk variants are noncoding, making it difficult to identify the impacted genes for further study and therapeutic development. To address this need, we used a multiomic approach and genome editing to identify and functionally characterize potential osteoarthritis risk genes. Computational analysis of genome-wide association studies and ChIP-seq data revealed that chondrocyte regulatory loci are enriched for osteoarthritis risk variants. We constructed a chondrocyte-specific regulatory network by mapping 3D chromatin structure and active enhancers in human chondrocytes. We then intersected these data with our previously collected RNA-seq dataset of chondrocytes responding to fibronectin fragment, a known osteoarthritis trigger. Integration of the 3 genomic datasets with recently reported osteoarthritis genome-wide association study variants revealed a refined set of putative causal osteoarthritis variants and their potential target genes. One of the putative target genes identified was SOCS2, which was connected to a putative causal variant by a 170-kb loop and is differentially regulated in response to fibronectin fragment. CRISPR-Cas9-mediated deletion of SOCS2 in primary human chondrocytes from 3 independent donors led to heightened expression of inflammatory markers after fibronectin fragment treatment. These data suggest that SOCS2 plays a role in resolving inflammation in response to cartilage matrix damage and provides a possible mechanistic explanation for its influence on osteoarthritis risk. In total, we identified 56 unique putative osteoarthritis risk genes for further research and potential therapeutic development.
Assuntos
Condrócitos , Osteoartrite , Humanos , Fibronectinas/genética , Fibronectinas/metabolismo , Estudo de Associação Genômica Ampla , Osteoartrite/genética , Osteoartrite/metabolismo , Cromatina/genética , Cromatina/metabolismoRESUMO
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.
Assuntos
Estudo de Associação Genômica Ampla , Transcriptoma , Bilirrubina , Carnitina , Glicerofosfolipídeos , Humanos , Masculino , Metabolômica , Locos de Características Quantitativas/genética , Membro 5 da Família 22 de Carreadores de Soluto/genética , Transcriptoma/genéticaRESUMO
Obesity perturbs central functions of human adipose tissue, centred on differentiation of preadipocytes to adipocytes, i.e., adipogenesis. The large environmental component of obesity makes it important to elucidate epigenetic regulatory factors impacting adipogenesis. Promoter Capture Hi-C (pCHi-C) has been used to identify chromosomal interactions between promoters and associated regulatory elements. However, long range interactions (LRIs) greater than 1 Mb are often filtered out of pCHi-C datasets, due to technical challenges and their low prevalence. To elucidate the unknown role of LRIs in adipogenesis, we investigated preadipocyte differentiation to adipocytes using pCHi-C and bulk and single nucleus RNA-seq data. We first show that LRIs are reproducible between biological replicates, and they increase >2-fold in frequency across adipogenesis. We further demonstrate that genomic loci containing LRIs are more epigenetically repressed than regions without LRIs, corresponding to lower gene expression in the LRI regions. Accordingly, as preadipocytes differentiate into adipocytes, LRI regions are more likely to contain repressed preadipocyte marker genes; whereas these same LRI regions are depleted of actively expressed adipocyte marker genes. Finally, we show that LRIs can be used to restrict multiple testing of the long-range cis-eQTL analysis to identify variants that regulate genes via LRIs. We exemplify this by identifying a putative long range cis regulatory mechanism at the LYPLAL1/TGFB2 obesity locus. In summary, we identify LRIs that mark repressed regions of the genome, and these interactions increase across adipogenesis, pinpointing developmental regions that need to be repressed in a cell-type specific way for adipogenesis to proceed.
Assuntos
Adipogenia , Metilação de DNA , Camundongos , Animais , Humanos , Adipogenia/genética , Células 3T3-L1 , Diferenciação Celular/genética , Obesidade/genética , Obesidade/metabolismo , Expressão GênicaRESUMO
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.