Next maSigPro: updating maSigPro bioconductor package for RNA-seq time series

Delta‐6 desaturase FADS2 is a tumor‐promoting factor in cholangiocarcinoma

Cholangiocarcinoma is a fatal disease with limited therapeutic options. We screened genes required for cholangiocarcinoma tumorigenicity and identified FADS2, a delta‐6 desaturase. FADS2 depletion reduced in vivo tumorigenicity and cell proliferation. In clinical samples, FADS2 was expressed in cancer cells but not in stromal cells. FADS2 inhibition also reduced the migration and sphere‐forming ability of cells and increased apoptotic cell death and ferroptosis markers. Lipidome assay revealed that triglyceride and cholesterol ester levels were decreased in FADS2‐knockdown cells. The oxygen consumption ratio was also decreased in FADS2‐depleted cells. These data indicate that FADS2 depletion causes a reduction in lipid levels, resulting in decrease of energy production and attenuation of cancer cell malignancy.

A longitudinal evaluation of oxidative stress - mitochondrial dysfunction - ferroptosis genes in anthracycline-induced cardiotoxicity

BACKGROUND

Antineoplastic medications, including doxorubicin, idarubicin, and epirubicin, have been found to adversely affect the heart due to oxidative stress - mitochondrial dysfunction - ferroptosis (ORMFs), which act as contributing attributes to anthracycline-induced cardiotoxicity. To better understand this phenomenon, the time-resolved measurements of ORMFS genes were analyzed in this study.

METHODS

The effect of three anthracycline drugs on ORMFs genes was studied using a human 3D cardiac microtissue cell model. Transcriptome data was collected over 14 days at two doses (therapeutic and toxic). WGCNA identified key module-related genes, and functional enrichment analysis investigated the biological processes quantified by ssGSEA, such as immune cell infiltration and angiogenesis. Biopsies were collected from heart failure patients and control subjects. GSE59672 and GSE2965 were collected for validation. Molecular docking was used to identify anthracyclines's interaction with key genes.

RESULTS

The ORMFs genes were screened in vivo or in vitro. Using WGCNA, six co-expressed gene modules were grouped, with MEblue emerging as the most significant module. Eight key genes intersecting the blue module with the dynamic response genes were obtained: CD36, CDH5, CHI3L1, HBA2, HSD11B1, OGN, RPL8, and VWF. Compared with control samples, all key genes except RPL8 were down-regulated in vitro ANT treatment settings, and their expression levels varied over time. According to functional analyses, the key module-related genes were engaged in angiogenesis and the immune system pathways. In all ANT-treated settings, ssGSEA demonstrated a significant down-regulation of angiogenesis score and immune cell activity, including Activated CD4 T cell, Immature B cell, Memory B cell, Natural killer cell, Type 1 T helper cell, and Type 2 T helper cell. Molecular docking revealed that RPL8 and CHI3L1 show significant binding affinity for anthracyclines.

CONCLUSION

This study focuses on the dynamic characteristics of ORMFs genes in both human cardiac microtissues and cardiac biopsies from ANT-treated patients. It has been highlighted that ORMFs genes may contribute to immune infiltration and angiogenesis in cases of anthracycline-induced cardiotoxicity. A thorough understanding of these genes could potentially lead to improved diagnosis and treatment of the disease.

scMaSigPro: Differential Expression Analysis along Single-Cell Trajectories

MOTIVATION

Understanding the dynamics of gene expression across different cellular states is crucial for discerning the mechanisms underneath cellular differentiation. Genes that exhibit variation in mean expression as a function of Pseudotime and between branching trajectories are expected to govern cell fate decisions. We introduce scMaSigPro, a method for the identification of differential gene expression patterns along Pseudotime and branching paths simultaneously.

RESULTS

We assessed the performance of scMaSigPro using synthetic and public datasets. Our evaluation shows that scMaSigPro outperforms existing methods in controlling the False Positive Rate and is computationally efficient.

AVAILABILITY AND IMPLEMENTATION

scMaSigPro is available as a free R package (version 4.0 or higher) under the GPL(≥2) license on GitHub at 'github.com/BioBam/scMaSigPro' and archived with version 0.03 on Zenodo at 'zenodo.org/records/12568922'.

Development of a screening system of gene sets for estimating the time of early skeletal muscle injury based on second-generation sequencing technology

The present study is aimed to address the challenge of wound age estimation in forensic science by identifying reliable genetic markers using low-cost and high-precision second-generation sequencing technology. A total of 54 Sprague-Dawley rats were randomly assigned to a control group or injury groups, with injury groups being further divided into time points (4 h, 8 h, 12 h, 16 h, 20 h, 24 h, 28 h, and 32 h after injury, n = 6) to establish rat skeletal muscle contusion models. Gene expression data were obtained using second-generation sequencing technology, and differential gene expression analysis, weighted gene co-expression network analysis (WGCNA) and time-dependent expression trend analysis were performed. A total of six sets of biomarkers were obtained: differentially expressed genes at adjacent time points (127 genes), co-expressed genes most associated with wound age (213 genes), hub genes exhibiting time-dependent expression (264 genes), and sets of transcription factors (TF) corresponding to the above sets of genes (74, 87, and 99 genes, respectively). Then, random forest (RF), support vector machine (SVM) and multilayer perceptron (MLP), were constructed for wound age estimation from the above gene sets. The results estimated by transcription factors were all superior to the corresponding hub genes, with the transcription factor group of WGCNA performed the best, with average accuracy rates of 96% for three models' internal testing, and 91.7% for the highest external validation. This study demonstrates the advantages of the indicator screening system based on second-generation sequencing technology and transcription factor level for wound age estimation.

Time-course cross-species transcriptomics reveals conserved hepatotoxicity pathways induced by repeated administration of cyclosporine A

Cyclosporine A (CsA) has shown efficacy against immunity-related diseases despite its toxicity in various organs, including the liver, emphasizing the need to elucidate its underlying hepatotoxicity mechanism. This study aimed to capture the alterations in genome-wide expression over time and the subsequent perturbations of corresponding pathways across species. Six data from humans, mice, and rats, including animal liver tissue, human liver microtissues, and two liver cell lines exposed to CsA toxic dose, were used. The microtissue exposed to CsA for 10 d was analyzed to obtain dynamically differentially expressed genes (DEGs). Single-time points data at 1, 3, 5, 7, and 28 d of different species were used to provide additional evidence. Using liver microtissue-based longitudinal design, DEGs that were consistently up- or down-regulated over time were captured, and the well-known mechanism involved in CsA toxicity was elucidated. Thirty DEGs that consistently changed in longitudinal data were also altered in 28-d rat in-house data with concordant expression. Some genes (e.g. TUBB2A, PLIN2, APOB) showed good concordance with identified DEGs in 1-d and 7-d mouse data. Pathway analysis revealed up-regulations of protein processing, asparagine N-linked glycosylation, and cargo concentration in the endoplasmic reticulum. Furthermore, the down-regulations of pathways related to biological oxidations and metabolite and lipid metabolism were elucidated. These pathways were also enriched in single-time-point data and conserved across species, implying their biological significance and generalizability. Overall, the human organoids-based longitudinal design coupled with cross-species validation provides temporal molecular change tracking, aiding mechanistic elucidation and biologically relevant biomarker discovery.

Transcriptomic profiling of gastrointestinal tracts in dairy cattle during lactation reveals molecular adaptations for milk synthesis

During lactation, dairy cattle's digestive tract requires significant adaptations to meet the increased nutrient demands for milk production. As we attempt to improve milk-related traits through selective pressure, it is crucial to understand the biological functions of the epithelia of the rumen, small intestine, and colonic tissues in response to changes in physiological state driven by changes in nutrient demands for milk synthesis. In this study, we obtained a total of 108 transcriptome profiles from three tissues (epithelia of the colon, duodenum, and rumen) of five Holstein cows, spanning eight time points from the early, mid, late lactation periods to the dry period. On average 97.06% of reads were successfully mapped to the reference genome assembly ARS-UCD1.2. We analyzed 27,607 gene expression patterns at multiple periods, enabling direct comparisons within and among tissues during different lactation stages, including early and peak lactation. We identified 1645, 813, and 2187 stage-specific genes in the colon, duodenum, and rumen, respectively, which were enriched for common or specific biological functions among different tissues. Time series analysis categorized the expressed genes within each tissue into four clusters. Furthermore, when the three tissues were analyzed collectively, 36 clusters of similarly expressed genes were identified. By integrating other comprehensive approaches such as gene co-expression analyses, functional enrichment, and cell type deconvolution, we gained profound insights into cattle lactation, revealing tissue-specific characteristics of the gastrointestinal tract and shedding light on the intricate molecular adaptations involved in nutrient absorption, immune regulation, and cellular processes for milk synthesis during lactation.

Pathogenic differences of cynomolgus macaques after Taï Forest virus infection depend on the viral stock propagation

