Analysis of the Expression of Repetitive DNA Elements in Osteosarcoma

Extracellular Vesicles: An Emerging Clinical Opportunity in Musculoskeletal Disease

Extracellular vesicles (EVs) are important mediators of cell-to-cell communication in the extracellular space. These membranous nanoparticles carry various molecules, often referred to as "cargo," which are delivered to nearby target cells. In the past decade, developments in nanotechnology have allowed for various new laboratory techniques for the increased utilization of EVs in cellular and animal studies. Such techniques have evolved for the isolation, characterization, and delivery of EVs to biological tissues. This emerging technology has immense clinical potential for both diagnostic and therapeutic applications. Various EV cargo molecules, including DNA, RNA, and proteins, can act as pathological biomarkers. Furthermore, EVs derived from certain cell sources have shown therapeutic benefit in certain pathologies. In addition to their native therapeutic benefit, EVs can be engineered to carry and selectively deliver therapeutic agents. While EVs have gained increasing interest in various pathologies, few studies have compiled their clinical potential in musculoskeletal pathologies. To bridge this gap, we present an overview of EVs, introduce current laboratory preparation techniques, and outline the most recent literature regarding the potential therapeutic applications of EVs in musculoskeletal pathologies.

Comprehensive analysis to identify PUS7 as a prognostic biomarker from pan-cancer analysis to osteosarcoma validation

AIM

Pseudouridylation has demonstrated the potential to control the development of numerous malignancies. PUS7(Pseudouridine Synthase 7) is one of the pseudouridine synthases, but the literature on this enzyme is limited to several cancer types. Currently, no investigation has been performed on the systematic pan-cancer analysis concerning PUS7 role in cancer diagnosis and prognosis.

METHODS

Employing public databases, including The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression Project (GTEx), Human Protein Atlas (HPA), UALCAN and Tumor Immune Single-cell Hub (TISCH), this work investigated the PUS7 carcinogenesis in pan-cancer. Differential expression analysis, prognostic survival analysis and biological function were systematically performed. Furthermore, PUS7 potential as an osteosarcoma biomarker for diagnosis and prognosis was assessed in this study.

RESULTS

The findings indicated that PUS7 was overexpressed in the majority of malignancies. High PUS7 expression contributed to the poor prognosis among 11 cancer types, including Adrenocortical Cancer (ACC), Bladder Cancer (BLCA), Liver Cancer (LIHC), Kidney Papillary Cell Carcinoma (KIRP), Mesothelioma (MESO), Lower Grade Glioma (LGG), Kidney Chromophobe (KICH), Sarcoma (SARC), osteosarcoma (OS), Pancreatic Cancer (PAAD), and Thyroid Cancer (THCA). In addition, elevated PUS7 expression was linked to advanced TNM across multiple malignancies, including ACC, BLCA, KIRP, LIHC and PAAD. The function enrichment analysis revealed that PUS7 participates in E2F targets, G2M checkpoint, ribosome biogenesis, and rRNA metabolic process. Moreover, PUS7 is also a reliable biomarker and a potential therapeutic target for osteosarcoma.

CONCLUSIONS

In summary, PUS7 is a putative pan-cancer biomarker that reliably forecasts cancer patients' prognosis. In addition, this enzyme regulates the cell cycle, ribosome biogenesis, and rRNA metabolism. Most importantly, PUS7 possibly regulates osteosarcoma initiation and progression.

RANK/RANKL axis promotes migration, invasion, and metastasis of osteosarcoma via activating NF-κB pathway

Osteosarcoma (OS) is one of the most prevalent primary bone tumors with a high degree of metastasis and poor prognosis. Epithelial-to-mesenchymal transition (EMT) is a cellular mechanism that contributes to the invasion and metastasis of cancer cells, and OS cells have been reported to exhibit EMT-like characteristics. Our previous studies have shown that the interaction between tumor necrosis factor superfamily member 11 (TNFRSF11A; also known as RANK) and its ligand TNFSF11 (also known as RANKL) promotes the EMT process in breast cancer cells. However, whether the interaction between RANK and RANKL enhances aggressive behavior by inducing EMT in OS cells has not yet been elucidated. In this study, we showed that the interaction between RANK and RANKL increased the migration, invasion, and metastasis of OS cells by promoting EMT. Importantly, we clarified that the RANK/RANKL axis induces EMT by activating the nuclear factor-kappa B (NF-κB) pathway. Furthermore, the NF-κB inhibitor dimethyl fumarate (DMF) suppressed migration, invasion, and EMT in OS cells. Our results suggest that the RANK/RANKL axis may serve as a potential tumor marker and promising therapeutic target for OS metastasis. Furthermore, DMF may have clinical applications in the treatment of lung metastasis in patients with OS.

Angio-Long Noncoding RNA MALAT1 (rs3200401) and MIAT (rs1061540) Gene Variants in Ovarian Cancer

