PREDA: an R-package to identify regional variations in genomic data

Definition of a Multi-Omics Signature for Esophageal Adenocarcinoma Prognosis Prediction

Esophageal cancer is a highly lethal malignancy, representing 5% of all cancer-related deaths. The two main subtypes are esophageal squamous cell carcinoma (ESCC) and esophageal adenocarcinoma (EAC). While most research has focused on ESCC, few studies have analyzed EAC for transcriptional signatures linked to diagnosis or prognosis. In this study, we utilized single-cell RNA sequencing and bulk RNA sequencing to identify specific immune cell types that contribute to anti-tumor responses, as well as differentially expressed genes (DEGs). We have characterized transcriptional signatures, validated against a wide cohort of TCGA patients, that are capable of predicting clinical outcomes and the prognosis of EAC post-surgery with efficacy comparable to the currently accepted prognostic factors. In conclusion, our findings provide insights into the immune landscape and therapeutic targets of EAC, proposing novel immunological biomarkers for predicting prognosis, aiding in patient stratification for post-surgical outcomes, follow-up, and personalized adjuvant therapy decisions.

Maximum-scoring path sets on pangenome graphs of constant treewidth

We generalize a problem of finding maximum-scoring segment sets, previously studied by Csűrös (IEEE/ACM Transactions on Computational Biology and Bioinformatics, 2004, 1, 139-150), from sequences to graphs. Namely, given a vertex-weighted graph G and a non-negative startup penalty c, we can find a set of vertex-disjoint paths in G with maximum total score when each path's score is its vertices' total weight minus c. We call this new problem maximum-scoring path sets (MSPS). We present an algorithm that has a linear-time complexity for graphs with a constant treewidth. Generalization from sequences to graphs allows the algorithm to be used on pangenome graphs representing several related genomes and can be seen as a common abstraction for several biological problems on pangenomes, including searching for CpG islands, ChIP-seq data analysis, analysis of region enrichment for functional elements, or simple chaining problems.

Sequencing-based transcriptome analysis reveals diversification of immune response- and angiogenesis-related expression patterns of early-stage cervical carcinoma as compared with high-grade CIN

Background

Molecular diversity of virus-associated cervical cancer remains a relatively underexplored issue, and interrelations of immunologic and angiogenic features during the establishment of a particular landscape of the cervical cancer microenvironment are not well-characterized, especially for its earliest clinical stages, although this may provide insight into the mechanisms behind the differences in tumor aggressiveness, treatment responsiveness and prognosis. In this research, we were aimed at identifying transcriptomic landscapes of early-stage cervical carcinoma that differ substantially in their immune-related characteristics, patterns of signaling pathways and composition of the microenvironment in comparison with immediate precursor (intraepithelial) lesions.

Methods

We performed the Illumina platform-based RNA sequencing using a panel of fresh tissue samples that included human papillomavirus-positive cervical intraepithelial neoplastic lesions (CIN), invasive squamous carcinoma of the cervix of FIGO IA1-IIB stages, and morphologically normal epithelium. The derived transcriptomic profiles were bioinformatically analyzed and compared by patterns of signaling pathway activation, distribution of tumor-infiltrating cell populations, and genomic regions involved.

Result

According to hierarchical cluster analysis of the whole-transcriptome profiles, tissue samples were distributed between three groups, or gene expression patterns (the one comprising most pre-cancer cases and the other two encompassing mostly early-stage invasive cancer cases). Differentially expressed genes were retrieved in each intergroup pairwise comparison followed by Gene Ontology analysis. Gene set enrichment analysis of the two groups of tumor samples in comparison with the CIN group identified substantial differences in immunological and angiogenic properties between tumorous groups suggesting the development of different molecular phenotypes. Cell composition analysis confirmed the diverse changes in the abundancies of immune and non-immune populations and, accordingly, different impacts of the immune and stromal compartments on the tumor microenvironment in these two groups of tumors compared to CIN. Positional gene expression analysis demonstrated that the identified transcriptomic differences were linked to different chromosomal regions and co-localized with particular gene families implicated in immune regulation, inflammation, cell differentiation, and tumor invasion.