Taï Forest virus (TAFV) is a negative-sense RNA virus in the Filoviridae family. TAFV has caused only a single human infection, but several disease outbreaks in chimpanzees have been linked to this virus. Limited research has been done on this human-pathogenic virus. We sought to establish an animal model to assess TAFV disease progression and pathogenicity at our facility. We had access to two different viral stock preparations from different institutions, both originating from the single human case. Type I interferon receptor knockout mice were inoculated with TAFV stock 1 or stock 2 by the intraperitoneal route. Inoculation resulted in 100% survival with no disease regardless of viral stock preparation or infectious dose. Next, cynomolgus macaques were inoculated with TAFV stock 1 or stock 2. Inoculation with TAFV stock 1 resulted in 100% survival and robust TAFV glycoprotein-specific IgG responses including neutralizing antibodies. In contrast, macaques infected with TAFV stock 2 developed disease and were euthanized 8-11 days after infection exhibiting viremia, thrombocytopenia, and increased inflammatory mediators identified by transcriptional analysis. Histopathologic analysis of tissue samples collected at necropsy confirmed classic filovirus disease in numerous organs. Genomic differences in both stock preparations were mapped to several viral genes which may have contributed to disease severity. Taken together, we demonstrate that infection with the two TAFV stocks resulted in no disease in mice and opposing disease phenotypes in cynomolgus macaques, highlighting the impact of viral stock propagation on pathogenicity in animal models.

FAXC interacts with ANXA2 and SRC in mitochondria and promotes tumorigenesis in cholangiocarcinoma

Cholangiocarcinoma (CCA) is one of the most difficult malignancies to treat as the therapeutic options are limited. Although several driver genes have been identified, most remain unknown. In this study, we identified a failed axon connection homolog (FAXC), whose function is unknown in mammals, by analyzing serially passaged CCA xenograft models. Knockdown of FAXC reduced subcutaneous tumorigenicity in mice. FAXC was bound to annexin A2 (ANXA2) and c-SRC, which are tumor-promoting genes. The FAXC/ANXA2/c-SRC complex forms in the mitochondria. FAXC enhances SRC-dependent ANXA2 phosphorylation at tyrosine-24, and the C-terminal amino acid residues (351-375) of FAXC are required for ANXA2 phosphorylation. Transcriptome data from a xenografted CCA cell line revealed that FAXC correlated with epithelial-mesenchymal transition, hypoxia, and KRAS signaling genes. Collectively, these findings advance our understanding of CCA tumorigenesis and provide candidate therapeutic targets.

Novel biomarkers of preterm brain injury from blood transcriptome in sheep model of intrauterine asphyxia

BACKGROUND

Infants born preterm have a higher incidence of neurological deficits. A key step in finding effective treatments is to identify biomarkers that reliably predict outcome.

METHODS

Following umbilical cord occlusion (UCO) in pregnant sheep, whole fetal blood RNA was sequenced pre- and post-UCO, brain injury outcome was determined by battery of neuropathology scoring and the transcriptome signature correlated to the degree of brain injury. Additionally, we developed a novel analytical procedure to deduce cell blood composition over time.

RESULTS

Sixty-one genes were identified with significant altered expression after UCO. In pre-UCO blood, the level of three mRNAs (Trex2, Znf280b, novel miRNA) and in post-UCO, four mRNAs (Fam184a, Angptl2, novel lincRNA and an unknown protein-coding gene) were associated to brain injury (FDR < 0.01). Several of these mRNAs are related to inflammation and angiogenesis. Pathway analysis highlighted genes playing a role in perinatal death and growth failure. Results also indicate that several leukocyte populations undergo significant changes after UCO.

CONCLUSION

We have used a whole transcriptomic approach to uncover novel biomarkers in fetal blood that correlate to neuropathology in the preterm sheep brain. The current data forms a basis for future studies to investigate mechanisms of these mRNAs in the injury progression.

IMPACT

Trend analysis of genes following asphyxia reveal a group of genes associated with perinatal death and growth failure. Several pre-asphyxia transcripts were associated to brain injury severity suggesting genomic susceptibility to injury. Several post-asphyxia transcripts were correlated to brain injury severity, thus, serve as potential novel biomarkers of injury outcome. Successfully adaptation of cell profiling algorithms suggests significant changes in blood cell composition following asphyxia.

Specific plasma microRNA profiles could be potential non-invasive biomarkers for biochemical pregnancy loss following embryo transfer

BACKGROUND

Plasma microRNAs act as biomarkers for predicting and diagnosing diseases. Reliable non-invasive biomarkers for biochemical pregnancy loss have not been established. We aim to analyze the dynamic microRNA profiles during the peri-implantation period and investigate if plasma microRNAs could be non-invasive biomarkers predicting BPL.

METHODS

In this study, we collected plasma samples from patients undergoing embryo transfer (ET) on ET day (ET0), 11 days after ET (ET11), and 14 days after ET (ET14). Patients were divided into the NP (negative pregnancy), BPL (biochemical pregnancy loss), and CP (clinical pregnancy) groups according to serum hCG levels at day11~14 and ultrasound at day28~35 following ET. MicroRNA profiles at different time-points were detected by miRNA-sequencing. We analyzed plasma microRNA signatures for BPL at the peri-implantation stage, we characterized the dynamic microRNA changes during the implantation period, constructed a microRNA co-expression network, and established predictive models for BPL. Finally, the sequencing results were confirmed by Taqman RT-qPCR.

RESULTS

BPL patients have distinct plasma microRNA profiles compared to CP patients at multiple time-points during the peri-implantation period. Machine learning models revealed that plasma microRNAs could predict BPL. RT-qPCR confirmed that miR-181a-2-3p, miR-9-5p, miR-150-3p, miR-150-5p, and miR-98-5p, miR-363-3p were significantly differentially expressed between patients with different reproductive outcomes.

CONCLUSION

Our study highlights the non-invasive value of plasma microRNAs in predicting BPL.

Local age-dependent neuromodulation in Rhodnius prolixus antennae

Kissing bugs do not respond to host cues when recently molted and only exhibit robust host-seeking several days after ecdysis. Behavioral plasticity has peripheral correlates in antennal gene expression changes through the week after ecdysis. The mechanisms regulating these peripheral changes are still unknown, but neuropeptide, G-protein coupled receptor, nuclear receptor, and takeout genes likely modulate peripheral sensory physiology. We evaluated their expression in antennal transcriptomes along the first week postecdysis of Rhodnius prolixus 5th instar larvae. Besides, we performed clustering and co-expression analyses to reveal relationships between neuromodulatory (NM) and sensory genes. Significant changes in transcript abundance were detected for 50 NM genes. We identified 73 sensory-related and NM genes that were assigned to nine clusters. According to their expression patterns, clusters were classified into four groups: two including genes up or downregulated immediately after ecdysis; and two with genes with expression altered at day 2. Several NM genes together with sensory genes belong to the first group, suggesting functional interactions. Co-expression network analysis revealed a set of genes that seem to connect with sensory system maturation. Significant expression changes in NM components were described in the antennae of R. prolixus after ecdysis, suggesting that a local NM system acts on antennal physiology. These changes may modify the sensitivity of kissing bugs to host cues during this maturation interval.

Exploring the effects of calycosin on anthracycline-induced cardiotoxicity: a network pharmacology, molecular docking, and experimental study

Background

Chemotherapy with anthracyclines can cause cardiotoxicity, possibly leading to stopping treatment in some cancer patients. In cardio-oncology research, preventing and minimizing anthracycline-induced cardiotoxicity (AIC) is a hot issue. For the treatment of AIC, calycosin (CA), an isoflavone component in astragali radix (AR), has become a research focus. However, the elaborate mechanisms of calycosin treating AIC remain to be unrevealed.

Aim of the study

To explore the effects of CA on AIC through multiple dimensions concerning network pharmacology, molecular docking, and experimental evaluations.

Methods

The study evaluated calycosin's potential targets and mechanisms for treating AIC using network pharmacology and molecular docking. The candidate genes/targets of CA and AIC were screened using the online-available database. Protein-protein interactions (PPI) between the common targets were constructed using the STRING platform, and the results were then visualized using Cytoscape. Molecular docking was used to evaluate the strength of the binding force between CA and the common targets. The possible pharmacological mechanisms of CA were explained by pathway enrichment and GSEA. Subsequently, the candidate targets were identified in vitro experiments.

Results

Network pharmacology effectively discovered the CA's multitarget intervention in AIC, including TNF, ABCC1, TOP2A, ABCB1, and XDH. CA binds to the ATP-binding cassette subfamily B member 1(ABCB1) had the highest binding energy (-7.5 kcal/mol) according to the molecular docking analysis and was selected and visualized for subsequent analysis. In vitro experiments showed that ABCB1 exhibited significant time-curve changes under different doses of doxorubicin (DOX) compared with DMSO control experiments. The anti-AIC pharmacological mechanism of CA were revealed by highlighting the biological processes of oxidative stress (OR) and inflammation.

Conclusions

We employed a practicable bioinformatics method to connect network and molecular docking to determine the calycosin's therapeutic mechanism against AIC and identified some bioinformatics results in in vitro experiments. The results presented show that CA may represent an encouraging treatment for AIC.

A clustering procedure for three-way RNA sequencing data using data transformations and matrix-variate Gaussian mixture models

