A New ℓ0-Regularized Log-Linear Poisson Graphical Model with Applications to RNA Sequencing Data

2021 ◽  
Author(s):  
Caesar Z. Li ◽  
Eric S. Kawaguchi ◽  
Gang Li
Keyword(s):  
Rna Sequencing ◽  
Graphical Model ◽  
Sequencing Data ◽  
Log Linear

scholarly journals A relative comparison between Hidden Markov- and Log-Linear-based models for differential expression analysis in a real time course RNA sequencing data

2018 ◽  
Author(s):  
Fatemeh Gholizadeh ◽  
Zahra Salehi ◽  
Ali Mohammad banaei-Moghaddam ◽  
Abbas Rahimi Foroushani ◽  
Kaveh kavousi
Keyword(s):  
Real Time ◽  
Differential Expression ◽  
Rna Sequencing ◽  
Time Course ◽  
Hidden Markov ◽  
Differential Expression Analysis ◽  
Rna Seq ◽  
Sequencing Data ◽  
Normalization Methods ◽  
Log Linear

AbstractWith the advent of the Next Generation Sequencing technologies, RNA-seq has become known as an optimal approach for studying gene expression profiling. Particularly, time course RNA-seq differential expression analysis has been used in many studies to identify candidate genes. However, applying a statistical method to efficiently identify differentially expressed genes (DEGs) in time course studies is challenging due to inherent characteristics of such data including correlation and dependencies over time. Here we aim to relatively compare EBSeq-HMM, a Hidden Markov-based model, with multiDE, a Log-Linear-based model, in a real time course RNA sequencing data. In order to conduct the comparison, common DEGs detected by edgeR, DESeq2 and Voom (referred to as Benchmark DEGs) were utilized as a measure. Each of the two models were compared using different normalization methods. The findings revealed that multiDE identified more Benchmark DEGs and showed a higher agreement with them than EBSeq-HMM. Furthermore, multiDE and EBSeq-HMM displayed their best performance using TMM and Upper-Quartile normalization methods, respectively.

scholarly journals Graphic Model-based Gene Regulatory Network Reconstruction using RNA Sequencing Count Data

2019 ◽  
Vol 8 ◽  
pp. 1078-1085
Author(s):  
Liliana Lopez-Kleine ◽  
Cristian Andres Gonzalez-Prieto
Keyword(s):  
Rna Sequencing ◽  
Negative Binomial Distribution ◽  
Count Data ◽  
Binomial Distribution ◽  
Regulatory Networks ◽  
Graphical Model ◽  
Negative Binomial ◽  
Sequencing Data ◽  
Model Based ◽  
Gene Regulatory

Interactions between genes, such as regulations are best represented by gene regulatory networks (GRN). These are often constructed based on gene expression data. Few methods for the construction of GRN exist for RNA sequencing count data. One of the most used methods for microarray data is based on graphical Gaussian networks. Considering that count data have different distributions, a method assuming RNA sequencing counts distribute Poisson has been proposed recently. Nevertheless, it has been argued that the most likely distribution of RNA sequencing counts is not Poisson due to overdispersion. Therefore, the negative binomial distribution is much more likely. For this distribution, no model-based method for the construction of GRN has been proposed until now. Here, we present a graphical, model-based method for the construction of GRN assuming a negative binomial distribution of the RNA sequencing count data. The R code is available under request. We used the method proposed both on simulated RNA sequencing count data and on real data. The graph is showed, and its descriptive measurements were assessed. They were found some interesting biological conclusions. We confirm that using negative binomial distribution for fitting the model is suitable because RNA sequencing data present overdispersion.

scholarly journals HEDGEHOG/GLI Modulates the PRR11-SKA2 Bidirectional Transcription Unit in Lung Squamous Cell Carcinomas

Genes ◽  
2021 ◽  
Vol 12 (1) ◽  
pp. 120
Author(s):  
Yiyun Sun ◽  
Dandan Xu ◽  
Chundong Zhang ◽  
Yitao Wang ◽  
Lian Zhang ◽  
...  
Keyword(s):  
Rna Sequencing ◽  
Squamous Cell ◽  
Gene Pair ◽  
Gene List ◽  
Ectopic Expression ◽  
Lung Squamous Cell Carcinoma ◽  
Gene Set Enrichment Analysis ◽  
The Cancer Genome Atlas ◽  
Sequencing Data ◽  
Gene Set