The genotyping of long non-coding RNA (lncRNA)-related single-nucleotide polymorphisms (SNPs) could be associated with cancer risk and/or progression. This study aimed to analyze the angiogenesis-related lncRNAs MALAT1 (rs3200401) and MIAT (rs1061540) variants in patients with ovarian cancer (OC) using "Real-Time allelic discrimination polymerase chain reaction" in 182 formalin-fixed paraffin-embedded (FFPE) samples of benign, borderline, and primary malignant ovarian tissues. Differences in the genotype frequencies between low-grade ovarian epithelial tumors (benign/borderline) and malignant tumors and between high-grade malignant epithelial tumors and malignant epithelial tumors other than high-grade serous carcinomas were compared. Odds ratios (ORs)/95% confidence intervals were calculated as measures of the association strength. Additionally, associations of the genotypes with the available pathological data were analyzed. The heterozygosity of MALAT1 rs3200401 was the most common genotype (47.8%), followed by C/C (36.3%). Comparing the study groups, no significant differences were observed regarding this variant. In contrast, the malignant epithelial tumors had a higher frequency of the MIAT rs1061540 C/C genotype compared to the low-grade epithelial tumor cohorts (56.7% vs. 37.6, p = 0.031). The same genotype was significantly higher in high-grade serous carcinoma than its counterparts (69.4% vs. 43.8%, p = 0.038). Multivariate Cox regression analysis showed that the age at diagnosis was significantly associated with the risk of OC development. In contrast, the MIAT T/T genotype was associated with a low risk of malignant epithelial tumors under the homozygote comparison model (OR = 0.37 (0.16-0.83), p = 0.017). Also, MIAT T allele carriers were less likely to develop high-grade serous carcinoma under heterozygote (CT vs. CC; OR = 0.33 (0.12-0.88), p = 0.027) and homozygote (TT vs. CC; OR = 0.26 (0.07-0.90), p = 0.034) comparison models. In conclusion, our data provide novel evidence for a potential association between the lncRNA MIAT rs1061540 and the malignant condition of ovarian cancer, suggesting the involvement of such lncRNAs in OC development.

TGF-β and BMP signaling are associated with the transformation of glioblastoma to gliosarcoma and then osteosarcoma

Background

Gliosarcoma, an isocitrate dehydrogenase wildtype (IDH-WT) variant of glioblastoma, is defined by clonal biphasic differentiation into gliomatous and sarcomatous components. While the transformation from a glioblastoma to gliosarcoma is uncommon, the subsequent transformation to osteosarcoma is rare but may provide additional insights into the biology of these typically distinct cancers. We observed a patient initially diagnosed with glioblastoma, that differentiated into gliosarcoma at recurrence, and further evolved to osteosarcoma at the second relapse. Our objective was to characterize the molecular mechanisms of tumor progression associated with this phenotypic transformation.

Methods

Tumor samples were collected at all 3 stages of disease and RNA sequencing was performed to capture their transcriptomic profiles. Sequential clonal evolution was confirmed by the maintenance of an identical PTEN mutation throughout the tumor differentiation using the TSO500 gene panel. Publicly available datasets and the Nanostring nCounter technology were used to validate the results.

Results

The glioblastoma tumor from this patient possessed mixed features of all 3 TCGA-defined transcriptomic subtypes of an IDH-WT glioblastoma and a proportion of osteosarcoma signatures were upregulated in the original tumor. Analysis showed that enhanced transforming growth factor-β (TGF-β) and bone morphogenic protein signaling was associated with tumor transformation. Regulatory network analysis revealed that TGF-β family signaling committed the lineage tumor to osteogenesis by stimulating the expression of runt-related transcription factor 2 (RUNX2), a master regulator of bone formation.

Conclusions

This unusual clinical case provided an opportunity to explore the modulators of longitudinal sarcomatous transformation, potentially uncovering markers indicating predisposition to this change and identification of novel therapeutic targets.

Adapted Deep Ensemble Learning-Based Voting Classifier for Osteosarcoma Cancer Classification

The study utilizes osteosarcoma hematoxylin and the Eosin-stained image dataset, which is unevenly dispersed, and it raises concerns about the potential impact on the overall performance and reliability of any analyses or models derived from the dataset. In this study, a deep-learning-based convolution neural network (CNN) and adapted heterogeneous ensemble-learning-based voting classifier have been proposed to classify osteosarcoma. The proposed methods can also resolve the issue and develop unbiased learning models by introducing an evenly distributed training dataset. Data augmentation is employed to boost the generalization abilities. Six different pre-trained CNN models, namely MobileNetV1, Mo-bileNetV2, ResNetV250, InceptionV2, EfficientNetV2B0, and NasNetMobile, are applied and evaluated in frozen and fine-tuned-based phases. In addition, a novel CNN model and adapted heterogeneous ensemble-learning-based voting classifier developed from the proposed CNN model, fine-tuned NasNetMobile model, and fine-tuned Efficient-NetV2B0 model are also introduced to classify osteosarcoma. The proposed CNN model outperforms other pre-trained models. The Kappa score obtained from the proposed CNN model is 93.09%. Notably, the proposed voting classifier attains the highest Kappa score of 96.50% and outperforms all other models. The findings of this study have practical implications in telemedicine, mobile healthcare systems, and as a supportive tool for medical professionals.

Construction and validation of an oxidative-stress-related risk model for predicting the prognosis of osteosarcoma

BACKGROUND

Osteosarcoma is the most common bone malignancy in teenagers, and warrants effective measures for diagnosis and prognosis. Oxidative stress (OS) is the key driver of several cancers and other diseases.

METHODS

The TARGET-osteosarcoma database was employed as the training cohort and GSE21257 and GSE39055 was applied for external validation. The patients were classified into the high- and low-risk groups based on the median risk score of each sample. ESTIMATE and CIBERSORT were applied for the evaluation of tumor microenvironment immune infiltration. GSE162454 of single-cell sequencing was employed for analyzing OS-related genes.

RESULTS

Based on the gene expression and clinical data of 86 osteosarcoma patients in the TARGET database, we identified eight OS-related genes, including MAP3K5, G6PD, HMOX1, ATF4, ACADVL, MAPK1, MAPK10, and INS. In both the training and validation sets, the overall survival of patients in the high-risk group was significantly worse than that in the low-risk group. The ESTIMATE algorithm revealed that patients in the high-risk group had higher tumor purity but lower immune score and stromal score. In addition, the CIBERSORT algorithm showed that the M0 and M2 macrophages were the predominant infiltrating cells in osteosarcoma. Based on the expression analysis of immune checkpoint, CD274(PDL1), CXCL12, BTN3A1, LAG3, and IL10 were identified as potential immune therapy targets. Analysis of the single cell sequencing data also revealed the expression patterns of OS-related genes in different cell types.

