Oxford Nanopore MinION Direct RNA-Seq for Systems Biology

The Translatome Map: RNC-Seq vs. Ribo-Seq for Profiling of HBE, A549, and MCF-7 Cell Lines

Gene expression is a tightly regulated process that involves multiple layers of control, including transcriptional, post-transcriptional, and translational regulation. To gain a comprehensive understanding of gene expression dynamics and its functional implications, it is crucial to compare translatomic, transcriptomic, and proteomic data. The two most common analysis methods, Ribo-seq and RNC-Seq, were used to analyze the translatome of the same sample, whose datasets were downloaded from the TranslatomeDB database. The resulting translatome maps obtained for three cell lines (HBE, A549, and MCF-7) using these two methods were comparatively analyzed. The two methods of translatome analysis were shown to provide comparable results and can be used interchangeably. The obtained mRNA translation patterns were annotated in the transcriptome and proteome context for the same sample, which may become the basis for the reconstruction of the molecular mechanisms of pathological process development in the future.

Transcriptome-wide N6-methyladenosine profiling reveals growth-defense trade-offs in the response of rice to brown planthopper (Nilaparvata lugens) infestation

BACKGROUND

N6-Methyladenosine (m6A) is a common messenger RNA (mRNA) modification that affects various physiological processes in stress responses. However, the role of m6A modifications in plants responses to herbivore stress remains unclear.

RESULTS

Here, we found that an infestation of brown planthopper (Nilaparvata lugens) female adults enhanced the resistance of rice to N. lugens. The m6A methylome analysis of N. lugens-infested and uninfested rice samples was performed to explore the interaction between rice and N. lugens. The m6A methylation mainly occurred in genes that were actively expressed in rice following N. lugens infestation, while an analysis of the whole-genomic mRNA distribution of m6A showed that N. lugens infestation caused an overall decrease in the number of m6A methylation sites across the chromosomes. The m6A methylation of genes involved in the m6A modification machinery and several defense-related phytohormones (jasmonic acid and salicylic acid) pathways was increased in N. lugens-infested rice compared to that in uninfested rice. In contrast, m6A modification levels of growth-related phytohormone (auxin and gibberellin) biosynthesis-related genes were significantly attenuated during N. lugens infestation, accompanied by the down-regulated expression of these transcripts, indicating that rice growth was restricted during N. lugens attack to rapidly optimize resource allocation for plant defense. Integrative analysis of the differential patterns of m6A methylation and the corresponding transcripts showed a positive correlation between m6A methylation and transcriptional regulation.

CONCLUSION

The m6A modification is an important strategy for regulating the expression of genes involved in rice defense and growth during rice-N. lugens interactions. These findings provide new ideas for formulating strategies to control herbivorous pests. © 2024 Society of Chemical Industry.

Exploring in vitro modeling in hepatocarcinogenesis research: morphological and molecular features and similarities to the corresponding human disease

The hepatocellular carcinoma (HCC) features a remarkable epidemiological burden, ranking as the third most lethal cancer worldwide. As the HCC-related molecular and cellular complexity unfolds as the disease progresses, the use of a myriad of in vitro models available is mandatory in translational preclinical research setups. In this review paper, we will compile cutting-edge information on the in vitro bioassays for HCC research, (A) emphasizing their morphological and molecular parallels with human HCC; (B) delineating the advantages and limitations of their application; and (C) offering perspectives on their prospective applications. While bidimensional (2D) (co) culture setups provide a rapid low-cost strategy for metabolism and drug screening investigations, tridimensional (3D) (co) culture bioassays - including patient-derived protocols as organoids and precision cut slices - surpass some of the 2D strategies limitations, mimicking the complex microarchitecture and cellular and non-cellular microenvironment observed in human HCC. 3D models have become invaluable tools to unveil HCC pathophysiology and targeted therapy. In both setups, the recapitulation of HCC in different etiologies/backgrounds (i.e., viral, fibrosis, and fatty liver) may be considered as a fundamental guide for obtaining translational findings. Therefore, a "multimodel" approach - encompassing the advantages of different in vitro bioassays - is encouraged to circumvent "model-biased" outcomes in preclinical HCC research.