We previously demonstrated that proline-rich protein 11 (PRR11) and spindle and kinetochore associated 2 (SKA2) constituted a head-to-head gene pair driven by a prototypical bidirectional promoter. This gene pair synergistically promoted the development of non-small cell lung cancer. However, the signaling pathways leading to the ectopic expression of this gene pair remains obscure. In the present study, we first analyzed the lung squamous cell carcinoma (LSCC) relevant RNA sequencing data from The Cancer Genome Atlas (TCGA) database using the correlation analysis of gene expression and gene set enrichment analysis (GSEA), which revealed that the PRR11-SKA2 correlated gene list highly resembled the Hedgehog (Hh) pathway activation-related gene set. Subsequently, GLI1/2 inhibitor GANT-61 or GLI1/2-siRNA inhibited the Hh pathway of LSCC cells, concomitantly decreasing the expression levels of PRR11 and SKA2. Furthermore, the mRNA expression profile of LSCC cells treated with GANT-61 was detected using RNA sequencing, displaying 397 differentially expressed genes (203 upregulated genes and 194 downregulated genes). Out of them, one gene set, including BIRC5, NCAPG, CCNB2, and BUB1, was involved in cell division and interacted with both PRR11 and SKA2. These genes were verified as the downregulated genes via RT-PCR and their high expression significantly correlated with the shorter overall survival of LSCC patients. Taken together, our results indicate that GLI1/2 mediates the expression of the PRR11-SKA2-centric gene set that serves as an unfavorable prognostic indicator for LSCC patients, potentializing new combinatorial diagnostic and therapeutic strategies in LSCC.

RNA Sequencing Data from Human Intracranial Aneurysm Tissue Reveals a Complex Inflammatory Environment Associated with Rupture

2021 ◽  
Author(s):  
Vincent M. Tutino ◽  
Haley R. Zebraski ◽  
Hamidreza Rajabzadeh-Oghaz ◽  
Lee Chaves ◽  
Adam A. Dmytriw ◽  
...  
Keyword(s):  
Intracranial Aneurysm ◽  
Rna Sequencing ◽  
Sequencing Data

scholarly journals In-depth transcriptomic analysis of human retina reveals molecular mechanisms underlying diabetic retinopathy

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Kolja Becker ◽  
Holger Klein ◽  
Eric Simon ◽  
Coralie Viollet ◽  
Christian Haslinger ◽  
...  
Keyword(s):  
Diabetic Retinopathy ◽  
Rna Sequencing ◽  
Molecular Mechanisms ◽  
Vision Loss ◽  
Ganglion Cells ◽  
Expression Profiles ◽  
Cell Types ◽  
Sequencing Data ◽  
Disease Stages ◽  
Post Transcriptional Regulation

AbstractDiabetic Retinopathy (DR) is among the major global causes for vision loss. With the rise in diabetes prevalence, an increase in DR incidence is expected. Current understanding of both the molecular etiology and pathways involved in the initiation and progression of DR is limited. Via RNA-Sequencing, we analyzed mRNA and miRNA expression profiles of 80 human post-mortem retinal samples from 43 patients diagnosed with various stages of DR. We found differentially expressed transcripts to be predominantly associated with late stage DR and pathways such as hippo and gap junction signaling. A multivariate regression model identified transcripts with progressive changes throughout disease stages, which in turn displayed significant overlap with sphingolipid and cGMP–PKG signaling. Combined analysis of miRNA and mRNA expression further uncovered disease-relevant miRNA/mRNA associations as potential mechanisms of post-transcriptional regulation. Finally, integrating human retinal single cell RNA-Sequencing data revealed a continuous loss of retinal ganglion cells, and Müller cell mediated changes in histidine and β-alanine signaling. While previously considered primarily a vascular disease, attention in DR has shifted to additional mechanisms and cell-types. Our findings offer an unprecedented and unbiased insight into molecular pathways and cell-specific changes in the development of DR, and provide potential avenues for future therapeutic intervention.

scholarly journals COVID-19 Severity Potentially Modulated by Cardiovascular-Disease-Associated Immune Dysregulation

Viruses ◽  
2021 ◽  
Vol 13 (6) ◽  
pp. 1018
Author(s):  
Abby C. Lee ◽  
Grant Castaneda ◽  
Wei Tse Li ◽  
Chengyu Chen ◽  
Neil Shende ◽  
...  
Keyword(s):  
Rna Sequencing ◽  
Mononuclear Cells ◽  
Human Peripheral Blood ◽  
Expression Profiles ◽  
Immune Dysregulation ◽  
Left Ventricular ◽  
Recurrent Event ◽  
Healthy Controls ◽  
Sequencing Data ◽  
Peripheral Blood Mononuclear

Patients with underlying cardiovascular conditions are particularly vulnerable to severe COVID-19. In this project, we aimed to characterize similarities in dysregulated immune pathways between COVID-19 patients and patients with cardiomyopathy, venous thromboembolism (VTE), or coronary artery disease (CAD). We hypothesized that these similarly dysregulated pathways may be critical to how cardiovascular diseases (CVDs) exacerbate COVID-19. To evaluate immune dysregulation in different diseases, we used four separate datasets, including RNA-sequencing data from human left ventricular cardiac muscle samples of patients with dilated or ischemic cardiomyopathy and healthy controls; RNA-sequencing data of whole blood samples from patients with single or recurrent event VTE and healthy controls; RNA-sequencing data of human peripheral blood mononuclear cells (PBMCs) from patients with and without obstructive CAD; and RNA-sequencing data of platelets from COVID-19 subjects and healthy controls. We found similar immune dysregulation profiles between patients with CVDs and COVID-19 patients. Interestingly, cardiomyopathy patients display the most similar immune landscape to COVID-19 patients. Additionally, COVID-19 patients experience greater upregulation of cytokine- and inflammasome-related genes than patients with CVDs. In all, patients with CVDs have a significant overlap of cytokine- and inflammasome-related gene expression profiles with that of COVID-19 patients, possibly explaining their greater vulnerability to severe COVID-19.

