Non-coding RNA Genes Research Articles

Re-evaluating the impact of alternative RNA splicing on proteomic diversity.

Alternative splicing (AS) constitutes a mechanism by which protein-coding genes and long non-coding RNA (lncRNA) genes produce more than a single mature transcript. From plants to humans, AS is a powerful process that increases transcriptome complexity. Importantly, splice variants produced from AS can potentially encode for distinct protein isoforms which can lose or gain specific domains and, hence, differ in their functional properties. Advances in proteomics have shown that the proteome is indeed diverse due to the presence of numerous protein isoforms. For the past decades, with the help of advanced high-throughput technologies, numerous alternatively spliced transcripts have been identified. However, the low detection rate of protein isoforms in proteomic studies raised debatable questions on whether AS contributes to proteomic diversity and on how many AS events are really functional. We propose here to assess and discuss the impact of AS on proteomic complexity in the light of the technological progress, updated genome annotation, and current scientific knowledge.

Combining multi-omics approaches to prioritize the variant-regulated functional long non-coding RNAs in autism spectrum disorder.

Rising evidence has indicated that long non-coding RNA (lncRNA) may play an essential role in the development of autism spectrum disorder (ASD). However, identifying the lncRNAs associated with ASD and the risk loci on them remains a major challenge. This study aims to identify potential causative variants and explore the related mechanisms. By leveraging differential expression analysis, WGCNA analysis and cis-expression quantitative analysis, our study mined functional SNPs with the regulated long non-coding RNA genes in brain tissues. We recruited 611 ASD children and 645 healthy children in the case-control study. Total 68 different expressed lncRNAs were validated by calculating the brain tissue-specific expression using RNA-seq data. By the WGCNA method, 9 functional lncRNAs classified as e-lncRNA were found to interact with 976 ASD risk genes. Furthermore, we mined functional SNPs regulated long non-coding RNAs in brain tissues. We analyzed the association between candidate SNPs and ASD risks in Chinese children, which showed BDNF-AS rs1565228 allele G to C reduced the risk of ASD (OR = 0.81, 95%CI: 0.66-0.98). Further bioinformatics analysis showed that the variant rs1565228 C>G with the low binding affinity of transcription factor SRF caused the decreased expression of lncRNA BDNF-AS. Our study revealed that rs2295412 in the non-coding sequence of the lncRNA gene region was significantly associated with the risk of ASD. These findings suggested that the SNPs in the non-coding region of lncRNA may play important roles in the genetic susceptibility of ASD, which may facilitate the early screening of ASD.

Full-Length RNA Sequencing Provides Insights into Goldfish Evolution under Artificial Selection.

Goldfish Carassius auratus is an ideal model for exploring fish morphology evolution. Although genes underlying several ornamental traits have been identified, little is known about the effects of artificial selection on embryo gene expression. In the present study, hybrid transcriptome sequencing was conducted to reveal gene expression profiles of Celestial-Eye (CE) and Ryukin (RK) goldfish embryos. Full-length transcriptome sequencing on the PacBio platform identified 54,218 and 54,106 transcript isoforms in CE and RK goldfish, respectively. Of particular note was that thousands of alternative splicing (AS) and alternative polyadenylation (APA) events were identified in both goldfish breeds, and most of them were inter-breed specific. RT-PCR and Sanger sequencing showed that most of the predicted AS and APA were correct. Moreover, abundant long non-coding RNA and fusion genes were detected, and again most of them were inter-breed specific. Through RNA-seq, we detected thousands of differentially expressed genes (DEGs) in each embryonic stage between the two goldfish breeds. KEGG enrichment analysis on DEGs showed extensive differences between CE and RK goldfish in gene expression. Taken together, our results demonstrated that artificial selection has led to far-reaching influences on goldfish gene expression, which probably laid the genetic basis for hundreds of goldfish variations.

Influence of the Cultivation Conditions of the Glioblastoma Neurosphere on the Expression of MALAT1 and LINCROR Long Non-coding RNA Genes

