Search

Obesity is a socially significant disease that is characterized by a disproportionate accumulation of fat. It is also associated with chronic inflammation, cancer, diabetes, and other comorbidities. Investigating biomarkers and pathological processes linked to obesity is especially vital for young individuals, given their increased potential for lifestyle modifications. By comparing the genetic, proteomic, and metabolomic profiles of individuals categorized as underweight, normal, overweight, and obese, we aimed to determine which omics layer most accurately reflects the phenotypic changes in an organism that result from obesity. We profiled blood plasma samples by employing three omics methodologies. The untargeted GC×GC–MS metabolomics approach identified 313 metabolites. To augment the metabolomic dataset, we integrated a label-free HPLC–MS/MS proteomics method, leading to the identification of 708 proteins. The genomic layer encompassed the genotyping of 647,250 SNPs. Utilizing omics data, we trained sparse Partial Least Squares models to predict body mass index. Molecular features exhibiting frequently non-zero coefficients were selected as potential biomarkers, and we further explored enriched biological pathways. Proteomics was the most effective in single-omics analyses, with a median absolute error (MAE) of 5.44 ± 0.31 kg/m2, incorporating an average of 24 proteins per model. Metabolomics showed slightly lower performance (MAE = 6.06 ± 0.33 kg/m2), followed by genomics (MAE = 6.20 ± 0.34 kg/m2). As expected, multiomic models demonstrated better accuracy, particularly the combination of proteomics and metabolomics (MAE = 4.77 ± 0.33 kg/m2), while including genomics data did not enhance the results. This manuscript is the first multiomics study of obesity in a gender-balanced cohort of young adults profiled by genomic, proteomic, and metabolomic methods. The comprehensive approach provides novel insights into the molecular mechanisms of obesity, opening avenues for more targeted interventions.

Recently, the concept of a mass spectrometric blood metabogram was introduced, which allows the analysis of the blood metabolome in terms of the time, cost, and reproducibility of clinical laboratory tests. It was demonstrated that the components of the metabogram are related groups of the blood metabolites associated with humoral regulation; the metabolism of lipids, carbohydrates, and amines; lipid intake into the organism; and liver function, thereby providing clinically relevant information. The purpose of this work was to evaluate the relevance of using the metabogram in a disease. To do this, the metabogram was used to analyze patients with various degrees of metabolic alterations associated with obesity. The study involved 20 healthy individuals, 20 overweight individuals, and 60 individuals with class 1, 2, or 3 obesity. The results showed that the metabogram revealed obesity-associated metabolic alterations, including changes in the blood levels of steroids, amino acids, fatty acids, and phospholipids, which are consistent with the available scientific data to date. Therefore, the metabogram allows testing of metabolically unhealthy overweight or obese patients, providing both a general overview of their metabolic alterations and detailing their individual characteristics. It was concluded that the metabogram is an accurate and clinically applicable test for assessing an individual’s metabolic status in disease.

Although modern biology is now in the post-genomic era with vastly increased access to high-quality data, the set of human genes with a known function remains far from complete. This is especially true for hundreds of mitochondria-associated genes, which are under-characterized and lack clear functional annotation. However, with the advent of multi-omics profiling methods coupled with systems biology algorithms, the cellular role of many such genes can be elucidated. Here, we report genes and pathways associated with TOMM34, Translocase of Outer Mitochondrial Membrane, which plays role in the mitochondrial protein import as a part of cytosolic complex together with Hsp70/Hsp90 and is upregulated in various cancers. We identified genes, proteins, and metabolites altered in TOMM34-/- HepG2 cells. To our knowledge, this is the first attempt to study the functional capacity of TOMM34 using a multi-omics strategy. We demonstrate that TOMM34 affects various processes including oxidative phosphorylation, citric acid cycle, metabolism of purine, and several amino acids. Besides the analysis of already known pathways, we utilized de novo network enrichment algorithm to extract novel perturbed subnetworks, thus obtaining evidence that TOMM34 potentially plays role in several other cellular processes, including NOTCH-, MAPK-, and STAT3-signaling. Collectively, our findings provide new insights into TOMM34’s cellular functions.

Both biological and technical variations can discredit the reliability of obtained data in omics studies. In this technical note, we investigated the effect of prolonged cultivation of the HepG2 hepatoma cell line on its metabolomic profile. Using the GC × GC-MS approach, we determined the degree of metabolic variability across HepG2 cells cultured in uniform conditions for 0, 5, 10, 15, and 20 days. Post-processing of obtained data revealed substantial changes in relative abundances of 110 metabolites among HepG2 samples under investigation. Our findings have implications for interpreting metabolomic results obtained from immortal cells, especially in longitudinal studies. There are still plenty of unanswered questions regarding metabolomics variability and many potential areas for future targeted and panoramic research. However, we suggest that the metabolome of cell lines is unstable and may undergo significant transformation over time, even if the culture conditions remain the same. Considering metabolomics variability on a relatively long-term basis, careful experimentation with particular attention to control samples is required to ensure reproducibility and relevance of the research results when testing both fundamentally and practically significant hypotheses.