RNA sequencing of time-course experiments results in three-way count data where the dimensions are the genes, the time points and the biological units. Clustering RNA-seq data allows to extract groups of co-expressed genes over time. After standardisation, the normalised counts of individual genes across time points and biological units have similar properties as compositional data. We propose the following procedure to suitably cluster three-way RNA-seq data: (1) pre-process the RNA-seq data by calculating the normalised expression profiles, (2) transform the data using the additive log ratio transform to map the composition in the D-part Aitchison simplex to a D - 1 -dimensional Euclidean vector, (3) cluster the transformed RNA-seq data using matrix-variate Gaussian mixture models and (4) assess the quality of the overall cluster solution and of individual clusters based on cluster separation in the transformed space using density-based silhouette information and on compactness of the cluster in the original space using cluster maps as a suitable visualisation. The proposed procedure is illustrated on RNA-seq data from fission yeast and results are also compared to an analogous two-way approach after flattening out the biological units.

Primary cilia promote the differentiation of human neurons through the WNT signaling pathway

BACKGROUND

Primary cilia emanate from most human cell types, including neurons. Cilia are important for communicating with the cell's immediate environment: signal reception and transduction to/from the ciliated cell. Deregulation of ciliary signaling can lead to ciliopathies and certain neurodevelopmental disorders. In the developing brain cilia play well-documented roles for the expansion of the neural progenitor cell pool, while information about the roles of cilia during post-mitotic neuron differentiation and maturation is scarce.

RESULTS

We employed ciliated Lund Human Mesencephalic (LUHMES) cells in time course experiments to assess the impact of ciliary signaling on neuron differentiation. By comparing ciliated and non-ciliated neuronal precursor cells and neurons in wild type and in RFX2 -/- mutant neurons with altered cilia, we discovered an early-differentiation "ciliary time window" during which transient cilia promote axon outgrowth, branching and arborization. Experiments in neurons with IFT88 and IFT172 ciliary gene knockdowns, leading to shorter cilia, confirm these results. Cilia promote neuron differentiation by tipping WNT signaling toward the non-canonical pathway, in turn activating WNT pathway output genes implicated in cyto-architectural changes.

CONCLUSIONS

We provide a mechanistic entry point into when and how ciliary signaling coordinates, promotes and translates into anatomical changes. We hypothesize that ciliary alterations causing neuron differentiation defects may result in "mild" impairments of brain development, possibly underpinning certain aspects of neurodevelopmental disorders.

miRNA profiling of chicken follicles during follicular development

MicroRNAs (miRNAs) play a crucial role as transcription regulators in various aspects of follicular development, including steroidogenesis, ovulation, apoptosis, and gene regulation in poultry. However, there is a paucity of studies examining the specific impact of miRNAs on ovarian granulosa cells (GCs) across multiple grades in laying hens. Consequently, this study aims to investigate the roles of miRNAs in chicken GCs. By constructing miRNA expression profiles of GCs at 10 different time points, encompassing 4 pre-hierarchical, 5 preovulatory, and 1 postovulatory follicles stage, we identified highly expressed miRNAs involved in GC differentiation (miR-148a-3p, miR-143-3p), apoptosis (let7 family, miR-363-3p, miR-30c-5p, etc.), and autophagy (miR-128-3p, miR-21-5p). Furthermore, we discovered 48 developmentally dynamic miRNAs (DDMs) that target 295 dynamic differentially expressed genes (DDGs) associated with follicular development and selection (such as oocyte meiosis, progesterone-mediated oocyte maturation, Wnt signaling pathway, TGF-β signaling pathway) as well as follicular regression (including autophagy and cellular senescence). These findings contribute to a more comprehensive understanding of the intricate mechanisms underlying follicle recruitment, selection, and degeneration, aiming to enhance poultry's reproductive capacity.

CAP-RNAseq: an integrated pipeline for functional annotation and prioritization of co-expression clusters

Cluster analysis is one of the most widely used exploratory methods for visualization and grouping of gene expression patterns across multiple samples or treatment groups. Although several existing online tools can annotate clusters with functional terms, there is no all-in-one webserver to effectively prioritize genes/clusters using gene essentiality as well as congruency of mRNA-protein expression. Hence, we developed CAP-RNAseq that makes possible (1) upload and clustering of bulk RNA-seq data followed by identification, annotation and network visualization of all or selected clusters; and (2) prioritization using DepMap gene essentiality and/or dependency scores as well as the degree of correlation between mRNA and protein levels of genes within an expression cluster. In addition, CAP-RNAseq has an integrated primer design tool for the prioritized genes. Herein, we showed using comparisons with the existing tools and multiple case studies that CAP-RNAseq can uniquely aid in the discovery of co-expression clusters enriched with essential genes and prioritization of novel biomarker genes that exhibit high correlations between their mRNA and protein expression levels. CAP-RNAseq is applicable to RNA-seq data from different contexts including cancer and available at http://konulabapps.bilkent.edu.tr:3838/CAPRNAseq/ and the docker image is downloadable from https://hub.docker.com/r/konulab/caprnaseq.

Modeling transcriptional regulation of the cell cycle using a novel cybernetic-inspired approach

Quantitative understanding of cellular processes, such as cell cycle and differentiation, is impeded by various forms of complexity ranging from myriad molecular players and their multilevel regulatory interactions, cellular evolution with multiple intermediate stages, lack of elucidation of cause-effect relationships among the many system players, and the computational complexity associated with the profusion of variables and parameters. In this paper, we present a modeling framework based on the cybernetic concept that biological regulation is inspired by objectives embedding rational strategies for dimension reduction, process stage specification through the system dynamics, and innovative causal association of regulatory events with the ability to predict the evolution of the dynamical system. The elementary step of the modeling strategy involves stage-specific objective functions that are computationally determined from experiments, augmented with dynamical network computations involving endpoint objective functions, mutual information, change-point detection, and maximal clique centrality. We demonstrate the power of the method through application to the mammalian cell cycle, which involves thousands of biomolecules engaged in signaling, transcription, and regulation. Starting with a fine-grained transcriptional description obtained from RNA sequencing measurements, we develop an initial model, which is then dynamically modeled using the cybernetic-inspired method, based on the strategies described above. The cybernetic-inspired method is able to distill the most significant interactions from a multitude of possibilities. In addition to capturing the complexity of regulatory processes in a mechanistically causal and stage-specific manner, we identify the functional network modules, including novel cell cycle stages. Our model is able to predict future cell cycles consistent with experimental measurements. We posit that this innovative framework has the promise to extend to the dynamics of other biological processes, with a potential to provide novel mechanistic insights.

Supergene evolution via gain of auto-regulation

The development of complex phenotypes requires the coordinated action of many genes across space and time, yet many species have evolved the ability to develop multiple discrete, alternate phenotypes1-3. Such polymorphisms are often controlled by supergenes, sets of tightly-linked mutations in one or more loci that function together to produce a complex phenotype4. Although theories of supergene evolution are well-established, the mutations that cause functional differences between supergene alleles remain essentially unknown. doublesex is the master regulator of insect sexual differentiation but functions as a supergene in multiple Papilio swallowtail butterflies, where divergent dsx alleles control development of discrete non-mimetic or mimetic female wing color patterns5-7. Here we demonstrate that the functional elements of the mimetic allele in Papilio alphenor are six new cis-regulatory elements (CREs) spread across 150 kb that are bound by DSX itself. Our findings provide experimental support to classic supergene theory and suggest that the evolution of auto-regulation may provide a simple route to supergene origination and to the co-option of pleiotropic genes into new developmental roles.

Dynamic profiles of lncRNAs reveal a functional natural antisense RNA that regulates the development of Schistosoma japonicum

Schistosomes are flatworm parasites that undergo a complex life cycle involving two hosts. The regulation of the parasite's developmental processes relies on both coding RNAs and non-coding RNAs. However, the roles of non-coding RNAs, including long non-coding RNAs (lncRNAs) in schistosomes remain largely unexplored. Here we conduct advanced RNA sequencing on male and female S. japonicum during their pairing and reproductive development, resulting in the identification of nearly 8,000 lncRNAs. This extensive dataset enables us to construct a comprehensive co-expression network of lncRNAs and mRNAs, shedding light on their interactions during the crucial reproductive stages within the mammalian host. Importantly, we have also revealed a specific lncRNA, LNC3385, which appears to play a critical role in the survival and reproduction of the parasite. These findings not only enhance our understanding of the dynamic nature of lncRNAs during the reproductive phase of schistosomes but also highlight LNC3385 as a potential therapeutic target for combating schistosomiasis.

Transcriptional profiling of immune responses in NHPs after low-dose, VSV-based vaccination against Marburg virus

