Demystifying non-coding GWAS variants: an overview of computational tools and methods

Using Genomics to Develop Personalized Cardiovascular Treatments

Advances in genomic technologies have significantly enhanced our understanding of both monogenic and polygenic etiologies of cardiovascular disease. In this review, we explore how the utilization of genomic information is bringing personalized medicine approaches to the forefront of cardiovascular disease management. We describe how genomic data can resolve diagnostic uncertainty, support cascade screening, and inform treatment strategies. We discuss how genome-wide association studies have identified thousands of genetic variants associated with polygenic cardiovascular diseases, and how integrating these insights into polygenic risk scores can enhance personalized risk prediction beyond traditional clinical algorithms. We detail how pharmacogenomics approaches leverage genotype information to guide drug selection and mitigate adverse events. Finally, we present the paradigm-shifting approach of gene therapy, which holds the promise of being a curative intervention for cardiovascular conditions.

Identification of Candidate Genes and eQTLs Related to Porcine Reproductive Function

Expression quantitative trait locus (eQTL) mapping is an effective tool for identifying genetic variations that regulate gene expression. An increasing number of studies suggested that SNPs associated with complex traits in farm animals are considered as expression quantitative trait loci. Identifying eQTLs associated with gene expression levels in the endometrium helps to unravel the regulatory mechanisms of genes related to reproductive functions in this tissue and provides molecular markers for the genetic improvement of high-fertility sow breeding. In this study, 218 RNA-seq data from pig endometrial tissue were used for eQTL analysis to identify genetic variants regulating gene expression. Additionally, weighted gene co-expression network analysis (WGCNA) was performed to identify hub genes involved in reproductive functions. The eQTL analysis identified 34,876 significant cis-eQTLs regulating the expression of 5632 genes (FDR ≤ 0.05), and 90 hub genes were identified by WGCNA analysis. By integrating eQTL and WGCNA results, 14 candidate genes and 16 fine-mapped cis-eQTLs were identified, including FRK, ARMC3, SLC35F3, TMEM72, FFAR4, SOWAHA, PSPH, FMO5, HPN, FUT2, RAP1GAP, C6orf52, SEL1L3, and CLGN, which were involved in the physiological processes of reproduction in sows through hormone regulation, cell adhesion, and amino acid and lipid metabolism. These eQTLs regulate the high expression of candidate genes in the endometrium, thereby affecting reproductive-related physiological functions. These findings enhance our understanding of the genetic basis of reproductive traits and provide valuable genetic markers for marker-assisted selection (MAS), which can be applied to improve sow fecundity and optimize breeding strategies for high reproductive performance.

Implications of gene × environment interactions in post-traumatic stress disorder risk and treatment

Exposure to traumatic stress is common in the general population. Variation in the brain's molecular encoding of stress potentially contributes to the heterogeneous clinical outcomes in response to traumatic experiences. For instance, only a minority of those exposed to trauma will develop post-traumatic stress disorder (PTSD). Risk for PTSD is at least partially heritable, with a growing number of genetic factors identified through GWAS. A major limitation of genetic studies is that they capture only the genetic component of risk, whereas PTSD by definition requires an environmental traumatic exposure. Furthermore, the extent, timing, and type of trauma affects susceptibility. Here, we discuss the molecular mechanisms of PTSD risk together with gene × environment interactions, with a focus on how either might inform genetic screening for individuals at high risk for disease, reveal biological mechanisms that might one day yield novel therapeutics, and impact best clinical practices even today. To close, we discuss the interaction of trauma with sex, gender, and race, with a focus on the implications for treatment. Altogether, we suggest that predicting, preventing, and treating PTSD will require integrating both genotypic and environmental information.

The Interplay of Genetic Predisposition, Circadian Misalignment, and Metabolic Regulation in Obesity

PURPOSE OF REVIEW

This review explores the complex interplay between genetic predispositions to obesity, circadian rhythms, metabolic regulation, and sleep. It highlights how genetic factors underlying obesity exacerbate metabolic dysfunction through circadian misalignment and examines promising interventions to mitigate these effects.

