Genome-wide identification of lncRNAs as novel prognosis biomarkers of glioma

Exosomal long non-coding RNAs in glioblastoma

Glioblastoma multiforme (GBM) is the most prevalent primary tumor found in the central nervous system, accounting for 70% of all adult brain tumors. The median overall survival rate is one year post-diagnosis with treatment, and only four months without treatment. Current GBM diagnostic methods, such as magnetic resonance imaging (MRI), surgery, and brain biopsies, have limitations. These include difficulty distinguishing between tumor recurrence and post-surgical necrotic regions, and operative risks associated with obtaining histological samples through direct surgery or biopsies. Consequently, there is a need for rapid, inexpensive, and minimally invasive techniques for early diagnosis and improved subsequent treatment. Research has shown that tumor-derived exosomes containing various long non-coding RNAs (lncRNAs) play critical regulatory roles in immunomodulation, cancer metastasis, cancer development, and drug resistance in GBM. They regulate genes that enhance cancer growth and progression and alter the expression of several key signaling pathways. Due to the specificity and sensitivity of exosomal lncRNAs, they have the potential to be used as biomarkers for early diagnosis and prognosis, as well as to monitor a patient's response to chemotherapy for GBM. In this review, we discuss the role of exosomal lncRNAs in the pathogenesis of GBM and their potential clinical applications for early diagnosis.

LncRNA as potential biomarker and therapeutic target in glioma

Glioma is the most frequent type of malignant tumor in the central nervous system, accounting for about 80% of primary malignant brain tumors, usually with a poor prognosis. A number of studies have been conducted on the molecular abnormalities in glioma to further understand its pathogenesis, and it has been found that lncRNAs (long non-coding RNA) play a key role in angiogenesis, tumor growth, infiltration and metastasis of glioma. Since specific lncRNAs have an aberrant expression in brain tissue, cerebrospinal fluid as well as peripheral circulation of glioma patients, they are considered to be potential biomarkers. This review focuses on the biological characteristics of lncRNA and its value as a biomarker for glioma diagnosis and prognosis. Moreover, in view of the role of lncRNAs in glioma proliferation and chemoradiotherapy resistance, we discussed the feasibility for lncRNAs as therapeutic targets. Finally, the persisting deficiencies and future prospects of using lncRNAs as clinical biomarkers and therapeutic targets were concluded.

Latest updates on cellular and molecular biomarkers of gliomas

Gliomas are the most common central nervous system malignancies, compromising almost 80% of all brain tumors and is associated with significant mortality. The classification of gliomas has shifted from basic histological perspective to one that is based on molecular biomarkers. Treatment of this type of tumors consists currently of surgery, chemotherapy and radiation therapy. During the past years, there was a limited development of effective glioma diagnostics and therapeutics due to multiple factors including the presence of blood-brain barrier and the heterogeneity of this type of tumors. Currently, it is necessary to highlight the advantage of molecular diagnosis of gliomas to develop patient targeted therapies based on multiple oncogenic pathway. In this review, we will evaluate the development of cellular and molecular biomarkers for the diagnosis of gliomas and the impact of these diagnostic tools for better tailored and targeted therapies.

Mechanisms of Long Non-Coding RNAs in Biological Characteristics and Aerobic Glycolysis of Glioma

Glioma is the most common and aggressive tumor of the central nervous system. The uncontrolled proliferation, cellular heterogeneity, and diffusive capacity of glioma cells contribute to a very poor prognosis of patients with high grade glioma. Compared to normal cells, cancer cells exhibit a higher rate of glucose uptake, which is accompanied with the metabolic switch from oxidative phosphorylation to aerobic glycolysis. The metabolic reprogramming of cancer cell supports excessive cell proliferation, which are frequently mediated by the activation of oncogenes or the perturbations of tumor suppressor genes. Recently, a growing body of evidence has started to reveal that long noncoding RNAs (lncRNAs) are implicated in a wide spectrum of biological processes in glioma, including malignant phenotypes and aerobic glycolysis. However, the mechanisms of diverse lncRNAs in the initiation and progression of gliomas remain to be fully unveiled. In this review, we summarized the diverse roles of lncRNAs in shaping the biological features and aerobic glycolysis of glioma. The thorough understanding of lncRNAs in glioma biology provides opportunities for developing diagnostic biomarkers and novel therapeutic strategies targeting gliomas.

