Methods for statistical fine-mapping and their applications to auto-immune diseases

From classical approaches to artificial intelligence, old and new tools for PDAC risk stratification and prediction

Pancreatic ductal adenocarcinoma (PDAC) is recognized as one of the most lethal malignancies, characterized by late-stage diagnosis and limited therapeutic options. Risk stratification has traditionally been performed using epidemiological studies and genetic analyses, through which key risk factors, including smoking, diabetes, chronic pancreatitis, and inherited predispositions, have been identified. However, the multifactorial nature of PDAC has often been insufficiently addressed by these methods, leading to limited precision in individualized risk assessments. Advances in artificial intelligence (AI) have been proposed as a transformative approach, allowing the integration of diverse datasets-spanning genetic, clinical, lifestyle, and imaging data into dynamic models capable of uncovering novel interactions and risk profiles. In this review, the evolution of PDAC risk stratification is explored, with classical epidemiological frameworks compared to AI-driven methodologies. Genetic insights, including genome-wide association studies and polygenic risk scores, are discussed, alongside AI models such as machine learning, radiomics, and deep learning. Strengths and limitations of these approaches are evaluated, with challenges in clinical translation, such as data scarcity, model interpretability, and external validation, addressed. Finally, future directions are proposed for combining classical and AI-driven methodologies to develop scalable, personalized predictive tools for PDAC, with the goal of improving early detection and patient outcomes.

Variant-to-function approaches for adipose tissue: Insights into cardiometabolic disorders

Genome-wide association studies (GWASs) have identified thousands of genetic loci associated with cardiometabolic disorders. However, the functional interpretation of these loci remains a daunting challenge. This is particularly true for adipose tissue, a critical organ in systemic metabolism and the pathogenesis of various cardiometabolic diseases. We discuss how variant-to-function (V2F) approaches are used to elucidate the mechanisms by which GWAS loci increase the risk of cardiometabolic disorders by directly influencing adipose tissue. We outline GWAS traits most likely to harbor adipose-related variants and summarize tools to pinpoint the putative causal variants, genes, and cell types for the associated loci. We explain how large-scale perturbation experiments, coupled with imaging and multi-omics, can be used to screen variants' effects on cellular phenotypes and how these phenotypes can be tied to physiological mechanisms. Lastly, we discuss the challenges and opportunities that lie ahead for V2F research and propose a roadmap for future studies.

Sequence-based genome-wide association study and fine-mapping in German Holstein reveal new quantitative trait loci for health traits

We conducted a large GWAS of 11 health traits belonging to 3 trait complexes: (1) metabolic diseases, (2) infectious and noninfectious feet and claw disorders, and (3) udder-related traits in 100,809 to 180,217 German Holstein cows to investigate the genetic architecture and underlying biological mechanisms behind these complex traits. The GWAS identified 12,306 genome-wide significant variants across 10 traits. The new association signals were inspected with a Bayesian fine-mapping approach, leading to the discovery of 159 novel variants with high potential for causality. Variants were in known and novel regions for the traits studied, leading to a list of 53 novel candidate genes. Our study represents the largest whole-genome sequence GWAS for health traits so far, hence ensuring the power to detect meaningful variants, especially when enhanced with fine-mapping.

Causal Variants in CHRNA1 and CHRNB1 Genes for Anti-acetylcholine Receptor Antibody Positive Myasthenia Gravis: Evidence from Bayesian Fine-Mapping and Genetic Association Study

Autoantibodies target the acetylcholine receptor (AChR) in 85% of myasthenia gravis (MG) patients. Genomic studies highlighted the association of genes encoding AChR subunits (CHRNA1 and CHRNB1) and MG in European populations. Additionally, Mendelian randomization revealed rs4151121 at the CHRNB1 locus as a potential causal variant. Here, we performed Bayesian fine-mapping of the CHRNA1 locus using GWAS summary statistics, a linkage disequilibrium matrix and functional annotations. The GWAS lead hit rs35274388 was identified as a causal variant overlapping with the promoter region (p < 0.01). Next, we performed a candidate gene study including 1038 participants from Serbia. Rs4151121 minor allele G was associated with late-onset MG (LOMG) (OR = 1.327, 95% CI = 1.084-1.625, p = 0.006, pperm = 0.007). Carriers of the rs4151121 GG and AG genotypes had an almost 1.5-fold increased risk of developing LOMG. A borderline association of the rs35274388 minor allele A with MG was observed (OR = 1.478, 95% CI = 1.009-2.166, p = 0.044, pperm = 0.060). Individuals with AA and GA genotypes also showed a nearly 1.5-fold higher risk of developing MG. In silico-identified causal variants at the CHRNA1 and CHRNB1 loci represent risk factors for MG in European populations, and to a greater extent for LOMG. Studies on non-European populations and functional research are needed to elucidate the role of AChR genes in the genetic architecture and development of MG.

