A network-based conditional genetic association analysis of the human metabolome

Workability of mRNA Sequencing for Predicting Protein Abundance

Transcriptomics methods (RNA-Seq, PCR) today are more routine and reproducible than proteomics methods, i.e., both mass spectrometry and immunochemical analysis. For this reason, most scientific studies are limited to assessing the level of mRNA content. At the same time, protein content (and its post-translational status) largely determines the cell's state and behavior. Such a forced extrapolation of conclusions from the transcriptome to the proteome often seems unjustified. The ratios of "transcript-protein" pairs can vary by several orders of magnitude for different genes. As a rule, the correlation coefficient between transcriptome-proteome levels for different tissues does not exceed 0.3-0.5. Several characteristics determine the ratio between the content of mRNA and protein: among them, the rate of movement of the ribosome along the mRNA and the number of free ribosomes in the cell, the availability of tRNA, the secondary structure, and the localization of the transcript. The technical features of the experimental methods also significantly influence the levels of the transcript and protein of the corresponding gene on the outcome of the comparison. Given the above biological features and the performance of experimental and bioinformatic approaches, one may develop various models to predict proteomic profiles based on transcriptomic data. This review is devoted to the ability of RNA sequencing methods for protein abundance prediction.

Development and validation of robust metabolism-related gene signature in the prognostic prediction of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the most common malignant tumours worldwide. Given metabolic reprogramming in tumours was a crucial hallmark, several studies have demonstrated its value in the diagnostics and surveillance of malignant tumours. The present study aimed to identify a cluster of metabolism-related genes to construct a prediction model for the prognosis of HCC. Multiple cohorts of HCC cases (466 cases) from public datasets were included in the present analysis. (GEO cohort) After identifying a list of metabolism-related genes associated with prognosis, a risk score based on metabolism-related genes was formulated via the LASSO-Cox and LASSO-pcvl algorithms. According to the risk score, patients were stratified into low- and high-risk groups, and further analysis and validation were accordingly conducted. The results revealed that high-risk patients had a significantly worse 5-year overall survival (OS) than low-risk patients in the GEO cohort. (30.0% vs. 57.8%; hazard ratio [HR], 0.411; 95% confidence interval [95% CI], 0.302-0.651; p < 0.001) This observation was confirmed in the external TCGA-LIHC cohort. (34.5% vs. 54.4%; HR 0.452; 95% CI, 0.299-0.681; p < 0.001) To promote the predictive ability of the model, risk score, age, gender and tumour stage were integrated into a nomogram. According to the results of receiver operating characteristic curves and decision curves analysis, the nomogram score possessed a superior predictive ability than conventional factors, which indicate that the risk score combined with clinicopathological features was able to achieve a robust prediction for OS and improve the individualized clinical decision making of HCC patients. In conclusion, the metabolic genes related to OS were identified and developed a metabolism-based predictive model for HCC. Through a series of bioinformatics and statistical analyses, the predictive ability of the model was approved.

Association of Physical Activity With Bioactive Lipids and Cardiovascular Events

BACKGROUND

To clarify the mechanisms underlying physical activity (PA)-related cardioprotection, we examined the association of PA with plasma bioactive lipids (BALs) and cardiovascular disease (CVD) events. We additionally performed genome-wide associations.

METHODS

PA-bioactive lipid associations were examined in VITAL (VITamin D and OmegA-3 TriaL)-clinical translational science center (REGISTRATION: URL: https://www.

CLINICALTRIALS

gov; Unique identifier: NCT01169259; N=1032) and validated in JUPITER (Justification for the Use of statins in Prevention: an Intervention Trial Evaluating Rosuvastatin)-NC (NCT00239681; N=589), using linear models adjusted for age, sex, race, low-density lipoprotein-cholesterol, total-C, and smoking. Significant BALs were carried over to examine associations with incident CVD in 2 nested CVD case-control studies: VITAL-CVD (741 case-control pairs) and JUPITER-CVD (415 case-control pairs; validation).

RESULTS

We detected 145 PA-bioactive lipid validated associations (false discovery rate <0.1). Annotations were found for 6 of these BALs: 12,13-diHOME, 9,10-diHOME, lysoPC(15:0), oxymorphone-3b-D-glucuronide, cortisone, and oleoyl-glycerol. Genetic analysis within JUPITER-NC showed associations of 32 PA-related BALs with 22 single-nucleotide polymorphisms. From PA-related BALs, 12 are associated with CVD.

