Housekeeping protein-coding genes interrogated with tissue and individual variations.

Housekeeping protein-coding genes are stably expressed genes in cells and tissues that are thought to be engaged in fundamental cellular biological functions. They are often utilized as normalization references in molecular biology research and are especially important in integrated bioinformatic investigations. Prior studies have examined human housekeeping protein-coding genes by analyzing various gene expression datasets. The inclusion of different tissue types significantly impacted the discovery of housekeeping genes. In this report, we investigated particularly individual human subject expression differences in protein-coding genes across different tissue types. We used GTEx V8 gene expression datasets obtained from more than 16,000 human normal tissue samples. Furthermore, the Gini index is utilized to investigate the expression variations of protein-coding genes between tissue and individual donor subjects. Housekeeping protein-coding genes found using Gini index profiles may vary depending on the tissue subtypes investigated, particularly given the diverse sample size collections across the GTEx tissue subtypes. We subsequently selected major tissues and identified subsets of housekeeping genes with stable expression levels among human donors within those tissues. In this work, we provide alternative sets of housekeeping protein-coding genes that show more consistent expression patterns in human subjects across major solid organs. Weblink: https://hpsv.ibms.sinica.edu.tw .

Synergistic Effect of Metformin and Lansoprazole Against Gastric Cancer through Growth Inhibition.

Cancer has been linked to metabolic disorders and diverse gene mutations. Metformin, which is widely used to treat type 2 diabetes, inhibits the growth of cancer cells in animal models. Here we investigated the effects of metformin on human gastric cancer cell lines. We also investigated the synergistic anticancer effect of metformin and proton pump inhibitors. Lansoprazole, a proton pump inhibitor, is effective for treating gastroesophageal reflux disease. Our results revealed that metformin and lansoprazole can significantly inhibit cancer cell growth in a dose-dependent manner by suppressing cell cycle progression and inducing apoptosis. Low concentrations of metformin and lansoprazole have a synergistic effect on AGS cell growth inhibition. In summary, our findings suggest a new and safe treatment protocol for treating stomach cancers.

TEx-MST: tissue expression profiles of MANE select transcripts.

Recently, a new reference transcript dataset [Matched Annotation from the NCBI and EMBL-EBI (MANE) select] was released by NCBI and EMBL-EBI to make available a new unified representative transcript for human protein-coding genes. While the main purpose of MANE project is to provide a harmonized gene and transcript information standard, there is no explicit tissue expression information about these MANE select transcripts. In this report, we tried to provide useful expression profiles of MANE select transcripts in various normal human tissues to allow further interrogation of their molecular modulations and functional significance. We obtained the new V9 transcript expression dataset from the Genotype-Tissue Expression (GTEx) web portal. This new GTEx dataset, based on a long-read sequencing platform, affords better assessment of the expression of alternative spliced transcripts. This tissue expression profiles of MANE select transcripts (TEx-MST) database not only provides the basic information of MANE select transcripts but also tissue expression profiles on alternative transcripts in protein-coding genes. Users can initiate the interrogation by gene symbol searches or by browsing the MANE genes with various criteria (such as genome locations or expression rankings). We further utilized the GENCODE biotype feature to identify the top-ranked protein-coding transcripts by choosing the most expressed protein-coding transcripts from GTEx datasets (both V8 and V9 datasets). In summary, there are 18 083 genes matched between MANE and GTEx. Among them, 13 245 MANE select transcripts matched with the top-ranked protein-coding transcripts in GTEx V9 dataset, which underlined the dominate expression of MANE select transcripts. This TEx-MST web bioinformatic database provides a visualized user interface for the normal tissue expression patterns of MANE select transcripts using the newly released GTEx dataset. Database URL: TEx-MST is available at https://texmst.ibms.sinica.edu.tw/.

Hepatitis B virus virion secretion is a CRM1-spike-mediated late event.

