Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences

Murine parainfluenza virus persists in lung innate immune cells sustaining chronic lung pathology

Common respiratory viruses, including the human parainfluenza viruses, threaten human health seasonally and associate with the development of chronic lung diseases. Evidence suggests that these viruses can persist, but the sources of viral products in vivo and their impact on chronic respiratory diseases remain unknown. Using the murine parainfluenza virus Sendai, we demonstrate that viral protein and RNA persist in lung macrophages, type 2 innate lymphoid cells (ILC2s) and dendritic cells long after the infectious virus is cleared. Cells containing persistent viral protein expressed Th2 inflammation-related transcriptomic signatures associated with the development of chronic lung diseases, including asthma. Lineage tracing demonstrated that distinct functional groups of cells contribute to the chronic pathology. Importantly, targeted ablation of infected cells significantly ameliorated chronic lung disease. Overall, we identified persistent infection of innate immune cells as a key factor in the progression from acute to chronic lung disease after infection with parainfluenza virus.

Dynamic changes in mRNA nucleocytoplasmic localization in the nitrate response of Arabidopsis roots

Nitrate is a nutrient and signal that regulates gene expression. The nitrate response has been extensively characterized at the organism, organ, and cell-type-specific levels, but intracellular mRNA dynamics remain unexplored. To characterize nuclear and cytoplasmic transcriptome dynamics in response to nitrate, we performed a time-course expression analysis after nitrate treatment in isolated nuclei, cytoplasm, and whole roots. We identified 402 differentially localized transcripts (DLTs) in response to nitrate treatment. Induced DLT genes showed rapid and transient recruitment of the RNA polymerase II, together with an increase in the mRNA turnover rates. DLTs code for genes involved in metabolic processes, localization, and response to stimulus indicating DLTs include genes with relevant functions for the nitrate response that have not been previously identified. Using single-molecule RNA FISH, we observed early nuclear accumulation of the NITRATE REDUCTASE 1 (NIA1) transcripts in their transcription sites. We found that transcription of NIA1, a gene showing delayed cytoplasmic accumulation, is rapidly and transiently activated; however, its transcripts become unstable when they reach the cytoplasm. Our study reveals the dynamic localization of mRNAs between the nucleus and cytoplasm as an emerging feature in the temporal control of gene expression in response to nitrate treatment in Arabidopsis roots.

Transcriptomic Analyses Predict Enhanced Metabolic Activity and Therapeutic Potential of mTOR Inhibitors in Acne-Prone Skin

Current acne therapies center on preventing new lesions in patients with acne. These therapies were historically found to be beneficial yet were chosen without knowledge of the specific changes in the skin that favor lesion development. A major challenge in developing new treatments is the incomplete understanding of nonlesional (NL), acne-prone skin's molecular characteristics. To address this, we compared RNA-sequencing data from NL skin of 49 patients with acne (denoted as NL acne [NLA]) with those from 19 healthy controls with no acne history. We found 77 differentially expressed genes in NLA (log fold change > 1; P < .05), including genes associated with innate immunity and epidermal barrier function. Notably, K RT 6C, K RT 16, S100A8, S100A9, and lactotransferrin were upregulated, and LCE4A, LCE6A, and CTSE were downregulated. Gene set enrichment analysis revealed that metabolic pathways were enriched in NLA skin, whereas keratinization was negatively enriched. To identify compounds that could shift the gene expression signature of NLA skin toward healthy control skin, we performed connectivity mapping with the Library of Integrated Network-Based Signatures. We identified 187 compounds, particularly mTOR inhibitors, that could potentially normalize the gene expression profile of acne-prone skin to that of healthy skin. Our findings indicate that NLA skin has distinct differences in epidermal differentiation, cellular metabolism, and innate immunity that may promote lesion formation and suggest that mTOR inhibitors could restore NLA skin toward a healthier state, potentially reversing the predisposition to lesion development.

Lineage tracing of nuclei in skeletal myofibers uncovers distinct transcripts and interplay between myonuclear populations

Multinucleated skeletal muscle cells need to acquire additional nuclei through fusion with activated skeletal muscle stem cells when responding to both developmental and adaptive growth stimuli. A fundamental question in skeletal muscle biology has been the reason underlying this need for new nuclei in cells that already harbor hundreds of nuclei. Here we utilize nuclear RNA-sequencing approaches and develop a lineage tracing strategy capable of defining the transcriptional state of recently fused nuclei and distinguishing this state from that of pre-existing nuclei. Our findings reveal the presence of conserved markers of newly fused nuclei both during development and after a hypertrophic stimulus in the adult. However, newly fused nuclei also exhibit divergent gene expression that is determined by the myogenic environment to which they fuse. Moreover, accrual of new nuclei through fusion is required for nuclei already resident in adult myofibers to mount a normal transcriptional response to a load-inducing stimulus. We propose a model of mutual regulation in the control of skeletal muscle development and adaptations, where newly fused and pre-existing myonuclear populations influence each other to maintain optimal functional growth.

Cell death induced by Lepeophtheirus salmonis labial gland protein 3 in salmonid fish leukocytes: A mechanism for disabling host immune responses

The salmon louse (Lepeophtheirus salmonis) is an ectoparasite feeding on mucus, skin, and blood of salmonids. On parasitised fish erosions and, at later lice stages, ulcerations appear at the louse feeding site. In susceptible species like Atlantic salmon (Salmo salar) with a limited rejection of lice, only a mild inflammatory response with minor influx of immune cells is seen at these lesions, as the salmon louse secrete proteins that can dampen immune responses. In a previous study, Lepeophtheirus salmonis labial gland protein 3 (LsLGP3) was suggested to dampen cellular responses, and the present study aimed at increasing our understanding of its mode of action. LsLGP3 was found to be secreted on to the host skin, and both in vivo and in vitro experiments were performed to elucidate its function. Histological analysis of the louse attachment site revealed an epidermal and dermal influx of mainly macrophages and granulocytes after 5 days post infestation. The immune cell influx was deeper in the dermis throughout the louse infestation, and LsLGP3 may be involved in dampening this response. Enriched populations of Atlantic salmon B-cells, T-cells, granulocytes, and monocytes were exposed to recombinant LsLGP3 (recLGP3) in vitro, resulting in a significant decrease in cell viability compared to non-exposed controls. An apoptotic cell morphology with "beads-on-a-string" like protrusions was seen in all leukocyte cell fractions after recLGP3 exposure, but not in erythrocytes or keratocytes. A decreased viability was also detected in pink salmon leucocytes, which was not in leucocytes from non-salmonid species. These functional insights suggest that LsLGP3 specifically induces apoptosis of salmonid leukocytes and is likely a key protein secreted by the lice that disables the Atlantic salmon ability to mount an adequate immune response towards the salmon louse. In vivo LsLGP3 knock down studies indicated that the effect is localised primarily at the lice feeding site, without affecting immune cells that are not situated adjacent to the lice-inflicted lesion. The findings from this study could significantly aid in the development of new immune based anti-salmon louse prophylactic measures and treatments.

Microbial dynamics and pulmonary immune responses in COVID-19 secondary bacterial pneumonia

Secondary bacterial pneumonia (2°BP) is associated with significant morbidity following respiratory viral infection, yet remains incompletely understood. In a prospective cohort of 112 critically ill adults intubated for COVID-19, we comparatively assess longitudinal airway microbiome dynamics and the pulmonary transcriptome of patients who developed 2°BP versus controls who did not. We find that 2°BP is significantly associated with both mortality and corticosteroid treatment. The pulmonary microbiome in 2°BP is characterized by increased bacterial RNA mass and dominance of culture-confirmed pathogens, detectable days prior to 2°BP clinical diagnosis, and frequently also present in nasal swabs. Assessment of the pulmonary transcriptome reveals suppressed TNFα signaling in patients with 2°BP, and sensitivity analyses suggest this finding is mediated by corticosteroid treatment. Further, we find that increased bacterial RNA mass correlates with reduced expression of innate and adaptive immunity genes in both 2°BP patients and controls. Taken together, our findings provide fresh insights into the microbial dynamics and host immune features of COVID-19-associated 2°BP, and suggest that suppressed immune signaling, potentially mediated by corticosteroid treatment, permits expansion of opportunistic bacterial pathogens.

Developmental-status-aware transcriptional decomposition establishes a cell state panorama of human cancers

BACKGROUND

Cancer cells evolve under unique functional adaptations that unlock transcriptional programs embedded in adult stem and progenitor-like cells for progression, metastasis, and therapeutic resistance. However, it remains challenging to quantify the stemness-aware cell state of a tumor based on its gene expression profile.

METHODS

We develop a developmental-status-aware transcriptional decomposition strategy using single-cell RNA-sequencing-derived tissue-specific fetal and adult cell signatures as anchors. We apply our method to various biological contexts, including developing human organs, adult human tissues, experimentally induced differentiation cultures, and bulk human tumors, to benchmark its performance and to reveal novel biology of entangled developmental signaling in oncogenic processes.

RESULTS