The importance of functional genomics studies in precision rheumatology

Rheumatic diseases, those that affect the musculoskeletal system, cause significant morbidity. Among risk factors of these diseases is a significant genetic component. Recent advances in high-throughput omics techniques now allow a comprehensive profiling of patients at a genetic level through genome-wide association studies. Without functional interpretation of variants identified through these studies, clinical insight remains limited. Strategies include statistical fine-mapping that refine the list of variants in loci associated with disease, whilst colocalization techniques attempt to attribute function to variants that overlap a genetically active chromatin annotation. Functional validation using genome editing techniques can be used to further refine genetic signals and identify key pathways in cell types relevant to rheumatic disease biology. Insight gained from the combination of genetic studies and functional validation can be used to improve precision medicine in rheumatic diseases by allowing risk prediction and drug repositioning.

Association between Cytotoxic T-Lymphocyte-Associated Antigen 4 (CTLA-4) Locus and Early-Onset Anti-acetylcholine Receptor-Positive Myasthenia Gravis in Serbian Patients

Genome-wide association studies (GWAS) have provided strong evidence that early- and late-onset MG have different genetic backgrounds. Recent in silico analysis based on GWAS results revealed rs231735 and rs231770 variants within CTLA-4 locus as possible MG causative genetic factors. We aimed to explore the association of rs231735 and rs231770 with MG in a representative cohort of Serbian patients. We conducted an age-, sex-, and ethnicity-matched case-control study. Using TaqMan allele discrimination assays, the frequency of rs231735 and rs231770 genetic variants was examined in 447 AChR-MG patients and 447 matched controls. There was no significant association of rs231735 and rs231770 with the entire MG cohort (P > 0.05). Nevertheless, when stratifying patients into early-onset (n = 183) and late-onset MG (n = 264), we found early-onset patients had a significantly lower frequency of the rs231735 allele T compared to controls (OR = 0.734, 95% CI = 0.575-0.938, p10e6 permutation < 0.05), and rs231735 genotype TT and rs231770 genotype TT had a protective effect on early-onset MG (OR = 0.548, 95% CI = 0.339-0.888, and OR = 0.563, 95% CI = 0.314-1.011, p10e6 permutation < 0.05). Consequently, we found that individuals with the rs231735-rs231770 haplotype GC had a higher risk for developing early-onset MG (OR = 1.360, P = 0.027, p10e6 permutation < 0.05). Our results suggest that CTLA-4 rs231735 and rs231770 may be risk factors only for patients with early-onset MG in Serbian population.

Genotype imputation methods for whole and complex genomic regions utilizing deep learning technology

The imputation of unmeasured genotypes is essential in human genetic research, particularly in enhancing the power of genome-wide association studies and conducting subsequent fine-mapping. Recently, several deep learning-based genotype imputation methods for genome-wide variants with the capability of learning complex linkage disequilibrium patterns have been developed. Additionally, deep learning-based imputation has been applied to a distinct genomic region known as the major histocompatibility complex, referred to as HLA imputation. Despite their various advantages, the current deep learning-based genotype imputation methods do have certain limitations and have not yet become standard. These limitations include the modest accuracy improvement over statistical and conventional machine learning-based methods. However, their benefits include other aspects, such as their "reference-free" nature, which ensures complete privacy protection, and their higher computational efficiency. Furthermore, the continuing evolution of deep learning technologies is expected to contribute to further improvements in prediction accuracy and usability in the future.

Genotype imputation in human genomic studies

