👤 Amadou Gaye

Type 2 diabetes (T2D) is a complex metabolic disorder driven by genetic and environmental factors. While genome-wide association studies (GWAS) have identified numerous T2D-associated variants, many remain functionally uncharacterized. Integration of GWAS with molecular phenotyping offers a path to revealing biological relevance. We investigated the influence of GWAS-variants, including sub-threshold T2D-associated variants (GWAS p-value ≤ 0.0001), on gene and protein expression to assign functional relevance. Genetic variants associated with T2D in the GWAS Catalog and present in our whole-genome sequencing (WGS) data were used to perform expression quantitative trait loci (eQTL) analysis in 242 whole-blood mRNA-sequenced samples. The same variants were used to perform protein quantitative trait loci (pQTL) analysis in a set of 362 plasma samples profiled on the Olink platform. For each analysis, the datasets were randomly split into discovery and validation subsets. Associations between variants and mRNA or protein levels were tested by multiple linear regression, and only QTLs that reached a false discovery rate adjusted p-value ≤ 0.05 in the discovery dataset and replicated in the validation dataset (p ≤ 0.05) with same direction of effect were carried forward. QTL-linked mRNAs and proteins were subsequently evaluated for their relationship with T2D status to connect them with T2D pathophysiology. We identified 1,291 eQTLs linked to 97 mRNAs and 1,273 pQTLs linked to 22 proteins. Among these, 10 mRNAs and 5 proteins were differentially expressed between non-diabetic and diabetic individuals. Notably, LPL, APOBR, APOM (lipid metabolism), NOTCH2, TREH (β-cell/endocrine regulation), and HLA-A, OAS3 (immune response) converged on three biological axes central to T2D pathophysiology. The directionality of molecular effects was consistent with known disease mechanisms, including insulin resistance (LPL, APOBR), β-cell stress (TREH, NOTCH2), and chronic inflammation (OAS3). Our findings indicate that variants falling below conventional GWAS significance thresholds can have demonstrable effects on gene expression and protein levels. This underscores the importance of prioritizing biological relevance alongside statistical significance, rather than relying solely on rigid p-value cutoffs. Show less

The analysis of rich catalogues of genetic variation from population-based sequencing provides an opportunity to screen for functional effects. Here we report a rare variant in APOC3 (rs138326449-A, minor allele frequency ~0.25% (UK)) associated with plasma triglyceride (TG) levels (-1.43 s.d. (s.e.=0.27 per minor allele (P-value=8.0 × 10(-8))) discovered in 3,202 individuals with low read-depth, whole-genome sequence. We replicate this in 12,831 participants from five additional samples of Northern and Southern European origin (-1.0 s.d. (s.e.=0.173), P-value=7.32 × 10(-9)). This is consistent with an effect between 0.5 and 1.5 mmol l(-1) dependent on population. We show that a single predicted splice donor variant is responsible for association signals and is independent of known common variants. Analyses suggest an independent relationship between rs138326449 and high-density lipoprotein (HDL) levels. This represents one of the first examples of a rare, large effect variant identified from whole-genome sequencing at a population scale. Show less