Prime-seq, efficient and powerful bulk RNA sequencing

Integrating nontargeted metabolomics and RNA sequencing of dexamethasone-treated and untreated asthmatic mice reveals changes of amino acids and aminoacyl-tRNA in group 2 innate lymphoid cells

Bronchial asthma is the most common multifactorial and heterogeneous disease in childhood. The glucocorticoid dexamethasone is a classic treatment for asthma. Research indicates that group 2 innate lymphoid cells (ILC2s) are crucial to the pathogenesis of asthma. However, few studies have focused on ILC2s metabolism and transcription. This study aims to establish an ovalbumin (OVA)-induced asthma model and a dexamethasone-treated asthma model to explore the regulation of lung ILC2s at the genetic and metabolic levels during the progression and remission of asthma, utilizing single-cell metabolomics and transcriptomics approaches. The results showed that ILC2s regulated the metabolic pathways and transcriptional levels of amino acids (such as arginine, proline, and histidine) and linoleic acid, as well as the metabolic biomarkers of arginine, urocanic acid, and linoleic acid in asthma. Additionally, the cytokine pathways and NF-γB pathways have been altered at the genetic level. At the same time, we revealed that dexamethasone regulates ILC2s amino acid and aminoacyl tRNA metabolism, as well as related genes, thereby alleviating asthma symptoms. Furthermore, we identified the genes Eno3 and Tap1, which are significantly associated with asthma. Quantitative real-time polymerase chain reaction (qRT-PCR) was performed to verify the accuracy of the RNA sequencing results. This study, for the first time, revealed the mechanistic changes of ILC2s in the development and treatment of asthma using multiomics techniques, laying a foundation for targeted therapies in asthma.

The glucocorticoid receptor elicited proliferative response in human erythropoiesis is BCL11A-dependent

Prior evidence indicates that the erythroid cellular response to glucocorticoids (GC) has developmental specificity, namely, that developmentally more advanced cells that are undergoing or have undergone fetal to adult globin switching are more responsive to GC-induced expansion. To investigate the molecular underpinnings of this, we focused on the major developmental globin regulator BCL11A. We compared: (1) levels of expression and nuclear content of BCL11A in adult erythroid cells upon GC stimulation; (2) response to GC of CD34+ cells from patients with BCL11A microdeletions and reduced BCL11A expression, and; (3) response to GC of 2 cellular models (HUDEP-2 and adult CD34+ cells) before and after reduction of BCL11A expression by shRNA. We observed that: (1) GC-expanded erythroid cells from a large cohort of blood donors displayed amplified expression and nuclear accumulation of BCL11A; (2) CD34 + cells from BCL11A microdeletion patients generated fewer erythroid cells when cultured with GC compared to their parents, while the erythroid expansion of the patients was similar to that of their parents in cultures without GC, and; (3) adult CD34+ cells and HUDEP-2 cells with shRNA-depleted expression of BCL11A exhibit reduced expansion in response to GC. In addition, RNA-seq profiling of shRNA-BCL11A CD34+ cells cultured with and without GC was similar (very few differentially expressed genes), while GC-specific responses (differential expression of GILZ and of numerous additional genes) were observed only in control cells with unperturbed BCL11A expression. These data indicate that BCL11A is an important participant in certain aspects of the stress pathway sustained by GC.

A deep learning method to integrate extracelluar miRNA with mRNA for cancer studies

MOTIVATION

Extracellular miRNAs (exmiRs) and intracellular mRNAs both can serve as promising biomarkers and therapeutic targets for various diseases. However, exmiR expression data is often noisy, and obtaining intracellular mRNA expression data usually involves intrusive procedures. To gain valuable insights into disease mechanisms, it is thus essential to improve the quality of exmiR expression data and develop noninvasive methods for assessing intracellular mRNA expression.

RESULTS

We developed CrossPred, a deep-learning multi-encoder model for the cross-prediction of exmiRs and mRNAs. Utilizing contrastive learning, we created a shared embedding space to integrate exmiRs and mRNAs. This shared embedding was then used to predict intracellular mRNA expression from noisy exmiR data and to predict exmiR expression from intracellular mRNA data. We evaluated CrossPred on three types of cancers and assessed its effectiveness in predicting the expression levels of exmiRs and mRNAs. CrossPred outperformed the baseline encoder-decoder model, exmiR or mRNA-based models, and variational autoencoder models. Moreover, the integration of exmiR and mRNA data uncovered important exmiRs and mRNAs associated with cancer. Our study offers new insights into the bidirectional relationship between mRNAs and exmiRs.

AVAILABILITY AND IMPLEMENTATION

The datasets and tool are available at https://doi.org/10.5281/zenodo.13891508.

Evidence for compensatory evolution within pleiotropic regulatory elements

