Perturbation Of Gene Expression Research Articles

Appraisal of Cinnamaldehyde Analogs as Dual-Acting Antibiofilm and Anthelmintic Agents.

Cinnamaldehyde has a broad range of biological activities, which include antibiofilm and anthelmintic activities. The ever-growing problem of drug resistance and limited treatment options have created an urgent demand for natural molecules with antibiofilm and anthelmintic properties. Hence, we hypothesized that molecules with a scaffold structurally similar to that of cinnamaldehyde might act as dual inhibitors against fungal biofilms and helminths. In this regard, eleven cinnamaldehyde analogs were tested to determine their effects on fungal Candida albicans biofilm and nematode Caenorhabditis elegans. α-Methyl and trans-4-methyl cinnamaldehydes efficiently inhibited C. albicans biofilm formation (>90% inhibition at 50 μg/mL) with minimum inhibitory concentrations (MICs) of ≥ 200 μg/mL and 4-bromo and 4-chloro cinnamaldehydes exhibited anthelmintic property at 20 μg/mL against C. elegans. α-Methyl and trans-4-methyl cinnamaldehydes inhibited hyphal growth and cell aggregation. Scanning electron microscopy was employed to determine the surface architecture of C. albicans biofilm and cuticle of C. elegans, and confocal laser scanning microscopy was used to determine biofilm characteristics. The perturbation in gene expression of C. albicans was investigated using qRT-PCR analysis and α-methyl and trans-4-methyl cinnamaldehydes exhibited down-regulation of ECE1, IFD6, RBT5, UCF1, and UME6 and up-regulation of CHT4 and YWP1. Additionally, molecular interaction of these two molecules with UCF1 and YWP1 were revealed by molecular docking simulation. Our observations collectively suggest α-methyl and trans-4-methyl cinnamaldehydes are potent biofilm inhibitors and that 4-bromo and 4-chloro cinnamaldehydes are anthelmintic agents. Efforts are required to determine the range of potential therapeutic applications of cinnamaldehyde analogs.

Assembly of Dynamic Gated and Cascaded Transient DNAzyme Networks.

The dynamic transient formation and depletion of G-quadruplexes regulate gene replication and transcription. This process was found to be related to various diseases such as cancer and premature aging. We report on the engineering of nucleic acid modules revealing dynamic, transient assembly and disassembly of G-quadruplex structures and G-quadruplex-based DNAzymes, gated transient processes, and cascaded dynamic transient reactions that involve G-quadruplex and DNAzyme structures. The dynamic transient processes are driven by functional DNA reaction modules activated by a fuel strand and guided toward dissipative operation by a nicking enzyme (Nt.BbvCI). The dynamic networks were further characterized by computational simulation of the experiments using kinetic models, allowing us to predict the dynamic performance of the networks under different auxiliary conditions applied to the systems. The systems reported herein could provide functional DNA machineries for the spatiotemporal control of G-quadruplex structures perturbing gene expression and thus provide a therapeutic means for related emergent diseases.

The flocculant Saccharomyces cerevisiae strain gains robustness via alteration of the cell wall hydrophobicity

When lignocellulosic biomass is utilized as a fermentative substrate to produce biochemicals, the existence of a yeast strain resistant to inhibitory chemical compounds (ICCs) released from the biomass becomes critical. To achieve the purpose, in this study, Saccharomyces yeast strains from a NBRC yeast culture collection were used for exploration and evaluated in two different media containing ICCs that mimic one another but resemble the hydrolysate of real biomass. Among them, S. cerevisiae F118 strain shows robustness upon the fermentation with unique flocculation trait that was strongly responsive to ICC stress. When this strain was cultured in the presence of ICCs, its cell wall hydrophobicity increased dramatically, and reduced significantly when the ICCs were depleted, demonstrating that cell-surface hydrophobicity can also act as an adaptive response to the ICCs. Cells from the strain with the highest cell-wall hydrophobicity displayed progressively stronger flocculation, indicating that the F118 strain is having unique robustness under ICC stress. Gene expression perturbation analysis revealed that mot3 gene encoding regulatory Mot3p from the F118 strain was expressed in response to the concentration of ICCs. This gene was found to control expression of ygp1 gene that encoding Ygp1p, one of cell wall proteins. Deep sequencing analysis revealed that the Mot3p of the F118 strain features a unique insertion and deletion of nucleotides that encode glutamine or asparagine residues, particularly in N-terminal domain, as determined by comparison to the Mot3p sequence from the S288c strain, which was employed as a control strain. Furthermore, the cell wall hydrophobicity of the S288c strain was greatly enhanced and became ICC-responsive after gene swapping with the mot3 gene from the F118 strain. The gene-swapped S288c strain fermented 6-fold faster than the wild-type strain, producing 14.5 g/L of ethanol from 30 g/L of glucose consumed within 24 h in a medium containing the ICCs. These such modifications to Mot3p in unique locations in its sequence have a potential to change the expression of a gene involved in cell wall hydrophobicity and boosted the flocculation response to ICC stress, allowing for the acquisition of extraordinary robustness.

