A genetics-led approach defines the drug target landscape of 30 immune-related traits

Revealing novel genomic insights and therapeutic targets for juvenile idiopathic arthritis through omics

BACKGROUND

The genetic architecture of JIA remains only partially comprehended. There is a clear imperative for continued endeavours to uncover insights into the underlying causes of JIA.

METHODS

This study encompassed a comprehensive spectrum of endeavours, including conducting a JIA genome-wide association study (GWAS) meta-analysis that incorporated data from 4550 JIA cases and 18 446 controls. We employed in silico and genome-editing approaches to prioritizing target genes. To investigate pleiotropic effects, we conducted phenome-wide association studies. Cell-type enrichment analyses were performed by integrating bulk and single-cell sequencing data. Finally, we delved into potential druggable targets for JIA.

RESULTS

Fourteen genome-wide significant non-HLA loci were identified, including four novel loci, each exhibiting pleiotropic associations with other autoimmune diseases or musculoskeletal traits. We uncovered strong genetic correlation between JIA and BMD traits at 52 genomic regions, including three GWAS loci for JIA. Candidate genes with immune functions were captured by in silico analyses at each novel locus, with additional findings identified through our experimental approach. Cell-type enrichment analysis revealed 21 specific immune cell types crucial for the affected organs in JIA, indicating their potential contribution to the disease. Finally, 24 known or candidate druggable target genes were prioritized.

CONCLUSIONS

Our identification of four novel JIA-associated genes, CD247, RHOH, COLEC10 and IRF8, broadens the novel potential drug repositioning opportunities. We established a new genetic link between COLEC10, TNFRSF11B and JIA/BMD. Additionally, the identification of RHOH underscores its role in positive thymocyte selection, thereby illuminating a critical facet of JIA's underlying biological mechanisms.

XGRm: A Web Server for Interpreting Mouse Summary-level Genomic Data

We introduce XGR-model (or XGRm), a web server made accessible at http://www.xgrm.pro, with the aim of meeting the increasing demand for effectively interpreting summary-level genomic data in model organisms. Currently, it hosts two enrichment analysers and two subnetwork analysers to support enrichment and subnetwork analyses for user-input mouse genomic data, whether gene-centric or genomic region-centric. The enrichment analysers identify ontology term enrichments for input genes (GElyser) or for genes linked from input genomic regions (RElyser). The subnetwork analysers rely on our previously established network algorithm to identify gene subnetworks from input gene-centric summary data (GSlyser) or from input region-centric summary data (RSlyser), leveraging network information about either functional interactions or pathway-derived interactions. Collectively, XGRm offers an all-in-one solution for gaining systems biology insights into summary-level genomic data in mice, underpinned by our commitment to regular updates as well as natural extensions to other model organisms.

Identifying genetic determinants of sarcopenia-related traits: a Mendelian randomization study of druggable genes

BACKGROUND

Sarcopenia, characterized by progressive muscle mass and function loss, particularly affects the elderly, and leads to severe consequences such as falls and mortality. Despite its prevalence, targeted pharmacotherapies for sarcopenia are lacking. Utilizing large-sample genome-wide association studies (GWAS) data is crucial for cost-effective drug discovery.

METHODS

Herein, we conducted four studies to understand the putative causal effects of genetic components on muscle mass and function. Study 1 employed a two-sample Mendelian randomization (MR) on 15,944 potential druggable genes, investigating their potential causality with muscle quantity and quality in a European population (N up to 461,089). Study 2 validated MR results through sensitivity analyses and colocalization analyses. Study 3 extended validation across other European cohorts, and study 4 conducted quantitative in vivo verification.

RESULTS

MR analysis revealed significant causality between four genes (BLOC-1 related complex subunit 7, BORCS7; peptidase m20 domain containing 1, PM20D1; nuclear casein kinase and cyclin dependent kinase substrate 1, NUCKS1 and ubiquinol-cytochrome c reductase complex assembly factor 1, UQCC1) and muscle mass and function (p-values range 5.98 × 10-6 to 9.26 × 10-55). To be specific, BORCS7 and UQCC1 negatively regulated muscle quantity and quality, whereas enhancing PM20D1 and NUCKS1 expression showed promise in promoting muscle mass and function. Causal relationships remained robust across sensitivity analyses, with UQCC1 exhibiting notable colocalization effects (PP·H4 93.4 % to 95.8 %). Further validation and in vivo replication verified the potential causality between these genes and muscle mass as well as function.

CONCLUSIONS

Our druggable genome-wide MR analysis identifies BORCS7, PM20D1, NUCKS1, and UQCC1 as causally associated with muscle mass and function. These findings offer insights into the genetic basis of sarcopenia, paving the way for these genes to become promising drug targets in mitigating this debilitating condition.

From genetic associations to genes: methods, applications, and challenges

Genome-wide association studies (GWASs) have identified numerous genetic loci associated with human traits and diseases. However, pinpointing the causal genes remains a challenge, which impedes the translation of GWAS findings into biological insights and medical applications. In this review, we provide an in-depth overview of the methods and technologies used for prioritizing genes from GWAS loci, including gene-based association tests, integrative analysis of GWAS and molecular quantitative trait loci (xQTL) data, linking GWAS variants to target genes through enhancer-gene connection maps, and network-based prioritization. We also outline strategies for generating context-dependent xQTL data and their applications in gene prioritization. We further highlight the potential of gene prioritization in drug repurposing. Lastly, we discuss future challenges and opportunities in this field.

Instrumental variable and colocalization analyses identify endotrophin and HTRA1 as potential therapeutic targets for coronary artery disease

Coronary artery disease (CAD) remains a leading cause of disease burden globally, and there is a persistent need for new therapeutic targets. Instrumental variable (IV) and genetic colocalization analyses can help identify novel therapeutic targets for human disease by nominating causal genes in genome-wide association study (GWAS) loci. We conducted cis-IV analyses for 20,125 genes and 1,746 plasma proteins with CAD using molecular trait quantitative trait loci variant (QTLs) data from three different studies. 19 proteins and 119 genes were significantly associated with CAD risk by IV analyses and demonstrated evidence of genetic colocalization. Notably, our analyses validated well-established targets such as PCSK9 and ANGPTL4 while also identifying HTRA1 and endotrophin (a cleavage product of COL6A3) as proteins whose levels are causally associated with CAD risk. Further experimental studies are needed to confirm the causal role of the genes and proteins identified through our multiomic cis-IV analyses on human disease.

Genetic and Epigenetic Landscape for Drug Development in Polycystic Ovary Syndrome

The treatment of polycystic ovary syndrome (PCOS) faces challenges as all known treatments are merely symptomatic. The US Food and Drug Administration has not approved any drug specifically for treating PCOS. As the significance of genetics and epigenetics rises in drug development, their pivotal insights have greatly enhanced the efficacy and success of drug target discovery and validation, offering promise for guiding the advancement of PCOS treatments. In this context, we outline the genetic and epigenetic advancement in PCOS, which provide novel insights into the pathogenesis of this complex disease. We also delve into the prospective method for harnessing genetic and epigenetic strategies to identify potential drug targets and ensure target safety. Additionally, we shed light on the preliminary evidence and distinctive challenges associated with gene and epigenetic therapies in the context of PCOS.

Genetics of rheumatoid arthritis

In the past four decades, a plethora of genetic association studies have been carried out in cohorts of patients with rheumatoid arthritis. These studies have highlighted key aspects of disease pathogenesis and suggested causal mechanisms. In this review, we discuss major advances in our understanding of the genetic architecture of rheumatoid arthritis susceptibility, severity and treatment response and explain how genetics supports current models of disease pathogenesis and outcome. We outline future research directions, like Mendelian randomisation, and present a number of potential avenues for clinical translation, including risk and outcome prediction, patient stratification into treatment response groups and pharmacological applications.

