Genome analysis and pleiotropy assessment using causal networks with loss of function mutation and metabolomics

Broadcasters, receivers, functional groups of metabolites, and the link to heart failure by revealing metabolomic network connectivity

BACKGROUND AND OBJECTIVE

Blood-based small molecule metabolites offer easy accessibility and hold significant potential for insights into health processes, the impact of lifestyle, and genetic variation on disease, enabling precise risk prevention. In a prospective study with records of heart failure (HF) incidence, we present metabolite profiling data from individuals without HF at baseline.

METHODS

We uncovered the interconnectivity of metabolites using data-driven and causal networks augmented with polygenic factors. Exploring the networks, we identified metabolite broadcasters, receivers, mediators, and subnetworks corresponding to functional classes of metabolites, and provided insights into the link between metabolomic architecture and regulation in health. We incorporated the network structure into the identification of metabolites associated with HF to control the effect of confounding metabolites.

RESULTS

We identified metabolites associated with higher and lower risk of HF incidence, such as glycine, ureidopropionic and glycocholic acids, and LPC 18:2. These associations were not confounded by the other metabolites due to uncovering the connectivity among metabolites and adjusting each association for the confounding metabolites. Examples of our findings include the direct influence of asparagine on glycine, both of which were inversely associated with HF. These two metabolites were influenced by polygenic factors and only essential amino acids, which are not synthesized in the human body and are obtained directly from the diet.

CONCLUSION

Metabolites may play a critical role in linking genetic background and lifestyle factors to HF incidence. Revealing the underlying connectivity of metabolites associated with HF strengthens the findings and facilitates studying complex conditions like HF.

WITHDRAWN: Broadcasters, receivers, functional groups of metabolites and the link to heart failure using polygenic factors

The full text of this preprint has been withdrawn, as it was submitted in error. Therefore, the authors do not wish this work to be cited as a reference. Questions should be directed to the corresponding author.

WITHDRAWN: Broadcasters, receivers, functional groups of metabolites and the link to heart failure using polygenic factors

The full text of this preprint has been withdrawn, as it was submitted in error. Therefore, the authors do not wish this work to be cited as a reference. Questions should be directed to the corresponding author.

Inference of causal metabolite networks in the presence of invalid instrumental variables with GWAS summary data

We propose structural equation models (SEMs) as a general framework to infer causal networks for metabolites and other complex traits. Traditionally SEMs are used only for individual-level data under the assumption that all instrumental variables (IVs) are valid. To overcome these limitations, we propose both one- and two-sample approaches for causal network inference based on SEMs that can: (1) perform causal analysis and discover causal relationships among multiple traits; (2) account for the possible presence of some invalid IVs; (3) allow for data analysis using only genome-wide association studies (GWAS) summary statistics when individual-level data are not available; (4) consider the possibility of bidirectional relationships between traits. Our method employs a simple stepwise selection to identify invalid IVs, thus avoiding false positives while possibly increasing true discoveries based on two-stage least squares (2SLS). We use both real GWAS data and simulated data to demonstrate the superior performance of our method over the standard 2SLS/SEMs. For real data analysis, our proposed approach is applied to a human blood metabolite GWAS summary data set to uncover putative causal relationships among the metabolites; we also identify some metabolites (putative) causal to Alzheimer's disease (AD), which, along with the inferred causal metabolite network, suggest some possible pathways of metabolites involved in AD.

Broadcasters, receivers, functional groups of metabolites and the link to heart failure progression using polygenic factors

In a prospective study with records of heart failure (HF) incidence, we present metabolite profiling data from individuals without HF at baseline. We uncovered the interconnectivity of metabolites using data-driven and causal networks augmented with polygenic factors. Exploring the networks, we identified metabolite broadcasters, receivers, mediators, and subnetworks corresponding to functional classes of metabolites, and provided insights into the link between metabolomic architecture and regulation in health. We incorporated the network structure into the identification of metabolites associated with HF to control the effect of confounding metabolites. We identified metabolites associated with higher or lower risk of HF incidence, the associations that were not confounded by the other metabolites, such as glycine, ureidopropionic and glycocholic acids, and LPC 18:2. We revealed the underlying relationships of the findings. For example, asparagine directly influenced glycine, and both were inversely associated with HF. These two metabolites were influenced by polygenic factors and only essential amino acids which are not synthesized in the human body and come directly from the diet. Metabolites may play a critical role in linking genetic background and lifestyle factors to HF progression. Revealing the underlying connectivity of metabolites associated with HF strengthens the findings and facilitates a mechanistic understanding of HF progression.

