RNA-Seq: Improving Our Understanding of Retinal Biology and Disease

125I Radioactive Particles Drive Protective Autophagy in Hepatocellular Carcinoma by Upregulating ATG9B

BACKGROUND AND AIMS

¹²⁵I radioactive particles implantation have demonstrated efficacy in eradicating hepatocellular carcinoma (HCC). However, progressive resistance of HCC to ¹²⁵I radioactive particles has limited its wide clinical application.

METHODS

We investigated the cellular responses to ¹²⁵I radioactive particles treatment and autophagy-related 9B (ATG9B) silencing in HCC cell lines and Hep3B xenografted tumor model using Cell Counting Kit-8 reagent, western blotting, immunofluorescence, flow cytometry, transmission electron microscopy and immunohistochemistry.

RESULTS

In this study, we demonstrated that ¹²⁵I radioactive particles induced cell apoptosis and protective autophagy of HCC in vitro and in vivo. Inhibition of autophagy enhanced the radiosensitivity of HCC to ¹²⁵I radioactive particles. Moreover, ¹²⁵I radioactive particles induced autophagy by upregulating ATG9B, with increased expression level of LC3B and decreased expression level of p62. Furthermore, ATG9B silencing downregulated LC3B expression and upregulated p62 expression and enhanced radiosensitivity of HCC to ¹²⁵I radioactive particles in vitro and in vivo.

CONCLUSIONS

Inhibition of ATG9B enhanced the antitumor effects of ¹²⁵I particle radiation against HCC in vitro and in vivo. Our findings suggest that ¹²⁵I particle radiation plus chloroquine or/and the ATG9B inhibitor may be a novel therapeutic strategy for HCC.

Transcriptomic Responses of Human Retinal Vascular Endothelial Cells to Inflammatory Cytokines

Purpose

Molecular profiling of human retinal endothelial cells provides opportunities to understand the roles of this cell population in maintenance of the blood-ocular barrier, and its involvements in diverse retinal vasculopathies. We aimed to generate a transcriptome of human retinal endothelial cells in the unstimulated state, and following treatment with inflammatory cytokines linked to cell dysfunction.

Methods

Endothelial cells were isolated from retinae of five human cadaveric donors, and treated for 60 minutes and 24 hours with interleukin-1β or tumor necrosis factor-α, or exposed to medium alone for the same intervals. Expression of intercellular adhesion molecule-1 was measured by RT-qPCR to confirm cytokine-induced activation of the cells. RNA was sequenced on the Illumina NovaSeq 6000 platform. Reads were aligned to the human GRCh38 genome, and reads that aligned to Ensembl-annotated genes were counted. Quality control of sequencing was performed with FastQC, and sequences were classified by Kraken.

Results

A human retinal endothelial cell RNA-sequencing dataset with mean of 99% reads aligned to the human genome was produced as raw RNA sequence data (FASTQ files) and processed read data (XLSX files). Multidimensional scaling analysis showed a strong donor effect, which was readily controlled by ComBat.

Conclusions

Our dataset may be useful for human retinal endothelial cell transcriptomic assemblies, functional gene annotating and/or gene expression and enrichment analyses, as well as cross-dataset harmonization.

Translational Relevance

The molecular profile of the human retinal endothelium is a source of candidate biologic targets for retinal vasculopathies.

Transcriptome analysis reveals the mechanism of stromal cell-derived factor-1 and exendin-4 synergistically promoted periodontal ligament stem cells osteogenic differentiation

