Aberrant Expression of Pseudogene-Derived lncRNAs as an Alternative Mechanism of Cancer Gene Regulation in Lung Adenocarcinoma

Enhanced bovine genome annotation through integration of transcriptomics and epi-transcriptomics datasets facilitates genomic biology

BACKGROUND

The accurate identification of the functional elements in the bovine genome is a fundamental requirement for high-quality analysis of data informing both genome biology and genomic selection. Functional annotation of the bovine genome was performed to identify a more complete catalog of transcript isoforms across bovine tissues.

RESULTS

A total of 160,820 unique transcripts (50% protein coding) representing 34,882 unique genes (60% protein coding) were identified across tissues. Among them, 118,563 transcripts (73% of the total) were structurally validated by independent datasets (PacBio isoform sequencing data, Oxford Nanopore Technologies sequencing data, de novo assembled transcripts from RNA sequencing data) and comparison with Ensembl and NCBI gene sets. In addition, all transcripts were supported by extensive data from different technologies such as whole transcriptome termini site sequencing, RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression, chromatin immunoprecipitation sequencing, and assay for transposase-accessible chromatin using sequencing. A large proportion of identified transcripts (69%) were unannotated, of which 86% were produced by annotated genes and 14% by unannotated genes. A median of two 5' untranslated regions were expressed per gene. Around 50% of protein-coding genes in each tissue were bifunctional and transcribed both coding and noncoding isoforms. Furthermore, we identified 3,744 genes that functioned as noncoding genes in fetal tissues but as protein-coding genes in adult tissues. Our new bovine genome annotation extended more than 11,000 annotated gene borders compared to Ensembl or NCBI annotations. The resulting bovine transcriptome was integrated with publicly available quantitative trait loci data to study tissue-tissue interconnection involved in different traits and construct the first bovine trait similarity network.

CONCLUSIONS

These validated results show significant improvement over current bovine genome annotations.

Methionine-producing tumor micro(be) environment fuels growth of solid tumors

BACKGROUND

Recent studies have uncovered the near-ubiquitous presence of microbes in solid tumors of diverse origins. Previous literature has shown the impact of specific bacterial species on the progression of cancer. We propose that local microbial dysbiosis enables certain cancer phenotypes through provisioning of essential metabolites directly to tumor cells.

METHODS

16S rDNA sequencing of 75 patient lung samples revealed the lung tumor microbiome specifically enriched for bacteria capable of producing methionine. Wild-type (WT) and methionine auxotrophic (metA mutant) E. coli cells were used to condition cell culture media and the proliferation of lung adenocarcinoma (LUAD) cells were measured using SYTO60 staining. Further, colony forming assay, Annexin V Staining, BrdU, AlamarBlue, western blot, qPCR, LINE microarray and subcutaneous injection with methionine modulated feed were used to analyze cellular proliferation, cell-cycle, cell death, methylation potential, and xenograft formation under methionine restriction. Moreover, C14-labeled glucose was used to illustrate the interplay between tumor cells and bacteria.

RESULTS/DISCUSSION

Our results show bacteria found locally within the tumor microenvironment are enriched for methionine synthetic pathways, while having reduced S-adenosylmethionine metabolizing pathways. As methionine is one of nine essential amino acids that mammals are unable to synthesize de novo, we investigated a potentially novel function for the microbiome, supplying essential nutrients, such as methionine, to cancer cells. We demonstrate that LUAD cells can utilize methionine generated by bacteria to rescue phenotypes that would otherwise be inhibited due to nutrient restriction. In addition to this, with WT and metA mutant E. coli, we saw a selective advantage for bacteria with an intact methionine synthetic pathway to survive under the conditions induced by LUAD cells. These results would suggest that there is a potential bi-directional cross-talk between the local microbiome and adjacent tumor cells. In this study, we focused on methionine as one of the critical molecules, but we also hypothesize that additional bacterial metabolites may also be utilized by LUAD. Indeed, our radiolabeling data suggest that other biomolecules are shared between cancer cells and bacteria. Thus, modulating the local microbiome may have an indirect effect on tumor development, progression, and metastasis.

Downregulated Expression of CLEC9A as Novel Biomarkers for Lung Adenocarcinoma

Purpose

Abnormal CLEC9A expression is concerned with carcinogenesis. However, the role of CLEC9A in lung adenocarcinoma (LUAD) remains unknown. The goal of this study was to reveal the role of CLEC9A in LUAD based on bioinformatics and cellular functional experiments.

Materials and methods