RECENT FINDINGS

Genome-wide association Studies (GWAS) have identified numerous Single Nucleotide Polymorphisms (SNPs) associated with obesity traits, attributing 40-75% heritability to body mass index (BMI). These findings illuminate critical links between genetic obesity, circadian clocks, and metabolic processes. SNPs in clock-related genes influence metabolic pathways, with disruptions in circadian rhythms-driven by poor sleep hygiene or erratic eating patterns-amplifying metabolic dysfunction. Circadian clocks, synchronized with the 24-h light-dark cycle, regulate key metabolic activities, including glucose metabolism, lipid storage, and energy utilization. Genetic mutations or external disruptions, such as irregular sleep or eating habits, can destabilize circadian rhythms, promoting weight gain and metabolic disorders. Circadian misalignment in individuals with genetic predispositions to obesity disrupts the release of key metabolic hormones, such as leptin and insulin, impairing hunger regulation and fat storage. Interventions like time-restricted feeding (TRF) and structured physical activity offer promising strategies to restore circadian harmony, improve metabolic health, and mitigate obesity-related risks.

Polygenic Score for the Prediction of Postoperative Nausea and Vomiting: A Retrospective Derivation and Validation Cohort Study

BACKGROUND

Postoperative nausea and vomiting (PONV) is a key driver of unplanned admission and patient satisfaction after surgery. Because traditional risk factors do not completely explain variability in risk, this study hypothesized that genetics may contribute to the overall risk for this complication. The objective of this research is to perform a genome-wide association study of PONV, derive a polygenic risk score for PONV, assess associations between the risk score and PONV in a validation cohort, and compare any genetic contributions to known clinical risks for PONV.

METHODS

Surgeries with integrated genetic and perioperative data performed under general anesthesia at Michigan Medicine (Ann Arbor, Michigan) and Vanderbilt University Medical Center (Nashville, Tennessee) were studied. PONV was defined as nausea or emesis occurring and documented in the postanesthesia care unit. In the discovery phase, genome-wide association studies were performed on each genetic cohort, and the results were meta-analyzed. Next, the polygenic phase assessed whether a polygenic score, derived from genome-wide association study in a derivation cohort from Vanderbilt University Medical Center, improved prediction within a validation cohort from Michigan Medicine, as quantified by discrimination (c-statistic) and net reclassification index.

RESULTS

Of 64,523 total patients, 5,703 developed PONV (8.8%). The study identified 46 genetic variants exceeding the threshold of P < 1 × 10-5, occurring with minor allele frequency greater than 1%, and demonstrating concordant effects in both cohorts. Standardized polygenic score was associated with PONV in a basic model, controlling for age and sex (adjusted odds ratio, 1.027 per SD increase in overall genetic risk; 95% CI, 1.001 to 1.053; P = 0.044), a model based on known clinical risks (adjusted odds ratio, 1.029; 95% CI, 1.003 to 1.055; P = 0.030), and a full clinical regression, controlling for 21 demographic, surgical, and anesthetic factors, (adjusted odds ratio, 1.029; 95% CI, 1.002 to 1.056; P = 0.033). The addition of polygenic score improved overall discrimination in models based on known clinical risk factors (c-statistic, 0.616 compared to 0.613; P = 0.028) and improved net reclassification of 4.6% of cases.

CONCLUSIONS

Standardized polygenic risk was associated with PONV in all three of the study's models, but the genetic influence was smaller than exerted by clinical risk factors. Specifically, a patient with a polygenic risk score greater than 1 SD above the mean has 2 to 3% greater odds of developing PONV when compared to the baseline population, which is at least an order of magnitude smaller than the increase associated with having prior PONV or motion sickness (55%), having a history of migraines (17%), or being female (83%) and is not clinically significant. Furthermore, the use of a polygenic risk score does not meaningfully improve discrimination compared to clinical risk factors and is not clinically useful.