Stromal cell-derived factor-1 (SDF-1) and Exendin-4 (EX-4) play beneficial roles in promoting periodontal ligament stem cells (PDLSCs) osteogenic differentiation, while the detailed mechanism has not been clarified. In this study, we aimed to evaluate the biological mechanism of SDF-1 and EX-4 alone or synergistic application in regulating PDLSCs differentiation by RNA-sequencing (RNA-seq). A total of 110, 116 and 109 differentially expressed genes (DEGs) were generated in osteogenic medium induced PDLSCs treated by SDF-1, EX-4, and SDF-1+EX-4, respectively. The DEGs in SDF-1 group were enriched in signal transduction related signaling pathways; the DEGs in EX-4 group were enriched in metabolism and biosynthesis-related pathways; and the DEGs generated in SDF-1+EX-4 group were mainly enriched in RNA polymerase II transcription, cell differentiation, chromatin organization, protein phosphorylation pathways. Based on Venn analysis, a total of 37 specific DEGs were identified in SDF-1+EX-4 group, which were mainly enriched in negative regulation of autophagy and cellular component disassembly signaling pathways. Short time-series expression miner (STEM) analysis grouped all expressed genes of PDLSCs into 49 clusters according to the dynamic expression patterns and 25 genes, including NRSN2, CHD9, TUBA1A, distributed in 10 gene clusters in SDF-1+EX-4 treated PDLSCs were significantly up-regulated compared with the SDF-1 and EX-4 alone groups. The gene set enrichment analysis indicated that SDF-1 could amplify the role of EX-4 in regulating varied signaling pathways, such as type II diabetes mellitus and insulin signaling pathways; while EX-4 could aggravate the effect of SDF-1 on PDLSCs biological roles via regulating primary immunodeficiency, tight junction signaling pathways. In summary, our study confirmed that SDF-1 and EX-4 combined application could enhance PDLSCs biological activity and promote PDLSCs osteogenic differentiation by regulating the metabolism, biosynthesis and immune-related signaling pathways.

Transcriptomics in Trichoderma reesei

Transcriptomics is a powerful technique to study gene expression. The main purpose of transcriptome studies in the filamentous fungus Trichoderma reesei is the analysis of differentially expressed genes as a transcriptional response of the genome to different environmental stimuli or physiological conditions such as sugar availability, nitrogen metabolism, pH response, and oxidative stress, among others. Here we describe the full protocol of RNA sequencing methodology from RNA isolation to data analysis in order to access the T. reesei transcriptome.

RNA sequencing analysis of the human retina and associated ocular tissues

The retina is a stratified layer of sensory neurons lining the posterior portion of the eye. In humans, fine detail and color vision are enabled by the macula, a central region of the retina dense in cone photoreceptors (PRs). Achromatic low light and peripheral vision are facilitated by rod PRs found with increasing density outside the macula in the peripheral retina. The outer retina is nourished by choroidal blood flow regulated by a single layer of intervening retinal pigment epithelial (RPE) cells. Existing human retinal transcriptome projects have been critical for studying aspects of retinal development and disease however, there are currently no publicly available data sets accurately describing the aging human central retina, peripheral retina, and supporting RPE/choroid. Here we used Illumina RNA sequencing (RNA-seq) analysis to characterize the mRNA transcriptome of rod and cone PR-enriched human retina as well as supporting macular RPE/choroid tissue. These data will be valuable to the vision research community for characterizing global changes in gene expression in clinically relevant ocular tissues.

Impacts of ciliary neurotrophic factor on the retinal transcriptome in a mouse model of photoreceptor degeneration

Ciliary neurotrophic factor (CNTF) has been tested in clinical trials for human retinal degeneration due to its potent neuroprotective effects in various animal models. To decipher CNTF-triggered molecular events in the degenerating retina, we performed high-throughput RNA sequencing analyses using the Rds/Prph2 (P216L) transgenic mouse as a preclinical model for retinitis pigmentosa. In the absence of CNTF treatment, transcriptome alterations were detected at the onset of rod degeneration compared with wild type mice, including reduction of key photoreceptor transcription factors Crx, Nrl, and rod phototransduction genes. Short-term CNTF treatments caused further declines of photoreceptor transcription factors accompanied by marked decreases of both rod- and cone-specific gene expression. In addition, CNTF triggered acute elevation of transcripts in the innate immune system and growth factor signaling. These immune responses were sustained after long-term CNTF exposures that also affected neuronal transmission and metabolism. Comparisons of transcriptomes also uncovered common pathways shared with other retinal degeneration models. Cross referencing bulk RNA-seq with single-cell RNA-seq data revealed the CNTF responsive cell types, including Müller glia, rod and cone photoreceptors, and bipolar cells. Together, these results demonstrate the influence of exogenous CNTF on the retinal transcriptome landscape and illuminate likely CNTF impacts in degenerating human retinas.