Glioblastoma multiforme (GBM) is the most aggressive malignant brain tumor. One of the reasons for the resistance of GBM to treatment is the extreme heterogeneity of the tumor and, in particular, the presence of cancer stem cells (CSCs) in the population of glioblastoma cells. In this work, we investigated the effect of conditions that reduce the proportion of CSCs in the GBM cell population on the levels of long noncoding RNAs (lincROR and MALAT1) involved in the formation of the phenotype of glioblastoma cancer stem cells. We have shown that culturing under conditions that cause a decrease in cell stemness (when fetal bovine serum is added to the culture medium) affected the content of these transcripts: in the cells of most of the analyzed lines, a decrease in the level of the positive stemness regulator lincROR and an increase in the content of MALAT1 were noted.

Identification of COL1A1/2 Mutations and Fusions With Noncoding RNA Genes in Bizarre Parosteal Osteochondromatous Proliferation (Nora Lesion)

Bizarre parosteal osteochondromatous proliferation (BPOP) (Nora lesion) is a benign bone surface lesion, which most commonly occurs in the digits of young patients and has a high rate of recurrence. Histologically, it is composed of a mixture of disorganized bone, cartilage, and spindle cells in variable proportions and characterized by amorphous “blue bone” mineralization. Recurrent chromosomal abnormalities, including t(1;17)(q32-42;q21-23) and inv(7)(q21.1-22q31.3-32), have been reported in BPOP. However, the exact genes involved in the rearrangements remain unknown. In this study, we analyzed 8 BPOP cases affecting the fingers, toe, ulna, radius, and fibula of 5 female and 3 male patients, aged 5 to 68 years. RNA sequencing of 5 cases identified genetic fusions between COL1A2 and LINC-PINT in 3 cases and COL1A1::MIR29B2CHG fusion in 1, both validated using fluorescence in situ hybridization and reverse transcription (RT)-PCR. The remaining fusion-negative case harbored 3 COL1A1 mutations as revealed by whole-exome sequencing and confirmed using Sanger sequencing. All these genetic alterations were predicted to cause frameshift and/or truncation of COL1A1/2. The chromosomal locations of COL1A2 (7q21.3), LINC-PINT (7q32.3), COL1A1 (17q21.33), and MIR29B2CHG (1q32.2) were consistent with the breakpoints identified in the previous cytogenetic studies. Subsequent screening of 3 BPOPs using fluorescence in situ hybridization identified 1 additional case each with COL1A1 or COL1A2 rearrangement. Our findings are consistent with reported chromosomal abnormalities and implicate the disruption of type I collagen, and perhaps of either noncoding RNA gene as a tumor suppressor, in the tumorigenesis of BPOP. The prevalence and tumorigenic mechanisms of these COL1A1/2 alterations in BPOP require further investigation.

Whole-genome sequencing, annotation, and biological characterization of a novel Siphoviridae phage against multi-drug resistant Propionibacterium acne.

Antibiotics-resistant Propionibacterium acne (P. acne) causes severe acne vulgaris, serious public health, and psychological threat. A new lytic bacteriophage (phage), φPaP11-13, infecting P. acne, was isolated from the sewage management center of Xinqiao Hospital. It can form transparent plaque with diameters of 1.0 ~ 5.0 mm on the double-layer agar plate, indicating a robust lytic ability against its host. Transmission electron microscopy (TEM) showed that φPaP11-13 belonged to the Siphoviridae family (head diameter 60 ± 4.5 nm, tail length 170 ± 6.4 nm, tail width 14 ± 2.4 nm). The one-step growth curve showed the incubation period was 5 h, and the burst size was 26 PFU (plaque-forming unit)/cell. Moreover, it exhibited tolerance over a broad range of pH and temperature ranges but was utterly inactivated by ultraviolet (UV) irradiation for 1 h. The whole-genome sequencing results revealed φPaP11-13 had a linear dsDNA with 29,648 bp length. The G/C content was 54.08%. Non-coding RNA genes and virulence factors were not found. Forty five open reading frames (ORFs) were identified after online annotation. This study reports a novel P. acne phage φPaP11-13, which has a robust lytic ability, no virulence factors, and good stability. The characterization and genomic analysis of φPaP11-13 will develop our understanding of phage biology and diversity and provide a potential arsenal for controlling antibiotics-resistant P. acne-induced severe acne vulgaris.

