Docker4Circ: A Framework for the Reproducible Characterization of circRNAs from RNA-Seq Data.

Recent improvements in cost-effectiveness of high-throughput technologies has allowed RNA sequencing of total transcriptomes suitable for evaluating the expression and regulation of circRNAs, a relatively novel class of transcript isoforms with suggested roles in transcriptional and post-transcriptional gene expression regulation, as well as their possible use as biomarkers, due to their deregulation in various human diseases. A limited number of integrated workflows exists for prediction, characterization, and differential expression analysis of circRNAs, none of them complying with computational reproducibility requirements. We developed Docker4Circ for the complete analysis of circRNAs from RNA-Seq data. Docker4Circ runs a comprehensive analysis of circRNAs in human and model organisms, including: circRNAs prediction; classification and annotation using six public databases; back-splice sequence reconstruction; internal alternative splicing of circularizing exons; alignment-free circRNAs quantification from RNA-Seq reads; and differential expression analysis. Docker4Circ makes circRNAs analysis easier and more accessible thanks to: (i) its R interface; (ii) encapsulation of computational tasks into docker images; (iii) user-friendly Java GUI Interface availability; and (iv) no need of advanced bash scripting skills for correct use. Furthermore, Docker4Circ ensures a reproducible analysis since all its tasks are embedded into a docker image following the guidelines provided by Reproducible Bioinformatics Project.

Laser capture microdissection for transcriptomic profiles in human skin biopsies.

BACKGROUND: The acquisition of reliable tissue-specific RNA sequencing data from human skin biopsy represents a major advance in research. However, the complexity of the process of isolation of specific layers from fresh-frozen human specimen by laser capture microdissection, the abundant presence of skin nucleases and RNA instability remain relevant methodological challenges. We developed and optimized a protocol to extract RNA from layers of human skin biopsies and to provide satisfactory quality and amount of mRNA sequencing data. RESULTS: The protocol includes steps of collection, embedding, freezing, histological coloration and relative optimization to preserve RNA extracted from specific components of fresh-frozen human skin biopsy of 14 subjects. Optimization of the protocol includes a preservation step in RNALater® Solution, the control of specimen temperature, the use of RNase Inhibitors and the time reduction of the staining procedure. The quality of extracted RNA was measured using the percentage of fragments longer than 200 nucleotides (DV200), a more suitable measurement for successful library preparation than the RNA Integrity Number (RIN). RNA was then enriched using the TruSeq® RNA Access Library Prep Kit (Illumina®) and sequenced on HiSeq® 2500 platform (Illumina®). Quality control on RNA sequencing data was adequate to get reliable data for downstream analysis. CONCLUSIONS: The described implemented and optimized protocol can be used for generating transcriptomics data on skin tissues, and it is potentially applicable to other tissues. It can be extended to multicenter studies, due to the introduction of an initial step of preservation of the specimen that allowed the shipment of biological samples.

Are oncoantigens suitable targets for anti-tumour therapy?

When a vaccine-elicited immune response is directed against oncoantigens--proteins required for the neoplastic process--the chance that the tumour will evade the vaccine should be reduced. But how can these causal oncoantigens be identified? One approach is to find tumour-associated and microenvironment-associated oncoantigens required for progression from one tumour stage to the next by comparing gene signatures isolated from the different stages of tumour progression in cancer-prone transgenic mice. Mouse oncoantigens subsequently shown to be involved in human cancer can then be validated in mouse vaccination experiments. This provides the groundwork for the rational design of cancer vaccines for clinical trials.

Alternative splicing detection workflow needs a careful combination of sample prep and bioinformatics analysis.