Drug repurposing opportunities for chronic kidney disease

The development of targeted drugs for the early prevention and management of chronic kidney disease (CKD) is of great importance. However, the success rates and cost-effectiveness of traditional drug development approaches are extremely low. Utilizing large sample genome-wide association study data for drug repurposing has shown promise in many diseases but has not yet been explored in CKD. Herein, we investigated actionable druggable targets to improve renal function using large-scale Mendelian randomization and colocalization analyses. We combined two population-scale independent genetic datasets and validated findings with cell-type-dependent eQTL data of kidney tubular and glomerular samples. We ultimately prioritized two drug targets, opioid receptor-like 1 and F12, with potential genetic support for restoring renal function and subsequent treatment of CKD. Our findings explore the potential pathological mechanisms of CKD, bridge the gap between the molecular mechanisms of pathogenesis and clinical intervention, and provide new strategies in future clinical trials of CKD.

A Novel Macrophage Subpopulation Conveys Increased Genetic Risk of Coronary Artery Disease

BACKGROUND

Coronary artery disease (CAD), the leading cause of death worldwide, is influenced by both environmental and genetic factors. Although over 250 genetic risk loci have been identified through genome-wide association studies, the specific causal variants and their regulatory mechanisms are still largely unknown, particularly in disease-relevant cell types such as macrophages.

METHODS

We utilized single-cell RNA-seq and single-cell multiomics approaches in primary human monocyte-derived macrophages to explore the transcriptional regulatory network involved in a critical pathogenic event of coronary atherosclerosis-the formation of lipid-laden foam cells. The relative genetic contribution to CAD was assessed by partitioning disease heritability across different macrophage subpopulations. Meta-analysis of single-cell RNA-seq data sets from 38 human atherosclerotic samples was conducted to provide high-resolution cross-referencing to macrophage subpopulations in vivo.

RESULTS

We identified 18 782 cis-regulatory elements by jointly profiling the gene expression and chromatin accessibility of >5000 macrophages. Integration with CAD genome-wide association study data prioritized 121 CAD-related genetic variants and 56 candidate causal genes. We showed that CAD heritability was not uniformly distributed and was particularly enriched in the gene programs of a novel CD52-hi lipid-handling macrophage subpopulation. These CD52-hi macrophages displayed significantly less lipoprotein accumulation and were also found in human atherosclerotic plaques. We investigated the cis-regulatory effect of a risk variant rs10488763 on FDX1, implicating the recruitment of AP-1 and C/EBP-β in the causal mechanisms at this locus.

CONCLUSIONS

Our results provide genetic evidence of the divergent roles of macrophage subsets in atherogenesis and highlight lipid-handling macrophages as a key subpopulation through which genetic variants operate to influence disease. These findings provide an unbiased framework for functional fine-mapping of genome-wide association study results using single-cell multiomics and offer new insights into the genotype-environment interactions underlying atherosclerotic disease.

Regulatory landscape enrichment analysis (RLEA): a computational toolkit for non-coding variant enrichment and cell type prioritization

BACKGROUND

As genomic studies continue to implicate non-coding sequences in disease, testing the roles of these variants requires insights into the cell type(s) in which they are likely to be mediating their effects. Prior methods for associating non-coding variants with cell types have involved approaches using linkage disequilibrium or ontological associations, incurring significant processing requirements. GaiaAssociation is a freely available, open-source software that enables thousands of genomic loci implicated in a phenotype to be tested for enrichment at regulatory loci of multiple cell types in minutes, permitting insights into the cell type(s) mediating the studied phenotype.

RESULTS

In this work, we present Regulatory Landscape Enrichment Analysis (RLEA) by GaiaAssociation and demonstrate its capability to test the enrichment of 12,133 variants across the cis-regulatory regions of 44 cell types. This analysis was completed in 134.0 ± 2.3 s, highlighting the efficient processing provided by GaiaAssociation. The intuitive interface requires only four inputs, offers a collection of customizable functions, and visualizes variant enrichment in cell-type regulatory regions through a heatmap matrix. GaiaAssociation is available on PyPi for download as a command line tool or Python package and the source code can also be installed from GitHub at https://github.com/GreallyLab/gaiaAssociation .

CONCLUSIONS

GaiaAssociation is a novel package that provides an intuitive and efficient resource to understand the enrichment of non-coding variants across the cis-regulatory regions of different cells, empowering studies seeking to identify disease-mediating cell types.

A cross-disease, pleiotropy-driven approach for therapeutic target prioritization and evaluation

Cross-disease genome-wide association studies (GWASs) unveil pleiotropic loci, mostly situated within the non-coding genome, each of which exerts pleiotropic effects across multiple diseases. However, the challenge "W-H-W" (namely, whether, how, and in which specific diseases pleiotropy can inform clinical therapeutics) calls for effective and integrative approaches and tools. We here introduce a pleiotropy-driven approach specifically designed for therapeutic target prioritization and evaluation from cross-disease GWAS summary data, with its validity demonstrated through applications to two systems of disorders (neuropsychiatric and inflammatory). We illustrate its improved performance in recovering clinical proof-of-concept therapeutic targets. Importantly, it identifies specific diseases where pleiotropy informs clinical therapeutics. Furthermore, we illustrate its versatility in accomplishing advanced tasks, including pathway crosstalk identification and downstream crosstalk-based analyses. To conclude, our integrated solution helps bridge the gap between pleiotropy studies and therapeutics discovery.

Proteome-wide Mendelian randomization identifies therapeutic targets for ankylosing spondylitis

Background

Ankylosing Spondylitis (AS) is a chronic inflammatory disorder which can lead to considerable pain and disability. Mendelian randomization (MR) has been extensively applied for repurposing licensed drugs and uncovering new therapeutic targets. Our objective is to pinpoint innovative therapeutic protein targets for AS and assess the potential adverse effects of druggable proteins.

Methods

We conducted a comprehensive proteome-wide MR study to assess the causal relationships between plasma proteins and the risk of AS. The plasma proteins were sourced from the UK Biobank Pharma Proteomics Project (UKB-PPP) database, encompassing GWAS data for 2,940 plasma proteins. Additionally, GWAS data for AS were extracted from the R9 version of the Finnish database, including 2,860 patients and 270,964 controls. The colocalization analysis was executed to identify shared causal variants between plasma proteins and AS. Finally, we examined the potential adverse effects of druggable proteins for AS therapy by conducting a phenome-wide association study (PheWAS) utilizing the extensive Finnish database in version R9, encompassing 2,272 phenotypes categorized into 46 groups.

Results

The findings revealed a positive genetic association between the predicted plasma levels of six proteins and an elevated risk of AS, while two proteins exhibited an inverse association with AS risk (P fdr < 0.05). Among these eight plasma proteins, colocalization analysis identified AIF1, TNF, FKBPL, AGER, ALDH5A1, and ACOT13 as shared variation with AS(PPH3+PPH4>0.8), suggesting that they represent potential direct targets for AS intervention. Further phenotype-wide association studies have shown some potential side effects of these six targets (P fdr < 0.05).

Conclusion

Our investigation examined the causal connections between six plasma proteins and AS, providing a comprehensive understanding of potential therapeutic targets.

Integration of human organoids single-cell transcriptomic profiles and human genetics repurposes critical cell type-specific drug targets for severe COVID-19