Phenotypic plasticity of myeloid cells in glioblastoma development, progression, and therapeutics

Glioblastoma (GBM) is the most common and malignant type of intracranial tumors with poor prognosis. Accumulating evidence suggests that phenotypic alterations of infiltrating myeloid cells in the tumor microenvironment are important for GBM progression. Conventional tumor immunotherapy commonly targets T-cells, while innate immunity as a therapeutic target is an emerging field. Targeting infiltrating myeloid cells that induce immune suppression in the TME provides a novel direction to improve the prognosis of patients with GBM. The factors released by tumor cells recruit myeloid cells into tumor bed and reprogram infiltrating myeloid cells into immunostimulatory/immunosuppressive phenotypes. Reciprocally, infiltrating myeloid cells, especially microglia/macrophages, regulate GBM progression and affect therapeutic efficacy. Herein, we revisited biological characteristics and functions of infiltrating myeloid cells and discussed the recent advances in immunotherapies targeting infiltrating myeloid cells in GBM. With an evolving understanding of the complex interactions between infiltrating myeloid cells and tumor cells in the tumor microenvironment, we will expand novel immunotherapeutic regimens targeting infiltrating myeloid cells in GBM treatment and improve the outcomes of GBM patients.

Identification and validation of a novel eight mutant-derived long non-coding RNAs signature as a prognostic biomarker for genome instability in low-grade glioma

Long non-coding RNAs (lncRNAs) comprise an integral part of the eukaryotic transcriptome. Alongside proteins, lncRNAs modulate lncRNA-based gene signatures of unstable transcripts, play a crucial role as antisense lncRNAs to control intracellular homeostasis and are implicated in tumorigenesis. However, the role of genomic instability-associated lncRNAs in low-grade gliomas (LGG) has not been fully explored. In this study, lncRNAs expression and somatic mutation profiles in low-grade glioma genome were used to identify eight novel mutant-derived genomic instability-associated lncRNAs including H19, FLG-AS1, AC091932.1, AC064875.1, AL138767.3, AC010273.2, AC131097.4 and ISX-AS1. Patients from the LGG gene mutagenome atlas were grouped into training and validation sets to test the performance of the signature. The genomic instability-associated lncRNAs signature (GILncSig) was then validated using multiple external cohorts. A total of 59 novel genomic instability-associated lncRNAs in LGG were used for least absolute shrinkage and selection operator (Lasso), single and multifactor Cox regression analysis using the training set. Furthermore, the independent predictive role of risk features in the training and validation sets were evaluated through survival analysis, receiver operating feature analysis and construction of a nomogram. Patients with IDH1 mutation status were grouped into two different risk groups based on the GILncSig score. The low-risk group showed a relatively higher rate of IDH1 mutations compared with patients in the high-risk group. Furthermore, patients in the low-risk group had better prognosis compared with patients in the high-risk group. In summary, this study reports a reliable prognostic prediction signature and provides a basis for further investigation of the role of lncRNAs on genomic instability. In addition, lncRNAs in the signature can be used as new targets for treatment of LGG.

LINC00511 knockdown suppresses glioma cell malignant progression through miR-15a-5p/AEBP1 axis

BACKGROUND

A strong relationship between long intergenic non-protein coding RNA 511 (LINC00511) and glioma has been previously reported but the mechanism of LINC00511 in glioma is yet to be determined. This study examined the mechanism of LINC00511 in glioma.

METHODS

The expression of LINC00511 in glioma was determined by bioinformatics analysis and real-time quantitative PCR (RT-qPCR) analysis. The target relationship between genes was predicted by starBase, TargetScan, and was verified by dual-luciferase. Subsequently, siRNA targeting LINC00511 (siLINC00511) and miR-15a-5p mimic were transfected into glioma cells to examine the effect on biological characteristics using cell counting kit-8, clone formation, flow cytometry, wound-healing, and transwell. MiR-15a-5p inhibitor and AEBP1 were used for in vitro rescue experiments, and tumorigenesis assay and immunohistochemical assays were performed for in vivo experiments. Epithelial-mesenchymal transition (EMT) and p65 phosphorylation were examined by Western blot.

RESULTS