Infection with Marburg virus (MARV), the causative agent of Marburg virus disease (MVD), results in haemorrhagic disease and high case fatality rates (>40%) in humans. Despite its public health relevance, there are no licensed vaccines or therapeutics to prevent or treat MVD. A vesicular stomatitis virus (VSV)-based vaccine expressing the MARV glycoprotein (VSV-MARV) is currently in clinical development. Previously, a single 10 million PFU dose of VSV-MARV administered 1-5 weeks before lethal MARV challenge conferred uniform protection in nonhuman primates (NHPs), demonstrating fast-acting potential. Additionally, our group recently demonstrated that even a low dose VSV-MARV (1000 PFU) protected NHPs when given 7 days before MARV challenge. In this study, we longitudinally profiled the transcriptional responses of NHPs vaccinated with this low dose of VSV-MARV either 14 or 7 days before lethal MARV challenge. NHPs vaccinated 14 days before challenge presented with transcriptional changes consistent with an antiviral response before challenge. Limited gene expression changes were observed in the group vaccinated 7 days before challenge. After challenge, genes related to lymphocyte-mediated immunity were only observed in the group vaccinated 14 days before challenge, indicating that the length of time between vaccination and challenge influenced gene expression. Our results indicate that a low dose VSV-MARV elicits distinct immune responses that correlate with protection against MVD. A low dose of VSV-MARV should be evaluated in clinical rails as it may be an option to deliver beneficial public health outcomes to more people in the event of future outbreaks.

Regulatory T cells shield muscle mitochondria from interferon-γ-mediated damage to promote the beneficial effects of exercise

Exercise enhances physical performance and reduces the risk of many disorders such as cardiovascular disease, type 2 diabetes, dementia, and cancer. Exercise characteristically incites an inflammatory response, notably in skeletal muscles. Although some effector mechanisms have been identified, regulatory elements activated in response to exercise remain obscure. Here, we have addressed the roles of Foxp3+CD4+ regulatory T cells (Tregs) in the healthful activities of exercise via immunologic, transcriptomic, histologic, metabolic, and biochemical analyses of acute and chronic exercise models in mice. Exercise rapidly induced expansion of the muscle Treg compartment, thereby guarding against overexuberant production of interferon-γ and consequent metabolic disruptions, particularly mitochondrial aberrancies. The performance-enhancing effects of exercise training were dampened in the absence of Tregs. Thus, exercise is a natural Treg booster with therapeutic potential in disease and aging contexts.

Differential gene expression during floral transition in pineapple

Pineapple (Ananas comosus var. comosus) and ornamental bromeliads are commercially induced to flower by treatment with ethylene or its analogs. The apex is transformed from a vegetative to a floral meristem and shows morphological changes in 8 to 10 days, with flowers developing 8 to 10 weeks later. During eight sampling stages ranging from 6 h to 8 days after treatment, 7961 genes were found to exhibit differential expression (DE) after the application of ethylene. In the first 3 days after treatment, there was little change in ethylene synthesis or in the early stages of the ethylene response. Subsequently, three ethylene response transcription factors (ERTF) were up-regulated and the potential gene targets were predicted to be the positive flowering regulator CONSTANS-like 3 (CO), a WUSCHEL gene, two APETALA1/FRUITFULL (AP1/FUL) genes, an epidermal patterning gene, and a jasmonic acid synthesis gene. We confirm that pineapple has lost the flowering repressor FLOWERING LOCUS C. At the initial stages, the SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) was not significantly involved in this transition. Another WUSCHEL gene and a PHD homeobox transcription factor, though not apparent direct targets of ERTF, were up-regulated within a day of treatment, their predicted targets being the up-regulated CO, auxin response factors, SQUAMOSA, and histone H3 genes with suppression of abscisic acid response genes. The FLOWERING LOCUS T (FT), TERMINAL FLOWER (TFL), AGAMOUS-like APETELAR (AP2), and SEPETALA (SEP) increased rapidly within 2 to 3 days after ethylene treatment. Two FT genes were up-regulated at the apex and not at the leaf bases after treatment, suggesting that transport did not occur. These results indicated that the ethylene response in pineapple and possibly most bromeliads act directly to promote the vegetative to flower transition via APETALA1/FRUITFULL (AP1/FUL) and its interaction with SPL, FT, TFL, SEP, and AP2. A model based on AP2/ERTF DE and predicted DE target genes was developed to give focus to future research. The identified candidate genes are potential targets for genetic manipulation to determine their molecular role in flower transition.

The flavour of grape colour: anthocyanin content tunes aroma precursor composition by altering the berry microenvironment

Anthocyaninless (white) instead of black/red (coloured) fruits develop in grapevine cultivars without functional VviMYBA1 and VviMYBA2 genes, and this conditions the colour of wines that can be produced. To evaluate whether this genetic variation has additional consequences on fruit ripening and composition, we performed comparisons of microenvironment, transcriptomics, and metabolomics of developing grapes between near-isogenic white- and black-berried somatic variants of Garnacha and Tempranillo cultivars. Berry temperature was as much as 3.5 ºC lower in white- compared to black-berried Tempranillo. An RNA-seq study combined with targeted and untargeted metabolomics revealed that ripening fruits of white-berried variants were characterized by the up-regulation of photosynthesis-related and other light-responsive genes and by their higher accumulation of specific terpene aroma precursors, fatty acid-derived aldehyde volatiles, and phenylpropanoid precursor amino acids. MYBA1-MYBA2 function proved essential for flavonol trihydroxylation in black-berried somatic variants, which were also characterized by enhanced expression of pathogen defence genes in the berry skin and increased accumulation of C6-derived alcohol and ester volatiles and γ-aminobutyric acid. Collectively, our results indicate that anthocyanin depletion has side-effects on grape composition by altering the internal microenvironment of the berry and the partitioning of the phenylpropanoid pathway. Our findings show how fruit colour can condition other fruit features, such as flavour potential and stress homeostasis.

Comparative analyses of dynamic transcriptome profiles highlight key response genes and dominant isoforms for muscle development and growth in chicken

BACKGROUND

Modern breeding strategies have resulted in significant differences in muscle mass between indigenous chicken and specialized broiler. However, the molecular regulatory mechanisms that underlie these differences remain elusive. The aim of this study was to identify key genes and regulatory mechanisms underlying differences in breast muscle development between indigenous chicken and specialized broiler.

RESULTS

Two time-series RNA-sequencing profiles of breast muscles were generated from commercial Arbor Acres (AA) broiler (fast-growing) and Chinese indigenous Lushi blue-shelled-egg (LS) chicken (slow-growing) at embryonic days 10, 14, and 18, and post-hatching day 1 and weeks 1, 3, and 5. Principal component analysis of the transcriptome profiles showed that the top four principal components accounted for more than 80% of the total variance in each breed. The developmental axes between the AA and LS chicken overlapped at the embryonic stages but gradually separated at the adult stages. Integrative investigation of differentially-expressed transcripts contained in the top four principal components identified 44 genes that formed a molecular network associated with differences in breast muscle mass between the two breeds. In addition, alternative splicing analysis revealed that genes with multiple isoforms always had one dominant transcript that exhibited a significantly higher expression level than the others. Among the 44 genes, the TNFRSF6B gene, a mediator of signal transduction pathways and cell proliferation, harbored two alternative splicing isoforms, TNFRSF6B-X1 and TNFRSF6B-X2. TNFRSF6B-X1 was the dominant isoform in both breeds before the age of one week. A switching event of the dominant isoform occurred at one week of age, resulting in TNFRSF6B-X2 being the dominant isoform in AA broiler, whereas TNFRSF6B-X1 remained the dominant isoform in LS chicken. Gain-of-function assays demonstrated that both isoforms promoted the proliferation of chicken primary myoblasts, but only TNFRSF6B-X2 augmented the differentiation and intracellular protein content of chicken primary myoblasts.

CONCLUSIONS

For the first time, we identified several key genes and dominant isoforms that may be responsible for differences in muscle mass between slow-growing indigenous chicken and fast-growing commercial broiler. These findings provide new insights into the regulatory mechanisms underlying breast muscle development in chicken.

Distinct defects in early innate and late adaptive immune responses typify impaired fracture healing in diet-induced obesity

Bone fractures, the most common musculoskeletal injuries, heal through three main phases: inflammatory, repair, and remodeling. Around 10% of fracture patients suffer from impaired healing that requires surgical intervention, a huge burden on the healthcare system. The rate of impaired healing increases with metabolic diseases such as obesity-associated hyperglycemia/type 2 diabetes (T2D), an increasing concern given the growing incidence of obesity/T2D. Immune cells play pivotal roles in fracture healing, and obesity/T2D is associated with defective immune-cell functions. However, there is a gap in knowledge regarding the stoichiometry of immune cells that populate the callus and how that population changes during different phases of healing. Here, we used complementary global and single-cell techniques to characterize the repertoire of immune cells in the fracture callus and to identify populations specifically enriched in the fracture callus relative to the unfractured bone or bone marrow. Our analyses identified two clear waves of immune-cell infiltration into the callus: the first wave occurs during the early inflammatory phase of fracture healing, while the second takes place during the late repair/early remodeling phase, which is consistent with previous publications. Comprehensive analysis of each wave revealed that innate immune cells were activated during the early inflammatory phase, but in later phases they returned to homeostatic numbers and activation levels. Of the innate immune cells, distinct subsets of activated dendritic cells were particularly enriched in the inflammatory healing hematoma. In contrast to innate cells, lymphocytes, including B and T cells, were enriched and activated in the callus primarily during the late repair phase. The Diet-Induced Obesity (DIO) mouse, an established model of obesity-associated hyperglycemia and insulin resistance, suffers from multiple healing defects. Our data demonstrate that DIO mice exhibit dysregulated innate immune responses during the inflammatory phase, and defects in all lymphocyte compartments during the late repair phase. Taken together, our data characterize, for the first time, immune populations that are enriched/activated in the callus during two distinct phases of fracture healing and identify defects in the healing-associated immune response in DIO mice, which will facilitate future development of immunomodulatory therapeutics for impaired fracture healing.