Data available from The Cancer Genome Atlas (TCGA) were employed to study CLEC9A expression and mutations in LUAD. Expression and alterations of CLEC9A were analyzed using UALCAN and cBioPortal, respectively. Kaplan-Meier analysis was used to analyze the effect of CLEC9A on the survival of LUAD. Protein-protein interaction (PPI) network was built using GeneMANIA analysis. The similar genes of CLEC9A were obtained using GEPIA analysis, while co-expression genes correlated with CLEC9A were identified using LinkedOmics analysis. The effects of CLEC9A expression on immune cell infiltration was assessed. The effect of CLEC9A on the proliferation, apoptosis, cell cycle distribution, and invasion of human LUAD cells was detected in the LUAD cell line.

Results

CLEC9A was downregulated and the CLEC9A gene was often altered in LUAD. The survival of LUAD patients was correlated with the expression level of CLEC9A. The similar genes of CLEC9A were linked to functional networks involving positive regulation of interleukin-12 production, plasma membrane and CD40 receptor binding, primary immunodeficiency, intestinal immune network for IgA production, and cell adhesion molecules pathways. Cell cycle, apoptosis, EMT, and RAS/MAPK were significantly enriched pathways in positive and negative correlation genes with CLEC9A. A difference in the immune infiltration level of immune cell between the high and low CLEC9A expression groups was observed. Somatic cell copy number alternations (CNAs) of the CLEC9A, including arm-level gain and arm-level deletion, observably changed the infiltration levels of B cells, CD4+ T cells, macrophages, and neutrophils in LUAD. Except for LAG3, the expression of CD274, CTLA4, PDCD1, and TIGIT was positively correlated with the expression level of CLEC9A. After transfection, overexpression and knockdown of CLEC9A could affect the proliferation, apoptosis, cell cycle distribution, and invasion of LUAD cells.

Conclusion

CLEC9A is associated with prognosis and tumor immune microenvironment of LUAD, suggesting that CLEC9A may be considered as a novel biomarker for LUAD.

The transcriptional stress response and its implications in cancer treatment

Cancer cells require high levels of transcription to survive and maintain their cancerous phenotype. For several years, global transcription inhibitors have been used in the treatment of cancer. However, recent advances in understanding the functioning of the basal transcription machinery and the discovery of new drugs that affect the components of this machinery have generated a new boom in the use of this type of drugs to treat cancer. Inhibiting transcription at the global level in the cell generates a stress situation in which the cancer cell responds by overexpressing hundreds of genes in response to this transcriptional stress. Many of these over-transcribed genes encode factors that may be involved in the selection of cells resistant to the treatment and with a greater degree of malignancy. In this study, we reviewed various examples of substances that inhibit global transcription, as well as their targets, that have a high potential to be used against cancer. We also analysed what kinds of genes are overexpressed in the response to transcriptional stress by different substances and finally we discuss what types of studies are necessary to understand this type of stress response to have more tools to fight cancer.

Pseudogene legumain promotes thyroid carcinoma progression via the microRNA-495/autophagy pathway

The pseudogene legumain (LGMN) has been reported to regulate cancer cell biology. However, the role of LGMN in thyroid carcinoma remains unknown. Herein, Cell Counting Kit 8 and Transwell assays were performed to evaluate cellular proliferation and invasion capacity, respectively. In addition, a tube formation assay was performed to assess HUVEC angiogenesis. The results showed that LGMN depletion attenuated cellular proliferation, invasion and tube formation ability, and that LGMN expression was dysregulated in thyroid carcinoma tumors. Furthermore, patients with high LGMN expression levels exhibited a lower overall survival rate than those with low expression levels. LGMN and microRNA (miR)-495 modulated the expression levels of autophagy-related gene 3 (ATG3) and p62. Finally, ATG3 overexpression rescued the LGMN-regulated thyroid carcinoma phenotype. In conclusion, LGMN was found to promote thyroid carcinoma progression via the miR-495/autophagy axis, thus providing novel insights for understanding the pathogenesis of thyroid carcinoma.

Differential methylation of G-protein coupled receptor signaling genes in gastrointestinal neuroendocrine tumors

Neuroendocrine tumors (NETs) of the small intestine undergo large chromosomal and methylation changes. The objective of this study was to identify methylation differences in NETs and consider how the differentially methylated genes may impact patient survival. Genome-wide methylation and chromosomal copy number variation (CNV) of NETs from the small intestine and appendix were measured. Tumors were divided into three molecular subtypes according to CNV results: chromosome 18 loss (18LOH), Multiple CNV, and No CNV. Comparison of 18LOH tumors with MultiCNV and NoCNV tumors identified 901 differentially methylated genes. Genes from the G-protein coupled receptor (GPCR) pathways are statistically overrepresented in the differentially methylated genes. One of the highlighted genes from the GPCR pathway is somatostatin (SST), a clinical target for NETs. Patient survival based on low versus high methylation in all samples identified four significant genes (p < 0.05) OR2S2, SMILR, RNU6-653P, and AC010543.1. Within the 18LOH molecular subtype tumors, survival differences were identified in high versus low methylation of 24 genes. The most significant is TRHR (p < 0.01), a GPCR with multiple FDA-approved drugs. By separating NETs into different molecular subtypes based on chromosomal changes, we find that multiple GPCRs and their ligands appear to be regulated through methylation and correlated with survival. These results suggest opportunities for better treatment strategies for NETs based on molecular features.