Genome-wide analysis of oxylipins and oxylipin profiles in a pediatric population

Background

Oxylipins are inflammatory biomarkers derived from omega-3 and-6 fatty acids implicated in inflammatory diseases but have not been studied in a genome-wide association study (GWAS). The aim of this study was to identify genetic loci associated with oxylipins and oxylipin profiles to identify biologic pathways and therapeutic targets for oxylipins.

Methods

We conducted a GWAS of plasma oxylipins in 316 participants in the Diabetes Autoimmunity Study in the Young (DAISY). DNA samples were genotyped using the TEDDY-T1D Exome array, and additional variants were imputed using the Trans-Omics for Precision Medicine (TOPMed) multi-ancestry reference panel. Principal components analysis of 36 plasma oxylipins was used to capture oxylipin profiles. PC1 represented linoleic acid (LA)- and alpha-linolenic acid (ALA)-related oxylipins, and PC2 represented arachidonic acid (ARA)-related oxylipins. Oxylipin PC1, PC2, and the top five loading oxylipins from each PC were used as outcomes in the GWAS (genome-wide significance: p < 5×10-8).

Results

The SNP rs143070873 was associated with (p < 5×10-8) the LA-related oxylipin 9-HODE, and rs6444933 (downstream of CLDN11) was associated with the LA-related oxylipin 13 S-HODE. A locus between MIR1302-7 and LOC100131146, rs10118380 and an intronic variant in TRPM3 were associated with the ARA-related oxylipin 11-HETE. These loci are involved in inflammatory signaling cascades and interact with PLA2, an initial step to oxylipin biosynthesis.

Conclusion

Genetic loci involved in inflammation and oxylipin metabolism are associated with oxylipin levels.

Genetic influences on human blood metabolites in the Japanese population

An increase in ethnic diversity in genetic studies has the potential to provide unprecedented insights into how genetic variations influence human phenotypes. In this study, we conducted a quantitative trait locus (QTL) analysis of 121 metabolites measured using gas chromatography-mass spectrometry with plasma samples from 4,888 Japanese individuals. We found 60 metabolite-gene associations, of which 13 have not been previously reported. Meta-analyses with another Japanese and a European study identified six and two additional unreported loci, respectively. Genetic variants influencing metabolite levels were more enriched in protein-coding regions than in the regulatory regions while being associated with the risk of various diseases. Finally, we identified a signature of strong negative selection for uric acid ( S ˆ  = -1.53, p = 6.2 × 10^-18). Our study expanded the knowledge of genetic influences on human blood metabolites, providing valuable insights into their physiological, pathological, and selective properties.

Deep neural networks with knockoff features identify nonlinear causal relations and estimate effect sizes in complex biological systems

BACKGROUND

Learning the causal structure helps identify risk factors, disease mechanisms, and candidate therapeutics for complex diseases. However, although complex biological systems are characterized by nonlinear associations, existing bioinformatic methods of causal inference cannot identify the nonlinear relationships and estimate their effect size.

RESULTS

To overcome these limitations, we developed the first computational method that explicitly learns nonlinear causal relations and estimates the effect size using a deep neural network approach coupled with the knockoff framework, named causal directed acyclic graphs using deep learning variable selection (DAG-deepVASE). Using simulation data of diverse scenarios and identifying known and novel causal relations in molecular and clinical data of various diseases, we demonstrated that DAG-deepVASE consistently outperforms existing methods in identifying true and known causal relations. In the analyses, we also illustrate how identifying nonlinear causal relations and estimating their effect size help understand the complex disease pathobiology, which is not possible using other methods.

CONCLUSIONS

With these advantages, the application of DAG-deepVASE can help identify driver genes and therapeutic agents in biomedical studies and clinical trials.