Homoeologs in Allopolyploids: Navigating Redundancy as Both an Evolutionary Opportunity and a Technical Challenge-A Transcriptomics Perspective

Allopolyploidy in plants involves the merging of two or more distinct parental genomes into a single nucleus, a significant evolutionary process in the plant kingdom. Transcriptomic analysis provides invaluable insights into allopolyploid plants by elucidating the fate of duplicated genes, revealing evolutionary novelties and uncovering their environmental adaptations. By examining gene expression profiles, scientists can discern how duplicated genes have evolved to acquire new functions or regulatory roles. This process often leads to the development of novel traits and adaptive strategies that allopolyploid plants leverage to thrive in diverse ecological niches. Understanding these molecular mechanisms not only enhances our appreciation of the genetic complexity underlying allopolyploidy but also underscores their importance in agriculture and ecosystem resilience. However, transcriptome profiling is challenging due to genomic redundancy, which is further complicated by the presence of multiple chromosomes sets and the variations among homoeologs and allelic genes. Prior to transcriptome analysis, sub-genome phasing and homoeology inference are essential for obtaining a comprehensive view of gene expression. This review aims to clarify the terminology in this field, identify the most challenging aspects of transcriptome analysis, explain their inherent difficulties, and suggest reliable analytic strategies. Furthermore, bulk RNA-seq is highlighted as a primary method for studying allopolyploid gene expression, focusing on critical steps like read mapping and normalization in differential gene expression analysis. This approach effectively captures gene expression from both parental genomes, facilitating a comprehensive analysis of their combined profiles. Its sensitivity in detecting low-abundance transcripts allows for subtle differences between parental genomes to be identified, crucial for understanding regulatory dynamics and gene expression balance in allopolyploids.

The dynamics of N6-methyladenine RNA modification in resistant and susceptible rice varieties responding to rice stem borer damage

N6-methyladenosine (m6A) is the most prevalent modification in cellular RNA which orchestrates diverse physiological and pathological processes during stress response. However, the differential m6A modifications that cope with herbivore stress in resistant and susceptible crop varieties remain unclear. Here, we found that rice stem borer (RSB) larvae grew better on indica rice (e.g., MH63, IR64, Nanjing 11) than on japonica rice varieties (e.g., Nipponbare, Zhonghua 11, Xiushui 11). Then, transcriptome-wide m6A profiling of representative resistant (Nipponbare) and susceptible (MH63) rice varieties were performed using a nanopore direct RNA sequencing approach, to reveal variety-specific m6A modifications against RSB. Upon RSB infestation, m6A methylation occurred in actively expressed genes in Nipponbare and MH63, but the number of methylation sites decreased across rice chromosomes. Integrative analysis showed that m6A methylation levels were closely associated with transcriptional regulation. Genes involved in herbivorous resistance related to mitogen-activated protein kinase, jasmonic acid (JA), and terpenoid biosynthesis pathways, as well as JA-mediated trypsin protease inhibitors, were heavily methylated by m6A, and their expression was more pronounced in RSB-infested Nipponbare than in RSB-infested MH63, which may have contributed to RSB resistance in Nipponbare. Therefore, dynamics of m6A modifications act as the main regulatory strategy for expression of genes involved in plant-insect interactions, which is attributed to differential responses of resistant and susceptible rice varieties to RSB infestation. These findings could contribute to developing molecular breeding strategies for controlling herbivorous pests.

Non-Targeted RNA Sequencing: Towards the Development of Universal Clinical Diagnosis Methods for Human and Veterinary Infectious Diseases

Metagenomics offers the potential to replace and simplify classical methods used in the clinical diagnosis of human and veterinary infectious diseases. Metagenomics boasts a high pathogen discovery rate and high specificity, advantages absent in most classical approaches. However, its widespread adoption in clinical settings is still pending, with a slow transition from research to routine use. While longer turnaround times and higher costs were once concerns, these issues are currently being addressed by automation, better chemistries, improved sequencing platforms, better databases, and automated bioinformatics analysis. However, many technical options and steps, each producing highly variable outcomes, have reduced the technology's operational value, discouraging its implementation in diagnostic labs. We present a case for utilizing non-targeted RNA sequencing (NT-RNA-seq) as an ideal metagenomics method for the detection of infectious disease-causing agents in humans and animals. Additionally, to create operational value, we propose to identify best practices for the "core" of steps that are invariably shared among many human and veterinary protocols. Reference materials, sequencing procedures, and bioinformatics standards should accelerate the validation processes necessary for the widespread adoption of this technology. Best practices could be determined through "implementation research" by a consortium of interested institutions working on common samples.