CONCLUSIONS

We identified a PA-related bioactive lipidome profile out of which 12 BALs also had opposite associations with incident CVD events.

A protocol for recruiting and analyzing the disease-oriented Russian disc degeneration study (RuDDS) biobank for functional omics studies of lumbar disc degeneration

Lumbar intervertebral disc degeneration (DD) disease is one of the main risk factors for low back pain and a leading cause of population absenteeism and disability worldwide. Despite a variety of biological studies, lumbar DD is not yet fully understood, partially because there are only few studies that use systematic and integrative approaches. This urges the need for studies that integrate different omics (including genomics and transcriptomics) measured on samples within a single cohort. This protocol describes a disease-oriented Russian disc degeneration study (RuDDS) biobank recruitment and analyses aimed to facilitate further omics studies of lumbar DD integrating genomic, transcriptomic and glycomic data. A total of 1,100 participants aged over 18 with available lumbar MRI scans, medical histories and biological material (whole blood, plasma and intervertebral disc tissue samples from surgically treated patients) will be enrolled during the three-year period from two Russian clinical centers. Whole blood, plasma and disc tissue specimens will be used for genotyping with genome-wide SNP-arrays, glycome profiling and RNA sequencing, respectively. Omics data will be further used for a genome-wide association study of lumbar DD with in silico functional annotation, analysis of plasma glycome and lumbar DD disease interactions and transcriptomic data analysis including an investigation of differential expression patterns associated with lumbar DD disease. Statistical tests applied in each of the analyses will meet the standard criteria specific to the attributed study field. In a long term, the results of the study will expand fundamental knowledge about lumbar DD development and contribute to the elaboration of novel personalized approaches for disease prediction and therapy. Additionally to the lumbar disc degeneration study, a RuDDS cohort could be used for other genetic studies, as it will have unique omics data.

Trial registration numberNCT04600544.

Defining Blood Plasma and Serum Metabolome by GC-MS

Metabolomics uses advanced analytical chemistry methods to analyze metabolites in biological samples. The most intensively studied samples are blood and its liquid components: plasma and serum. Armed with advanced equipment and progressive software solutions, the scientific community has shown that small molecules’ roles in living systems are not limited to traditional “building blocks” or “just fuel” for cellular energy. As a result, the conclusions based on studying the metabolome are finding practical reflection in molecular medicine and a better understanding of fundamental biochemical processes in living systems. This review is not a detailed protocol of metabolomic analysis. However, it should support the reader with information about the achievements in the whole process of metabolic exploration of human plasma and serum using mass spectrometry combined with gas chromatography.

Multivariate genome-wide analysis of immunoglobulin G N-glycosylation identifies new loci pleiotropic with immune function

The N-glycosylation of immunoglobulin G (IgG) affects its structure and function. It has been demonstrated that IgG N-glycosylation patterns are inherited as complex quantitative traits. Genome-wide association studies identified loci harboring genes encoding enzymes directly involved in protein glycosylation as well as loci likely to be involved in regulation of glycosylation biochemical pathways. Many of these loci could be linked to immune functions and risk of inflammatory and autoimmune diseases. The aim of the present study was to discover and replicate new loci associated with IgG N-glycosylation and to investigate possible pleiotropic effects of these loci onto immune function and the risk of inflammatory and autoimmune diseases. We conducted a multivariate genome-wide association analysis of 23 IgG N-glycosylation traits measured in 8090 individuals of European ancestry. The discovery stage was followed up by replication in 3147 people and in silico functional analysis. Our study increased the total number of replicated loci from 22 to 29. For the discovered loci, we suggest a number of genes potentially involved in the control of IgG N-glycosylation. Among the new loci, two (near RNF168 and TNFRSF13B) were previously implicated in rare immune deficiencies and were associated with levels of circulating immunoglobulins. For one new locus (near AP5B1/OVOL1), we demonstrated a potential pleiotropic effect on the risk of asthma. Our findings underline an important link between IgG N-glycosylation and immune function and provide new clues to understanding their interplay.

Varicose veins of lower extremities: Insights from the first large-scale genetic study