Pleiotropy, measured as expression breadth across tissues, is one of the best predictors for protein sequence and expression conservation. In this study, we investigated its effect on the evolution of cis-regulatory elements (CREs). To this end, we carefully reanalyzed the Epigenomics Roadmap data for nine fetal tissues, assigning a measure of pleiotropic degree to nearly half a million CREs. To assess the functional conservation of CREs, we generated ATAC-seq and RNA-seq data from humans and macaques. We found that more pleiotropic CREs exhibit greater conservation in accessibility, and the mRNA expression levels of the associated genes are more conserved. This trend of higher conservation for higher degrees of pleiotropy persists when analyzing the transcription factor binding repertoire. In contrast, simple DNA sequence conservation of orthologous sites between species tends to be even lower for pleiotropic CREs than for species-specific CREs. Combining various lines of evidence, we propose that the lack of sequence conservation in functionally conserved pleiotropic CREs is owing to within-element compensatory evolution. In summary, our findings suggest that pleiotropy is also a good predictor for the functional conservation of CREs, even though this is not reflected in the sequence conservation of pleiotropic CREs.

Dissecting the Impact of Maternal Androgen Exposure on Developmental Programming through Targeting the Androgen Receptor

Women with polycystic ovary syndrome (PCOS) exhibit sustained elevation in circulating androgens during pregnancy, an independent risk factor linked to pregnancy complications and adverse outcomes in offspring. Yet, further studies are required to understand the effects of elevated androgens on cell type-specific placental dysfunction and fetal development. Therefore, a PCOS-like mouse model induced by continuous androgen exposure is examined. The PCOS-mice exhibited impaired placental and embryonic development, resulting in mid-gestation lethality. Co-treatment with the androgen receptor blocker, flutamide, prevents these phenotypes including germ cell specification. Comprehensive profiling of the placenta by whole-genome bisulfite and RNA sequencing shows a reduced proportion of trophoblast precursors, possibly due to the downregulation of Cdx2 expression. Reduced expression of Gcm1, Synb, and Prl3b1 is associated with reduced syncytiotrophoblasts and sinusoidal trophoblast giant cells, impairs placental labyrinth formation. Importantly, human trophoblast organoids exposed to androgens exhibit analogous changes, showing impaired trophoblast differentiation as a key feature in PCOS-related pregnancy complications. These findings provide new insights into the potential cellular targets for future treatments.

siqRNA-seq is a spike-in-independent technique for quantitative mapping of mRNA landscape

BACKGROUND

RNA sequencing (RNA-seq) is widely used for gene expression profiling and quantification. Quantitative RNA sequencing usually requires cell counting and spike-in, which is not always applicable to many samples. Here, we present a novel quantitative RNA sequencing method independent of spike-ins or cell counting, named siqRNA-seq, which can be used to quantitatively profile gene expression by utilizing gDNA as an internal control. Single-stranded library preparation used in siqRNA-seq profiles gDNA and cDNA with equal efficiency.

RESULTS

To quantify mRNA expression levels, siqRNA-seq constructs libraries for total nucleic acid to establish a model for expression quantification. Compared to Relative Quantification RNA-seq, siqRNA-seq is technically reliable and reproducible for expression profiling but also can sequence reads from gDNA which can be used as an internal reference for accurate expression quantification. Applying siqRNA-seq to investigate the effects of actinomycin D on gene expression in HEK293T cells, we show the advantages of siqRNA-seq in accurately identifying differentially expressed genes between samples with distinct global mRNA levels. Furthermore, we analyzed factors influencing the downward trend of gene expression regulated by ActD using siqRNA-seq and found that mRNA with m6A modification exhibited a faster decay rate compared to mRNA without m6A modification. Additionally, applying this technique to the quantitative analysis of seven tumor cell lines revealed a high degree of diversity in total mRNA expression among tumor cell lines.

CONCLUSIONS

Collectively, siqRNA-seq is a spike-in independent quantitative RNA sequencing method, which creatively uses gDNA as an internal reference to absolutely quantify gene expression. We consider that siqRNA-seq provides a convenient and versatile method to quantitatively profile the mRNA landscape in various samples.

Direct neuronal reprogramming of mouse astrocytes is associated with multiscale epigenome remodeling and requires Yy1

Direct neuronal reprogramming is a promising approach to regenerate neurons from local glial cells. However, mechanisms of epigenome remodeling and co-factors facilitating this process are unclear. In this study, we combined single-cell multiomics with genome-wide profiling of three-dimensional nuclear architecture and DNA methylation in mouse astrocyte-to-neuron reprogramming mediated by Neurogenin2 (Ngn2) and its phosphorylation-resistant form (PmutNgn2), respectively. We show that Ngn2 drives multilayered chromatin remodeling at dynamic enhancer-gene interaction sites. PmutNgn2 leads to higher reprogramming efficiency and enhances epigenetic remodeling associated with neuronal maturation. However, the differences in binding sites or downstream gene activation cannot fully explain this effect. Instead, we identified Yy1, a transcriptional co-factor recruited by direct interaction with Ngn2 to its target sites. Upon deletion of Yy1, activation of neuronal enhancers, genes and ultimately reprogramming are impaired without affecting Ngn2 binding. Thus, our work highlights the key role of interactors of proneural factors in direct neuronal reprogramming.

Generation and characterization of inducible KRAB-dCas9 iPSCs from primates for cross-species CRISPRi

