YBX1-interacting small RNAs and RUNX2 can be blocked in primary bone cancer using CADD522.

Primary bone cancer (PBC) comprises several subtypes each underpinned by distinctive genetic drivers. This driver diversity produces novel morphological features and clinical behaviour that serendipitously makes PBC an excellent metastasis model. Here, we report that some transfer RNA-derived small RNAs termed tRNA fragments (tRFs) perform as a constitutive tumour suppressor mechanism by blunting a potential pro-metastatic protein-RNA interaction. This mechanism is reduced in PBC progression with a gradual loss of tRNAGlyTCC cleavage into 5' end tRF-GlyTCC when comparing low-grade, intermediate-grade and high-grade patient tumours. We detected recurrent activation of miR-140 leading to upregulated RUNX2 expression in high-grade patient tumours. Both tRF-GlyTCC and RUNX2 share a sequence motif in their 3' ends that matches the YBX1 recognition site known to stabilise pro-metastatic mRNAs. Investigating some aspects of this interaction network, gain- and loss-of-function experiments using small RNA mimics and antisense LNAs, respectively, showed that ectopic tRF-GlyTCC reduced RUNX2 expression and dispersed 3D micromass architecture in vitro. iCLIP sequencing revealed YBX1 physical binding to the 3' UTR of RUNX2. The interaction between YBX1, tRF-GlyTCC and RUNX2 led to the development of the RUNX2 inhibitor CADD522 as a PBC treatment. CADD522 assessment in vitro revealed significant effects on PBC cell behaviour. In xenograft mouse models, CADD522 as a single agent without surgery significantly reduced tumour volume, increased overall and metastasis-free survival and reduced cancer-induced bone disease. Our results provide insight into PBC molecular abnormalities that have led to the identification of new targets and a new therapeutic.

Transcriptome analysis of long noncoding RNAs reveals their potential roles in anthracycline-induced cardiotoxicity.

Aims: Anthracyclines (ANTs) are essential chemotherapeutic agents; however, their adverse effects can lead to heart failure in cancer survivors. While long non-coding RNAs (lncRNAs) have become new players in cellular processes, there is limited knowledge on lncRNA expression related to anthracyclines-induced cardiotoxicity. This study investigates the lncRNA profiles in human cardiac microtissues exposed to 3 popular ANTs, namely doxorubicin, epirubicin, and idarubicin, as well as in heart biopsies from ANT-treated patients. Methods and results: The in vitro microtissues were exposed to each ANT at 2 doses over 2 weeks; the transcriptome data was collected at 7 time points. The human biopsies were collected from heart failure patients who underwent ANT treatment and control subjects. Over 100 lncRNAs were differentially expressed in each in vitro ANT treatment condition compared to control samples; 16 of them were differentially expressed across all ANT-treated conditions. The lncRNA databases and literature revealed insight on how these lncRNAs relate to heart failure and cellular functions. For instance, H19 and RMRP are involved in heart failure progression, while BDNF-AS is a cardiomyocyte damage-associated gene; SNHG7 is a cardiac hypertrophy regulator. PCAT19 can promote the miR-182/PDK4 axis and modulate p53 expression, whereas SNHG29 can regulate the Wnt/ß-catenin signaling pathway via the miR-223-3p/CTNND1 axis. Other lncRNAs, which were only differentially expressed in particular ANT-treated conditions, are also involved in cardiomyocyte damage and heart failure disease. The alterations of these lncRNA expressions in the in vitro cardiac tissue were also affirmed by similar changes in the human biopsies. Conclusion: This study revealed several lncRNAs that can be potential biomarkers or targets for further ANT-induced cardiotoxicity investigation, according to the transcriptome in both human cardiac microtissues expose to ANTs as well as in heart biopies form ANT-treated patients. Especially, H19 lncRNA showed its contribution to on-target toxicity, in which it is involved in both chemoresistance and cardiotoxic mechanism.