BACKGROUND: RNA-Seq provides remarkable power in the area of biomarkers discovery and disease characterization. Two crucial steps that affect RNA-Seq experiment results are Library Sample Preparation (LSP) and Bioinformatics Analysis (BA). This work describes an evaluation of the combined effect of LSP methods and BA tools in the detection of splice variants. RESULTS: Different LSPs (TruSeq unstranded/stranded, ScriptSeq, NuGEN) allowed the detection of a large common set of splice variants. However, each LSP also detected a small set of unique transcripts that are characterized by a low coverage and/or FPKM. This effect was particularly evident using the low input RNA NuGEN v2 protocol. A benchmark dataset, in which synthetic reads as well as reads generated from standard (Illumina TruSeq 100) and low input (NuGEN) LSPs were spiked-in was used to evaluate the effect of LSP on the statistical detection of alternative splicing events (AltDE). Statistical detection of AltDE was done using as prototypes for splice variant-quantification Cuffdiff2 and RSEM-EBSeq. As prototype for exon-level analysis DEXSeq was used. Exon-level analysis performed slightly better than splice variant-quantification approaches, although at most only 50% of the spiked-in transcripts was detected. The performances of both splice variant-quantification and exon-level analysis improved when raising the number of input reads. CONCLUSION: Data, derived from NuGEN v2, were not the ideal input for AltDE, especially when the exon-level approach was used. We observed that both splice variant-quantification and exon-level analysis performances were strongly dependent on the number of input reads. Moreover, the ribosomal RNA depletion protocol was less sensitive in detecting splicing variants, possibly due to the significant percentage of the reads mapping to non-coding transcripts.

Prdm5 regulates collagen gene transcription by association with RNA polymerase II in developing bone.

PRDM family members are transcriptional regulators involved in tissue specific differentiation. PRDM5 has been reported to predominantly repress transcription, but a characterization of its molecular functions in a relevant biological context is lacking. We demonstrate here that Prdm5 is highly expressed in developing bones; and, by genome-wide mapping of Prdm5 occupancy in pre-osteoblastic cells, we uncover a novel and unique role for Prdm5 in targeting all mouse collagen genes as well as several SLRP proteoglycan genes. In particular, we show that Prdm5 controls both Collagen I transcription and fibrillogenesis by binding inside the Col1a1 gene body and maintaining RNA polymerase II occupancy. In vivo, Prdm5 loss results in delayed ossification involving a pronounced impairment in the assembly of fibrillar collagens. Collectively, our results define a novel role for Prdm5 in sustaining the transcriptional program necessary to the proper assembly of osteoblastic extracellular matrix.

Additional copies of 1q negatively impact the outcome of multiple myeloma patients and induce transcriptomic deregulation in malignant plasma cells.

Additional copies of chromosome 1 long arm (1q) are frequently found in multiple myeloma (MM) and predict high-risk disease. Available data suggest a different outcome and biology of patients with amplification (Amp1q, ≥4 copies of 1q) vs. gain (Gain1q, 3 copies of 1q) of 1q. We evaluated the impact of Amp1q/Gain1q on the outcome of newly diagnosed MM patients enrolled in the FORTE trial (NCT02203643). Among 400 patients with available 1q data, 52 (13%) had Amp1q and 129 (32%) Gain1q. After a median follow-up of 62 months, median progression-free survival (PFS) was 21.2 months in the Amp1q group, 54.9 months in Gain1q, and not reached (NR) in Normal 1q. PFS was significantly hampered by the presence of Amp1q (HR 3.34 vs. Normal 1q, P < 0.0001; HR 1.99 vs. Gain1q, P = 0.0008). Patients with Gain1q had also a significantly shorter PFS compared with Normal 1q (HR 1.68, P = 0.0031). Concomitant poor prognostic factors or the failure to achieve MRD negativity predicted a median PFS < 12 months in Amp1q patients. Carfilzomib-lenalidomide-dexamethasone plus autologous stem cell transplantation treatment improved the adverse effect of Gain1q but not Amp1q. Transcriptomic data showed that additional 1q copies were associated with deregulation in apoptosis signaling, p38 MAPK signaling, and Myc-related genes.

Circulating hsa-miR-5096 predicts 18F-FDG PET/CT positivity and modulates somatostatin receptor 2 expression: a novel miR-based assay for pancreatic neuroendocrine tumors.