Our strategy successfully captures complex dynamics in developmental tissue bulks, reveals remarkable cellular heterogeneity in adult tissues, and resolves the ambiguity of cell identities in in vitro transformations. Applying it to large patient cohorts of bulk RNA-seq, we identify clinically relevant cell-of-origin patterns and observe that decomposed fetal cell signals significantly increase in tumors versus normal tissues and metastases versus primary tumors. Across cancer types, the inferred fetal-state strength outperforms published stemness indices in predicting patient survival and confers substantially improved predictive power for therapeutic responses.

CONCLUSIONS

Our study not only provides a general approach to quantifying developmental-status-aware cell states of bulk samples but also constructs an information-rich, biologically interpretable, cell-state panorama of human cancers, enabling diverse translational applications.

DNA-protein quasi-mapping for rapid differential gene expression analysis in non-model organisms

BACKGROUND

Conventional differential gene expression analysis pipelines for non-model organisms require computationally expensive transcriptome assembly. We recently proposed an alternative strategy of directly aligning RNA-seq reads to a protein database, and demonstrated drastic improvements in speed, memory usage, and accuracy in identifying differentially expressed genes.

RESULT

Here we report a further speed-up by replacing DNA-protein alignment by quasi-mapping, making our pipeline > 1000× faster than assembly-based approach, and still more accurate. We also compare quasi-mapping to other mapping techniques, and show that it is faster but at the cost of sensitivity.

CONCLUSION

We provide a quick-and-dirty differential gene expression analysis pipeline for non-model organisms without a reference transcriptome, which directly quasi-maps RNA-seq reads to a reference protein database, avoiding computationally expensive transcriptome assembly.

A transcriptomic hourglass in brown algae

Complex multicellularity has emerged independently across a few eukaryotic lineages and is often associated with the rise of elaborate, tightly coordinated developmental processes1,2. How multicellularity and development are interconnected in evolution is a major question in biology. The hourglass model of embryonic evolution depicts how developmental processes are conserved during evolution, and predicts morphological and molecular divergence in early and late embryogenesis, bridged by a conserved mid-embryonic (phylotypic) period linked to the formation of the basic body plan3,4. Initially found in animal embryos5-8, molecular hourglass patterns have recently been proposed for land plants and fungi9,10. However, whether the hourglass pattern is an intrinsic feature of all complex multicellular eukaryotes remains unknown. Here we tested the presence of a molecular hourglass in the brown algae, a eukaryotic lineage that has evolved multicellularity independently from animals, fungi and plants1,11,12. By exploring transcriptome evolution patterns of brown algae with distinct morphological complexities, we uncovered an hourglass pattern during embryogenesis in morphologically complex species. Filamentous algae without canonical embryogenesis display transcriptome conservation in multicellular stages of the life cycle, whereas unicellular stages are more rapidly evolving. Our findings suggest that transcriptome conservation in brown algae is associated with cell differentiation stages, but is not necessarily linked to embryogenesis. Together with previous work in animals, plants and fungi, we provide further evidence for the generality of a developmental hourglass pattern across complex multicellular eukaryotes.

Vascular endothelial growth factor receptor 2 as a potential host target for the inhibition of enterovirus replication

Because host kinases are key regulators of multiple signaling pathways in response to viral infections, we previously screened a kinase inhibitor library using rhabdomyosarcoma cells and human intestinal organoids in parallel to identify potent inhibitors against EV-A71 infection. We found that Rho-associated coiled-coil-containing protein kinase (Rock) inhibitor efficiently suppressed the EV-A71 replication and further revealed Rock1 as a novel EV-A71 host factor. In this study, subsequent analysis found that a variety of vascular endothelial growth factor receptor (VEGFR) inhibitors also had potent antiviral effects. Among the hits, Pazopanib, with a selectivity index as high as 254, which was even higher than that of Pirodavir, a potent broad-spectrum picornavirus inhibitor targeting viral capsid protein VP1, was selected for further analysis. We demonstrated that Pazopanib not only efficiently suppressed the replication of EV-A71 in a dose-dependent manner, but also exhibited broad-spectrum anti-enterovirus activity. Mechanistically, Pazopanib probably induces alterations in host cells, thereby impeding viral genome replication and transcription. Notably, VEGFR2 knockdown and overexpression suppressed and facilitated EV-A71 replication, respectively, indicating that VEGFR2 is a novel host dependency factor for EV-A71 replication. Transcriptome analysis further proved that VEGFR2 potentially plays a crucial role in combating EV-A71 infection through the TSAd-Src-PI3K-Akt pathway. These findings expand the range of potential antiviral candidates of anti-enterovirus therapeutics and suggest that VEGFR2 may be a key host factor involved in EV-A71 replication, making it a potential target for the development of anti-enterovirus therapeutics.

IMPORTANCE

As the first clinical case was identified in the United States, EV-A71, a significant neurotropic enterovirus, has been a common cause of hand, foot, and mouth disease (HFMD) in infants and young children. Developing an effective antiviral agent for EV-A71 and other human enteroviruses is crucial, as these viral pathogens consistently cause outbreaks in humans. In this study, we demonstrated that multiple inhibitors against VEGFRs effectively reduced EV-A71 replication, with Pazopanib emerging as the top candidate. Furthermore, Pazopanib also attenuated the replication of other enteroviruses, including CVA10, CVB1, EV-D70, and HRV-A, displaying broad-spectrum anti-enterovirus activity. Given that Pazopanib targets various VEGFRs, we narrowed the focus to VEGFR2 using knockdown and overexpression experiments. Transcriptomic analysis suggests that Pazopanib's potential downstream targets involve the TSAd-Src-PI3K-Akt pathway. Our work may contribute to identifying targets for antiviral inhibitors and advancing treatments for human enterovirus infections.

Activation of ATF3 via the integrated stress response pathway regulates innate immune response to restrict Zika virus

Zika virus (ZIKV) is a re-emerging mosquito-borne flavivirus that can have devastating health consequences. The developmental and neurological effects of a ZIKV infection arise in part from the virus triggering cellular stress pathways and perturbing transcriptional programs. To date, the underlying mechanisms of transcriptional control directing viral restriction and virus-host interaction are understudied. Activating Transcription Factor 3 (ATF3) is a stress-induced transcriptional effector that modulates the expression of genes involved in a myriad of cellular processes, including inflammation and antiviral responses, to restore cellular homeostasis. While ATF3 is known to be upregulated during ZIKV infection, the mode by which ATF3 is activated, and the specific role of ATF3 during ZIKV infection is unknown. In this study, we show via inhibitor and RNA interference approaches that ZIKV infection initiates the integrated stress response pathway to activate ATF4 which in turn induces ATF3 expression. Additionally, by using CRISPR-Cas9 system to delete ATF3, we found that ATF3 acts to limit ZIKV gene expression in A549 cells. We also determined that ATF3 enhances the expression of antiviral genes such as STAT1 and other components in the innate immunity pathway to induce an ATF3-dependent anti-ZIKV response. Our study reveals crosstalk between the integrated stress response and innate immune response pathways and highlights an important role for ATF3 in establishing an antiviral effect during ZIKV infection.

IMPORTANCE

Zika virus (ZIKV) is a re-emerging mosquito-borne flavivirus that co-opts cellular mechanisms to support viral processes that can reprogram the host transcriptional profile. Such viral-directed transcriptional changes and the pro- or anti-viral outcomes remain understudied. We previously showed that ATF3, a stress-induced transcription factor, is significantly upregulated in ZIKV-infected mammalian cells, along with other cellular and immune response genes. We now define the intracellular pathway responsible for ATF3 activation and elucidate the impact of ATF3 expression on ZIKV infection. We show that during ZIKV infection, the integrated stress response pathway stimulates ATF3 which enhances the innate immune response to antagonize ZIKV infection. This study establishes a link between viral-induced stress response and transcriptional regulation of host defense pathways and thus expands our knowledge of virus-mediated transcriptional mechanisms and transcriptional control of interferon-stimulated genes during ZIKV infection.

Progesterone increases metabolism via the pentose phosphate pathway in bovine uterine epithelial cells

BACKGROUND

During early pregnancy, glucose is essential for the uterine epithelium and the developing embryo. In cows, progesterone increases the secretion of glucose into the uterine lumen. The uterine epithelium can convert glucose to fructose, but other fates of glucose in the uterine epithelium have been sparsely investigated. Therefore, our objective was to investigate how progesterone influences glucose metabolism in immortalized bovine uterine epithelial (BUTE) cells.

METHODS

BUTE cells were grown to 80% confluence and treated with vehicle (DMSO) or 10 µM progesterone for 24 h. Cells were collected and analyzed. Immunohistochemistry was performed on endometrial samples collected from the bovine endometrium on days 1 and 11 of the reproductive cycle.

RESULTS

Progesterone treatment increased glucose consumption of BUTE cells. RNAseq identified 3,072 genes regulated by progesterone. KEGG analysis indicated that progesterone altered genes associated with metabolic pathways and glutathione metabolism. Manually examining genes unique to specific glucose metabolic pathways identified an increase in the rate-limiting enzyme in the pentose phosphate pathway-glucose-6-phosphate dehydrogenase. Functionally, a major product of the pentose phosphate pathway is NADPH, and progesterone treatment increased NADPH levels in BUTE cells. In cows, immunohistochemistry confirmed that glucose-6-phosphate dehydrogenase levels were higher in the uterine epithelium in the luteal phase when progesterone concentrations are high.

CONCLUSIONS