Conclusions

Overall, detection of different transcriptomic patterns of invasive cervical carcinoma at its earliest stages supports the diverse impacts of immune response- and angiogenesis-related mechanisms on the onset of tumor invasion and progression. This may provide new options for broadening the applicability and increasing the efficiency of target anti-angiogenic and immune-based therapy of virus-associated cervical carcinoma.

3D microenvironment attenuates simulated microgravity-mediated changes in T cell transcriptome

Human space travel and exploration are of interest to both the industrial and scientific community. However, there are many adverse effects of spaceflight on human physiology. In particular, there is a lack of understanding of the extent to which microgravity affects the immune system. T cells, key players of the adaptive immune system and long-term immunity, are present not only in blood circulation but also reside within the tissue. As of yet, studies investigating the effects of microgravity on T cells are limited to peripheral blood or traditional 2D cell culture that recapitulates circulating blood. To better mimic interstitial tissue, 3D cell culture has been well established for physiologically and pathologically relevant models. In this work, we utilize 2D cell culture and 3D collagen matrices to gain an understanding of how simulated microgravity, using a random positioning machine, affects both circulating and tissue-resident T cells. T cells were studied in both resting and activated stages. We found that 3D cell culture attenuates the effects of simulated microgravity on the T cells transcriptome and nuclear irregularities compared to 2D cell culture. Interestingly, simulated microgravity appears to have less effect on activated T cells compared to those in the resting stage. Overall, our work provides novel insights into the effects of simulated microgravity on circulating and tissue-resident T cells which could provide benefits for the health of space travellers.

Duchenne muscular dystrophy gene expression is an independent prognostic marker for IDH mutant low-grade glioma

Alterations in the expression of the Duchenne muscular dystrophy (DMD) gene have been associated with the development, progression and survival outcomes of numerous cancers including tumours of the central nervous system. We undertook a detailed bioinformatic analysis of low-grade glioma (LGG) bulk RNAseq data to characterise the association between DMD expression and LGG survival outcomes. High DMD expression was significantly associated with poor survival in LGG with a difference in median overall survival between high and low DMD groups of over 7 years (P = < 0.0001). In a multivariate model, DMD expression remained significant (P = 0.02) and was an independent prognostic marker for LGG. The effect of DMD expression on overall survival was only apparent in isocitrate dehydrogenase (IDH) mutant cases where non-1p/19q co-deleted LGG patients could be further stratified into high/low DMD groups. Patients in the high DMD group had a median overall survival time almost halve that of the low DMD group. The expression of the individual DMD gene products Dp71, Dp71ab and Dp427m were also significantly associated with overall survival in LGG which have differential biological effects relevant to the pathogenesis of LGG. Differential gene expression and pathway analysis identifies dysregulated biological processes relating to ribosome biogenesis, synaptic signalling, neurodevelopment, morphogenesis and immune pathways. Genes spanning almost the entirety of chromosome 1p are upregulated in patients with high overall DMD, Dp71 and Dp427m expression which worsens survival outcomes for these patients. We confirmed dystrophin protein is variably expressed in LGG tumour tissue by immunohistochemistry and, overall, demonstrate that DMD expression has potential utility as an independent prognostic marker which can further stratify IDH mutant LGG to identify those at risk of poor survival. This knowledge may improve risk stratification and management of LGG.

Sex differences in skeletal muscle revealed through fiber type, capillarity, and transcriptomics profiling in mice