Gastro-entero-pancreatic neuroendocrine tumors (GEP-NETs) are rare diseases encompassing pancreatic (PanNETs) and ileal NETs (SINETs), characterized by heterogeneous somatostatin receptors (SSTRs) expression. Treatments for inoperable GEP-NETs are limited, and SSTR-targeted Peptide Receptor Radionuclide Therapy (PRRT) achieves variable responses. Prognostic biomarkers for the management of GEP-NET patients are required. 18F-FDG uptake is a prognostic indicator of aggressiveness in GEP-NETs. This study aims to identify circulating and measurable prognostic miRNAs associated with 18F-FDG-PET/CT status, higher risk and lower response to PRRT. Methods: Whole miRNOme NGS profiling was conducted on plasma samples obtained from well-differentiated advanced, metastatic, inoperable G1, G2 and G3 GEP-NET patients enrolled in the non-randomized LUX (NCT02736500) and LUNET (NCT02489604) clinical trials prior to PRRT (screening set, n= 24). Differential expression analysis was performed between 18F-FDG positive (n=12) and negative (n=12) patients. Validation was conducted by Real Time quantitative PCR in two distinct well-differentiated GEP-NET validation cohorts, considering the primary site of origin (PanNETs n=38 and SINETs n=30). The Cox regression was applied to assess independent clinical parameters and imaging for progression-free survival (PFS) in PanNETs. In situ RNA hybridization combined with immunohistochemistry was performed to simultaneously detect miR and protein expression in the same tissue specimens. This novel semi-automated miR-protein protocol was applied in PanNET FFPE specimens (n=9). In vitro functional experiments were performed in PanNET models. Results: While no miRNAs emerged to be deregulated in SINETs, hsa-miR-5096, hsa-let-7i-3p and hsa-miR-4311 were found to correlate with 18F-FDG-PET/CT in PanNETs (p-value:<0.005). Statistical analysis has shown that, hsa-miR-5096 can predict 6-month PFS (p-value:<0.001) and 12-month Overall Survival upon PRRT treatment (p-value:<0.05), as well as identify 18F-FDG-PET/CT positive PanNETs with worse prognosis after PRRT (p-value:<0.005). In addition, hsa-miR-5096 inversely correlated with both SSTR2 expression in PanNET tissue and with the 68Gallium-DOTATOC captation values (p-value:<0.05), and accordingly it was able to decrease SSTR2 when ectopically expressed in PanNET cells (p-value:<0.01). Conclusions: hsa-miR-5096 well performs as a biomarker for 18F-FDG-PET/CT and as independent predictor of PFS. Moreover, exosome-mediated delivery of hsa-miR-5096 may promote SSTR2 heterogeneity and thus resistance to PRRT.

Circulating miRNome of Trachemys scripta after elective gonadectomy under general anesthesia.

Post-surgical management is an important issue in veterinary medicine, requiring biomarkers with high sensitivity and specificity for timely and effective treatment. Emerging evidence suggests that miRNAs are promising stress- and pain-related markers. The aims were to profile the circulating miRNA signature in plasma of turtles (Trachemys scripta) and point out potential candidate biomarkers to assess the status of the animal. The plasma of female turtles underwent surgical gonadectomy were collected 24 h pre-surgery, and 2.5 h and 36 h post-surgery. The expression of miRNAs was profiled by Next Generation Sequencing and the dysregulated miRNAs were validated using RT-qPCR. The diagnostic value of miRNAs was calculated by ROC curves. The results showed that 14 miRNAs were differentially expressed over time. RT-qPCR validation highlighted that 2-miR-499-3p and miR-203-5p-out of 8 miRNAs tested were effectively modulated. The Area Under the Curve (AUC) of miR-203-5p was fair (AUC 0.7934) in discriminating pre- and 36 h post-surgery samples and poor for other time points; the AUC of miR-499-3p was excellent (AUC 0.944) in discriminating pre-surgery and 2.5 h post-surgery samples, and fair in discriminating pre-surgery and 36 h post-surgery (AUC 0.7292) and 2.5 h and 36 h post-surgery (AUC 0.7569) samples. In conclusion, we demonstrated for the first time that miRNAs profile changes in plasma of turtles underwent surgical oophorectomy and identified miR-203-5p and miR-499-3p as potential candidate biomarkers to assess animals' status. Further studies are necessary to confirm their diagnostic value and to investigate functional and mechanistic networks to improve our understanding of the biological processes.

Identification of TENM4 as a Novel Cancer Stem Cell-Associated Molecule and Potential Target in Triple Negative Breast Cancer.