Human organoids recapitulate the cell type diversity and function of their primary organs holding tremendous potentials for basic and translational research. Advances in single-cell RNA sequencing (scRNA-seq) technology and genome-wide association study (GWAS) have accelerated the biological and therapeutic interpretation of trait-relevant cell types or states. Here, we constructed a computational framework to integrate atlas-level organoid scRNA-seq data, GWAS summary statistics, expression quantitative trait loci, and gene-drug interaction data for distinguishing critical cell populations and drug targets relevant to coronavirus disease 2019 (COVID-19) severity. We found that 39 cell types across eight kinds of organoids were significantly associated with COVID-19 outcomes. Notably, subset of lung mesenchymal stem cells increased proximity with fibroblasts predisposed to repair COVID-19-damaged lung tissue. Brain endothelial cell subset exhibited significant associations with severe COVID-19, and this cell subset showed a notable increase in cell-to-cell interactions with other brain cell types, including microglia. We repurposed 33 druggable genes, including IFNAR2, TYK2, and VIPR2, and their interacting drugs for COVID-19 in a cell-type-specific manner. Overall, our results showcase that host genetic determinants have cellular-specific contribution to COVID-19 severity, and identification of cell type-specific drug targets may facilitate to develop effective therapeutics for treating severe COVID-19 and its complications.

Priority index for critical Covid-19 identifies clinically actionable targets and drugs

While genome-wide studies have identified genomic loci in hosts associated with life-threatening Covid-19 (critical Covid-19), the challenge of resolving these loci hinders further identification of clinically actionable targets and drugs. Building upon our previous success, we here present a priority index solution designed to address this challenge, generating the target and drug resource that consists of two indexes: the target index and the drug index. The primary purpose of the target index is to identify clinically actionable targets by prioritising genes associated with Covid-19. We illustrate the validity of the target index by demonstrating its ability to identify pre-existing Covid-19 phase-III drug targets, with the majority of these targets being found at the leading prioritisation (leading targets). These leading targets have their evolutionary origins in Amniota ('four-leg vertebrates') and are predominantly involved in cytokine-cytokine receptor interactions and JAK-STAT signaling. The drug index highlights opportunities for repurposing clinically approved JAK-STAT inhibitors, either individually or in combination. This proposed strategic focus on the JAK-STAT pathway is supported by the active pursuit of therapeutic agents targeting this pathway in ongoing phase-II/III clinical trials for Covid-19.

Identifying a panel of nine genes as novel specific model in endometriosis noninvasive diagnosis

OBJECTIVE

To study biomarkers to develop a novel diagnosis model for endometriosis and validate it using clinical samples.

DESIGN

We used publicly available data sets and weighted gene coexpression network analysis to identify differentially expressed genes. Ten machine learning algorithms were used to develop an integrative model for predicting endometriosis. The accuracy and robustness of the model were validated using data sets and clinical samples.

SETTING

Department of Obstetrics and Gynecology, Tangdu Hospital, Air Force Medical University, Xi'an, Shaanxi, China.

PATIENT(S)

The study included clinical patients between the ages of 20 and 40 years who required laparoscopic surgery and who had not undergone hormone therapy within the previous 3 months. All the healthy individuals had given birth to a child at least once in their lives. Patients with inflammatory conditions, malignant diseases, immune diseases, myoma, or adenomyosis were excluded. Paraffin blocks of the samples were collected (case, n = 5; control, n = 5). Blood samples of 58 individuals were collected (case, n = 28; control, n = 30).

INTERVENTION(S)

None.

MAIN OUTCOME MEASURE(S)

The areas under the receiver operator characteristic curve of our diagnostic model were measured for data sets and clinical samples. Multiplex immunohistochemical staining and real-time quantitative polymerase chain reaction assays were used for the validation of the model from tissue slides and peripheral blood samples.

RESULT(S)

A nine-gene panel endometriosis messenger RNA score (EMScore), was constructed to distinguish the patients with endometriosis from healthy individuals using algorithms. The EMScore accurately predicted endometriosis, and the areas under the receiver operator characteristic curve of our diagnostic model were 0.920, and 0.942 for tissue and blood samples, respectively. Moreover, the EMScore outperformed other acknowledged signatures for predicting endometriosis across seven clinical cohorts. Overall, the EMScore constitutes a sensitive and specific noninvasive diagnostic method for endometriosis.

CONCLUSION(S)

We developed the EMScore, a novel model that can aid in the diagnosis of endometriosis using peripheral blood samples. This study will contribute to the development of improved clinical noninvasive and sensitive diagnostic tools for endometriosis. These nine genes might be potential target molecules for treating endometriosis.

Cross-population applications of genomics to understand the risk of multifactorial traits involving inflammation and immunity

The interplay between genetic and environmental factors plays a significant role in interindividual variation in immune and inflammatory responses. The availability of high-throughput low-cost genotyping and next-generation sequencing has revolutionized our ability to identify human genetic variation and understand how this varies within and between populations, and the relationship with disease. In this review, we explore the potential of genomics for patient benefit, specifically in the diagnosis, prognosis and treatment of inflammatory and immune-related diseases. We summarize the knowledge arising from genetic and functional genomic approaches, and the opportunity for personalized medicine. The review covers applications in infectious diseases, rare immunodeficiencies and autoimmune diseases, illustrating advances in diagnosis and understanding risk including use of polygenic risk scores. We further explore the application for patient stratification and drug target prioritization. The review highlights a key challenge to the field arising from the lack of sufficient representation of genetically diverse populations in genomic studies. This currently limits the clinical utility of genetic-based diagnostic and risk-based applications in non-Caucasian populations. We highlight current genome projects, initiatives and biobanks from diverse populations and how this is being used to improve healthcare globally by improving our understanding of genetic susceptibility to diseases and regional pathogens such as malaria and tuberculosis. Future directions and opportunities for personalized medicine and wider application of genomics in health care are described, for the benefit of individual patients and populations worldwide.

Novel insight into the etiology of ischemic stroke gained by integrative multiome-wide association study

Stroke, characterized by sudden neurological deficits, is the second leading cause of death worldwide. Although genome-wide association studies (GWAS) have successfully identified many genomic regions associated with ischemic stroke (IS), the genes underlying risk and their regulatory mechanisms remain elusive. Here, we integrate a large-scale GWAS (N = 1 296 908) for IS together with molecular QTLs data, including mRNA, splicing, enhancer RNA (eRNA), and protein expression data from up to 50 tissues (total N = 11 588). We identify 136 genes/eRNA/proteins associated with IS risk across 60 independent genomic regions and find IS risk is most enriched for eQTLs in arterial and brain-related tissues. Focusing on IS-relevant tissues, we prioritize 9 genes/proteins using probabilistic fine-mapping TWAS analyses. In addition, we discover that blood cell traits, particularly reticulocyte cells, have shared genetic contributions with IS using TWAS-based pheWAS and genetic correlation analysis. Lastly, we integrate our findings with a large-scale pharmacological database and identify a secondary bile acid, deoxycholic acid, as a potential therapeutic component. Our work highlights IS risk genes/splicing-sites/enhancer activity/proteins with their phenotypic consequences using relevant tissues as well as identify potential therapeutic candidates for IS.

Development of a human genetics-guided priority score for 19,365 genes and 399 drug indications

Studies have shown that drug targets with human genetic support are more likely to succeed in clinical trials. Hence, a tool integrating genetic evidence to prioritize drug target genes is beneficial for drug discovery. We built a genetic priority score (GPS) by integrating eight genetic features with drug indications from the Open Targets and SIDER databases. The top 0.83%, 0.28% and 0.19% of the GPS conferred a 5.3-, 9.9- and 11.0-fold increased effect of having an indication, respectively. In addition, we observed that targets in the top 0.28% of the score were 1.7-, 3.7- and 8.8-fold more likely to advance from phase I to phases II, III and IV, respectively. Complementary to the GPS, we incorporated the direction of genetic effect and drug mechanism into a directional version of the score called the GPS with direction of effect. We applied our method to 19,365 protein-coding genes and 399 drug indications and made all results available through a web portal.