BACKGROUND: Hepatitis B virus (HBV) is a major human pathogen worldwide. To date, there is no curative treatment for chronic hepatitis B. The mechanism of virion secretion remains to be investigated. Previously, we found that nuclear export of HBc particles can be facilitated via two CRM1-specific nuclear export signals (NES) at the spike tip. METHODS: In this study, we used site-directed mutagenesis at the CRM1 NES, as well as treatment with CRM1 inhibitors at a low concentration, or CRM1-specific shRNA knockdown, in HBV-producing cell culture, and measured the secretion of various HBV viral and subviral particles via a native agarose gel electrophoresis assay. Separated HBV particles were characterized by Western blot analysis, and their genomic DNA contents were measured by Southern blot analysis. Secreted extracellular particles were compared with intracellular HBc capsids for DNA synthesis and capsid formation. Virion secretion and the in vivo interactions among HBc capsids, CRM1 and microtubules, were examined by proximity ligation assay, immunofluorescence microscopy, and nocodazole treatment. RESULTS: We report here that the tip of spike of HBV core (HBc) particles (capsids) contains a complex sensor for secretion of both HBV virions and naked capsids. HBV virion secretion is closely associated with HBc nuclear export in a CRM1-dependent manner. At the conformationally flexible spike tips of HBc particles, NES motifs overlap extensively with motifs important for secretion of HBV virions and naked capsids. CONCLUSIONS: We provided experimental evidence that virions and naked capsids can egress via two distinct, yet overlapping, pathways. Unlike the secretion of naked capsids, HBV virion secretion is highly CRM1- and microtubule-dependent. CRM1 is well known for its involvement in nuclear transport in literature. To our knowledge, this is the first report that CRM1 is required for virion secretion. CRM1 inhibitors could be a promising therapeutic candidate for chronic HBV patients in clinical medicine.

Dominant transcript expression profiles of human protein-coding genes interrogated with GTEx dataset.

The discovery and quantification of mRNA transcripts using short-read next-generation sequencing (NGS) data is a complicated task. There are far more alternative mRNA transcripts expressed by human genes than can be identified from NGS transcriptome data and various bioinformatic pipelines, while the numbers of annotated human protein-coding genes has gradually declined in recent years. It is essential to learn more about the thorough tissue expression profiles of alternative transcripts in order to obtain their molecular modulations and actual functional significance. In this report, we present a bioinformatic database for interrogating the representative tissue of human protein-coding transcripts. The database allows researchers to visually explore the top-ranked transcript expression profiles in particular tissue types. Most transcripts of protein-coding genes were found to have certain tissue expression patterns. This observation demonstrated that many alternative transcripts were particularly modulated in different cell types. This user-friendly tool visually represents transcript expression profiles in a tissue-specific manner. Identification of tissue specific protein-coding genes and transcripts is a substantial advance towards interpreting their biological functions and further functional genomics studies.

Circular RNAs Modulate Cancer Hallmark and Molecular Pathways to Support Cancer Progression and Metastasis.

Circular RNAs (circRNAs) are noncoding products of backsplicing of pre-mRNAs which have been established to possess potent biological functions. Dysregulated circRNA expression has been linked to diseases including different types of cancer. Cancer progression is known to result from the dysregulation of several molecular mechanisms responsible for the maintenance of cellular and tissue homeostasis. The dysregulation of these processes is defined as cancer hallmarks, and the molecular pathways implicated in them are regarded as the targets of therapeutic interference. In this review, we summarize the literature on the investigation of circRNAs implicated in cancer hallmark molecular signaling. First, we present general information on the properties of circRNAs, such as their biogenesis and degradation mechanisms, as well as their basic molecular functions. Subsequently, we summarize the roles of circRNAs in the framework of each cancer hallmark and finally discuss the potential as therapeutic targets.

Strategic Decoy Peptides Interfere with MSI1/AGO2 Interaction to Elicit Tumor Suppression Effects.