Comparisons of molecular phenotypes across primates provide unique information to understand human biology and evolution, and single-cell RNA-seq CRISPR interference (CRISPRi) screens are a powerful approach to analyze them. Here, we generate and validate three human, three gorilla, and two cynomolgus iPS cell lines that carry a dox-inducible KRAB-dCas9 construct at the AAVS1 locus. We show that despite variable expression levels of KRAB-dCas9 among lines, comparable downregulation of target genes and comparable phenotypic effects are observed in a single-cell RNA-seq CRISPRi screen. Hence, we provide valuable resources for performing and further extending CRISPRi in human and non-human primates.

Multiomic analyses uncover immunological signatures in acute and chronic coronary syndromes

Acute and chronic coronary syndromes (ACS and CCS) are leading causes of mortality. Inflammation is considered a key pathogenic driver of these diseases, but the underlying immune states and their clinical implications remain poorly understood. Multiomic factor analysis (MOFA) allows unsupervised data exploration across multiple data types, identifying major axes of variation and associating these with underlying molecular processes. We hypothesized that applying MOFA to multiomic data obtained from blood might uncover hidden sources of variance and provide pathophysiological insights linked to clinical needs. Here we compile a longitudinal multiomic dataset of the systemic immune landscape in both ACS and CCS (n = 62 patients in total, n = 15 women and n = 47 men) and validate this in an external cohort (n = 55 patients in total, n = 11 women and n = 44 men). MOFA reveals multicellular immune signatures characterized by distinct monocyte, natural killer and T cell substates and immune-communication pathways that explain a large proportion of inter-patient variance. We also identify specific factors that reflect disease state or associate with treatment outcome in ACS as measured using left ventricular ejection fraction. Hence, this study provides proof-of-concept evidence for the ability of MOFA to uncover multicellular immune programs in cardiovascular disease, opening new directions for mechanistic, biomarker and therapeutic studies.

Direct neuronal reprogramming of NDUFS4 patient cells identifies the unfolded protein response as a novel general reprogramming hurdle

Mitochondria account for essential cellular pathways, from ATP production to nucleotide metabolism, and their deficits lead to neurological disorders and contribute to the onset of age-related diseases. Direct neuronal reprogramming aims at replacing neurons lost in such conditions, but very little is known about the impact of mitochondrial dysfunction on the direct reprogramming of human cells. Here, we explore the effects of mitochondrial dysfunction on the neuronal reprogramming of induced pluripotent stem cell (iPSC)-derived astrocytes carrying mutations in the NDUFS4 gene, important for Complex I and associated with Leigh syndrome. This led to the identification of the unfolded protein response as a major hurdle in the direct neuronal conversion of not only astrocytes and fibroblasts from patients but also control human astrocytes and fibroblasts. Its transient inhibition potently improves reprogramming by influencing the mitochondria-endoplasmic-reticulum-stress-mediated pathways. Taken together, disease modeling using patient cells unraveled novel general hurdles and ways to overcome these in human astrocyte-to-neuron reprogramming.

Appropriate glycemic management protects the germline but not the uterine environment in hyperglycemia

Emerging evidence indicates that parental diseases can impact the health of subsequent generations through epigenetic inheritance. Recently, it was shown that maternal diabetes alters the metaphase II oocyte transcriptome, causing metabolic dysfunction in offspring. However, type 1 diabetes (T1D) mouse models frequently utilized in previous studies may be subject to several confounding factors due to severe hyperglycemia. This limits clinical translatability given improvements in glycemic control for T1D subjects. Here, we optimize a T1D mouse model to investigate the effects of appropriately managed maternal glycemic levels on oocytes and intrauterine development. We show that diabetic mice with appropriate glycemic control exhibit better long-term health, including maintenance of the oocyte transcriptome and chromatin accessibility. We further show that human oocytes undergoing in vitro maturation challenged with mildly increased levels of glucose, reflecting appropriate glycemic management, also retain their transcriptome. However, fetal growth and placental function are affected in mice despite appropriate glycemic control, suggesting the uterine environment rather than the germline as a pathological factor in developmental programming in appropriately managed diabetes.

Peripheral priming induces plastic transcriptomic and proteomic responses in circulating neutrophils required for pathogen containment

Neutrophils rapidly respond to inflammation and infection, but to which degree their functional trajectories after mobilization from the bone marrow are shaped within the circulation remains vague. Experimental limitations have so far hampered neutrophil research in human disease. Here, using innovative fixation and single-cell-based toolsets, we profile human and murine neutrophil transcriptomes and proteomes during steady state and bacterial infection. We find that peripheral priming of circulating neutrophils leads to dynamic shifts dominated by conserved up-regulation of antimicrobial genes across neutrophil substates, facilitating pathogen containment. We show the TLR4/NF-κB signaling-dependent up-regulation of canonical neutrophil activation markers like CD177/NB-1 during acute inflammation, resulting in functional shifts in vivo. Blocking de novo RNA synthesis in circulating neutrophils abrogates these plastic shifts and prevents the adaptation of antibacterial neutrophil programs by up-regulation of distinct effector molecules upon infection. These data underline transcriptional plasticity as a relevant mechanism of functional neutrophil reprogramming during acute infection to foster bacterial containment within the circulation.