CONCLUSIONS

An OS-related prognostic model can accurately provide the prognosis of osteosarcoma patients, and may help identify suitable candidates for immunotherapy.

ac4C acetylation regulates mRNA stability and translation efficiency in osteosarcoma

Objective

N4-acetylcytidine (ac4C) acetylation can promote target gene expression through improved mRNA stability. To explore the role of ac4C acetylation in osteosarcoma, U2OS and MG63 cell lines were treated with the N-acetyltransferase 10 (NAT10) inhibitor Remodelin. Reverse transcription-polymerase chain reaction (RT-PCR) and Western blot were used to test the gene and protein expression efficiency.

Methods

The proliferation rate of osteosarcoma cells was measured by a cell counting kit-8 (CCK8) assay. The cell cycle and apoptosis were analyzed by flow cytometry. The invasiveness of osteosarcoma cells was detected by a transwell invasion assay. The ac4C acetylation of target genes was screened by acetylated RNA immunoprecipitation and sequencing (acRIP-seq).

Results

We found that when osteosarcoma cells were treated with Remodelin at the optimal concentration, their NAT10 expression and the cell proliferation was inhibited, the cells in the G1 phase increased (P < 0.05) but those in the S phase decreased, the apoptotic cells in the early and late stages increased, and the cells invasiveness decreased (P < 0.05).

Conclusions

The farnesyltransferase subunit beta gene (FNTB) was identified by acRIP-seq as one of the target genes of ac4C acetylation and was further verified by RT-PCR and Western blot analyses. Remodelin was demonstrated to reduce the stability and protein translation efficiency of target gene mRNA in osteosarcoma cells. In conclusion, inhibition of ac4C acetylation in osteosarcoma can block proliferation and metastasis as well as promote apoptosis and cell cycle arrest. Ac4C acetylation contributes to the stability and protein translation efficiency of the downstream target gene mRNA.

Chondrosarcoma Resistance to Radiation Therapy: Origins and Potential Therapeutic Solutions

Chondrosarcoma is a malignant cartilaginous tumor that is particularly chemoresistant and radioresistant to X-rays. The first line of treatment is surgery, though this is almost impossible in some specific locations. Such resistances can be explained by the particular composition of the tumor, which develops within a dense cartilaginous matrix, producing a resistant area where the oxygen tension is very low. This microenvironment forces the cells to adapt and dedifferentiate into cancer stem cells, which are described to be more resistant to conventional treatments. One of the main avenues considered to treat this type of tumor is hadrontherapy, in particular for its ballistic properties but also its greater biological effectiveness against tumor cells. In this review, we describe the different forms of chondrosarcoma resistance and how hadrontherapy, combined with other treatments involving targeted inhibitors, could help to better treat high-grade chondrosarcoma.

HERVs and Cancer-A Comprehensive Review of the Relationship of Human Endogenous Retroviruses and Human Cancers

Genomic instability and genetic mutations can lead to exhibition of several cancer hallmarks in affected cells such as sustained proliferative signaling, evasion of growth suppression, activated invasion, deregulation of cellular energetics, and avoidance of immune destruction. Similar biological changes have been observed to be a result of pathogenic viruses and, in some cases, have been linked to virus-induced cancers. Human endogenous retroviruses (HERVs), once external pathogens, now occupy more than 8% of the human genome, representing the merge of genomic and external factors. In this review, we outline all reported effects of HERVs on cancer development and discuss the HERV targets most suitable for cancer treatments as well as ongoing clinical trials for HERV-targeting drugs. We reviewed all currently available reports of the effects of HERVs on human cancers including solid tumors, lymphomas, and leukemias. Our review highlights the central roles of HERV genes, such as gag, env, pol, np9, and rec in immune regulation, checkpoint blockade, cell differentiation, cell fusion, proliferation, metastasis, and cell transformation. In addition, we summarize the involvement of HERV long terminal repeat (LTR) regions in transcriptional regulation, creation of fusion proteins, expression of long non-coding RNAs (lncRNAs), and promotion of genome instability through recombination.

Comprehensive analysis of potential cellular communication networks in advanced osteosarcoma using single-cell RNA sequencing data

Osteosarcoma (OS) is a common bone cancer in children and adolescents, and metastasis and recurrence are the major causes of poor treatment outcomes. A better understanding of the tumor microenvironment is required to develop an effective treatment for OS. In this paper, a single-cell RNA sequencing dataset was taken to a systematic genetic analysis, and potential signaling pathways linked with osteosarcoma development were explored. Our findings revealed 25 clusters across 11 osteosarcoma tissues, with 11 cell types including “Chondroblastic cells”, “Osteoblastic cells”, “Myeloid cells”, “Pericytes”, “Fibroblasts”, “Proliferating osteoblastic cells”, “Osteoclasts”, “TILs”, “Endothelial cells”, “Mesenchymal stem cells”, and “Myoblasts”. The results of Cell communication analysis showed 17 potential cellular communication networks including “COLLAGEN signaling pathway network”, “CD99 signaling pathway network”, “PTN signaling pathway network”, “MIF signaling pathway network”, “SPP1 signaling pathway network”, “FN1 signaling pathway network”, “LAMININ signaling pathway network”, “FGF signaling pathway network”, “VEGF signaling pathway network”, “GALECTIN signaling pathway network”, “PERIOSTIN signaling pathway network”, “VISFATIN signaling pathway network”, “ITGB2 signaling pathway network”, “NOTCH signaling pathway network”, “IGF signaling pathway network”, “VWF signaling pathway network”, “PDGF signaling pathway network”. This research may provide novel insights into the pathophysiology of OS’s molecular processes.