Activity-Dependent Non-Coding RNA MAPK Interactome of the Human Epileptic Brain.

The human brain has evolved to have extraordinary capabilities, enabling complex behaviors. The uniqueness of the human brain is increasingly posited to be due in part to the functions of primate-specific, including human-specific, long non-coding RNA (lncRNA) genes, systemically less conserved than protein-coding genes in evolution. Patients who have surgery for drug-resistant epilepsy are subjected to extensive electrical recordings of the brain tissue that is subsequently removed in order to treat their epilepsy. Precise localization of brain tissues with distinct electrical properties offers a rare opportunity to explore the effects of brain activity on gene expression. Here, we identified 231 co-regulated, activity-dependent lncRNAs within the human MAPK signaling cascade. Six lncRNAs, four of which were antisense to known protein-coding genes, were further examined because of their high expression and potential impact on the disease phenotype. Using a model of repeated depolarizations in human neuronal-like cells (Sh-SY5Y), we show that five out of six lncRNAs were electrical activity-dependent, with three of four antisense lncRNAs having reciprocal expression patterns relative to their protein-coding gene partners. Some were directly regulated by MAPK signaling, while others effectively downregulated the expression of the protein-coding genes encoded on the opposite strands of their genomic loci. These lncRNAs, therefore, likely contribute to highly evolved and primate-specific human brain regulatory functions that could be therapeutically modulated to treat epilepsy.

Study of LncRNA BANCR Expression in Tumor Tissues and Adjacent Normal Tissues in Gastric Cancer Patients.

Long non-coding RNAs (lncRNAs) have emerged as crucial regulators in various biological processes, including cancer development and progression. This study aimed to investigate the expression differences of the BRAF-activated non-coding RNA (BANCR) gene in GC tissues compared to adjacent normal tissues. The potential diagnostic significance of BANCR in GC was explored, with the aim of improving diagnostic and therapeutic approaches for this global health burden. Tissue samples from 100 gastric cancer (GC) patients were collected, and BANCR expression was analyzed using quantitative real-time PCR. Correlations between BANCR expression and clinicopathological features were assessed, and its biomarker potential was evaluated. In individuals diagnosed with GC, the expression of BANCR was notably elevated in tumor tissues compared to adjacent normal tissues (P < 0.0001). However, the analysis of gene expression data did not demonstrate any statistically significant correlation between elevated BANCR expression and clinicopathological features. According to the ROC analysis, BANCR demonstrated an AUC of 0.6733 (P < 0.0001), with a sensitivity of 73% and a specificity of 45%. However, further evaluation is required to determine its potential as a biomarker (CI 95% = 0.5992 to 0.7473). The observed upregulation of BANCR in GC tissues implies its potential involvement as an oncogenic lncRNA in GC patients. Furthermore, BANCR may serve as a promising biomarker for identification and treatment of GC.

An increase in a long noncoding RNA ANRIL in peripheral plasma is an indicator of stable angina

IntroductionStable angina develops during physical activity or stress, and it is typically an aspect of Coronary Heart Disease (CHD) that can lead to arrhythmia, heart failure and even sudden death. ANRIL, an Antisense Noncoding RNA gene in the INK4 Locus, is associated with multiple disorders including CHD; however, expressional levels of ANRIL in between patients with stable angina and myocardial infarction, one of the acute coronary syndrome, have not been clarified yet. MethodsThe authors enrolled 62 patients with myocardial infarction and 59 with stable angina before primary percutaneous coronary intervention, as well as 48 healthy volunteers. Their peripheral blood was collected for analysis of ANRIL and cardiac troponin I, a traditional diagnostic index of CHD by real-time PCR. ResultsThe data showed that ANRIL is a better diagnostic indicator than cardiac troponin I in patients with stable angina and that the levels of ANRIL are higher in patients with stable angina than those with the myocardial infarction. DiscussionThe levels of ANRIL in peripheral plasma could be used as a good biomarker for stable angina.