Systematic review of differential methylation in rare ophthalmic diseases

Rare ophthalmic diseases have a devastating impact on a patient's vision and consequently negatively affect their independence, ability to work and overall quality of life. Methylation is an important emerging biomarker of disease and may improve understanding of rare ophthalmic disorders. This systematic review sought to identify and evaluate literature on methylation and rare ophthalmic disease. MEDLINE, EMBASE, PubMed, Cochrane Database of Systematic Reviews and grey literature resources were searched for publications prior to 20 August 2019. Articles written in English which featured key terms such as 'methylation' and rare ophthalmic diseases were included. Titles, abstracts, keywords and full texts of publications were screened, as well as reference lists for reverse citations and Web of Science 'cited reference search' for forward citation searching. Study characteristics were extracted, and methodological rigour appraised using a standardised template. Fourteen articles were selected for full inclusion. Rare ophthalmic conditions include congenital fibrosis of extraocular muscles, retinitis pigmentosa, Fuchs endothelial corneal dystrophy, granular corneal dystrophy, choroideraemia, brittle cornea syndrome, retinopathy of prematurity, keratoconus and congenital cataracts. Outcomes include identification of methylation as contributor to disease and identification of potential novel therapeutic targets. The studies included were heterogeneous with no scope for meta-analysis following review; a narrative synthesis was undertaken. Differential methylation has been identified in a small number of rare ophthalmic diseases and few studies have been performed to date. Further multiomic research will improve understanding of rare eye diseases and hopefully lead to improved provision of diagnostic/prognostic biomarkers, and help identify novel therapeutic targets.

The Molecular Basis of Human Anophthalmia and Microphthalmia

Human eye development is coordinated through an extensive network of genetic signalling pathways. Disruption of key regulatory genes in the early stages of eye development can result in aborted eye formation, resulting in an absent eye (anophthalmia) or a small underdeveloped eye (microphthalmia) phenotype. Anophthalmia and microphthalmia (AM) are part of the same clinical spectrum and have high genetic heterogeneity, with >90 identified associated genes. By understanding the roles of these genes in development, including their temporal expression, the phenotypic variation associated with AM can be better understood, improving diagnosis and management. This review describes the genetic and structural basis of eye development, focusing on the function of key genes known to be associated with AM. In addition, we highlight some promising avenues of research involving multiomic approaches and disease modelling with induced pluripotent stem cell (iPSC) technology, which will aid in developing novel therapies.

Human papillomavirus 16E6/E7 activates autophagy via Atg9B and LAMP1 in cervical cancer cells

BACKGROUNDS

Although the role of high-risk human papillomavirus (HPV) E6 and E7 in cellular malignant transformation has been elucidated， the function of both genes in cellular homeostasis is still unknown. Autophagy functions in maintenance of cellular homeostasis play a key role in the initiation and development of cancer and infectious disease.

METHODS

Cervical cancer cell lines SiHa and CaSki were utilized in this study.

RESULTS

We found that HPV 16E6/E7 (16E6/E7) downregulation inhibited autophagy, and consequently suppressed cell proliferation and promoted early apoptosis. Transcriptome sequencing demonstrated that Atg9B and LAMP1 were downregulated in 16E6/E7 knockdown cells. Gene function experiments revealed that 16E6/E7 downregulation depressed Atg9B and LAMP1, and Atg9B and LAMP1 overexpression compensated, at least partially, autophagy blockage induced by 16E6/E7 knockdown. Immunoprecipitation assay showed that 16E7 interacted with Atg9B and dual-luciferase reporter system revealed that 16E6 most likely regulated -1750 to -2000 nt in Atg9B and -1800 to -2000 nt in LAMP1 promoter region.