The three-dimensional structure of monomers and homodimers of CYP102A1/WT (wild-type) proteins and their A83F and A83I mutant forms was predicted using the AlphaFold2 (AF2) and AlphaFold Multimer (AFMultimer) programs, which were compared with the rate constants of hydroxylation reactions of these enzyme forms to determine the efficiency of intra- and interprotein electron transport in the CYP102A1 hydroxylase system. The electron transfer rate constants (ket), which determine the rate of indole hydroxylation by the CYP102A1 system, were calculated based on the distances (R) between donor-acceptor prosthetic groups (PG) FAD→FMN→HEME of these proteins using factor β, which describes an exponential decay from R the speed of electron transport (ET) according to the tunnelling mechanism. It was shown that the structure of monomers in the homodimer, calculated using the AlpfaFold Multimer program, is in good agreement with the experimental structures of globular domains (HEME-, FMN-, and FAD-domains) in CYP102A1/WT obtained by X-ray structural analysis, and the structure of isolated monomers predicted in AF2 does not coincide with the structure of monomers in the homodimer, although a high level of similarity in individual domains remains. The structures of monomers and homodimers of A83F and A83I mutants were also calculated, and their structures were compared with the wild-type protein. Significant differences in the structure of all isolated monomers with respect to the structures of monomers in homodimers were also found for them, and at the same time, insignificant differences were revealed for all homodimers. Comparative analysis for CYP102A1/WT between the calculated intra- and interprotein distances FAD→FMN→HEME and the rate constants of hydroxylation in these proteins showed that the distance between prosthetic groups both in the monomer and in the dimer allows the implementation of electron transfer between PGs, which is consistent with experimental literature data about kcat. For the mutant form of monomer A83I, an increase in the distance between PGs was obtained, which can restrict electron transportation compared to WT; however, for the dimer of this protein, a decrease in the distance between PGs was observed compared to the WT form, which can lead to an increase in the electron transfer rate constant and, accordingly, kcat. For the monomer and homodimer of the A83F mutant, the calculations showed an increase in the distance between the PGs compared to the WT form, which should have led to a decrease in the electron transfer rate, but at the same time, for the homodimer, the approach of the aromatic group F262 with heme can speed up transportation for this form and, accordingly, the rate of hydroxylation.

Long-read direct RNA sequencing developed by Oxford Nanopore Technologies (ONT) is quickly gaining popularity for transcriptome studies, while fast turnaround time and low cost make it an attractive instrument for clinical applications. There is a growing interest to utilize transcriptome data to unravel activated biological processes responsible for disease progression and response to therapies. This trend is of particular interest for precision medicine which aims at single-patient analysis. Here we evaluated whether gene abundances measured by MinION direct RNA sequencing are suited to produce robust estimates of pathway activation for single sample scoring methods. We performed multiple RNA-seq analyses for a single sample that originated from the HepG2 cell line, namely five ONT replicates, and three replicates using Illumina NovaSeq. Two pathway scoring methods were employed—ssGSEA and singscore. We estimated the ONT performance in terms of detected protein-coding genes and average pairwise correlation between pathway activation scores using an exhaustive computational scheme for all combinations of replicates. In brief, we found that at least two ONT replicates are required to obtain reproducible pathway scores for both algorithms. We hope that our findings may be of interest to researchers planning their ONT direct RNA-seq experiments.

The Human Proteome Project is moving into the next phase of creating and/or reconsidering the functional annotations of proteins using the chromosome-centric paradigm. This challenge cannot be solved exclusively using automated means, but rather requires human intelligence for interpreting the combined data. To foster the integration between human cognition and post-genome array a number of specific tools were recently developed, among them CAPER, GenomewidePDB, and The Proteome Browser (TPB). For the purpose of tackling the task of protein functional annotating the Gene-Centric Content Management System (GenoCMS) was expanded with new features. The goal was to enable bioinformaticans to develop self-made applications and to position these applets within the generalized informational canvas supported by GenoCMS. We report the results of GenoCMS-enabled integration of the concordant informational flows in the chromosome-centric framework of the human chromosome 18 project. The workflow described in the article can be scaled to other human chromosomes, and also supplemented with new tracks created by the user. The GenoCMS is an example of a project-oriented informational system, which are important for public data sharing.