Alterations in 3D chromatin organization contribute to tumorigenesis of EGFR-amplified glioblastoma.

Background: EGFR amplification and/or mutation are found in more than half of the cases with glioblastoma. Yet, the role of chromatin interactions and its regulation of gene expression in EGFR-amplified glioblastoma remains unclear. Methods: In this study, we explored alterations in 3D chromatin organization of EGFR-amplified glioblastoma and its subsequent impact by performing a comparative analysis of Hi-C, RNA-seq, and whole-genome sequencing (WGS) on EGFR-amplified glioblastoma-derived A172 and normal astrocytes (HA1800 cell line). Results: A172 cells showed an elevated chromatin relaxation, and unexpected entanglement of chromosome regions. A genome-wide landscape of switched compartments and differentially expressed genes between HA1800 and A172 cell lines demonstrated that compartment activation reshaped chromatin accessibility and activated tumorigenesis-related genes. Topological associating domain (TAD) analysis revealed that altered TAD domains in A172 also contribute to oncogene activation and tumor repressor deactivation. Interestingly, glioblastoma-derived A172 cells showed a different chromatin loop contact propensity. Genes in tumorigenesis-associated signaling pathways were significantly enriched at the anchor loci of altered chromatin loops. Oncogene activation and tumor repressor deactivation were associated with chromatin loop alteration. Structure variations (SVs) had a dramatic impact on the chromatin conformation of EGFR-amplified glioblastoma-derived tumor cells. Moreover, our results revealed that 7p11.2 duplication activated EGFR expression in EGFR-amplified glioblastoma via neo-TAD formation and novel enhancer-promoter interaction emergence between LINC01446 and EGFR. Conclusions: The disordered 3D genomic map and multi-omics data of EGFR-amplified glioblastoma provide a resource for future interrogation of the relationship between chromatin interactions and transcriptome in tumorigenesis.

Transcriptomic Profiling of Plaque Psoriasis and Cutaneous T-Cell Subsets during Treatment with Secukinumab.

The IL-17A inhibitor secukinumab is efficacious for the treatment of psoriasis. To better understand its mechanism of action, we investigated its impact on psoriatic lesions from 15 patients with moderate-to-severe plaque psoriasis undergoing secukinumab treatment. We characterized the longitudinal transcriptomic changes of whole lesional skin tissue as well as cutaneous CD4+ and CD8+ T effector cells and CD4+ T regulatory cells across 12 weeks of treatment. Secukinumab was clinically effective and reduced disease-associated overexpression of IL17A , IL17F, IL23A, IL23R, and IFNG in whole tissue as soon as 2 weeks after initiation of treatment. IL17A overexpression in T-cell subsets, primarily CD8+ T cells, was also reduced. Although secukinumab treatment resolved 89‒97% of psoriasis-associated expression differences in bulk tissue and T-cell subsets by week 12 of treatment, we observed expression differences involved in IFN signaling and metallothionein synthesis that remained unresolved at this time point as well as potential treatment-associated expression differences involved in IL-15 signaling. These changes were accompanied by shifts in broader immune cell composition on the basis of deconvolution of RNA-sequencing data. In conclusion, our study reveals several phenotypic and cellular changes within the lesion that underlie clinical improvement from secukinumab.

Mechanistic insights into the digestion of complex dietary fibre by the rumen microbiota using combinatorial high-resolution glycomics and transcriptomic analyses.

