Finemap-MiXeR: A variational Bayesian approach for genetic finemapping.

Genome-wide association studies (GWAS) implicate broad genomic loci containing clusters of highly correlated genetic variants. Finemapping techniques can select and prioritize variants within each GWAS locus which are more likely to have a functional influence on the trait. Here, we present a novel method, Finemap-MiXeR, for finemapping causal variants from GWAS summary statistics, controlling for correlation among variants due to linkage disequilibrium. Our method is based on a variational Bayesian approach and direct optimization of the Evidence Lower Bound (ELBO) of the likelihood function derived from the MiXeR model. After obtaining the analytical expression for ELBO's gradient, we apply Adaptive Moment Estimation (ADAM) algorithm for optimization, allowing us to obtain the posterior causal probability of each variant. Using these posterior causal probabilities, we validated Finemap-MiXeR across a wide range of scenarios using both synthetic data, and real data on height from the UK Biobank. Comparison of Finemap-MiXeR with two existing methods, FINEMAP and SuSiE RSS, demonstrated similar or improved accuracy. Furthermore, our method is computationally efficient in several aspects. For example, unlike many other methods in the literature, its computational complexity does not increase with the number of true causal variants in a locus and it does not require any matrix inversion operation. The mathematical framework of Finemap-MiXeR is flexible and may also be applied to other problems including cross-trait and cross-ancestry finemapping.

COSGAP: COntainerized Statistical Genetics Analysis Pipelines.

Summary: The collection and analysis of sensitive data in large-scale consortia for statistical genetics is hampered by multiple challenges, due to their non-shareable nature. Time-consuming issues in installing software frequently arise due to different operating systems, software dependencies, and limited internet access. For federated analysis across sites, it can be challenging to resolve different problems, including format requirements, data wrangling, setting up analysis on high-performance computing (HPC) facilities, etc. Easier, more standardized, automated protocols and pipelines can be solutions to overcome these issues. We have developed one such solution for statistical genetic data analysis using software container technologies. This solution, named COSGAP: "COntainerized Statistical Genetics Analysis Pipelines," consists of already established software tools placed into Singularity containers, alongside corresponding code and instructions on how to perform statistical genetic analyses, such as genome-wide association studies, polygenic scoring, LD score regression, Gaussian Mixture Models, and gene-set analysis. Using provided helper scripts written in Python, users can obtain auto-generated scripts to conduct the desired analysis either on HPC facilities or on a personal computer. COSGAP is actively being applied by users from different countries and projects to conduct genetic data analyses without spending much effort on software installation, converting data formats, and other technical requirements. Availability and implementation: COSGAP is freely available on GitHub (https://github.com/comorment/containers) under the GPLv3 license.

A Breast Cancer Polygenic Risk Score Is Feasible for Risk Stratification in the Norwegian Population.

BACKGROUND: Statistical associations of numerous single nucleotide polymorphisms with breast cancer (BC) have been identified in genome-wide association studies (GWAS). Recent evidence suggests that a Polygenic Risk Score (PRS) can be a useful risk stratification instrument for a BC screening strategy, and a PRS test has been developed for clinical use. The performance of the PRS is yet unknown in the Norwegian population. AIM: To evaluate the performance of PRS models for BC in a Norwegian dataset. METHODS: We investigated a sample of 1053 BC cases and 7094 controls from different regions of Norway. PRS values were calculated using four PRS models, and their performance was evaluated by the area under the curve (AUC) and the odds ratio (OR). The effect of the PRS on the age of onset of BC was determined by a Cox regression model, and the lifetime absolute risk of developing BC was calculated using the iCare tool. RESULTS: The best performing PRS model included 3820 SNPs, which yielded an AUC = 0.625 and an OR = 1.567 per one standard deviation increase. The PRS values of the samples correlate with an increased risk of BC, with a hazard ratio of 1.494 per one standard deviation increase (95% confidence interval of 1.406-1.588). The individuals in the highest decile of the PRS have at least twice the risk of developing BC compared to the individuals with a median PRS. The results in this study with Norwegian samples are coherent with the findings in the study conducted using Estonian and UK Biobank samples. CONCLUSION: The previously validated PRS models have a similar observed accuracy in the Norwegian data as in the UK and Estonian populations. A PRS provides a meaningful association with the age of onset of BC and lifetime risk. Therefore, as suggested in Estonia, a PRS may also be integrated into the screening strategy for BC in Norway.

Molecular Communication for Equilibrium State Estimation in Biochemical Processes on a Lab-on-a-Chip.

A basic problem in molecular biology is to estimate equilibrium states of biochemical processes. To this end, advanced spectroscopy methods have been developed in order to estimate chemical concentrations in situ or in vivo. However, such spectroscopy methods can require special conditions that do not allow direct observation of the biochemical process. A natural means of resolving this problem is to transmit chemical signals to another location within a lab-on-a-chip device; that is, employing molecular communication in order to perform spectroscopy in a different location. In this paper, we develop such a signaling strategy and estimation algorithms for equilibrium states of a biochemical process. In two biologically-inspired models, we then study via simulation the tradeoff between the rate of obtaining spectroscopy measurements and the estimation error, providing insights into requirements of spectroscopy devices for high-throughput biological assays.

A General Analytical Approximation to Impulse Response of 3-D Microfluidic Channels in Molecular Communication.

In this paper, the impulse response for a 3-D microfluidic channel in the presence of Poiseuille flow is obtained by solving the diffusion equation in radial coordinates. Using the radial distribution, the axial distribution is then approximated accordingly. Since Poiseuille flow velocity changes with radial position, molecules have different axial properties for different radial distributions. We, therefore, present a piecewise function for the axial distribution of the molecules in the channel considering this radial distribution. We lay evidence for our theoretical derivations for impulse response of the microfluidic channel and radial distribution of molecules through comparing them using various Monte Carlo simulations. Finally, the communication performance of the channel is examined.