The glucocorticoid receptor elicited proliferative response in human erythropoiesis is BCL11A-dependent

Prior evidence indicates that the erythroid cellular response to glucocorticoids (GC) has developmental specificity, namely, that developmentally more advanced cells that are undergoing or have undergone fetal to adult globin switching are more responsive to GC-induced expansion. To investigate the molecular underpinnings of this, we focused on the major developmental globin regulator BCL11A. We compared: a) levels of expression and nuclear content of BCL11A in adult erythroid cells upon GC stimulation; b) response to GC of CD34+ cells from patients with BCL11A microdeletions and reduced BCL11A expression, and; c) response to GC of two cellular models (HUDEP-2 and adult CD34+ cells) before and after reduction of BCL11A expression by shRNA. We observed that: a) GC-expanded erythroid cells from a large cohort of blood donors displayed amplified expression and nuclear accumulation of BCL11A; b) CD34+ cells from BCL11A microdeletion patients generated fewer erythroid cells when cultured with GC compared to their parents, while the erythroid expansion of the patients was similar to that of their parents in cultures without GC, and; c) adult CD34+ cells and HUDEP-2 cells with shRNA-depleted expression of BCL11A exhibit reduced expansion in response to GC. In addition, RNA-seq profiling of shRNA-BCL11A CD34+ cells cultured with and without GC was similar (very few differentially expressed genes), while GC-specific responses (differential expression of GILZ and of numerous additional genes) were observed only in controls cells with unperturbed BCL11A expression. These data indicate that BCL11A is an important participant of certain aspects of the stress pathway sustained by GC.

Integrated analysis reveals crosstalk between pyroptosis and immune regulation in renal fibrosis

Purpose

The pathogenesis of renal fibrosis (RF) involves intricate interactions between profibrotic processes and immune responses. This study aimed to explore the potential involvement of the pyroptosis signaling pathway in immune microenvironment regulation within the context of RF. Through comprehensive bioinformatics analysis and experimental validation, we investigated the influence of pyroptosis on the immune landscape in RF.

Methods

We obtained RNA-seq datasets from Gene Expression Omnibus (GEO) databases and identified Pyroptosis-Associated Regulators (PARs) through literature reviews. Systematic evaluation of alterations in 27 PARs was performed in RF and normal kidney samples, followed by relevant functional analyses. Unsupervised cluster analysis revealed distinct pyroptosis modification patterns. Using single-sample gene set enrichment analysis (ssGSEA), we examined the correlation between pyroptosis and immune infiltration. Hub regulators were identified via weighted gene coexpression network analysis (WGCNA) and further validated in a single-cell RNA-seq dataset. We also established a unilateral ureteral obstruction-induced RF mouse model to verify the expression of key regulators at the mRNA and protein levels.

Results

Our comprehensive analysis revealed altered expression of 19 PARs in RF samples compared to normal samples. Five hub regulators, namely PYCARD, CASP1, AIM2, NOD2, and CASP9, exhibited potential as biomarkers for RF. Based on these regulators, a classifier capable of distinguishing normal samples from RF samples was developed. Furthermore, we identified correlations between immune features and PARs expression, with PYCARD positively associated with regulatory T cells abundance in fibrotic tissues. Unsupervised clustering of RF samples yielded two distinct subtypes (Subtype A and Subtype B), with Subtype B characterized by active immune responses against RF. Subsequent WGCNA analysis identified PYCARD, CASP1, and NOD2 as hub PARs in the pyroptosis modification patterns. Single-cell level validation confirmed PYCARD expression in myofibroblasts, implicating its significance in the stress response of myofibroblasts to injury. In vivo experimental validation further demonstrated elevated PYCARD expression in RF, accompanied by infiltration of Foxp3+ regulatory T cells.

Conclusions

Our findings suggest that pyroptosis plays a pivotal role in orchestrating the immune microenvironment of RF. This study provides valuable insights into the pathogenesis of RF and highlights potential targets for future therapeutic interventions.

Targeting a cell-specific microRNA repressor of CXCR4 ameliorates atherosclerosis in mice

The CXC chemokine receptor 4 (CXCR4) in endothelial cells (ECs) and vascular smooth muscle cells (VSMCs) is crucial for vascular integrity. The atheroprotective functions of CXCR4 in vascular cells may be counteracted by atherogenic functions in other nonvascular cell types. Thus, strategies for cell-specifically augmenting CXCR4 function in vascular cells are crucial if this receptor is to be useful as a therapeutic target in treating atherosclerosis and other vascular disorders. Here, we identified miR-206-3p as a vascular-specific CXCR4 repressor and exploited a target-site blocker (CXCR4-TSB) that disrupted the interaction of miR-206-3p with CXCR4 in vitro and in vivo. In vitro, CXCR4-TSB enhanced CXCR4 expression in human and murine ECs and VSMCs to modulate cell viability, proliferation, and migration. Systemic administration of CXCR4-TSB in Apoe-deficient mice enhanced Cxcr4 expression in ECs and VSMCs in the walls of blood vessels, reduced vascular permeability and monocyte adhesion to endothelium, and attenuated the development of diet-induced atherosclerosis. CXCR4-TSB also increased CXCR4 expression in B cells, corroborating its atheroprotective role in this cell type. Analyses of human atherosclerotic plaque specimens revealed a decrease in CXCR4 and an increase in miR-206-3p expression in advanced compared with early lesions, supporting a role for the miR-206-3p-CXCR4 interaction in human disease. Disrupting the miR-206-3p-CXCR4 interaction in a cell-specific manner with target-site blockers is a potential therapeutic approach that could be used to treat atherosclerosis and other vascular diseases.