Imputation is a method that supplies missing information about genetic variants that could not be directly genotyped with DNA microarrays or low-coverage sequencing. Imputation plays a critical role in genome-wide association studies (GWAS). It leads to a significant increase in the number of studied variants, which improves the resolution of the method and enhances the comparability of data obtained in different cohorts and/or by using different technologies, which is important for conducting meta-analyses. When performing imputation, genotype information from the study sample, in which only part of the genetic variants are known, is complemented using the standard (reference) sample, which has more complete genotype data (most often the results of whole-genome sequencing). Imputation has become an integral part of human genomic research due to the benefits it provides and the increasing availability of imputation tools and reference sample data. This review focuses on imputation in human genomic research. The first section of the review provides a description of technologies for obtaining information about human genotypes and characteristics of these types of data. The second section describes the imputation methodology, lists the stages of its implementation and the corresponding programs, provides a description of the most popular reference panels and methods for assessing the quality of imputation. The review concludes with examples of the use of imputation in genomic studies of samples from Russia. This review shows the importance of imputation, provides information on how to carry it out, and systematizes the results of its application using Russian samples.

The goldmine of GWAS summary statistics: a systematic review of methods and tools

Genome-wide association studies (GWAS) have revolutionized our understanding of the genetic architecture of complex traits and diseases. GWAS summary statistics have become essential tools for various genetic analyses, including meta-analysis, fine-mapping, and risk prediction. However, the increasing number of GWAS summary statistics and the diversity of software tools available for their analysis can make it challenging for researchers to select the most appropriate tools for their specific needs. This systematic review aims to provide a comprehensive overview of the currently available software tools and databases for GWAS summary statistics analysis. We conducted a comprehensive literature search to identify relevant software tools and databases. We categorized the tools and databases by their functionality, including data management, quality control, single-trait analysis, and multiple-trait analysis. We also compared the tools and databases based on their features, limitations, and user-friendliness. Our review identified a total of 305 functioning software tools and databases dedicated to GWAS summary statistics, each with unique strengths and limitations. We provide descriptions of the key features of each tool and database, including their input/output formats, data types, and computational requirements. We also discuss the overall usability and applicability of each tool for different research scenarios. This comprehensive review will serve as a valuable resource for researchers who are interested in using GWAS summary statistics to investigate the genetic basis of complex traits and diseases. By providing a detailed overview of the available tools and databases, we aim to facilitate informed tool selection and maximize the effectiveness of GWAS summary statistics analysis.

Identifying drug targets for schizophrenia through gene prioritization

Background

Schizophrenia genome-wide association studies (GWASes) have identified >250 significant loci and prioritized >100 disease-related genes. However, gene prioritization efforts have mostly been restricted to locus-based methods that ignore information from the rest of the genome.

Methods

To more accurately characterize genes involved in schizophrenia etiology, we applied a combination of highly-predictive tools to a published GWAS of 67,390 schizophrenia cases and 94,015 controls. We combined both locus-based methods (fine-mapped coding variants, distance to GWAS signals) and genome-wide methods (PoPS, MAGMA, ultra-rare coding variant burden tests). To validate our findings, we compared them with previous prioritization efforts, known neurodevelopmental genes, and results from the PsyOPS tool.

Results

We prioritized 62 schizophrenia genes, 41 of which were also highlighted by our validation methods. In addition to DRD2, the principal target of antipsychotics, we prioritized 9 genes that are targeted by approved or investigational drugs. These included drugs targeting glutamatergic receptors (GRIN2A and GRM3), calcium channels (CACNA1C and CACNB2), and GABAB receptor (GABBR2). These also included genes in loci that are shared with an addiction GWAS (e.g. PDE4B and VRK2).

Conclusions

We curated a high-quality list of 62 genes that likely play a role in the development of schizophrenia. Developing or repurposing drugs that target these genes may lead to a new generation of schizophrenia therapies. Rodent models of addiction more closely resemble the human disorder than rodent models of schizophrenia. As such, genes prioritized for both disorders could be explored in rodent addiction models, potentially facilitating drug development.

Relationship between HLA genetic variations, COVID-19 vaccine antibody response, and risk of breakthrough outcomes

The rapid global distribution of COVID-19 vaccines, with over a billion doses administered, has been unprecedented. However, in comparison to most identified clinical determinants, the implications of individual genetic factors on antibody responses post-COVID-19 vaccination for breakthrough outcomes remain elusive. Here, we conducted a population-based study including 357,806 vaccinated participants with high-resolution HLA genotyping data, and a subset of 175,000 with antibody serology test results. We confirmed prior findings that single nucleotide polymorphisms associated with antibody response are predominantly located in the Major Histocompatibility Complex region, with the expansive HLA-DQB1*06 gene alleles linked to improved antibody responses. However, our results did not support the claim that this mutation alone can significantly reduce COVID-19 risk in the general population. In addition, we discovered and validated six HLA alleles (A*03:01, C*16:01, DQA1*01:02, DQA1*01:01, DRB3*01:01, and DPB1*10:01) that independently influence antibody responses and demonstrated a combined effect across HLA genes on the risk of breakthrough COVID-19 outcomes. Lastly, we estimated that COVID-19 vaccine-induced antibody positivity provides approximately 20% protection against infection and 50% protection against severity. These findings have immediate implications for functional studies on HLA molecules and can inform future personalised vaccination strategies.

