The study of differential gene expression patterns through RNA-Seq comprises a routine task in the daily lives of molecular bioscientists, who produce vast amounts of data requiring proper management and analysis. Despite widespread use, there are still no widely accepted golden standards for the normalization and statistical analysis of RNA-Seq data, and critical biases, such as gene lengths and problems in the detection of certain types of molecules, remain largely unaddressed. Stimulated by these unmet needs and the lack of in-depth research into the potential of combinatorial methods to enhance the analysis of differential gene expression, we had previously introduced the PANDORA P-value combination algorithm while presenting evidence for PANDORA’s superior performance in optimizing the tradeoff between precision and sensitivity. In this article, we present the next generation of the algorithm along with a more in-depth investigation of its capabilities to effectively analyze RNA-Seq data. In particular, we show that PANDORA-reported lists of differentially expressed genes are unaffected by biases introduced by different normalization methods, while, at the same time, they comprise a reliable input option for downstream pathway analysis. Additionally, PANDORA outperforms other methods in detecting differential expression patterns in certain transcript types, including long non-coding RNAs.

Autotaxin (ATX; ENPP2) is a secreted lysophospholipase D catalyzing the extracellular production of lysophosphatidic acid (LPA), a pleiotropic signaling phospholipid. Genetic and pharmacologic studies have previously established a pathologic role for ATX and LPA signaling in pulmonary injury, inflammation, and fibrosis. Here, increased ENPP2 mRNA levels were detected in immune cells from nasopharyngeal swab samples of COVID-19 patients, and increased ATX serum levels were found in severe COVID-19 patients. ATX serum levels correlated with the corresponding increased serum levels of IL-6 and endothelial damage biomarkers, suggesting an interplay of the ATX/LPA axis with hyperinflammation and the associated vascular dysfunction in COVID-19. Accordingly, dexamethasone (Dex) treatment of mechanically ventilated patients reduced ATX levels, as shown in two independent cohorts, indicating that the therapeutic benefits of Dex include the suppression of ATX. Moreover, large scale analysis of multiple single cell RNA sequencing datasets revealed the expression landscape of ENPP2 in COVID-19 and further suggested a role for ATX in the homeostasis of dendritic cells, which exhibit both numerical and functional deficits in COVID-19. Therefore, ATX has likely a multifunctional role in COVID-19 pathogenesis, suggesting that its pharmacological targeting might represent an additional therapeutic option, both during and after hospitalization.

The activation and accumulation of lung fibroblasts resulting in aberrant deposition of extracellular matrix components, is a pathogenic hallmark of Idiopathic Pulmonary Fibrosis, a lethal and incurable disease. In this report, increased expression of TKS5, a scaffold protein essential for the formation of podosomes, was detected in the lung tissue of Idiopathic Pulmonary Fibrosis patients and bleomycin-treated mice. Τhe profibrotic milieu is found to induceTKS5expression and the formation of prominent podosome rosettes in lung fibroblasts, that are retained ex vivo, culminating in increased extracellular matrix invasion.Tks5^+/-mice are found resistant to bleomycin-induced pulmonary fibrosis, largely attributed to diminished podosome formation in fibroblasts and decreased extracellular matrix invasion. As computationally predicted, inhibition of src kinase is shown to potently attenuate podosome formation in lung fibroblasts and extracellular matrix invasion, and bleomycin-induced pulmonary fibrosis, suggesting pharmacological targeting of podosomes as a very promising therapeutic option in pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is characterized by exuberant deposition of extracellular matrix components, leading to the deterioration of lung architecture and respiratory functions. Profibrotic mechanisms are controlled by multiple regulatory molecules, including MAPKs, in turn regulated by multiple phosphorylation cascades. MAP3K8 is an MAPK kinase kinase suggested to pleiotropically regulate multiple pathogenic pathways in the context of inflammation and cancer; however, a possible role in the pathogenesis of IPF has not been investigated. In this report, MAP3K8 mRNA levels were found decreased in the lungs of IPF patients and of mice upon bleomycin-induced pulmonary fibrosis. Ubiquitous genetic deletion of Map3k8 in mice exacerbated the modeled disease, whereas bone marrow transfer experiments indicated that although MAP3K8 regulatory functions are active in both hematopoietic and nonhematopoietic cells, Map3k8 in hematopoietic cells has a more dominant role. Macrophage-specific deletion of Map3k8 was further found to be sufficient for disease exacerbation thus confirming a major role for macrophages in pulmonary fibrotic responses and suggesting a main role for Map3k8 in the homeostasis of their effector functions in the lung. Map3k8 deficiency was further shown to be associated with decreased Cox-2 expression, followed by a decrease in PGE₂ production in the lung; accordingly, exogenous administration of PGE₂ reduced inflammation and reversed the exacerbated fibrotic profile of Map3k8^−/− mice. Therefore, MAP3K8 has a central role in the regulation of inflammatory responses and Cox-2–mediated PGE₂ production in the lung, and the attenuation of its expression is integral to pulmonary fibrosis development. Key Points Pulmonary fibrosis entails a transcriptional downregulation of Map3k8 expression. Map3k8 suppress pulmonary fibrosis via Cox-2–mediated PGE₂ production.