Mural cell-derived chemokines provide a protective niche to safeguard vascular macrophages and limit chronic inflammation

Maladaptive, non-resolving inflammation contributes to chronic inflammatory diseases such as atherosclerosis. Because macrophages remove necrotic cells, defective macrophage programs can promote chronic inflammation with persistent tissue injury. Here, we investigated the mechanisms sustaining vascular macrophages. Intravital imaging revealed a spatiotemporal macrophage niche across vascular beds alongside mural cells (MCs)-pericytes and smooth muscle cells. Single-cell transcriptomics, co-culture, and genetic deletion experiments revealed MC-derived expression of the chemokines CCL2 and MIF, which actively preserved macrophage survival and their homeostatic functions. In atherosclerosis, this positioned macrophages in viable plaque areas, away from the necrotic core, and maintained a homeostatic macrophage phenotype. Disruption of this MC-macrophage unit via MC-specific deletion of these chemokines triggered detrimental macrophage relocalizing, exacerbated plaque necrosis, inflammation, and atheroprogression. In line, CCL2 inhibition at advanced stages of atherosclerosis showed detrimental effects. This work presents a MC-driven safeguard toward maintaining the homeostatic vascular macrophage niche.

CSF3R T618I Collaborates With RUNX1-RUNX1T1 to Expand Hematopoietic Progenitors and Sensitizes to GLI Inhibition

Activating colony-stimulating factor-3 receptor gene (CSF3R) mutations are recurrent in acute myeloid leukemia (AML) with t(8;21) translocation. However, the nature of oncogenic collaboration between alterations of CSF3R and the t(8;21) associated RUNX1-RUNX1T1 fusion remains unclear. In CD34+ hematopoietic stem and progenitor cells from healthy donors, double oncogene expression led to a clonal advantage, increased self-renewal potential, and blast-like morphology and distinct immunophenotype. Gene expression profiling revealed hedgehog signaling as a potential mechanism, with upregulation of GLI2 constituting a putative pharmacological target. Both primary hematopoietic cells and the t(8;21) positive AML cell line SKNO-1 showed increased sensitivity to the GLI inhibitor GANT61 when expressing CSF3R T618I. Our findings suggest that during leukemogenesis, the RUNX1-RUNXT1 fusion and CSF3R mutation act in a synergistic manner to alter hedgehog signaling, which can be exploited therapeutically.

3' RNA-seq is superior to standard RNA-seq in cases of sparse data but inferior at identifying toxicity pathways in a model organism

Introduction: The application of RNA-sequencing has led to numerous breakthroughs related to investigating gene expression levels in complex biological systems. Among these are knowledge of how organisms, such as the vertebrate model organism zebrafish (Danio rerio), respond to toxicant exposure. Recently, the development of 3' RNA-seq has allowed for the determination of gene expression levels with a fraction of the required reads compared to standard RNA-seq. While 3' RNA-seq has many advantages, a comparison to standard RNA-seq has not been performed in the context of whole organism toxicity and sparse data. Methods and results: Here, we examined samples from zebrafish exposed to perfluorobutane sulfonamide (FBSA) with either 3' or standard RNA-seq to determine the advantages of each with regards to the identification of functionally enriched pathways. We found that 3' and standard RNA-seq showed specific advantages when focusing on annotated or unannotated regions of the genome. We also found that standard RNA-seq identified more differentially expressed genes (DEGs), but that this advantage disappeared under conditions of sparse data. We also found that standard RNA-seq had a significant advantage in identifying functionally enriched pathways via analysis of DEG lists but that this advantage was minimal when identifying pathways via gene set enrichment analysis of all genes. Conclusions: These results show that each approach has experimental conditions where they may be advantageous. Our observations can help guide others in the choice of 3' RNA-seq vs standard RNA sequencing to query gene expression levels in a range of biological systems.

WDR5 is required for DUX4 expression and its pathological effects in FSHD muscular dystrophy