Progesterone increased glucose-6-phosphate dehydrogenase expression and metabolism via the pentose phosphate pathway in the bovine uterine epithelium. This metabolism could provide substrates for cell proliferation, molecules to be secreted into the uterine lumen, or maintain reduction/oxidation balance in the uterine epithelium.

Stabilizing selection and adaptation shape cis and trans gene expression variation in C. elegans

An outstanding question in the evolution of gene expression is the relative influence of neutral processes versus natural selection, including adaptive change driven by directional selection as well as stabilizing selection, which may include compensatory dynamics. These forces shape patterns of gene expression variation within and between species, including the regulatory mechanisms governing expression in cis and trans. In this study, we interrogate intraspecific gene expression variation among seven wild C. elegans strains, with varying degrees of genomic divergence from the reference strain N2, leveraging this system's unique advantages to comprehensively evaluate gene expression evolution. By capturing allele-specific and between-strain changes in expression, we characterize the regulatory architecture and inheritance mode of gene expression variation within C. elegans and assess their relationship to nucleotide diversity, genome evolutionary history, gene essentiality, and other biological factors. We conclude that stabilizing selection is a dominant influence in maintaining expression phenotypes within the species, and the discovery that genes with higher overall expression tend to exhibit fewer expression differences supports this conclusion, as do widespread instances of cis differences compensated in trans. Moreover, analyses of human expression data replicate our finding that higher expression genes have less variable expression. We also observe evidence for directional selection driving expression divergence, and that expression divergence accelerates with increasing genomic divergence. To provide community access to the data from this first analysis of allele-specific expression in C. elegans, we introduce an interactive web application, where users can submit gene-specific queries to view expression, regulatory pattern, inheritance mode, and other information: https://wildworm.biosci.gatech.edu/ase/.

A Plasmodium falciparum MORC protein complex modulates epigenetic control of gene expression through interaction with heterochromatin

Dynamic control of gene expression is critical for blood stage development of malaria parasites. Here, we used multi-omic analyses to investigate transcriptional regulation by the chromatin-associated microrchidia protein, MORC, during asexual blood stage development of the human malaria parasite Plasmodium falciparum. We show that PfMORC (PF3D7_1468100) interacts with a suite of nuclear proteins, including APETALA2 (ApiAP2) transcription factors (PfAP2-G5, PfAP2-O5, PfAP2-I, PF3D7_0420300, PF3D7_0613800, PF3D7_1107800, and PF3D7_1239200), a DNA helicase DS60 (PF3D7_1227100), and other chromatin remodelers (PfCHD1 and PfEELM2). Transcriptomic analysis of PfMORCHA-glmS knockdown parasites revealed 163 differentially expressed genes belonging to hypervariable multigene families, along with upregulation of genes mostly involved in host cell invasion. In vivo genome-wide chromatin occupancy analysis during both trophozoite and schizont stages of development demonstrates that PfMORC is recruited to repressed, multigene families, including the var genes in subtelomeric chromosomal regions. Collectively, we find that PfMORC is found in chromatin complexes that play a role in the epigenetic control of asexual blood stage transcriptional regulation and chromatin organization.

Integrative Multiomics in the Lung Reveals a Protective Role of Asporin in Pulmonary Arterial Hypertension

BACKGROUND

Integrative multiomics can elucidate pulmonary arterial hypertension (PAH) pathobiology, but procuring human PAH lung samples is rare.

METHODS

We leveraged transcriptomic profiling and deep phenotyping of the largest multicenter PAH lung biobank to date (96 disease and 52 control) by integration with clinicopathologic data, genome-wide association studies, Bayesian regulatory networks, single-cell transcriptomics, and pharmacotranscriptomics.

RESULTS

We identified 2 potentially protective gene network modules associated with vascular cells, and we validated ASPN, coding for asporin, as a key hub gene that is upregulated as a compensatory response to counteract PAH. We found that asporin is upregulated in lungs and plasma of multiple independent PAH cohorts and correlates with reduced PAH severity. We show that asporin inhibits proliferation and transforming growth factor-β/phosphorylated SMAD2/3 signaling in pulmonary artery smooth muscle cells from PAH lungs. We demonstrate in Sugen-hypoxia rats that ASPN knockdown exacerbated PAH and recombinant asporin attenuated PAH.

CONCLUSIONS

Our integrative systems biology approach to dissect the PAH lung transcriptome uncovered asporin as a novel protective target with therapeutic potential in PAH.

Systems analysis of miR-199a/b-5p and multiple miR-199a/b-5p targets during chondrogenesis

Changes in chondrocyte gene expression can contribute to the development of osteoarthritis (OA), and so recognition of the regulative processes during chondrogenesis can lead to a better understanding of OA. microRNAs (miRNAs) are key regulators of gene expression in chondrocytes/OA, and we have used a combined experimental, bioinformatic, and systems biology approach to explore the multiple miRNA-mRNA interactions that regulate chondrogenesis. A longitudinal chondrogenesis bioinformatic analysis identified paralogues miR-199a-5p and miR-199b-5p as pro-chondrogenic regulators. Experimental work in human cells demonstrated alteration of miR-199a-5p or miR-199b-5p expression led to significant inverse modulation of key chondrogenic genes and extracellular matrix production. miR-199a/b-5p targets FZD6, ITGA3 and CAV1 were identified by inhibition experiments and verified as direct targets by luciferase assay. The experimental work was used to generate and parameterise a multi-miRNA 14-day chondrogenesis kinetic model to be used as a repository for the experimental work and as a resource for further investigation of this system. This is the first multi-miRNA model of a chondrogenesis-based system, and highlights the complex relationships between regulatory miRNAs, and their target mRNAs.

miR-210 is essential to retinal homeostasis in fruit flies and mice

BACKGROUND

miR-210 is one of the most evolutionarily conserved microRNAs. It is known to be involved in several physiological and pathological processes, including response to hypoxia, angiogenesis, cardiovascular diseases and cancer. Recently, new roles of this microRNA are emerging in the context of eye and visual system homeostasis. Recent studies in Drosophila melanogaster unveiled that the absence of miR-210 leads to a progressive retinal degeneration characterized by the accumulation of lipid droplets and disruptions in lipid metabolism. However, the possible conservation of miR-210 knock-out effect in the mammalian retina has yet to be explored.

RESULTS

We further investigated lipid anabolism and catabolism in miR-210 knock-out (KO) flies, uncovering significant alterations in gene expression within these pathways. Additionally, we characterized the retinal morphology of flies overexpressing (OE) miR-210, which was not affected by the increased levels of the microRNA. For the first time, we also characterized the retinal morphology of miR-210 KO and OE mice. Similar to flies, miR-210 OE did not affect retinal homeostasis, whereas miR-210 KO mice exhibited photoreceptor degeneration. To explore other potential parallels between miR-210 KO models in flies and mice, we examined lipid metabolism, circadian behaviour, and retinal transcriptome in mice, but found no similarities. Specifically, RNA-seq confirmed the lack of involvement of lipid metabolism in the mice's pathological phenotype, revealing that the differentially expressed genes were predominantly associated with chloride channel activity and extracellular matrix homeostasis. Simultaneously, transcriptome analysis of miR-210 KO fly brains indicated that the observed alterations extend beyond the eye and may be linked to neuronal deficiencies in signal detection and transduction.

CONCLUSIONS

We provide the first morphological characterization of the retina of miR-210 KO and OE mice, investigating the role of this microRNA in mammalian retinal physiology and exploring potential parallels with phenotypes observed in fly models. Although the lack of similarities in lipid metabolism, circadian behaviour, and retinal transcriptome in mice suggests divergent mechanisms of retinal degeneration between the two species, transcriptome analysis of miR-210 KO fly brains indicates the potential existence of a shared upstream mechanism contributing to retinal degeneration in both flies and mammals.

Functional screen for mediators of onco-mRNA translation specificity

Oncogenic protein dosage is tightly regulated to enable cancer cells to adapt and survive. Whether this is regulated at the level of translational control and the key factors in cis and trans remain unknown. The Myc oncogene is a central paradigm of an exquisitely regulated oncogene and a major driver of pancreatic ductal adenocarcinoma (PDAC). Using a functional genome-wide CRISPRi screen in PDAC cells, we identified activators of selective MYC translation through its 5' untranslated region (5'UTR) and validated four RNA binding proteins (RBPs), including epitranscriptome modifiers. Among these RBPs, our top hit was RBM42, which is highly expressed in PDAC and predicts poor survival. Combining polysome sequencing and CLIP-seq analyses, we find that RBM42 binds and selectively regulates the translation of MYC and a precise, yet vital suite of pro-oncogenic transcripts, including JUN and EGFR . Mechanistically, employing IP-mass spectrometry analysis, we find that RMB42 is a novel ribosome-associated protein (RAP). Using DMS-Seq and mutagenesis analysis, we show that RBM42 directly binds and remodels the MYC 5'UTR RNA structure, facilitating the formation of the translation pre-initiation complex. Importantly, RBM42 is necessary for human PDAC cell growth and fitness and PDAC tumorigenesis in xenograft mouse models in a Myc-dependent manner in vivo . In PDAC patient samples, RBM42 expression is correlated with Myc protein levels and transcriptional activity. This work transforms our understanding of the translational code in cancer and offers a new therapeutic opening to target the expression of oncogenes.