LINC00511 was predicted and verified to be up-regulated in glioma. SiLINC00511 suppressed cell viability, proliferation, migration and invasion, accelerated apoptosis of glioma cells. Mechanically, siLINC00511 promoted E-cadherin expression but suppressed N-cadherin and Snail expressions. MiR-15a-5p bound to LINC00511, and miR-15a-5p inhibitor partially reversed the effect and regulation of siLINC00511 on glioma cells. AEBP1, a target gene of miR-15a-5p, could activate p65 phosphorylation to promote EMT protein expression and partially reverse the inhibitory effect of miR-15a-5p mimic on the malignant phenotype of glioma cells. SiLINC00511 inhibited tumor growth, down-regulated miR-15a-5p expression and up-regulated AEBP1 and Ki67 expressions in vivo.

CONCLUSION

LINC00511 knockdown inhibits glioma cell progression via miR-15a-5p/AEBP1 axis.

LINC00908 Promotes Diffuse Large B-Cell Lymphoma Development by Down-Regulating miR-671-5p

Introduction

Emerging evidence has revealed that long noncoding RNA (lncRNA) play important role in almost all kinds of human cancers. LINC00908 has been reported to be involved in the development of prostate cancer, colorectal cancer and gastric cancer which was functioned as an oncogene. However, the potential biology role and molecular mechanism of LINC00908 in diffuse large B-cell lymphoma are still unclear.

Methods

LINC00908 and miR-671-5p expression were evaluated in DLBCL tissues and cell lines using RT-qPCR. CCK-8 and transwell assay were used to analyze the in vitro role of LINC00908 in DLBCL progression. The xenograft model was used to explore the in vivo role of LINC00908 in DLBCL growth. The physical interaction between LINC00908 and miR-671-5p was confirmed using bioinformatics analysis and a dual luciferase assay, RIP and RNA pull down.

Results

The expression of LINC00908 was markedly up-regulated in diffuse large B-cell lymphoma tissues and cell lines, and the decreased expression of LINC00908 significantly inhibited diffuse large B-cell lymphoma cell proliferation and invasion. Then, we revealed that LINC00908 directly interacted with miR-671-5p, which was down-regulated in diffuse large B-cell lymphoma cells and highly expressed with LINC00908 knockdown. Moreover, luciferase reporter assays and RNA immunoprecipitation (RIP) assay further proved that miR-671-5p is a direct target of LINC00908 in diffuse large B-cell lymphoma cells. Rescue experiments were also performed, and we confirmed that LINC00908 acts as an oncogene role in diffuse large B-cell lymphoma through miR-671-5p. Finally, the influence of LINC00908 silence significantly inhibited diffuse large B-cell lymphoma growth in vivo.

Conclusion

LINC00908 promotes malignancy of diffuse large B-cell lymphoma through regulating miR-671-5p.

The role of M6A modification in the regulation of tumor-related lncRNAs

N⁶-methyladenosine (m⁶A) is the most abundant modification in eukaryotic cells, and it regulates RNA transcription, processing, splicing, degradation, and translation. Long non-coding RNAs (lncRNAs), as transcriptional products with no or limited protein coding ability more than 200 nt in length, play an important role in epigenetic modification, mRNA transcription, splicing, stability, translation, and other biological functions. Extensive studies have shown that both m⁶A modification and lncRNAs are involved in the pathogenesis of various diseases, such as kinds of cancers, heart failure, Alzheimer’s disease, periodontitis, human abdominal aortic aneurysm, and obesity. To date, m⁶A modification has been identified as an important biological function in enrichment and regulation of lncRNAs. In this review, we summarize the role of m⁶A modification in the regulation and function of tumor-related lncRNAs. Moreover, we discuss the potential applications and possible future directions in the field.

Silencing of Long Non-coding RNA ENST00000606790.1 Inhibits the Malignant Behaviors of Papillary Thyroid Carcinoma through the PI3K/AKT Pathway

PURPOSE

This study aimed to investigate the role and mechanism of lncRNA ENST00000606790.1 (ENST) in promoting the progression of papillary thyroid carcinoma (PTC).

METHODS

The expression of ENST in human PTC and normal para-cancerous thyroid (NPTC) tissues or cell lines was determined by RT-qPCR. Cell growth was determined by CCK8 assay. Cell colony formation was determined by cell colony formation assay. Cell cycle analysis was performed by staining cells with PI (Propidium Iodide). Cell invasion was assessed by transwell assay. Protein expression was examined by western-blot. siRNA was constructed to inhibit the expression of ENST. 740-Y-P was used to activate PI3K. The correlation between ENST expression and clinical outcomes was analyzed.