Facioscapulohumeral muscular dystrophy (FSHD) is one of the most prevalent neuromuscular disorders. The disease is linked to copy number reduction and/or epigenetic alterations of the D4Z4 macrosatellite on chromosome 4q35 and associated with aberrant gain of expression of the transcription factor DUX4, which triggers a pro-apoptotic transcriptional program leading to muscle wasting. As today, no cure or therapeutic option is available to FSHD patients. Given its centrality in FSHD, blocking DUX4 expression with small molecule drugs is an attractive option. We previously showed that the long non protein-coding RNA DBE-T is required for aberrant DUX4 expression in FSHD. Using affinity purification followed by proteomics, here we identified the chromatin remodeling protein WDR5 as a novel DBE-T interactor and a key player required for the biological activity of the lncRNA. We found that WDR5 is required for the expression of DUX4 and its targets in primary FSHD muscle cells. Moreover, targeting WDR5 rescues both cell viability and myogenic differentiation of FSHD patient cells. Notably, comparable results were obtained by pharmacological inhibition of WDR5. Importantly, WDR5 targeting was safe to healthy donor muscle cells. Our results support a pivotal role of WDR5 in the activation of DUX4 expression identifying a druggable target for an innovative therapeutic approach for FSHD.

Nucleocapsid-specific T cell responses associate with control of SARS-CoV-2 in the upper airways before seroconversion

Despite intensive research since the emergence of SARS-CoV-2, it has remained unclear precisely which components of the early immune response protect against the development of severe COVID-19. Here, we perform a comprehensive immunogenetic and virologic analysis of nasopharyngeal and peripheral blood samples obtained during the acute phase of infection with SARS-CoV-2. We find that soluble and transcriptional markers of systemic inflammation peak during the first week after symptom onset and correlate directly with upper airways viral loads (UA-VLs), whereas the contemporaneous frequencies of circulating viral nucleocapsid (NC)-specific CD4⁺ and CD8⁺ T cells correlate inversely with various inflammatory markers and UA-VLs. In addition, we show that high frequencies of activated CD4⁺ and CD8⁺ T cells are present in acutely infected nasopharyngeal tissue, many of which express genes encoding various effector molecules, such as cytotoxic proteins and IFN-γ. The presence of IFNG mRNA-expressing CD4⁺ and CD8⁺ T cells in the infected epithelium is further linked with common patterns of gene expression among virus-susceptible target cells and better local control of SARS-CoV-2. Collectively, these results identify an immune correlate of protection against SARS-CoV-2, which could inform the development of more effective vaccines to combat the acute and chronic illnesses attributable to COVID-19.

Transcriptome-Powered Pluripotent Stem Cell Differentiation for Regenerative Medicine

Pluripotent stem cells are endless sources for in vitro engineering human tissues for regenerative medicine. Extensive studies have demonstrated that transcription factors are the key to stem cell lineage commitment and differentiation efficacy. As the transcription factor profile varies depending on the cell type, global transcriptome analysis through RNA sequencing (RNAseq) has been a powerful tool for measuring and characterizing the success of stem cell differentiation. RNAseq has been utilized to comprehend how gene expression changes as cells differentiate and provide a guide to inducing cellular differentiation based on promoting the expression of specific genes. It has also been utilized to determine the specific cell type. This review highlights RNAseq techniques, tools for RNAseq data interpretation, RNAseq data analytic methods and their utilities, and transcriptomics-enabled human stem cell differentiation. In addition, the review outlines the potential benefits of the transcriptomics-aided discovery of intrinsic factors influencing stem cell lineage commitment, transcriptomics applied to disease physiology studies using patients' induced pluripotent stem cell (iPSC)-derived cells for regenerative medicine, and the future outlook on the technology and its implementation.

Regulatory and coding sequences of TRNP1 co-evolve with brain size and cortical folding in mammals

Brain size and cortical folding have increased and decreased recurrently during mammalian evolution. Identifying genetic elements whose sequence or functional properties co-evolve with these traits can provide unique information on evolutionary and developmental mechanisms. A good candidate for such a comparative approach is TRNP1, as it controls proliferation of neural progenitors in mice and ferrets. Here, we investigate the contribution of both regulatory and coding sequences of TRNP1 to brain size and cortical folding in over 30 mammals. We find that the rate of TRNP1 protein evolution (ω) significantly correlates with brain size, slightly less with cortical folding and much less with body size. This brain correlation is stronger than for >95% of random control proteins. This co-evolution is likely affecting TRNP1 activity, as we find that TRNP1 from species with larger brains and more cortical folding induce higher proliferation rates in neural stem cells. Furthermore, we compare the activity of putative cis-regulatory elements (CREs) of TRNP1 in a massively parallel reporter assay and identify one CRE that likely co-evolves with cortical folding in Old World monkeys and apes. Our analyses indicate that coding and regulatory changes that increased TRNP1 activity were positively selected either as a cause or a consequence of increases in brain size and cortical folding. They also provide an example how phylogenetic approaches can inform biological mechanisms, especially when combined with molecular phenotypes across several species.

UHMK1 is a novel splicing regulatory kinase

