The vast amount of RNA-seq data deposited in Gene Expression Omnibus (GEO) and Sequence Read Archive (SRA) is still a grossly underutilized resource for biomedical research. To remove technical roadblocks for reusing these data, we have developed a web-application GREIN (GEO RNA-seq Experiments Interactive Navigator) which provides user-friendly interfaces to manipulate and analyze GEO RNA-seq data. GREIN is powered by the back-end computational pipeline for uniform processing of RNA-seq data and the large number (>6,500) of already processed datasets. The front-end user interfaces provide a wealth of user-analytics options including sub-setting and downloading processed data, interactive visualization, statistical power analyses, construction of differential gene expression signatures and their comprehensive functional characterization, and connectivity analysis with LINCS L1000 data. The combination of the massive amount of back-end data and front-end analytics options driven by user-friendly interfaces makes GREIN a unique open-source resource for re-using GEO RNA-seq data. GREIN is accessible at:https://shiny.ilincs.org/grein, the source code at:https://github.com/uc-bd2k/grein, and the Docker container at:https://hub.docker.com/r/ucbd2k/grein.

The Library of Integrated Network-based Cellular Signatures (LINCS) program is a national consortium funded by the NIH to generate a diverse and extensive reference library of cell-based perturbation-response signatures, along with novel data analytics tools to improve our understanding of human diseases at the systems level. In contrast to other large-scale data generation efforts, LINCS Data and Signature Generation Centers (DSGCs) employ a wide range of assay technologies cataloging diverse cellular responses. Integration of, and unified access to LINCS data has therefore been particularly challenging. The Big Data to Knowledge (BD2K) LINCS Data Coordination and Integration Center (DCIC) has developed data standards specifications, data processing pipelines, and a suite of end-user software tools to integrate and annotate LINCS-generated data, to make LINCS signatures searchable and usable for different types of users. Here, we describe the LINCS Data Portal (LDP) (http://lincsportal.ccs.miami.edu/), a unified web interface to access datasets generated by the LINCS DSGCs, and its underlying database, LINCS Data Registry (LDR). LINCS data served on the LDP contains extensive metadata and curated annotations. We highlight the features of the LDP user interface that is designed to enable search, browsing, exploration, download and analysis of LINCS data and related curated content.

Rhinovirus (HRV) is associated with the large majority of virus-induced asthma exacerbations in children and young adults, but the mechanisms remain poorly defined. Asthmatics and non-asthmatic controls were inoculated with HRV-A16, and nasal epithelial samples were obtained 7 days before, 36 hours after, and 7 days after viral inoculation. RNA was extracted and subjected to RNA-seq analysis. At baseline, 57 genes were differentially expressed between asthmatics and controls, and the asthmatics had decreased expression of viral replication inhibitors and increased expression of genes involved in inflammation. At 36 hours (before the emergence of peak symptoms), 1329 genes were significantly altered from baseline in the asthmatics compared to 62 genes in the controls. At this time point, asthmatics lacked an increase in IL-10 signaling observed in the controls. At 7 days following HRV inoculation, 222 genes were significantly dysregulated in the asthmatics, whereas only 4 genes were dysregulated among controls. At this time point, the controls but not asthmatics demonstrated upregulation of SPINK5. As judged by the magnitude and persistence of dysregulated genes, asthmatics have a substantially different host response to HRV-A16 infection compared with non-asthmatic controls. Gene expression differences illuminate biologically plausible mechanisms that contribute to a better understanding of the pathogenesis of HRV-induced asthma exacerbations.

One of the main objectives of many biological studies is to explore differential gene expression profiles between samples. Genes are referred to as differentially expressed (DE) if the read counts change across treatments or conditions systematically. Poisson and negative binomial (NB) regressions are widely used methods for non-over-dispersed (NOD) and over-dispersed (OD) count data respectively. However, in the presence of excessive number of zeros, these methods need adjustments. In this paper, we consider a zero-inflated Poisson mixed effects model (ZIPMM) and zero-inflated negative binomial mixed effects model (ZINBMM) to address excessive zero counts in the NOD and OD RNA-seq data respectively in the presence of random effects. We apply these methods to both simulated and real RNA-seq datasets. The ZIPMM and ZINBMM perform better on both simulated and real datasets.

Single-cell RNA sequencing (scRNA-seq) provides a comprehensive understanding of cellular complexity; however, its requirement for fresh or frozen samples limits its flexibility. To overcome this limitation to effectively leverage clinical samples, Chromium Fixed RNA Profiling on formalin-fixed paraffin-embedded (FFPE) tissue blocks (scFFPE-seq) was developed to perform single-nucleus RNA sequencing from nuclei isolated from FFPE. In this study, we utilized fresh tissue samples from colon, ileum, and skin to assess the viability of scFFPE-seq compared to these fresh samples. We were able to recover unique cell types from challenging FFPE tissues and validated scFFPE-seq findings through Hematoxylin and Eosin (H&E) images. The results demonstrated that scFFPE-seq effectively captured the single-cell transcriptome in FFPE tissues, obtaining comparable cell abundance, cell type annotation, and pathway characterization to those in fresh tissues. Overall, the study presents strong evidence of the potential of scFFPE-seq to enhance scientific knowledge by enabling the generation of high-quality, sensitive single-nucleus RNA-seq data from preserved tissue samples. This technique unlocks the vast archives of FFPE samples for extensive retrospective genomic studies.