Peptide drugs that target protein-protein interactions have attracted mounting research efforts towards clinical developments over the past decades. Increasing reports have indicated that expression of Musashi 1 (MSI1) is tightly correlated to high grade of cancers as well as enrichment of cancer stem cells. Treatment failure in malignant tumors glioblastoma multiform (GBM) had also been correlated to CSC-regulating properties of MSI1. It is thus imperative to develop new therapeutics that could effectively improve current regimens used in clinics. MSI1 and AGO2 are two emerging oncogenic molecules that both contribute to GBM tumorigenesis through mRNA regulation of targets involved in apoptosis and cell cycle. In this study, we designed peptide arrays covering the C-terminus of MSI1 and identified two peptides (Pep#11 and Pep#26) that could specifically interfere with the binding with AGO2. Our Biacore analyses ascertained binding between the identified peptides and AGO2. Recombinant reporter system Gaussian luciferase and fluorescent bioconjugate techniques were employed to determine biological functions and pharmacokinetic characteristics of these two peptides. Our data suggested that Pep#11 and Pep#26 could function as decoy peptides by mimicking the interaction function of MSI1 with its binding partner AGO2 in vitro and in vivo. Further experiments using GMB animal models corroborated the ability of Pep#11 and Pep#26 in disrupting MSI1/AGO2 interaction and consequently anti-tumorigenicity and prolonged survival rates. These striking therapeutic efficacies orchestrated by the synthetic peptides were attributed to the decoy function to C-terminal MSI1, especially in malignant brain tumors and glioblastoma.

Musashi-1 Regulates MIF1-Mediated M2 Macrophage Polarization in Promoting Glioblastoma Progression.

Glioblastoma (GBM) is the most malignant brain tumor which is characterized by high proliferation and migration capacity. The poor survival rate has been attributed to limitations of the current standard therapies. The search for novel biological targets that can effectively hamper tumor progression remains extremely challenging. Previous studies indicated that tumor-associated macrophages (TAMs) are the abundant elements in the tumor microenvironment that are closely implicated in glioma progression and tumor pathogenesis. M2 type TAMs are immunosuppressive and promote GBM proliferation. RNA-binding protein Musashi-1 (MSI1) has recently been identified as a marker of neural stem/progenitor cells, and its high expression has been shown to correlate with the growth of GBM. Nevertheless, the relationship between MSI1 and TAMs in GBM is still unknown. Thus, in our present study, we aimed to investigate the molecular interplay between MSI1 and TAMs in contributing to GBM tumorigenesis. Our data revealed that the secretion of macrophage inhibitory factor 1 (MIF1) is significantly upregulated by MSI1 overexpression in vitro. Importantly, M2 surface markers of THP-1-derived macrophages were induced by recombinant MIF1 and reduced by using MIF1 inhibitor (S,R)-3-(4-hHydroxyphenyl)-4,5-dihydro-5-isoxazole acetic acid (ISO-1). Furthermore, GBM tumor model data suggested that the tumor growth, MIF1 expression and M2 macrophage population were significantly downregulated when MSI1 expression was silenced in vivo. Collectively, our findings identified a novel role of MSI1 in the secretion of MIF1 and the consequent polarization of macrophages into the M2 phenotype in promoting GBM tumor progression.

Chemical and Biochemical Aspects of Molecular Hydrogen in Treating Kawasaki Disease and COVID-19.

Kawasaki disease (KD) is a systemic vasculitis and is the most commonly acquired heart disease among children in many countries, which was first reported 50 years ago in Japan. The 2019 coronavirus disease (COVID-19, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)) has been a pandemic in most of the world since 2020, and since late 2019 in China. Kawasaki-like disease caused by COVID-19 shares some symptoms with KD, referred to as multisystem inflammatory syndrome in children, and has been reported in the United States, Italy, France, England, and other areas of Europe, with an almost 6-10 times or more increase compared with previous years of KD prevalence. Hydrogen gas is a stable and efficient antioxidant, which has a positive effect on oxidative damage, inflammation, cell apoptosis, and abnormal blood vessel inflammation. This review reports the chemical and biochemical aspects of hydrogen gas inhalation in treating KD and COVID-19.

