A Weighted SNP Correlation Network Method for Estimating Polygenic Risk Scores

Association of whole-person eigen-polygenic risk scores with Alzheimer's disease

Late-Onset Alzheimer's Disease (LOAD) is a heterogeneous neurodegenerative disorder with complex etiology and high heritability. Its multifactorial risk profile and large portions of unexplained heritability suggest the involvement of yet unidentified genetic risk factors. Here we describe the "whole person" genetic risk landscape of polygenic risk scores for 2218 traits in 2044 elderly individuals and test if novel eigen-PRSs derived from clustered subnetworks of single-trait PRSs can improve the prediction of LOAD diagnosis, rates of cognitive decline, and canonical LOAD neuropathology. Network analyses revealed distinct clusters of PRSs with clinical and biological interpretability. Novel eigen-PRSs (ePRS) from these clusters significantly improved LOAD-related phenotypes prediction over current state-of-the-art LOAD PRS models. Notably, an ePRS representing clusters of traits related to cholesterol levels was able to improve variance explained in a model of the brain-wide beta-amyloid burden by 1.7% (likelihood ratio test P = 9.02 × 10-7). All associations of ePRS with LOAD phenotypes were eliminated by the removal of APOE-proximal loci. However, our association analysis identified modules characterized by PRSs of high cholesterol and LOAD. We believe this is due to the influence of the APOE region from both PRSs. We found significantly higher mean SNP effects for LOAD in the intersecting APOE region SNPs. Combining genetic risk factors for vascular traits and dementia could improve current single-trait PRS models of LOAD, enhancing the use of PRS in risk stratification. Our results are catalogued for the scientific community, to aid in generating new hypotheses based on our maps of clustered PRSs and associations with LOAD-related phenotypes.

Unraveling the Variability of Human Satiation: Implications for Precision Obesity Management

Satiation is the physiologic process that regulates meal size and termination, and it is quantified by the calories consumed to reach satiation. Given its role in energy intake, changes in satiation contribute to obesity's pathogenesis. Our study employed a protocolized approach to study the components of food intake regulation including a standardized breakfast, a gastric emptying study, appetite sensation testing, and a satiation measurement by an ad libitummeal test. These studies revealed that satiation is highly variable among individuals, and while baseline characteristics, anthropometrics, body composition and hormones, contribute to this variability, these factors do not fully account for it. To address this gap, we explored the role of a germline polygenic risk score, which demonstrated a robust association with satiation. Furthermore, we developed a machine-learning-assisted gene risk score to predict satiation and leveraged this prediction to anticipate responses to anti-obesity medications. Our findings underscore the significance of satiation, its inherent variability, and the potential of a genetic risk score to forecast it, ultimately allowing us to predict responses to different anti-obesity interventions.

Cancer Risk Score Prediction Based on a Single-Nucleotide Polymorphism Network

Objectives

Genome-wide association studies (GWAS) are performed to study the associations between genetic variants with respect to certain phenotypic traits such as cancer. However, the method that is commonly used in GWAS assumes that certain traits are solely affected by a single mutation. We propose a network analysis method, in which we generate association networks of single-nucleotide polymorphisms (SNPs) that can differentiate case and control groups. We hypothesize that certain phenotypic traits are attributable to mutations in groups of associated SNPs.

Methods

We propose a method based on a network analysis framework to study SNP-SNP interactions related to cancer incidence. We employed logistic regression to measure the significance of all SNP pairs from GWAS for the incidence of colorectal cancer and computed a cancer risk score based on the generated SNP networks.

Results

We demonstrated our method in a dataset from a case-control study of colorectal cancer in the South Sulawesi population. From the GWAS results, 20,094 pairs of 200 SNPs were created. We obtained one cluster containing four pairs of five SNPs that passed the filtering threshold based on their p-values. A locus on chromosome 12 (12:54410007) was found to be strongly connected to the four variants on chromosome 1. A polygenic risk score was computed from the five SNPs, and a significant difference in colorectal cancer risk was obtained between the case and control groups.

Conclusions

Our results demonstrate the applicability of our method to understand SNP-SNP interactions and compute risk scores for various types of cancer.

Integration of Imaging Genomics Data for the Study of Alzheimer's Disease Using Joint-Connectivity-Based Sparse Nonnegative Matrix Factorization