RESULTS

ENST was significantly up-regulated in PTC tissues or PTC cell lines (PTC and IHH4 cell lines), compared to NPTC tissues or normal cell lines, respectively. High expression of ENST was strongly correlated to lymph node metastasis and tumor size at diagnosis. Silencing of ENST significantly inhibited cell growth and colony formation, arrested the cell cycle at G2/M phase, upregulated the expression of CHK1, downregulated the expression of CDC25C, and inhibited cell invasion. Silencing of ENST significantly down-regulated the expression of PI3K, p-PI3K, AKT, and p-AKT in IHH4 cells. Furthermore, treatment with the PI3K activator  740-Y-P partially abolished the effect of silencing of ENST on PTC cells.

CONCLUSIONS

Overall, our results demonstrated that ENST can promote PTC progression by activating the PI3K/AKT signaling pathway, suggesting that ENST can serve as a potential biomarker and new therapeutic target for patients with PTC.

Identification and Verification on Prognostic Index of Lower-Grade Glioma Immune-Related LncRNAs

Previous studies have shown that the prognosis of patients with lower-grade glioma (LGG) is closely related to the infiltration of immune cells and the expression of long non-coding RNAs (lncRNAs). In this paper, we applied single-sample gene set enrichment analysis (ssGSEA) algorithm to evaluate the expression level of immune genes from tumor tissues in The Cancer Genome Atlas (TCGA) database, and divided patients into the high immune group and the low immune group, which were separately analyzed for differential expression. Venn analysis was taken to select 36 immune-related lncRNAs. To construct a prognostic model of LGG based on immune-related lncRNAs, we divided patients into a training set and a verification set at a ratio of 2:1. Univariate Cox regression and the Least Absolute Shrinkage and Selection Operator (LASSO) regression were performed to select 11 immune-related lncRNAs associated with the prognosis of LGG, and based on these selected lncRNAs, the risk scoring model was constructed. Through Kaplan-Meier analysis, the overall survival (OS) of patients in the high-risk group was significantly lower than that of the low-risk group. Then, established a nomogram including age, gender, neoplasm histologic grade, and risk score. Meanwhile, the predictive performance of the model was evaluated by calculating the C-index, drawing the calibration chart, the clinical decision curve as well as the Receiver Operating Characteristic (ROC) curve. Similar results were obtained by utilizing the validation data to verify the above consequences. Based on the TIMER database, the correlation analysis showed that the 11 immune-related lncRNAs risk score of LGG were in connection with the infiltration of the subtypes of immune cells. Subsequently, we performed enrichment analysis, whose results showed that these immune-related lncRNAs played important roles in the progress of LGG. In conclusion, these 11 immune-related lncRNAs have the potential to predict the prognosis of patients with LGG, which may play a key role in the development of LGG.

Individualized Prediction of Survival by a 10-Long Non-coding RNA-Based Prognostic Model for Patients With Breast Cancer

Deregulations of long non-coding RNAs (lncRNAs) have been implicated in the progression of breast cancer (BC). However, the prognostic values of those lncRNAs in BC remain elusive. This study aimed at constructing a lncRNA-based prognostic model to improve the clinical management of BC. Systematic investigation of lncRNA expression profiles and clinical data from The Cancer Genome Atlas (TCGA) database were utilized to establish a 10-lncRNA signature. The prognostic signature efficiently discriminated patients with significantly different prognosis regardless of intrinsic molecular subtypes and tumor-node-metastasis (TNM) stage. A combined model was constructed by multivariate Cox proportional hazards regression (CPHR) analysis, which combined the lncRNA-based signature with certain clinical risk factors (TNM stage, age, and human epidermal growth factor receptor 2 status). This model predicted a survival probability that closely corresponds to the actual survival probability. With respect to the entire set, the time-dependent receiver-operating characteristic curves revealed that the area under the curve of this model was the highest than any of the clinical risk factors. Moreover, functional enrichment analysis indicated that the molecular signature was mainly involved in DNA replication, which was firmly related to BC tumorigenesis. Consistent with the discovery, the knockdown of LHX1-DT, one of the 10 prognostic lncRNAs, attenuated the proliferation of BC cells in vitro and in vivo. Taken together, our study constructed a novel 10-lncRNA signature for prediction prognosis, and the signature-based model could provide new insight into accurate management of BC patients.