Top-ranked expressed gene transcripts of human protein-coding genes investigated with GTEx dataset.

With considerable accumulation of RNA-Seq transcriptome data, we have extended our understanding about protein-coding gene transcript compositions. However, alternatively compounded patterns of human protein-coding gene transcripts would complicate gene expression data processing and interpretation. It is essential to exhaustively interrogate complex mRNA isoforms of protein-coding genes with an unified data resource. In order to investigate representative mRNA transcript isoforms to be utilized as transcriptome analysis references, we utilized GTEx data to establish a top-ranked transcript isoform expression data resource for human protein-coding genes. Distinctive tissue specific expression profiles and modulations could be observed for individual top-ranked transcripts of protein-coding genes. Protein-coding transcripts or genes do occupy much higher expression fraction in transcriptome data. In addition, top-ranked transcripts are the dominantly expressed ones in various normal tissues. Intriguingly, some of the top-ranked transcripts are noncoding splicing isoforms, which imply diverse gene regulation mechanisms. Comprehensive investigation on the tissue expression patterns of top-ranked transcript isoforms is crucial. Thus, we established a web tool to examine top-ranked transcript isoforms in various human normal tissue types, which provides concise transcript information and easy-to-use graphical user interfaces. Investigation of top-ranked transcript isoforms would contribute understanding on the functional significance of distinctive alternatively spliced transcript isoforms.

miR-TV: an interactive microRNA Target Viewer for microRNA and target gene expression interrogation for human cancer studies.

MicroRNAs (miRNAs) have been identified in many organisms, and they are essential for gene expression regulation in many critical cellular processes. The expression levels of these genes and miRNAs are closely associated with the progression of diseases such as cancers. Furthermore, survival analysis is a significant indicator for evaluating the criticality of these cellular processes in cancer progression. We established a web tool, miRNA Target Viewer (miR-TV), which integrates 5p-arm and 3p-arm miRNA expression profiles, mRNA target gene expression levels in healthy and cancer populations, and clinical data of cancer patients and their survival information. The developed miR-TV obtained miRNA-seq, mRNA-seq and clinical data from the Cancer Genome Atlas and potential miRNA target gene predictions from miRDB, targetScan and miRanda. The data presentation was implemented using the D3 javascript toolkit. The D3 toolkit is frequently used to provide an easy-to-use interactive interface. Our miR-TV provides a user-friendly and interactive interface, which can be beneficial for biomedical researchers to freely interrogate miRNA expression information and their potential target genes. We believe that such a data visualization bioinformatics tool is excellent for obtaining information from massive biological data. Database URL: http://mirtv.ibms.sinica.edu.tw.

Musashi-1 promotes stress-induced tumor progression through recruitment of AGO2.

Carcinomatous progression and recurrence are the main therapeutic challenges frequently faced by patients with refractory tumors. However, the underlined molecular mechanism remains obscure. Methods: We found Musashi-1 (MSI1) transported into cytosol under stress condition by confocal microscopy and cell fractionation. Argonaute 2 (AGO2) was then identified as a cytosolic binding partner of MSI1 by Mass Spectrametry, immunoprecipitation, and recombinant protein pull-down assay. We used RNA-IP to determine the MSI1/AGO2 associated regions on downstream target mRNAs. Finally, we overexpressed C-terminus of MSI1 to disrupt endogenous MSI1/AGO2 interaction and confirm it effects on tmor progression. Results: Malignant tumors exhibit elevated level of cytosolic Musashi-1 (MSI1), which translocates into cytosol in response to stress and promote tumor progression. Cytosolic MSI1 forms a complex with AGO2 and stabilize or destabilize its target mRNAs by respectively binding to their 3´ untranslated region or coding domain sequence. Both MSI1 translocation and MSI1/AGO2 binding are essential for promoting tumor progression. Blocking MSI1 shuttling by either chemical inhibition or point mutation attenuates the growth of GBM-xenografts in mice. Importantly, overexpression of the C-terminus of MSI1 disrupts endogenous MSI1/AGO2 interaction and effectively reduces stress-induced tumor progression. Conclusion: Our findings highlight novel molecular functions of MSI1 during stress-induced carcinomatous recurrence, and suggest a new therapeutic strategy for refractory malignancies by targeting MSI1 translocation and its interaction with AGOs.