CONCLUSIONS

Our findings verified that 16E6/E7 activated autophagy via accelerating autophagosome formation and degradation, and Atg9B and LAMP1 were involved in the process of 16E6/E7 modulating autophagy, suggesting that targeting autophagy may be a potential approach in cervical cancer therapeutics.

Digital gene expression analyses of mammary glands from meat ewes naturally infected with clinical mastitis

Clinical mastitis in sheep has gravely restrained production performance for a long time. Knowledge of mechanisms of its pathogenesis and resistance in meat sheep mammary gland with clinical mastitis are not yet understood, especially for clinical mastitis caused by natural infection. In this work, RNA-sequencing was firstly used to screen the differentially expressed genes (DEGs) in clinical mastitic mammary tissues (CMMTs) when compared with healthy mammary tissues (HMTs) from meat sheep flocks. We identified 420 DEGs including 316 upregulated and 104 downregulated genes in CMMTs. Gene ontology annotation revealed these DEGs were mainly engaged in immune response and inflammation response. Pathway enrichment showed they were primarily enriched in pathways relevant to inflammation, immune response and metabolism. Alternative splicing analysis showed most common differential splicing genes in CMMTs and HMTs were implicated in immune response. Immunostaining for three immune response-related proteins encoded by DEGs were mainly observed in mammary epithelium from both CMMTs and HMTs, and their positive signals were more intensive in CMMTs than those in HMTs. These findings provide experimental basis and reference for further researching the molecular genetic mechanisms, particularly immune defence mechanisms, of sheep mammary gland during clinical mastitis.

Protocol for a scoping review of multi-omic analysis for rare diseases

INTRODUCTION

The development of next generation sequencing technology has enabled cost-efficient, large scale, multiple 'omic' analysis, including epigenomic, genomic, metabolomic, phenomic, proteomic and transcriptomic research. These integrated approaches hold significant promise for rare disease research, with the potential to aid biomarker discovery, improve our understanding of disease pathogenesis and identify novel therapeutic targets. In this paper we outline a systematic approach for a scoping review designed to evaluate what primary research has been performed to date on multi-omics and rare disease.

METHODS AND ANALYSIS

This protocol was designed using the Joanna Briggs Institute methodology for scoping reviews. Databases to be searched will include: MEDLINE, EMBASE, PubMed, Web of Science, Scopus and Google Scholar for primary studies relevant to the key terms 'multi-omics' and 'rare disease', published prior to 30th December 2018. Grey literature databases GreyLit and OpenGrey will also be searched, as well as reverse citation screening of relevant articles and forward citation searching using Web of Science Cited Reference Search Tool. Data extraction will be performed using customised forms and a narrative synthesis of the results will be presented.

ETHICS AND DISSEMINATION

As a secondary analysis study with no primary data generated, this scoping review does not require ethical approval. We anticipate this review will highlight a gap in rare disease research and provide direction for novel research. The completed review will be submitted for publication in peer-reviewed journals and presented at relevant conferences discussing rare disease research and/or molecular strategies.

Identification of the miRNA-mRNA regulatory network of antler growth centers

Antler growth is a unique event compared to other growth and development processes in mammals. Antlers grow extremely fast during the rapid growth stage when growth rate peaks at 2 cm per day. Antler growth is driven by a specific endochondral ossification process in the growth center that is in the distal region of the antler tip. In this study, we used state-of-art RNA-seq technology to analyze the expression profiles of mRNAs and miRNAs during antler growth. Our results indicated that the expression levels of multiple genes involved in chondrogenesis and endochondral ossification, including Fn1, Sox9, Col2a1, Acan, Col9a1, Col11a1, Hapln1, Wwp2, Fgfr3, Comp, Sp7 and Ihh, were significantly increased at the rapid growth stage. Our results also indicated that there were multiple differentially expressed miRNAs interacting with differentially expressed genes with opposite expression patterns. Furthermore, some of the miRNAs, including miR-3072-5p, miR-1600, miR-34-5p, miR-6889-5p and miR-6729-5p, simultaneously interacted with and controlled multiple genes involved in the process of chondrogenesis and endochondral ossification. Therefore, we established a miRNA-mRNA regulatory network by identifying miRNAs and their target genes that were differentially expressed in the antler growth centers by comparing the rapid growth stage and the initial growth stage.