Identification of the potential biomarkers in patients with glioma: a weighted gene co-expression network analysis

Glioma is the most common brain tumor with high mortality. However, there are still challenges for the timely and accurate diagnosis and effective treatment of the tumor. One hundred and twenty-one samples with grades II, III and IV from the Gene Expression Omnibus database were used to construct gene co-expression networks to identify hub modules closely related to glioma grade, and performed pathway enrichment analysis on genes from significant modules. In gene co-expression network constructed by 2345 differentially expressed genes from 121 gene expression profiles for glioma, we identified the black and blue modules that associated with grading. The module preservation analysis based on 118 samples indicates that the two modules were replicable. Enrichment analysis showed that the extracellular matrix genes were enriched for blue module, while cell division genes were enriched for black module. According to survival analysis, 21 hub genes were significantly up-regulated and one gene was significantly down-regulated. What’s more, IKBIP, SEC24D, and FAM46A are the genes with little attention among the 22 hub genes. In this study, IKBIP, SEC24D, and FAM46A related to glioma were mentioned for the first time to the current knowledge, which might provide a new idea for us to study the disease in the future. IKBIP, SEC24D and FAM46A among the 22 hub genes identified that are related to the malignancy degree of glioma might be used as new biomarkers to improve the diagnosis, treatment and prognosis of glioma.

HOTAIRM1, an enhancer lncRNA, promotes glioma proliferation by regulating long-range chromatin interactions within HOXA cluster genes

The long noncoding RNA HOTAIRM1 reportedly plays important roles in acute myeloid leukemia, gastric cancer and colorectal cancer. Here, we analyzed potential function of HOTAIRM1 in glioma and asked whether it participates in long-range chromatin interactions. We monitored expression of HOTAIRM1 in glioma tissues and correlated levels with patient survival using the TCGA dataset. HOTAIRM1 was highly expressed in glioma tissue, with high levels associated with shortened patient survival time. We then suppressed HOTAIRM1 activity in the human glioblastoma U251 line using CRISPR-cas9 to knock in a truncating polyA fragment. Reporter analysis of these and control cells confirmed that the HOTAIRM1 locus serves as an active enhancer. We then performed Capture-C analysis to identify target genes of that locus and applied RNA antisense purification to assess chromatin interactions between the HOTAIRM1 locus and HOXA cluster genes. HOTAIRM1 knockdown in glioma cells decreased proliferation and reduced expression of HOXA cluster genes. HOTAIRM1 regulates long-range interactions between the HOTAIRM1 locus and HOXA genes. Our work suggests a new mechanism by which HOTAIRM1 regulates glioma progression by regulating high-order chromatin structure and could suggest novel therapeutic targets to treat an intractable cancer.

Identification of Novel lncRNA Markers in Glioblastoma Multiforme and Their Clinical Significance: A Study Based on Multiple Sequencing Data

Background

Long non-coding RNAs (lncRNAs) have been verified to have a vital role in the progression of glioblastoma multiforme (GBM). Our research was about to identify the potential lncRNAs which was closely associated with the pathogenesis and prognosis of glioblastoma multiforme.

Methods

All RNA sequence profiling data from patients with GBM were obtained from The Genotype-Tissue Expression (GTEx) and The Cancer Genome Atlas (TCGA). Differently expressed genes identified from GBM and control samples were used to construct competing endogenous RNA (ceRNA) network and perform corresponding functional enrichment analysis. Univariate Cox regression followed by lasso regression and multivariate Cox was used to validate independent lncRNA factors and construct a risk prediction model. Quantitative polymerase chain reaction (qPCR) was performed to verify the expression levels of potential lncRNA biomarkers in human GBM clinical specimens. A gene set enrichment analysis (GSEA) was subsequently conducted to explore potential signaling pathways in which critical lncRNAs may be involved. Moreover, nomogram plot was applied based on our prediction model and significant clinical covariates to visualize the prognosis of GBM patients.

Results

A total of 2023 differentially expressed genes (DEGs) including 56 lncRNAs, 1587 message RNAs (mRNAs) and 380 other RNAs were included. Based on predictive databases, 16lncRNAs, 32 microRNAs (miRNAs) and 99 mRNAs were used to construct a ceRNA network. Moreover, we performed a novel risk prediction model with 5 potential prognostic lncRNAs, in which 4 of them were newly identified in GBM, to predict the prognosis of GBM patients. Finally, a nomogram plot was constructed to illustrate the potential relationship between the prognosis of GBM and our risk prediction model and significant clinical covariates.