Overlapping protein-coding genes in human genome and their coincidental expression in tissues.

The completion of human genome sequences and the advancement of next-generation sequencing technologies have engendered a clear understanding of all human genes. Overlapping genes are usually observed in compact genomes, such as those of bacteria and viruses. Notably, overlapping protein-coding genes do exist in human genome sequences. Accordingly, we used the current Ensembl gene annotations to identify overlapping human protein-coding genes. We analysed 19,200 well-annotated protein-coding genes and determined that 4,951 protein-coding genes overlapped with their adjacent genes. Approximately a quarter of all human protein-coding genes were overlapping genes. We observed different clusters of overlapping protein-coding genes, ranging from two genes (paired overlapping genes) to 22 genes. We also divided the paired overlapping protein-coding gene groups into four subtypes. We found that the divergent overlapping gene subtype had a stronger expression association than did the subtypes of 5'-tandem overlapping and 3'-tandem overlapping genes. The majority of paired overlapping genes exhibited comparable coincidental tissue expression profiles; however, a few overlapping gene pairs displayed distinctive tissue expression association patterns. In summary, we have carefully examined the genomic features and distributions about human overlapping protein-coding genes and found coincidental expression in tissues for most overlapping protein-coding genes.

MicroRNA let-7-TGFBR3 signalling regulates cardiomyocyte apoptosis after infarction.

BACKGROUND: Myocardial infarction (MI) is a life-threatening disease, often leading to heart failure. Defining therapeutic targets at an early time point is important to prevent heart failure. METHODS: MicroRNA screening was performed at early time points after MI using paired samples isolated from the infarcted and remote myocardium of pigs. We also examined the microRNA expression in plasma of MI patients and pigs. For mechanistic studies, AAV9-mediated microRNA knockdown and overexpression were administrated in mice undergoing MI. FINDINGS: MicroRNAs let-7a and let-7f were significantly downregulated in the infarct area within 24 h post-MI in pigs. We also observed a reduction of let-7a and let-7f in plasma of MI patients and pigs. Inhibition of let-7 exacerbated cardiomyocyte apoptosis, induced a cardiac hypertrophic phenotype, and resulted in worsened left ventricular ejection fraction. In contrast, ectopic let-7 overexpression significantly reduced those phenotypes and improved heart function. We then identified TGFBR3 as a target of let-7, and found that induction of Tgfbr3 in cardiomyocytes caused apoptosis, likely through p38 MAPK activation. Finally, we showed that the plasma TGFBR3 level was elevated after MI in plasma of MI patients and pigs. INTERPRETATION: Together, we conclude that the let-7-Tgfbr3-p38 MAPK signalling plays an important role in cardiomyocyte apoptosis after MI. Furthermore, microRNA let-7 and Tgfbr3 may serve as therapeutic targets and biomarkers for myocardial damage. FUND: Ministry of Science and Technology, National Health Research Institutes, Academia Sinica Program for Translational Innovation of Biopharmaceutical Development-Technology Supporting Platform Axis, Thematic Research Program and the Summit Research Program, Taiwan.

Covariate-adjusted heatmaps for visualizing biological data via correlation decomposition.