Aberrant Expressions of Co-stimulatory and Co-inhibitory Molecules in Autoimmune Diseases

Co-signaling molecules include co-stimulatory and co-inhibitory molecules and play important roles in modulating immune responses. The roles of co-signaling molecules in autoimmune diseases have not been clearly defined. We assessed the expressions of co-stimulatory and co-inhibitory molecules in autoimmune diseases through a bioinformatics-based study. By using datasets of whole-genome transcriptome, the expressions of 54 co-stimulatory or co-inhibitory genes in common autoimmune diseases were analyzed using Robust rank aggregation (RRA) method. Nineteen array datasets and 6 RNA-seq datasets were included in the RRA discovery study and RRA validation study, respectively. Significant genes were further validated in several autoimmune diseases including Graves' disease (GD). RRA discovery study suggested that CD160 was the most significant gene aberrantly expressed in autoimmune diseases (Adjusted P = 5.9E-12), followed by CD58 (Adjusted P = 5.7E-06) and CD244 (Adjusted P = 9.5E-05). RRA validation study also identified CD160 as the most significant gene aberrantly expressed in autoimmune diseases (Adjusted P = 5.9E-09). We further found that the aberrant expression of CD160 was statistically significant in multiple autoimmune diseases including GD (P < 0.05), and CD160 had a moderate role in diagnosing those autoimmune diseases. Flow cytometry confirmed that CD160 was differentially expressed on the surface of CD8⁺ T cells between GD patients and healthy controls (P = 0.002), which proved the aberrant expression of CD160 in GD at the protein level. This study suggests that CD160 is the most significant co-signaling gene aberrantly expressed in autoimmune diseases. Treatment strategy targeting CD160-related pathway may be promising for the therapy of autoimmune diseases.

A two-component signal transduction system contributes to the virulence of Riemerella anatipestifer

Similar to other studies of bacterial pathogens, current studies of the pathogenesis of Riemerella anatipestifer (RA) are focused mainly on in vitro culture conditions. To elucidate further the pathogenesis of RA in vivo, bacterial RNA was extracted from overnight tryptic soy broth cultures (in vitro) and from the blood of infected ducks (in vivo) for comparative RNA sequencing analysis. In total, 682 upregulated genes were identified in vivo. Among the upregulated genes, a signal transduction response regulator (ArsR) and a signal transduction histidine kinase (SthK) were predicted to be located on the same operon. A mutant was constructed by deletion of both of these genes. Duck infection tests showed that genes ArsR and SthK were related to the virulence of the pathogen in vivo. Differentially expressed genes identified by comparison of in vitro and in vivo conditions provided an insight into the physiological process of RA infection and provided an opportunity to identify additional virulence factors.

Functional networks of aging markers in the glomeruli of IgA nephropathy: a new therapeutic opportunity