EDITOR’S PERSPECTIVE

Genetic studies through the lens of gene networks

Understanding the genetic basis of complex traits is a longstanding challenge in the field of genomics. Genome-wide association studies (GWAS) have identified thousands of variant-trait associations, but most of these variants are located in non-coding regions, making the link to biological function elusive. While traditional approaches, such as transcriptome-wide association studies (TWAS), have advanced our understanding by linking genetic variants to gene expression, they often overlook gene-gene interactions. Here, we review current approaches to integrate different molecular data, leveraging machine learning methods to identify gene modules based on co-expression and functional relationships. These integrative approaches, like PhenoPLIER, combine TWAS and drug-induced transcriptional profiles to effectively capture biologically meaningful gene networks. This integration provides a context-specific understanding of disease processes while highlighting both core and peripheral genes. These insights pave the way for novel therapeutic targets and enhance the interpretability of genetic studies in personalized medicine.

Machine Learning Reveals the Diversity of Human 3D Chromatin Contact Patterns

Understanding variation in chromatin contact patterns across diverse humans is critical for interpreting noncoding variants and their effects on gene expression and phenotypes. However, experimental determination of chromatin contact patterns across large samples is prohibitively expensive. To overcome this challenge, we develop and validate a machine learning method to quantify the variation in 3D chromatin contacts at 2 kilobase resolution from genome sequence alone. We apply this approach to thousands of human genomes from the 1000 Genomes Project and the inferred hominin ancestral genome. While patterns of 3D contact divergence genome wide are qualitatively similar to patterns of sequence divergence, we find substantial differences in 3D divergence and sequence divergence in local 1 megabase genomic windows. In particular, we identify 392 windows with significantly greater 3D divergence than expected from sequence. Moreover, for 31% of genomic windows, a single individual has a rare divergent 3D contact map pattern. Using in silico mutagenesis, we find that most single nucleotide sequence changes do not result in changes to 3D chromatin contacts. However, in windows with substantial 3D divergence just one or a few variants can lead to divergent 3D chromatin contacts without the individuals carrying those variants having high sequence divergence. In summary, inferring 3D chromatin contact maps across human populations reveals variable contact patterns. We anticipate that these genetically diverse maps of 3D chromatin contact will provide a reference for future work on the function and evolution of 3D chromatin contact variation across human populations.

Translation of genome-wide association study: from genomic signals to biological insights

Since the turn of the 21st century, genome-wide association study (GWAS) have successfully identified genetic signals associated with a myriad of common complex traits and diseases. As we transition from establishing robust genetic associations with diverse phenotypes, the central challenge is now focused on characterizing the underlying functional mechanisms driving these signals. Previous GWAS efforts have revealed multiple variants, each conferring relatively subtle susceptibility, collectively contributing to the pathogenesis of various common diseases. Such variants can further exhibit associations with multiple other traits and differ across ancestries, plus disentangling causal variants from non-causal due to linkage disequilibrium complexities can lead to challenges in drawing direct biological conclusions. Combined with cellular context considerations, such challenges can reduce the capacity to definitively elucidate the biological significance of GWAS signals, limiting the potential to define mechanistic insights. This review will detail current and anticipated approaches for functional interpretation of GWAS signals, both in terms of characterizing the underlying causal variants and the corresponding effector genes.

Integrative computational analyses implicate regulatory genomic elements contributing to spina bifida

Purpose

Spina bifida (SB) arises from complex genetic interactions that converge to interfere with neural tube closure. Understanding the precise patterns conferring SB risk requires a deep exploration of the genomic networks and molecular pathways that govern neurulation. This study aims to delineate genome-wide regulatory signatures underlying SB pathophysiology.

Methods

An untargeted, genome-wide approach was used to interrogate regulatory regions for rare single-nucleotide and copy-number variants (rSNVs and rCNVs, respectively) predicted to affect gene expression, comparing results from SB patients with healthy controls. Qualifying variants were subjected to a deep learning prioritization framework to identify the most functionally relevant variants, as well as the likely target genes affected by these rare regulatory variants.