Triple-negative breast cancer (TNBC) is insensitive to endocrine and Her2-directed therapies, making the development of TNBC-targeted therapies an unmet medical need. Since patients with TNBC frequently show a quicker relapse and metastatic progression compared to other breast cancer subtypes, we hypothesized that cancer stem cells (CSC) could have a role in TNBC. To identify putative TNBC CSC-associated targets, we compared the gene expression profiles of CSC-enriched tumorspheres and their parental cells grown as monolayer. Among the up-regulated genes coding for cell membrane-associated proteins, we selected Teneurin 4 (TENM4), involved in cell differentiation and deregulated in tumors of different histotypes, as the object for this study. Meta-analysis of breast cancer datasets shows that TENM4 mRNA is up-regulated in invasive carcinoma specimens compared to normal breast and that high expression of TENM4 correlates with a shorter relapse-free survival in TNBC patients. TENM4 silencing in mammary cancer cells significantly impaired tumorsphere-forming ability, migratory capacity and Focal Adhesion Kinase (FAK) phosphorylation. Moreover, we found higher levels of TENM4 in plasma from tumor-bearing mice and TNBC patients compared to the healthy controls. Overall, our results indicate that TENM4 may act as a novel biomarker and target for the treatment of TNBC.

AP-2alpha regulates migration of GN-11 neurons via a specific genetic programme involving the Axl receptor tyrosine kinase.

BACKGROUND: Neuronal migration is a crucial process that allows neurons to reach their correct target location to allow the nervous system to function properly. AP-2alpha is a transcription factor essential for neural crest cell migration and its mutation results in apoptosis within this cell population, as demonstrated by genetic models. RESULTS: We down-modulated AP-2alpha expression in GN-11 neurons by RNA interference and observe reduced neuron migration following the activation of a specific genetic programme including the Adhesion Related Kinase (Axl) gene. We prove that Axl is able to coordinate migration per se and by ChIP and promoter analysis we observe that its transcription is directly driven by AP-2alpha via the binding to one or more functional AP-2alpha binding sites present in its regulatory region. Analysis of migration in AP-2alpha null mouse embryo fibroblasts also reveals an essential role for AP-2alpha in cell movement via the activation of a distinct genetic programme. CONCLUSION: We show that AP-2alpha plays an essential role in cell movement via the activation of cell-specific genetic programmes. Moreover, we demonstrate that the AP-2alpha regulated gene Axl is an essential player in GN-11 neuron migration.

Inflammation and breast cancer. Inflammatory component of mammary carcinogenesis in ErbB2 transgenic mice.

This review addresses genes differentially expressed in the mammary gland transcriptome during the progression of mammary carcinogenesis in BALB/c mice that are transgenic for the rat neu (ERBB2, or HER-2/neu) oncogene (BALB-neuT664V-E mice). The Ingenuity knowledge database was used to characterize four functional association networks whose hub genes are directly linked to inflammation (specifically, the genes encoding IL-1beta, tumour necrosis factor, interferon-gamma, and monocyte chemoattractant protein-1/CC chemokine ligand-2) and are increasingly expressed during such progression. In silico meta-analysis in a human breast cancer dataset suggests that proinflammatory activation in the mammary glands of these mice reflects a general pattern of human breast cancer.

Immunotargeting of Antigen xCT Attenuates Stem-like Cell Behavior and Metastatic Progression in Breast Cancer.

Resistance to therapy and lack of curative treatments for metastatic breast cancer suggest that current therapies may be missing the subpopulation of chemoresistant and radioresistant cancer stem cells (CSC). The ultimate success of any treatment may well rest on CSC eradication, but specific anti-CSC therapies are still limited. A comparison of the transcriptional profiles of murine Her2(+) breast tumor TUBO cells and their derived CSC-enriched tumorspheres has identified xCT, the functional subunit of the cystine/glutamate antiporter system xc(-), as a surface protein that is upregulated specifically in tumorspheres. We validated this finding by cytofluorimetric analysis and immunofluorescence in TUBO-derived tumorspheres and in a panel of mouse and human triple negative breast cancer cell-derived tumorspheres. We further show that downregulation of xCT impaired tumorsphere generation and altered CSC intracellular redox balance in vitro, suggesting that xCT plays a functional role in CSC biology. DNA vaccination based immunotargeting of xCT in mice challenged with syngeneic tumorsphere-derived cells delayed established subcutaneous tumor growth and strongly impaired pulmonary metastasis formation by generating anti-xCT antibodies able to alter CSC self-renewal and redox balance. Finally, anti-xCT vaccination increased CSC chemosensitivity to doxorubicin in vivo, indicating that xCT immunotargeting may be an effective adjuvant to chemotherapy.