Genetic and Gene Expression Resources for Osteoporosis and Bone Biology Research

PURPOSE OF REVIEW

The integration of data from multiple genomic assays from humans and non-human model organisms is an effective approach to identify genes involved in skeletal fragility and fracture risk due to osteoporosis and other conditions. This review summarizes genome-wide genetic variation and gene expression data resources relevant to the discovery of genes contributing to skeletal fragility and fracture risk.

RECENT FINDINGS

Genome-wide association studies (GWAS) of osteoporosis-related traits are summarized, in addition to gene expression in bone tissues in humans and non-human organisms, with a focus on rodent models related to skeletal fragility and fracture risk. Gene discovery approaches using these genomic data resources are described. We also describe the Musculoskeletal Knowledge Portal (MSKKP) that integrates much of the available genomic data relevant to fracture risk. The available genomic resources provide a wealth of knowledge and can be analyzed to identify genes related to fracture risk. Genomic resources that would fill particular scientific gaps are discussed.

Gibson T, McCain S, Stone BC, Gill T, Pickett BE. Pathway2Targets: an open-source pathway-based approach to repurpose therapeutic drugs and prioritize human targets

Background

Recent efforts to repurpose existing drugs to different indications have been accompanied by a number of computational methods, which incorporate protein-protein interaction networks and signaling pathways, to aid with prioritizing existing targets and/or drugs. However, many of these existing methods are focused on integrating additional data that are only available for a small subset of diseases or conditions.

Methods

We have designed and implemented a new R-based open-source target prioritization and repurposing method that integrates both canonical intracellular signaling information from five public pathway databases and target information from public sources including OpenTargets.org. The Pathway2Targets algorithm takes a list of significant pathways as input, then retrieves and integrates public data for all targets within those pathways for a given condition. It also incorporates a weighting scheme that is customizable by the user to support a variety of use cases including target prioritization, drug repurposing, and identifying novel targets that are biologically relevant for a different indication.

Results

As a proof of concept, we applied this algorithm to a public colorectal cancer RNA-sequencing dataset with 144 case and control samples. Our analysis identified 430 targets and ~700 unique drugs based on differential gene expression and signaling pathway enrichment. We found that our highest-ranked predicted targets were significantly enriched in targets with FDA-approved therapeutics for colorectal cancer (p-value < 0.025) that included EGFR, VEGFA, and PTGS2. Interestingly, there was no statistically significant enrichment of targets for other cancers in this same list suggesting high specificity of the results. We also adjusted the weighting scheme to prioritize more novel targets for CRC. This second analysis revealed epidermal growth factor receptor (EGFR), phosphoinositide-3-kinase (PI3K), and two mitogen-activated protein kinases (MAPK14 and MAPK3). These observations suggest that our open-source method with a customizable weighting scheme can accurately prioritize targets that are specific and relevant to the disease or condition of interest, as well as targets that are at earlier stages of development. We anticipate that this method will complement other approaches to repurpose drugs for a variety of indications, which can contribute to the improvement of the quality of life and overall health of such patients.

Prioritization of therapeutic targets for dyslipidemia using integrative multi-omics and multi-trait analysis

Drug targets with genetic support are several-fold more likely to succeed in clinical trials. We introduce a genetic-driven approach based on causal inferences that can inform drug target prioritization, repurposing, and adverse effects of using lipid-lowering agents. Given that a multi-trait approach increases the power to detect meaningful variants/genes, we conduct multi-omics and multi-trait analyses, followed by network connectivity investigations, and prioritize 30 potential therapeutic targets for dyslipidemia, including SORT1, PSRC1, CELSR2, PCSK9, HMGCR, APOB, GRN, HFE2, FJX1, C1QTNF1, and SLC5A8. 20% (6/30) of prioritized targets from our hypothesis-free drug target search are either approved or under investigation for dyslipidemia. The prioritized targets are 22-fold higher in likelihood of being approved or under investigation in clinical trials than genome-wide association study (GWAS)-curated targets. Our results demonstrate that the genetic-driven approach used in this study is a promising strategy for prioritizing targets while informing about the potential adverse effects and repurposing opportunities.

Priority index for asthma (PIA): In silico discovery of shared and distinct drug targets for adult- and childhood-onset disease

Asthma is a chronic disease that is caused by a combination of genetic risks and environmental triggers and can affect both adults and children. Genome-wide association studies have revealed partly distinct genetic architectures for its two age-of-onset subtypes (namely, adult-onset and childhood-onset). We reason that identifying shared and distinct drug targets between these subtypes may inform the development of subtype-specific therapeutic strategies. In attempting this, we here introduce Priority Index for Asthma or PIA, a genetics-led and network-driven drug target prioritisation tool for asthma. We demonstrate the validity of the tool in improving drug target prioritisation for asthma compared to the status quo methods, as well as in capturing the underlying etiology and existing therapeutics for the disease. We also illustrate how PIA can be used to prioritise drug targets for adult- and childhood-onset asthma, as well as to identify shared and distinct pathway crosstalk genes. Shared crosstalk genes are mostly involved in JAK-STAT signaling, with clinical evidence supporting that targeting this pathway may be a promising drug repurposing opportunity for both subtypes. Crosstalk genes specific to childhood-onset asthma are enriched for PI3K-AKT-mTOR signaling, and we identify genes that are already targeted by licensed medications as repurposed drug candidates for this subtype. We make all our results accessible and reproducible at http://www.genetictargets.com/PIA. Collectively, our study has significant implications for asthma computational medicine research and can guide the future development of subtype-specific therapeutic strategies for the disease.

The dcGO Domain-Centric Ontology Database in 2023: New Website and Extended Annotations for Protein Structural Domains

Protein structural domains have been less studied than full-length proteins in terms of ontology annotations. The dcGO database has filled this gap by providing mappings from protein domains to ontologies. The dcGO update in 2023 extends annotations for protein domains of multiple definitions (SCOP, Pfam, and InterPro) with commonly used ontologies that are categorised into functions, phenotypes, diseases, drugs, pathways, regulators, and hallmarks. This update adds new dimensions to the utility of both ontology and protein domain resources. A newly designed website at http://www.protdomainonto.pro/dcGO offers a more centralised and user-friendly way to access the dcGO database, with enhanced faceted search returning term- and domain-specific information pages. Users can navigate both ontology terms and annotated domains through improved ontology hierarchy browsing. A newly added facility enables domain-based ontology enrichment analysis.

Gene, cell type, and drug prioritization analysis suggest genetic basis for the utility of diuretics in treating Alzheimer disease

We introduce a user-friendly tool for risk gene, cell type, and drug prioritization for complex traits: GCDPipe. It uses gene-level GWAS-derived data and gene expression data to train a model for the identification of disease risk genes and relevant cell types. Gene prioritization information is then coupled with known drug target data to search for applicable drug agents based on their estimated functional effects on the identified risk genes. We illustrate the utility of our approach in different settings: identification of the cell types, implicated in disease pathogenesis, was tested in inflammatory bowel disease (IBD) and Alzheimer disease (AD); gene target and drug prioritization was tested in IBD and schizophrenia. The analysis of phenotypes with known disease-affected cell types and/or existing drug candidates shows that GCDPipe is an effective tool to unify genetic risk factors with cellular context and known drug targets. Next, analysis of the AD data with GCDPipe suggested that gene targets of diuretics, as an Anatomical Therapeutic Chemical drug subgroup, are significantly enriched among the genes prioritized by GCDPipe, indicating their possible effect on the course of the disease.

