Colocalization of GWAS and eQTL signals at loci with multiple signals identifies additional candidate genes for body fat distribution

Druggable genome-wide Mendelian randomization for identifying the role of integrated stress response in therapeutic targets of bipolar disorder

For bipolar disorder (BD), the inconsistency of treatment guidelines and the long phases of pharmacological adjustment remain major challenges. BD is known to be comorbid with many medical and psychiatric conditions and they may share inflammatory and stress-related aetiologies, which could give rise to this association. The integrated stress response (ISR) responds to various stress conditions that lead to alterations in cellular homeostasis. However, as a causative mechanism underlying cognitive deficits and neurodegeneration in a broad range of brain disorders, the impact of ISR on BD is understudied. Mendelian randomization has been widely used to repurpose licensed drugs and discover novel therapeutic targets. Thus, we aimed to identify novel therapeutic targets for BD and analyze their pathophysiological mechanisms, using the summary data-based Mendelian Randomization (SMR) and Bayesian colocalization (COLOC) methods to integrate the summary-level data of the GWAS on BD and the expression quantitative trait locus (eQTL) study in blood. We utilized the GWAS data including 41,917 BD cases and 371,549 controls from the Psychiatric Genomics Consortium and the eQTL data from 31,684 participants of predominantly European ancestry from the eQTLGen consortium. The SMR analysis identified the EIF2B5 gene that was associated with BD due to no linkage but pleiotropy or causality. The COLOC analysis strongly suggested that EIF2B5 and the trait of BD were affected by shared causal variants, and thus were colocalized. Utilizing data in EpiGraphDB we find other putative causal BD genes (EIF2AK4 and GSK3B) to prioritize potential alternative drug targets.

Single nucleus RNA-sequencing integrated into risk variant colocalization discovers 17 cell-type-specific abdominal obesity genes for metabolic dysfunction-associated steatotic liver disease

BACKGROUND

Abdominal obesity increases the risk for non-alcoholic fatty liver disease (NAFLD), now known as metabolic dysfunction-associated steatotic liver disease (MASLD).

METHODS

To elucidate the directional cell-type level biological mechanisms underlying the association between abdominal obesity and MASLD, we integrated adipose and liver single nucleus RNA-sequencing and bulk cis-expression quantitative trait locus (eQTL) data with the UK Biobank genome-wide association study (GWAS) data using colocalization. Then we used colocalized cis-eQTL variants as instrumental variables in Mendelian randomization (MR) analyses, followed by functional validation experiments on the target genes of the cis-eQTL variants.

FINDINGS

We identified 17 colocalized abdominal obesity GWAS variants, regulating 17 adipose cell-type marker genes. Incorporating these 17 variants into MR discovers a putative tissue-of-origin, cell-type-aware causal effect of abdominal obesity on MASLD consistently with multiple MR methods without significant evidence for pleiotropy or heterogeneity. Single cell data confirm the adipocyte-enriched mean expression of the 17 genes. Our cellular experiments across human adipogenesis identify risk variant -specific epigenetic and transcriptional mechanisms. Knocking down two of the 17 genes, PPP2R5A and SH3PXD2B, shows a marked decrease in adipocyte lipidation and significantly alters adipocyte function and adipogenesis regulator genes, including DGAT2, LPL, ADIPOQ, PPARG, and SREBF1. Furthermore, the 17 genes capture a characteristic MASLD expression signature in subcutaneous adipose tissue.

INTERPRETATION

Overall, we discover a significant cell-type level effect of abdominal obesity on MASLD and trace its biological effect to adipogenesis.

FUNDING