Autotaxin (ATX), encoded by the ctonucleotide pyrophosphatase/phosphodiesterase 2 (ENPP2) gene, is a key enzyme in lysophosphatidic acid (LPA) synthesis. We have recently described ENPP2 methylation profiles in health and multiple malignancies and demonstrated correlation to its aberrant expression. Here we focus on breast cancer (BrCa), analyzing in silico publicly available BrCa methylome datasets, to identify differentially methylated CpGs (DMCs) and correlate them with expression. Numerous DMCs were identified between BrCa and healthy breast tissues in the gene body and promoter-associated regions (PA). PA DMCs were upregulated in BrCa tissues in relation to normal, in metastatic BrCa in relation to primary, and in stage I BrCa in relation to normal, and this was correlated to decreased mRNA expression. The first exon DMC was also investigated in circulating cell free DNA (ccfDNA) isolated by BrCa patients; methylation was increased in BrCa in relation to ccfDNA from healthy individuals, confirming in silico results. It also differed between patient groups and was correlated to the presence of multiple metastatic sites. Our data indicate that promoter methylation of ENPP2 arrests its transcription in BrCa and introduce first exon methylation as a putative biomarker for diagnosis and monitoring which can be assessed in liquid biopsy.

Autotaxin (ATX) encoded by Ectonucleotide Pyrophosphatase/Phosphodiesterase 2 (ENPP2) is a key enzyme in Lysophosphatidic Acid (LPA) synthesis implicated in cancer. Although its aberrant expression has been reported, ENPP2 methylation profiles in health and malignancy are not described. We examined in silico the methylation of ENPP2 analyzing publicly available methylome datasets, to identify Differentially Methylated CpGs (DMCs) which were then correlated with expression at gene and isoform levels. Significance indication was set to be FDR corrected p-value < 0.05. Healthy tissues presented methylation in all gene body CGs and lower levels in Promoter Associated (PA) regions, whereas in the majority of the tumors examined (HCC, melanoma, CRC, LC and PC) the methylation pattern was reversed. DMCs identified in the promoter were located in sites recognized by multiple transcription factors, suggesting involvement in gene expression. Alterations in methylation were correlated to an aggressive phenotype in cancer cell lines. In prostate and lung adenocarcinomas, increased methylation of PA CGs was correlated to decreased ENPP2 mRNA expression and to poor prognosis parameters. Collectively, our results corroborate that methylation is an active level of ATX expression regulation in cancer. Our study provides an extended description of the methylation status of ENPP2 in health and cancer and points out specific DMCs of value as prognostic biomarkers.