Chromosome-level genome assembly and population genetic analysis of a near-threatened rosewood species (Dalbergia cultrata Pierre Graham ex Benth) provide insights into its evolutionary and cold stress responses

Dalbergia cultrata Pierre Graham ex Benth (D. cultrata) is a precious rosewood tree species that grows in the tropical and subtropical regions of Asia. In this study, we used PacBio long-reading sequencing technology and Hi-C assistance to sequence and assemble the reference genome of D. cultrata. We generated 171.47 Gb PacBio long reads and 72.43 Gb Hi-C data and yielded an assembly of 10 pseudochromosomes with a total size of 690.99 Mb and Scaffold N50 of 65.76 Mb. The analysis of specific genes revealed that the triterpenoids represented by lupeol may play an important role in D. cultrata's potential medicinal value. Using the new reference genome, we analyzed the resequencing of 19 Dalbergia accessions and found that D. cultrata and D. cochinchinensis have the latest genetic relationship. Transcriptome sequencing of D. cultrata leaves grown under cold stress revealed that MYB transcription factor and E3 ubiquitin ligase may be playing an important role in the cold response of D. cultrata. Genome resources and identified genetic variation, especially those genes related to the biosynthesis of phytochemicals and cold stress response, will be helpful for the introduction, domestication, utilization, and further breeding of Dalbergia species.

Optimization of universal allogeneic CAR-T cells combining CRISPR and transposon-based technologies for treatment of acute myeloid leukemia

Despite the potential of CAR-T therapies for hematological malignancies, their efficacy in patients with relapse and refractory Acute Myeloid Leukemia has been limited. The aim of our study has been to develop and manufacture a CAR-T cell product that addresses some of the current limitations. We initially compared the phenotype of T cells from AML patients and healthy young and elderly controls. This analysis showed that T cells from AML patients displayed a predominantly effector phenotype, with increased expression of activation (CD69 and HLA-DR) and exhaustion markers (PD1 and LAG3), in contrast to the enriched memory phenotype observed in healthy donors. This differentiated and more exhausted phenotype was also observed, and corroborated by transcriptomic analyses, in CAR-T cells from AML patients engineered with an optimized CAR construct targeting CD33, resulting in a decreased in vivo antitumoral efficacy evaluated in xenograft AML models. To overcome some of these limitations we have combined CRISPR-based genome editing technologies with virus-free gene-transfer strategies using Sleeping Beauty transposons, to generate CAR-T cells depleted of HLA-I and TCR complexes (HLA-IKO/TCRKO CAR-T cells) for allogeneic approaches. Our optimized protocol allows one-step generation of edited CAR-T cells that show a similar phenotypic profile to non-edited CAR-T cells, with equivalent in vitro and in vivo antitumoral efficacy. Moreover, genomic analysis of edited CAR-T cells revealed a safe integration profile of the vector, with no preferences for specific genomic regions, with highly specific editing of the HLA-I and TCR, without significant off-target sites. Finally, the production of edited CAR-T cells at a larger scale allowed the generation and selection of enough HLA-IKO/TCRKO CAR-T cells that would be compatible with clinical applications. In summary, our results demonstrate that CAR-T cells from AML patients, although functional, present phenotypic and functional features that could compromise their antitumoral efficacy, compared to CAR-T cells from healthy donors. The combination of CRISPR technologies with transposon-based delivery strategies allows the generation of HLA-IKO/TCRKO CAR-T cells, compatible with allogeneic approaches, that would represent a promising option for AML treatment.

Impaired resolution of blood transcriptomes through tuberculosis treatment with diabetes comorbidity

BACKGROUND

People with diabetes are more likely to develop tuberculosis (TB) and to have poor TB-treatment outcomes than those without. We previously showed that blood transcriptomes in people with TB-diabetes (TB-DM) co-morbidity have excessive inflammatory and reduced interferon responses at diagnosis. It is unknown whether this persists through treatment and contributes to the adverse outcomes.

METHODS

Pulmonary TB patients recruited in South Africa, Indonesia and Romania were classified as having TB-DM, TB with prediabetes, TB-related hyperglycaemia or TB-only, based on glycated haemoglobin concentration at TB diagnosis and after 6 months of TB treatment. Gene expression in blood at diagnosis and intervals throughout treatment was measured by unbiased RNA-Seq and targeted Multiplex Ligation-dependent Probe Amplification. Transcriptomic data were analysed by longitudinal mixed-model regression to identify whether genes were differentially expressed between clinical groups through time. Predictive models of TB-treatment response across groups were developed and cross-tested.

RESULTS

Gene expression differed between TB and TB-DM patients at diagnosis and was modulated by TB treatment in all clinical groups but to different extents, such that differences remained in TB-DM relative to TB-only throughout. Expression of some genes increased through TB treatment, whereas others decreased: some were persistently more highly expressed in TB-DM and others in TB-only patients. Genes involved in innate immune responses, anti-microbial immunity and inflammation were significantly upregulated in people with TB-DM throughout treatment. The overall pattern of change was similar across clinical groups irrespective of diabetes status, permitting models predictive of TB treatment to be developed.

CONCLUSIONS

Exacerbated transcriptome changes in TB-DM take longer to resolve during TB treatment, meaning they remain different from those in uncomplicated TB after treatment completion. This may indicate a prolonged inflammatory response in TB-DM, requiring prolonged treatment or host-directed therapy for complete cure. Development of transcriptome-based biomarker signatures of TB-treatment response should include people with diabetes for use across populations.

Acute and Long-Term Consequences of Co-opted doublesex on the Development of Mimetic Butterfly Color Patterns

Novel phenotypes are increasingly recognized to have evolved by co-option of conserved genes into new developmental contexts, yet the process by which co-opted genes modify existing developmental programs remains obscure. Here, we provide insight into this process by characterizing the role of co-opted doublesex in butterfly wing color pattern development. dsx is the master regulator of insect sex differentiation but has been co-opted to control the switch between discrete nonmimetic and mimetic patterns in Papilio alphenor and its relatives through the evolution of novel mimetic alleles. We found dynamic spatial and temporal expression pattern differences between mimetic and nonmimetic butterflies throughout wing development. A mimetic color pattern program is switched on by a pulse of dsx expression in early pupal development that causes acute and long-term differential gene expression, particularly in Wnt and Hedgehog signaling pathways. RNAi suggested opposing, novel roles for these pathways in mimetic pattern development. Importantly, Dsx co-option caused Engrailed, a primary target of Hedgehog signaling, to gain a novel expression domain early in pupal wing development that is propagated through mid-pupal development to specify novel mimetic patterns despite becoming decoupled from Dsx expression itself. Altogether, our findings provide multiple views into how co-opted genes can both cause and elicit changes to conserved networks and pathways to result in development of novel, adaptive phenotypes.

Single-dose VSV-based vaccine protects cynomolgus macaques from disease after Taï Forest virus infection

Taï Forest virus (TAFV) is a lesser-known ebolavirus that causes lethal infections in chimpanzees and is responsible for a single human case. Limited research has been done on this human pathogen; however, with the recent emergence of filoviruses in West Africa, further investigation and countermeasure development against this virus is warranted.We developed a vesicular stomatitis virus (VSV)-based vaccine expressing the TAFV glycoprotein as the viral antigen and assessed it for protective efficacy in nonhuman primates (NHPs). Following a single high-dose vaccination, NHPs developed antigen-specific binding and neutralizing antibodies as well as modest T cell responses. Importantly, all vaccinated NHPs were uniformly protected from disease after lethal TAFV challenge while the naïve control group succumbed to the disease. Histopathologic lesions consistent with filovirus disease were present in control NHPs but were not observed in vaccinated NHPs. Transcriptional analysis of whole blood samples obtained after vaccination and challenge was performed to gain insight into molecular underpinnings conferring protection. Differentially expressed genes (DEG) detected 7 days post-vaccination were enriched to processes associated with innate immunity and antiviral responses. Only a small number of DEG was detected in vaccinated NHPs post-challenge while over 1,000 DEG were detected in control NHPs at end-stage disease which mapped to gene ontology terms indicative of defense responses and inflammation. Taken together, this data demonstrates the effective single-dose protection of the VSV-TAFV vaccine, and its potential for use in outbreaks.

An optogenetic-phosphoproteomic study reveals dynamic Akt1 signaling profiles in endothelial cells