Skeletal muscle anatomy and physiology are sexually dimorphic but molecular underpinnings and muscle-specificity are not well-established. Variances in metabolic health, fitness level, sedentary behavior, genetics, and age make it difficult to discern inherent sex effects in humans. Therefore, mice under well-controlled conditions were used to determine female and male (n = 19/sex) skeletal muscle fiber type/size and capillarity in superficial and deep gastrocnemius (GA-s, GA-d), soleus (SOL), extensor digitorum longus (EDL), and plantaris (PLT), and transcriptome patterns were also determined (GA, SOL). Summed muscle weight strongly correlated with lean body mass (r2  = 0.67, p < 0.0001, both sexes). Other phenotypes were muscle-specific: e.g., capillarity (higher density, male GA-s), myofiber size (higher, male EDL), and fiber type (higher, lower type I and type II prevalences, respectively, in female SOL). There were broad differences in transcriptomics, with >6000 (GA) and >4000 (SOL) mRNAs differentially-expressed by sex; only a minority of these were shared across GA and SOL. Pathway analyses revealed differences in ribosome biology, transcription, and RNA processing. Curation of sexually dimorphic muscle transcripts shared in GA and SOL, and literature datasets from mice and humans, identified 11 genes that we propose are canonical to innate sex differences in muscle: Xist, Kdm6a, Grb10, Oas2, Rps4x (higher, females) and Ddx3y, Kdm5d, Irx3, Wwp1, Aldh1a1, Cd24a (higher, males). These genes and those with the highest "sex-biased" expression in our study do not contain estrogen-response elements (exception, Greb1), but a subset are proposed to be regulated through androgen response elements. We hypothesize that innate muscle sexual dimorphism in mice and humans is triggered and then maintained by classic X inactivation (Xist, females) and Y activation (Ddx3y, males), with coincident engagement of X encoded (Kdm6a) and Y encoded (Kdm5d) demethylase epigenetic regulators that are complemented by modulation at some regions of the genome that respond to androgen.

TRAIP modulates the IGFBP3/AKT pathway to enhance the invasion and proliferation of osteosarcoma by promoting KANK1 degradation

Osteosarcoma is one of the most common primary malignancies in bones and is characterized by high metastatic rates. Circulating tumor cells (CTCs) derived from solid tumors can give rise to metastatic lesions, increasing the risk of death in patients with cancer. Here, we used bioinformatics tools to compare the gene expression between CTCs and metastatic lesions in osteosarcoma to identify novel molecular mechanisms underlying osteosarcoma metastasis. We identified TRAIP as a key differentially expressed gene with prognostic significance in osteosarcoma. We demonstrated that TRAIP regulated the proliferation and invasion of osteosarcoma cells. In addition, we found that TRAIP promoted KANK1 polyubiquitination and subsequent degradation, downregulating IGFBP3 and activating the AKT pathway in osteosarcoma cells. These results support the critical role of the TRAIP/KANK1/IGFBP3/AKT signaling axis in osteosarcoma progression and suggest that TRAIP may represent a promising therapeutic target for osteosarcoma.

Analysis of MCM Proteins' Role as a Potential Target of Statins in Patients with Acute Type A Aortic Dissection through Bioinformatics

Acute aortic dissection is one of the most severe vascular diseases. The molecular mechanisms of aortic expansion and dissection are unclear. Clinical studies have found that statins play a protective role in aortic dissection development and therapy; however, the mechanism of statins' effects on the aorta is unknown. The Gene Expression Omnibus (GEO) dataset GSE52093, GSE2450and GSE8686 were analyzed, and genes expressed differentially between aortic dissection samples and normal samples were determined using the Networkanalyst and iDEP tools. Weight gene correlation network analysis (WGCNA), functional annotation, pathway enrichment analysis, and the analysis of the regional variations of genomic features were then performed. We found that the minichromosome maintenance proteins (MCMs), a family of proteins targeted by statins, were upregulated in dissected aortic wall tissues and play a central role in cell-cycle and mitosis regulation in aortic dissection patients. Our results indicate a potential molecular target and mechanism for statins' effects in patients with acute type A aortic dissection.

Oncogenic transcriptomic profile is sustained in the liver after the eradication of the hepatitis C virus