Imaging genetics reveals the connection between microscopic genetics and macroscopic imaging, enabling the identification of disease biomarkers. In this work, we make full use of prior knowledge that has significant reference value for investigating the correlation between the brain and genetics to explore more biologically substantial biomarkers. In this paper, we propose joint-connectivity-based sparse nonnegative matrix factorization (JCB-SNMF). The algorithm simultaneously projects structural magnetic resonance imaging (sMRI), single-nucleotide polymorphism sites (SNPs), and gene expression data onto a common feature space, where heterogeneous variables with large coefficients in the same projection direction form a common module. In addition, the connectivity information for each region of the brain and genetic data are added as prior knowledge to identify regions of interest (ROIs), SNPs, and gene-related risks related to Alzheimer's disease (AD) patients. GraphNet regularization increases the anti-noise performance of the algorithm and the biological interpretability of the results. The simulation results show that compared with other NMF-based algorithms (JNMF, JSNMNMF), JCB-SNMF has better anti-noise performance and can identify and predict biomarkers closely related to AD from significant modules. By constructing a protein-protein interaction (PPI) network, we identified SF3B1, RPS20, and RBM14 as potential biomarkers of AD. We also found some significant SNP-ROI and gene-ROI pairs. Among them, two SNPs rs4472239 and rs11918049 and three genes KLHL8, ZC3H11A, and OSGEPL1 may have effects on the gray matter volume of multiple brain regions. This model provides a new way to further integrate multimodal impact genetic data to identify complex disease association patterns.

Machine learning based disease prediction from genotype data

Using results from genome-wide association studies for understanding complex traits is a current challenge. Here we review how genotype data can be used with different machine learning (ML) methods to predict phenotype occurrence and severity from genotype data. We discuss common feature encoding schemes and how studies handle the often small number of samples compared to the huge number of variants. We compare which ML methods are being applied, including recent results using deep neural networks. Further, we review the application of methods for feature explanation and interpretation.

Dementia key gene identification with multi-layered SNP-gene-disease network

MOTIVATION

Recently, various approaches for diagnosing and treating dementia have received significant attention, especially in identifying key genes that are crucial for dementia. If the mutations of such key genes could be tracked, it would be possible to predict the time of onset of dementia and significantly aid in developing drugs to treat dementia. However, gene finding involves tremendous cost, time and effort. To alleviate these problems, research on utilizing computational biology to decrease the search space of candidate genes is actively conducted.

In this study, we propose a framework in which diseases, genes and single-nucleotide polymorphisms are represented by a layered network, and key genes are predicted by a machine learning algorithm. The algorithm utilizes a network-based semi-supervised learning model that can be applied to layered data structures.

RESULTS

The proposed method was applied to a dataset extracted from public databases related to diseases and genes with data collected from 186 patients. A portion of key genes obtained using the proposed method was verified in silico through PubMed literature, and the remaining genes were left as possible candidate genes.

AVAILABILITY AND IMPLEMENTATION

The code for the framework will be available at http://www.alphaminers.net/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Exploring the complexity of soybean (Glycine max) transcriptional regulation using global gene co-expression networks

MAIN CONCLUSION

We report a soybean gene co-expression network built with data from 1284 RNA-Seq experiments, which was used to identify important regulators, modules and to elucidate the fates of gene duplicates. Soybean (Glycine max (L.) Merr.) is one of the most important crops worldwide, constituting a major source of protein and edible oil. Gene co-expression networks (GCN) have been extensively used to study transcriptional regulation and evolution of genes and genomes. Here, we report a soybean GCN using 1284 publicly available RNA-Seq samples from 15 distinct tissues. We found modules that are differentially regulated in specific tissues, comprising processes such as photosynthesis, gluconeogenesis, lignin metabolism, and response to biotic stress. We identified transcription factors among intramodular hubs, which probably integrate different pathways and shape the transcriptional landscape in different conditions. The top hubs for each module tend to encode proteins with critical roles, such as succinate dehydrogenase and RNA polymerase subunits. Importantly, gene essentiality was strongly correlated with degree centrality and essential hubs were enriched in genes involved in nucleic acids metabolism and regulation of cell replication. Using a guilt-by-association approach, we predicted functions for 93 of 106 hubs without functional description in soybean. Most of the duplicated genes had different transcriptional profiles, supporting their functional divergence, although paralogs originating from whole-genome duplications (WGD) are more often preserved in the same module than those from other mechanisms. Together, our results highlight the importance of GCN analysis in unraveling key functional aspects of the soybean genome, in particular those associated with hub genes and WGD events.