Varicose veins of lower extremities (VVs) are a common multifactorial vascular disease. Genetic factors underlying VVs development remain largely unknown. Here we report the first large-scale study of VVs performed on a freely available genetic data of 408,455 European-ancestry individuals. We identified the 12 reliably associated loci that explain 13% of the SNP-based heritability, and prioritized the most likely causal genes CASZ1, PIEZO1, PPP3R1, EBF1, STIM2, HFE, GATA2, NFATC2, and SOX9. VVs-associated variants within these loci exhibited pleiotropic effects on several phenotypes including blood pressure/hypertension and blood cell traits. Gene set enrichment analysis revealed gene categories related to abnormal vasculogenesis. Genetic correlation analysis confirmed known epidemiological associations between VVs and deep venous thrombosis, weight, rough labor, and standing job, and found a genetic overlap with multiple traits that have not been previously suspected to share common genetic background with VVs. These traits included educational attainment, fluid intelligence and prospective memory scores, walking pace (negative correlation with VVs), smoking, height, number of operations, pain, and gonarthrosis (positive correlation with VVs). Finally, Mendelian randomization analysis provided evidence for causal effects of plasma levels of MICB and CD209 proteins, and anthropometric traits such as waist and hip circumference, height, weight, and both fat and fat-free mass. Our results provide novel insight into both VVs genetics and etiology. The revealed genes and proteins can be considered as good candidates for follow-up functional studies and might be of interest as potential drug targets.

Varicose veins of lower extremities (VVs) affect about 30% of adults in developed countries and cause both cosmetic and health problems. A strong body of evidence indicates that heredity plays an important role in the etiology of this condition. However, genetic basis of VVs remains poorly understood. Here, we present the results of the first large-scale genetic study for VVs. We identified genes which are the most likely involved in VVs pathogenesis. We show that VVs are correlated at a genetic level with numerous traits and phenotypes, including those already known from prior epidemiological studies (deep venous thrombosis, body mass index, standing job, etc.) as well as with those that have not been suspected to share common genetic background with VVs (fluid intelligence and prospective memory scores, smoking, walking pace, pain all over the body, and other traits). Finally, using genetic variants as instruments, we demonstrate direct causal effects of the traits related to anthropometry, such as height and weight, and plasma levels of immune-related proteins MICB and CD209. Our study provides novel insight into both VVs genetics and etiology. The revealed genes (CASZ1, PIEZO1, PPP3R1, EBF1, STIM2, HFE, GATA2, NFATC2, and SOX9) and proteins (MICB and CD209) can be considered as good candidates for follow-up functional studies and might be of interest as potential drug targets.

A network-based conditional genetic association analysis of the human metabolome

Background

Genome-wide association studies have identified hundreds of loci that influence a wide variety of complex human traits; however, little is known regarding the biological mechanism of action of these loci. The recent accumulation of functional genomics (“omics”), including metabolomics data, has created new opportunities for studying the functional role of specific changes in the genome. Functional genomic data are characterized by their high dimensionality, the presence of (strong) statistical dependency between traits, and, potentially, complex genetic control. Therefore, the analysis of such data requires specific statistical genetics methods.

Results

To facilitate our understanding of the genetic control of omics phenotypes, we propose a trait-centered, network-based conditional genetic association (cGAS) approach for identifying the direct effects of genetic variants on omics-based traits. For each trait of interest, we selected from a biological network a set of other traits to be used as covariates in the cGAS. The network can be reconstructed either from biological pathway databases (a mechanistic approach) or directly from the data, using a Gaussian graphical model applied to the metabolome (a data-driven approach). We derived mathematical expressions that allow comparison of the power of univariate analyses with conditional genetic association analyses. We then tested our approach using data from a population-based Cooperative Health Research in the region of Augsburg (KORA) study (n = 1,784 subjects, 1.7 million single-nucleotide polymorphisms) with measured data for 151 metabolites.

Conclusions

We found that compared to single-trait analysis, performing a genetic association analysis that includes biologically relevant covariates can either gain or lose power, depending on specific pleiotropic scenarios, for which we provide empirical examples. In the context of analyzed metabolomics data, the mechanistic network approach had more power compared to the data-driven approach. Nevertheless, we believe that our analysis shows that neither a prior-knowledge-only approach nor a phenotypic-data-only approach is optimal, and we discuss possibilities for improvement.