From classical mendelian randomization to causal networks for systematic integration of multi-omics

The number of studies with information at multiple biological levels of granularity, such as genomics, proteomics, and metabolomics, is increasing each year, and a biomedical questaion is how to systematically integrate these data to discover new biological mechanisms that have the potential to elucidate the processes of health and disease. Causal frameworks, such as Mendelian randomization (MR), provide a foundation to begin integrating data for new biological discoveries. Despite the growing number of MR applications in a wide variety of biomedical studies, there are few approaches for the systematic analysis of omic data. The large number and diverse types of molecular components involved in complex diseases interact through complex networks, and classical MR approaches targeting individual components do not consider the underlying relationships. In contrast, causal network models established in the principles of MR offer significant improvements to the classical MR framework for understanding omic data. Integration of these mostly distinct branches of statistics is a recent development, and we here review the current progress. To set the stage for causal network models, we review some recent progress in the classical MR framework. We then explain how to transition from the classical MR framework to causal networks. We discuss the identification of causal networks and evaluate the underlying assumptions. We also introduce some tests for sensitivity analysis and stability assessment of causal networks. We then review practical details to perform real data analysis and identify causal networks and highlight some of the utility of causal networks. The utilities with validated novel findings reveal the full potential of causal networks as a systems approach that will become necessary to integrate large-scale omic data.

Whole Genome Association Study of the Plasma Metabolome Identifies Metabolites Linked to Cardiometabolic Disease in Black Individuals

Integrating genetic information with metabolomics has provided new insights into genes affecting human metabolism. However, gene-metabolite integration has been primarily studied in individuals of European Ancestry, limiting the opportunity to leverage genomic diversity for discovery. In addition, these analyses have principally involved known metabolites, with the majority of the profiled peaks left unannotated. Here, we perform a whole genome association study of 2,291 metabolite peaks (known and unknown features) in 2,466 Black individuals from the Jackson Heart Study. We identify 519 locus-metabolite associations for 427 metabolite peaks and validate our findings in two multi-ethnic cohorts. A significant proportion of these associations are in ancestry specific alleles including findings in APOE, TTR and CD36. We leverage tandem mass spectrometry to annotate unknown metabolites, providing new insight into hereditary diseases including transthyretin amyloidosis and sickle cell disease. Our integrative omics approach leverages genomic diversity to provide novel insights into diverse cardiometabolic diseases.

In Silico Pleiotropy Analysis in KEGG Signaling Networks Using a Boolean Network Model

Pleiotropy, which refers to the ability of different mutations on the same gene to cause different pathological effects in human genetic diseases, is important in understanding system-level biological diseases. Although some biological experiments have been proposed, still little is known about pleiotropy on gene–gene dynamics, since most previous studies have been based on correlation analysis. Therefore, a new perspective is needed to investigate pleiotropy in terms of gene–gene dynamical characteristics. To quantify pleiotropy in terms of network dynamics, we propose a measure called in silico Pleiotropic Scores (sPS), which represents how much a gene is affected against a pair of different types of mutations on a Boolean network model. We found that our model can identify more candidate pleiotropic genes that are not known to be pleiotropic than the experimental database. In addition, we found that many types of functionally important genes tend to have higher sPS values than other genes; in other words, they are more pleiotropic. We investigated the relations of sPS with the structural properties in the signaling network and found that there are highly positive relations to degree, feedback loops, and centrality measures. This implies that the structural characteristics are principles to identify new pleiotropic genes. Finally, we found some biological evidence showing that sPS analysis is relevant to the real pleiotropic data and can be considered a novel candidate for pleiotropic gene research. Taken together, our results can be used to understand the dynamics pleiotropic characteristics in complex biological systems in terms of gene–phenotype relations.

Whole-exome sequencing identifies rare genetic variants associated with human plasma metabolites