There is a knowledge gap regarding the factors that impede the ruminal digestion of plant cell walls or if rumen microbiota possess the functional activities to overcome these constraints. Innovative experimental methods were adopted to provide a high-resolution understanding of plant cell wall chemistries, identify higher-order structures that resist microbial digestion, and determine how they interact with the functional activities of the rumen microbiota. We characterized the total tract indigestible residue (TTIR) from cattle fed a low-quality straw diet using two comparative glycomic approaches: ELISA-based glycome profiling and total cell wall glycosidic linkage analysis. We successfully detected numerous and diverse cell wall glycan epitopes in barley straw (BS) and TTIR and determined their relative abundance pre- and post-total tract digestion. Of these, xyloglucans and heteroxylans were of higher abundance in TTIR. To determine if the rumen microbiota can further saccharify the residual plant polysaccharides within TTIR, rumen microbiota from cattle fed a diet containing BS were incubated with BS and TTIR ex vivo in batch cultures. Transcripts coding for carbohydrate-active enzymes (CAZymes) were identified and characterized for their contribution to cell wall digestion based on glycomic analyses, comparative gene expression profiles, and associated CAZyme families. High-resolution phylogenetic fingerprinting of these sequences encoded CAZymes with activities predicted to cleave the primary linkages within heteroxylan and arabinan. This experimental platform provides unprecedented precision in the understanding of forage structure and digestibility, which can be extended to other feed-host systems and inform next-generation solutions to improve the performance of ruminants fed low-quality forages.

Integrative RNA profiling of TBEV-infected neurons and astrocytes reveals potential pathogenic effectors.

Tick-borne encephalitis virus (TBEV), the most medically relevant tick-transmitted flavivirus in Eurasia, targets the host central nervous system and frequently causes severe encephalitis. The severity of TBEV-induced neuropathogenesis is highly cell-type specific and the exact mechanism responsible for such differences has not been fully described yet. Thus, we performed a comprehensive analysis of alterations in host poly-(A)/miRNA/lncRNA expression upon TBEV infection in vitro in human primary neurons (high cytopathic effect) and astrocytes (low cytopathic effect). Infection with severe but not mild TBEV strain resulted in a high neuronal death rate. In comparison, infection with either of TBEV strains in human astrocytes did not. Differential expression and splicing analyses with an in silico prediction of miRNA/mRNA/lncRNA/vd-sRNA networks found significant changes in inflammatory and immune response pathways, nervous system development and regulation of mitosis in TBEV Hypr-infected neurons. Candidate mechanisms responsible for the aforementioned phenomena include specific regulation of host mRNA levels via differentially expressed miRNAs/lncRNAs or vd-sRNAs mimicking endogenous miRNAs and virus-driven modulation of host pre-mRNA splicing. We suggest that these factors are responsible for the observed differences in the virulence manifestation of both TBEV strains in different cell lines. This work brings the first complex overview of alterations in the transcriptome of human astrocytes and neurons during the infection by two TBEV strains of different virulence. The resulting data could serve as a starting point for further studies dealing with the mechanism of TBEV-host interactions and the related processes of TBEV pathogenesis.

The P2RX7B splice variant modulates osteosarcoma cell behaviour and metastatic properties.

Background: Osteosarcoma (OS) is the most common type of primary bone cancer affecting children and adolescents. OS has a high propensity to spread meaning the disease is often incurable and fatal. There have been no improvements in survival rates for decades. This highlights an urgent need for the development of novel therapeutic strategies. Here, we report in vitro and in vivo data that demonstrates the role of purinergic signalling, specifically, the B isoform of the purinergic receptor P2RX7 (P2RX7B), in OS progression and metastasis. Methods: TE85 and MNNG-HOS OS cells were transfected with P2RX7B. These cell lines were then characterised and assessed for proliferation, cell adhesion, migration and invasion in vitro. We used these cells to perform both paratibial and tail vein injected mouse studies where the primary tumour, bone and lungs were analysed. We used RNA-seq to identify responsive pathways relating to P2RX7B. Results: Our data shows that P2RX7B expression confers a survival advantage in TE85 + P2RX7B and MNNG-HOS + P2RX7B human OS cell lines in vitro that is minimised following treatment with A740003, a specific P2RX7 antagonist. P2RX7B expression reduced cell adhesion and P2RX7B activation promoted invasion and migration in vitro, demonstrating a metastatic phenotype. Using an in vivo OS xenograft model, MNNG-HOS + P2RX7B tumours exhibited cancer-associated ectopic bone formation that was abrogated with A740003 treatment. A pro-metastatic phenotype was further demonstrated in vivo as expression of P2RX7B in primary tumour cells increased the propensity of tumour cells to metastasise to the lungs. RNA-seq identified a novel gene axis, FN1/LOX/PDGFB/IGFBP3/BMP4, downregulated in response to A740003 treatment. Conclusion: Our data illustrates a role for P2RX7B in OS tumour growth, progression and metastasis. We show that P2RX7B is a future therapeutic target in human OS.