Estimating gene-level false discovery probability improves eQTL statistical fine-mapping precision

Statistical fine-mapping prioritizes putative causal variants from a large number of candidate variants, and is widely used in expression quantitative loci (eQTLs) studies. In eQTL fine-mapping, the existence of causal variants for gene expression is not guaranteed, since the genetic heritability of gene expression explained by nearby (cis-) variants is limited. Here we introduce a refined fine-mapping algorithm, named Knockoff-Finemap combination (KFc). KFc estimates the probability that the causal variant(s) exist in the cis-window of a gene through construction of knockoff genotypes (i.e. a set of synthetic genotypes that resembles the original genotypes), and uses it to adjust the posterior inclusion probabilities (PIPs). Utilizing simulated gene expression data, we show that KFc results in calibrated PIP distribution with improved precision. When applied to gene expression data of 465 genotyped samples from the Japan COVID-19 Task Force (JCTF), KFc resulted in significant enrichment of a functional score as well as reporter assay hits in the top PIP bins. When combined with functional priors derived from an external fine-mapping study (GTEx), KFc resulted in a significantly higher proportion of hematopoietic trait putative causal variants in the top PIP bins. Our work presents improvements in the precision of a major fine-mapping algorithm.

Coexistence of Multiple Functional Variants and Genes Underlies Genetic Risk Locus 11p11.2 of Alzheimer's Disease

BACKGROUND

Genome-wide association studies have identified dozens of genetic risk loci for Alzheimer's disease (AD), yet the underlying causal variants and biological mechanisms remain elusive, especially for loci with complex linkage disequilibrium and regulation.

METHODS

To fully untangle the causal signal at a single locus, we performed a functional genomic study of 11p11.2 (the CELF1/SPI1 locus). Genome-wide association study signals at 11p11.2 were integrated with datasets of histone modification, open chromatin, and transcription factor binding to distill potentially functional variants (fVars). Their allelic regulatory activities were confirmed by allele imbalance, reporter assays, and base editing. Expressional quantitative trait loci and chromatin interaction data were incorporated to assign target genes to fVars. The relevance of these genes to AD was assessed by convergent functional genomics using bulk brain and single-cell transcriptomic, epigenomic, and proteomic datasets of patients with AD and control individuals, followed by cellular assays.

RESULTS

We found that 24 potential fVars, rather than a single variant, were responsible for the risk of 11p11.2. These fVars modulated transcription factor binding and regulated multiple genes by long-range chromatin interactions. Besides SPI1, convergent evidence indicated that 6 target genes (MTCH2, ACP2, NDUFS3, PSMC3, C1QTNF4, and MADD) of fVars were likely to be involved in AD development. Disruption of each gene led to cellular amyloid-β and phosphorylated tau changes, supporting the existence of multiple likely causal genes at 11p11.2.

CONCLUSIONS

Multiple variants and genes at 11p11.2 may contribute to AD risk. This finding provides new insights into the mechanistic and therapeutic challenges of AD.

Identification of genetic variants associated with diabetic kidney disease in multiple Korean cohorts via a genome-wide association study mega-analysis

BACKGROUND

The pathogenesis of diabetic kidney disease (DKD) is complex, involving metabolic and hemodynamic factors. Although DKD has been established as a heritable disorder and several genetic studies have been conducted, the identification of unique genetic variants for DKD is limited by its multiplex classification based on the phenotypes of diabetes mellitus (DM) and chronic kidney disease (CKD). Thus, we aimed to identify the genetic variants related to DKD that differentiate it from type 2 DM and CKD.

METHODS