Metabolite levels measured in the human population are endophenotypes for biological processes. We combined sequencing data for 3,924 (whole-exome sequencing, WES, discovery) and 2,805 (whole-genome sequencing, WGS, replication) donors from a prospective cohort of blood donors in England. We used multiple approaches to select and aggregate rare genetic variants (minor allele frequency [MAF] < 0.1%) in protein-coding regions and tested their associations with 995 metabolites measured in plasma by using ultra-high-performance liquid chromatography-tandem mass spectrometry. We identified 40 novel associations implicating rare coding variants (27 genes and 38 metabolites), of which 28 (15 genes and 28 metabolites) were replicated. We developed algorithms to prioritize putative driver variants at each locus and used mediation and Mendelian randomization analyses to test directionality at associations of metabolite and protein levels at the ACY1 locus. Overall, 66% of reported associations implicate gene targets of approved drugs or bioactive drug-like compounds, contributing to drug targets' validating efforts.

Clinical metabolomics for inborn errors of metabolism

Metabolism is a highly regulated process that provides nutrients to cells and essential building blocks for the synthesis of protein, DNA and other macromolecules. In healthy biological systems, metabolism maintains a steady state in which the concentrations of metabolites are relatively constant yet are subject to metabolic demands and environmental stimuli. Rare genetic disorders, such as inborn errors of metabolism (IEM), cause defects in regulatory enzymes or proteins leading to metabolic pathway disruption and metabolite accumulation or deficiency. Traditionally, the laboratory diagnosis of IEMs has been limited to analytical methods that target specific metabolites such as amino acids and acyl carnitines. This approach is effective as a screening method for the most common IEM disorders but lacks the comprehensive coverage of metabolites that is necessary to identify rare disorders that present with nonspecific clinical symptoms. Fortunately, advancements in technology and data analytics has introduced a new field of study called metabolomics which has allowed scientists to perform comprehensive metabolite profiling of biological systems to provide insight into mechanism of action and gene function. Since metabolomics seeks to measure all small molecule metabolites in a biological specimen, it provides an innovative approach to evaluating disease in patients with rare genetic disorders. In this review we provide insight into the appropriate application of metabolomics in clinical settings. We discuss the advantages and limitations of the method and provide details related to the technology, data analytics and statistical modeling required for metabolomic profiling of patients with IEMs.

Multivariate analysis of NMR-based metabolomic data

Nuclear magnetic resonance (NMR) spectroscopy allows for simultaneous detection of a wide range of metabolites and lipids. As metabolites act together in complex metabolic networks, they are often highly correlated, and optimal biological insight is achieved when using methods that take the correlation into account. For this reason, latent-variable-based methods, such as principal component analysis and partial least-squares discriminant analysis, are widely used in metabolomic studies. However, with increasing availability of larger population cohorts, and a shift from analysis of spectral data to using quantified metabolite levels, both more traditional statistical approaches and alternative machine learning methods have become more widely used. This review aims at providing an overview of the current state-of-the-art multivariate methods for the analysis of NMR-based metabolomic data as well as alternative methods, highlighting their strengths and limitations.

Multi-Omic Approaches to Identify Genetic Factors in Metabolic Syndrome

Metabolic syndrome (MetS) is a highly heritable disease and a major public health burden worldwide. MetS diagnosis criteria are met by the simultaneous presence of any three of the following: high triglycerides, low HDL/high LDL cholesterol, insulin resistance, hypertension, and central obesity. These diseases act synergistically in people suffering from MetS and dramatically increase risk of morbidity and mortality due to stroke and cardiovascular disease, as well as certain cancers. Each of these component features is itself a complex disease, as is MetS. As a genetically complex disease, genetic risk factors for MetS are numerous, but not very powerful individually, often requiring specific environmental stressors for the disease to manifest. When taken together, all sequence variants that contribute to MetS disease risk explain only a fraction of the heritable variance, suggesting additional, novel loci have yet to be discovered. In this article, we will give a brief overview on the genetic concepts needed to interpret genome-wide association studies (GWAS) and quantitative trait locus (QTL) data, summarize the state of the field of MetS physiological genomics, and to introduce tools and resources that can be used by the physiologist to integrate genomics into their own research on MetS and any of its component features. There is a wealth of phenotypic and molecular data in animal models and humans that can be leveraged as outlined in this article. Integrating these multi-omic QTL data for complex diseases such as MetS provides a means to unravel the pathways and mechanisms leading to complex disease and promise for novel treatments. © 2022 American Physiological Society. Compr Physiol 12:1-40, 2022.