The role of the long non-coding RNA gene PVT1 polymorphism in the development of multiple sclerosis

Long non-coding RNA locus PVT1 produces multiple regulatory RNA molecules via alternative splicing, and therefore is a promising object for the development of targeted therapies of different disorders. Performing associationanalysis of genetic variants rs4410871 andrs759648 at the PVT1 locus with multiple sclerosis (MS) in the ethnic group of Russians, Tatars, and Bashkirs from the Volga-Ural region of Russia, we established an association between rs759648*С and MS in the group of Tatars (OR=1.42, PBonf0.046), confirmed by the results of meta-analysisin the three ethnic groups (OR=1.28, P = 0.015). Haplotypic analysis has revealed an association of the rs759648*С/rs4410871*С haplotype with MS in the total study group (χ2=6.083, Рperm =0.012), whilers4410871*C/C + rs759648* C/C combination conferred an increased risk of MS (OR = 2.71, Pperm = 0.027) in this groupaccording to the data yielded by the multilocus analysis. The results of our analysis indicate the significance of the rs759648 polymorphism located near the DNA sequence producing miR-1208. Further studies will help to elucidate the molecular mechanisms of the involvement of PVT1 gene in the etiopathogenesis of MS.

The Genetic and Epigenetic Footprint in Idiopathic Pulmonary Fibrosis and Familial Pulmonary Fibrosis: A State-of-the-Art Review.

Idiopathic pulmonary fibrosis (IPF) is a rare disease of the lung with a largely unknown etiology and a poor prognosis. Intriguingly, forms of familial pulmonary fibrosis (FPF) have long been known and linked to specific genetic mutations. There is little evidence of the possible role of genetics in the etiology of sporadic IPF. We carried out a non-systematic, narrative literature review aimed at describing the main known genetic and epigenetic mechanisms that are involved in the pathogenesis and prognosis of IPF and FPF. In this review, we highlighted the mutations in classical genes associated with FPF, including those encoding for telomerases (TERT, TERC, PARN, RTEL1), which are also found in about 10-20% of cases of sporadic IPF. In addition to the Mendelian forms, mutations in the genes encoding for the surfactant proteins (SFTPC, SFTPA1, SFTPA2, ABCA3) and polymorphisms of genes for the mucin MUC5B and the Toll-interacting protein TOLLIP are other pathways favoring the fibrogenesis that have been thoroughly explored. Moreover, great attention has been paid to the main epigenetic alterations (DNA methylation, histone modification and non-coding RNA gene silencing) that are emerging to play a role in fibrogenesis. Finally, a gaze on the shared mechanisms between cancer and fibrogenesis, and future perspectives on the genetics of pulmonary fibrosis have been analyzed.

Complete Genome Sequence of Antibiotic-Producing Bacillus velezensis H208, Isolated from Ginger Rhizosphere in Laifeng County, China.

The genome of an antibiotic-producing bacterium, Bacillus velezensis H208, was sequenced. Strain H208 was isolated from ginger rhizosphere in Laifeng County, China. The genome consisted of 3,929,792 bp, with a GC content of 46.5%, and contained 3,773 protein-coding genes and 118 noncoding RNA genes.

Increased copy number of imprinted genes in the chromosomal region 20q11-q13.32 is associated with resistance to antitumor agents in cancer cell lines