Association of RDoC dimensions with post-mortem brain transcriptional profiles in Alzheimer's Disease

INTRODUCTION

Neuropsychiatric symptoms are common in people with Alzheimer's disease (AD) across all severity stages. Their heterogeneous presentation and variable temporal association with cognitive decline suggest shared and distinct biological mechanisms. We hypothesized that specific patterns of gene expression associate with distinct NIMH Research Domain Criteria (RDoC) domains in AD.

METHODS

Post-mortem bulk RNAseq on the insula and anterior cingulate cortex from 60 brain donors representing the spectrum of canonical AD neuropathology combined with natural language processing approaches based on the RDoC Clinical Domains.

RESULTS

Distinct sets of >100 genes (p FDR <0.05) were specifically associated with at least one clinical domain (Cognitive, Social, Negative, Positive, Arousal). In addition, dysregulation of immune response pathways was shared across domains and brain regions.

DISCUSSION

Our findings provide evidence for distinct transcriptional profiles associated with RDoC domains suggesting that each dimension is characterized by specific sets of genes providing insight into the underlying mechanisms.

Analysis of Differential Gene Expression under Acute Lead or Mercury Exposure in Larval Zebrafish Using RNA-Seq

This study was first conducted to investigate the effects of acute lead exposure on developing zebrafish embryos or larvae from 24 to 120 h post-fertilization (hpf). Our data showed that treatment with 50-200 μM lead significantly affected larval survivability and morphology compared to the respective control. Second, we chose 120 hpf larvae treated with 12.5 μM lead for RNA sequencing due to its exposure level being sufficient to produce toxic effects with minimum death and lead bioaccumulation in developing zebrafish. A total of 137.45 million raw reads were obtained, and more than 86% of clean data were mapped to the zebrafish reference genome. Differential expression profiles generated 116 up- and 34 down-regulated genes upon lead exposure. The most enriched GO terms for representative DEGs were ion transport and lipid metabolism. Third, a comparison with the dataset of mercury-regulated gene expression identified 94 genes (64 up-regulated and 30 down-regulated) for exposure specific to lead, as well as 422 genes (338 up-regulated and 84 down-regulated) for exposure specific to mercury. In addition, 56 genes were co-regulated by micromolar mercury and lead treatment, and the expression of thirteen genes, including mt2, ctssb.1, prdx1, txn, sqrdl, tmprss13a, socs3a, trpv6, abcb6a, gsr, hbz, fads2, and zgc:92590 were validated by qRT-PCR. These genes were mainly associated with metal ion binding, proteolysis, antioxidant activity, signal transduction, calcium ion or oxygen transport, the fatty acid biosynthetic process, and protein metabolism. Taken together, these findings help better understand the genome-wide responses of developing zebrafish to lead or mercury and provide potential biomarkers for acute exposure to toxic metals.

Butterfly eggs prime anti-herbivore defense in an annual but not perennial Arabidopsis species

MAIN CONCLUSION

Unlike Arabidopsis thaliana, defenses of Arabidopsis lyrata against Pieris brassicae larval feeding are not primable by P. brassicae eggs. Thus, egg primability of plant anti-herbivore defenses is not phylogenetically conserved in the genus Arabidopsis. While plant anti-herbivore defenses of the annual species Arabidopsis thaliana were shown to be primable by Pieris brassicae eggs, the primability of the phylogenetically closely related perennial Arabidopsis lyrata has not yet been investigated. Previous studies revealed that closely related wild Brassicaceae plant species, the annual Brassica nigra and the perennial B. oleracea, exhibit an egg-primable defense trait, even though they have different life spans. Here, we tested whether P. brassicae eggs prime anti-herbivore defenses of the perennial A. lyrata. We exposed A. lyrata to P. brassicae eggs and larval feeding and assessed their primability by (i) determining the biomass of P. brassicae larvae after feeding on plants with and without prior P. brassicae egg deposition and (ii) investigating the plant transcriptomic response after egg deposition and/or larval feeding. For comparison, these studies were also conducted with A. thaliana. Consistent with previous findings, A. thaliana's response to prior P. brassicae egg deposition negatively affected conspecific larvae feeding upon A. thaliana. However, this was not observed in A. lyrata. Arabidopsis thaliana responded to P. brassicae eggs with strong transcriptional reprogramming, whereas A. lyrata responses to eggs were negligible. In response to larval feeding, A. lyrata exhibited a greater transcriptome change compared to A. thaliana. Among the strongly feeding-induced A. lyrata genes were those that are egg-primed in feeding-induced A. thaliana, i.e., CAX3, PR1, PR5, and PDF1.4. These results suggest that A. lyrata has evolved a robust feeding response that is independent from prior egg exposure.

Comparative transcriptomic analysis of articular cartilage of post-traumatic osteoarthritis models

Animal models of post-traumatic osteoarthritis (PTOA) recapitulate the pathological changes observed in human PTOA. Here, skeletally mature C57Bl6 mice were subjected to either rapid-onset non-surgical mechanical rupture of the anterior cruciate ligament (ACL) or to surgical destabilisation of the medial meniscus (DMM). Transcriptome profiling of micro-dissected cartilage at day 7 or day 42 following ACL or DMM procedure, respectively, showed that the two models were comparable and highly correlative. Gene ontology (GO) enrichment analysis identified similarly enriched pathways that were overrepresented by anabolic terms. To address the transcriptome changes more completely in the ACL model, we also performed small RNA sequencing, describing the first microRNA profile of this model. miR-199-5p was amongst the most abundant, yet differentially expressed, microRNAs, and its inhibition in primary human chondrocytes led to a transcriptome response that was comparable to that observed in both human 'OA damaged vs intact cartilage' and murine DMM cartilage datasets. We also experimentally verified CELSR1, GIT1, ECE1 and SOS2 as novel miR-199-5p targets. Together, these data support the use of the ACL rupture model as a non-invasive companion to the DMM model.

LsRTDv1, a reference transcript dataset for accurate transcript-specific expression analysis in lettuce

Accurate quantification of gene and transcript-specific expression, with the underlying knowledge of precise transcript isoforms, is crucial to understanding many biological processes. Analysis of RNA sequencing data has benefited from the development of alignment-free algorithms which enhance the precision and speed of expression analysis. However, such algorithms require a reference transcriptome. Here we generate a reference transcript dataset (LsRTDv1) for lettuce (cv. Saladin), combining long- and short-read sequencing with publicly available transcriptome annotations, and filtering to keep only transcripts with high-confidence splice junctions and transcriptional start and end sites. LsRTDv1 identifies novel genes (mostly long non-coding RNAs) and increases the number of transcript isoforms per gene in the lettuce genome from 1.4 to 2.7. We show that LsRTDv1 significantly increases the mapping rate of RNA-seq data from a lettuce time-series experiment (mock- and Botrytis cinerea-inoculated) and enables detection of genes that are differentially alternatively spliced in response to infection as well as transcript-specific expression changes. LsRTDv1 is a valuable resource for investigation of transcriptional and alternative splicing regulation in lettuce.

Transcriptional and methylation outcomes of didehydro-cortistatin A use in HIV-1-infected CD4+ T cells

Ongoing viral transcription from the reservoir of HIV-1 infected long-lived memory CD4+ T cells presents a barrier to cure and associates with poorer health outcomes for people living with HIV, including chronic immune activation and inflammation. We previously reported that didehydro-cortistatin A (dCA), an HIV-1 Tat inhibitor, blocks HIV-1 transcription. Here, we examine the impact of dCA on host immune CD4+ T-cell transcriptional and epigenetic states. We performed a comprehensive analysis of genome-wide transcriptomic and DNA methylation profiles upon long-term dCA treatment of primary human memory CD4+ T cells. dCA prompted specific transcriptional and DNA methylation changes in cell cycle, histone, interferon-response, and T-cell lineage transcription factor genes, through inhibition of both HIV-1 and Mediator kinases. These alterations establish a tolerogenic Treg/Th2 phenotype, reducing viral gene expression and mitigating inflammation in primary CD4+ T cells during HIV-1 infection. In addition, dCA suppresses the expression of lineage-defining transcription factors for Th17 and Th1 cells, critical HIV-1 targets, and reservoirs. dCA's benefits thus extend beyond viral transcription inhibition, modulating the immune cell landscape to limit HIV-1 acquisition and inflammatory environment linked to HIV infection.

Smoking-informed methylation and expression QTLs in human brain and colocalization with smoking-associated genetic loci