The serine/threonine kinase AKT is a central node in cell signaling. While aberrant AKT activation underlies the development of a variety of human diseases, how different patterns of AKT-dependent phosphorylation dictate downstream signaling and phenotypic outcomes remains largely enigmatic. Herein, we perform a systems-level analysis that integrates methodological advances in optogenetics, mass spectrometry-based phosphoproteomics, and bioinformatics to elucidate how different intensity, duration, and pattern of Akt1 stimulation lead to distinct temporal phosphorylation profiles in vascular endothelial cells. Through the analysis of ~35,000 phosphorylation sites across multiple conditions precisely controlled by light stimulation, we identify a series of signaling circuits activated downstream of Akt1 and interrogate how Akt1 signaling integrates with growth factor signaling in endothelial cells. Furthermore, our results categorize kinase substrates that are preferably activated by oscillating, transient, and sustained Akt1 signals. We validate a list of phosphorylation sites that covaried with Akt1 phosphorylation across experimental conditions as potential Akt1 substrates. Our resulting dataset provides a rich resource for future studies on AKT signaling and dynamics.

Phase 1/2a clinical trial in ALS with ropinirole, a drug candidate identified by iPSC drug discovery

iPSC-based drug discovery led to a phase 1/2a trial of ropinirole in ALS. 20 participants with sporadic ALS received ropinirole or placebo for 24 weeks in the double-blind period to evaluate safety, tolerability, and therapeutic effects. Adverse events were similar in both groups. During the double-blind period, muscle strength and daily activity were maintained, but a decline in the ALSFRS-R, which assesses the functional status of ALS patients, was not different from that in the placebo group. However, in the open-label extension period, the ropinirole group showed significant suppression of ALSFRS-R decline and an additional 27.9 weeks of disease-progression-free survival. iPSC-derived motor neurons from participants showed dopamine D2 receptor expression and a potential involvement of the SREBP2-cholesterol pathway in therapeutic effects. Lipid peroxide represents a clinical surrogate marker to assess disease progression and drug efficacy. Limitations include small sample sizes and high attrition rates in the open-label extension period, requiring further validation.

Spatio-temporal transcriptome dynamics coordinate rapid transition of core crop functions in 'lactating' pigeon

Pigeons (Columba livia) are among a select few avian species that have developed a specialized reproductive mode wherein the parents produce a 'milk' in their crop to feed newborn squabs. Nonetheless, the transcriptomic dynamics and role in the rapid transition of core crop functions during 'lactation' remain largely unexplored. Here, we generated a de novo pigeon genome assembly to construct a high resolution spatio-temporal transcriptomic landscape of the crop epithelium across the entire breeding stage. This multi-omics analysis identified a set of 'lactation'-related genes involved in lipid and protein metabolism, which contribute to the rapid functional transitions in the crop. Analysis of in situ high-throughput chromatin conformation capture (Hi-C) sequencing revealed extensive reorganization of promoter-enhancer interactions linked to the dynamic expression of these 'lactation'-related genes between stages. Moreover, their expression is spatially localized in specific epithelial layers, and can be correlated with phenotypic changes in the crop. These results illustrate the preferential de novo synthesis of 'milk' lipids and proteins in the crop, and provides candidate enhancer loci for further investigation of the regulatory elements controlling pigeon 'lactation'.

Actionable cancer vulnerability due to translational arrest, p53 aggregation and ribosome biogenesis stress evoked by the disulfiram metabolite CuET

Drug repurposing is a versatile strategy to improve current therapies. Disulfiram has long been used in the treatment of alcohol dependency and multiple clinical trials to evaluate its clinical value in oncology are ongoing. We have recently reported that the disulfiram metabolite diethyldithiocarbamate, when combined with copper (CuET), targets the NPL4 adapter of the p97VCP segregase to suppress the growth of a spectrum of cancer cell lines and xenograft models in vivo. CuET induces proteotoxic stress and genotoxic effects, however important issues concerning the full range of the CuET-evoked tumor cell phenotypes, their temporal order, and mechanistic basis have remained largely unexplored. Here, we have addressed these outstanding questions and show that in diverse human cancer cell models, CuET causes a very early translational arrest through the integrated stress response (ISR), later followed by features of nucleolar stress. Furthermore, we report that CuET entraps p53 in NPL4-rich aggregates leading to elevated p53 protein and its functional inhibition, consistent with the possibility of CuET-triggered cell death being p53-independent. Our transcriptomics profiling revealed activation of pro-survival adaptive pathways of ribosomal biogenesis (RiBi) and autophagy upon prolonged exposure to CuET, indicating potential feedback responses to CuET treatment. The latter concept was validated here by simultaneous pharmacological inhibition of RiBi and/or autophagy that further enhanced CuET's tumor cytotoxicity, using both cell culture and zebrafish in vivo preclinical models. Overall, these findings expand the mechanistic repertoire of CuET's anti-cancer activity, inform about the temporal order of responses and identify an unorthodox new mechanism of targeting p53. Our results are discussed in light of cancer-associated endogenous stresses as exploitable tumor vulnerabilities and may inspire future clinical applications of CuET in oncology, including combinatorial treatments and focus on potential advantages of using certain validated drug metabolites, rather than old, approved drugs with their, often complex, metabolic profiles.

Robust classification of wound healing stages in both mice and humans for acute and burn wounds based on transcriptomic data

BACKGROUND

Wound healing involves careful coordination among various cell types carrying out unique or even multifaceted functions. The abstraction of this complex dynamic process into four primary wound stages is essential to the study of wound care for timing treatment and tracking wound progression. For example, a treatment that may promote healing in the inflammatory stage may prove detrimental in the proliferative stage. Additionally, the time scale of individual responses varies widely across and within the same species. Therefore, a robust method to assess wound stages can help advance translational work from animals to humans.

RESULTS

In this work, we present a data-driven model that robustly identifies the dominant wound healing stage using transcriptomic data from biopsies gathered from mouse and human wounds, both burn and surgical. A training dataset composed of publicly available transcriptomic arrays is used to derive 58 shared genes that are commonly differentially expressed. They are divided into 5 clusters based on temporal gene expression dynamics. The clusters represent a 5-dimensional parametric space containing the wound healing trajectory. We then create a mathematical classification algorithm in the 5-dimensional space and demonstrate that it can distinguish between the four stages of wound healing: hemostasis, inflammation, proliferation, and remodeling.

CONCLUSIONS

In this work, we present an algorithm for wound stage detection based on gene expression. This work suggests that there are universal characteristics of gene expression in wound healing stages despite the seeming disparities across species and wounds. Our algorithm performs well for human and mouse wounds of both burn and surgical types. The algorithm has the potential to serve as a diagnostic tool that can advance precision wound care by providing a way of tracking wound healing progression with more accuracy and finer temporal resolution compared to visual indicators. This increases the potential for preventive action.

The gut microbiota promotes distal tissue regeneration via RORγ+ regulatory T cell emissaries

Specific microbial signals induce the differentiation of a distinct pool of RORγ+ regulatory T (Treg) cells crucial for intestinal homeostasis. We discovered highly analogous populations of microbiota-dependent Treg cells that promoted tissue regeneration at extra-gut sites, notably acutely injured skeletal muscle and fatty liver. Inflammatory meditators elicited by tissue damage combined with MHC-class-II-dependent T cell activation to drive the accumulation of gut-derived RORγ+ Treg cells in injured muscle, wherein they regulated the dynamics and tenor of early inflammation and helped balance the proliferation vs. differentiation of local stem cells. Reining in IL-17A-producing T cells was a major mechanism underlying the rheostatic functions of RORγ+ Treg cells in compromised tissues. Our findings highlight the importance of gut-trained Treg cell emissaries in controlling the response to sterile injury of non-mucosal tissues.

Progression-associated molecular changes in basal/squamous and sarcomatoid bladder carcinogenesis

The aggressive basal/squamous (Ba/Sq) bladder cancer (BLCA) subtype is often diagnosed at the muscle-invasive stage and can progress to the sarcomatoid variant. Identification of molecular changes occurring during progression from non-muscle-invasive BLCA (NMIBC) to Ba/Sq muscle-invasive BLCA (MIBC) is thus challenging in human disease. We used the N-butyl-N-(4-hydroxybutyl)-nitrosamine (BBN) mouse model of Ba/Sq MIBC to study longitudinally the molecular changes leading to the Ba/Sq phenotype and to the sarcomatoid variant using IHC and microdissection followed by RNA-seq at all stages of progression. A shift to the Ba/Sq phenotype started in early progression stages. Pathway analysis of gene clusters with coordinated expression changes revealed Shh signaling loss and a shift from fatty acid metabolism to glycolysis. An upregulated cluster, appearing early in carcinogenesis, showed relevance to human disease, identifying NMIBC patients at risk of progression. Similar to the human counterpart, sarcomatoid BBN tumors displayed a Ba/Sq phenotype and epithelial-mesenchymal transition (EMT) features. An EGFR/FGFR1 signaling switch occurred with sarcomatoid dedifferentiation and correlated with EMT. BLCA cell lines with high EMT were the most sensitive to FGFR1 knockout and resistant to EGFR knockout. Taken together, these findings provide insights into the underlying biology of Ba/Sq BLCA progression and sarcomatoid dedifferentiation with potential clinical implications. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