Multi-layered genetic approaches to identify approved drug targets

Drugs targeting genes linked to disease via evidence from human genetics have increased odds of approval. Approaches to prioritize such genes include genome-wide association studies (GWASs), rare variant burden tests in exome sequencing studies (Exome), or integration of a GWAS with expression/protein quantitative trait loci (eQTL/pQTL-GWAS). Here, we compare gene-prioritization approaches on 30 clinically relevant traits and benchmark their ability to recover drug targets. Across traits, prioritized genes were enriched for drug targets with odds ratios (ORs) of 2.17, 2.04, 1.81, and 1.31 for the GWAS, eQTL-GWAS, Exome, and pQTL-GWAS methods, respectively. Adjusting for differences in testable genes and sample sizes, GWAS outperforms e/pQTL-GWAS, but not the Exome approach. Furthermore, performance increased through gene network diffusion, although the node degree, being the best predictor (OR = 8.7), revealed strong bias in literature-curated networks. In conclusion, we systematically assessed strategies to prioritize drug target genes, highlighting the promises and pitfalls of current approaches.

Molecular profiling reveals potential targets in cholangiocarcinoma

BACKGROUND

Cholangiocarcinoma (CCA) is a devastating malignancy and has a very poor prognosis if tumors spread outside the liver. Understanding the molecular mechanisms underlying the CCA progression will likely yield therapeutic approaches toward treating this deadly disease.

AIM

To determine the molecular pathogenesis in CCA progression.

METHODS

In silico analysis, in vitro cell culture, CCA transgenic animals, histological, and molecular assays were adopted to determine the molecular pathogenesis.

RESULTS

The transcriptomic data of human CCA samples were retrieved from The Cancer Genome Atlas (TGCA, CHOL), European Bioinformatics Institute (EBI, GAD00001001076), and Gene Expression Omnibus (GEO, GSE107943) databases. Using Gene set enrichment analysis, the cell cycle and Notch related pathways were demonstrated to be significantly activated in CCA in TCGA and GEO datasets. We, through differentially expressed genes, found several cell cycle and notch associated genes were significantly up-regulated in cancer tissues when compared with the non-cancerous control samples. The associated genes, via quantitative real-time PCR and western blotting assays, were further examined in normal human cholangiocytes, CCA cell lines, mouse normal bile ducts, and mouse CCA tumors established by specifically depleting P53 and expressing KrasG12D mutation in the liver. Consistently, we validated that the cell cycle and Notch pathways are up-regulated in CCA cell lines and mouse CCA tumors. Interestingly, targeting cell cycle and notch pathways using small molecules also exhibited significant beneficial effects in controlling tumor malignancy. More importantly, we demonstrated that several cell cycle and Notch associated genes are significantly associated with poor overall survival and disease-free survival using the Log-Rank test.

CONCLUSION

In summary, our study comprehensively analyzed the gene expression pattern of CCA samples using publicly available datasets and identified the cell cycle and Notch pathways are potential therapeutic targets in this deadly disease.

Genome-wide association analyses define pathogenic signaling pathways and prioritize drug targets for IgA nephropathy

IgA nephropathy (IgAN) is a progressive form of kidney disease defined by glomerular deposition of IgA. Here we performed a genome-wide association study of 10,146 kidney-biopsy-diagnosed IgAN cases and 28,751 controls across 17 international cohorts. We defined 30 genome-wide significant risk loci explaining 11% of disease risk. A total of 16 loci were new, including TNFSF4/TNFSF18, REL, CD28, PF4V1, LY86, LYN, ANXA3, TNFSF8/TNFSF15, REEP3, ZMIZ1, OVOL1/RELA, ETS1, IGH, IRF8, TNFRSF13B and FCAR. The risk loci were enriched in gene orthologs causing abnormal IgA levels when genetically manipulated in mice. We also observed a positive genetic correlation between IgAN and serum IgA levels. High polygenic score for IgAN was associated with earlier onset of kidney failure. In a comprehensive functional annotation analysis of candidate causal genes, we observed convergence of biological candidates on a common set of inflammatory signaling pathways and cytokine ligand-receptor pairs, prioritizing potential new drug targets.

A human proteogenomic-cellular framework identifies KIF5A as a modulator of astrocyte process integrity with relevance to ALS

Genome-wide association studies identified several disease-causing mutations in neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS). However, the contribution of genetic variants to pathway disturbances and their cell type-specific variations, especially in glia, is poorly understood. We integrated ALS GWAS-linked gene networks with human astrocyte-specific multi-omics datasets to elucidate pathognomonic signatures. It predicts that KIF5A, a motor protein kinesin-1 heavy-chain isoform, previously detected only in neurons, can also potentiate disease pathways in astrocytes. Using postmortem tissue and super-resolution structured illumination microscopy in cell-based perturbation platforms, we provide evidence that KIF5A is present in astrocyte processes and its deficiency disrupts structural integrity and mitochondrial transport. We show that this may underly cytoskeletal and trafficking changes in SOD1 ALS astrocytes characterised by low KIF5A levels, which can be rescued by c-Jun N-terminal Kinase-1 (JNK1), a kinesin transport regulator. Altogether, our pipeline reveals a mechanism controlling astrocyte process integrity, a pre-requisite for synapse maintenance and suggests a targetable loss-of-function in ALS.

Comprehensive epigenomic profiling reveals the extent of disease-specific chromatin states and informs target discovery in ankylosing spondylitis

Ankylosing spondylitis (AS) is a common, highly heritable inflammatory arthritis characterized by enthesitis of the spine and sacroiliac joints. Genome-wide association studies (GWASs) have revealed more than 100 genetic associations whose functional effects remain largely unresolved. Here, we present a comprehensive transcriptomic and epigenomic map of disease-relevant blood immune cell subsets from AS patients and healthy controls. We find that, while CD14⁺ monocytes and CD4⁺ and CD8⁺ T cells show disease-specific differences at the RNA level, epigenomic differences are only apparent upon multi-omics integration. The latter reveals enrichment at disease-associated loci in monocytes. We link putative functional SNPs to genes using high-resolution Capture-C at 10 loci, including PTGER4 and ETS1, and show how disease-specific functional genomic data can be integrated with GWASs to enhance therapeutic target discovery. This study combines epigenetic and transcriptional analysis with GWASs to identify disease-relevant cell types and gene regulation of likely pathogenic relevance and prioritize drug targets.

AI/ML in Precision Medicine: A Look Beyond the Hype

Artificial Intelligence (AI) and Machine Learning (ML) are making headlines in medical research, especially in drug discovery, digital imaging, disease diagnostics, genetic testing, and optimal care pathway (personalized care). However, the potential uses and benefits of AI/ML applications need to be distinguished from hype. In the 2022 American Statistical Association Biopharmaceutical Section Regulatory-Industry Statistical Workshop, we convened a panel of experts from FDA and industry to talk about the challenges of successfully applying AI/ML in precision medicine and how to overcome those challenges. This paper provides a summary and expansion on the topics discussed in the panel: the application of AI/ML, bias, and data quality.

A publication-wide association study (PWAS), historical language models to prioritise novel therapeutic drug targets