Results

This ensemble of computational tools identified rSNVs in specific transcription factor binding sites (TFBSs) that distinguish SB cases from controls. rSNV enrichment was found in specific TFBSs, especially CCCTC-binding factor binding sites. These TFBSs were subjected to a deep learning prioritization framework to identify the most functionally relevant variants, as well as the likely target genes affected by these rSNVs. The functional pathways or modules implicated by these regulated genes serve protein transport, cilia assembly, and central nervous system development. Moreover, the detected rare copy-number variants in SB cases are positioned to disrupt gene regulatory networks and alter 3-dimensional genomic architectures, including brain-specific enhancers and topologically associated domain boundaries of relevant cell types.

Conclusion

Our study provides a resource for identifying and interpreting genomic regulatory DNA variant contributions to human SB genetic predisposition.

Implementing Whole Genome Sequencing (WGS) in Clinical Practice: Advantages, Challenges, and Future Perspectives

The integration of whole genome sequencing (WGS) into all aspects of modern medicine represents the next step in the evolution of healthcare. Using this technology, scientists and physicians can observe the entire human genome comprehensively, generating a plethora of new sequencing data. Modern computational analysis entails advanced algorithms for variant detection, as well as complex models for classification. Data science and machine learning play a crucial role in the processing and interpretation of results, using enormous databases and statistics to discover new and support current genotype-phenotype correlations. In clinical practice, this technology has greatly enabled the development of personalized medicine, approaching each patient individually and in accordance with their genetic and biochemical profile. The most propulsive areas include rare disease genomics, oncogenomics, pharmacogenomics, neonatal screening, and infectious disease genomics. Another crucial application of WGS lies in the field of multi-omics, working towards the complete integration of human biomolecular data. Further technological development of sequencing technologies has led to the birth of third and fourth-generation sequencing, which include long-read sequencing, single-cell genomics, and nanopore sequencing. These technologies, alongside their continued implementation into medical research and practice, show great promise for the future of the field of medicine.

Comparative analysis of models in predicting the effects of SNPs on TF-DNA binding using large-scale in vitro and in vivo data

Non-coding variants associated with complex traits can alter the motifs of transcription factor (TF)-deoxyribonucleic acid binding. Although many computational models have been developed to predict the effects of non-coding variants on TF binding, their predictive power lacks systematic evaluation. Here we have evaluated 14 different models built on position weight matrices (PWMs), support vector machines, ordinary least squares and deep neural networks (DNNs), using large-scale in vitro (i.e. SNP-SELEX) and in vivo (i.e. allele-specific binding, ASB) TF binding data. Our results show that the accuracy of each model in predicting SNP effects in vitro significantly exceeds that achieved in vivo. For in vitro variant impact prediction, kmer/gkm-based machine learning methods (deltaSVM_HT-SELEX, QBiC-Pred) trained on in vitro datasets exhibit the best performance. For in vivo ASB variant prediction, DNN-based multitask models (DeepSEA, Sei, Enformer) trained on the ChIP-seq dataset exhibit relatively superior performance. Among the PWM-based methods, tRap demonstrates better performance in both in vitro and in vivo evaluations. In addition, we find that TF classes such as basic leucine zipper factors could be predicted more accurately, whereas those such as C2H2 zinc finger factors are predicted less accurately, aligning with the evolutionary conservation of these TF classes. We also underscore the significance of non-sequence factors such as cis-regulatory element type, TF expression, interactions and post-translational modifications in influencing the in vivo predictive performance of TFs. Our research provides valuable insights into selecting prioritization methods for non-coding variants and further optimizing such models.

cisDynet: An integrated platform for modeling gene-regulatory dynamics and networks