LiveG4ID-Seq for Profiling the Dynamic Landscape of Chromatin G-Quadruplexes During Cell Cycle in Living Cells

G-quadruplex (G4) structures exist in the single-stranded DNA of chromatin and regulate genome function. However, the native chromatin G4 landscape in living cells has yet to be fully characterized. Herein, a genetic-encoded live-cell G4 identifier probe (LiveG4ID) is constructed and its cellular localization, biocompatibility, and G4-binding specificity is evaluated. By coupling LiveG4ID with cleavage under targets and tagmentation (CUT&Tag), LiveG4ID-seq, a method for mapping native chromatin G4 landscape in living cells with high accuracy is established. Compared to the conventional G4 CUT&Tag method, LiveG4ID-seq can identify more chromatin G4 signals and have a higher ratio of true positive signals. Using LiveG4ID-seq, the dynamic landscape of chromatin G4 structures during the cell cycle is profiled. It is discovered that chromatin G4 structures are prevalent in the promoter regions of cell cycle-specific genes, even in the early M phase when the chromatin is condensed. These data demonstrate the capacity of LiveG4ID-seq to profile a more accurate G4 landscape in living cells and promote future studies on chromatin G4 structures.

Integrated transcriptome and proteome analysis reveals posttranscriptional regulation of ribosomal genes in human brain organoids

During development of the human cerebral cortex, multipotent neural progenitors generate excitatory neurons and glial cells. Investigations of the transcriptome and epigenome have revealed important gene regulatory networks underlying this crucial developmental event. However, the posttranscriptional control of gene expression and protein abundance during human corticogenesis remains poorly understood. We addressed this issue by using human telencephalic brain organoids grown using a dual reporter cell line to isolate neural progenitors and neurons and performed cell class and developmental stage-specific transcriptome and proteome analysis. Integrating the two datasets revealed modules of gene expression during human corticogenesis. Investigation of one such module uncovered mTOR-mediated regulation of translation of the 5'TOP element-enriched translation machinery in early progenitor cells. We show that in early progenitors partial inhibition of the translation of ribosomal genes prevents precocious translation of differentiation markers. Overall, our multiomics approach proposes novel posttranscriptional regulatory mechanisms crucial for the fidelity of cortical development.

mRNA translational specialization by RBPMS presets the competence for cardiac commitment in hESCs

The blueprints of developing organs are preset at the early stages of embryogenesis. Transcriptional and epigenetic mechanisms are proposed to preset developmental trajectories. However, we reveal that the competence for the future cardiac fate of human embryonic stem cells (hESCs) is preset in pluripotency by a specialized mRNA translation circuit controlled by RBPMS. RBPMS is recruited to active ribosomes in hESCs to control the translation of essential factors needed for cardiac commitment program, including Wingless/Integrated (WNT) signaling. Consequently, RBPMS loss specifically and severely impedes cardiac mesoderm specification, leading to patterning and morphogenetic defects in human cardiac organoids. Mechanistically, RBPMS specializes mRNA translation, selectively via 3'UTR binding and globally by promoting translation initiation. Accordingly, RBPMS loss causes translation initiation defects highlighted by aberrant retention of the EIF3 complex and depletion of EIF5A from mRNAs, thereby abrogating ribosome recruitment. We demonstrate how future fate trajectories are programmed during embryogenesis by specialized mRNA translation.

Mechanism of pulmonary plague biphasic syndrome: inhibition or activation of NF-κB signaling pathway

Background: Pneumonic plague is a fatal respiratory disease caused by Yersinia pestis. Time-course transcriptome analysis on the mechanism of pneumonic plague biphasic syndrome is lacking in the literature. Materials & methods: This study documented the disease course through bacterial load, histopathology, cytokine levels and flow cytometry. RNA-sequencing technology was used to investigate the global transcriptome profile of lung tissue in mice infected with Y. pestis. Results: Inflammation-related genes were significantly upregulated at 48 h post-infection, while genes related to cell adhesion and cytoskeletal structure were downregulated. Conclusion: NOD-like receptor and TNF signaling pathways play a plausible role in pneumonic plague biphasic syndrome and lung injury by controlling the activation and inhibition of the NF-κB signaling pathway.

Modeling transcriptional regulation of the cell cycle using a novel cybernetic-inspired approach

Quantitative understanding of cellular processes, such as cell cycle and differentiation, is impeded by various forms of complexity ranging from myriad molecular players and their multilevel regulatory interactions, cellular evolution with multiple intermediate stages, lack of elucidation of cause-effect relationships among the many system players, and the computational complexity associated with the profusion of variables and parameters. In this paper, we present an elegant modeling framework based on the cybernetic concept that biological regulation is inspired by objectives embedding entirely novel strategies for dimension reduction, process stage specification through the system dynamics, and innovative causal association of regulatory events with the ability to predict the evolution of the dynamical system. The elementary step of the modeling strategy involves stage-specific objective functions that are computationally-determined from experiments, augmented with dynamical network computations involving end point objective functions, mutual information, change point detection, and maximal clique centrality. We demonstrate the power of the method through application to the mammalian cell cycle, which involves thousands of biomolecules engaged in signaling, transcription, and regulation. Starting with a fine-grained transcriptional description obtained from RNA sequencing measurements, we develop an initial model, which is then dynamically modeled using the cybernetic-inspired method (CIM), utilizing the strategies described above. The CIM is able to distill the most significant interactions from a multitude of possibilities. In addition to capturing the complexity of regulatory processes in a mechanistically causal and stage-specific manner, we identify the functional network modules, including novel cell cycle stages. Our model is able to predict future cell cycles consistent with experimental measurements. We posit that this state-of-the-art framework has the promise to extend to the dynamics of other biological processes, with a potential to provide novel mechanistic insights.

STATEMENT OF SIGNIFICANCE

Cellular processes like cell cycle are overly complex, involving multiple players interacting at multiple levels, and explicit modeling of such systems is challenging. The availability of longitudinal RNA measurements provides an opportunity to "reverse-engineer" for novel regulatory models. We develop a novel framework, inspired using goal-oriented cybernetic model, to implicitly model transcriptional regulation by constraining the system using inferred temporal goals. A preliminary causal network based on information-theory is used as a starting point, and our framework is used to distill the network to temporally-based networks containing essential molecular players. The strength of this approach is its ability to dynamically model the RNA temporal measurements. The approach developed paves the way for inferring regulatory processes in many complex cellular processes.

Developmental dynamic transcriptome and systematic analysis reveal the major genes underlying isoflavone accumulation in soybean

INTRODUCTION

Soy isoflavone, a class of polyphenolic compounds exclusively occurred in legumes, is an important bioactive compound for both plants and human beings. The outline of isoflavones biosynthesis pathway has been drawn up basically in the previous research. However, research on the subject has been mostly restricted to investigate the static regulation of isoflavone content in soybean, rather than characterize its dynamic variation and modulation network in developing seeds.

METHODS

In this study, by using eight recombinant inbred lines (RIL), the contents of six isoflavone components in the different stages of developing soybean seeds were determined to characterize the dynamic variation of isoflavones, and the isoflavones accumulation pattern at physiological level was investigated. Meanwhile, we integrated and analyzed the whole genome expression profile of four lines and 42 meta-transcriptome data, based on the multiple algorithms.

RESULTS

This study: 1) obtained 4 molecular modules strongly correlated with isoflavone accumulation; 2) identified 28 novel major genes that could affect the accumulation of isoflavones in developing seeds free from the limitation of environments; 3) discussed the dynamic molecular patterns regulating isoflavones accumulation in developing seed; 4) expanded the isoflavone biosynthesis pathway.

DISCUSSION

The results not only promote the understandings on the biosynthesis and regulation of isoflavones at physiological and molecular level, but also facilitate to breed elite soybean cultivars with high isoflavone contents.

Gene Co-Expression Network Analysis Reveals the Hub Genes and Key Pathways Associated with Resistance to Salmonella Enteritidis Colonization in Chicken

Salmonella negatively impacts the poultry industry and threatens animals' and humans' health. The gastrointestinal microbiota and its metabolites can modulate the host's physiology and immune system. Recent research demonstrated the role of commensal bacteria and short-chain fatty acids (SCFAs) in developing resistance to Salmonella infection and colonization. However, the complex interactions among chicken, Salmonella, host-microbiome, and microbial metabolites remain unelucidated. Therefore, this study aimed to explore these complex interactions by identifying the driver and hub genes highly correlated with factors that confer resistance to Salmonella. Differential gene expression (DEGs) and dynamic developmental genes (DDGs) analyses and weighted gene co-expression network analysis (WGCNA) were performed using transcriptome data from the cecum of Salmonella Enteritidis-infected chicken at 7 and 21 days after infection. Furthermore, we identified the driver and hub genes associated with important traits such as the heterophil/lymphocyte (H/L) ratio, body weight post-infection, bacterial load, propionate and valerate cecal contents, and Firmicutes, Bacteroidetes, and Proteobacteria cecal relative abundance. Among the multiple genes detected in this study, EXFABP, S100A9/12, CEMIP, FKBP5, MAVS, FAM168B, HESX1, EMC6, and others were found as potential candidate gene and transcript (co-) factors for resistance to Salmonella infection. In addition, we found that the PPAR and oxidative phosphorylation (OXPHOS) metabolic pathways were also involved in the host's immune response/defense against Salmonella colonization at the earlier and later stage post-infection, respectively. This study provides a valuable resource of transcriptome profiles from chicken cecum at the earlier and later stage post-infection and mechanistic understanding of the complex interactions among chicken, Salmonella, host-microbiome, and associated metabolites.