Most biomedical knowledge is published as text, making it challenging to analyse using traditional statistical methods. In contrast, machine-interpretable data primarily comes from structured property databases, which represent only a fraction of the knowledge present in the biomedical literature. Crucial insights and inferences can be drawn from these publications by the scientific community. We trained language models on literature from different time periods to evaluate their ranking of prospective gene-disease associations and protein-protein interactions. Using 28 distinct historical text corpora of abstracts published between 1995 and 2022, we trained independent Word2Vec models to prioritise associations that were likely to be reported in future years. This study demonstrates that biomedical knowledge can be encoded as word embeddings without the need for human labelling or supervision. Language models effectively capture drug discovery concepts such as clinical tractability, disease associations, and biochemical pathways. Additionally, these models can prioritise hypotheses years before their initial reporting. Our findings underscore the potential for extracting yet-to-be-discovered relationships through data-driven approaches, leading to generalised biomedical literature mining for potential therapeutic drug targets. The Publication-Wide Association Study (PWAS) enables the prioritisation of under-explored targets and provides a scalable system for accelerating early-stage target ranking, irrespective of the specific disease of interest.

Novel insight into the etiology of ischemic stroke gained by integrative transcriptome-wide association study

Stroke, characterized by sudden neurological deficits, is the second leading cause of death worldwide. Although genome-wide association studies (GWAS) have successfully identified many genomic regions associated with ischemic stroke (IS), the genes underlying risk and their regulatory mechanisms remain elusive. Here, we integrate a large-scale GWAS (N=1,296,908) for IS together with mRNA, splicing, enhancer RNA (eRNA) and protein expression data (N=11,588) from 50 tissues. We identify 136 genes/eRNA/proteins associated with IS risk across 54 independent genomic regions and find IS risk is most enriched for eQTLs in arterial and brain-related tissues. Focusing on IS-relevant tissues, we prioritize 9 genes/proteins using probabilistic fine-mapping TWAS analyses. In addition, we discover that blood cell traits, particularly reticulocyte cells, have shared genetic contributions with IS using TWAS-based pheWAS and genetic correlation analysis. Lastly, we integrate our findings with a large-scale pharmacological database and identify a secondary bile acid, deoxycholic acid, as a potential therapeutic component. Our work highlights IS risk genes/splicing-sites/enhancer activity/proteins with their phenotypic consequences using relevant tissues as well as identify potential therapeutic candidates for IS.

The long noncoding RNA THBS1-AS1 promotes cardiac fibroblast activation in cardiac fibrosis by regulating TGFBR1

Cardiac fibrosis is associated with an adverse prognosis in cardiovascular disease that results in a decreased cardiac compliance and, ultimately, heart failure. Recent studies have identified the role of long noncoding RNA (lncRNA) in cardiac fibrosis. However, the functions of many lncRNAs in cardiac fibrosis remain to be characterized. Through a whole-transcriptome sequencing and bioinformatics analysis on a mouse model of pressure overload-induced cardiac fibrosis, we screened a key lncRNA termed thrombospondin 1 antisense 1 (THBS1-AS1), which was positively associated with cardiac fibrosis. In vitro functional studies demonstrated that the silencing of THBS1-AS1 ameliorated TGF-β1 effects on cardiac fibroblast (CF) activation, and the overexpression of THBS1-AS1 displayed the opposite effect. A mechanistic study revealed that THBS1-AS1 could sponge miR-221/222 to regulate the expression of TGFBR1. Moreover, under TGF-β1 stimulation, the forced expression of miR-221/222 or the knockdown TGFBR1 significantly reversed the THBS1-AS1 overexpression induced by further CF activation. In vivo, specific knockdown of THBS1-AS1 in activated CFs significantly alleviated transverse aorta constriction-induced (TAC-induced) cardiac fibrosis in mice. Finally, we demonstrated that the human THBS1-AS1 can also affect the activation of CFs by regulating TGFBR1. In conclusion, this study reveals that lncRNA THBS1-AS1 is a potentially novel regulator of cardiac fibrosis and may serve as a target for the treatment of cardiac fibrosis.

Network expansion of genetic associations defines a pleiotropy map of human cell biology

Interacting proteins tend to have similar functions, influencing the same organismal traits. Interaction networks can be used to expand the list of candidate trait-associated genes from genome-wide association studies. Here, we performed network-based expansion of trait-associated genes for 1,002 human traits showing that this recovers known disease genes or drug targets. The similarity of network expansion scores identifies groups of traits likely to share an underlying genetic and biological process. We identified 73 pleiotropic gene modules linked to multiple traits, enriched in genes involved in processes such as protein ubiquitination and RNA processing. In contrast to gene deletion studies, pleiotropy as defined here captures specifically multicellular-related processes. We show examples of modules linked to human diseases enriched in genes with known pathogenic variants that can be used to map targets of approved drugs for repurposing. Finally, we illustrate the use of network expansion scores to study genes at inflammatory bowel disease genome-wide association study loci, and implicate inflammatory bowel disease-relevant genes with strong functional and genetic support.

Using human genetics to improve safety assessment of therapeutics

Human genetics research has discovered thousands of proteins associated with complex and rare diseases. Genome-wide association studies (GWAS) and studies of Mendelian disease have resulted in an increased understanding of the role of gene function and regulation in human conditions. Although the application of human genetics has been explored primarily as a method to identify potential drug targets and support their relevance to disease in humans, there is increasing interest in using genetic data to identify potential safety liabilities of modulating a given target. Human genetic variants can be used as a model to anticipate the effect of lifelong modulation of therapeutic targets and identify the potential risk for on-target adverse events. This approach is particularly useful for non-clinical safety evaluation of novel therapeutics that lack pharmacologically relevant animal models and can contribute to the intrinsic safety profile of a drug target. This Review illustrates applications of human genetics to safety studies during drug discovery and development, including assessing the potential for on- and off-target associated adverse events, carcinogenicity risk assessment, and guiding translational safety study designs and monitoring strategies. A summary of available human genetic resources and recommended best practices is provided. The challenges and future perspectives of translating human genetic information to identify risks for potential drug effects in preclinical and clinical development are discussed.

Inhibition of long noncoding RNA Gm41724 alleviates pressure overload-induced cardiac fibrosis by regulating lamina-associated polypeptide 2α

Cardiac fibrosis is a pathological process underlying myocardial remodeling and is characterized by excessive deposition of the myocardial extracellular matrix. Long noncoding RNAs (lncRNAs) have emerged as critical regulators of various biological processes. In this study, we investigated the role of a novel lncRNA, Gm41724, in cardiac fibrosis induced by pressure overload. High-throughput whole transcriptome sequencing analysis was performed to detect differentially expressed lncRNAs in cardiac fibroblasts (CFs) with or without TGF-β1 treatment. Differential expression analysis and gene set enrichment analysis identified Gm41724 as a potential molecule targeting fibrosis. Gm41724 positively regulated the activation of CFs induced by TGF-β1 and pressure overload. Knocking down Gm41724 could inhibit the differentiation of CFs into myofibroblasts and alleviate cardiac fibrosis induced by pressure overload. Mechanistically, comprehensive identification of RNA-binding proteins by mass spectrometry (CHIRP-MS) and RNA immunoprecipitation (RIP) assay combined with other methods of molecular biological revealed the important role of Gm41724 binding to lamina-associated polypeptide 2α (lap2α) for the activation of CFs. Further mechanistic studies indicated that the regulator of G protein signaling 4 (Rgs4), as the downstream effector of Gm41724/lap2α, regulated CFs activation. Our results implicated the involvement of Gm41724 in cardiac fibrosis induced by pressure overload and it is expected to be a promising target for anti-fibrotic therapy.

When cancer drug resistance meets metabolomics (bulk, single-cell and/or spatial): Progress, potential, and perspective