IgA nephropathy(IgAN) is the most common primary glomerular disease in China. Primary infections always occur before IgAN. However, the pathology of IgAN is still unclear. Previously we found that LL37, a protein secreted by senescent cells, was specific for the progression of IgAN, and also played a role in the neutrophil function. So we hypothesized that the infiltration of neutrophils, inflammation factors, and aging markers, which were modulated by functional networks, induced the immune response and renal injury. RNA-Sequencing (RNA-seq) can be used to study the whole transcriptome and detect splicing variants that are expressed in a specific cell type or tissue. We separate glomerulus from the renal biopsy tissues. After RNA extraction, the sequences were analyzed with Illumina HiSeq 2000/2500. 381 genes with differential expression between the IgAN patients and the healthy controls were identified. Only PLAU, JUN, and FOS were related to DNA damage, telomere dysfunction-induced aging markers, neutrophil function and IgA nephropathy. The networks showed the possibility of these genes being connected. We conclude that DNA damage and telomere dysfunction could play important roles in IgA nephropathy. In addition, neutrophils are also important factors in this disease. The networks of these markers showed the mechanism pathways that are involved in the duration of the occurrence and progression of IgA nephropathy and might be a new therapeutic opportunity for disease treatment.

Correlation between RNA-Seq and microarrays results using TCGA data

RNA sequencing (RNA-Seq) and microarray are two of the most commonly used high-throughput technologies for transcriptome profiling; however, they both have their own inherent strengths and limitations. This research aims to analyze the correlation between microarrays and RNA-Seq detection of transcripts in the same tissue sample to explore the reproducibility between the techniques. Using data of RNA-Seq v2 and three different microarrays provided by The Cancer Genome Atlas, 11,120 genes of 111 lung squamous cell carcinoma samples were simultaneously detected by the four methods. Then we analyzed the Pearson correlation between microarrays and RNA-Seq. Finally, in the six comparison results, 9984 (89.8%) genes, irrespective of which two methods were used, simultaneously showed the existence of correlation, whereas only 83 (0.1%) genes proved to have no significant correlation in either comparison. In addition, the comparisons between 3266 (29.3%) genes showed high correlation (R≥0.8) in all six comparisons, only for 1643 (14.8%) genes correlation were not as high in either comparison. Meanwhile, transcripts with extreme high or low expression levels were more highly discrepant across the methods. In conclusion, we found that, for most transcripts, the results obtained by RNA-Seq and microarrays were highly reproducible.

Identification of a RNA-Seq Based 8-Long Non-Coding RNA Signature Predicting Survival in Esophageal Cancer

Background

Accumulating evidence suggests the involvement of long non-coding RNAs (lncRNAs) as oncogenic or tumor suppressive regulators in the development of various cancers. In the present study, we aimed to identify a lncRNA signature based on RNA sequencing (RNA-seq) data to predict survival in esophageal cancer.

Material/Methods

The RNA-seq lncRNA expression data and clinical information were downloaded from The Cancer Genome Atlas (TCGA) database. Differentially expressed lncRNAs were screened out between esophageal cancer and normal tissues. Univariate and multivariate Cox regression analysis were performed to establish a lncRNA-related prognostic model. Receiver operating characteristic (ROC) analysis was conducted to test the sensitivity and specificity of the model. GO (gene ontology) functional and KEGG pathway enrichment analyses were performed for mRNAs co-expressed with the lncRNAs to explore the potential functions of the prognostic lncRNAs.

Results

A total of 265 differentially expressed lncRNAs were identified between esophageal cancer and normal tissues. After univariate and multivariate Cox regression analysis, eight lncRNAs (GS1-600G8.5, LINC00365, CTD-2357A8.3, RP11-705O24.1, LINC01554, RP1-90J4.1, RP11-327J17.1, and LINC00176) were finally screened out to establish a predictive model by which patients could be classified into high-risk and low-risk groups with significantly different overall survival. Further analysis indicated independent prognostic capability of the 8-lncRNA signature from other clinicopathological factors. ROC curve analysis demonstrated good performance of the 8-lncRNA signature. Functional enrichment analysis showed that the prognostic lncRNAs were mainly associated with esophageal cancer related biological processes such as regulation of glucose metabolic process and amino acid and lipids metabolism.

Conclusions

Our study developed a novel candidate model providing additional and more powerful prognostic information beyond conventional clinicopathological factors for survival prediction of esophageal cancer patients. Moreover, it also brings us new insights into the molecular mechanisms underlying esophageal cancer.