scholarly journals IMMU-27. SINGLE CELL RNA-SEQUENCING IDENTIFIES NOVEL BONE MARROW DERIVED MYELOID CELLS IN GLIOBLASTOMA ASSOCIATED WITH TUMOR AGGRESSION

2020 ◽  
Vol 22 (Supplement_2) ◽  
pp. ii110-ii110
Author(s):  
Christina Jackson ◽  
Christopher Cherry ◽  
Sadhana Bom ◽  
Hao Zhang ◽  
John Choi ◽  
...  
Keyword(s):  
Bone Marrow ◽  
Single Cell ◽  
Tumor Cells ◽  
Rna Sequencing ◽  
Metabolic Pathways ◽  
Myeloid Cells ◽  
Tumor Grade ◽  
Sequencing Data ◽  
Single Cell Rna Sequencing ◽  
Two Populations

Abstract BACKGROUND Glioma associated myeloid cells (GAMs) can be induced to adopt an immunosuppressive phenotype that can lead to inhibition of anti-tumor responses in glioblastoma (GBM). Understanding the composition and phenotypes of GAMs is essential to modulating the myeloid compartment as a therapeutic adjunct to improve anti-tumor immune response. METHODS We performed single-cell RNA-sequencing (sc-RNAseq) of 435,400 myeloid and tumor cells to identify transcriptomic and phenotypic differences in GAMs across glioma grades. We further correlated the heterogeneity of the GAM landscape with tumor cell transcriptomics to investigate interactions between GAMs and tumor cells. RESULTS sc-RNAseq revealed a diverse landscape of myeloid-lineage cells in gliomas with an increase in preponderance of bone marrow derived myeloid cells (BMDMs) with increasing tumor grade. We identified two populations of BMDMs unique to GBMs; Mac-1and Mac-2. Mac-1 demonstrates upregulation of immature myeloid gene signature and altered metabolic pathways. Mac-2 is characterized by expression of scavenger receptor MARCO. Pseudotime and RNA velocity analysis revealed the ability of Mac-1 to transition and differentiate to Mac-2 and other GAM subtypes. We further found that the presence of these two populations of BMDMs are associated with the presence of tumor cells with stem cell and mesenchymal features. Bulk RNA-sequencing data demonstrates that gene signatures of these populations are associated with worse survival in GBM. CONCLUSION We used sc-RNAseq to identify a novel population of immature BMDMs that is associated with higher glioma grades. This population exhibited altered metabolic pathways and stem-like potentials to differentiate into other GAM populations including GAMs with upregulation of immunosuppressive pathways. Our results elucidate unique interactions between BMDMs and GBM tumor cells that potentially drives GBM progression and the more aggressive mesenchymal subtype. Our discovery of these novel BMDMs have implications in new therapeutic targets in improving the efficacy of immune-based therapies in GBM.

scholarly journals Differential expression of TgMIC1 in isolates of Chinese 1 Toxoplasma with different virulence

2021 ◽  
Vol 14 (1) ◽  
Author(s):  
Yang Wang ◽  
Chengjian Han ◽  
Rongsheng zhou ◽  
Jinjin Zhu ◽  
Famin Zhang ◽  
...  
Keyword(s):  
Differential Expression ◽  
Rna Sequencing ◽  
Polyclonal Antibody ◽  
Protein Level ◽  
Polyclonal Antibodies ◽  
Mrna Levels ◽  
Quantitative Real Time Pcr ◽  
Sequencing Data ◽  
Predominant Genotype ◽  
High Virulence

Abstract Background The predominant genotype of Toxoplasma in China is the Chinese 1 (ToxoDB#9) lineage. TgCtwh3 and TgCtwh6 are two representative strains of Chinese 1, exhibiting high and low virulence to mice, respectively. Little is known regarding the virulence mechanism of this non-classical genotype. Our previous RNA sequencing data revealed differential mRNA levels of TgMIC1 in TgCtwh3 and TgCtwh6. We aim to further confirm the differential expression of TgMIC1 and its significance in this atypical genotype. Methods Quantitative real-time PCR was used to verify the RNA sequencing data; then, polyclonal antibodies against TgMIC1 were prepared and identified. Moreover, the invasion and proliferation of the parasite in HFF cells were observed after treatment with TgMIC1 polyclonal antibody or not. Results The data showed that the protein level of TgMIC1 was significantly higher in high-virulence strain TgCtwh3 than in low-virulence strain TgCtwh6 and that the invasion and proliferation of TgCtwh3 were inhibited by TgMIC1 polyclonal antibody. Conclusion Differential expression of TgMIC1 in TgCtwh3 and TgCtwh6 may explain, at least partly, the virulence mechanism of this atypical genotype.

Sign in / Sign up

Export Citation Format

Share Document