Resistance to drug treatment is a critical barrier in cancer therapy. There is an unmet need to explore cancer hallmarks that can be targeted to overcome this resistance for therapeutic gain. Over time, metabolic reprogramming has been recognised as one hallmark that can be used to prevent therapeutic resistance. With the advent of metabolomics, targeting metabolic alterations in cancer cells and host patients represents an emerging therapeutic strategy for overcoming cancer drug resistance. Driven by technological and methodological advances in mass spectrometry imaging, spatial metabolomics involves the profiling of all the metabolites (metabolomics) so that the spatial information is captured bona fide within the sample. Spatial metabolomics offers an opportunity to demonstrate the drug-resistant tumor profile with metabolic heterogeneity, and also poses a data-mining challenge to reveal meaningful insights from high-dimensional spatial information. In this review, we discuss the latest progress, with the focus on currently available bulk, single-cell and spatial metabolomics technologies and their successful applications in pre-clinical and translational studies on cancer drug resistance. We provide a summary of metabolic mechanisms underlying cancer drug resistance from different aspects; these include the Warburg effect, altered amino acid/lipid/drug metabolism, generation of drug-resistant cancer stem cells, and immunosuppressive metabolism. Furthermore, we propose solutions describing how to overcome cancer drug resistance; these include early detection during cancer initiation, monitoring of clinical drug response, novel anticancer drug and target metabolism, immunotherapy, and the emergence of spatial metabolomics. We conclude by describing the perspectives on how spatial omics approaches (integrating spatial metabolomics) could be further developed to improve the management of drug resistance in cancer patients.

Single-cell-led drug repurposing for Alzheimer's disease

Alzheimer's disease is the most common form of dementia. Notwithstanding the huge investments in drug development, only one disease-modifying treatment has been recently approved. Here we present a single-cell-led systems biology pipeline for the identification of drug repurposing candidates. Using single-cell RNA sequencing data of brain tissues from patients with Alzheimer's disease, genome-wide association study results, and multiple gene annotation resources, we built a multi-cellular Alzheimer's disease molecular network that we leveraged for gaining cell-specific insights into Alzheimer's disease pathophysiology and for the identification of drug repurposing candidates. Our computational approach pointed out 54 candidate drugs, mainly targeting MAPK and IGF1R signaling pathways, which could be further evaluated for their potential as Alzheimer's disease therapy.

Network analysis of genome-wide association studies for drug target prioritisation

Over the past decades, genome-wide association studies (GWAS) have led to a dramatic expansion of genetic variants implicated with human traits and diseases. These advances are expected to result in new drug targets but the identification of causal genes and the cell biology underlying human diseases from GWAS remains challenging. Here, we review protein interaction network-based methods to analyse GWAS data. These approaches can rank candidate drug targets at GWAS-associated loci or among interactors of disease genes without direct genetic support. These methods identify the cell biology affected in common across diseases, offering opportunities for drug repurposing, as well as be combined with expression data to identify focal tissues and cell types. Going forward, we expect that these methods will further improve from advances in the characterisation of context specific interaction networks and the joint analysis of rare and common genetic signals.

SLC11A1 as a stratification indicator for immunotherapy or chemotherapy in patients with glioma

Background

Glioma is a fatal tumor originating from the brain, which accounts for most intracranial malignancies. Currently, Immunotherapy has turned into a novel and promising treatment in glioma patients. however, there are still few effective biomarkers to mirror the reaction to immunotherapy in patients with glioma. Therefore, we intended to elucidate the evaluable efficacy of SLC11A1 in glioma patients.

Methods

In this study, samples from Shanghai General Hospital and data from TCGA, GEO, CGGA datasets were used to investigate and validate the relationship between SLC11A1 and the progression of glioma. We evaluated the predictive value of SLC11A1 on the prognosis of glioma with cox regression analysis. Then the relationship between immune infiltration and SLC11A1 was also analyzed. Ultimately, we performed the prediction on the immunotherapeutic response and therapeutic drugs according to the expression of SLC11A1.

Results

Expression of SLC11A1 increased with progression and predicted unfavorable prognosis for glioma patients. The hazard ratio for SLC11A1 expression was 2.33 with 95% CI (1.92-2.58) (P < 0.001) in cox analysis. And based on expression, we found SLC11A1 stratified glioma patients into subgroups with different immune activation statuses. Moreover, we observed that patients with higher SLC11A1 levels companied with better immunotherapeutic response, while those with lower SLC11A1 levels may respond better to temozolomide.

Conclusion

This study provided evidence that SLC11A1 was a novel prognostic marker and immunotherapy response indicator for gliomas. In some cases, SLC11A1 could be an effective marker for identifying patients who might benefit from immunotherapy or chemotherapy.

QTLbase2: an enhanced catalog of human quantitative trait loci on extensive molecular phenotypes