Smoking is a leading cause of preventable morbidity and mortality. Smoking is heritable, and genome-wide association studies (GWASs) of smoking behaviors have identified hundreds of significant loci. Most GWAS-identified variants are noncoding with unknown neurobiological effects. We used genome-wide genotype, DNA methylation, and RNA sequencing data in postmortem human nucleus accumbens (NAc) to identify cis-methylation/expression quantitative trait loci (meQTLs/eQTLs), investigate variant-by-cigarette smoking interactions across the genome, and overlay QTL evidence at smoking GWAS-identified loci to evaluate their regulatory potential. Active smokers (N = 52) and nonsmokers (N = 171) were defined based on cotinine biomarker levels and next-of-kin reporting. We simultaneously tested variant and variant-by-smoking interaction effects on methylation and expression, separately, adjusting for biological and technical covariates and correcting for multiple testing using a two-stage procedure. We found >2 million significant meQTL variants (padj < 0.05) corresponding to 41,695 unique CpGs. Results were largely driven by main effects, and five meQTLs, mapping to NUDT12, FAM53B, RNF39, and ADRA1B, showed a significant interaction with smoking. We found 57,683 significant eQTL variants for 958 unique eGenes (padj < 0.05) and no smoking interactions. Colocalization analyses identified loci with smoking-associated GWAS variants that overlapped meQTLs/eQTLs, suggesting that these heritable factors may influence smoking behaviors through functional effects on methylation/expression. One locus containing MUSTN1 and ITIH4 colocalized across all data types (GWAS, meQTL, and eQTL). In this first genome-wide meQTL map in the human NAc, the enriched overlap with smoking GWAS-identified genetic loci provides evidence that gene regulation in the brain helps explain the neurobiology of smoking behaviors.

ETENLNC: An end to end lncRNA identification and analysis framework to facilitate construction of known and novel lncRNA regulatory networks

Long non-coding RNAs (lncRNAs) play crucial roles in the regulation of gene expression and maintenance of genomic integrity through various interactions with DNA, RNA, and proteins. The availability of large-scale sequence data from various high-throughput platforms has opened possibilities to identify, predict, and functionally annotate lncRNAs. As a result, there is a growing demand for an integrative computational framework capable of identifying known lncRNAs, predicting novel lncRNAs, and inferring the downstream regulatory interactions of lncRNAs at the genome-scale. We present ETENLNC (End-To-End-Novel-Long-NonCoding), a user-friendly, integrative, open-source, scalable, and modular computational framework for identifying and analyzing lncRNAs from raw RNA-Seq data. ETENLNC employs six stringent filtration steps to identify novel lncRNAs, performs differential expression analysis of mRNA and lncRNA transcripts, and predicts regulatory interactions between lncRNAs, mRNAs, miRNAs, and proteins. We benchmarked ETENLNC against six existing tools and optimized it for desktop workstations and high-performance computing environments using data from three different species. ETENLNC is freely available on GitHub: https://github.com/EvolOMICS-TU/ETENLNC.

DifferentialRegulation: a Bayesian hierarchical approach to identify differentially regulated genes

Although transcriptomics data is typically used to analyze mature spliced mRNA, recent attention has focused on jointly investigating spliced and unspliced (or precursor-) mRNA, which can be used to study gene regulation and changes in gene expression production. Nonetheless, most methods for spliced/unspliced inference (such as RNA velocity tools) focus on individual samples, and rarely allow comparisons between groups of samples (e.g. healthy vs. diseased). Furthermore, this kind of inference is challenging, because spliced and unspliced mRNA abundance is characterized by a high degree of quantification uncertainty, due to the prevalence of multi-mapping reads, ie reads compatible with multiple transcripts (or genes), and/or with both their spliced and unspliced versions. Here, we present DifferentialRegulation, a Bayesian hierarchical method to discover changes between experimental conditions with respect to the relative abundance of unspliced mRNA (over the total mRNA). We model the quantification uncertainty via a latent variable approach, where reads are allocated to their gene/transcript of origin, and to the respective splice version. We designed several benchmarks where our approach shows good performance, in terms of sensitivity and error control, vs. state-of-the-art competitors. Importantly, our tool is flexible, and works with both bulk and single-cell RNA-sequencing data. DifferentialRegulation is distributed as a Bioconductor R package.

Spatial transcriptomic validation of a biomimetic model of fibrosis enables re-evaluation of a therapeutic antibody targeting LOXL2

Matrix stiffening by lysyl oxidase-like 2 (LOXL2)-mediated collagen cross-linking is proposed as a core feedforward mechanism that promotes fibrogenesis. Failure in clinical trials of simtuzumab (the humanized version of AB0023, a monoclonal antibody against human LOXL2) suggested that targeting LOXL2 may not have disease relevance; however, target engagement was not directly evaluated. We compare the spatial transcriptome of active human lung fibrogenesis sites with different human cell culture models to identify a disease-relevant model. Within the selected model, we then evaluate AB0023, identifying that it does not inhibit collagen cross-linking or reduce tissue stiffness, nor does it inhibit LOXL2 catalytic activity. In contrast, it does potently inhibit angiogenesis consistent with an alternative, non-enzymatic mechanism of action. Thus, AB0023 is anti-angiogenic but does not inhibit LOXL2 catalytic activity, collagen cross-linking, or tissue stiffening. These findings have implications for the interpretation of the lack of efficacy of simtuzumab in clinical trials of fibrotic diseases.

Concurrent analysis of genome and transcriptome in one single cell

Thus far, multiple techniques for single cell analysis have been developed, yet we lack a relatively simple tool to assess DNA and RNA from the same cell at whole-transcriptome and whole-genome depths. Here we present an updated method for physical separation of cytoplasmic RNA from the nuclei, which allows for simultaneous studies of DNA and RNA from the same single cell. The method consists of three steps-(1) immobilization of a single cell on solid substrate, (2) hypotonic lysis of immobilized single cell, and (3) separation of cytosol containing aqueous phase and immobilized nucleus. We found that DNA and RNA extracted from single cell using our approach is suitable for downstream sequencing-based applications. We demonstrated that the coverage of transcriptome and genome sequencing data obtained after DNA/RNA separation is similar to that observed without separation. We also showed that the separation procedure does not create any noticeable bias in observed mutational load or mutation spectra. Thus, our method can serve as a tool for simultaneous complex analysis of the genome and transcriptome, providing necessary information on the relationship between somatic mutations and the regulation of gene expression.

Dysregulated FOXO1 activity drives skeletal muscle intrinsic dysfunction in amyotrophic lateral sclerosis

Amyotrophic Lateral Sclerosis (ALS) is a multisystemic neurodegenerative disorder, with accumulating evidence indicating metabolic disruptions in the skeletal muscle preceding disease symptoms, rather than them manifesting as a secondary consequence of motor neuron (MN) degeneration. Hence, energy homeostasis is deeply implicated in the complex physiopathology of ALS and skeletal muscle has emerged as a key therapeutic target. Here, we describe intrinsic abnormalities in ALS skeletal muscle, both in patient-derived muscle cells and in muscle cell lines with genetic knockdown of genes related to familial ALS, such as TARDBP (TDP-43) and FUS. We found a functional impairment of myogenesis that parallels defects of glucose oxidation in ALS muscle cells. We identified FOXO1 transcription factor as a key mediator of these metabolic and functional features in ALS muscle, via gene expression profiling and biochemical surveys in TDP-43 and FUS-silenced muscle progenitors. Strikingly, inhibition of FOXO1 mitigated the impaired myogenesis in both the genetically modified and the primary ALS myoblasts. In addition, specific in vivo conditional knockdown of TDP-43 or FUS orthologs (TBPH or caz) in Drosophila muscle precursor cells resulted in decreased innervation and profound dysfunction of motor nerve terminals and neuromuscular synapses, accompanied by motor abnormalities and reduced lifespan. Remarkably, these phenotypes were partially corrected by foxo inhibition, bolstering the potential pharmacological management of muscle intrinsic abnormalities associated with ALS. The findings demonstrate an intrinsic muscle dysfunction in ALS, which can be modulated by targeting FOXO factors, paving the way for novel therapeutic approaches that focus on the skeletal muscle as complementary target tissue.

PDGFRα signaling regulates Srsf3 transcript binding to affect PI3K signaling and endosomal trafficking

Signaling through the platelet-derived growth factor receptor alpha (PDGFRα) plays a critical role in craniofacial development, as mutations in PDGFRA are associated with cleft lip/palate in humans and Pdgfra mutant mouse models display varying degrees of facial clefting. Phosphatidylinositol 3-kinase (PI3K)/Akt is the primary effector of PDGFRα signaling during skeletal development in the mouse. We previously demonstrated that Akt phosphorylates the RNA-binding protein serine/arginine-rich splicing factor 3 (Srsf3) downstream of PI3K-mediated PDGFRα signaling in mouse embryonic palatal mesenchyme (MEPM) cells, leading to its nuclear translocation. We further showed that ablation of Srsf3 in the murine neural crest lineage results in severe midline facial clefting, due to defects in proliferation and survival of cranial neural crest cells, and widespread alternative RNA splicing (AS) changes. Here, we sought to determine the molecular mechanisms by which Srsf3 activity is regulated downstream of PDGFRα signaling to control AS of transcripts necessary for craniofacial development. We demonstrated via enhanced UV-crosslinking and immunoprecipitation (eCLIP) of MEPM cells that PDGF-AA stimulation leads to preferential binding of Srsf3 to exons and loss of binding to canonical Srsf3 CA-rich motifs. Through the analysis of complementary RNA-seq data, we showed that Srsf3 activity results in the preferential inclusion of exons with increased GC content and lower intron to exon length ratio. Moreover, we found that the subset of transcripts that are bound by Srsf3 and undergo AS upon PDGFRα signaling commonly encode regulators of PI3K signaling and early endosomal trafficking. Functional validation studies further confirmed that Srsf3 activity downstream of PDGFRα signaling leads to retention of the receptor in early endosomes and increases in downstream PI3K-mediated Akt signaling. Taken together, our findings reveal that growth factor-mediated phosphorylation of an RNA-binding protein underlies gene expression regulation necessary for mammalian craniofacial development.