Hepatocellular carcinoma (HCC) developing after hepatitis C virus (HCV) eradication is a serious clinical concern. However, molecular basis for the hepatocarcinogenesis after sustained virologic response (SVR) remains unclear. In this study, we aimed to unveil the transcriptomic profile of post-SVR liver tissues and explore the molecules associated with post-SVR carcinogenesis. We analysed 90 RNA sequencing datasets, consisting of non-cancerous liver tissues including 20 post-SVR, 40 HCV-positive and 7 normal livers, along with Huh7 cell line specimens before and after HCV infection and eradication. Comparative analysis demonstrated that cell cycle- and mitochondrial function-associated pathways were altered only in HCV-positive non-cancerous liver tissues, whereas some cancer-related pathways were up-regulated in the non-cancerous liver tissues of both post-SVR and HCV-positive cases. The persistent up-regulation of carcinogenesis-associated gene clusters after viral clearance was reconfirmed through in vitro experiments, of which, CYR61, associated with liver fibrosis and carcinogenesis in several cancer types, was the top enriched gene and co-expressed with cell proliferation-associated gene modules. To evaluate whether this molecule could be a predictor of hepatocarcinogenesis after cure of HCV infection, we also examined 127 sera from independent HCV-positive cohorts treated with direct-acting antivirals (DAAs), including 60 post-SVR-HCC patients, and found that the elevated serum Cyr61 was significantly associated with early carcinogenesis after receiving DAA therapy. In conclusion, some oncogenic transcriptomic profiles are sustained in liver tissues after HCV eradication, which might be a molecular basis for the liver cancer development even after viral clearance. Among them, up-regulated CYR61 could be a possible biomarker for post-SVR-HCC.

WoPPER: Web server for Position Related data analysis of gene Expression in Prokaryotes

The structural and conformational organization of chromosomes is crucial for gene expression regulation in eukaryotes and prokaryotes as well. Up to date, gene expression data generated using either microarray or RNA-sequencing are available for many bacterial genomes. However, differential gene expression is usually investigated with methods considering each gene independently, thus not taking into account the physical localization of genes along a bacterial chromosome. Here, we present WoPPER, a web tool integrating gene expression and genomic annotations to identify differentially expressed chromosomal regions in bacteria. RNA-sequencing or microarray-based gene expression data are provided as input, along with gene annotations. The user can select genomic annotations from an internal database including 2780 bacterial strains, or provide custom genomic annotations. The analysis produces as output the lists of positionally related genes showing a coordinated trend of differential expression. Graphical representations, including a circular plot of the analyzed chromosome, allow intuitive browsing of the results. The analysis procedure is based on our previously published R-package PREDA. The release of this tool is timely and relevant for the scientific community, as WoPPER will fill an existing gap in prokaryotic gene expression data analysis and visualization tools. WoPPER is open to all users and can be reached at the following URL: https://WoPPER.ba.itb.cnr.it

iDEP: an integrated web application for differential expression and pathway analysis of RNA-Seq data

BACKGROUND

RNA-seq is widely used for transcriptomic profiling, but the bioinformatics analysis of resultant data can be time-consuming and challenging, especially for biologists. We aim to streamline the bioinformatic analyses of gene-level data by developing a user-friendly, interactive web application for exploratory data analysis, differential expression, and pathway analysis.

RESULTS

iDEP (integrated Differential Expression and Pathway analysis) seamlessly connects 63 R/Bioconductor packages, 2 web services, and comprehensive annotation and pathway databases for 220 plant and animal species. The workflow can be reproduced by downloading customized R code and related pathway files. As an example, we analyzed an RNA-Seq dataset of lung fibroblasts with Hoxa1 knockdown and revealed the possible roles of SP1 and E2F1 and their target genes, including microRNAs, in blocking G1/S transition. In another example, our analysis shows that in mouse B cells without functional p53, ionizing radiation activates the MYC pathway and its downstream genes involved in cell proliferation, ribosome biogenesis, and non-coding RNA metabolism. In wildtype B cells, radiation induces p53-mediated apoptosis and DNA repair while suppressing the target genes of MYC and E2F1, and leads to growth and cell cycle arrest. iDEP helps unveil the multifaceted functions of p53 and the possible involvement of several microRNAs such as miR-92a, miR-504, and miR-30a. In both examples, we validated known molecular pathways and generated novel, testable hypotheses.