Bulk RNA-seq Assessment of Murine Spleen Using a Portable MinION Sequencing Device

Comprehensive profiling of the transcriptome reflects the dynamic changes in complex regulator interactions between genes and proteins in the pathological processes of autoimmune disease. To enable real-time adaptive sampling and run rapid sequencing with maximal expected readings of 50 Gb of data with accessibility to clinicians, basic researchers, and even students, a portable and affordable sequencing device, MinION, was employed in our laboratory for both basic and clinical studies. Here, the workflow of bulk RNA-seq in murine spleen using MinION is introduced. The methodology of both laboratory library preparation and the establishment of a bioinformatic pipeline are included in this chapter.

Identifying N6-Methyladenosine Sites in HepG2 Cell Lines Using Oxford Nanopore Technology

RNA modifications, particularly N6-methyladenosine (m6A), are pivotal regulators of RNA functionality and cellular processes. We analyzed m6A modifications by employing Oxford Nanopore technology and the m6Anet algorithm, focusing on the HepG2 cell line. We identified 3968 potential m6A modification sites in 2851 transcripts, corresponding to 1396 genes. A gene functional analysis revealed the active involvement of m6A-modified genes in ubiquitination, transcription regulation, and protein folding processes, aligning with the known role of m6A modifications in histone ubiquitination in cancer. To ensure data robustness, we assessed reproducibility across technical replicates. This study underscores the importance of evaluating algorithmic reproducibility, especially in supervised learning. Furthermore, we examined correlations between transcriptomic, translatomic, and proteomic levels. A strong transcriptomic-translatomic correlation was observed. In conclusion, our study deepens our understanding of m6A modifications' multifaceted impacts on cellular processes and underscores the importance of addressing reproducibility concerns in analytical approaches.

Special Issue on Differential Gene Expression and Coexpression

The most common approach in transcriptomics (RNA-seq and microarrays) is differential gene expression analysis (DGEA) [...].

Comparative Transcriptomic Analysis of Three Common Liver Cell Lines

BACKGROUND

Comparative transcriptomic analysis is a powerful approach for investigating the molecular mechanisms underlying various physiological and pathological processes, including liver disease. The liver is a vital organ with diverse functions, including metabolism and detoxification. In vitro models of liver cells, such as HepG2, Huh7, and Hep3B, have been widely used to study liver biology and pathology. However, there is limited information on the heterogeneity of these cell lines at the transcriptomic level.

OBJECTIVE

This study aimed to conduct a comparative transcriptomic analysis of three common liver cell lines (HepG2, Huh7, and Hep3B) using publicly available RNA-sequencing data. In addition, we compared these cell lines to primary hepatocytes, cells isolated directly from liver tissue and considered the gold standard for studying liver function and disease.

METHODS

Our study included sequencing data with the following criteria: total number of reads over 20,000,000, average read length of over 60 base pairs, Illumina sequencing, and non-treated cells. The data for the three cell lines were compiled: HepG2 (97 samples), Huh7 (39 samples), and Hep3B (16 samples). We performed differential gene expression analysis using the DESeq2 package, principal component analysis, hierarchical clustering on principal components, and correlation analysis to explore the heterogeneity within each cell line.

RESULTS

We identified numerous genes and pathways differentially expressed between HepG2, Huh7, and Hep3B, such as oxidative phosphorylation, cholesterol metabolism, and DNA damage. We report that the expression levels of important genes differ significantly between primary hepatocytes and liver cell lines.

CONCLUSION

Our study provides new insights into the transcriptional heterogeneity of commonly used liver cell lines and highlights the importance of considering specific cell line. Consequently, transferring results without considering the heterogeneity of cell lines is impractical and may lead to inaccurate or distorted conclusions.

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.