Transcriptomic profiles in major depressive disorder: the role of immunometabolic and cell-cycle-related pathways in depression with different levels of inflammation

Transcriptomic profiles are important indicators for molecular mechanisms and pathways involved in major depressive disorder (MDD) and its different phenotypes, such as immunometabolic depression. We performed whole-transcriptome and pathway analyses on 139 individuals from the observational, case-control, BIOmarkers in DEPression (BIODEP) study, 105 with MDD and 34 controls. We divided MDD participants based on levels of inflammation, as measured by serum high-sensitivity C-reactive protein (CRP), in n = 39 'not inflamed' (CRP < 1 mg/L), n = 31 with 'elevated CRP' (1-3 mg/L), and n = 35 with 'low-grade inflammation' (>3 mg/L). We performed whole-blood RNA sequencing using Illumina NextSeq 550 and statistical analyses with the Deseq2 package for R statistics (RUV-corrected) and subsequent pathway analyses with Ingenuity Pathway Analysis. Immunometabolic pathways were activated in individuals with CRP > 1 mg/L, although surprisingly the CRP 1-3 group showed stronger immune activation than the CRP > 3 group. The main pathways identified in the comparison between CRP < 1 group and controls were cell-cycle-related, which may be protective against immunometabolic abnormalities in this 'non-inflamed' depressed group. We further divided MDD participants based on exposure and response to antidepressants (n = 47 non-responders, n = 37 responders, and n = 22 unmedicated), and identified specific immunomodulatory and neuroprotective pathways in responders (especially vs. non-responders), which could be relevant to treatment response. In further subgroup analyses, we found that the specific transcriptional profile of responders is independent of CRP levels, and that the inhibition of cell-cycle-related pathways in MDD with CRP < 1 mg/L is present only in those who are currently depressed, and not in the responders. The present study demonstrates immunometabolic and cell-cycle-related transcriptomic pathways associated with MDD and different (CRP-based and treatment-based) MDD phenotypes, while shedding light on potential molecular mechanisms that could prevent or facilitate an individual's trajectory toward immunometabolic depression and/or treatment-non-responsive depression. The recognition and integration of these mechanisms will facilitate a precision-medicine approach in MDD.

Transposable elements alter gene expression and may impact response to cisplatin therapy in ovarian cancer

Cisplatin is widely employed for cancer treatment; therefore, understanding resistance to this drug is critical for therapeutic practice. While studies have delved into differential gene expression in the context of cisplatin resistance, findings remain somewhat scant. We performed a comprehensive investigation of transposable elements (TEs) expression and their impact in host genes in two cisplatin-treated ovarian cancer cell lines. RNA-seq, ATAC-seq, and in-depth bioinformatics analysis were used to compare cisplatin-sensitive and -resistant ovarian cancer cell lines. Our results reveal that cisplatin therapy alters not only the expression of protein-coding genes, but also key TEs, including LINE1, Alu, and endogenous retroviruses, in both cisplatin-sensitive and -resistant cell lines. By co-expressing with downstream genes or by creating chimeric transcripts with host genes at their insertion sites, these TEs seem to control the expression of protein-coding genes, including tumor-related genes. Our model uncovers TEs influencing the expression of cancer genes and cancer pathways. Collectively, our findings indicate that TE alterations associated with cisplatin treatment occur in critical cancer genes and cellular pathways synergically. This research highlights the importance of considering the entire spectrum of transcribed elements in the genome, especially TE expression, for a complete understanding of complex models like cancer response to treatment.

Stable and robust Xi and Y transcriptomes drive cell-type-specific autosomal and Xa responses in vivo and in vitro in four human cell types

Recent in vitro studies of human sex chromosome aneuploidy showed that the Xi ("inactive" X) and Y chromosomes broadly modulate autosomal and Xa ("active" X) gene expression. We tested these findings in vivo. Linear modeling of CD4+ T cells and monocytes from individuals with one to three X chromosomes and zero to two Y chromosomes revealed 82 sex-chromosomal and 344 autosomal genes whose expression changed significantly with Xi and/or Y dosage in vivo. Changes in sex-chromosomal expression were remarkably constant in vivo and in vitro; autosomal responses to Xi and/or Y dosage were largely cell-type specific (∼2.6-fold more variation than sex-chromosomal responses). Targets of the sex-chromosomal transcription factors ZFX and ZFY accounted for a significant fraction of these autosomal responses both in vivo and in vitro. We conclude that the human Xi and Y transcriptomes are surprisingly robust and stable, yet they modulate autosomal and Xa genes in a cell-type-specific fashion.

Identification of diagnostic candidates in Mendelian disorders using an RNA sequencing-centric approach

BACKGROUND

RNA sequencing (RNA-seq) is increasingly being used as a complementary tool to DNA sequencing in diagnostics where DNA analysis has been uninformative. RNA-seq enables the identification of aberrant splicing and aberrant gene expression, improving the interpretation of variants of unknown significance (VUSs), and provides the opportunity to scan the transcriptome for aberrant splicing and expression in relevant genes that may be the cause of a patient's phenotype. This work aims to investigate the feasibility of generating new diagnostic candidates in patients without a previously reported VUS using an RNA-seq-centric approach.

METHODS

We systematically assessed the transcriptomic profiles of 86 patients with suspected Mendelian disorders, 38 of whom had no candidate sequence variant, using RNA from blood samples. Each VUS was visually inspected to search for splicing abnormalities. Once aberrant splicing was identified in cases with VUS, multiple open-source alternative splicing tools were used to investigate if they would identify what was observed in IGV. Expression outliers were detected using OUTRIDER. Diagnoses in cases without a VUS were explored using two separate strategies.

RESULTS

RNA-seq allowed us to assess 71% of VUSs, detecting aberrant splicing in 14/48 patients with a VUS. We identified four new diagnoses by detecting novel aberrant splicing events in patients with no candidate sequence variants from prior DNA testing (n = 32) or where the candidate VUS did not affect splicing (n = 23). An additional diagnosis was made through the detection of skewed X-inactivation.

CONCLUSION

This work demonstrates the utility of an RNA-centric approach in identifying novel diagnoses in patients without candidate VUSs. It underscores the utility of blood-based RNA analysis in improving diagnostic yields and highlights optimal approaches for such analyses.

Evaluating metagenomics and targeted approaches for diagnosis and surveillance of viruses

BACKGROUND

Metagenomics is a powerful approach for the detection of unknown and novel pathogens. Workflows based on Illumina short-read sequencing are becoming established in diagnostic laboratories. However, high sequencing depth requirements, long turnaround times, and limited sensitivity hinder broader adoption. We investigated whether we could overcome these limitations using protocols based on untargeted sequencing with Oxford Nanopore Technologies (ONT), which offers real-time data acquisition and analysis, or a targeted panel approach, which allows the selective sequencing of known pathogens and could improve sensitivity.

METHODS

We evaluated detection of viruses with readily available untargeted metagenomic workflows using Illumina and ONT, and an Illumina-based enrichment approach using the Twist Bioscience Comprehensive Viral Research Panel (CVRP), which targets 3153 viruses. We tested samples consisting of a dilution series of a six-virus mock community in a human DNA/RNA background, designed to resemble clinical specimens with low microbial abundance and high host content. Protocols were designed to retain the host transcriptome, since this could help confirm the absence of infectious agents. We further compared the performance of commonly used taxonomic classifiers.

RESULTS

Capture with the Twist CVRP increased sensitivity by at least 10-100-fold over untargeted sequencing, making it suitable for the detection of low viral loads (60 genome copies per ml (gc/ml)), but additional methods may be needed in a diagnostic setting to detect untargeted organisms. While untargeted ONT had good sensitivity at high viral loads (60,000 gc/ml), at lower viral loads (600-6000 gc/ml), longer and more costly sequencing runs would be required to achieve sensitivities comparable to the untargeted Illumina protocol. Untargeted ONT provided better specificity than untargeted Illumina sequencing. However, the application of robust thresholds standardized results between taxonomic classifiers. Host gene expression analysis is optimal with untargeted Illumina sequencing but possible with both the CVRP and ONT.

CONCLUSIONS

Metagenomics has the potential to become standard-of-care in diagnostics and is a powerful tool for the discovery of emerging pathogens. Untargeted Illumina and ONT metagenomics and capture with the Twist CVRP have different advantages with respect to sensitivity, specificity, turnaround time and cost, and the optimal method will depend on the clinical context.

ATM phosphorylation of CD98HC increases antiporter membrane localization and prevents chronic toxic glutamate accumulation in Ataxia telangiectasia

Ataxia telangiectasia (A-T) is a rare genetic disorder characterized by neurological defects, immunodeficiency, cancer predisposition, radiosensitivity, decreased blood vessel integrity, and diabetes. ATM, the protein mutated in A-T, responds to DNA damage and oxidative stress, but its functional relationship to the progressive clinical manifestation of A-T is not understood. CD98HC chaperones cystine/glutamate (x c - ) and cationic/neutral amino acid (y + L) antiporters to the cell membrane, and CD98HC phosphorylation by ATM accelerates membrane localization to acutely increase amino acid transport. Loss of ATM impacts tissues reliant on SLC family antiporters relevant to A-T phenotypes, such as endothelial cells (telangiectasia) and pancreatic α-cells (fatty liver and diabetes) with toxic glutamate accumulation. Bypassing the antiporters restores intracellular metabolic balance both in ATM-deficient cells and mouse models. These findings provide new insight into the long-known benefits of N-acetyl cysteine to A-T cells beyond oxidative stress through removing excess glutamate by production of glutathione.