The U2AF Homology Motif Kinase 1 (UHMK1) is the only kinase that contains the U2AF homology motif, a common protein interaction domain among splicing factors. Through this motif, UHMK1 interacts with the splicing factors SF1 and SF3B1, known to participate in the 3' splice site recognition during the early steps of spliceosome assembly. Although UHMK1 phosphorylates these splicing factors in vitro, the involvement of UHMK1 in RNA processing has not previously been demonstrated. Here, we identify novel putative substrates of this kinase and evaluate UHMK1 contribution to overall gene expression and splicing, by integrating global phosphoproteomics, RNA-seq, and bioinformatics approaches. Upon UHMK1 modulation, 163 unique phosphosites were differentially phosphorylated in 117 proteins, of which 106 are novel potential substrates of this kinase. Gene Ontology analysis showed enrichment of terms previously associated with UHMK1 function, such as mRNA splicing, cell cycle, cell division, and microtubule organization. The majority of the annotated RNA-related proteins are components of the spliceosome but are also involved in several steps of gene expression. Comprehensive analysis of splicing showed that UHMK1 affected over 270 alternative splicing events. Moreover, splicing reporter assay further supported UHMK1 function on splicing. Overall, RNA-seq data demonstrated that UHMK1 knockdown had a minor impact on transcript expression and pointed to UHMK1 function in epithelial-mesenchymal transition. Functional assays demonstrated that UHMK1 modulation affects proliferation, colony formation, and migration. Taken together, our data implicate UHMK1 as a splicing regulatory kinase, connecting protein regulation through phosphorylation and gene expression in key cellular processes.

High-throughput transcriptional profiling of perturbations by Panax ginseng saponins and Panax notoginseng saponins using TCM-seq

Panax ginseng (PG) and Panax notoginseng (PN) are highly valuable Chinese medicines (CM). Although both CMs have similar active constituents, their clinical applications are clearly different. Over the past decade, RNA sequencing (RNA-seq) analysis has been employed to investigate the molecular mechanisms of extracts or monomers. However, owing to the limited number of samples in standard RNA-seq, few studies have systematically compared the effects of PG and PN spanning multiple conditions at the transcriptomic level. Here, we developed an approach that simultaneously profiles transcriptome changes for multiplexed samples using RNA-seq (TCM-seq), a high-throughput, low-cost workflow to molecularly evaluate CM perturbations. A species-mixing experiment was conducted to illustrate the accuracy of sample multiplexing in TCM-seq. Transcriptomes from repeated samples were used to verify the robustness of TCM-seq. We then focused on the primary active components, Panax notoginseng saponins (PNS) and Panax ginseng saponins (PGS) extracted from PN and PG, respectively. We also characterized the transcriptome changes of 10 cell lines, treated with four different doses of PNS and PGS, using TCM-seq to compare the differences in their perturbing effects on genes, functional pathways, gene modules, and molecular networks. The results of transcriptional data analysis showed that the transcriptional patterns of various cell lines were significantly distinct. PGS exhibited a stronger regulatory effect on genes involved in cardiovascular disease, whereas PNS resulted in a greater coagulation effect on vascular endothelial cells. This study proposes a paradigm to comprehensively explore the differences in mechanisms of action between CMs based on transcriptome readouts.

GPR55 in B cells limits atherosclerosis development and regulates plasma cell maturation

Dissecting the pathways regulating the adaptive immune response in atherosclerosis is of particular therapeutic interest. Here we report that the lipid G-protein coupled receptor GPR55 is highly expressed by splenic plasma cells (PC), upregulated in mouse spleens during atherogenesis and human unstable or ruptured compared to stable plaques. Gpr55-deficient mice developed larger atherosclerotic plaques with increased necrotic core size compared to their corresponding controls. Lack of GPR55 hyperactivated B cells, disturbed PC maturation and resulted in immunoglobulin (Ig)G overproduction. B cell-specific Gpr55 depletion or adoptive transfer of Gpr55-deficient B cells was sufficient to promote plaque development and elevated IgG titers. In vitro, the endogenous GPR55 ligand lysophsophatidylinositol (LPI) enhanced PC proliferation, whereas GPR55 antagonism blocked PC maturation and increased their mitochondrial content. Collectively, these discoveries provide previously undefined evidence for GPR55 in B cells as a key modulator of the adaptive immune response in atherosclerosis.

Machine learning alternative to systems biology should not solely depend on data

In recent years, artificial intelligence (AI)/machine learning has emerged as a plausible alternative to systems biology for the elucidation of biological phenomena and in attaining specified design objective in synthetic biology. Although considered highly disruptive with numerous notable successes so far, we seek to bring attention to both the fundamental and practical pitfalls of their usage, especially in illuminating emergent behaviors from chaotic or stochastic systems in biology. Without deliberating on their suitability and the required data qualities and pre-processing approaches beforehand, the research and development community could experience similar 'AI winters' that had plagued other fields. Instead, we anticipate the integration or combination of the two approaches, where appropriate, moving forward.

PARP14 is a novel target in STAT6 mutant follicular lymphoma