Motivation: Heatmap is a popular visualization technique in biology and related fields. In this study, we extend heatmaps within the framework of matrix visualization (MV) by incorporating a covariate adjustment process through the estimation of conditional correlations. MV can explore the embedded information structure of high-dimensional large-scale datasets effectively without dimension reduction. The benefit of the proposed covariate-adjusted heatmap is in the exploration of conditional association structures among the subjects or variables that cannot be done with conventional MV. Results: For adjustment of a discrete covariate, the conditional correlation is estimated by the within and between analysis. This procedure decomposes a correlation matrix into the within- and between-component matrices. The contribution of the covariate effects can then be assessed through the relative structure of the between-component to the original correlation matrix while the within-component acts as a residual. When a covariate is of continuous nature, the conditional correlation is equivalent to the partial correlation under the assumption of a joint normal distribution. A test is then employed to identify the variable pairs which possess the most significant differences at varying levels of correlation before and after a covariate adjustment. In addition, a z-score significance map is constructed to visualize these results. A simulation and three biological datasets are employed to illustrate the power and versatility of our proposed method. Availability and implementation: GAP is available to readers and is free to non-commercial applications. The installation instructions, the user's manual, and the detailed tutorials can be found at http://gap.stat.sinica.edu.tw/Software/GAP. Supplementary information: Supplementary Data are available at Bioinformatics online.

A Comprehensive Analysis of Transcript-Supported De Novo Genes in Saccharomyces sensu stricto Yeasts.

Novel genes arising from random DNA sequences (de novo genes) have been suggested to be widespread in the genomes of different organisms. However, our knowledge about the origin and evolution of de novo genes is still limited. To systematically understand the general features of de novo genes, we established a robust pipeline to analyze >20,000 transcript-supported coding sequences (CDSs) from the budding yeast Saccharomyces cerevisiae. Our analysis pipeline combined phylogeny, synteny, and sequence alignment information to identify possible orthologs across 20 Saccharomycetaceae yeasts and discovered 4,340 S. cerevisiae-specific de novo genes and 8,871 S. sensu stricto-specific de novo genes. We further combine information on CDS positions and transcript structures to show that >65% of de novo genes arose from transcript isoforms of ancient genes, especially in the upstream and internal regions of ancient genes. Fourteen identified de novo genes with high transcript levels were chosen to verify their protein expressions. Ten of them, including eight transcript isoform-associated CDSs, showed translation signals and five proteins exhibited specific cytosolic localizations. Our results suggest that de novo genes frequently arise in the S. sensu stricto complex and have the potential to be quickly integrated into ancient cellular network.

Identification, chromosomal arrangements and expression analyses of the evolutionarily conserved prmt1 gene in chicken in comparison with its vertebrate paralogue prmt8.

Nine protein arginine methyltransferases (PRMTs) are conserved in mammals and fish. Among these, PRMT1 is the major type I PRMT for asymmetric dimethylarginine (ADMA) formation and is the most conserved and widely distributed one. Two chicken prmt1 splicing variants were assembled and confirmed by RT-PCR experiments. However, only two scaffolds containing single separate prmt1 exon with high GC contents are present in the current chicken genome assembly. Besides, prmt1 exons are scattered in separate small scaffolds in most avian species. Complete prmt1 gene has only been predicted from two falcon species with few neighboring genes. Crocodilians are considered close to the common ancestor shared by crocodilians and birds. The gene arrangements around prmt1 in American alligator are different from that in birds but are largely conserved in human. Orthologues of genes in a large segment of human chromosomal 19 around PRMT1 are missing or not assigned to the current chicken chromosomes. In comparison, prmt8, the prmt1 paralogue, is on chicken chromosome 1 with the gene arrangements downstream of prmt8 highly conserved in birds, crocodilians, and human. However, the ones upstream vary greatly in birds. Biochemically, we found that though prmt1 transcripts were detected, limited or none PRMT1 protein was present in chicken tissues. Moreover, a much higher level of PRMT8 protein was detected in chicken brain than in mouse brain. While PRMT8 is brain specific in other vertebrate species studied, low level of PRMT8 was present in chicken but not mouse liver and muscle. We also showed that the ADMA level in chicken was similar to that in mouse. This study provides the critical information of chicken PRMT1 and PRMT8 for future analyses of the function of protein arginine methyltransferases in birds.