ER procollagen storage defect without coupled unfolded protein response drives precocious arthritis

Collagenopathies are a group of clinically diverse disorders caused by defects in collagen folding and secretion. For example, mutations in the gene encoding collagen type-II, the primary collagen in cartilage, can lead to diverse chondrodysplasias. One example is the Gly1170Ser substitution in procollagen-II, which causes precocious osteoarthritis. Here, we biochemically and mechanistically characterize an induced pluripotent stem cell-based cartilage model of this disease, including both hetero- and homozygous genotypes. We show that Gly1170Ser procollagen-II is notably slow to fold and secrete. Instead, procollagen-II accumulates intracellularly, consistent with an endoplasmic reticulum (ER) storage disorder. Likely owing to the unique features of the collagen triple helix, this accumulation is not recognized by the unfolded protein response. Gly1170Ser procollagen-II interacts to a greater extent than wild-type with specific ER proteostasis network components, consistent with its slow folding. These findings provide mechanistic elucidation into the etiology of this disease. Moreover, the easily expandable cartilage model will enable rapid testing of therapeutic strategies to restore proteostasis in the collagenopathies.

Lipid droplets and small extracellular vesicles: More than two independent entities

Despite increasing knowledge about small extracellular vesicle (sEV) composition and functions in cell-cell communication, the mechanism behind their biogenesis remains unclear. Here, we reveal for the first time that sEV biogenesis and release into the microenvironment are tightly connected with another important organelle, Lipid Droplets (LDs). The correlation was observed in several human cancer cell lines as well as patient-derived colorectal cancer stem cells (CR-CSCs). Our results demonstrated that external stimuli such as radiation, pH, hypoxia or lipid-interfering drugs, known to affect the number of LDs/cell, similarly influenced sEV secretion. Importantly, through multiple omics data, at both mRNA and protein levels, we revealed RAB5C as a potential important molecular player behind this organelle connection. Altogether, the potential to fine-tune sEV biogenesis by targeting LDs could significantly impact the amount, cargos and properties of these sEVs, opening new clinical perspectives.

Single-cell transcriptomic analysis of retinal immune regulation and blood-retinal barrier function during experimental autoimmune uveitis

Uveitis is characterised by breakdown of the blood-retinal barrier (BRB), allowing infiltration of immune cells that mediate intraocular inflammation, which can lead to irreversible damage of the neuroretina and the loss of sight. Treatment of uveitis relies heavily on corticosteroids and systemic immunosuppression due to limited understanding of disease pathogenesis. We performed single-cell RNA-sequencing of retinas, as well as bulk RNA-sequencing of retinal pigment epithelial (RPE) cells from mice with experimental autoimmune uveitis (EAU) versus healthy control. This revealed that the Th1/Th17-driven disease induced strong gene expression changes in response to inflammation in rods, cones, Müller glia and RPE. In particular, Müller glia and RPE cells were found to upregulate expression of chemokines, complement factors, leukocyte adhesion molecules and MHC class II, thus highlighting their contributions to immune cell recruitment and antigen presentation at the inner and outer BRB, respectively. Additionally, ligand-receptor interaction analysis with CellPhoneDB revealed key interactions between Müller glia and T cell / natural killer cell subsets via chemokines, galectin-9 to P4HB/TIM-3, PD-L1 to PD-1, and nectin-2/3 to TIGIT signalling axes. Our findings elucidate mechanisms contributing to breakdown of retinal immune privilege during uveitis and identify novel targets for therapeutic interventions.

Robustness and variability in Caenorhabditis elegans dauer gene expression

Both plasticity and robustness are pervasive features of developmental programs. The dauer in Caenorhabditis elegans is an arrested, hypometabolic alternative to the third larval stage of the nematode. Dauers undergo dramatic tissue remodeling and extensive physiological, metabolic, behavioral, and gene expression changes compared to conspecifics that continue development and can be induced by several adverse environments or genetic mutations that act as independent and parallel inputs into the larval developmental program. Therefore, dauer induction is an example of phenotypic plasticity. However, whether gene expression in dauer larvae induced to arrest development by different genetic or environmental triggers is invariant or varies depending on their route into dauer has not been examined. By using RNA-sequencing to characterize gene expression in different types of dauer larvae and computing the variance and concordance within Gene Ontologies (GO) and gene expression networks, we find that the expression patterns within most pathways are strongly correlated between dauer larvae, suggestive of transcriptional robustness. However, gene expression within specific defense pathways, pathways regulating some morphological traits, and several metabolic pathways differ between the dauer larvae. We speculate that the transcriptional robustness of core dauer pathways allows for the buffering of variation in the expression of genes involved in adaptation, allowing the dauers induced by different stimuli to survive in and exploit different niches.

Non-canonical Wnt signaling triggered by WNT2B drives adrenal aldosterone production

The steroid hormone aldosterone, produced by the zona glomerulosa (zG) of the adrenal gland, is a master regulator of plasma electrolytes and blood pressure. While aldosterone control by the renin-angiotensin system is well understood, other key regulatory factors have remained elusive. Here, we replicated a prior association between a non-coding variant in WNT2B and an increased risk of primary aldosteronism, a prevalent and debilitating disease caused by excessive aldosterone production. We further show that in both mice and humans, WNT2B is expressed in the mesenchymal capsule surrounding the adrenal cortex, in close proximity to the zG. Global loss of Wnt2b in the mouse results in a dysmorphic and hypocellular zG, with impaired aldosterone production. Similarly, humans harboring WNT2B loss-of-function mutations develop a novel form of Familial Hyperreninemic Hypoaldosteronism, designated here as Type 4. Additionally, we demonstrate that WNT2B signals by activating the non-canonical Wnt/planar cell polarity pathway. Our findings identify WNT2B as a key regulator of zG function and aldosterone production with important clinical implications.

The role of heterochronic gene expression and regulatory architecture in early developmental divergence

New developmental programs can evolve through adaptive changes to gene expression. The annelid Streblospio benedicti has a developmental dimorphism, which provides a unique intraspecific framework for understanding the earliest genetic changes that take place during developmental divergence. Using comparative RNAseq through ontogeny, we find that only a small proportion of genes are differentially expressed at any time, despite major differences in larval development and life history. These genes shift expression profiles across morphs by either turning off any expression in one morph or changing the timing or amount of gene expression. We directly connect the contributions of these mechanisms to differences in developmental processes. We examine F1 offspring - using reciprocal crosses - to determine maternal mRNA inheritance and the regulatory architecture of gene expression. These results highlight the importance of both novel gene expression and heterochronic shifts in developmental evolution, as well as the trans-acting regulatory factors in initiating divergence.

Nitrate and ammonium, the yin and yang of nitrogen uptake: a time-course transcriptomic study in rice

Nitrogen is an essential nutrient for plants and a major determinant of plant growth and crop yield. Plants acquire nitrogen mainly in the form of nitrate and ammonium. Both nitrogen sources affect plant responses and signaling pathways in a different way, but these signaling pathways interact, complicating the study of nitrogen responses. Extensive transcriptome analyses and the construction of gene regulatory networks, mainly in response to nitrate, have significantly advanced our understanding of nitrogen signaling and responses in model plants and crops. In this study, we aimed to generate a more comprehensive gene regulatory network for the major crop, rice, by incorporating the interactions between ammonium and nitrate. To achieve this, we assessed transcriptome changes in rice roots and shoots over an extensive time course under single or combined applications of the two nitrogen sources. This dataset enabled us to construct a holistic co-expression network and identify potential key regulators of nitrogen responses. Next to known transcription factors, we identified multiple new candidates, including the transcription factors OsRLI and OsEIL1, which we demonstrated to induce the primary nitrate-responsive genes OsNRT1.1b and OsNIR1. Our network thus serves as a valuable resource to obtain novel insights in nitrogen signaling.

Suppression of the host antiviral response by non-infectious varicella zoster virus extracellular vesicles

Varicella zoster virus (VZV) reactivates from ganglionic sensory neurons to produce herpes zoster (shingles) in a unilateral dermatomal distribution, typically in the thoracic region. Reactivation not only heightens the risk of stroke and other neurological complications but also increases susceptibility to co-infections with various viral and bacterial pathogens at sites distant from the original infection. The mechanism by which VZV results in complications remote from the initial foci remains unclear. Small extracellular vesicles (sEVs) are membranous signaling structures that can deliver proteins and nucleic acids to modify the function of distal cells and tissues during normal physiological conditions. Although viruses have been documented to exploit the sEV machinery to propagate infection, the role of non-infectious sEVs released from VZV-infected neurons in viral spread and disease has not been studied. Using multi-omic approaches, we characterized the content of sEVs released from VZV-infected human sensory neurons (VZV sEVs). One viral protein was detected (immediate-early 62), as well as numerous immunosuppressive and vascular disease-associated host proteins and miRNAs that were absent in sEVs from uninfected neurons. Notably, VZV sEVs are non-infectious yet transcriptionally altered primary human cells, suppressing the antiviral type 1 interferon response and promoting neuroinvasion of a secondary pathogen in vivo. These results challenge our understanding of VZV infection, proposing that the virus may contribute to distant pathologies through non-infectious sEVs beyond the primary infection site. Furthermore, this study provides a previously undescribed immune-evasion mechanism induced by VZV that highlights the significance of non-infectious sEVs in early VZV pathogenesis.

