Exploiting Multi-Omics Profiling and Systems Biology to Investigate Functions of TOMM34.

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.

Exploring Dynamic Metabolome of the HepG2 Cell Line: Rise and Fall.

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.

Oxford Nanopore MinION Direct RNA-Seq for Systems Biology.

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.

Heterotypic Multicellular Spheroids as Experimental and Preclinical Models of Sprouting Angiogenesis.

Sprouting angiogenesis is the common response of live tissues to physiological and pathological angiogenic stimuli. Its accurate evaluation is of utmost importance for basic research and practical medicine and pharmacology and requires adequate experimental models. A variety of assays for angiogenesis were developed, none of them perfect. In vitro approaches are generally less physiologically relevant due to the omission of essential components regulating the process. However, only in vitro models can be entirely non-xenogeneic. The limitations of the in vitro angiogenesis assays can be partially overcome using 3D models mimicking tissue O2 and nutrient gradients, the influence of the extracellular matrix (ECM), and enabling cell-cell interactions. Here we present a review of the existing models of sprouting angiogenesis that are based on the use of endothelial cells (ECs) co-cultured with perivascular or other stromal cells. This approach provides an excellent in vitro platform for further decoding of the cellular and molecular mechanisms of sprouting angiogenesis under conditions close to the in vivo conditions, as well as for preclinical drug testing and preclinical research in tissue engineering and regenerative medicine.

Magnetic levitational bioassembly of 3D tissue construct in space.

Magnetic levitational bioassembly of three-dimensional (3D) tissue constructs represents a rapidly emerging scaffold- and label-free approach and alternative conceptual advance in tissue engineering. The magnetic bioassembler has been designed, developed, and certified for life space research. To the best of our knowledge, 3D tissue constructs have been biofabricated for the first time in space under microgravity from tissue spheroids consisting of human chondrocytes. Bioassembly and sequential tissue spheroid fusion presented a good agreement with developed predictive mathematical models and computer simulations. Tissue constructs demonstrated good viability and advanced stages of tissue spheroid fusion process. Thus, our data strongly suggest that scaffold-free formative biofabrication using magnetic fields is a feasible alternative to traditional scaffold-based approaches, hinting a new perspective avenue of research that could significantly advance tissue engineering. Magnetic levitational bioassembly in space can also advance space life science and space regenerative medicine.

Fabrication of calcium phosphate 3D scaffolds for bone repair using magnetic levitational assembly.

The calcium phosphate particles can be used as building blocks for fabrication of 3D scaffolds intended for bone tissue engineering. This work presents for the first time a rapid creation of 3D scaffolds using magnetic levitation of calcium phosphate particles. Namely, tricalcium phosphate particles of equal size and certain porosity are used, which undergo the process of recrystallization after magnetic levitational assembly of the scaffold to ensure stitching of the scaffold. Label-free levitational assembly is achieved by using a custom-designed magnetic system in the presence of gadolinium salts, which allows the levitation of calcium phosphate particles. Chemical transformation of tricalcium- to octacalcium phosphate under the condition of magnetic levitation in non-homogeneous magnetic field is also demonstrated. This approach allows obtaining rapidly the octacalcium phosphate phase in the final 3D product, which is biocompatible.

Application of selected reaction monitoring and parallel reaction monitoring for investigation of HL-60 cell line differentiation.

Targeted mass spectrometry represents a powerful tool for investigation of biological processes. The convenient approach of selected reaction monitoring using stable isotope-labeled peptide standard (SIS) is widely applied for protein quantification. Along with this method, high-resolution parallel reaction monitoring has been increasingly used for protein targeted analysis. Here we applied two targeted approaches (selected reaction monitoring with SIS and label-free parallel reaction monitoring) to investigate expression of 11 proteins during all-trans retinoic acid-induced differentiation of HL-60 cells. In our experiments, we have determined the proteins expression ratio at 3, 24, 48, and 96 h after all-trans retinoic acid treatment in comparison with 0 h, respectively. Expression profiles of four proteins (VAV1, PRAM1, LYN, and CEBPB) were highly correlated ( r > 0.75) and FGR expression was detected on proteome level starting from 24 h by both techniques. For prominent differences (fold change ≥ 2) label-free parallel reaction monitoring is not inferior to selected reaction monitoring with isotopically labeled peptide standards. Differentially expressed proteins, that have been determined in our study, can be considered as potential drug targets for acute myeloid leukemia (AML) treatment.

State of the Art of Chromosome 18-Centric HPP in 2016: Transcriptome and Proteome Profiling of Liver Tissue and HepG2 Cells.

A gene-centric approach was applied for a large-scale study of expression products of a single chromosome. Transcriptome profiling of liver tissue and HepG2 cell line was independently performed using two RNA-Seq platforms (SOLiD and Illumina) and also by Droplet Digital PCR (ddPCR) and quantitative RT-PCR. Proteome profiling was performed using shotgun LC-MS/MS as well as selected reaction monitoring with stable isotope-labeled standards (SRM/SIS) for liver tissue and HepG2 cells. On the basis of SRM/SIS measurements, protein copy numbers were estimated for the Chromosome 18 (Chr 18) encoded proteins in the selected types of biological material. These values were compared with expression levels of corresponding mRNA. As a result, we obtained information about 158 and 142 transcripts for HepG2 cell line and liver tissue, respectively. SRM/SIS measurements and shotgun LC-MS/MS allowed us to detect 91 Chr 18-encoded proteins in total, while an intersection between the HepG2 cell line and liver tissue proteomes was ∼66%. In total, there were 16 proteins specifically observed in HepG2 cell line, while 15 proteins were found solely in the liver tissue. Comparison between proteome and transcriptome revealed a poor correlation (R2 ≈ 0.1) between corresponding mRNA and protein expression levels. The SRM and shotgun data sets (obtained during 2015-2016) are available in PASSEL (PASS00697) and ProteomeExchange/PRIDE (PXD004407). All measurements were also uploaded into the in-house Chr 18 Knowledgebase at http://kb18.ru/protein/matrix/416126 .