Repeat expression is linked to patient survival and exhibits single nucleotide variation in pancreatic cancer revealing LTR70:r.879A>G

Despite an overwhelming number of cancer literature reporting the links between patient survival and the expression levels of genes or mutations/single nucleotide variations (SNVs) on them, there is only limited information on repeat elements, which make at least half the human genome. Here, we analysed RNA-seq data obtained from primary pancreatic cancer tissues of 51 patients and revealed that two transposons, HERVI-int and X6A_LINE, showed an upregulation trend in the patients who lived shorter, along with 56 other potential repeats which were linked to survival. We also detected expressed single nucleotide variations (SNVs) on repeats, among which LTR70:r.879A>G stands out with the effect of its presence on this particular repeat's expression levels and a significant link to overall patient survival. Interestingly, the expression of LTR70:r.879A>G correlated with different cancer genes in comparison to its reference version highlighting the involvement of BRAF and Fumerate Hydratase with this expressed SNV. This is one of the first studies revealing possible links between repeat expression and survival in cancer and it warrants further research in this avenue.

Single-Cell Profiling of Tumor Microenvironment Heterogeneity in Osteosarcoma Identifies a Highly Invasive Subcluster for Predicting Prognosis

Osteosarcoma is the most common malignant bone tumor in adolescents, and metastasis is the key reason for treatment failure and poor prognosis. Once metastasis occurs, the 5-year survival rate is only approximately 20%, and assessing and predicting the risk of osteosarcoma metastasis are still difficult tasks. In this study, cellular communication between tumor cells and nontumor cells was identified through comprehensive analysis of osteosarcoma single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data, illustrating the complex regulatory network in the osteosarcoma microenvironment. In line with the heterogeneity of osteosarcoma, we found subpopulations of osteosarcoma cells that highly expressed COL6A1, COL6A3 and MIF and were closely associated with lung metastasis. Then, BCDEG, a reliable risk regression model that could accurately assess the metastasis risk and prognosis of patients, was established, providing a new strategy for the diagnosis and treatment of osteosarcoma.

Pericentromeric Non-Coding DNA Transcription Is Associated with Niche Impairment in Patients with Ineffective or Partially Effective Multiple Myeloma Treatment

Mesenchymal stromal cells (MSC) ‘educated’ by tumor cells are an essential component of the multiple myeloma (MM) tumor microenvironment (TME) involved in tumor progression. Transcription of tandemly repeated (TR) non-coding DNA is often activated in many tumors and is required for tumor progression and cancer cells genome reorganization. The aim of the work was to study functional properties including the TR DNA transcription profile of MSC from the hematopoietic niche of treated MM patients. Healthy donors (HD) and patients after bortezomib-based treatment (with partial or complete response, PoCR, and non-responders, NR) were enrolled in the study. Their trephine biopsies were examined histologically to evaluate the hematopoietic niche. MSC cultures obtained from the biopsies were used for evaluation of the proliferation rate, osteogenic differentiation, presence of tumor MSC markers, resistance to bortezomib, and pericentromeric TR DNA transcription level. The MSC ‘education’ by multiple myeloma cells was mimicked in co-culture experiments with or without bortezomib. The TR DNA transcription profile was accessed. The histological examination revealed the persistence of the tumor microenvironment (especially of the vasculature) in treated patients. In co-culture experiments, MSC of bortezomib-treated patients were more resistant to bortezomib and protected cancer MM cells of the RPMI8226 cell line more effectively than HD-MSC did. The MSC obtained from PoCR and NR samples differed in their functional properties (proliferation capacity, osteogenic potential, and cancer-associated fibroblasts markers). Transcriptome analysis revealed activation of the TR transcription in cells of non-hematopoietic origin from NR patients’ bone marrow. The pericentromeric TR DNA of HS2/HS3 families was among the most upregulated in stromal MSC but not in cancer cells. The highest level of transcription was observed in NR-MSC. Transcription of HS2/HS3 was not detected in healthy donors MSC unless they were co-cultured with MM cancer cells and acquired cancer-associated phenotype. Treatment with TNFα downregulated HS2/HS3 transcription in MSC and upregulated in MM cells. Our results suggest that the hematopoietic niche retains the cancer-associated phenotype after treatment. Pericentromeric non-coding DNA transcription is associated with the MSC cancer-associated phenotype in patients with ineffective or partially effective multiple myeloma treatment.

FNC: An Advanced Anticancer Therapeutic or Just an Underdog?

Azvudine (FNC) is a novel cytidine analogue that has both antiviral and anticancer activities. This minireview focuses on its underlying molecular mechanisms of suppressing viral life cycle and cancer cell growth and discusses applications of this nucleoside drug for advanced therapy of tumors and malignant blood diseases. FNC inhibits positive-stand RNA viruses, like HCV, EV, SARS-COV-2, HBV, and retroviruses, including HIV, by suppressing their RNA-dependent polymerase enzymes. It may also inhibit such enzyme (reverse transcriptase) in the human retrotransposons, including human endogenous retroviruses (HERVs). As the activation of retrotransposons can be the major factor of ongoing cancer genome instability and consequently higher aggressiveness of tumors, FNC has a potential to increase the efficacy of multiple anticancer therapies. Furthermore, FNC also showed other aspects of anticancer activity by inhibiting adhesion, migration, invasion, and proliferation of malignant cells. It was also reported to be involved in cell cycle arrest and apoptosis, thereby inhibiting the progression of cancer through different pathways. To the date, the grounds of FNC effects on cancer cells are not fully understood and hence additional studies are needed for better understanding molecular mechanisms of its anticancer activities to support its medical use in oncology.