The Transition from Cancer "omics" to "epi-omics" through Next- and Third-Generation Sequencing

Deciphering cancer etiopathogenesis has proven to be an especially challenging task since the mechanisms that drive tumor development and progression are far from simple. An astonishing amount of research has revealed a wide spectrum of defects, including genomic abnormalities, epigenomic alterations, disturbance of gene transcription, as well as post-translational protein modifications, which cooperatively promote carcinogenesis. These findings suggest that the adoption of a multidimensional approach can provide a much more precise and comprehensive picture of the tumor landscape, hence serving as a powerful tool in cancer research and precision oncology. The introduction of next- and third-generation sequencing technologies paved the way for the decoding of genetic information and the elucidation of cancer-related cellular compounds and mechanisms. In the present review, we discuss the current and emerging applications of both generations of sequencing technologies, also referred to as massive parallel sequencing (MPS), in the fields of cancer genomics, transcriptomics and proteomics, as well as in the progressing realms of epi-omics. Finally, we provide a brief insight into the expanding scope of sequencing applications in personalized cancer medicine and pharmacogenomics.

A cohesin-associated gene score may predict immune checkpoint blockade in hepatocellular carcinoma

Stromal antigen 1 (STAG1), a component of cohesion, is overexpressed in various cancers, but it is unclear whether it has a role in the transcriptional regulation of hepatocellular carcinoma (HCC). To test this hypothesis, here, we screened global HCC datasets and performed multiscale embedded gene co-expression network analysis to identify the potential functional modules of differentially expressed STAG1 co-expressed genes. The putative transcriptional targets of STAG1 were identified using chromatin immunoprecipitation followed by high-throughput DNA sequencing. The cohesin-associated gene score (CAGS) was quantified using the The Cancer Genome Atlas HCC cohort and single-sample gene set enrichment analysis. Distinct cohesin-associated gene patterns were identified by calculating the euclidean distance of each patient. We assessed the potential ability of the CAGS in predicting immune checkpoint blockade (ICB) treatment response using IMvigor210 and GSE78220 cohorts. STAG1 was upregulated in 3313 HCC tissue samples compared with 2692 normal liver tissue samples (standard mean difference = 0.54). A total of three cohesin-associated gene patterns were identified, where cluster 2 had a high TP53 mutated rate and a poor survival outcome. Low CAGS predicted a significant survival advantage but presaged poor immunotherapy response. Differentially expressed STAG1 co-expression genes were enriched in the mitotic cell cycle, lymphocyte activation, and blood vessel development. PDS5A and PDGFRA were predicted as the downstream transcriptional targets of STAG1. In summary, STAG1 is significantly upregulated in global HCC tissue samples and may participate in blood vessel development and the mitotic cell cycle. A cohesin-associated gene scoring system may have potential to predict the ICB response.

Improved transcriptome assembly using a hybrid of long and short reads with StringTie

Short-read RNA sequencing and long-read RNA sequencing each have their strengths and weaknesses for transcriptome assembly. While short reads are highly accurate, they are rarely able to span multiple exons. Long-read technology can capture full-length transcripts, but its relatively high error rate often leads to mis-identified splice sites. Here we present a new release of StringTie that performs hybrid-read assembly. By taking advantage of the strengths of both long and short reads, hybrid-read assembly with StringTie is more accurate than long-read only or short-read only assembly, and on some datasets it can more than double the number of correctly assembled transcripts, while obtaining substantially higher precision than the long-read data assembly alone. Here we demonstrate the improved accuracy on simulated data and real data from Arabidopsis thaliana, Mus musculus, and human. We also show that hybrid-read assembly is more accurate than correcting long reads prior to assembly while also being substantially faster. StringTie is freely available as open source software at https://github.com/gpertea/stringtie.

[Splice variants of mRNA of cytochrome P450 genes: analysis by the nanopore sequencing method in human liver tissue and HepG2 cell line]