Virtual-Experimental 2DE Approach in Chromosome-Centric Human Proteome Project.

To obtain more information about human proteome, especially about proteoforms (protein species) coded by 18th chromosome, we separated proteins from human cancer cell line (HepG2) by two-dimensional gel electrophoresis (2DE). Initially, proteins in major spots were identified by MALDI-MS peptide mass fingerprinting. According to parameters (pI/Mw) of identified proteins the gel was calibrated. Using this calibrated gel, a virtual 2D map of proteoforms coded by Chromosome 18 was constructed. Next, the produced gel was divided into 96 sections with determined coordinates. Each section was cut, shredded, and treated by trypsin according to mass-spectrometry protocol. After protein identification by shotgun mass spectrometry using ESI LC-MS/MS, a list of 20 462 proteoforms (product of 3774 genes) was generated. Among them, 165 proteoforms are representing 39 genes of 18th chromosome. The 3D graphs showing the distribution of different proteoforms from the same gene in 2D map were generated. This is a first step in creation of 2DE-based knowledge database of proteins coded by 18th chromosome.

Combination of virtual and experimental 2DE together with ESI LC-MS/MS gives a clearer view about proteomes of human cells and plasma.

Virtual and experimental 2DE coupled with ESI LC-MS/MS was introduced to obtain better representation of the information about human proteome. The proteins from HEPG2 cells and human blood plasma were run by 2DE. After staining and protein spot identification by MALDI-TOF MS, the protein maps were generated. The experimental physicochemical parameters (pI/Mw) of the proteoforms further detected by ESI LC-MS/MS in these spots were obtained. Next, the theoretical pI and Mw of identified proteins were calculated using program Compute pI/Mw (http://web.expasy.org/compute_pi/pi_tool-doc.html). Accordingly, the relationship between theoretical and experimental parameters was analyzed, and the correlation plots were built. Additionally, virtual/experimental information about different protein species/proteoforms from the same genes was extracted. As it was revealed from the plots, the major proteoforms detected in HepG2 cell line have pI/Mw parameters similar to theoretical values. In opposite, the minor protein species have mainly very different from theoretical pI and Mw parameters. A similar situation was observed in plasma in much higher degree. It means that minor protein species are heavily modified in cell and even more in plasma proteome.

RNA-Seq gene expression profiling of HepG2 cells: the influence of experimental factors and comparison with liver tissue.

BACKGROUND: Human hepatoma HepG2 cells are used as an in vitro model of the human liver. High-throughput transcriptomic sequencing is an advanced approach for assessing the functional state of a tissue or cell type. However, the influence of experimental factors, such as the sample preparation method and inter-laboratory variation, on the transcriptomic profile has not been evaluated. RESULTS: The whole-transcriptome sequencing of HepG2 cells was performed using the SOLiD platform and validated using droplet digital PCR. The gene expression profile was compared to the results obtained with the same sequencing method in another laboratory and using another sample preparation method. We also compared the transcriptomic profile HepG2 cells with that of liver tissue. Comparison of the gene expression profiles between the HepG2 cell line and liver tissue revealed the highest variation, followed by HepG2 cells submitted to two different sample preparation protocols. The lowest variation was observed between HepG2 cells prepared by two different laboratories using the same protocol. The enrichment analysis of the genes that were differentially expressed between HepG2 cells and liver tissue mainly revealed the cancer-associated gene signature of HepG2 cells and the activation of the response to chemical stimuli in the liver tissue. The HepG2 transcriptome obtained with the SOLiD platform was highly correlated with the published transcriptome obtained with the Illumina and Helicos platforms, with moderate correspondence to microarrays. CONCLUSIONS: In the present study, we assessed the influence of experimental factors on the HepG2 transcriptome and identified differences in gene expression between the HepG2 cell line and liver cells. These findings will facilitate robust experimental design in the fields of pharmacology and toxicology. Our results were supported by a comparative analysis with previous HepG2 gene expression studies.

Chromosome 18 transcriptoproteome of liver tissue and HepG2 cells and targeted proteome mapping in depleted plasma: update 2013.

We report the results obtained in 2012-2013 by the Russian Consortium for the Chromosome-centric Human Proteome Project (C-HPP). The main scope of this work was the transcriptome profiling of genes on human chromosome 18 (Chr 18), as well as their encoded proteome, from three types of biomaterials: liver tissue, the hepatocellular carcinoma-derived cell line HepG2, and blood plasma. The transcriptome profiling for liver tissue was independently performed using two RNaseq platforms (SOLiD and Illumina) and also by droplet digital PCR (ddPCR) and quantitative RT-PCR. The proteome profiling of Chr 18 was accomplished by quantitatively measuring protein copy numbers in the three types of biomaterial (the lowest protein concentration measured was 10(-13) M) using selected reaction monitoring (SRM). In total, protein copy numbers were estimated for 228 master proteins, including quantitative data on 164 proteins in plasma, 171 in the HepG2 cell line, and 186 in liver tissue. Most proteins were present in plasma at 10(8) copies/µL, while the median abundance was 10(4) and 10(5) protein copies per cell in HepG2 cells and liver tissue, respectively. In summary, for liver tissue and HepG2 cells a "transcriptoproteome" was produced that reflects the relationship between transcript and protein copy numbers of the genes on Chr 18. The quantitative data acquired by RNaseq, PCR, and SRM were uploaded into the "Update_2013" data set of our knowledgebase (www.kb18.ru) and investigated for linear correlations.