Joint-Connectivity-Based Sparse Canonical Correlation Analysis of Imaging Genetics for Detecting Biomarkers of Parkinson's Disease

Imaging genetics is a method used to detect associations between imaging and genetic variables. Some researchers have used sparse canonical correlation analysis (SCCA) for imaging genetics. This study was conducted to improve the efficiency and interpretability of SCCA. We propose a connectivity-based penalty for incorporating biological prior information. Our proposed approach, named joint connectivity-based SCCA (JCB-SCCA), includes the proposed penalty and can handle multi-modal neuroimaging datasets. Different neuroimaging techniques provide distinct information on the brain and have been used to investigate various neurological disorders, including Parkinson's disease (PD). We applied our algorithm to simulated and real imaging genetics datasets for performance evaluation. Our algorithm was able to select important features in a more robust manner compared with other multivariate methods. The algorithm revealed promising features of single-nucleotide polymorphisms and brain regions related to PD by using a real imaging genetic dataset. The proposed imaging genetics model can be used to improve clinical diagnosis in the form of novel potential biomarkers. We hope to apply our algorithm to cohorts such as Alzheimer's patients or healthy subjects to determine the generalizability of our algorithm.

Data Integration in Poplar: 'Omics Layers and Integration Strategies

Populus trichocarpa is an important biofuel feedstock that has been the target of extensive research and is emerging as a model organism for plants, especially woody perennials. This research has generated several large ‘omics datasets. However, only few studies in Populus have attempted to integrate various data types. This review will summarize various ‘omics data layers, focusing on their application in Populus species. Subsequently, network and signal processing techniques for the integration and analysis of these data types will be discussed, with particular reference to examples in Populus.

Multi-Phenotype Association Decomposition: Unraveling Complex Gene-Phenotype Relationships

Various patterns of multi-phenotype associations (MPAs) exist in the results of Genome Wide Association Studies (GWAS) involving different topologies of single nucleotide polymorphism (SNP)-phenotype associations. These can provide interesting information about the different impacts of a gene on closely related phenotypes or disparate phenotypes (pleiotropy). In this work we present MPA Decomposition, a new network-based approach which decomposes the results of a multi-phenotype GWAS study into three bipartite networks, which, when used together, unravel the multi-phenotype signatures of genes on a genome-wide scale. The decomposition involves the construction of a phenotype powerset space, and subsequent mapping of genes into this new space. Clustering of genes in this powerset space groups genes based on their detailed MPA signatures. We show that this method allows us to find multiple different MPA and pleiotropic signatures within individual genes and to classify and cluster genes based on these SNP-phenotype association topologies. We demonstrate the use of this approach on a GWAS analysis of a large population of 882 Populus trichocarpa genotypes using untargeted metabolomics phenotypes. This method should prove invaluable in the interpretation of large GWAS datasets and aid in future synthetic biology efforts designed to optimize phenotypes of interest.

On the Influence of Structural Connectivity on the Correlation Patterns and Network Synchronization

Since brain structural connectivity is the foundation of its functionality, in order to understand brain abilities, studying the relation between structural and functional connectivity is essential. Several approaches have been applied to measure the role of the structural connectivity in the emergent correlation/synchronization patterns. In this study, we investigates the cross-correlation and synchronization sensitivity to coupling strength between neural regions for different topological networks. We model the neural populations by a neural mass model that express an oscillatory dynamic. The results highlight that coupling between neural ensembles leads to various cross-correlation patterns and local synchrony even on an ordered network. Moreover, as the network departs from an ordered organization to a small-world architecture, correlation patterns, and synchronization dynamics change. Interestingly, at a certain range of the synaptic strength, by fixing the structural conditions, different organized patterns are seen at the different input signals. This variety switches to a bifurcation region by increasing the synaptic strength. We show that topological variations is a major factor of synchronization behavior and lead to alterations in correlated local clusters. We found the coupling strength (between cortical areas) to be especially important at conversions of correlation and synchronization states. Since correlation patterns generate functional connections and transitions of functional connectivity have been related to cognitive operations, these diverse correlation patterns may be considered as different dynamical states corresponding to various cognitive tasks.