Seminal fluid gene expression and reproductive fitness in Drosophila melanogaster

BackgroundThe rapid evolution of seminal fluid proteins (SFPs) has been suggested to be driven by adaptations to postcopulatory sexual selection (e.g. sperm competition). However, we have recently shown that most SFPs evolve rapidly under relaxed selective pressures. Given the role of SFPs in competition for fertilization phenotypes, like the ability to transfer and store sperm and the modulation of female receptivity and ovulation, the prevalence of selectively relaxed SFPs appears as a conundrum. One possible explanation is that selection on SFPs might be relaxed in terms of protein amino acid content, but adjustments of expression are essential for post-mating function. Interestingly, there is a general lack of systematic implementation of gene expression perturbation assays to monitor their effect on phenotypes related to sperm competition.ResultsWe successfully manipulated the expression of 16 SFP encoding genes using tissue-specific knockdowns (KDs) and determined the effect of these genes’ perturbation on three important post-mating phenotypes: female refractoriness to remating, defensive (P1), and offensive (P2) sperm competitive abilities in Drosophila melanogaster. Our analyses show that KDs of tested SFP genes do not affect female refractoriness to remating and P2, however, most gene KDs significantly decreased P1. Moreover, KDs of SFP genes that are selectively constrained in terms of protein-coding sequence evolution have lower P1 than KDs of genes evolving under relaxed selection.ConclusionsOur results suggest a more predominant role, than previously acknowledged, of variation in gene expression than coding sequence changes on sperm competitive ability in D. melanogaster.

Potential toxicity of inorganic ions in particulate matter: Ion permeation in lung and disruption of cell metabolism

Exposure to ambient particulate matter (PM) is associated with adverse health effects. Yet, due to the complexity of its chemical composition, the molecular effects of PM exposure and the mechanism of PM-mediated toxicity remain largely unknown. Here, we show that water-soluble inorganics such as nitrate and sulfate ions, rather than PM itself, rapidly penetrate the lung surfactant barrier to the alveolar region and perturb gene expression in the lungs. Through high-throughput sequencing of lung adenocarcinoma cells, we find that exposure to nitrate and sulfate ions activates the cholesterol biosynthetic metabolism and induces the expression of genes related to tumorigenesis. Transcriptome analysis of mouse lungs exposed to nitrate/sulfate aerosols reveals interferon gamma-associated immune response. Interestingly, we find that exposure to a nitrate/sulfate mixture leads to a unique gene expression pattern that is not observed when nitrate or sulfate is treated alone. Our work suggests that the water-soluble ions are a potential source of PM-mediated toxicity and provides a roadmap to unveil the molecular mechanism of health hazards from PM exposure.

Directed mutational scanning reveals a balance between acidic and hydrophobic residues in strong human activation domains.

SUMMARYAcidic activation domains are intrinsically disordered regions of the transcription factors that bind coactivators. The intrinsic disorder and low evolutionary conservation of activation domains have made it difficult to identify the sequence features that control activity. To address this problem, we designed thousands of variants in seven acidic activation domains and measured their activities with a high-throughput assay in human cell culture. We found that strong activation domain activity requires a balance between the number of acidic residues and aromatic and leucine residues. These findings motivated a predictor of acidic activation domains that scans the human proteome for clusters of aromatic and leucine residues embedded in regions of high acidity. This predictor identifies known activation domains and accurately predicts previously unidentified ones. Our results support a flexible acidic exposure model of activation domains in which the acidic residues solubilize hydrophobic motifs so that they can interact with coactivators. A record of this paper’s transparent peer review process is included in the supplemental information.