Urine miR-21-5p as a potential non-invasive biomarker for gastric cancer.

Many reports have implicated that microRNAs involve in cancer development and progression, such as miR-155 in breast cancers and miR-196 in gastric cancers. Furthermore, microRNAs are more stable than typical protein-coding gene mRNAs in varieties of clinical samples including body fluids. This suggests that they are potentially valuable biomarkers for cancer monitoring. In this study, we have used urine samples of gastric cancer patients to demonstrate the feasibility of urine microRNAs for gastric cancer detection. Urine samples of gastric cancer patients were extracted for total RNA, which were examined for the expression of miR-21-5p using quantitative stem-loop PCR. Our results demonstrated that miR-21-5p could be detected in small amounts of urine samples with good stability, and the expression levels of miR-21-5p were reduced following surgical removal of gastric cancer tissues. These results implicate that urine miR-21-5p could be utilized as a novel non-invasive biomarker of gastric cancer detection and monitoring.

Visual Display of 5p-arm and 3p-arm miRNA Expression with a Mobile Application.

MicroRNAs (miRNAs) play important roles in human cancers. In previous studies, we have demonstrated that both 5p-arm and 3p-arm of mature miRNAs could be expressed from the same precursor and we further interrogated the 5p-arm and 3p-arm miRNA expression with a comprehensive arm feature annotation list. To assist biologists to visualize the differential 5p-arm and 3p-arm miRNA expression patterns, we utilized a user-friendly mobile App to display. The Cancer Genome Atlas (TCGA) miRNA-Seq expression information. We have collected over 4,500 miRNA-Seq datasets from 15 TCGA cancer types and further processed them with the 5p-arm and 3p-arm annotation analysis pipeline. In order to be displayed with the RNA-Seq Viewer App, annotated 5p-arm and 3p-arm miRNA expression information and miRNA gene loci information were converted into SQLite tables. In this distinct application, for any given miRNA gene, 5p-arm miRNA is illustrated on the top of chromosome ideogram and 3p-arm miRNA is illustrated on the bottom of chromosome ideogram. Users can then easily interrogate the differentially 5p-arm/3p-arm expressed miRNAs with their mobile devices. This study demonstrates the feasibility and utility of RNA-Seq Viewer App in addition to mRNA-Seq data visualization.

Identification of lncRNA functions in lung cancer based on associated protein-protein interaction modules.

Long non-coding RNAs (lncRNAs) have been found to play important roles in various biological processes; however, many of their functions remain unclear. In this study, we present a novel approach to identify the lncRNA-associated protein-protein interaction (PPI) modules and ascertain their functions in human lung squamous cell carcinoma. We collected lncRNA and mRNA expression profiles of lung squamous cell carcinoma from The Cancer Genome Atlas. To identify the lncRNA-associated PPI modules, lncRNA-mRNA co-expression networks were first constructed based on the mutual ranks of expression correlations. Next, we examined whether the co-expressed mRNAs of a specific lncRNA were closely connected by PPIs. For this, a significantly connected mRNA set was considered to be the lncRNA-associated PPI module. Finally, the prospective functions of a lncRNA was inferred using Gene Ontology enrichment analysis on the associated module. We found that lncRNA-associated PPI modules were subtype-dependent and each subtype had unique molecular mechanisms. In addition, antisense lncRNAs and sense genes tended to be functionally associated. Our results might provide new directions for understanding lncRNA regulations in lung cancer. The analysis pipeline was implemented in a web tool, available at http://lncin.ym.edu.tw/.