A pooled CRISPR/AsCpf1 screen using paired gRNAs to induce genomic deletions in Chinese hamster ovary cells.

Chinese hamster ovary (CHO) cells are the most widely used host for the expression of therapeutic proteins. Recently, significant progress has been made due to advances in genome sequence and annotation quality to unravel the black box CHO. Nevertheless, in many cases the link between genotype and phenotype in the context of suspension cultivated production cell lines is still not fully understood. While frameshift approaches targeting coding genes are frequently used, the non-coding regions of the genome have received less attention with respect to such functional annotation. Importantly, for non-coding regions frameshift knock-out strategies are not feasible. In this study, we developed a CRISPR-mediated screening approach that performs full deletions of genomic regions to enable the functional study of both the translated and untranslated genome. An in silico pipeline for the computational high-throughput design of paired guide RNAs (pgRNAs) directing CRISPR/AsCpf1 was established and used to generate a library tackling process-related genes and long non-coding RNAs. Next generation sequencing analysis of the plasmid library revealed a sufficient, but highly variable pgRNA composition. Recombinase-mediated cassette exchange was applied for pgRNA library integration rather than viral transduction to ensure single copy representation of pgRNAs per cell. After transient AsCpf1 expression, cells were cultivated over two sequential batches to identify pgRNAs which massively affected growth and survival. By comparing pgRNA abundance, depleted candidates were identified and individually validated to verify their effect.

ceRNA network development and tumour-infiltrating immune cell analysis of metastatic breast cancer to bone.

PURPOSE: Advanced breast cancer commonly metastasises to bone; however, the molecular mechanisms underlying the affinity for breast cancer cells to bone remains unclear. Thus, we developed nomograms based on a competing endogenous RNA (ceRNA) network and analysed tumour-infiltrating immune cells to elucidate the molecular pathways that may predict prognosis in patients with breast cancer. METHODS: We obtained the RNA expression profile of 1091 primary breast cancer samples included in The Cancer Genome Atlas database, 58 of which were from patients with bone metastasis. We analysed the differential RNA expression patterns between breast cancer with and without bone metastasis and developed a ceRNA network. Cibersort was employed to differentiate between immune cell types based on tumour transcripts. Nomograms were then established based on the ceRNA network and immune cell analysis. The value of prognostic factors was evaluated by Kaplan-Meier survival analysis and a Cox proportional risk model. RESULTS: We found significant differences in long non-coding RNAs (lncRNAs), 18 microRNAs (miRNAs), and 20 messenger RNAs (mRNAs) between breast cancer with and without bone metastasis, which were used to construct a ceRNA network. We found that the protein-coding genes GJB3, CAMMV, PTPRZ1, and FBN3 were significantly differentially expressed by Kaplan-Meier analysis. We also observed significant differences in the abundance of plasma cell and follicular helper T cell populations between the two groups. In addition, the proportion of mast cells, gamma delta T cells, and plasma cells differed depending on disease location and stage. Our analysis showed that a high proportion of follicular helper T cells and a low proportion of eosinophils promoted survival and that DLX6-AS1, Wnt6, and GABBR2 expression may be associated with bone metastasis in breast cancer. CONCLUSIONS: We developed a bioinformatic tool for exploring the molecular mechanisms of bone metastasis in patients with breast cancer and identified factors that may predict the occurrence of bone metastasis.