Hsp90 regulates the tumorigenic function of tyrosine protein kinase in osteosarcoma

Despite recent advances in diagnosis and treatment, osteosarcoma remains as the most common bone cancer in children and is associated with poor prognosis. Growing evidence has supported dysregulation of threonine and tyrosine protein kinase (TTK) expression as a hallmark of multiple cancers, however, its function in osteosarcoma remains to be elucidated. In the present study, we found that TTK was frequently overexpressed in osteosarcoma and associated with increased tumour growth and progression. Moreover, using both in vitro and in vivo assays, we provided evidence that TTK level was regulated by a molecular chaperone, heat shock protein 90 (Hsp90). Hsp90 directly interacted with TTK and prevents proteasome-dependent TTK degradation, leading to the accumulation of TTK in osteosarcoma cells. Elevated TTK promoted cancer cell proliferation and survival by activating cell-cycle progression and inhibiting apoptosis. Consistently, depletion of TTK by Hsp90 inhibition induced cell-cycle arrest, generated aneuploidy and eventually resulted in apoptotic cancer cell death. Together, our study revealed an important Hsp90-TTK regulatory axis in osteosarcoma cells to promote cancer cell growth and survival. These findings expand our knowledge on osteosarcoma pathogenesis and offer novel therapeutic options for clinical practice.

Differential enrichment of H3K9me3 at annotated satellite DNA repeats in human cell lines and during fetal development in mouse

BACKGROUND

Trimethylation of histone H3 on lysine 9 (H3K9me3) at satellite DNA sequences has been primarily studied at (peri)centromeric regions, where its level shows differences associated with various processes such as development and malignant transformation. However, the dynamics of H3K9me3 at distal satellite DNA repeats has not been thoroughly investigated.

RESULTS

We exploit the sets of publicly available data derived from chromatin immunoprecipitation combined with massively parallel DNA sequencing (ChIP-Seq), produced by the The Encyclopedia of DNA Elements (ENCODE) project, to analyze H3K9me3 at assembled satellite DNA repeats in genomes of human cell lines and during mouse fetal development. We show that annotated satellite elements are generally enriched for H3K9me3, but its level in cancer cell lines is on average lower than in normal cell lines. We find 407 satellite DNA instances with differential H3K9me3 enrichment between cancer and normal cells including a large 115-kb cluster of GSATII elements on chromosome 12. Differentially enriched regions are not limited to satellite DNA instances, but instead encompass a wider region of flanking sequences. We found no correlation between the levels of H3K9me3 and noncoding RNA at corresponding satellite DNA loci. The analysis of data derived from multiple tissues identified 864 instances of satellite DNA sequences in the mouse reference genome that are differentially enriched between fetal developmental stages.

CONCLUSIONS

Our study reveals significant differences in H3K9me3 level at a subset of satellite repeats between biological states and as such contributes to understanding of the role of satellite DNA repeats in epigenetic regulation during development and carcinogenesis.

Identification of abnormally high expression of POGZ as a new biomarker associated with a poor prognosis in osteosarcoma

Osteosarcoma (OS) is the most prevalent malignant bone tumor in children and young adults. There is an urgent need for a novel biomarker related to the prognosis of OS. We performed a meta-analysis incorporating six independent datasets and performed a survival analysis with one independent dataset GSE21257 in the GEO database for gene screening. The results revealed that one potential biomarker related to OS survival, POGZ was the most significantly upregulated gene. We also verified that the POGZ was overexpressed in clinical samples. The survival analysis revealed that POGZ is associated with a poor prognosis in OS. Moreover, flow cytometry analysis of isolated OS cells demonstrated that OS cells were arrested in the G1 phase after POGZ knockdown. The RNA-seq results indicated that POGZ was co-expressed with CCNE1 and CCNB1. Pathway analysis showed that genes associated with high expression levels of POGZ were related to the cell cycle pathway. A cell model was constructed to detect the effects of POGZ. After POGZ knockdown, OS cell proliferation, invasion and migration were all decreased. Therefore, POGZ is an important gene for evaluating the prognosis of OS patients and is a potential therapeutic target.

Immunohistochemical Characterisation of GLUT1, MMP3 and NRF2 in Osteosarcoma

Osteosarcoma (OSA) is an aggressive bone malignancy. Unlike many other malignancies, OSA outcomes have not improved in recent decades. One challenge to the development of better diagnostic and therapeutic methods for OSA has been the lack of well characterized experimental model systems. Spontaneous OSA in dogs provides a good model for the disease seen in people and also remains an important veterinary clinical challenge. We recently used RNA sequencing and qRT-PCR to provide a detailed molecular characterization of OSA relative to non-malignant bone in dogs. We identified differential mRNA expression of the solute carrier family 2 member 1 (SLC2A1/GLUT1), matrix metallopeptidase 3 (MMP3) and nuclear factor erythroid 2–related factor 2 (NFE2L2/NRF2) genes in canine OSA tissue in comparison to paired non-tumor tissue. Our present work characterizes protein expression of GLUT1, MMP3 and NRF2 using immunohistochemistry. As these proteins affect key processes such as Wnt activation, heme biosynthesis, glucose transport, understanding their expression and the enriched pathways and gene ontologies enables us to further understand the potential molecular pathways and mechanisms involved in OSA. This study further supports spontaneous OSA in dogs as a model system to inform the development of new methods to diagnose and treat OSA in both dogs and people.