A hypoxic signature marks tumors formed by disseminated tumor cells in the BALB-neuT mammary cancer model.

Metastasis is the final stage of cancer progression. Some evidence indicates that tumor cell dissemination occurs early in the natural history of cancer progression. Disseminated tumor cells (DTC) have been described in the bone marrow (BM) of cancer patients as well as in experimental models, where they correlate with later development of metastasis. However, little is known about the tumorigenic features of DTC obtained at different time points along tumor progression. Here, we found that early DTC isolated from BM of 15-17 week-old Her2/neu transgenic (BALB-neuT) mice were not tumorigenic in immunodeficient mice. In contrast, DTC-derived tumors were easily detectable when late DTC obtained from 19-22 week-old BALB-neuT mice were injected. Angiogenesis, which contributes to regulate tumor dormancy, appeared dispensable to reactivate late DTC, although it accelerated growth of secondary DTC tumors. Compared with parental mammary tumors, gene expression profiling disclosed a distinctive transcriptional signature of late DTC tumors which was enriched for hypoxia-related transcripts and was maintained in ex-vivo cell culture. Altogether, these findings highlight a different tumorigenic potential of early and late DTC in the BALB-neuT model and describe a HIF-1α-related transcriptional signature in DTC tumors, which may render DTC angiogenesis-competent, when placed in a favourable environment.

A novel infection- and inflammation-associated molecular signature in peripheral blood of myasthenia gravis patients.

Myasthenia gravis (MG) is a T-cell dependent autoimmune disorder of the neuromuscular junction, characterised by muscle weakness and fatigability. Autoimmunity is thought to initiate in the thymus of acetylcholine receptor (AChR)-positive MG patients; however, the molecular mechanisms linking intra-thymic MG pathogenesis with autoreactivity via the circulation to the muscle target organ are poorly understood. Using whole-transcriptome sequencing, we compared the transcriptional profile of peripheral blood mononuclear cells from AChR-early onset MG (AChR-EOMG) patients with healthy controls: 178 coding transcripts and 229 long non-coding RNAs, including 11 pre-miRNAs, were differentially expressed. Among the 178 coding transcripts, 128 were annotated of which 17% were associated with the 'infectious disease' functional category and 46% with 'inflammatory disease' and 'inflammatory response-associated' categories. Validation of selected transcripts by qPCR indicated that of the infectious disease-related transcripts, ETF1, NFKB2, PLK3, and PPP1R15A were upregulated, whereas CLC and IL4 were downregulated in AChR-EOMG patients; in the 'inflammatory' categories, ABCA1, FUS, and RELB were upregulated, suggesting a contribution of these molecules to immunological dysfunctions in MG. Data selection and validation were also based on predicted microRNA-mRNA interactions. We found that miR-612, miR-3654, and miR-3651 were increased, whereas miR-612-putative AKAp12 and HRH4 targets and the miR-3651-putative CRISP3 target were downregulated in AChR-EOMG, also suggesting altered immunoregulation. Our findings reveal a novel peripheral molecular signature in AChR-EOMG, reflecting a critical involvement of inflammatory- and infectious disease-related immune responses in disease pathogenesis.

Novel insights into Notum and glypicans regulation in colorectal cancer.

The connection between colorectal cancer (CRC) and Wnt signaling pathway activation is well known, but full elucidation of the underlying regulation of the Wnt/ß-catenin pathway and its biological functions in CRC pathogenesis is still needed. Here, the azoxymethane/dextran sulfate sodium salt (AOM/DSS) murine model has been used as an experimental platform able to mimic human sporadic CRC development with predictable timing. We performed genome-wide expression profiling of AOM/DSS-induced tumors and normal colon mucosa to identify potential novel CRC biomarkers. Remarkably, the enhanced expression of Notum, a conserved feedback antagonist of Wnt, was observed in tumors along with alterations in Glypican-1 and Glypican-3 levels. These findings were confirmed in a set of human CRC samples. Here, we provide the first demonstration of significant changes in Notum and glypicans gene expression during CRC development and present evidence to suggest them as potential new biomarkers of CRC pathogenesis.