The analysis of cytochrome P450 transcripts was carried out by the nanopore sequencing in liver tissue samples of three donors and HepG2 line cells. It has been demonstrated that direct mRNA sequencing with a MinION nanopore sequencer (Oxford Nanopore Technologies) allows one to obtained quantitative profiles for transcripts (and their splice variants) of cytochrome P450 superfamily genes encoding isoforms involved in metabolism of the large (~80%) part of drugs. The splice variant profiles substantially differ for donors. The cytochrome P450 gene expression at the transcript level is significantly weaker in cells of the HepG2 line compared with that in the normal liver tissue. This limits the capability of the direct mRNA nanopore sequencing for studying alternative splicing of cytochrome P450 transcripts in HepG2 cells. Both quantitative and qualitative profiles of the cytochrome P450 gene expression at the transcript level are notably differ in human liver tissue and HepG2 cells.

The Curious Case of the HepG2 Cell Line: 40 Years of Expertise

Liver cancer is the third leading cause of cancer death worldwide. Representing such a dramatic impact on our lives, liver cancer is a significant public health concern. Sustainable and reliable methods for preventing and treating liver cancer require fundamental research on its molecular mechanisms. Cell lines are treated as in vitro equivalents of tumor tissues, making them a must-have for basic research on the nature of cancer. According to recent discoveries, certified cell lines retain most genetic properties of the original tumor and mimic its microenvironment. On the other hand, modern technologies allowing the deepest level of detail in omics landscapes have shown significant differences even between samples of the same cell line due to cross- and mycoplasma infection. This and other observations suggest that, in some cases, cell cultures are not suitable as cancer models, with limited predictive value for the effectiveness of new treatments. HepG2 is a popular hepatic cell line. It is used in a wide range of studies, from the oncogenesis to the cytotoxicity of substances on the liver. In this regard, we set out to collect up-to-date information on the HepG2 cell line to assess whether the level of heterogeneity of the cell line allows in vitro biomedical studies as a model with guaranteed production and quality.

The Curious Case of the HepG2 Cell Line: 40 Years of Expertise

Liver cancer is the third leading cause of cancer death worldwide. Representing such a dramatic impact on our lives, liver cancer is a significant public health concern. Sustainable and reliable methods for preventing and treating liver cancer require fundamental research on its molecular mechanisms. Cell lines are treated as in vitro equivalents of tumor tissues, making them a must-have for basic research on the nature of cancer. According to recent discoveries, certified cell lines retain most genetic properties of the original tumor and mimic its microenvironment. On the other hand, modern technologies allowing the deepest level of detail in omics landscapes have shown significant differences even between samples of the same cell line due to cross- and mycoplasma infection. This and other observations suggest that, in some cases, cell cultures are not suitable as cancer models, with limited predictive value for the effectiveness of new treatments. HepG2 is a popular hepatic cell line. It is used in a wide range of studies, from the oncogenesis to the cytotoxicity of substances on the liver. In this regard, we set out to collect up-to-date information on the HepG2 cell line to assess whether the level of heterogeneity of the cell line allows in vitro biomedical studies as a model with guaranteed production and quality.

The Curious Case of the HepG2 Cell Line: 40 Years of Expertise

Liver cancer is the third leading cause of cancer death worldwide. Representing such a dramatic impact on our lives, liver cancer is a significant public health concern. Sustainable and reliable methods for preventing and treating liver cancer require fundamental research on its molecular mechanisms. Cell lines are treated as in vitro equivalents of tumor tissues, making them a must-have for basic research on the nature of cancer. According to recent discoveries, certified cell lines retain most genetic properties of the original tumor and mimic its microenvironment. On the other hand, modern technologies allowing the deepest level of detail in omics landscapes have shown significant differences even between samples of the same cell line due to cross- and mycoplasma infection. This and other observations suggest that, in some cases, cell cultures are not suitable as cancer models, with limited predictive value for the effectiveness of new treatments. HepG2 is a popular hepatic cell line. It is used in a wide range of studies, from the oncogenesis to the cytotoxicity of substances on the liver. In this regard, we set out to collect up-to-date information on the HepG2 cell line to assess whether the level of heterogeneity of the cell line allows in vitro biomedical studies as a model with guaranteed production and quality.