Prospects Of Non-Coding Elements In Genomic Dna Based Gene Therapy

Gene therapy has made significant development since the commencement of the first clinical trials a few decades ago and has remained a dynamic area of research regardless of obstacles such as immune response and insertional mutagenesis. Progression in various technologies like next-generation sequencing (NGS) and nanotechnology has established the importance of non-coding segments of a genome, thereby taking gene therapy to the next level. In this review, we have summarized the importance of non-coding elements, highlighting the advantages of using full-length genomic DNA loci (gDNA) compared to complementary DNA (cDNA) or minigene, currently used in gene therapy. The focus of this review is to provide an overview of the advances and the future of potential use of gDNA loci in gene therapy, expanding the therapeutic repertoire in molecular medicine.

Non-Coding RNAs in Cancer Diagnosis and Therapy: Focus on Lung Cancer

Over the last several decades, clinical evaluation and treatment of lung cancers have largely improved with the classification of genetic drivers of the disease, such as EGFR, ALK, and ROS1. There are numerous regulatory factors that exert cellular control over key oncogenic pathways involved in lung cancers. In particular, non-coding RNAs (ncRNAs) have a diversity of regulatory roles in lung cancers such that they have been shown to be involved in inducing proliferation, suppressing apoptotic pathways, increasing metastatic potential of cancer cells, and acquiring drug resistance. The dysregulation of various ncRNAs in human cancers has prompted preclinical studies examining the therapeutic potential of restoring and/or inhibiting these ncRNAs. Furthermore, ncRNAs demonstrate tissue-specific expression in addition to high stability within biological fluids. This makes them excellent candidates as cancer biomarkers. This review aims to discuss the relevance of ncRNAs in cancer pathology, diagnosis, and therapy, with a focus on lung cancer.

MutCombinator: identification of mutated peptides allowing combinatorial mutations using nucleotide-based graph search

Motivation

Proteogenomics has proven its utility by integrating genomics and proteomics. Typical approaches use data from next-generation sequencing to infer proteins expressed. A sample-specific protein sequence database is often adopted to identify novel peptides from matched mass spectrometry-based proteomics; nevertheless, there is no software that can practically identify all possible forms of mutated peptides suggested by various genomic information sources.

Results

We propose MutCombinator, which enables us to practically identify mutated peptides from tandem mass spectra allowing combinatorial mutations during the database search. It uses an upgraded version of a variant graph, keeping track of frame information. The variant graph is indexed by nine nucleotides for fast access. Using MutCombinator, we could identify more mutated peptides than previous methods, because combinations of point mutations are considered and also because it can be practically applied together with a large mutation database such as COSMIC. Furthermore, MutCombinator supports in-frame search for coding regions and three-frame search for non-coding regions.

Availability and implementation

https://prix.hanyang.ac.kr/download/mutcombinator.jsp.

Supplementary information

Supplementary data are available at Bioinformatics online.

Deregulation of a Cis-Acting lncRNA in Non-small Cell Lung Cancer May Control HMGA1 Expression

Background

Long non-coding RNAs (lncRNAs) have long been implicated in cancer-associated phenotypes. Recently, a class of lncRNAs, known as cis-acting, have been shown to regulate the expression of neighboring protein-coding genes and may represent undiscovered therapeutic action points. The chromatin architecture modification gene HMGA1 has recently been described to be aberrantly expressed in lung adenocarcinoma (LUAD). However, the mechanisms mediating the expression of HMGA1 in LUAD remain unknown. Here we investigate the deregulation of a putative cis-acting lncRNA in LUAD, and its effect on the oncogene HMGA1.

Methods

LncRNA expression was determined from RNA-sequencing data of tumor and matched non-malignant tissues from 36 LUAD patients. Transcripts with significantly deregulated expression were identified and validated in a secondary LUAD RNA-seq dataset (TCGA). SiRNA-mediated knockdown of a candidate cis-acting lncRNA was performed in BEAS-2B cells. Quantitative real-time PCR was used to observe the effects of lncRNA knockdown on the expression of HMGA1.

Results

We identified the lncRNA RP11.513I15.6, which we refer to as HMGA1-lnc, neighboring HMGA1 to be significantly downregulated in both LUAD cohorts. Conversely, we found HMGA1 significantly overexpressed in LUAD and anticorrelated with HMGA1-lnc. In vitro experiments demonstrated siRNA-mediated inhibition of HMGA1-lnc in immortalized non-malignant lung epithelial cells resulted in a significant increase in HMGA1 gene expression.