Genomic and proteomic analyses of Prdm5 reveal interactions with insulator binding proteins in embryonic stem cells.

PRDM proteins belong to the SET domain protein family, which is involved in the regulation of gene expression. Although few PRDM members possess histone methyltransferase activity, the molecular mechanisms by which the other members exert transcriptional regulation remain to be delineated. In this study, we find that Prdm5 is highly expressed in mouse embryonic stem (mES) cells and exploit this cellular system to characterize molecular functions of Prdm5. By combining proteomics and next-generation sequencing technologies, we identify Prdm5 interaction partners and genomic occupancy. We demonstrate that although Prdm5 is dispensable for mES cell maintenance, it directly targets genomic regions involved in early embryonic development and affects the expression of a subset of developmental regulators during cell differentiation. Importantly, Prdm5 interacts with Ctcf, cohesin, and TFIIIC and cooccupies genomic loci. In summary, our data indicate how Prdm5 modulates transcription by interacting with factors involved in genome organization in mouse embryonic stem cells.

Cell-to-cell signaling influences the fate of prostate cancer stem cells and their potential to generate more aggressive tumors.

An increasing number of malignancies has been shown to be initiated and propelled by small subpopulations of cancer stem cells (CSC). However, whether tumor aggressiveness is driven by CSC and by what extent this property may be relevant within the tumor mass is still unsettled. To address this issue, we isolated a rare tumor cell population on the basis of its CD44(+)CD24(-) phenotype from the human androgen-independent prostate carcinoma cell line DU145 and established its CSC properties. The behavior of selected CSC was investigated with respect to the bulk DU145 cells. The injection of CSC in nude mice generated highly vascularized tumors infiltrating the adjacent tissues, showing high density of neuroendocrine cells and expressing low levels of E-cadherin and ß-catenin as well as high levels of vimentin. On the contrary, when a comparable number of unsorted DU145 cells were injected the resulting tumors were less aggressive. To investigate the different features of tumors in vivo, the influence of differentiated tumor cells on CSC was examined in vitro by growing CSC in the absence or presence of conditioned medium from DU145 cells. CSC grown in permissive conditions differentiated into cell populations with features similar to those of cells held in aggressive tumors generated from CSC injection. Differently, conditioned medium induced CSC to differentiate into a cell phenotype comparable to cells of scarcely aggressive tumors originated from bulk DU145 cell injection. These findings show for the first time that CSC are able to generate differentiated cells expressing either highly or scarcely aggressive phenotype, thus influencing prostate cancer progression. The fate of CSC was determined by signals released from tumor environment. Moreover, using microarray analysis we selected some molecules which could be involved in this cell-to-cell signaling, hypothesizing their potential value for prognostic or therapeutic applications.

Oncoantigens as anti-tumor vaccination targets: the chance of a lucky strike?

Neoplastic transformation is a multistage process and distinct gene products of specific cell regulatory pathways are involved at each stage. Identification of genes overexpressed at a specific stage provides an unprecedented opportunity to address the immune system against antigens with a driving role in tumor progression (oncoantigens). The ERBB2 oncogene is a prototype of deregulated oncogenic protein kinase membrane receptors. Mice transgenic for rat ERBB2 (BALB-neuT mice) were used in this study to identify an additional set of oncoantigens expressed at defined stages by most breast carcinomas to be used as alternatives to ERBB2-driven vaccination. To address this question, we integrated the transcription data generated by comparing preneoplastic lesions and neoplasia in BALB-neuT mice with a meta-analysis on transcription profiles generated from normal and breast tumor human specimens. Forty-six putative oncoantigens identified and prioritized according to their expression on the cell membrane or in the extra cellular space, cytoplasm and nucleus were chosen for preclinical investigation as vaccination targets.