Data from differentially expressed proteins in platelet components associated with adverse transfusion reactions.

The presented dataset was used for the study focused on the search for differentially expressed proteins in blood platelet components (PCs) associated with adverse transfusion reactions (ATRs). Pellets of ATR platelet components and their controls were subjected to high-throughput proteomics analysis using a Q Exactive high-resolution tandem mass spectrometer. The data reported here constitutes an extension of "Differential protein expression of blood platelet components associated with adverse transfusion reactions" article Aloui et al., 2018. The reported data herein have been deposited into the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifiers PXD003510 for the pooled platelet components (PPCs) and PXD008886 for the apheresis platelet components (SDA-PCs) associated with ATRs.

Radiogenomic Analysis of F-18-Fluorodeoxyglucose Positron Emission Tomography and Gene Expression Data Elucidates the Epidemiological Complexity of Colorectal Cancer Landscape.

PURPOSE: Transcriptomic profiling has enabled the neater genomic characterization of several cancers, among them colorectal cancer (CRC), through the derivation of genes with enhanced causal role and informative gene sets. However, the identification of small-sized gene signatures, which can serve as potential biomarkers in CRC, remains challenging, mainly due to the great genetic heterogeneity of the disease. METHODS: We developed and exploited an analytical framework for the integrative analysis of CRC datasets, encompassing transcriptomic data and positron emission tomography (PET) measurements. Profiling data comprised two microarray datasets, pertaining biopsy specimen from 30 untreated patients with primary CRC, coupled by their F-18-Fluorodeoxyglucose (FDG) PET values, using tracer kinetic analysis measurements. The computational framework incorporates algorithms for semantic processing, multivariate analysis, data mining and dimensionality reduction. RESULTS: Transcriptomic and PET data feature sets, were evaluated for their discrimination performance between primary colorectal adenocarcinomas and adjacent normal mucosa. A composite signature was derived, pertaining 12 features: 7 genes and 5 PET variables. This compact signature manifests superior performance in classification accuracy, through the integration of gene expression and PET data. CONCLUSIONS: This work represents an effort for the integrative, multilayered, signature-oriented analysis of CRC, in the context of radio-genomics, inferring a composite signature with promising results for patient stratification.

Microarray data and gene expression statistics for Saccharomyces cerevisiae exposed to simulated asbestos mine drainage.

Here we describe microarray expression data (raw and normalized), experimental metadata, and gene-level data with expression statistics from Saccharomyces cerevisiae exposed to simulated asbestos mine drainage from the Vermont Asbestos Group (VAG) Mine on Belvidere Mountain in northern Vermont, USA. For nearly 100 years (between the late 1890s and 1993), chrysotile asbestos fibers were extracted from serpentinized ultramafic rock at the VAG Mine for use in construction and manufacturing industries. Studies have shown that water courses and streambeds nearby have become contaminated with asbestos mine tailings runoff, including elevated levels of magnesium, nickel, chromium, and arsenic, elevated pH, and chrysotile asbestos-laden mine tailings, due to leaching and gradual erosion of massive piles of mine waste covering approximately 9 km2. We exposed yeast to simulated VAG Mine tailings leachate to help gain insight on how eukaryotic cells exposed to VAG Mine drainage may respond in the mine environment. Affymetrix GeneChip® Yeast Genome 2.0 Arrays were utilized to assess gene expression after 24-h exposure to simulated VAG Mine tailings runoff. The chemistry of mine-tailings leachate, mine-tailings leachate plus yeast extract peptone dextrose media, and control yeast extract peptone dextrose media is also reported. To our knowledge this is the first dataset to assess global gene expression patterns in a eukaryotic model system simulating asbestos mine tailings runoff exposure. Raw and normalized gene expression data are accessible through the National Center for Biotechnology Information Gene Expression Omnibus (NCBI GEO) Database Series GSE89875 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE89875).