BackgroundParent of origin-specific allelic expression of imprinted genes is epigenetically controlled. In cancer, imprinted genes undergo both genomic and epigenomic alterations, including frequent copy number changes. We investigated whether copy number loss or gain of imprinted genes in cancer cell lines is associated with response to chemotherapy treatment.ResultsWe analyzed 198 human imprinted genes including protein-coding genes and noncoding RNA genes using data from tumor cell lines from the Cancer Cell Line Encyclopedia and Genomics of Drug Sensitivity in Cancer datasets. We examined whether copy number of the imprinted genes in 35 different genome locations was associated with response to cancer drug treatment. We also analyzed associations of pretreatment expression and DNA methylation of imprinted genes with drug response. Higher copy number of BLCAP, GNAS, NNAT, GNAS-AS1, HM13, MIR296, MIR298, and PSIMCT-1 in the chromosomal region 20q11-q13.32 was associated with resistance to multiple antitumor agents. Increased expression of BLCAP and HM13 was also associated with drug resistance, whereas higher methylation of gene regions of BLCAP, NNAT, SGK2, and GNAS was associated with drug sensitivity. While expression and methylation of imprinted genes in several other chromosomal regions was also associated with drug response and many imprinted genes in different chromosomal locations showed a considerable copy number variation, only imprinted genes at 20q11-q13.32 had a consistent association of their copy number with drug response. Copy number values among the imprinted genes in the 20q11-q13.32 region were strongly correlated. They were also correlated with the copy number of cancer-related non-imprinted genes MYBL2, AURKA, and ZNF217 in that chromosomal region. Expression of genes at 20q11-q13.32 was associated with ex vivo drug response in primary tumor samples from the Beat AML 1.0 acute myeloid leukemia patient cohort. Association of the increased copy number of the 20q11-q13.32 region with drug resistance may be complex and could involve multiple genes.ConclusionsCopy number of imprinted and non-imprinted genes in the chromosomal region 20q11-q13.32 was associated with cancer drug resistance. The genes in this chromosomal region may have a modulating effect on tumor response to chemotherapy.

Structural basis of SNAPc-dependent snRNA transcription initiation by RNA polymerase II

RNA polymerase II (Pol II) carries out transcription of both protein-coding and non-coding genes. Whereas Pol II initiation at protein-coding genes has been studied in detail, Pol II initiation at non-coding genes, such as small nuclear RNA (snRNA) genes, is less well understood at the structural level. Here, we study Pol II initiation at snRNA gene promoters and show that the snRNA-activating protein complex (SNAPc) enables DNA opening and transcription initiation independent of TFIIE and TFIIH in vitro. We then resolve cryo-EM structures of the SNAPc-containing Pol IIpre-initiation complex (PIC) assembled on U1 and U5 snRNA promoters. The core of SNAPc binds two turns of DNA and recognizes the snRNA promoter-specific proximal sequence element (PSE), located upstream of the TATA box-binding protein TBP. Two extensions of SNAPc, called wing-1 and wing-2, bind TFIIA and TFIIB, respectively, explaining how SNAPc directs Pol II to snRNA promoters. Comparison of structures of closed and open promoter complexes elucidates TFIIH-independent DNA opening. These results provide the structural basis of Pol II initiation at non-coding RNA gene promoters.

Estimation of Early Postmortem Interval from Long Noncoding RNA Gene Expression in the Incised Cutaneous Wound: An Experimental Study.

The assessment of alteration of postmortem RNA expression has forensic significance in estimating postmortem interval. To evaluate wound healing progression and the effect of different postmortem intervals, histopathological changes, immunohistochemical matrix metalloproteinase-9 (MMP-9) expression, and long noncoding fatty acid oxidation (lncFAO), RNA expression was assessed in the incised cutaneous wound model. A full-thickness cutaneous wound was inflicted on 75 rats. All 15 rats were sacrificed at different post-infliction intervals (0, 2, 4, 8 and 10 days), and the cutaneous wounds (n = 5) were excised at different postmortem intervals (0, 5, and 24 h after euthanasia). The maximal inflammatory healing stage was detected at day 4 post-infliction, while near complete healing, thick mature collagen deposition was detected at day 10 post-infliction. LncFAO expression was significantly over-expressed with increasing wound age. MMP-9 was detectable on injury day with continuous elevation until 8 days post-wounding, which later decreased. Although histopathological and immunohistochemical examinations within 24 h postmortem did not show any remarkable changes, lncFAO RNA expression showed a significant negative correlation with hours passed since death. The combined use of histopathological changes, immunohistochemical expression of MMP-9, and molecular expression of lncFAO could be appropriate in wound dating verification. Among these factors, lncFAO could be a reliable indicator in postmortem interval estimation.

Improving Infinium MethylationEPIC data processing: re-annotation of enhancers and long noncoding RNA genes and benchmarking of normalization methods