Conclusion

Our results suggest that HMGA1-lnc is a novel cis-acting lncRNA that negatively regulates HMGA1 gene expression in lung cells. Further characterization of this regulatory mechanism may advance our understanding of the maintenance of lung cancer phenotypes and uncover a novel therapeutic intervention point for tumors driven by HMGA1.

Pseudogene KRT17P3 drives cisplatin resistance of human NSCLC cells by modulating miR-497-5p/mTOR

Chemoresistance is a major obstacle in non–small cell lung cancer (NSCLC) treatment. The pseudogene keratin 17 pseudogene 3 (KRT17P3) has been previously shown to be upregulated in lung cancer tissues of patients with cisplatin resistance. In the present study, RT‐qPCR was performed to evaluate KRT17P3 levels in plasma samples collected from 30 cisplatin‐resistant and 32 cisplatin‐sensitive patients. We found that the plasma level of KRT17P3 is upregulated in cisplatin‐resistant patients, and the increased expression of plasma KRT17P3 is associated with poor chemotherapy response. Functional studies demonstrated that KRT17P3 overexpression in cultured NSCLC cells increases cell viability and decreases apoptosis upon cisplatin treatment in vitro and in vivo, while KRT17P3 knockdown has the opposite effect. Mechanistically, bioinformatics analysis, RNA immunoprecipitation, and dual luciferase reporter assay indicated that KRT17P3 acts as a molecular sponge for miR‐497‐5p and relieves the binding of miR‐497‐5p to its target gene mTOR. Rescue experiments validated the functional interaction between KRT17P3, miR‐497‐5p, and mTOR. Taken together, our findings indicate that KRT17P3/miR‐497‐5p/mTOR regulates the chemosensitivity of NSCLC, suggesting a potential therapeutic target for cisplatin‐resistant NSCLC patients. KRT17P3 may be a potential peripheral blood marker of NSCLC patients resistant to cisplatin.

Assessment of long non-coding RNA expression reveals novel mediators of the lung tumour immune response

The tumour immune microenvironment is a crucial mediator of lung tumourigenesis, and characterizing the immune landscape of patient tumours may guide immunotherapy treatment regimens and uncover novel intervention points. We sought to identify the landscape of tumour-infiltrating immune cells in the context of long non-coding RNA (lncRNAs), known regulators of gene expression. We examined the lncRNA profiles of lung adenocarcinoma (LUAD) tumours by interrogating RNA sequencing data from microdissected and non-microdissected samples (BCCRC and TCGA). Subsequently, analysis of single-cell RNA sequencing data from lung tumours and flow-sorted healthy peripheral blood mononuclear cells identified lncRNAs in immune cells, highlighting their biological and prognostic relevance. We discovered lncRNA expression patterns indicative of regulatory relationships with immune-related protein-coding genes, including the relationship between AC008750.1 and NKG7 in NK cells. Activation of NK cells in vitro was sufficient to induce AC008750.1 expression. Finally, siRNA-mediated knockdown of AC008750.1 significantly impaired both the expression of NKG7 and the anti-tumour capacity of NK cells. We present an atlas of cancer-cell extrinsic immune cell-expressed lncRNAs, in vitro evidence for a functional role of lncRNAs in anti-tumour immune activity, which upon further exploration may reveal novel clinical utility as markers of immune infiltration.

Long noncoding RNA ARHGAP27P1 inhibits gastric cancer cell proliferation and cell cycle progression through epigenetically regulating p15 and p16

Long noncoding RNAs (lncRNAs) have emerged as important regulators in the development and progression of gastric cancer (GC). ARHGAP27P1 is a pseudogene-derived lncRNA, and it has been found to be associated with GC in our preliminary study, but this association has not been studied further. Herein, we confirmed that ARHGAP27P1 was significantly downregulated in GC tissues, plasma and cells. Low expression of ARHGAP27P1 was closely associated with advanced TNM stage, increased invasion depth and lymphatic metastasis. Low ARHGAP27P1 expression also predicted a poor prognosis in GC patients. Functionally, overexpression of ARHGAP27P1 inhibited proliferation, invasion, and migration in GC cells, while silencing of ARHGAP27P1 showed the opposite effects. Mechanistic investigations showed that ARHGAP27P1 had a key role in G0/G1 arrest. We further demonstrated that ARHGAP27P1 was associated with Jumonji-domain containing 3 (JMJD3) and that this association was required for the demethylation of H3K27me3, thereby epigenetically activating expression of p15, p16 and p57. Moreover, knockdown of JMJD3, p15, or p16 consistently reversed the inhibitory effects of ARHGAP27P1 in cell proliferation and cell cycle progression. Taken together, these results suggest that lncRNA ARHGAP27P1, as a novel cell cycle regulator, may serve as a potential target for GC prevention and treatment in human GC.