Chromatin accessibility sequencing has been widely used for uncovering genetic regulatory mechanisms and inferring gene regulatory networks. However, effectively integrating large-scale chromatin accessibility datasets has posed a significant challenge. This is due to the lack of a comprehensive end-to-end solution, as many existing tools primarily emphasize data preprocessing and overlook downstream analyses. To bridge this gap, we have introduced cisDynet, a holistic solution that combines streamlined data preprocessing using Snakemake and R functions with advanced downstream analysis capabilities. cisDynet excels in conventional data analyses, encompassing peak statistics, peak annotation, differential analysis, motif enrichment analysis, and more. Additionally, it allows to perform sophisticated data exploration, such as tissue-specific peak identification, time course data modeling, integration of RNA-seq data to establish peak-to-gene associations, constructing regulatory networks, and conducting enrichment analysis of genome-wide association study (GWAS) variants. As a proof of concept, we applied cisDynet to reanalyze comprehensive ATAC-seq datasets across various tissues from the Encyclopedia of DNA Elements (ENCODE) project. The analysis successfully delineated tissue-specific open chromatin regions (OCRs), established connections between OCRs and target genes, and effectively linked these discoveries with 1861 GWAS variants. Furthermore, cisDynet was instrumental in dissecting the time course open chromatin data of mouse embryonic development, revealing the dynamic behavior of OCRs over developmental stages and identifying key transcription factors governing differentiation trajectories. In summary, cisDynet offers researchers a user-friendly solution that minimizes the need for extensive coding, ensures the reproducibility of results, and greatly simplifies the exploration of epigenomic data.

Genome-wide multi-trait analysis of irritable bowel syndrome and related mental conditions identifies 38 new independent variants

BACKGROUND

Irritable bowel syndrome (IBS) is a chronic disorder of gut-brain interaction frequently accompanied by mental conditions, including depression and anxiety. Despite showing substantial heritability and being partly determined by a genetic component, the genetic underpinnings explaining the high rates of comorbidity remain largely unclear and there are no conclusive data on the temporal relationship between them. Exploring the overlapping genetic architecture between IBS and mental conditions may help to identify novel genetic loci and biological mechanisms underlying IBS and causal relationships between them.

METHODS

We quantified the genetic overlap between IBS, neuroticism, depression and anxiety, conducted a multi-trait genome-wide association study (GWAS) considering these traits and investigated causal relationships between them by using the largest GWAS to date.

RESULTS

IBS showed to be a highly polygenic disorder with extensive genetic sharing with mental conditions. Multi-trait analysis of IBS and neuroticism, depression and anxiety identified 42 genome-wide significant variants for IBS, of which 38 are novel. Fine-mapping risk loci highlighted 289 genes enriched in genes upregulated during early embryonic brain development and gene-sets related with psychiatric, digestive and autoimmune disorders. IBS-associated genes were enriched for target genes of anti-inflammatory and antirheumatic drugs, anesthetics and opioid dependence pharmacological treatment. Mendelian-randomization analysis accounting for correlated pleiotropy identified bidirectional causal effects between IBS and neuroticism and depression and causal effects of the genetic liability of IBS on anxiety.

CONCLUSIONS

These findings provide evidence of the polygenic architecture of IBS, identify novel genome-wide significant variants for IBS and extend previous knowledge on the genetic overlap and relationship between gastrointestinal and mental disorders.

Multiplexed functional genomic assays to decipher the noncoding genome

Linkage disequilibrium and the incomplete regulatory annotation of the noncoding genome complicates the identification of functional noncoding genetic variants and their causal association with disease. Current computational methods for variant prioritization have limited predictive value, necessitating the application of highly parallelized experimental assays to efficiently identify functional noncoding variation. Here, we summarize two distinct approaches, massively parallel reporter assays and CRISPR-based pooled screens and describe their flexible implementation to characterize human noncoding genetic variation at unprecedented scale. Each approach provides unique advantages and limitations, highlighting the importance of multimodal methodological integration. These multiplexed assays of variant effects are undoubtedly poised to play a key role in the experimental characterization of noncoding genetic risk, informing our understanding of the underlying mechanisms of disease-associated loci and the development of more robust predictive classification algorithms.