ABSTRACT Illumina Infinium DNA Methylation (5mC) arrays are a popular technology for low-cost, high-throughput, genome-scale measurement of 5mC distribution, especially in cancer and other complex diseases. After the success of its HumanMethylation450 array (450k), Illumina released the MethylationEPIC array (850k) featuring increased coverage of enhancers. Despite the widespread use of 850k, analysis of the corresponding data remains suboptimal: it still relies mostly on Illumina’s default annotation, which underestimates enhancerss and long noncoding RNAs. Results: We have thus developed an approach, based on the ENCODE and LNCipedia databases, which greatly improves upon Illumina’s default annotation of enhancers and long noncoding transcripts. We compared the re-annotated 850k with both 450k and reduced-representation bisulphite sequencing (RRBS), another high-throughput 5mC profiling technology. We found 850k to cover at least three times as many enhancers and long noncoding RNAs as either 450k or RRBS. We further investigated the reproducibility of the three technologies, applying various normalization methods to the 850k data. Most of these methods reduced variability to a level below that of RRBS data. We then used 850k with our new annotation and normalization to profile 5mC changes in breast cancer biopsies. 850k highlighted aberrant enhancer methylation as the predominant feature, in agreement with previous reports. Our study provides an updated processing approach for 850k data, based on refined probe annotation and normalization, allowing for improved analysis of methylation at enhancers and long noncoding RNA genes. Our findings will help to further advance understanding of the DNA methylome in health and disease.

Agarose-Degrading Characteristics of a Deep-Sea Bacterium Vibrio Natriegens WPAGA4 and Its Cold-Adapted GH50 Agarase Aga3420.

Up until now, the characterizations of GH50 agarases from Vibrio species have rarely been reported compared to GH16 agarases. In this study, a deep-sea strain, WPAGA4, was isolated and identified as Vibrio natriegens due to the maximum similarity of its 16S rRNA gene sequence, the values of its average nucleotide identity, and through digital DNA-DNA hybridization. Two circular chromosomes in V. natriegens WPAGA4 were assembled. A total of 4561 coding genes, 37 rRNA, 131 tRNA, and 59 other non-coding RNA genes were predicted in the genome of V. natriegens WPAGA4. An agarase gene belonging to the GH50 family was annotated in the genome sequence and expressed in E. coli cells. The optimum temperature and pH of the recombinant Aga3420 (rAga3420) were 40 °C and 7.0, respectively. Neoagarobiose (NA2) was the only product during the degradation process of agarose by rAga3420. rAga3420 had a favorable stability following incubation at 10-30 °C for 50 min. The Km, Vmax, and kcat values of rAga3420 were 2.8 mg/mL, 78.1 U/mg, and 376.9 s-1, respectively. rAga3420 displayed cold-adapted properties as 59.7% and 41.2% of the relative activity remained at 10 3 °C and 0 °C, respectively. This property ensured V. natriegens WPAGA4 could degrade and metabolize the agarose in cold deep-sea environments and enables rAga3420 to be an appropriate industrial enzyme for NA2 production, with industrial potential in medical and cosmetic fields.

Population specific methylome remodeling in high and low elevation populations of Indian west Himalayan Arabidopsis thaliana in response to elevated CO2

Unravelling the plants responses towards rising atmospheric [CO 2 ] have largely been explored in genetical and morpho-physiological contexts. However, epigenetic factors including DNA methylation which play critical roles in adaptation are largely unexplored. Here we investigated the methylome, transcriptome and morpho-physiological responses of the two Arabidopsis thaliana populations evolved at high (3400 m amsl) and low elevation (700 m amsl) zones to elevated [CO 2 ] (e[CO 2 ]). We show that depending on the origin of the population, there were local level variations in DNA methylation when exposed to e[CO 2 ], but global methyl cytosine (mC) content remained unchanged. Further, there was loss of methylation and more protein coding genes were differentially methylated than transposable elements and non-coding RNA genes, in both the populations but more so in the low elevation one. The differentially methylated genes of the two populations belonged to distinct functional categories. More genes related to methylation machinery were down-regulated in the high elevation population than the low elevation one. Although there was no correlation between methylation and gene expression at the global scale but a few genes exhibited methylation dependent expression level. Finally the hyper- and hypo-methylated status of a few genes due to e[CO 2 ] treatment were validated. Overall, our data suggested the two populations responded differently towards e[CO 2 ] with respect to methylome remodeling, phenotypic and molecular plasticity. However, methylome remodeling and molecular plasticity were more prominent in the low elevation population. Understanding the evolution of epigenetic response towards e[CO 2 ] may help in future crop improvement strategies. • Elevated CO 2 affects DNA methylation remodeling in Arabidopsis thaliana depending on origin of the population. • Methylome remodeling and molecular plasticity of the low elevation population were more than the higher elevation one. • More protein coding genes were differentially methylated in response to elevated CO 2. • Elevated CO 2 induced global gene expression and methylation is not co-linear.