PKI: A bioinformatics method of quantifying the importance of nodes in gene regulatory network via a pseudo knockout index

BACKGROUND

Gene regulatory network (GRN) is a model that characterizes the complex relationships between genes and thereby provides an informatics environment to measure the importance of nodes. The evaluation of important nodes in a GRN can effectively refer to their functional implications severing as key players in particular biological processes, such as master regulator and driver gene. Currently, it is mainly based on network topological parameters and focuses only on evaluating a single node individually. However, genes and products play their functions by interacting with each other. It is worth noting that the effects of gene combinations in GRN are not simply additive. Key combinations discovery is of significance in revealing gene sets with important functions. Recently, with the development of single-cell RNA-sequencing (scRNA-seq) technology, we can quantify gene expression profiles of individual cells that provide the potential to identify crucial nodes in gene regulations regarding specific condition, e.g., stem cell differentiation.

RESULTS

In this paper, we propose a bioinformatics method, called Pseudo Knockout Importance (PKI), to quantify the importance of node and node sets in a specific GRN structure using time-course scRNA-seq data. First, we construct ordinary differential equations to approach the gene regulations during cell differentiation. Then we design gene pseudo knockout experiments and define PKI score evaluation criteria based on the coefficient of determination. The importance of nodes can be described as the influence on the ODE system of removing variables. For key gene combinations, PKI is derived as a combinatorial optimization problem of quantifying the in silico gene knockout effects.

CONCLUSIONS

Here, we focus our analyses on the specific GRN of embryonic stem cells with time series gene expression profile. To verify the effectiveness and advantage of PKI method, we compare its node importance rankings with other twelve kinds of centrality-based methods, such as degree and Latora closeness. For key node combinations, we compare the results with the method based on minimum dominant set. Moreover, the famous combinations of transcription factors in induced pluripotent stem cell are also employed to verify the vital gene combinations identified by PKI. These results demonstrate the reliability and superiority of the proposed method.

Parallelized multidimensional analytic framework applied to mammary epithelial cells uncovers regulatory principles in EMT

A proper understanding of disease etiology will require longitudinal systems-scale reconstruction of the multitiered architecture of eukaryotic signaling. Here we combine state-of-the-art data acquisition platforms and bioinformatics tools to devise PAMAF, a workflow that simultaneously examines twelve omics modalities, i.e., protein abundance from whole-cells, nucleus, exosomes, secretome and membrane; N-glycosylation, phosphorylation; metabolites; mRNA, miRNA; and, in parallel, single-cell transcriptomes. We apply PAMAF in an established in vitro model of TGFβ-induced epithelial to mesenchymal transition (EMT) to quantify >61,000 molecules from 12 omics and 10 timepoints over 12 days. Bioinformatics analysis of this EMT-ExMap resource allowed us to identify; -topological coupling between omics, -four distinct cell states during EMT, -omics-specific kinetic paths, -stage-specific multi-omics characteristics, -distinct regulatory classes of genes, -ligand-receptor mediated intercellular crosstalk by integrating scRNAseq and subcellular proteomics, and -combinatorial drug targets (e.g., Hedgehog signaling and CAMK-II) to inhibit EMT, which we validate using a 3D mammary duct-on-a-chip platform. Overall, this study provides a resource on TGFβ signaling and EMT.

Toward a hypothesis-free understanding of how phosphorylation dynamically impacts protein turnover

The turnover measurement of proteins and proteoforms has been largely facilitated by workflows coupling metabolic labeling with mass spectrometry (MS), including dynamic stable isotope labeling by amino acids in cell culture (dynamic SILAC) or pulsed SILAC (pSILAC). Very recent studies including ours have integrated themeasurement of post-translational modifications (PTMs) at the proteome level (i.e., phosphoproteomics) with pSILAC experiments in steady state systems, exploring the link between PTMs and turnover at the proteome-scale. An open question in the field is how to exactly interpret these complex datasets in a biological perspective. Here, we present a novel pSILAC phosphoproteomic dataset which was obtained during a dynamic process of cell starvation using data-independent acquisition MS (DIA-MS). To provide an unbiased "hypothesis-free" analysis framework, we developed a strategy to interrogate how phosphorylation dynamically impacts protein turnover across the time series data. With this strategy, we discovered a complex relationship between phosphorylation and protein turnover that was previously underexplored. Our results further revealed a link between phosphorylation stoichiometry with the turnover of phosphorylated peptidoforms. Moreover, our results suggested that phosphoproteomic turnover diversity cannot directly explain the abundance regulation of phosphorylation during cell starvation, underscoring the importance of future studies addressing PTM site-resolved protein turnover.

Circadian rhythm disruption is associated with skeletal muscle dysfunction within the blind Mexican Cavefish

Circadian control of physiology and metabolism is pervasive throughout nature, with circadian disruption contributing to premature aging, neurodegenerative disease, and type 2 diabetes (Musiek et al. 2016; Panda, 2016). It has become increasingly clear that peripheral tissues, such as skeletal muscle, possess cell-autonomous clocks crucial for metabolic homeostasis (Gabriel et al. 2021). In fact, disruption of the skeletal muscle circadian rhythm results in insulin resistance, sarcomere disorganization, and muscle weakness in both vertebrates and non-vertebrates - indicating that maintenance of a functional muscle circadian rhythm provides an adaptive advantage. We and others have found that cavefish possess a disrupted central circadian rhythm and, interestingly, a skeletal muscle phenotype strikingly similar to circadian knock-out mutants; namely, muscle loss, muscle weakness, and insulin resistance (Olsen et al. 2022; Riddle et al. 2018; Mack et al. 2021). However, whether the cavefish muscle phenotype results from muscle-specific circadian disruption remains untested. To this point, we investigated genome-wide, circadian-regulated gene expression within the skeletal muscle of the Astyanax mexicanus - comprised of the river-dwelling surface fish and troglobitic cavefish - providing novel insights into the evolutionary consequence of circadian disruption on skeletal muscle physiology.

Disabling Uncompetitive Inhibition of Oncogenic IDH Mutations Drives Acquired Resistance

Mutations in IDH genes occur frequently in acute myeloid leukemia (AML) and other human cancers to generate the oncometabolite R-2HG. Allosteric inhibition of mutant IDH suppresses R-2HG production in a subset of patients with AML; however, acquired resistance emerges as a new challenge, and the underlying mechanisms remain incompletely understood. Here we establish isogenic leukemia cells containing common IDH oncogenic mutations by CRISPR base editing. By mutational scanning of IDH single amino acid variants in base-edited cells, we describe a repertoire of IDH second-site mutations responsible for therapy resistance through disabling uncompetitive enzyme inhibition. Recurrent mutations at NADPH binding sites within IDH heterodimers act in cis or trans to prevent the formation of stable enzyme-inhibitor complexes, restore R-2HG production in the presence of inhibitors, and drive therapy resistance in IDH-mutant AML cells and patients. We therefore uncover a new class of pathogenic mutations and mechanisms for acquired resistance to targeted cancer therapies.

SIGNIFICANCE

Comprehensive scanning of IDH single amino acid variants in base-edited leukemia cells uncovers recurrent mutations conferring resistance to IDH inhibition through disabling NADPH-dependent uncompetitive inhibition. Together with targeted sequencing, structural, and functional studies, we identify a new class of pathogenic mutations and mechanisms for acquired resistance to IDH-targeting cancer therapies. This article is highlighted in the In This Issue feature, p. 1.

MicroRNA 3' ends shorten during adolescent brain maturation

MicroRNA (miRNA) dysregulation is well-documented in psychiatric disease, but miRNA dynamics remain poorly understood during adolescent and early adult brain maturation, when symptoms often first appear. Here, we use RNA sequencing to examine miRNAs and their mRNA targets in cortex and hippocampus from early-, mid-, and late-adolescent and adult mice. Furthermore, we use quantitative proteomics by tandem mass tag mass spectrometry (TMT-MS) to examine protein dynamics in cortex from the same subjects. We found that ~25% of miRNAs' 3' ends shorten with age due to increased 3' trimming and decreased U tailing. Particularly, shorter but functionally competent isoforms (isomiRs) of miR-338-3p increase up to 10-fold during adolescence and only in brain. MiRNAs that undergo 3' shortening exhibit stronger negative correlations with targets that decrease with age and stronger positive correlations with targets that increase with age, than miRNAs with stable 3' ends. Increased 3' shortening with age was also observed in available mouse and human miRNA-seq data sets, and stronger correlations between miRNAs that undergo shortening and their mRNA targets were observed in two of the three available data sets. We conclude that age-associated miRNA 3' shortening is a well-conserved feature of postnatal brain maturation.