We conducted a large-scale genome-wide association study mega-analysis, combining Korean multi-cohorts using multinomial logistic regression. A total of 33,879 patients were classified into four groups-normal, DM without CKD, CKD without DM, and DKD-and were further analyzed to identify novel single-nucleotide polymorphisms (SNPs) associated with DKD. Additionally, fine-mapping analysis was conducted to investigate whether the variants of interest contribute to a trait. Conditional analyses adjusting for the effect of type 1 DM (T1D)-associated HLA variants were also performed to remove confounding factors of genetic association with T1D. Moreover, analysis of expression quantitative trait loci (eQTL) was performed using the Genotype-Tissue Expression project. Differentially expressed genes (DEGs) were analyzed using the Gene Expression Omnibus database (GSE30529). The significant eQTL DEGs were used to explore the predicted interaction networks using search tools for the retrieval of interacting genes and proteins.

RESULTS

We identified three novel SNPs [rs3128852 (P = 8.21×10^-25), rs117744700 (P = 8.28×10^-10), and rs28366355 (P = 2.04×10^-8)] associated with DKD. Moreover, the fine-mapping study validated the causal relationship between rs3128852 and DKD. rs3128852 is an eQTL for TRIM27 in whole blood tissues and HLA-A in adipose-subcutaneous tissues. rs28366355 is an eQTL for HLA-group genes present in most tissues.

CONCLUSIONS

We successfully identified SNPs (rs3128852, rs117744700, and rs28366355) associated with DKD and verified the causal association between rs3128852 and DKD. According to the in silico analysis, TRIM27 and HLA-A can define DKD pathophysiology and are associated with immune response and autophagy. However, further research is necessary to understand the mechanism of immunity and autophagy in the pathophysiology of DKD and to prevent and treat DKD.

SysBiolPGWAS: simplifying post-GWAS analysis through the use of computational technologies and integration of diverse omics datasets

MOTIVATION

Post-genome-wide association studies (pGWAS) analysis is designed to decipher the functional consequences of significant single-nucleotide polymorphisms (SNPs) in the era of GWAS. This can be translated into research insights and clinical benefits such as the effectiveness of strategies for disease screening, treatment and prevention. However, the setup of pGWAS (pGWAS) tools can be quite complicated, and it mostly requires big data. The challenge however is, scientists are required to have sufficient experience with several of these technically complex and complicated tools in order to complete the pGWAS analysis.

RESULTS

We present SysBiolPGWAS, a pGWAS web application that provides a comprehensive functionality for biologists and non-bioinformaticians to conduct several pGWAS analyses to overcome the above challenges. It provides unique functionalities for analysis involving multi-omics datasets and visualization using various bioinformatics tools. SysBiolPGWAS provides access to individual pGWAS tools and a novel custom pGWAS pipeline that integrates several individual pGWAS tools and data. The SysBiolPGWAS app was developed to be a one-stop shop for pGWAS analysis. It targets researchers in the area of the human genome and performs its analysis mainly in the autosomal chromosomes.

AVAILABILITY AND IMPLEMENTATION

SysBiolPGWAS web app was developed using JavaScript/TypeScript web frameworks and is available at: https://spgwas.waslitbre.org/. All codes are available in this GitHub repository https://github.com/covenant-university-bioinformatics.

3D genome organization links non-coding disease-associated variants to genes

Genome sequencing has revealed over 300 million genetic variations in human populations. Over 90% of variants are single nucleotide polymorphisms (SNPs), the remainder include short deletions or insertions, and small numbers of structural variants. Hundreds of thousands of these variants have been associated with specific phenotypic traits and diseases through genome wide association studies which link significant differences in variant frequencies with specific phenotypes among large groups of individuals. Only 5% of disease-associated SNPs are located in gene coding sequences, with the potential to disrupt gene expression or alter of the function of encoded proteins. The remaining 95% of disease-associated SNPs are located in non-coding DNA sequences which make up 98% of the genome. The role of non-coding, disease-associated SNPs, many of which are located at considerable distances from any gene, was at first a mystery until the discovery that gene promoters regularly interact with distal regulatory elements to control gene expression. Disease-associated SNPs are enriched at the millions of gene regulatory elements that are dispersed throughout the non-coding sequences of the genome, suggesting they function as gene regulation variants. Assigning specific regulatory elements to the genes they control is not straightforward since they can be millions of base pairs apart. In this review we describe how understanding 3D genome organization can identify specific interactions between gene promoters and distal regulatory elements and how 3D genomics can link disease-associated SNPs to their target genes. Understanding which gene or genes contribute to a specific disease is the first step in designing rational therapeutic interventions.