Spermine synthase (SMS) serves as a prognostic biomarker in head and neck squamous cell carcinoma: a bioinformatics analysis

Head and neck squamous cell carcinoma (HNSC) is an aggressive type of cancer that lacks early detection, and therefore, has a low 5-year survival rate. The spermine synthase (SMS) gene has been shown to be associated with Snyder-Robinson syndrome and poor prognosis of multiple cancers; however, its regulatory role in HNSC has never been investigated. Therefore, we explored the potential predictive value of SMS in HNSC. We explored the association between SMS expression and clinicopathological parameters of HNSC patients by using data from The Cancer Genome Atlas datasets (TCGA). The prognostic value of SMS was evaluated using the Kaplan-Meier plotter, Gene Expression Profiling Interactive Analysis (GEPIA) 2 and univariate and multivariate Cox regression analyses. We further used gene set enrichment analysis (GESA) to investigate the potential roles of SMS in HNSC prognosis and Tumor Immunity Estimation Resource 2.0 (TIMER2.0) to analyze the correlation between immune cell infiltration and SMS expression. Finally, starBase was used to screen out prognosis-associated non-coding RNA genes to constructed the competing endogenous RNA (ceRNA) network. Co-expression and survival analyses were used to identify the ceRNA network's effect on HNSC prognosis. We found that SMS expression was increased in HNSC compared with normal tissues (P<0.05). In addition, SMS expression was associated with tumor grade (P=0.006), N stage (P=0.001), and prognosis. Survival analysis revealed that high expression of SMS showed worse overall survival (OS) (HR =1.4, P=0.01) and worse disease-free survival (DFS) (HR =1.5, P=0.014). Multivariate Cox analysis further supported the prognostic value of SMS in HNSC (HR =1.006636, P=0.0056). GESA showed that SMS was involved in metabolism- and immune-related pathways. The immune infiltration analyses results showed a decrease in the landscape of immune cell infiltration with high SMS expression and SMS deletion in HNSC. Finally, a ceRNA network (SMS/hsa-miR-23b-3p/KTN1-AS1 and VPS9D1-AS axis) was constructed based on the co-expression and survival analyses in HNSC. Our findings first revealed that SMS functioned as a potential prognostic biomarker and provide insights into the molecular mechanisms of its function in HNSC. The use of SMS may be powerful for determining worse prognosis HNSC patients.

Epigenetic control of chromosome-associated lncRNA genes essential for replication and stability

ASARs are long noncoding RNA genes that control replication timing of entire human chromosomes in cis. The three known ASAR genes are located on human chromosomes 6 and 15, and are essential for chromosome integrity. To identify ASARs on all human chromosomes we utilize a set of distinctive ASAR characteristics that allow for the identification of hundreds of autosomal loci with epigenetically controlled, allele-restricted behavior in expression and replication timing of coding and noncoding genes, and is distinct from genomic imprinting. Disruption of noncoding RNA genes at five of five tested loci result in chromosome-wide delayed replication and chromosomal instability, validating their ASAR activity. In addition to the three known essential cis-acting chromosomal loci, origins, centromeres, and telomeres, we propose that all mammalian chromosomes also contain “Inactivation/Stability Centers” that display allele-restricted epigenetic regulation of protein coding and noncoding ASAR genes that are essential for replication and stability of each chromosome.