Endogenous retrovirus envelope as a tumor-associated immunotherapeutic target in murine osteosarcoma

Osteosarcoma remains one of the deadliest cancers in pediatrics and young adults. We administered two types of immunotherapies, oncolytic virotherapy and immune checkpoint inhibition, to two murine osteosarcoma models and observed divergent results. Mice bearing F420 showed no response, whereas those with K7M2 showed prolonged survival in response to combination therapy. K7M2 had higher expression of immune-related genes and higher baseline immune cell infiltrates, but there were no significant differences in tumor mutational burden or predicted MHC class I binding of nonsynonymous mutations. Instead, we found several mouse endogenous retrovirus sequences highly expressed in K7M2 compared with F420. T cell tetramer staining for one of them, gp70, was detected in mice with K7M2 but not F420, suggesting that endogenous retrovirus proteins are targets for the anti-tumor immune reaction. Given prior observations of endogenous retrovirus expression in human osteosarcomas, our findings may be translatable to human disease.

Expression Quantitative Trait Loci (eQTLs) Associated with Retrotransposons Demonstrate their Modulatory Effect on the Transcriptome

Transposable elements (TEs) are repetitive elements that belong to a variety of functional classes and have an important role in shaping genome evolution. Around 50% of the human genome contains TEs, and they have been termed the "dark matter" of the genome because relatively little is known about their function. While TEs have been shown to participate in aberrant gene regulation and the pathogenesis of diseases, only a few studies have explored the systemic effect of TEs on gene expression. In the present study, we analysed whole genome sequences and blood whole transcriptome data from 570 individuals within the Parkinson's Progressive Markers Initiative (PPMI) cohort to identify expression quantitative trait loci (eQTL) regulating genome-wide gene expression associated with TEs. We identified 2132 reference TEs that were polymorphic for their presence or absence in our study cohort. The presence or absence of the TE element could change the expression of the gene or gene clusters from zero to tens of thousands of copies of RNA. The main finding is that many TEs possess very strong regulatory effects, and they have the potential to modulate large genetic networks with hundreds of target genes over the genome. We illustrate the plethora of regulatory mechanisms using examples of their action at the HLA gene cluster and data showing different TEs' convergence to modulate WFS1 gene expression. In conclusion, the presence or absence of polymorphisms of TEs has an eminent genome-wide regulatory function with large effect size at the level of the whole transcriptome. The role of TEs in explaining, in part, the missing heritability for complex traits is convincing and should be considered.

Identification of Gene as Predictive Biomarkers for the Occurrence and Recurrence of Osteosarcoma

Purpose

Osteosarcoma is the most common malignant bone cancer affecting adolescents and young adults. This study aimed to screen potential diagnostic and therapeutic markers for osteosarcoma.

Methods

Differential expression analysis between osteosarcoma and control was performed in GSE99671, the differentially expressed genes (DEGs) were subjected to co-expression analysis. Enrichment analysis was employed to identify the biological functions and KEGG signaling pathways of module genes. In addition, a differential analysis was also performed between recurrent and non-recurrent osteosarcoma samples in GSE39055, and enrichment analysis was performed for DEGs. Further, Kaplan–Meier curve analysis was performed on the module genes, and receiver operating characteristic (ROC) curve was drawn. Comparison of the module with the highest correlation to osteosarcoma identified key genes. Cox regression model was utilized to identify the predictive ability of key genes for the prognosis of osteosarcoma.

Results

A total of 13 co-expression modules were identified from 4871 DEGs of GSE99671, module 1 had the highest positive correlation with osteosarcoma. Module genes were mainly enriched in autophagy and macrophage migration functions. A total of 1126 DEGs were obtained from GSE39055, significantly involved in neutrophil mediated immunity. Screening of genes with area under the ROC curve (AUC) values greater than 0.73 in both GSE99671 and GSE39055 identified 5 key genes when compared with genes from module 1. The nomogram results showed that ATF5, CHCHD8, ENOPH1, and LOC286367 might predict 5-year or 8-year survival time of osteosarcoma patients. The Cox model results confirmed that the signals of ATF5, CHCHD8, and LOC286367 were robust, and it may be used in the diagnosis, treatment, and prognosis of osteosarcoma.

Conclusion

We found that ATF5, CHCHD8, and LOC286367 can effectively identify osteosarcoma tumorigenesis and even recurrence status. This is helpful for early diagnosis and treatment, improving the clinical treatment of patients with osteosarcoma.

Human Endogenous Retroviruses in Glioblastoma Multiforme

Glioblastoma multiforme (GBM) is the most aggressive and deadly brain tumor. It is primarily diagnosed in the elderly and has a 5-year survival rate of less than 6% even with the most aggressive therapies. The lack of biomarkers has made the development of immunotherapy for GBM challenging. Human endogenous retroviruses (HERVs) are a group of viruses with long terminal repeat (LTR) elements, which are believed to be relics from ancient viral infections. Recent studies have found that those repetitive elements play important roles in regulating various biological processes. The differentially expressed LTR elements from HERVs are potential biomarkers for immunotherapy to treat GBM. However, the understanding of the LTR element expression in GBM is greatly lacking.

METHODS

We obtained 1077.4 GB of sequencing data from public databases. These data were generated from 111 GBM tissue studies, 30 GBM cell lines studies, and 45 normal brain tissues studies. We analyzed repetitive elements that were differentially expressed in GBM and normal brain samples.

RESULTS

We found that 48 LTR elements were differentially expressed (p-value < 0.05) between GBM and normal brain tissues, of which 46 were HERV elements. Among these 46 elements, 34 significantly changed HERVs belong to the ERV1 superfamily. Furthermore, 43 out of the 46 differentially expressed HERV elements were upregulated.