NIH grants R01HG010505, R01DK132775, and R01HL170604; the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (Grant No. 802825), Academy of Finland (Grants Nos. 333021), the Finnish Foundation for Cardiovascular Research the Sigrid Jusélius Foundation and the Jane and Aatos Erkko Foundation; American Association for the Study of Liver Diseases (AASLD) Advanced Transplant Hepatology award and NIH/NIDDK (P30DK41301) Pilot and Feasibility award; NIH/NIEHS F32 award (F32ES034668); Finnish Diabetes Research Foundation, Kuopio University Hospital Project grant (EVO/VTR grants 2005-2021), the Academy of Finland grant (Contract no. 138006); Academy of Finland (Grant Nos 335443, 314383, 272376 and 266286), Sigrid Jusélius Foundation, Finnish Medical Foundation, Finnish Diabetes Research Foundation, Novo Nordisk Foundation (#NNF20OC0060547, NNF17OC0027232, NNF10OC1013354) and Government Research Funds to Helsinki University Hospital; Orion Research Foundation, Maud Kuistila Foundation, Finish Medical Foundation, and University of Helsinki.

Systems genetics analysis of human body fat distribution genes identifies adipocyte processes

Excess abdominal fat is a sexually dimorphic risk factor for cardio-metabolic disease and is approximated by the waist-to-hip ratio adjusted for body mass index (WHRadjBMI). Whereas this trait is highly heritable, few causal genes are known. We aimed to identify novel drivers of WHRadjBMI using systems genetics. We used two independent cohorts of adipose tissue gene expression and constructed sex- and depot-specific Bayesian networks to model gene-gene interactions from 8,492 genes. Using key driver analysis, we identified genes that, in silico and putatively in vitro, regulate many others. 51-119 key drivers in each network were replicated in both cohorts. In other cell types, 23 of these genes are found in crucial adipocyte pathways: Wnt signaling or mitochondrial function. We overexpressed or down-regulated seven key driver genes in human subcutaneous pre-adipocytes. Key driver genes ANAPC2 and RSPO1 inhibited adipogenesis, whereas PSME3 increased adipogenesis. RSPO1 increased Wnt signaling activity. In differentiated adipocytes, MIGA1 and UBR1 down-regulation led to mitochondrial dysfunction. These five genes regulate adipocyte function, and we hypothesize that they regulate fat distribution.

Mendelian randomization reveals apolipoprotein B shortens healthspan and possibly increases risk for Alzheimer's disease

Apolipoprotein B-100 (APOB) is a component of fat- and cholesterol-transporting molecules in the bloodstream. It is the main lipoprotein in low-density lipoprotein cholesterol (LDL) and has been implicated in conditions that end healthspan (the interval between birth and onset of chronic disease). However, APOB's direct relationship with healthspan remains uncertain. With Mendelian randomization, we show that higher levels of APOB and LDL shorten healthspan in humans. Multivariable Mendelian randomization of APOB and LDL on healthspan suggests that the predominant trait accounting for the relationship is APOB. In addition, we provide preliminary evidence that APOB increases risk for Alzheimer's disease, a condition that ends healthspan. If these relationships are causal, they suggest that interventions to improve healthspan in aging populations could include strategies targeting APOB. Ultimately, given that more than 44 million people currently suffer from Alzheimer's disease worldwide, such interventions are needed.

Molecular quantitative trait loci in reproductive tissues impact male fertility in cattle

Breeding bulls are well suited to investigate inherited variation in male fertility because they are genotyped and their reproductive success is monitored through semen analyses and thousands of artificial inseminations. However, functional data from relevant tissues are lacking in cattle, which prevents fine-mapping fertility-associated genomic regions. Here, we characterize gene expression and splicing variation in testis, epididymis, and vas deferens transcriptomes of 118 mature bulls and conduct association tests between 414,667 molecular phenotypes and 21,501,032 genome-wide variants to identify 41,156 regulatory loci. We show broad consensus in tissue-specific and tissue-enriched gene expression between the three bovine tissues and their human and murine counterparts. Expression- and splicing-mediating variants are more than three times as frequent in testis than epididymis and vas deferens, highlighting the transcriptional complexity of testis. Finally, we identify genes (WDR19, SPATA16, KCTD19, ZDHHC1) and molecular phenotypes that are associated with quantitative variation in male fertility through transcriptome-wide association and colocalization analyses.

Adipose tissue eQTL meta-analysis reveals the contribution of allelic heterogeneity to gene expression regulation and cardiometabolic traits

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.

Genetics and epigenetics in the obesity phenotyping scenario

Obesity is a common complex trait that elevates the risk for various diseases, including type 2 diabetes and cardiovascular disease. A combination of environmental and genetic factors influences the pathogenesis of obesity. Advances in genomic technologies have driven the identification of multiple genetic loci associated with this disease, ranging from studying severe onset cases to investigating common multifactorial polygenic forms. Additionally, findings from epigenetic analyses of modifications to the genome that do not involve changes to the underlying DNA sequence have emerged as key signatures in the development of obesity. Such modifications can mediate the effects of environmental factors, including diet and lifestyle, on gene expression and clinical presentation. This review outlines what is known about the genetic and epigenetic contributors to obesity susceptibility, along with the albeit limited therapeutic options currently available. Furthermore, we delineate the potential mechanisms of actions through which epigenetic changes can mediate environmental influences and the related opportunities they present for future interventions in the management of obesity.

Assessing efficiency of fine-mapping obesity-associated variants through leveraging ancestry architecture and functional annotation using PAGE and UKBB cohorts

Inadequate representation of non-European ancestry populations in genome-wide association studies (GWAS) has limited opportunities to isolate functional variants. Fine-mapping in multi-ancestry populations should improve the efficiency of prioritizing variants for functional interrogation. To evaluate this hypothesis, we leveraged ancestry architecture to perform comparative GWAS and fine-mapping of obesity-related phenotypes in European ancestry populations from the UK Biobank (UKBB) and multi-ancestry samples from the Population Architecture for Genetic Epidemiology (PAGE) consortium with comparable sample sizes. In the investigated regions with genome-wide significant associations for obesity-related traits, fine-mapping in our ancestrally diverse sample led to 95% and 99% credible sets (CS) with fewer variants than in the European ancestry sample. Lead fine-mapped variants in PAGE regions had higher average coding scores, and higher average posterior probabilities for causality compared to UKBB. Importantly, 99% CS in PAGE loci contained strong expression quantitative trait loci (eQTLs) in adipose tissues or harbored more variants in tighter linkage disequilibrium (LD) with eQTLs. Leveraging ancestrally diverse populations with heterogeneous ancestry architectures, coupled with functional annotation, increased fine-mapping efficiency and performance, and reduced the set of candidate variants for consideration for future functional studies. Significant overlap in genetic causal variants across populations suggests generalizability of genetic mechanisms underpinning obesity-related traits across populations.

Systems genetics analysis of human body fat distribution genes identifies Wnt signaling and mitochondrial activity in adipocytes

BACKGROUND

Excess fat in the abdomen is a sexually dimorphic risk factor for cardio-metabolic disease. The relative storage between abdominal and lower-body subcutaneous adipose tissue depots is approximated by the waist-to-hip ratio adjusted for body mass index (WHRadjBMI). Genome-wide association studies (GWAS) identified 346 loci near 495 genes associated with WHRadjBMI. Most of these genes have unknown roles in fat distribution, but many are expressed and putatively act in adipose tissue. We aimed to identify novel sex- and depot-specific drivers of WHRadjBMI using a systems genetics approach.

METHODS

We used two independent cohorts of adipose tissue gene expression with 362 - 444 males and 147 - 219 females, primarily of European ancestry. We constructed sex- and depot- specific Bayesian networks to model the gene-gene interactions from 8,492 adipose tissue genes. Key driver analysis identified genes that, in silico and putatively in vitro, regulate many others, including the 495 WHRadjBMI GWAS genes. Key driver gene function was determined by perturbing their expression in human subcutaneous pre-adipocytes using lenti-virus or siRNA.

RESULTS

51 - 119 key drivers in each network were replicated in both cohorts. We used single-cell expression data to select replicated key drivers expressed in adipocyte precursors and mature adipocytes, prioritized genes which have not been previously studied in adipose tissue, and used public human and mouse data to nominate 53 novel key driver genes (10 - 21 from each network) that may regulate fat distribution by altering adipocyte function. In other cell types, 23 of these genes are found in crucial adipocyte pathways: Wnt signaling or mitochondrial function. We selected seven genes whose expression is highly correlated with WHRadjBMI to further study their effects on adipogenesis/Wnt signaling (ANAPC2, PSME3, RSPO1, TYRO3) or mitochondrial function (C1QTNF3, MIGA1, PSME3, UBR1).Adipogenesis was inhibited in cells overexpressing ANAPC2 and RSPO1 compared to controls. RSPO1 results are consistent with a positive correlation between gene expression in the subcutaneous depot and WHRadjBMI, therefore lower relative storage in the subcutaneous depot. RSPO1 inhibited adipogenesis by increasing β-catenin activation and Wnt-related transcription, thus repressing PPARG and CEBPA. PSME3 overexpression led to more adipogenesis than controls. In differentiated adipocytes, MIGA1 and UBR1 downregulation led to mitochondrial dysfunction, with lower oxygen consumption than controls; MIGA1 knockdown also lowered UCP1 expression.

SUMMARY

ANAPC2, MIGA1, PSME3, RSPO1, and UBR1 affect adipocyte function and may drive body fat distribution.

Identifying novel regulatory effects for clinically relevant genes through the study of the Greek population

BACKGROUND

Expression quantitative trait loci (eQTL) studies provide insights into regulatory mechanisms underlying disease risk. Expanding studies of gene regulation to underexplored populations and to medically relevant tissues offers potential to reveal yet unknown regulatory variants and to better understand disease mechanisms. Here, we performed eQTL mapping in subcutaneous (S) and visceral (V) adipose tissue from 106 Greek individuals (Greek Metabolic study, GM) and compared our findings to those from the Genotype-Tissue Expression (GTEx) resource.

RESULTS

We identified 1,930 and 1,515 eGenes in S and V respectively, over 13% of which are not observed in GTEx adipose tissue, and that do not arise due to different ancestry. We report additional context-specific regulatory effects in genes of clinical interest (e.g. oncogene ST7) and in genes regulating responses to environmental stimuli (e.g. MIR21, SNX33). We suggest that a fraction of the reported differences across populations is due to environmental effects on gene expression, driving context-specific eQTLs, and suggest that environmental effects can determine the penetrance of disease variants thus shaping disease risk. We report that over half of GM eQTLs colocalize with GWAS SNPs and of these colocalizations 41% are not detected in GTEx. We also highlight the clinical relevance of S adipose tissue by revealing that inflammatory processes are upregulated in individuals with obesity, not only in V, but also in S tissue.

CONCLUSIONS

By focusing on an understudied population, our results provide further candidate genes for investigation regarding their role in adipose tissue biology and their contribution to disease risk and pathogenesis.

Functional characterization of human genomic variation linked to polygenic diseases

The burden of human disease lies predominantly in polygenic diseases. Since the early 2000s, genome-wide association studies (GWAS) have identified genetic variants and loci associated with complex traits. These have ranged from variants in coding sequences to mutations in regulatory regions, such as promoters and enhancers, as well as mutations affecting mediators of mRNA stability and other downstream regulators, such as 5' and 3'-untranslated regions (UTRs), long noncoding RNA (lncRNA), and miRNA. Recent research advances in genetics have utilized a combination of computational techniques, high-throughput in vitro and in vivo screening modalities, and precise genome editing to impute the function of diverse classes of genetic variants identified through GWAS. In this review, we highlight the vastness of genomic variants associated with polygenic disease risk and address recent advances in how genetic tools can be used to functionally characterize them.

Small open reading frames: a comparative genetics approach to validation

Open reading frames (ORFs) with fewer than 100 codons are generally not annotated in genomes, although bona fide genes of that size are known. Newer biochemical studies have suggested that thousands of small protein-coding ORFs (smORFs) may exist in the human genome, but the true number and the biological significance of the micropeptides they encode remain uncertain. Here, we used a comparative genomics approach to identify high-confidence smORFs that are likely protein-coding. We identified 3,326 high-confidence smORFs using constraint within human populations and evolutionary conservation as additional lines of evidence. Next, we validated that, as a group, our high-confidence smORFs are conserved at the amino-acid level rather than merely residing in highly conserved non-coding regions. Finally, we found that high-confidence smORFs are enriched among disease-associated variants from GWAS. Overall, our results highlight that smORF-encoded peptides likely have important functional roles in human disease.

Genetic Regulation of SMC Gene Expression and Splicing Predict Causal CAD Genes

BACKGROUND

Coronary artery disease (CAD) is the leading cause of death worldwide. Recent meta-analyses of genome-wide association studies have identified over 175 loci associated with CAD. The majority of these loci are in noncoding regions and are predicted to regulate gene expression. Given that vascular smooth muscle cells (SMCs) play critical roles in the development and progression of CAD, we aimed to identify the subset of the CAD loci associated with the regulation of transcription in distinct SMC phenotypes.

METHODS

We measured gene expression in SMCs isolated from the ascending aortas of 151 heart transplant donors of various genetic ancestries in quiescent or proliferative conditions and calculated the association of their expression and splicing with ~6.3 million imputed single-nucleotide polymorphism markers across the genome.

RESULTS

We identified 4910 expression and 4412 splicing quantitative trait loci (sQTLs) representing regions of the genome associated with transcript abundance and splicing. A total of 3660 expression quantitative trait loci (eQTLs) had not been observed in the publicly available Genotype-Tissue Expression dataset. Further, 29 and 880 eQTLs were SMC-specific and sex-biased, respectively. We made these results available for public query on a user-friendly website. To identify the effector transcript(s) regulated by CAD loci, we used 4 distinct colocalization approaches. We identified 84 eQTL and 164 sQTL that colocalized with CAD loci, highlighting the importance of genetic regulation of mRNA splicing as a molecular mechanism for CAD genetic risk. Notably, 20% and 35% of the eQTLs were unique to quiescent or proliferative SMCs, respectively. One CAD locus colocalized with a sex-specific eQTL (TERF2IP), and another locus colocalized with SMC-specific eQTL (ALKBH8). The most significantly associated CAD locus, 9p21, was an sQTL for the long noncoding RNA CDKN2B-AS1, also known as ANRIL, in proliferative SMCs.

CONCLUSIONS

Collectively, our results provide evidence for the molecular mechanisms of genetic susceptibility to CAD in distinct SMC phenotypes.

Effect of all-but-one conditional analysis for eQTL isolation in peripheral blood

Expression quantitative trait locus detection has become increasingly important for understanding how noncoding variants contribute to disease susceptibility and complex traits. The major challenges in expression quantitative trait locus fine-mapping and causal variant discovery relate to the impact of linkage disequilibrium on signals due to one or multiple functional variants that lie within a credible set. We perform expression quantitative trait locus fine-mapping using the all-but-one approach, conditioning each signal on all others detected in an interval, on the Consortium for the Architecture of Gene Expression cohorts of microarray-based peripheral blood gene expression in 2,138 European-ancestry human adults. We contrast these results with traditional forward stepwise conditional analysis and a Bayesian localization method. All-but-one conditioning significantly modifies effect-size estimates for 51% of 2,351 expression quantitative trait locus peaks, but only modestly affects credible set size and location. On the other hand, both conditioning approaches result in unexpectedly low overlap with Bayesian credible sets, with just 57% peak concordance and between 50% and 70% SNP sharing, leading us to caution against the assumption that any one localization method is superior to another. We also cross reference our results with ATAC-seq data, cell-type-specific expression quantitative trait locus, and activity-by-contact-enhancers, leading to the proposal of a 5-tier approach to further reduce credible set sizes and prioritize likely causal variants for all known inflammatory bowel disease risk loci active in immune cells.

Cellular Genome-wide Association Study Identifies Common Genetic Variation Influencing Lithium-Induced Neural Progenitor Proliferation

BACKGROUND

Bipolar disorder is a highly heritable neuropsychiatric condition affecting more than 1% of the human population. Lithium salts are commonly prescribed as a mood stabilizer for individuals with bipolar disorder. Lithium is clinically effective in approximately half of treated individuals, and their genetic backgrounds are known to influence treatment outcomes. While the mechanism of lithium's therapeutic action is unclear, it stimulates adult neural progenitor cell proliferation, similar to some antidepressant drugs.

METHODS

To identify common genetic variants that modulate lithium-induced proliferation, we conducted an EdU incorporation assay in a library of 80 genotyped human neural progenitor cells treated with lithium. These data were used to perform a genome-wide association study to identify common genetic variants that influence lithium-induced neural progenitor cell proliferation. We manipulated the expression of a putatively causal gene using CRISPRi/a (clustered regularly interspaced short palindromic repeats interference/activation) constructs to experimentally verify lithium-induced proliferation effects.

RESULTS

We identified a locus on chr3p21.1 associated with lithium-induced proliferation. This locus is also associated with bipolar disorder risk, schizophrenia risk, and interindividual differences in intelligence. We identified a single gene, GNL3, whose expression temporally increased in an allele-specific fashion following lithium treatment. Experimentally increasing the expression of GNL3 led to increased proliferation under baseline conditions, while experimentally decreasing GNL3 expression suppressed lithium-induced proliferation.

CONCLUSIONS

Our experiments reveal that common genetic variation modulates lithium-induced neural progenitor proliferation and that GNL3 expression is necessary for the full proliferation-stimulating effects of lithium. These results suggest that performing genome-wide associations in genetically diverse human cell lines is a useful approach to discover context-specific pharmacogenomic effects.

Cis-eQTL Analysis and Functional Validation of Candidate Genes for Carcass Yield Traits in Beef Cattle

Carcass yield traits are of considerable economic importance for farm animals, which act as a major contributor to the world’s food supply. Genome-wide association studies (GWASs) have identified many genetic variants associated with carcass yield traits in beef cattle. However, their functions are not effectively illustrated. In this study, we performed an integrative analysis of gene-based GWAS with expression quantitative trait locus (eQTL) analysis to detect candidate genes for carcass yield traits and validate their effects on bovine skeletal muscle satellite cells (BSCs). The gene-based GWAS and cis-eQTL analysis revealed 1780 GWAS and 1538 cis-expression genes. Among them, we identified 153 shared genes that may play important roles in carcass yield traits. Notably, the identified cis-eQTLs of PON3 and PRIM2 were significantly (p < 0.001) enriched in previous GWAS loci for carcass traits. Furthermore, overexpression of PON3 and PRIM2 promoted the BSCs’ proliferation, increased the expression of MYOD and downregulated the expression of MYOG, which indicated that these genes may inhibit myogenic differentiation. In contrast, PON3 and PRIM2 were significantly downregulated during the differentiation of BSCs. These findings suggested that PON3 and PRIM2 may promote the proliferation of BSCs and inhibit them in the pre-differentiation stage. Our results further contribute to the understanding of the molecular mechanisms of carcass yield traits in beef cattle.

Association between telomere length and insomnia: A mendelian randomization and colocalization study

BACKGROUND

Previous studies have suggested a potential association between sleep and telomere length (TL), but its genetic basis remains unclear. In this study, we aimed to explore the genetic correlation and potential causal association between TL and insomnia.

METHODS

The genome-wide association study (GWAS) datasets of TL and insomnia-related traits were used, including insomnia, snoring, daytime dozing and napping. Based on the polygenic risk scores (PRS) of TL, linear regression and linkage disequilibrium score (LDSC) regression were used to preliminarily explore the association between TL and insomnia parameters in the UK Biobank cohort. Then, we investigated the causal association between TL and insomnia by mendelian randomization (MR) analysis and colocalization analysis.

RESULTS

In the UK Biobank cohort, the association between TL and insomnia was observed in the female samples (t = 2.968, P = 3.00 × 10^-3). LDSC detected a genetic correlation between short TL and insomnia (Rg = -9.27 × 10^-2, P = 8.00 × 10^-4). We found no evidence supporting significant causal association between insomnia and TL in IVW method (b = -5.95 × 10^-3, P = 0.57), with horizontal pleiotropy and heterogeneity tests indicating the validity of our MR study. Finally, rs12638862 was classified as colocalized by COLOC (PP4 = 0.99), and TERC may be involved in regulating the association between insomnia and TL.

CONCLUSIONS

Our study found no evidence for causal association between insomnia and TL in individuals of European ancestry. We detected a candidate gene associated with both insomnia and TL, providing novel clues for understanding the roles of this association.

Colocalization of Gene Expression and DNA Methylation with Genetic Risk Variants Supports Functional Roles of MUC5B and DSP in Idiopathic Pulmonary Fibrosis

Rationale: Common genetic variants have been associated with idiopathic pulmonary fibrosis (IPF). Objectives: To determine functional relevance of the 10 IPF-associated common genetic variants we previously identified. Methods: We performed expression quantitative trait loci (eQTL) and methylation quantitative trait loci (mQTL) mapping, followed by co-localization of eQTL and mQTL with genetic association signals and functional validation by luciferase reporter assays. Illumina multi-ethnic genotyping arrays, mRNA sequencing, and Illumina 850k methylation arrays were performed on lung tissue of participants with IPF (234 RNA and 345 DNA samples) and non-diseased controls (188 RNA and 202 DNA samples). Measurements and Main Results: Focusing on genetic variants within 10 IPF-associated genetic loci, we identified 27 eQTLs in controls and 24 eQTLs in cases (false-discovery-rate-adjusted P < 0.05). Among these signals, we identified associations of lead variants rs35705950 with expression of MUC5B and rs2076295 with expression of DSP in both cases and controls. mQTL analysis identified CpGs in gene bodies of MUC5B (cg17589883) and DSP (cg08964675) associated with the lead variants in these two loci. We also demonstrated strong co-localization of eQTL/mQTL and genetic signal in MUC5B (rs35705950) and DSP (rs2076295). Functional validation of the mQTL in MUC5B using luciferase reporter assays demonstrates that the CpG resides within a putative internal repressor element. Conclusions: We have established a relationship of the common IPF genetic risk variants rs35705950 and rs2076295 with respective changes in MUC5B and DSP expression and methylation. These results provide additional evidence that both MUC5B and DSP are involved in the etiology of IPF.

Identification of genetic loci simultaneously associated with multiple cardiometabolic traits

BACKGROUND AND AIMS

Cardiometabolic disorders (CMD) arise from a constellation of features such as increased adiposity, hyperlipidemia, hypertension and compromised glucose control. Many genetic loci have shown associations with individual CMD-related traits, but no investigations have focused on simultaneously identifying loci showing associations across all domains. We therefore sought to identify loci associated with risk across seven continuous CMD-related traits.

METHODS AND RESULTS

We conducted separate genome-wide association studies (GWAS) for systolic and diastolic blood pressure (SBP/DBP), hemoglobin A1c (HbA1c), low- and high- density lipoprotein cholesterol (LDL-C/HDL-C), waist-to-hip-ratio (WHR), and triglycerides (TGs) in the UK Biobank (N = 356,574-456,823). Multiple loci reached genome-wide levels of significance (N = 145-333) for each trait, but only four loci (in/near VEGFA, GRB14-COBLL1, KLF14, and RGS19-OPRL1) were associated with risk across all seven traits (P < 5 × 10^-8). We sought replication of these four loci in an independent set of seven trait-specific GWAS meta-analyses. GRB14-COBLL1 showed the most consistent replication, revealing nominally significant associations (P < 0.05) with all traits except DBP.

CONCLUSIONS

Our analyses suggest that very few loci are associated in the same direction of risk with traits representing the full spectrum of CMD features. We identified four such loci, and an understanding of the pathways between these loci and CMD risk may eventually identify factors that can be used to identify pathologic disturbances that represent broadly beneficial therapeutic targets.

Discovery of novel eGFR-associated multiple independent signals using a quasi-adaptive method

A decreased estimated glomerular filtration rate (eGFR) leading to chronic kidney disease is a significant public health problem. Kidney function is a heritable trait, and recent application of genome-wide association studies (GWAS) successfully identified multiple eGFR-associated genetic loci. To increase statistical power for detecting independent associations in GWAS loci, we improved our recently developed quasi-adaptive method estimating SNP-specific alpha levels for the conditional analysis, and applied it to the GWAS meta-analysis results of eGFR among 783,978 European-ancestry individuals. Among known eGFR loci, we revealed 19 new independent association signals that were subsequently replicated in the United Kingdom Biobank (n = 408,608). These associations have remained undetected by conditional analysis using the established conservative genome-wide significance level of 5 × 10^-8. Functional characterization of known index SNPs and novel independent signals using colocalization of conditional eGFR association results and gene expression in cis across 51 human tissues identified two potentially causal genes across kidney tissues: TSPAN33 and TFDP2, and three candidate genes across other tissues: SLC22A2, LRP2, and CDKN1C. These colocalizations were not identified in the original GWAS. By applying our improved quasi-adaptive method, we successfully identified additional genetic variants associated with eGFR. Considering these signals in colocalization analyses can increase the precision of revealing potentially functional genes of GWAS loci.

Multi-Omic Approaches to Identify Genetic Factors in Metabolic Syndrome

Metabolic syndrome (MetS) is a highly heritable disease and a major public health burden worldwide. MetS diagnosis criteria are met by the simultaneous presence of any three of the following: high triglycerides, low HDL/high LDL cholesterol, insulin resistance, hypertension, and central obesity. These diseases act synergistically in people suffering from MetS and dramatically increase risk of morbidity and mortality due to stroke and cardiovascular disease, as well as certain cancers. Each of these component features is itself a complex disease, as is MetS. As a genetically complex disease, genetic risk factors for MetS are numerous, but not very powerful individually, often requiring specific environmental stressors for the disease to manifest. When taken together, all sequence variants that contribute to MetS disease risk explain only a fraction of the heritable variance, suggesting additional, novel loci have yet to be discovered. In this article, we will give a brief overview on the genetic concepts needed to interpret genome-wide association studies (GWAS) and quantitative trait locus (QTL) data, summarize the state of the field of MetS physiological genomics, and to introduce tools and resources that can be used by the physiologist to integrate genomics into their own research on MetS and any of its component features. There is a wealth of phenotypic and molecular data in animal models and humans that can be leveraged as outlined in this article. Integrating these multi-omic QTL data for complex diseases such as MetS provides a means to unravel the pathways and mechanisms leading to complex disease and promise for novel treatments. © 2022 American Physiological Society. Compr Physiol 12:1-40, 2022.

Integrating genome-wide association and transcriptome prediction model identifies novel target genes for osteoporosis

In this study, we integrated large-scale GWAS summary data and used the predicted transcriptome-wide association study method to discover novel genes associated with osteoporosis. We identified 204 candidate genes, which provide novel clues for understanding the genetic mechanism of osteoporosis and indicate potential therapeutic targets.

INTRODUCTION

Osteoporosis is a highly polygenetic disease characterized by low bone mass and deterioration of the bone microarchitecture. Our objective was to discover novel candidate genes associated with osteoporosis.

METHODS

To identify potential causal genes of the associated loci, we investigated trait-gene expression associations using the transcriptome-wide association study (TWAS) method. This method directly imputes gene expression effects from genome-wide association study (GWAS) data using a statistical prediction model trained on GTEx reference transcriptome data. We then performed a colocalization analysis to evaluate the posterior probability of biological patterns: associations characterized by a single causal variant or multiple distinct causal variants. Finally, a functional enrichment analysis of gene sets was performed using the VarElect and CluePedia tools, which assess the causal relationships between genes and a disease and search for potential gene's functional pathways. The osteoporosis-associated genes were further confirmed based on the differentially expressed genes profiled from mRNA expression data of bone tissue.

RESULTS

Our analysis identified 204 candidate genes, including 154 genes that have been previously associated with osteoporosis, 50 genes that have not been previously discovered. A biological function analysis found that 20 of the candidate genes were directly associated with osteoporosis. Further analysis of multiple gene expression profiles showed that 15 genes were differentially expressed in patients with osteoporosis. Among these, SLC11A2, MAP2K5, NFATC4, and HSP90B1 were enriched in four pathways, namely, mineral absorption pathway, MAPK signaling pathway, Wnt signaling pathway, and PI3K-Akt signaling pathway, which indicates a causal relationship with the occurrence of osteoporosis.

CONCLUSIONS

We demonstrated that transcriptome fine-mapping identifies more osteoporosis-related genes and provides key insight into the development of novel targeted therapeutics for the treatment of osteoporosis.

Gaining insight into metabolic diseases from human genetic discoveries

Human large-scale genetic association studies have identified sequence variations at thousands of genetic risk loci that are more common in patients with diverse metabolic disease compared with healthy controls. While these genetic associations have been replicated in multiple large cohorts and sometimes can explain up to 50% of heritability, the molecular and cellular mechanisms affected by common genetic variation associated with metabolic disease remains mostly unknown. A variety of new genome-wide data types, in conjunction with novel biostatistical and computational analytical methodologies and foundational experimental technologies, are paving the way for a principled approach to systematic variant-to-function (V2F) studies for metabolic diseases, turning associated regions into causal variants, cell types and states of action, effector genes, and cellular and physiological mechanisms. Identification of new target genes and cellular programs for metabolic risk loci will improve mechanistic understanding of disease biology and identification of novel therapeutic strategies.

FIVEx: an interactive eQTL browser across public datasets

SUMMARY

Expression quantitative trait loci (eQTLs) characterize the associations between genetic variation and gene expression to provide insights into tissue-specific gene regulation. Interactive visualization of tissue-specific eQTLs or splice QTLs (sQTLs) can facilitate our understanding of functional variants relevant to disease-related traits. However, combining the multi-dimensional nature of eQTLs/sQTLs into a concise and informative visualization is challenging. Existing QTL visualization tools provide useful ways to summarize the unprecedented scale of transcriptomic data but are not necessarily tailored to answer questions about the functional interpretations of trait-associated variants or other variants of interest. We developed FIVEx, an interactive eQTL/sQTL browser with an intuitive interface tailored to the functional interpretation of associated variants. It features the ability to navigate seamlessly between different data views while providing relevant tissue- and locus-specific information to offer users a better understanding of population-scale multi-tissue transcriptomic profiles. Our implementation of the FIVEx browser on the EBI eQTL catalogue, encompassing 16 publicly available RNA-seq studies, provides important insights for understanding potential tissue-specific regulatory mechanisms underlying trait-associated signals.

AVAILABILITY AND IMPLEMENTATION

A FIVEx instance visualizing EBI eQTL catalogue data can be found at https://fivex.sph.umich.edu. Its source code is open source under an MIT license at https://github.com/statgen/fivex.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Integration analysis of GWAS and expression quantitative trait loci to identify candidate genes and pathways for clozapine-related neutropaenia

AIMS

Little is known about the genetic basis of clozapine-related neutropaenia. This study aims to explore candidate genes and pathways involved in clozapine-related neutropaenia.

METHODS

This study conducted a two-stage integrative analysis of the summary statistics from the genome-wide association study (GWAS, n = 552) of the lowest absolute neutrophil count (ANC) during clozapine treatment and the summary data of the expressed quantitative trait locus (eQTL). First, we use the probabilistic Mendelian randomization (PMR-Egger) to identify genes whose expression is causally related to ANC, and then use Bayesian co-localization analysis to investigate whether there are shared causal variants between them [posterior probability for hypotheses 4 (PP.H4) > 0.80]. Finally, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were conducted to explore the pathways that may be associated with ANC during clozapine treatment.

RESULTS

PMR-Egger analysis identified 146 genes that may be causally associated with ANC after Bonferroni correction (P-value < 3.25e-6). Bayesian co-localization analysis identified six further genes whose gene expression shared common variants with ANC, including NT5E (PP.H4 = 0.96), GLDC (PP.H4 = 0.82), NUDT17 (PP.H4 = 0.88), MSH4 (PP.H4 = 0.88), PTER (PP.H4 = 0.89) and SERPINB6 (PP.H4 = 0.83). Enrichment analysis identified 52 GO terms and seven pathways associated with ANC, such as NAD metabolic process, drug catabolic process and glyoxylate and dicarboxylate metabolism.

CONCLUSION

This study identified multiple candidate genes and pathways that may be involved in clozapine-related neutropaenia, providing novel clues for the mechanism of clozapine-related neutropaenia.

Wnt Signaling: From Mesenchymal Cell Fate to Lipogenesis and Other Mature Adipocyte Functions

Wnt signaling is an ancient and evolutionarily conserved pathway with fundamental roles in the development of adipose tissues. Roles of this pathway in mesenchymal stem cell fate determination and differentiation have been extensively studied. Indeed, canonical Wnt signaling is a significant endogenous inhibitor of adipogenesis and promoter of other cell fates, including osteogenesis, chondrogenesis, and myogenesis. However, emerging genetic evidence in both humans and mice suggests central roles for Wnt signaling in body fat distribution, obesity, and metabolic dysfunction. Herein, we highlight recent studies that have begun to unravel the contributions of various Wnt pathway members to critical adipocyte functions, including carbohydrate and lipid metabolism. We further explore compelling evidence of complex and coordinated interactions between adipocytes and other cell types within adipose tissues, including stromal, immune, and endothelial cells. Given the evolutionary conservation and ubiquitous cellular distribution of this pathway, uncovering the contributions of Wnt signaling to cell metabolism has exciting implications for therapeutic intervention in widespread pathologic states, including obesity, diabetes, and cancers.

Genetics of Obesity: What We Have Learned Over Decades of Research

There is a genetic component to human obesity that accounts for 40% to 50% of the variability in body weight status but that is lower among normal weight individuals (about 30%) and substantially higher in the subpopulation of individuals with obesity and severe obesity (about 60%-80%). The appreciation that heritability varies across classes of BMI represents an important advance. After controlling for BMI, ectopic fat and fat distribution traits are characterized by heritability levels ranging from 30% to 55%. Defects in at least 15 genes are the cause of monogenic obesity cases, resulting mostly from deficiencies in the leptin-melanocortin signaling pathway. Approximately two-thirds of the BMI heritability can be imputed to common DNA variants, whereas low-frequency and rare variants explain the remaining fraction. Diminishing allele effect size is observed as the number of obesity-associated variants expands, with most BMI-increasing or -decreasing alleles contributing only a few grams or less to body weight. Obesity-promoting alleles exert minimal effects in normal weight individuals but have larger effects in individuals with a proneness to obesity, suggesting a higher penetrance; however, it is not known whether these larger effect sizes precede obesity or are caused by an obese state. The obesity genetic risk is conditioned by thousands of DNA variants that make genetically based obesity prevention and treatment a major challenge.

Cell-type-specific expression quantitative trait loci associated with Alzheimer disease in blood and brain tissue

Because regulation of gene expression is heritable and context-dependent, we investigated AD-related gene expression patterns in cell types in blood and brain. Cis-expression quantitative trait locus (eQTL) mapping was performed genome-wide in blood from 5257 Framingham Heart Study (FHS) participants and in brain donated by 475 Religious Orders Study/Memory & Aging Project (ROSMAP) participants. The association of gene expression with genotypes for all cis SNPs within 1 Mb of genes was evaluated using linear regression models for unrelated subjects and linear-mixed models for related subjects. Cell-type-specific eQTL (ct-eQTL) models included an interaction term for the expression of "proxy" genes that discriminate particular cell type. Ct-eQTL analysis identified 11,649 and 2533 additional significant gene-SNP eQTL pairs in brain and blood, respectively, that were not detected in generic eQTL analysis. Of note, 386 unique target eGenes of significant eQTLs shared between blood and brain were enriched in apoptosis and Wnt signaling pathways. Five of these shared genes are established AD loci. The potential importance and relevance to AD of significant results in myeloid cell types is supported by the observation that a large portion of GWS ct-eQTLs map within 1 Mb of established AD loci and 58% (23/40) of the most significant eGenes in these eQTLs have previously been implicated in AD. This study identified cell-type-specific expression patterns for established and potentially novel AD genes, found additional evidence for the role of myeloid cells in AD risk, and discovered potential novel blood and brain AD biomarkers that highlight the importance of cell-type-specific analysis.

Integrating GWAS and eQTL to predict genes and pathways for non-syndromic cleft lip with or without palate

OBJECTIVE

To explore susceptibility genes and pathways for non-syndromic cleft lip with or without cleft palate (NSCL/P).

MATERIALS AND METHODS

Two genome-wide association studies (GWAS) datasets, including 858 NSCL/P cases and 1,248 controls, were integrated with expression quantitative trait loci (eQTL) dataset identified by Genotype-Tissue Expression (GTEx) project in whole-blood samples. The expression of the candidate genes in mouse orofacial development was inquired from FaceBase. Protein-protein interaction (PPI) network was visualized to identify protein functions. Go and KEGG pathway analyses were performed to explore the underlying risk pathways.

RESULTS

A total of 233 eQTL single-nucleotide polymorphisms (SNPs) in 432 candidate genes were identified to be associated with the risk of NSCL/P. One hundred and eighty-three susceptible genes were expressed in mouse orofacial development according to FaceBase. PPI network analysis highlighted that these genes involved in ubiquitin-mediated proteolysis (KCTD7, ASB1, UBOX5, ANAPC4) and DNA synthesis (XRCC3, RFC3, KAT5, RHNO1) were associated with the risk of NSCL/P. GO and KEGG pathway analyses revealed that the fatty acid metabolism pathway (ACADL, HSD17B12, ACSL5, PPT1, MCAT) played an important role in the development of NSCL/P.

CONCLUSIONS

Our results identified novel susceptibility genes and pathways associated with the development of NSCL/P.

Phenome-wide Mendelian randomization mapping the influence of the plasma proteome on complex diseases

The human proteome is a major source of therapeutic targets. Recent genetic association analyses of the plasma proteome enable systematic evaluation of the causal consequences of variation in plasma protein levels. Here we estimated the effects of 1,002 proteins on 225 phenotypes using two-sample Mendelian randomization (MR) and colocalization. Of 413 associations supported by evidence from MR, 130 (31.5%) were not supported by results of colocalization analyses, suggesting that genetic confounding due to linkage disequilibrium is widespread in naïve phenome-wide association studies of proteins. Combining MR and colocalization evidence in cis-only analyses, we identified 111 putatively causal effects between 65 proteins and 52 disease-related phenotypes ( https://www.epigraphdb.org/pqtl/ ). Evaluation of data from historic drug development programs showed that target-indication pairs with MR and colocalization support were more likely to be approved, evidencing the value of this approach in identifying and prioritizing potential therapeutic targets.

Functional genomics, genetic risk profiling and cell phenotypes in neurodegenerative disease

Human genetics provides unbiased insights into the causes of human disease, which can be used to create a foundation for effective ways to more accurately diagnose patients, stratify patients for more successful clinical trials, discover and develop new therapies, and ultimately help patients choose the safest and most promising therapeutic option based on their risk profile. But the process for translating basic observations from human genetics studies into pathogenic disease mechanisms and treatments is laborious and complex, and this challenge has particularly slowed the development of interventions for neurodegenerative disease. In this review, we discuss the many steps in the process, the important considerations at each stage, and some of the latest tools and technologies that are available to help investigators translate insights from human genetics into diagnostic and therapeutic strategies that will lead to the sort of advances in clinical care that make a difference for patients.

Genome-wide pQTL analysis of protein expression regulatory networks in the human liver

BACKGROUND

Previous expression quantitative trait loci (eQTL) studies have identified thousands of genetic variants to be associated with gene expression at the mRNA level in the human liver. However, protein expression often correlates poorly with mRNA levels. Thus, protein quantitative trait loci (pQTL) study is required to identify genetic variants that regulate protein expression in human livers.

RESULTS

We conducted a genome-wide pQTL study in 287 normal human liver samples and identified 900 local pQTL variants and 4026 distant pQTL variants. We further discovered 53 genome hotspots of pQTL variants. Transcriptional region mapping analysis showed that 1133 pQTL variants are in transcriptional regulatory regions. Genomic region enrichment analysis of the identified pQTL variants revealed 804 potential regulatory interactions among 595 predicted regulators (e.g., non-coding RNAs) and 394 proteins. Moreover, pQTL variants and trait-variant integration analysis implied several novel mechanisms underlying the relationships between protein expression and liver diseases, such as alcohol dependence. Notably, over 2000 of the identified pQTL variants have not been reported in previous eQTL studies, suggesting extensive involvement of genetic polymorphisms in post-transcriptional regulation of protein expression in human livers.

CONCLUSIONS

We have partially established protein expression regulation networks in human livers and generated a wealth of pQTL data that could serve as a valuable resource for the scientific community.