The variable clinical course of follicular lymphoma (FL) is determined by the molecular heterogeneity of tumor cells and complex interactions within the tumor microenvironment (TME). IL-4 producing follicular helper T cells (TFH) are critical components of the FL TME. Binding of IL-4 to IL-4R on FL cells activates JAK/STAT signaling. We identified STAT6 mutations (STAT6MUT) in 13% of FL (N = 33/258), all clustered within the DNA binding domain. Gene expression data and immunohistochemistry showed upregulation of IL-4/STAT6 target genes in STAT6MUT FL, including CCL17, CCL22, and FCER2 (CD23). Functionally, STAT6MUT was gain-of-function by serial replating phenotype in pre-B CFU assays. Expression of STAT6MUT enhanced IL-4 induced FCER2/CD23, CCL17 and CCL22 expression and was associated with nuclear accumulation of pSTAT6. RNA sequencing identified PARP14 -a transcriptional switch and co-activator of STAT6- among the top differentially upregulated genes in IL-4 stimulated STAT6MUT lymphoma cells and in STAT6MUT primary FL cells. Quantitative chromatin immunoprecipitation (qChIP) demonstrated binding of STAT6MUT but not STAT6WT to the PARP14 promotor. Reporter assays showed increased IL-4 induced transactivation activity of STAT6MUT at the PARP14 promotor, suggesting a self-reinforcing regulatory circuit. Knock-down of PARP14 or PARP-inhibition abrogated the STAT6MUT gain-of-function phenotype. Thus, our results identify PARP14 as a novel therapeutic target in STAT6MUT FL.

Combi-seq for multiplexed transcriptome-based profiling of drug combinations using deterministic barcoding in single-cell droplets

Anti-cancer therapies often exhibit only short-term effects. Tumors typically develop drug resistance causing relapses that might be tackled with drug combinations. Identification of the right combination is challenging and would benefit from high-content, high-throughput combinatorial screens directly on patient biopsies. However, such screens require a large amount of material, normally not available from patients. To address these challenges, we present a scalable microfluidic workflow, called Combi-Seq, to screen hundreds of drug combinations in picoliter-size droplets using transcriptome changes as a readout for drug effects. We devise a deterministic combinatorial DNA barcoding approach to encode treatment conditions, enabling the gene expression-based readout of drug effects in a highly multiplexed fashion. We apply Combi-Seq to screen the effect of 420 drug combinations on the transcriptome of K562 cells using only ~250 single cell droplets per condition, to successfully predict synergistic and antagonistic drug pairs, as well as their pathway activities.

TET1 regulates gene expression and repression of endogenous retroviruses independent of DNA demethylation

DNA methylation (5-methylcytosine (5mC)) is critical for genome stability and transcriptional regulation in mammals. The discovery that ten-eleven translocation (TET) proteins catalyze the oxidation of 5mC to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC) revolutionized our perspective on the complexity and regulation of DNA modifications. However, to what extent the regulatory functions of TET1 can be attributed to its catalytic activity remains unclear. Here, we use genome engineering and quantitative multi-omics approaches to dissect the precise catalytic vs. non-catalytic functions of TET1 in murine embryonic stem cells (mESCs). Our study identifies TET1 as an essential interaction hub for multiple chromatin modifying complexes and a global regulator of histone modifications. Strikingly, we find that the majority of transcriptional regulation depends on non-catalytic functions of TET1. In particular, we show that TET1 is critical for the establishment of H3K9me3 and H4K20me3 at endogenous retroviral elements (ERVs) and their silencing that is independent of its canonical role in DNA demethylation. Furthermore, we provide evidence that this repression of ERVs depends on the interaction between TET1 and SIN3A. In summary, we demonstrate that the non-catalytic functions of TET1 are critical for regulation of gene expression and the silencing of endogenous retroviruses in mESCs.

Antibiotic-Induced Treatments Reveal Stress-Responsive Gene Expression in the Endangered Lichen Lobaria pulmonaria

Antibiotics are primarily found in the environment due to human activity, which has been reported to influence the structure of biotic communities and the ecological functions of soil and water ecosystems. Nonetheless, their effects in other terrestrial ecosystems have not been well studied. As a result of oxidative stress in organisms exposed to high levels of antibiotics, genotoxicity can lead to DNA damage and, potentially, cell death. In addition, in symbiotic organisms, removal of the associated microbiome by antibiotic treatment has been observed to have a big impact on the host, e.g., corals. The lung lichen Lobaria pulmonaria has more than 800 associated bacterial species, a microbiome which has been hypothesized to increase the lichen's fitness. We artificially exposed samples of L. pulmonaria to antibiotics and a stepwise temperature increase to determine the relative effects of antibiotic treatments vs. temperature on the mycobiont and photobiont gene expression and the viability and on the community structure of the lichen-associated bacteria. We found that the mycobiont and photobiont highly reacted to different antibiotics, independently of temperature exposure. We did not find major differences in bacterial community composition or alpha diversity between antibiotic treatments and controls. For these reasons, the upregulation of stress-related genes in antibiotic-treated samples could be caused by genotoxicity in L. pulmonaria and its photobiont caused by exposure to antibiotics, and the observed stress responses are reactions of the symbiotic partners to reduce damage to their cells. Our study is of great interest for the community of researchers studying symbiotic organisms as it represents one of the first steps to understanding gene expression in an endangered lichen in response to exposure to toxic environments, along with dynamics in its associated bacterial communities.