CONCLUSION

Our results indicate significant differential expression of many HERV LTR elements in GBM and normal brain tissues. Expression levels of these elements could be developed as biomarkers for GBM treatments.

At the dawn of the transcriptomic medicine

Progress in genomic analytical technologies has improved our possibilities to obtain information regarding DNA, RNA, and their dynamic changes that occur over time or in response to specific challenges. This information describes the blueprint for cells, tissues, and organisms and has fundamental importance for all living organisms. This review focuses on the technological challenges to analyze the transcriptome and what is the impact of transcriptomics on precision medicine. The transcriptome is a term that covers all RNA present in cells and a substantial part of it will never be translated into protein but is nevertheless functional in determining cell phenotype. Recent developments in transcriptomics have challenged the fundamentals of the central dogma of biology by providing evidence of pervasive transcription of the genome. Such massive transcriptional activity is challenging the definition of a gene and especially the term "pseudogene" that has now been demonstrated in many examples to be both transcribed and translated. We also review the common sources of biomaterials for transcriptomics and justify the suitability of whole blood RNA as the current optimal analyte for clinical transcriptomics. At the end of the review, a brief overview of the clinical implications of transcriptomics in clinical trial design and clinical diagnosis is given. Finally, we introduce the transcriptome as a target for modern drug development as a tool for extending our capacity for precision medicine in multiple diseases.

Extracellular vesicle-associated repetitive element DNAs as candidate osteosarcoma biomarkers

Osteosarcoma (OS) is the most common malignant bone tumor in children and young adults. Despite that high-risk factors have been identified, no test for early detection is available. This study aimed to identify circulating nucleic acid sequences associated with serum extracellular vesicle (EV) preparations at the time of OS diagnosis, as a step towards an OS early detection assay. Sequencing of small nucleic acids extracted from serum EV preparations revealed increased representation of diverse repetitive element sequences in OS patient versus control sera. Analysis of a validation cohort using qPCR of PEG-precipitated EV preparations revealed the over-representation of HSATI, HSATII, LINE1-P1, and Charlie 3 at the DNA but not RNA level, with receiver operating characteristic (ROC) area under the curve (AUC) ≥ 0.90. HSATI and HSATII DNAs co-purified with EVs prepared by precipitation and size exclusion chromatography but not by exosome immunocapture, indicative of packaging in a non-exosomal complex. The consistent over-representation of EV-associated repetitive element DNA sequences suggests their potential utility as biomarkers for OS and perhaps other cancers.

HSATII RNA is induced via a noncanonical ATM-regulated DNA damage response pathway and promotes tumor cell proliferation and movement

Pericentromeric human satellite II (HSATII) repeats are normally silent but can be actively transcribed in tumor cells, where increased HSATII copy number is associated with a poor prognosis in colon cancer, and in human cytomegalovirus (HCMV)-infected fibroblasts, where the RNA facilitates viral replication. Here, we report that HCMV infection or treatment of ARPE-19 diploid epithelial cells with DNA-damaging agents, etoposide or zeocin, induces HSATII RNA expression, and a kinase-independent function of ATM is required for the induction. Additionally, various breast cancer cell lines growing in adherent, two-dimensional cell culture express HSATII RNA at different levels, and levels are markedly increased when cells are infected with HCMV or treated with zeocin. High levels of HSATII RNA expression correlate with enhanced migration of breast cancer cells, and knockdown of HSATII RNA reduces cell migration and the rate of cell proliferation. Our investigation links high expression of HSATII RNA to the DNA damage response, centered on a noncanonical function of ATM, and demonstrates a role for the satellite RNA in tumor cell proliferation and movement.

Alternative splicing of leptin receptor overlapping transcript in osteosarcoma

IMPACT STATEMENT

Osteosarcoma (OS, also known as osteogenic sarcoma) is the most common primary malignancy of bone in children and adolescents. The molecular mechanisms of OS are extremely complicated and its molecular mediators remain to be elucidated. We sequenced total RNA from 18 OS bone samples (paired normal-tumor biopsies). We found statistically significant (FDR <0.05) 26 differentially expressed transcript variants of LEPROT gene with different expressions in normal and tumor samples. These findings contribute to the understanding of molecular mechanisms of OS development and provide encouragement to pursue further research.

Biomarker potential of repetitive-element transcriptome in lung cancer

Since repetitive elements (REs) account for nearly 53% of the human genome, profiling its transcription after an oncogenic change might help in the search for new biomarkers. Lung cancer was selected as target since it is the most frequent cause of cancer death. A bioinformatic workflow based on well-established bioinformatic tools (such as RepEnrich, RepBase, SAMTools, edgeR and DESeq2) has been developed to identify differentially expressed RNAs from REs. It was trained and tested with public RNA-seq data from matched sequencing of tumour and healthy lung tissues from the same patient to reveal differential expression within the RE transcriptome. Healthy lung tissues express a specific set of REs whose expression, after an oncogenic process, is strictly and specifically changed. Discrete sets of differentially expressed REs were found for lung adenocarcinoma, for small-cell lung cancer, and for both cancers. Differential expression affects more HERV-than LINE-derived REs and seems biased towards down-regulation in cancer cells. REs behaving consistently in all patients were tested in a different patient cohort to validate the proposed biomarkers. Down-regulation of AluYg6 and LTR18B was confirmed as potential lung cancer biomarkers, while up-regulation of HERVK11D-Int is specific for lung adenocarcinoma and up-regulation of UCON88 is specific for small cell lung cancer. Hence, the study of RE transcriptome might be considered another research target in cancer, making REs a promising source of lung cancer biomarkers.