Potential functional variants of KIAA genes are associated with breast cancer risk in a case control study

Background

KIAA genes identified in the Kazusa cDNA-sequencing project may play important roles in biological processes and are involved in carcinogenesis of many cancers. Genetic variants of KIAA genes are implicated in the abnormal expression of these genes and are linked to susceptibility of several human complex diseases.

Methods

The differentially expressed KIAA genes were screened and identified in The Cancer Genome Atlas (TCGA) database of breast cancer. A total of 48 variants located in the 28 KIAA genes were selected to investigate the associations between polymorphism and breast cancer in 1,032 cases and 1,063 cancer-free controls in a Chinese population.

Results

Two coding variants, which included a SNP rs2306369 in KIAA1109 and a SNP rs1205434 in KIAA1755, were identified to be associated with the incidences of breast cancer. Logistic regression analysis showed that the SNP rs2306369 G allele was associated with a decreased risk of breast cancer (additive model: OR =0.81, 95% CI: 0.66-0.99, P=0.038), whereas the SNP rs1205434 A allele was involved with a higher risk of breast cancer (additive model: OR =1.19, 95% CI: 1.02-1.38, P= 0.025). Further stratified analysis revealed that the SNP rs1205434 showed a significant difference for age at menarche strata (heterogeneity test P=0.009). Multiplicative interaction analysis indicated that there was positive multiplicative interaction between the SNP rs1205434 and menarche age (OR =1.09, 95% CI: 1.01-1.17, P=0.036). Additionally, expression quantitative trait loci analysis revealed that the SNP rs1205434 A allele could decrease the KIAA1755 expression in the Genotype-Tissue Expression (GTEx) database (P=0.002). The Kaplan-Meier plotter showed that breast cancer patients with high KIAA1755 expression have significantly better outcomes than those with low levels of expression (HR =0.84, 95% CI: 0.72-0.99, P=0.033).

Conclusions

The results indicate that the genetic variants (rs2306369 and rs1205434) in the coding region of KIAA1109 and KIAA1755 respectively may affect Chinese females' breast cancer susceptibility and act as potential predictive biomarkers for breast cancer.

Genomic-Metabolomic Associations Support the Role of LIPC and Glycerophospholipids in Age-Related Macular Degeneration

Purpose

Large-scale genome-wide association studies (GWAS) have reported important single nucleotide polymorphisms (SNPs) with significant associations with age-related macular degeneration (AMD). However, their role in disease development remains elusive. This study aimed to assess SNP–metabolite associations (i.e., metabolite quantitative trait loci [met-QTL]) and to provide insights into the biological mechanisms of AMD risk SNPs.

Design

Cross-sectional multicenter study (Boston, Massachusetts, and Coimbra, Portugal).

Participants

Patients with AMD (n = 388) and control participants (n = 98) without any vitreoretinal disease (> 50 years).

Methods

Age-related macular degeneration grading was performed using color fundus photographs according to the Age-Related Eye Disease Study classification scheme. Fasting blood samples were collected and evaluated with mass spectrometry for metabolomic profiling and Illumina OmniExpress for SNPs profiling. Analyses of met-QTL of endogenous metabolites were conducted using linear regression models adjusted for age, gender, smoking, 10 metabolite principal components (PCs), and 10 SNP PCs. Additionally, we analyzed the cumulative effect of AMD risk SNPs on plasma metabolites by generating genetic risk scores and assessing their associations with metabolites using linear regression models, accounting for the same covariates. Modeling was performed first for each cohort, and then combined by meta-analysis. Multiple comparisons were accounted for using the false discovery rate (FDR).

Main Outcome Measures

Plasma metabolite levels associated with AMD risk SNPs.

Results

After quality control, data for 544 plasma metabolites were included. Meta-analysis of data from all individuals (AMD patients and control participants) identified 28 significant met-QTL (β = 0.016–0.083; FDR q-value < 1.14 × 10^–2), which corresponded to 5 metabolites and 2 genes: ASPM and LIPC. Polymorphisms in the LIPC gene were associated with phosphatidylethanolamine metabolites, which are glycerophospholipids, and polymorphisms in the ASPM gene with branched-chain amino acids. Similar results were observed when considering only patients with AMD. Genetic risk score–metabolite associations further supported a global impact of AMD risk SNPs on the plasma metabolome.