Deciphering the fine-scale molecular mechanisms that shape the genetic effects at disease-associated loci from genome-wide association studies (GWAS) remains challenging. The key avenue is to identify the essential molecular phenotypes that mediate the causal variant and disease under particular biological conditions. Therefore, integrating GWAS signals with context-specific quantitative trait loci (QTLs) (such as different tissue/cell types, disease states, and perturbations) from extensive molecular phenotypes would present important strategies for full understanding of disease genetics. Via persistent curation and systematic data processing of large-scale human molecular trait QTLs (xQTLs), we updated our previous QTLbase database (now QTLbase2, http://mulinlab.org/qtlbase) to comprehensively analyze and visualize context-specific QTLs across 22 molecular phenotypes and over 95 tissue/cell types. Overall, the resource features the following major updates and novel functions: (i) 960 more genome-wide QTL summary statistics from 146 independent studies; (ii) new data for 10 previously uncompiled QTL types; (iii) variant query scope expanded to fit 195 QTL datasets based on whole-genome sequencing; (iv) supports filtering and comparison of QTLs for different biological conditions, such as stimulation types and disease states; (v) a new linkage disequilibrium viewer to facilitate variant prioritization across tissue/cell types and QTL types.

Integrative analyses of immune-related biomarkers and associated mechanisms in coronary heart disease

Various studies showed that the effect of immune activation is pro-atherogenic and coronary heart disease (CHD) should therefore be considered an autoimmune disease. This study aimed to identify potential immune-related biomarkers, pathways, and the potential regulatory networks underlying CHD. Differentially expressed genes (DEGs) between CHD and control samples were determined by analyzing GSE71226 and GSE9128. The overlapping differential expression immune-related genes (DE-IRGs) for CHD were identified by analyzing the ImmPort database and two GEO databases. A total of 384 DE-IRGs were identified. Subsequently, comprehensive enrichment analyses suggested that DE-IRGs were enriched in immune-related pathways, including autoimmune thyroid disease, the intestinal immune network for IGA production, and downstream signaling events of B cell receptors. The signature of DE-IRGs was validated using an external independent dataset GSE20681 (AUC = 0.875). Furthermore, we conducted protein–protein interaction network analysis and identified eight hub genes, which were most enriched in regulation of defense response, NF-κB signaling pathway, regulation of JNK cascade, and regulation of cytokine production. Moreover, networks of miRNAs-mRNAs and transcription factors (TFs)-mRNA underlying the integrated data were established, involving eight miRNAs and 76 TF-targeting hub genes. Ultimately, 17 SNPs in miRNA-mediated gene networks were identified. We screened potential immune-related genes in CHD and constructed miRNA-mRNA-TF and SNP-miRNA networks, which not only provide inspired insights into the occurrence and the molecular mechanisms of CHD but also lay a foundation for targeting potential biomarkers using immunotherapy and for understanding the molecular mechanisms of CHD.

Scalable approaches for functional analyses of whole-genome sequencing non-coding variants

Non-coding genetic variants outside of protein-coding genome regions play an important role in genetic and epigenetic regulation. It has become increasingly important to understand their roles, as non-coding variants often make up the majority of top findings of genome-wide association studies (GWAS). In addition, the growing popularity of disease-specific whole-genome sequencing (WGS) efforts expands the library of and offers unique opportunities for investigating both common and rare non-coding variants, which are typically not detected in more limited GWAS approaches. However, the sheer size and breadth of WGS data introduce additional challenges to predicting functional impacts in terms of data analysis and interpretation. This review focuses on the recent approaches developed for efficient, at-scale annotation and prioritization of non-coding variants uncovered in WGS analyses. In particular, we review the latest scalable annotation tools, databases and functional genomic resources for interpreting the variant findings from WGS based on both experimental data and in silico predictive annotations. We also review machine learning-based predictive models for variant scoring and prioritization. We conclude with a discussion of future research directions which will enhance the data and tools necessary for the effective functional analyses of variants identified by WGS to improve our understanding of disease etiology.

DrugMAP: molecular atlas and pharma-information of all drugs

The efficacy and safety of drugs are widely known to be determined by their interactions with multiple molecules of pharmacological importance, and it is therefore essential to systematically depict the molecular atlas and pharma-information of studied drugs. However, our understanding of such information is neither comprehensive nor precise, which necessitates the construction of a new database providing a network containing a large number of drugs and their interacting molecules. Here, a new database describing the molecular atlas and pharma-information of drugs (DrugMAP) was therefore constructed. It provides a comprehensive list of interacting molecules for >30 000 drugs/drug candidates, gives the differential expression patterns for >5000 interacting molecules among different disease sites, ADME (absorption, distribution, metabolism and excretion)-relevant organs and physiological tissues, and weaves a comprehensive and precise network containing >200 000 interactions among drugs and molecules. With the great efforts made to clarify the complex mechanism underlying drug pharmacokinetics and pharmacodynamics and rapidly emerging interests in artificial intelligence (AI)-based network analyses, DrugMAP is expected to become an indispensable supplement to existing databases to facilitate drug discovery. It is now fully and freely accessible at: https://idrblab.org/drugmap/.

A practical guideline of genomics-driven drug discovery in the era of global biobank meta-analysis

Genomics-driven drug discovery is indispensable for accelerating the development of novel therapeutic targets. However, the drug discovery framework based on evidence from genome-wide association studies (GWASs) has not been established, especially for cross-population GWAS meta-analysis. Here, we introduce a practical guideline for genomics-driven drug discovery for cross-population meta-analysis, as lessons from the Global Biobank Meta-analysis Initiative (GBMI). Our drug discovery framework encompassed three methodologies and was applied to the 13 common diseases targeted by GBMI (N _mean = 1,329,242). Individual methodologies complementarily prioritized drugs and drug targets, which were systematically validated by referring previously known drug-disease relationships. Integration of the three methodologies provided a comprehensive catalog of candidate drugs for repositioning, nominating promising drug candidates targeting the genes involved in the coagulation process for venous thromboembolism and the interleukin-4 and interleukin-13 signaling pathway for gout. Our study highlighted key factors for successful genomics-driven drug discovery using cross-population meta-analyses.

Identifying disease-critical cell types and cellular processes by integrating single-cell RNA-sequencing and human genetics

Genome-wide association studies provide a powerful means of identifying loci and genes contributing to disease, but in many cases, the related cell types/states through which genes confer disease risk remain unknown. Deciphering such relationships is important for identifying pathogenic processes and developing therapeutics. In the present study, we introduce sc-linker, a framework for integrating single-cell RNA-sequencing, epigenomic SNP-to-gene maps and genome-wide association study summary statistics to infer the underlying cell types and processes by which genetic variants influence disease. The inferred disease enrichments recapitulated known biology and highlighted notable cell-disease relationships, including γ-aminobutyric acid-ergic neurons in major depressive disorder, a disease-dependent M-cell program in ulcerative colitis and a disease-specific complement cascade process in multiple sclerosis. In autoimmune disease, both healthy and disease-dependent immune cell-type programs were associated, whereas only disease-dependent epithelial cell programs were prominent, suggesting a role in disease response rather than initiation. Our framework provides a powerful approach for identifying the cell types and cellular processes by which genetic variants influence disease.

Recent advances in elucidating the genetic basis of systemic sclerosis

PURPOSE OF REVIEW

Systemic sclerosis (SSc) is a complex autoimmune disorder that affects the connective tissue and causes severe vascular damage and fibrosis of the skin and internal organs. There are recent advances in the field that apply novel methods to high throughput genotype information of thousands of patients with SSc and provide promising results towards the use of genomic data to help SSc diagnosis and clinical care.

RECENT FINDINGS

This review addresses the development of the first SSc genomic risk score, which can contribute to differentiating SSc patients from healthy controls and other immune-mediated diseases. Moreover, we explore the implementation of data mining strategies on the results of genome-wide association studies to highlight subtype-specific HLA class II associations and a strong association of the HLA class I locus with SSc for the first time. Finally, the combination of genomic data with transcriptomics informed drug repurposing and genetic association studies in well characterized SSc patient cohorts identified markers of severe complications of the disease.

SUMMARY

Early diagnosis and clinical management of SSc and SSc-related complications are still challenging for rheumatologists. The development of predictive models and tools using genotype data may help to finally deliver personalized clinical care and treatment for patients with SSc in the near future.

SNP-to-gene linking strategies reveal contributions of enhancer-related and candidate master-regulator genes to autoimmune disease

We assess contributions to autoimmune disease of genes whose regulation is driven by enhancer regions (enhancer-related) and genes that regulate other genes in trans (candidate master-regulator). We link these genes to SNPs using several SNP-to-gene (S2G) strategies and apply heritability analyses to draw three conclusions about 11 autoimmune/blood-related diseases/traits. First, several characterizations of enhancer-related genes using functional genomics data are informative for autoimmune disease heritability after conditioning on a broad set of regulatory annotations. Second, candidate master-regulator genes defined using trans-eQTL in blood are also conditionally informative for autoimmune disease heritability. Third, integrating enhancer-related and master-regulator gene sets with protein-protein interaction (PPI) network information magnified their disease signal. The resulting PPI-enhancer gene score produced >2-fold stronger heritability signal and >2-fold stronger enrichment for drug targets, compared with the recently proposed enhancer domain score. In each case, functionally informed S2G strategies produced 4.1- to 13-fold stronger disease signals than conventional window-based strategies.

PiER: web-based facilities tailored for genetic target prioritisation harnessing human disease genetics, functional genomics and protein interactions

Integrative prioritisation promotes translational use of disease genetic findings in target discovery. I report 'PiER' (http://www.genetictargets.com/PiER), web-based facilities that support ab initio and real-time genetic target prioritisation through integrative use of human disease genetics, functional genomics and protein interactions. By design, the PiER features two facilities: elementary and combinatory. The elementary facility is designed to perform specific tasks, including three online tools: eV2CG, utilising functional genomics to link disease-associated variants (particularly located at the non-coding genome) to core genes likely responsible for genetic associations in disease; eCG2PG, using knowledge of protein interactions to 'network' core genes and additional peripheral genes, producing a ranked list of core and peripheral genes; and eCrosstalk, exploiting the information of pathway-derived interactions to identify highly-ranked genes mediating crosstalk between molecular pathways. Each of elementary tasks giving results is sequentially piped to the next one. By chaining together elementary tasks, the combinatory facility automates genetics-led and network-based integrative prioritisation for genetic targets at the gene level (cTGene) and at the crosstalk level (cTCrosstalk). Together with a tutorial-like booklet describing instructions on how to use, the PiER facilities meet multi-tasking needs to accelerate computational translational medicine that leverages human disease genetics and genomics for early-stage target discovery and drug repurposing.