CONCLUSIONS

Combining comprehensive analytic functionalities with massive annotation databases, iDEP ( http://ge-lab.org/idep/ ) enables biologists to easily translate transcriptomic and proteomic data into actionable insights.

Chromosome positioning in interphase nuclei of hematopoietic stem cell and myeloid precursor

Human myelopoiesis is an intriguing biological process during which multipotent stem cells limit their differentiation potential generating precursors that evolve into terminally differentiated cells. The differentiation process is correlated with differential gene expression and changes in nuclear architecture. In interphase, chromosomes are distinct entities known as chromosome territories and they show a radial localization that could result in a constrain of inter-homologous distance. This element plays a role in genome stability and gene expression. Here, we provide the first experimental evidence of 3D chromosomal arrangement considering two steps of human normal myelopoiesis. Specifically, multicolor 3D-FISH and 3D image analysis revealed that, in both normal human hematopoietic stem cells and myelod precursors CD14-, chromosomal position is correlated with gene density. However, we observed that inter-homologue distances are totally different during differentiation. This could be associated with differential gene expression that we found comparing the two cell types. Our results disclose an unprecedented framework relevant for deciphering the genomic mechanisms at the base of normal human myelopoiesis.

Integrative analysis of copy number and gene expression data suggests novel pathogenetic mechanisms in primary myelofibrosis

Primary myelofibrosis (PMF) is a Myeloproliferative Neoplasm (MPN) characterized by megakaryocyte hyperplasia, progressive bone marrow fibrosis, extramedullary hematopoiesis and transformation to Acute Myeloid Leukemia (AML). A number of phenotypic driver (JAK2, CALR, MPL) and additional subclonal mutations have been described in PMF, pointing to a complex genomic landscape. To discover novel genomic lesions that can contribute to disease phenotype and/or development, gene expression and copy number signals were integrated and several genomic abnormalities leading to a concordant alteration in gene expression levels were identified. In particular, copy number gain in the polyamine oxidase (PAOX) gene locus was accompanied by a coordinated transcriptional up-regulation in PMF patients. PAOX inhibition resulted in rapid cell death of PMF progenitor cells, while sparing normal cells, suggesting that PAOX inhibition could represent a therapeutic strategy to selectively target PMF cells without affecting normal hematopoietic cells' survival. Moreover, copy number loss in the chromatin modifier HMGXB4 gene correlates with a concomitant transcriptional down-regulation in PMF patients. Interestingly, silencing of HMGXB4 induces megakaryocyte differentiation, while inhibiting erythroid development, in human hematopoietic stem/progenitor cells. These results highlight a previously un-reported, yet potentially interesting role of HMGXB4 in the hematopoietic system and suggest that genomic and transcriptional imbalances of HMGXB4 could contribute to the aberrant expansion of the megakaryocytic lineage that characterizes PMF patients.

Prospective Biomarker Analysis of the Randomized CHER-LOB Study Evaluating the Dual Anti-HER2 Treatment With Trastuzumab and Lapatinib Plus Chemotherapy as Neoadjuvant Therapy for HER2-Positive Breast Cancer

BACKGROUND

The CHER-LOB randomized phase II study showed that the combination of lapatinib and trastuzumab plus chemotherapy increases the pathologic complete remission (pCR) rate compared with chemotherapy plus either trastuzumab or lapatinib. A biomarker program was prospectively planned to identify potential predictors of sensitivity to different treatments and to evaluate treatment effect on tumor biomarkers.

MATERIALS AND METHODS

Overall, 121 breast cancer patients positive for human epidermal growth factor 2 (HER2) were randomly assigned to neoadjuvant chemotherapy plus trastuzumab, lapatinib, or both trastuzumab and lapatinib. Pre- and post-treatment samples were centrally evaluated for HER2, p95-HER2, phosphorylated AKT (pAKT), phosphatase and tensin homolog, Ki67, apoptosis, and PIK3CA mutations. Fresh-frozen tissue samples were collected for genomic analyses.

RESULTS

A mutation in PIK3CA exon 20 or 9 was documented in 20% of cases. Overall, the pCR rates were similar in PIK3CA wild-type and PIK3CA-mutated patients (33.3% vs. 22.7%; p = .323). For patients receiving trastuzumab plus lapatinib, the probability of pCR was higher in PIK3CA wild-type tumors (48.4% vs. 12.5%; p = .06). Ki67, pAKT, and apoptosis measured on the residual disease were significantly reduced from baseline. The degree of Ki67 inhibition was significantly higher in patients receiving the dual anti-HER2 blockade. The integrated analysis of gene expression and copy number data demonstrated that a 50-gene signature specifically predicted the lapatinib-induced pCR.

CONCLUSION

PIK3CA mutations seem to identify patients who are less likely to benefit from dual anti-HER2 inhibition. p95-HER2 and markers of phosphoinositide 3-kinase pathway deregulation are not confirmed as markers of different sensitivity to trastuzumab or lapatinib.

IMPLICATIONS FOR PRACTICE

HER2 is currently the only validated marker to select breast cancer patients for anti-HER2 treatment; however, it is becoming evident that HER2-positive breast cancer is a heterogeneous disease. In addition, more and more new anti-HER2 treatments are becoming available. There is a need to identify markers of sensitivity to different treatments to move in the direction of treatment personalization. This study identified PIK3CA mutations as a potential predictive marker of resistance to dual anti-HER2 treatment that should be further studied in breast cancer.

Cancer gene prioritization by integrative analysis of mRNA expression and DNA copy number data: a comparative review

A variety of genome-wide profiling techniques are available to investigate complementary aspects of genome structure and function. Integrative analysis of heterogeneous data sources can reveal higher level interactions that cannot be detected based on individual observations. A standard integration task in cancer studies is to identify altered genomic regions that induce changes in the expression of the associated genes based on joint analysis of genome-wide gene expression and copy number profiling measurements. In this review, we highlight common approaches to genomic data integration and provide a transparent benchmarking procedure to quantitatively compare method performances in cancer gene prioritization. Algorithms, data sets and benchmarking results are available at http://intcomp.r-forge.r-project.org.

Comparative genomics and transcriptional profiles of Saccharopolyspora erythraea NRRL 2338 and a classically improved erythromycin over-producing strain

Background

The molecular mechanisms altered by the traditional mutation and screening approach during the improvement of antibiotic-producing microorganisms are still poorly understood although this information is essential to design rational strategies for industrial strain improvement. In this study, we applied comparative genomics to identify all genetic changes occurring during the development of an erythromycin overproducer obtained using the traditional mutate-and- screen method.

Results

Compared with the parental Saccharopolyspora erythraea NRRL 2338, the genome of the overproducing strain presents 117 deletion, 78 insertion and 12 transposition sites, with 71 insertion/deletion sites mapping within coding sequences (CDSs) and generating frame-shift mutations. Single nucleotide variations are present in 144 CDSs. Overall, the genomic variations affect 227 proteins of the overproducing strain and a considerable number of mutations alter genes of key enzymes in the central carbon and nitrogen metabolism and in the biosynthesis of secondary metabolites, resulting in the redirection of common precursors toward erythromycin biosynthesis. Interestingly, several mutations inactivate genes coding for proteins that play fundamental roles in basic transcription and translation machineries including the transcription anti-termination factor NusB and the transcription elongation factor Efp. These mutations, along with those affecting genes coding for pleiotropic or pathway-specific regulators, affect global expression profile as demonstrated by a comparative analysis of the parental and overproducer expression profiles. Genomic data, finally, suggest that the mutate-and-screen process might have been accelerated by mutations in DNA repair genes.

Conclusions

This study helps to clarify the mechanisms underlying antibiotic overproduction providing valuable information about new possible molecular targets for rationale strain improvement.