Conclusion

In this study, we identified 4 novel potential lncRNA biomarkers and constructed a prediction model of GBM prognosis. A simple-to-use nomogram was provided for further clinical application.

SNHG29 regulates miR-223-3p/CTNND1 axis to promote glioblastoma progression via Wnt/β-catenin signaling pathway

Background

Glioblastoma has been seen as the most common malignancy of brain tumor. Emerging reports has claimed that SNHG29 (LRRC75A-AS1) was involved in several biological processes via modulation of signaling pathway, and served as an malignant facilitatorin osteosarcoma. However, the specific role of SNHG29 in glioblastoma remains unknown.

Methods

RT-qPCR and microarray were operated to measure genes expression. Western blot was performed to examine protein expression. CCK-8 and colony formation assays were used to evaluate cell proliferation. Cell migration was tested by transwell assay. Nuclear-cytoplasmic fractionation was conducted to locate SNHG29. The binding capacity of miR-223-3p to SNHG29 or CTNND1 3′UTR was verified by RIP and luciferase reporter assay.

Results

SNHG29 presented high expression in glioblastoma to boost cell proliferation, migration and EMT process. In addition, miR-223-3p was validated to bind with SNHG29 after prediction and screening. Furthermore, miR-223-3p was proved to be a negative regulator for its target CTNND1. Then, the inhibition on cell proliferation, migration and EMT process resulted from SNHG29 knockdown was recovered by CTNND1 overexpression. At last, the inhibitive impacts on cell proliferation, migration and EMT process of CTNND1 deficiency was abrogated by LiCl.

Conclusions

In conclusion, SNHG29 regulates miR-223-3p/CTNND1 axis to promote glioblastoma progression via Wnt/β-catenin signaling pathway, offering a potential therapeutic point for glioblastoma patients.

The Long Noncoding RNA LINC00908 Facilitates Hepatocellular Carcinoma Progression Via Interaction With Sox-4

Background

The hepatocellular carcinoma (HCC) is a highly aggressive and common malignancy worldwide. Accumulating evidence has demonstrated a pivotal role of long noncoding RNAs (lncRNAs) in various tumors. However, the function of intergenic lncRNA LINC00908 is still unknown in HCC.

Methods

The RT-qPCR method was used to quantify the expression of LINC00908. Migration and viability assay were performed to evaluate the in vitro effect and xenograft tumor model was used to measure the in vivo effect. Immunoblot was used to identify the association of LINC00908 with Sox-4 and the stability of Sox-4.

Results

We found a novel lncRNA related to HCC. LINC00908 is highly expressed in tumorous tissues and cell lines compared with normal ones. High LINC00908 expression correlated with advanced TNM stages, tumor size and metastasis. LINC00908 promoted the migration and viability of HCC cells. The in vivo xenograft tumor growth and proliferation were also enhanced by LINC00908 overexpression and inhibited by LINC00908 silence. LINC00908 physically interacted with Sox-4, and the association between LINC00908 and Sox-4 increased the stability of Sox-4 by reducing proteasomal degradation.

Conclusion

Taken together, our current work has identified a novel lncRNA LINC00908 which is critically involved in HCC progression. The LINC00908-Sox-4 axis might provide a new and potential target for pharmaceutical therapies.

Microsyntenic Clusters Reveal Conservation of lncRNAs in Chordates Despite Absence of Sequence Conservation

Homologous long non-coding RNAs (lncRNAs) are elusive to identify by sequence similarity due to their fast-evolutionary rate. Here we develop LincOFinder, a pipeline that finds conserved intergenic lncRNAs (lincRNAs) between distant related species by means of microsynteny analyses. Using this tool, we have identified 16 bona fide homologous lincRNAs between the amphioxus and human genomes. We characterized and compared in amphioxus and Xenopus the expression domain of one of them, Hotairm1, located in the anterior part of the Hox cluster. In addition, we analyzed the function of this lincRNA in Xenopus, showing that its disruption produces a severe headless phenotype, most probably by interfering with the regulation of the Hox cluster. Our results strongly suggest that this lincRNA has probably been regulating the Hox cluster since the early origin of chordates. Our work pioneers the use of syntenic searches to identify non-coding genes over long evolutionary distances and helps to further understand lncRNA evolution.