MetaLINCS: an R package for meta-level analysis of LINCS L1000 drug signatures using stratified connectivity mapping.

Accessing the collection of perturbed gene expression profiles, such as the LINCS L1000 connectivity map, is usually performed at the individual dataset level, followed by a summary performed by counting individual hits for each perturbagen. With the metaLINCS R package, we present an alternative approach that combines rank correlation and gene set enrichment analysis to identify meta-level enrichment at the perturbagen level and, in the case of drugs, at the mechanism of action level. This significantly simplifies the interpretation and highlights overarching themes in the data. We demonstrate the functionality of the package and compare its performance against those of three currently used approaches. metaLINCS is released under GPL3 license. Source code and documentation are freely available on GitHub (https://github.com/bigomics/metaLINCS). Supplementary data are available at Bioinformatics Advances online.

Topologically switchable and gated transcription machinery.

Topological barriers control in nature the transcription machinery, thereby perturbing gene expression. Here we introduce synthetically designed DNA templates that include built-in topological barriers for switchable, triggered-controlled transcription of RNA aptamers. This is exemplified with the design of transcription templates that include reversible and switchable topological barriers consisting of a Sr2+-ion-stabilized G-quadruplex and its separation by kryptofix [2.2.2], KP, for the switchable transcription of the malachite green (MG) RNA aptamer, the T-A·T triplex barrier being separated by a fuel-strand for the cyclic triggered transcription of the 3,5-difluoro-4-hydroxybenzylidene imidazolinone (DFHBI)-binding aptamer, and the use of a photoactivated cis/trans azobenzene-modified nucleic acid barrier for the switchable “ON”/“OFF” transcription of the MG RNA aptamer. By applying a mixture of topologically triggered templates consisting of the photoresponsive barrier and the T-A·T triplex barrier, the gated transcription of the MG aptamer or the DFHBI-binding aptamer is demonstrated. In addition, a Sr2+-ion/KP topologically triggered DNA tetrahedra promoter-transcription scaffold, for the replication of the MG RNA aptamer, and T7 RNA polymerase are integrated into DNA-based carboxymethyl cellulose hydrogel microcapsules acting as cell-like assemblies. The switchable, reversible transcription of the MG RNA aptamer in a cell-like containment is introduced.

Measuring mRNA Decay with Roadblock-qPCR.

The control of mRNA stability is fundamental to gene regulation, and a deeper understanding of this post-transcriptional regulatory step can provide key insights into gene function. Measuring mRNA half-lives directly, however, is challenging. The most common strategies for evaluating mRNA stability and decay involve blocking general transcription and then measuring the decline in mRNA levels over time. The downside of these approaches, however, is that they severely impact cell function and viability, indirectly perturbing gene expression. Here, we describe Roadblock-qPCR, a simple method for measuring mRNA decay kinetics in living cells that is both economical and quick. Cells are first incubated with the nucleoside analog 4-thiouridine (4sU), which is readily incorporated into nascent mRNAs during transcription. RNA is then extracted and treated with N-ethylmaleimide (NEM), a sulfhydryl alkylating agent that selectively modifies 4sU, before proceeding to cDNA synthesis. Because the NEM-modified 4sU creates a chemical "roadblock" that interferes with reverse transcription, this treatment ultimately results in the depletion of the nascent 4sU-containing transcripts from the cDNA pool. As such, the decay rate of the non-4sU-labeled pre-existing mRNAs can be monitored by quantitative PCR (qPCR). In combination with spike-in standards, this approach can be used to efficiently and accurately measure the half-lives of endogenous mRNAs with a wide range of stabilities, while avoiding the artifacts of transcription shutoff strategies. © 2022 Wiley Periodicals LLC. Basic Protocol: Roadblock-qPCR Support Protocol: Synthesis of spike-in mRNA.

Transcriptomic signatures induced by the Ebola virus vaccine rVSVΔG-ZEBOV-GP in adult cohorts in Europe, Africa, and North America: a molecular biomarker study

SummaryBackgroundA recombinant vesicular stomatitis virus vector expressing the Zaire Ebola virus glycoprotein (rVSVΔG-ZEBOV-GP) vaccine has been reported as safe, immunogenic, and highly protective in a ring vaccination trial. We aimed to identify transcriptomic immune response biomarker signatures induced by vaccination and associated signatures with its immunogenicity and reactogenicity to better understand the potential mechanisms of action of the vaccine.Methods354 healthy adult volunteers were vaccinated in randomised, double-blind, placebo-controlled trials in Europe (Geneva, Switzerland [November, 2014, to January, 2015]) and North America (USA [Dec 5, 2014, to June 23, 2015]), and dose-escalation trials in Africa (Lambaréné, Gabon [November, 2014, to January, 2015], and Kilifi, Kenya [December, 2014, to January, 2015]) using different doses of the recombinant vesicular stomatitis virus vector expressing the Zaire Ebola virus glycoprotein (rVSVΔG-ZEBOV-GP; 3 × 105 to 1 × 108 plaque-forming units [pfu]). Longitudinal transcriptomic responses (days 0, 1, 2, 3, 7, 14, and 28) were measured in whole blood using a targeted gene expression profiling platform (dual-colour reverse-transcriptase multiplex ligation-dependent probe amplification) focusing on 144 immune-related genes. The effect of time and dose on transcriptomic response was also assessed. Logistic regression with lasso regularisation was applied to identify host signatures with optimal discriminatory capability of vaccination at day 1 or day 7 versus baseline, whereas random-effects models and recursive feature elimination combined with regularised logistic regression were used to associate signatures with immunogenicity and reactogenicity.FindingsOur results indicated that perturbation of gene expression peaked on day 1 and returned to baseline levels between day 7 and day 28. The magnitude of the response was dose-dependent, with vaccinees receiving a high dose (≥9 × 106 pfu) of rVSVΔG-ZEBOV-GP exhibiting the largest amplitude. The most differentially expressed genes that were significantly upregulated following vaccination consisted of type I and II interferon-related genes and myeloid cell-associated markers, whereas T cell, natural killer cell, and cytotoxicity-associated genes were downregulated. A gene signature associated with immunogenicity (common to all four cohorts) was identified correlating gene expression profiles with ZEBOV-GP antibody titres and a gene signatures associated with reactogenicity (Geneva cohort) was identified correlating gene expression profiles with an adverse event (ie, arthritis).InterpretationCollectively, our results identify and cross-validate immune-related transcriptomic signatures induced by rVSVΔG-ZEBOV-GP vaccination in four cohorts of adult participants from different genetic and geographical backgrounds. These signatures will aid in the rational development, testing, and evaluation of novel vaccines and will allow evaluation of the effect of host factors such as age, co-infection, and comorbidity on responses to vaccines.FundingInnovative Medicines Initiative 2 Joint Undertaking.

Elucidating cellular contributions to tau-mediated neurodegeneration using drosophila and single-cell transcriptomics.

Alzheimer's Dementia (AD) is neuropathologically defined by aggregations of amyloid-beta plaques and tau neurofibrillary tangles in the brain, with tau pathology being most strongly correlated with cognitive decline. Systems biology approaches profiling post-mortem human AD brain tissue have highlighted involvement of key processes such as innate immunity and translational regulation. However, interpretation of cross-sectional post-mortem gene expression data is confounded by co-morbidities and mixed brain pathologies at autopsy. Moreover, characterizing longitudinal, age-related changes in Alzheimer's dementia (AD) is of particular importance because clinical manifestation of AD is a decades-long process, and age is the strongest known risk factor. We previously completed RNA-sequencing and proteomics of several highly tractable, longitudinal, Drosophila models of tauopathy to characterize tau- and age-specific perturbations. However, cellular heterogeneity confounds interpretation of bulk tissue data. We thus leverage single cell RNA sequencing (scRNA-seq) to explore cell-specific contributions to tau- and age-mediated gene expression. Whole brains of adult Drosophila pan-neuronally expressing a mutant variant of human Tau (R406W) and controls were dissociated for scRNA-seq (10x Genomics Chromium) at 1-, 10-, and 20- days of age. Perturbations in cellular composition and gene expression were quantified. To address cross-species conservation of tau-mediated scRNAseq signatures, we leverage existing human AD brain scRNA-seq data. As expected, we highlight evidence for tau-mediated degeneration of Kenyon cells, a neuronal population in the Drosophila "mushroom body" that is required for learning and memory, along with a concomitant increase in glial cell numbers (e.g. gliosis). Gene set enrichment analysis reveal overrepresentation of acetylcholine receptor regulation pathways in vulnerable mushroom body neurons. Interestingly, regularized regression analyses highlight innate immune activity as predictive of tau-related neuronal reduction. Cluster-specific transcriptional signatures between flies and humans demonstrate strong correlations among inferred cell identities. Our results comprise a powerful single cell transcriptomic resource for studying tau-mediated disruption of gene expression and dynamic age-dependent changes at cellular resolution in a tractable model system that is amenable for high throughput functional validation.

Endoderm development requires centrioles to restrain p53-mediated apoptosis in the absence of ERK activity.

Centrioles comprise the heart of centrosomes, microtubule-organizing centers. To study the function of centrioles in lung and gut development, we genetically disrupted centrioles throughout the mouse endoderm. Surprisingly, removing centrioles from the endoderm did not disrupt intestinal growth or development but blocked lung branching. In the lung, acentriolar SOX2-expressing airway epithelial cells apoptosed. Loss of centrioles activated p53, and removing p53 restored survival of SOX2-expressing cells, lung branching, and mouse viability. To investigate how endodermal p53 activation specifically killed acentriolar SOX2-expressing cells, we assessed ERK, a prosurvival cue. ERK was active throughout the intestine and in the distal lung buds, correlating with tolerance to centriole loss. Pharmacologically inhibiting ERK activated apoptosis in acentriolar cells, revealing that ERK activity protects acentriolar cells from apoptosis. Therefore, centrioles are largely dispensable for endodermal growth and the spatial distribution of ERK activity in the endoderm shapes the developmental consequences of centriolar defects and p53 activation.

Abstract 13644: Depict Analysis Suggests Enrichment of Heart Valve Tissue Expression in African Rheumatic Heart Disease Patients From the RHDGen GWAS Study Cohort

Introduction: Rheumatic heart disease (RHD) is a sequela of rheumatic fever that results in permanent heart valve damage. While it remains the primary cause of cardiac surgery in Africa, its pathophysiology is poorly understood.Genome-wide association studies (GWASs) are ‘hypothesis-free’, genetic discovery tools for common complex diseases. Recently, our African GWAS of RHD susceptibility uncovered a novel locus (11q24.1; OR=1.63; 95% CI, 1.38-1.94; P=4.36x10 -8 ) which may explain the disproportionate burden in Africans.Here we further report secondary analyses investigating tissue type expression enrichment among suggestively associated loci, to elucidate the underlying biological processes involved in RHD susceptibility. Aims and Objectives: This study sought to identify tissue types enriched in gene-set analyses associated with RHD GWAS susceptibility data in Africans. Methods: A multicenter, ethnically matched case-control GWAS was performed in 2,548 RHD cases and 2,261 controls, from eight African countries. All 4,809 participants were assessed according to the 2012 World Heart Federation (WHF) criteria for echocardiographic diagnosis of RHD. After cleaning, the GWAS data were stratified by ethnicity and an association test was run via mixed linear models in GCTA that adjusted for gender, chip type and the first ten principal components. Using the resulting GWAS summary statistics suggestive loci (P<1x10 -5 ), we conducted secondary analyses in DEPICT (Data-driven Expression Prioritized Integration for Complex Traits) to test for association with RHD and its related traits and tissues, to validate the GWAS findings. Results: DEPICT indicated promising and biologically relevant results; the most enriched tissues were heart valves (p=2.47x10 -3 ; False Discovery Rate (FDR)=<0.20). Discussion and Conclusion: Our results indicate that genetic susceptibility to RHD may involve perturbed gene expression signatures characteristic of heart valves, known for their pathogenic involvement in RHD development. Future larger studies are warranted to identify the specific genes involved.

Potentiating antibiotic efficacy via perturbation of non-essential gene expression

Proliferation of multidrug-resistant (MDR) bacteria poses a threat to human health, requiring new strategies. Here we propose using fitness neutral gene expression perturbations to potentiate antibiotics. We systematically explored 270 gene knockout-antibiotic combinations in Escherichia coli, identifying 90 synergistic interactions. Identified gene targets were subsequently tested for antibiotic synergy on the transcriptomic level via multiplexed CRISPR-dCas9 and showed successful sensitization of E. coli without a separate fitness cost. These fitness neutral gene perturbations worked as co-therapies in reducing a Salmonella enterica intracellular infection in HeLa. Finally, these results informed the design of four antisense peptide nucleic acid (PNA) co-therapies, csgD, fnr, recA and acrA, against four MDR, clinically isolated bacteria. PNA combined with sub-minimal inhibitory concentrations of trimethoprim against two isolates of Klebsiella pneumoniae and E. coli showed three cases of re-sensitization with minimal fitness impacts. Our results highlight a promising approach for extending the utility of current antibiotics.

HI-NESS: a family of genetically encoded DNA labels based on a bacterial nucleoid-associated protein.

The interplay between three-dimensional chromosome organisation and genomic processes such as replication and transcription necessitates in vivo studies of chromosome dynamics. Fluorescent organic dyes are often used for chromosome labelling in vivo. The mode of binding of these dyes to DNA cause its distortion, elongation, and partial unwinding. The structural changes induce DNA damage and interfere with the binding dynamics of chromatin-associated proteins, consequently perturbing gene expression, genome replication, and cell cycle progression. We have developed a minimally-perturbing, genetically encoded fluorescent DNA label consisting of a (photo-switchable) fluorescent protein fused to the DNA-binding domain of H-NS — a bacterial nucleoid-associated protein. We show that this DNA label, abbreviated as HI-NESS (H-NS-based indicator for nucleic acid stainings), is minimally-perturbing to genomic processes and labels chromosomes in eukaryotic cells in culture, and in zebrafish embryos with preferential binding to AT-rich chromatin.

Identification of Pathway-Based Biomarkers with Crosstalk Analysis for Overall Survival Risk Prediction in Breast Cancer

Recently, many studies have investigated the role of gene-signature on the prognostic assessment of breast cancer (BC), however, the tumor heterogeneity and sequencing noise have limited the clinical usage of these models. Pathway-based approaches are more stable to the perturbation of certain gene expression. In this study, we constructed a prognostic classifier based on survival-related pathway crosstalk analysis. We estimated pathway’s deregulation scores (PDSs) for samples collected from public databases to select survival-related pathways. After pathway crosstalk analysis, we conducted K-means clustering analysis to cluster the patients into G1 and G2 subgroups. The survival outcome of the G2 subgroup was significantly worse than the G1 subgroup. Internal and external dataset exhibits high consistency with the training dataset. Significant differences were found between G2 and G1 subgroups on pathway activity, gene mutation, immune cell infiltration levels, and in particular immune cells/pathway’s activities were significantly negatively associated with BC patient’s outcomes. In conclusion, we established a novel classifier reflecting the overall survival risk of BC and successfully validated its clinical usage on multiple BC datasets, which could offer clinicians inspiration in formulating the clinical treatment plan.

Agent Repurposing for the Treatment of Advanced Stage Diffuse Large B-Cell Lymphoma Based on Gene Expression and Network Perturbation Analysis.

Over 50% of diffuse large B-cell lymphoma (DLBCL) patients are diagnosed at an advanced stage. Although there are a few therapeutic strategies for DLBCL, most of them are more effective in limited-stage cancer patients. The prognosis of patients with advanced-stage DLBCL is usually poor with frequent recurrence and metastasis. In this study, we aimed to identify gene expression and network differences between limited- and advanced-stage DLBCL patients, with the goal of identifying potential agents that could be used to relieve the severity of DLBCL. Specifically, RNA sequencing data of DLBCL patients at different clinical stages were collected from the cancer genome atlas (TCGA). Differentially expressed genes were identified using DESeq2, and then, weighted gene correlation network analysis (WGCNA) and differential module analysis were performed to find variations between different stages. In addition, important genes were extracted by key driver analysis, and potential agents for DLBCL were identified according to gene-expression perturbations and the Crowd Extracted Expression of Differential Signatures (CREEDS) drug signature database. As a result, 20 up-regulated and 73 down-regulated genes were identified and 79 gene co-expression modules were found using WGCNA, among which, the thistle1 module was highly related to the clinical stage of DLBCL. KEGG pathway and GO enrichment analyses of genes in the thistle1 module indicated that DLBCL progression was mainly related to the NOD-like receptor signaling pathway, neutrophil activation, secretory granule membrane, and carboxylic acid binding. A total of 47 key drivers were identified through key driver analysis with 11 up-regulated key driver genes and 36 down-regulated key diver genes in advanced-stage DLBCL patients. Five genes (MMP1, RAB6C, ACCSL, RGS21 and MOCOS) appeared as hub genes, being closely related to the occurrence and development of DLBCL. Finally, both differentially expressed genes and key driver genes were subjected to CREEDS analysis, and 10 potential agents were predicted to have the potential for application in advanced-stage DLBCL patients. In conclusion, we propose a novel pipeline to utilize perturbed gene-expression signatures during DLBCL progression for identifying agents, and we successfully utilized this approach to generate a list of promising compounds.

Resolving the intertwining of inflammation and fibrosis in human heart failure at single-cell level.

Inflammation and fibrosis are intertwined mechanisms fundamentally involved in heart failure. Detailed deciphering gene expression perturbations and cell-cell interactions of leukocytes and non-myocytes is required to understand cell-type-specific pathology in the failing human myocardium. To this end, we performed single-cell RNA sequencing and single T cell receptor sequencing of 200,615 cells in both human dilated cardiomyopathy (DCM) and ischemic cardiomyopathy (ICM) hearts. We sampled both lesion and mild-lesion tissues from each heart to sequentially capture cellular and molecular alterations to different extents of cardiac fibrosis. By which, left (lesion) and right ventricle (mild-lesion) for DCM hearts were harvest while infarcted (lesion) and non-infarcted area (mild-lesion) were dissected from ICM hearts. A novel transcription factor AEBP1 was identified as a crucial cardiac fibrosis regulator in ACTA2+ myofibroblasts. Within fibrotic myocardium, an infiltration of a considerable number of leukocytes was witnessed, especially cytotoxic and exhausted CD8+ T cells and pro-inflammatory CD4+ T cells. Furthermore, a subset of tissue-resident macrophage, CXCL8hiCCR2+HLA-DRhi macrophage was particularly identified in severely fibrotic area, which interacted with activated endothelial cell via DARC, that potentially facilitate leukocyte recruitment and infiltration in human heart failure.

Multi-species single-cell transcriptomic analysis of ocular compartment regulons

The retina is a widely profiled tissue in multiple species by single-cell RNA sequencing studies. However, integrative research of the retina across species is lacking. Here, we construct the first single-cell atlas of the human and porcine ocular compartments and study inter-species differences in the retina. In addition to that, we identify putative adult stem cells present in the iris tissue. We also create a disease map of genes involved in eye disorders across compartments of the eye. Furthermore, we probe the regulons of different cell populations, which include transcription factors and receptor-ligand interactions and reveal unique directional signalling between ocular cell types. In addition, we study conservation of regulons across vertebrates and zebrafish to identify common core factors. Here, we show perturbation of KLF7 gene expression during retinal ganglion cells differentiation and conclude that it plays a significant role in the maturation of retinal ganglion cells.

Repurposable drugs for SARS-CoV-2 and influenza sepsis with scRNA-seq data targeting post-transcription modifications.

Coronavirus disease 2019 (COVID-19) has impacted almost every part of human life worldwide, posing a massive threat to human health. The lack of time for new drug discovery and the urgent need for rapid disease control to reduce mortality have led to a search for quick and effective alternatives to novel therapeutics, for example drug repurposing. To identify potentially repurposable drugs, we employed a systematic approach to mine candidates from U.S. FDA-approved drugs and preclinical small-molecule compounds by integrating gene expression perturbation data for chemicals from the Library of Integrated Network-Based Cellular Signatures project with a publicly available single-cell RNA sequencing dataset from patients with mild and severe COVID-19 (GEO: GSE145926, public data available and accessed on 22 April 2020). We identified 281 FDA-approved drugs that have the potential to be effective against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, 16 of which are currently undergoing clinical trials to evaluate their efficacy against COVID-19. We experimentally tested and demonstrated the inhibitory effects of tyrphostin-AG-1478 and brefeldin-a, two chemical inhibitors of glycosylation (a post-translational modification) on the replication of the single-stranded ribonucleic acid (ssRNA) virus influenza A virus as well as on the transcription and translation of host cell cytokines and their regulators (IFNs and ISGs). In conclusion, we have identified and experimentally validated repurposable anti-SARS-CoV-2 and IAV drugs using a systems biology approach, which may have the potential for treating these viral infections and their complications (sepsis).