Hormonal regulation of visfatin gene in avian Leghorn male hepatoma (LMH) cells

Visfain has been extensively studied in mammals and has been shown to play an important role in obesity and insulin resistance. However, there is a paucity of information on visfatin regulation in non-mammalian species. After characterization of chicken visfatin gene, we undertook this study to determine its hormonal regulation in avian (non-mammalian) liver cells. Addition of 5 ng/mL TNFα, 100 ng/mL leptin, 1, 3, 10 or 100 ng/mL T3 for 24 h upregulated visfatin gene expression by 1.2, 1.8, 1.95, 1.75, 1.80, and 2.45 folds (P < .05), respectively, compared to untreated LMH cells. Administration of 10 ng/mL of orexin A significantly down regulated visfatin gene expression by 1.35 folds compared to control cells. In contrast, treatment with IL-6 or orexin B for 24 h did not influence visfatin mRNA abundance. These pro-inflammatory cytokines and obesity-related hormones modulate the expression of CRP, INSIG2, and nuclear orphan receptors. Hepatic CRP gene expression was significantly upregulated by IL-6, TNFα, orexin B, and T3 and down regulated by leptin and orexin A. LXR mRNA abundances were increased by orexin A, decreased by orexin B, and T3, and did not affected by IL6, TNFα, or leptin. The expression of FXR gene was induced by IL-6, leptin, and T3, but it was not influenced by TNFα, orexin A or B. CXR gene expression was up regulated by TNFα, leptin, orexin B, and T3, down regulated by 5 ng/mL orexin A, and did not affected by IL-6. INSIG2 mRNA levels were increased by TNFα (5 ng/mL), leptin (100 ng/mL), and T3 (1, 3, 10, and 100 ng/mL), decreased by orexin A, and remained unchanged with IL-6 or orexin B treatment. Together, this is the first report showing hormonal regulation of visfatin in avian hepatocyte cells and suggesting a potential role of CRP, INSIG2, and nuclear orphan receptor LXR, FXR, and CXR in mediating these hormonal effects.

Close to the Bedside: A Systematic Review of Endogenous Retroviruses and Their Impact in Oncology

BACKGROUND

Human endogenous retroviruses (HERVs) are genetic elements in the human genome, which resulted from ancient retroviral germline infections. HERVs have strong transcriptional promoters and enhancers that affect a cell's transcriptome. They also encode proteins that can exert effects in human cells. This review examines how our increased understanding of HERVs have led to their potential use as biomarkers and immunologic targets.

MATERIAL AND METHODS

PubMed/Medline, Embase, Web of Science, and Cochrane databases were used in a systematic search to identify all articles studying the potential impact of HERVs on surgical diseases. The search included studies that involved clinical patient samples in diseases including cancer, inflammatory conditions, and autoimmune disease. Articles focused on conditions not routinely managed by surgeons were excluded.

RESULTS

Eighty six articles met inclusion and quality criteria for this review and were included. Breast cancer and melanoma have robust evidence regarding the use of HERVs as potential tumor markers and immunologic targets. Reported evidence of the activity of HERVs in colorectal cancer, pancreatic cancer, hepatocellular cancer, prostate and ovarian cancer, germ cell tumors as well as idiopathic pulmonary hypertension, and the inflammatory response in burns was also reviewed.

CONCLUSIONS

Increasingly convincing evidence indicates that HERVs may play a role in solid organ malignancy and present important biomarkers or immunologic targets in multiple cancers. Innovative investigation of HERVs is a valuable focus of translational research and can deepen our understanding of cellular physiology and the effects of endogenous retroviruses on human biology. As strategies for treatment continue to focus on genome-based interventions, understanding the impact of endogenous retroviruses on human disease will be critical.

ERVmap analysis reveals genome-wide transcription of human endogenous retroviruses

Endogenous retrovirus (ERV) sequences make up a large fraction of our genome, yet little is understood about their function and biological relevance. Deep-sequencing data contain valuable information on a genome-wide scale. Yet, due to their highly repetitive nature, analysis of ERVs has been computationally challenging. We describe a bioinformatics tool called ERVmap to analyze transcription of unique sets of human ERVs in a range of cell types in health and disease settings. Our open-source code and accompanied web tool should facilitate researchers in all fields to study the expression patterns of ERVs in sequencing data and should lead to significant advancement in understanding the biological relevance of ERVs in health and disease.

Endogenous retroviruses (ERVs) are integrated retroviral elements that make up 8% of the human genome. However, the impact of ERVs on human health and disease is not well understood. While select ERVs have been implicated in diseases, including autoimmune disease and cancer, the lack of tools to analyze genome-wide, locus-specific expression of proviral autonomous ERVs has hampered the progress in the field. Here we describe a method called ERVmap, consisting of an annotated database of 3,220 human proviral ERVs and a pipeline that allows for locus-specific genome-wide identification of proviral ERVs that are transcribed based on RNA-sequencing data, and provide examples of the utility of this tool. Using ERVmap, we revealed cell-type–specific ERV expression patterns in commonly used cell lines as well as in primary cells. We identified 124 unique ERV loci that are significantly elevated in the peripheral blood mononuclear cells of patients with systemic lupus erythematosus that represent an IFN-independent signature. Finally, we identified additional tumor-associated ERVs that correlate with cytolytic activity represented by granzyme and perforin expression in breast cancer tissue samples. The open-source code of ERVmap and the accompanied web tool are made publicly available to quantify proviral ERVs in RNA-sequencing data with ease. Use of ERVmap across a range of diseases and experimental conditions has the potential to uncover novel disease-associated antigens and effectors involved in human health that is currently missed by focusing on protein-coding sequences.