IMPORTANCE

Varicella zoster virus (VZV) is a ubiquitous human virus that predominantly spreads by direct cell-cell contact and requires efficient and immediate host immune evasion strategies to spread. The mechanisms of immune evasion prior to virion entry have not been fully elucidated and represent a critical gap in our complete understanding of VZV pathogenesis. This study describes a previously unreported antiviral evasion strategy employed by VZV through the exploitation of the infected host cell's small extracellular vesicle (sEV) machinery. These findings suggest that non-infectious VZV sEVs could travel throughout the body, affecting cells remote from the site of infection and challenging the current understanding of VZV clinical disease, which has focused on local effects and direct infection. The significance of these sEVs in early VZV pathogenesis highlights the importance of further investigating their role in viral spread and secondary disease development to reduce systemic complications following VZV infections.

The frequency and function of nucleoprotein-specific CD8+ T cells are critical for heterosubtypic immunity against influenza virus infection

Cytotoxic T lymphocytes (CTLs) mediate host defense against viral and intracellular bacterial infections and tumors. However, the magnitude of CTL response and their function needed to confer heterosubtypic immunity against influenza virus infection are unknown. We addressed the role of CD8+ T cells in the absence of any cross-reactive antibody responses to influenza viral proteins using an adenoviral vector expressing a 9mer amino acid sequence recognized by CD8+ T cells. Our results indicate that both CD8+ T cell frequency and function are crucial for heterosubtypic immunity. Low morbidity, lower viral lung titers, low to minimal lung pathology, and better survival upon heterosubtypic virus challenge correlated with the increased frequency of NP-specific CTLs. NP-CD8+ T cells induced by differential infection doses displayed distinct RNA transcriptome profiles and functional properties. CD8+ T cells induced by a high dose of influenza virus secreted significantly higher levels of IFN-γ and exhibited higher levels of cytotoxic function. The mice that received NP-CD8+ T cells from the high-dose virus recipients through adoptive transfer had lower viral titers following viral challenge than those induced by the low dose of virus, suggesting differential cellular programming by antigen dose. Enhanced NP-CD8+ T-cell functions induced by a higher dose of influenza virus strongly correlated with the increased expression of cellular and metabolic genes, indicating a shift to a more glycolytic metabolic phenotype. These findings have implications for developing effective T cell vaccines against infectious diseases and cancer.

IMPORTANCE

Cytotoxic T lymphocytes (CTLs) are an important component of the adaptive immune system that clears virus-infected cells or tumor cells. Hence, developing next-generation vaccines that induce or recall CTL responses against cancer and infectious diseases is crucial. However, it is not clear if the frequency, function, or both are essential in conferring protection, as in the case of influenza. In this study, we demonstrate that both CTL frequency and function are crucial for providing heterosubtypic immunity to influenza by utilizing an Ad-viral vector expressing a CD8 epitope only to rule out the role of antibodies, single-cell RNA-seq analysis, as well as adoptive transfer experiments. Our findings have implications for developing T cell vaccines against infectious diseases and cancer.

Insulin-like growth factor 2 targets IGF1R signaling transduction to facilitate metastasis and imatinib resistance in gastrointestinal stromal tumors

BACKGROUND

Gastrointestinal stromal tumors (GISTs) are typical gastrointestinal tract neoplasms. Imatinib is the first-line therapy for GIST patients. Drug resistance limits the long-term effectiveness of imatinib. The regulatory effect of insulin-like growth factor 2 (IGF2) has been confirmed in various cancers and is related to resistance to chemotherapy and a worse prognosis.

AIM

To further investigate the mechanism of IGF2 specific to GISTs.

METHODS

IGF2 was screened and analyzed using Gene Expression Omnibus (GEO: GSE225819) data. After IGF2 knockdown or overexpression by transfection, the phenotypes (proliferation, migration, invasion, apoptosis) of GIST cells were characterized by cell counting kit 8, Transwell, and flow cytometry assays. We used western blotting to evaluate pathway-associated and epithelial-mesenchymal transition (EMT)-associated proteins. We injected transfected cells into nude mice to establish a tumor xenograft model and observed the occurrence and metastasis of GIST.

RESULTS

Data from the GEO indicated that IGF2 expression is high in GISTs, associated with liver metastasis, and closely related to drug resistance. GIST cells with high expression of IGF2 had increased proliferation and migration, invasiveness and EMT. Knockdown of IGF2 significantly inhibited those activities. In addition, OE-IGF2 promoted GIST metastasis in vivo in nude mice. IGF2 activated IGF1R signaling in GIST cells, and IGF2/IGF1R-mediated glycolysis was required for GIST with liver metastasis. GIST cells with IGF2 knockdown were sensitive to imatinib treatment when IGF2 overexpression significantly raised imatinib resistance. Moreover, 2-deoxy-D-glucose (a glycolysis inhibitor) treatment reversed IGF2 overexpression-mediated imatinib resistance in GISTs.

CONCLUSION

IGF2 targeting of IGF1R signaling inhibited metastasis and decreased imatinib resistance by driving glycolysis in GISTs.

Bipotential B-neutrophil progenitors are present in human and mouse bone marrow and emerge in the periphery upon stress hematopoiesis

Hematopoiesis is a tightly regulated process that gets skewed toward myelopoiesis. This restrains lymphopoiesis, but the role of lymphocytes in this process is not well defined. To unravel the intricacies of neutrophil responses in COVID-19, we performed bulk RNAseq on neutrophils from healthy controls and COVID-19 patients. Principal component analysis revealed distinguishing neutrophil gene expression alterations in COVID-19 patients. ICU and ward patients displayed substantial transcriptional changes, with ICU patients exhibiting a more pronounced response. Intriguingly, neutrophils from COVID-19 patients, notably ICU patients, exhibited an enrichment of immunoglobulin (Ig) and B cell lineage-associated genes, suggesting potential lineage plasticity. We validated our RNAseq findings in a larger cohort. Moreover, by reanalyzing single-cell RNA sequencing (scRNAseq) data on human bone marrow (BM) granulocytes, we identified the cluster of granulocyte-monocyte progenitors (GMP) enriched with Ig and B cell lineage-associated genes. These cells with lineage plasticity may serve as a resource depending on the host's needs during severe systemic infection. This distinct B cell subset may play a pivotal role in promoting myelopoiesis in response to infection. The scRNAseq analysis of BM neutrophils in infected mice further supported our observations in humans. Finally, our studies using an animal model of acute infection implicate IL-7/GM-CSF in influencing neutrophil and B cell dynamics. Elevated GM-CSF and reduced IL-7 receptor expression in COVID-19 patients imply altered hematopoiesis favoring myeloid cells over B cells. Our findings provide novel insights into the relationship between the B-neutrophil lineages during severe infection, hinting at potential implications for disease pathogenesis.

IMPORTANCE

This study investigates the dynamics of hematopoiesis in COVID-19, focusing on neutrophil responses. Through RNA sequencing of neutrophils from healthy controls and COVID-19 patients, distinct gene expression alterations are identified, particularly in ICU patients. Notably, neutrophils from COVID-19 patients, especially in the ICU, exhibit enrichment of immunoglobulin and B cell lineage-associated genes, suggesting potential lineage plasticity. Validation in a larger patient cohort and single-cell analysis of bone marrow granulocytes support the presence of granulocyte-monocyte progenitors with B cell lineage-associated genes. The findings propose a link between B-neutrophil lineages during severe infection, implicating a potential role for these cells in altered hematopoiesis favoring myeloid cells over B cells. Elevated GM-CSF and reduced IL-7 receptor expression in stress hematopoiesis suggest cytokine involvement in these dynamics, providing novel insights into disease pathogenesis.

Applications of Zebrafish Embryo Models to Predict Developmental Toxicity for Agrochemical Product Development

The development of safe crop protection products is a complex process that traditionally relies on intensive animal use for hazard identification. Methods that capture toxicity in early stages of agrochemical discovery programs enable a more efficient and sustainable product development pipeline. Here, we explored whether the zebrafish model can be leveraged to identify mammalian-relevant toxicity. We used transgenic zebrafish to assess developmental toxicity following exposures to known mammalian teratogens and captured larval morphological malformations, including bone and vascular perturbations. We further applied toxicogenomics to identify common biomarker signatures of teratogen exposure. The results show that the larval malformation assay predicted teratogenicity with 82.35% accuracy, 87.50% specificity, and 77.78% sensitivity. Similar and slightly lower accuracies were obtained with the vascular and bone assays, respectively. A set of 20 biomarkers were identified that efficiently segregated teratogenic chemicals from nonteratogens. In conclusion, zebrafish are valuable, robust, and cost-effective models for toxicity testing in the early stages of product development.