Conclusions

This study demonstrated that genomic–metabolomic associations can provide insights into the biological relevance of AMD risk SNPs. In particular, our results support that the LIPC gene and the glycerophospholipid metabolic pathway may play an important role in AMD, thus offering new potential therapeutic targets for this disease.

Steroid disulfates - Sulfation double trouble

Sulfation pathways have recently come into the focus of biomedical research. For steroid hormones and related compounds, sulfation represents an additional layer of regulation as sulfated steroids are more water-soluble and tend to be biologically less active. For steroid diols, an additional sulfation is possible, carried out by the same sulfotransferases that catalyze the first sulfation step. The steroid disulfates that are formed are the focus of this review. We discuss both their biochemical production as well as their putative biological function. Steroid disulfates have also been linked to various clinical conditions in numerous untargeted metabolomics studies. New analytical techniques exploring the biosynthetic routes of steroid disulfates have led to novel insights, changing our understanding of sulfation in human biology. They promise a bright future for research into sulfation pathways, hopefully too for the diagnosis and treatment of several associated diseases.

Integrating lipidomics and genomics: emerging tools to understand cardiovascular diseases

Cardiovascular diseases (CVDs) are the leading cause of mortality and morbidity worldwide leading to 31% of all global deaths. Early prediction and prevention could greatly reduce the enormous socio-economic burden posed by CVDs. Plasma lipids have been at the center stage of the prediction and prevention strategies for CVDs that have mostly relied on traditional lipids (total cholesterol, total triglycerides, HDL-C and LDL-C). The tremendous advancement in the field of lipidomics in last two decades has facilitated the research efforts to unravel the metabolic dysregulation in CVDs and their genetic determinants, enabling the understanding of pathophysiological mechanisms and identification of predictive biomarkers, beyond traditional lipids. This review presents an overview of the application of lipidomics in epidemiological and genetic studies and their contributions to the current understanding of the field. We review findings of these studies and discuss examples that demonstrates the potential of lipidomics in revealing new biology not captured by traditional lipids and lipoprotein measurements. The promising findings from these studies have raised new opportunities in the fields of personalized and predictive medicine for CVDs. The review further discusses prospects of integrating emerging genomics tools with the high-dimensional lipidome to move forward from the statistical associations towards biological understanding, therapeutic target development and risk prediction. We believe that integrating genomics with lipidome holds a great potential but further advancements in statistical and computational tools are needed to handle the high-dimensional and correlated lipidome.

Genetics and Not Shared Environment Explains Familial Resemblance in Adult Metabolomics Data

Metabolites are small molecules involved in cellular metabolism where they act as reaction substrates or products. The term 'metabolomics' refers to the comprehensive study of these molecules. The concentrations of metabolites in biological tissues are under genetic control, but this is limited by environmental factors such as diet. In adult mono- and dizygotic twin pairs, we estimated the contribution of genetic and shared environmental influences on metabolite levels by structural equation modeling and tested whether the familial resemblance for metabolite levels is mainly explained by genetic or by environmental factors that are shared by family members. Metabolites were measured across three platforms: two based on proton nuclear magnetic resonance techniques and one employing mass spectrometry. These three platforms comprised 237 single metabolic traits of several chemical classes. For the three platforms, metabolites were assessed in 1407, 1037 and 1116 twin pairs, respectively. We carried out power calculations to establish what percentage of shared environmental variance could be detected given these sample sizes. Our study did not find evidence for a systematic contribution of shared environment, defined as the influence of growing up together in the same household, on metabolites assessed in adulthood. Significant heritability was observed for nearly all 237 metabolites; significant contribution of the shared environment was limited to 6 metabolites. The top quartile of the heritability distribution was populated by 5 of the 11 investigated chemical classes. In this quartile, metabolites of the class lipoprotein were significantly overrepresented, whereas metabolites of classes glycerophospholipids and glycerolipids were significantly underrepresented.