Search

Spodoptera litura commonly known as the cutworm, is among the most destructive lepidopteran pests affecting over 120 plants species. The powerful destructive nature of this lepidopteran is attributable to its high reproductive capacity. The testicular fusion that occurs during metamorphosis from larvae to pupa in S.litura positively influences the reproductive success of the offspring. In contrast, Bombyx mori, the silkworm, retains separate testes throughout its life and does not undergo this fusion process. Microscopic examination reveals that during testicular fusion in S.litura, the peritoneal sheath becomes thinner and more translucent, whereas in B.mori, the analogous region thickens. The outer basement membrane in S.litura exhibits fractures, discontinuity, and uneven thickness accompanied by a significant presence of cellular secretions, large cell size, increased vesicles, liquid droplets, and a proliferation of rough endoplasmic reticulum and mitochondria. In contrast, the testicular peritoneal sheath of B.mori at comparable developmental stage exhibits minimal change. Comparative transcriptomic analysis of the testicular peritoneal sheath reveals a substantial difference in gene expression between the two species. The disparity in differential expressed genes (DEGs) is linked to an enrichment of numerous transcription factors, intracellular signaling pathways involving Ca 2+ and GTPase, as well as intracellular protein transport and signaling pathways. Meanwhile, structural proteins including actin, chitin-binding proteins, membrane protein fractions, cell adhesion, extracellular matrix proteins are predominantly identified. Moreover, the study highlights the enrichment of endopeptidases, serine proteases, proteolytic enzymes and matrix metalloproteins, which may play a role in the degradation of the outer membrane. Five transcription factors-Slforkhead, Slproline, Slcyclic, Slsilk, and SlD-ETS were identified, and their expression pattern were confirmed by qRT-PCR. they are candidates for participating in the regulation of testicular fusion. Our findings underscore significant morphological and trancriptomic variation in the testicular peritoneal sheath of S.litura compared to the silkworm, with substantial changes at the transcriptomic level coinciding with testicular fusion. The research provides valuable clues for understanding the complex mechanisms underlying this unique phenomenon in insects., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Background: Leukocyte- and platelet-rich fibrin (L-PRF) is an autologous platelet concentrate, prepared by centrifugation of blood and consisting of a dense fibrin network with incorporated leukocytes and platelets. This study aims to perform an in-depth analysis of the cells, growth factors, and transcriptome of L-PRF., Methods and Results: Fresh, 1 week and 2 weeks cultured human L-PRF membranes and liquid L-PRF glue were characterized on cellular and transcriptional level using flow cytometry (n = 4), single-cell RNA sequencing (n = 5) and RT-qPCR. Growth factor kinetics were investigated using ELISA (EGF, VEGF, PDGF-AB, TGF-β1, bFGF). L-PRF contained a large number of viable cells (fresh 97.14 ± 1.09%, 1 week cultured 93.57 ± 1.68%), mainly granulocytes in fresh samples (53.9 ± 19.86%) and T cells in cultured samples (84.7 ± 6.1%), confirmed with scRNA-seq. Monocytes differentiate to macrophages during 1 week incubation. Specifically arterial L-PRF membranes were found to release significant amounts of VEGF, EGF, PDGF-AB and TGF-β1., Conclusion: We characterized L-PRF using in vitro experiments, to obtain an insight in the composition of the material including a possible mechanistic role for tissue healing. This was the first study characterizing L-PRF at a combined cellular, proteomic, and transcriptional level., (© 2024. The Author(s).)

The domestication of crops, coupled with agroecosystem development, is associated with major environmental changes and provides an ideal model of phenotypic plasticity. Here, we examined 32 genotypes of three tetraploid wheat (Triticum turgidum L.) subspecies, wild emmer, emmer, and durum wheat, which are representative of the key stages in the domestication of tetraploid wheat. We developed a pipeline that integrates RNA-Seq data and population genomics to assess gene expression plasticity and identify selection signatures under diverse nitrogen availability conditions. Our analysis revealed differing gene expression responses to nitrogen availability across primary (wild emmer to emmer) and secondary (emmer to durum wheat) domestication. Notably, nitrogen triggered the expression of twice as many genes in durum wheat compared to that in emmer and wild emmer. Unique selection signatures were identified at each stage: primary domestication mainly influenced genes related to biotic interactions, whereas secondary domestication affected genes related to amino acid metabolism, in particular lysine. Selection signatures were found in differentially expressed genes (DEGs), notably those associated with nitrogen metabolism, such as the gene encoding glutamate dehydrogenase (GDH). Overall, our study highlights the pivotal role of nitrogen availability in the domestication and adaptive responses of a major food crop, with varying effects across different traits and growth conditions., Competing Interests: Conflict of interest statement. None declared., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists.)

Decapping is a crucial step in mRNA degradation in eucaryotes and requires the formation of a holoenzyme complex between the decapping enzyme DECAPPING 2 (DCP2) and the decapping enhancer DCP1. In Arabidopsis (Arabidopsis thaliana), DCP1-ASSOCIATED NYN ENDORIBONUCLEASE 1 (DNE1) is a direct protein partner of DCP1. The function of both DNE1 and decapping is necessary to maintain phyllotaxis, the regularity of organ emergence in the apex. In this study, we combined in vivo mRNA editing, RNA degradome sequencing, transcriptomics, and small RNA-omics to identify targets of DNE1 and study how DNE1 and DCP2 cooperate in controlling mRNA fate. Our data reveal that DNE1 mainly contacts and cleaves mRNAs in the coding sequence and has sequence cleavage preferences. DNE1 targets are also degraded through decapping, and both RNA degradation pathways influence the production of mRNA-derived small interfering RNAs. Finally, we detected mRNA features enriched in DNE1 targets including RNA G-quadruplexes and translated upstream open reading frames. Combining these four complementary high-throughput sequencing strategies greatly expands the range of DNE1 targets and allowed us to build a conceptual framework describing the influence of DNE1 and decapping on mRNA fate. These data will be crucial to unveil the specificity of DNE1 action and understand its importance for developmental patterning., Competing Interests: Conflict of interest statement. None declared., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

Understanding plant responses to individual stresses does not mean that we understand real-world situations, where stresses usually combine and interact. These interactions arise at different levels, from stress exposure to the molecular networks of the stress response. Here, we built an in-depth multiomic description of plant responses to mild water (W) and nitrogen (N) limitations, either individually or combined, among 5 genetically different Arabidopsis (Arabidopsis thaliana) accessions. We highlight the different dynamics in stress response through integrative traits such as rosette growth and the physiological status of the plants. We also used transcriptomic and metabolomic profiling during a stage when the plant response was stabilized to determine the wide diversity in stress-induced changes among accessions, highlighting the limited reality of a "universal" stress response. The main effect of the W × N interaction was an attenuation of the N-deficiency syndrome when combined with mild drought, but to a variable extent depending on the accession. Other traits subject to W × N interactions are often accession specific. Multiomic analyses identified a subset of transcript-metabolite clusters that are critical to stress responses but essentially variable according to the genotype factor. Including intraspecific diversity in our descriptions of plant stress response places our findings in perspective., Competing Interests: Conflict of interest statement. None declared., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

Background: NUP98 rearrangements (NUP98-r) are rare but overrepresented mutations in pediatric acute myeloid leukemia (AML) patients. NUP98-r is often associated with chemotherapy resistance and a particularly poor prognosis. Therefore, characterizing pediatric AML with NUP98-r to identify aberrations is critically important., Methods: Here, we retrospectively analyzed the clinicopathological features, genomic and transcriptomic landscapes, treatments, and outcomes of pediatric patients with AML., Results: Nine patients with NUP98-r mutations were identified in our cohort of 142 patients. Ten mutated genes were detected in patients with NUP98-r. The frequency of FLT3-ITD mutations differed significantly between the groups harboring NUP98-r and those without NUP98-r (P = 0.035). Unsupervised hierarchical clustering via RNA sequencing data from 21 AML patients revealed that NUP98-r samples clustered together, strongly suggesting a distinct subtype. Compared with that in the non-NUP98-r fusion and no fusion groups, CMAHP expression was significantly upregulated in the NUP98-r samples (P < 0.001 and P = 0.001, respectively). Multivariate Cox regression analyses demonstrated that patients harboring NUP98-r (P < 0.001) and WT1 mutations (P = 0.030) had worse relapse-free survival, and patients harboring NUP98-r (P < 0.008) presented lower overall survival., Conclusions: These investigations contribute to the understanding of the molecular characteristics, risk stratification, and prognostic evaluation of pediatric AML patients., (© 2024. The Author(s).)

Trehalose 6-phosphate (Tre6P) is an essential signal metabolite that regulates the level of sucrose, linking growth and development to the metabolic status. We hypothesized that Tre6P plays a role in mediating the regulation of gene expression by sucrose. To test this, we performed transcriptomic profiling on Arabidopsis (Arabidopsis thaliana) plants that expressed a bacterial TREHALOSE 6-PHOSPHATE SYNTHASE (TPS) under the control of an ethanol-inducible promoter. Induction led to a 4-fold rise in Tre6P levels, a concomitant decrease in sucrose, significant changes (FDR ≤ 0.05) of over 13,000 transcripts, and 2-fold or larger changes of over 5,000 transcripts. Comparison with nine published responses to sugar availability allowed some of these changes to be linked to the rise in Tre6P, while others were probably due to lower sucrose or other indirect effects. Changes linked to Tre6P included repression of photosynthesis-related gene expression and induction of many growth-related processes including ribosome biogenesis. About 500 starvation-related genes are known to be induced by SUCROSE-NON-FERMENTING-1-RELATED KINASE 1 (SnRK1). They were largely repressed by Tre6P in a manner consistent with SnRK1 inhibition by Tre6P. SnRK1 also represses many genes that are involved in biosynthesis and growth. These responded to Tre6P in a more complex manner, pointing toward Tre6P interacting with other C-signaling pathways. Additionally, elevated Tre6P modified the expression of genes encoding regulatory subunits of the SnRK1 complex and TPS class II and FCS-LIKE ZINC FINGER proteins that are thought to modulate SnRK1 function and genes involved in circadian, TARGET OF RAPAMYCIN, light, abscisic acid, and other hormone signaling., Competing Interests: Conflict of interest statement. None declared., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists.)

Citrus is one of the most important fruit crop genera in the world, but many Citrus species are vulnerable to cold stress. Ichang papeda (Citrus ichangensis), a cold-hardy citrus species, holds great potential for identifying valuable metabolites that are critical for cold tolerance in Citrus. However, the metabolic changes and underlying mechanisms that regulate Ichang papeda cold tolerance remain largely unknown. In this study, we compared the metabolomes and transcriptomes of Ichang papeda and HB pummelo (Citrus grandis "Hirado Buntan", a cold-sensitive species) to explore the critical metabolites and genes responsible for cold tolerance. Metabolomic analyses led to the identification of common and genotype-specific metabolites, consistent with transcriptomic alterations. Compared to HB pummelo under cold stress, Ichang papeda accumulated more sugars, flavonoids, and unsaturated fatty acids, which are well-characterized metabolites involved in stress responses. Interestingly, sphingosine and chlorogenic acid substantially accumulated only in Ichang papeda. Knockdown of CiSPT (C. ichangensis serine palmitoyltransferase) and CiHCT2 (C. ichangensis hydroxycinnamoyl-CoA: shikimate hydroxycinnamoyltransferase2), two genes involved in sphingosine and chlorogenic acid biosynthesis, dramatically decreased endogenous sphingosine and chlorogenic acid levels, respectively. This reduction in sphingosine and chlorogenic acid notably compromised the cold tolerance of Ichang papeda, whereas exogenous application of these metabolites increased plant cold tolerance. Taken together, our findings indicate that greater accumulation of a spectrum of metabolites, particularly sphingosine and chlorogenic acid, promotes cold tolerance in cold-tolerant citrus species. These findings broaden our understanding of plant metabolic alterations in response to cold stress and provide valuable targets that can be manipulated to improve Citrus cold tolerance., Competing Interests: Conflict of interest statement. The authors declare no conflicts of interest associated with this work., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

As a commonly used medicinal plant, the flavonoid metabolites of Blumea balsamifera and their association with genes are still elusive. In this study, the total flavonoid content (TFC), flavonoid metabolites and biosynthetic gene expression patterns of B. balsamifera after application of exogenous methyl jasmonate (MeJA) were scrutinized. The different concentrations of exogenous MeJA increased the TFC of B. balsamifera leaves after 48 h of exposure, and there was a positive correlation between TFC and the elicitor concentration. A total of 48 flavonoid metabolites, falling into 10 structural classes, were identified, among which flavones and flavanones were predominant. After screening candidate genes by transcriptome mining, the comprehensive analysis of gene expression level and TFC suggested that FLS and MYB may be key genes that regulate the TFC in B. balsamifera leaves under exogenous MeJA treatment. This study lays a foundation for elucidating flavonoids of B. balsamifera, and navigates the breeding of flavonoid-rich B. balsamifera varieties., (© 2024 Scandinavian Plant Physiology Society.)

Summary: Measurement of single-cell gene expression at different timepoints enables the study of cell development. However, due to the resource constraints and technical challenges associated with the single-cell experiments, researchers can only profile gene expression at discrete and sparsely sampled timepoints. This missing timepoint information impedes downstream cell developmental analyses. We propose scNODE, an end-to-end deep learning model that can predict in silico single-cell gene expression at unobserved timepoints. scNODE integrates a variational autoencoder with neural ordinary differential equations to predict gene expression using a continuous and nonlinear latent space. Importantly, we incorporate a dynamic regularization term to learn a latent space that is robust against distribution shifts when predicting single-cell gene expression at unobserved timepoints. Our evaluations on three real-world scRNA-seq datasets show that scNODE achieves higher predictive performance than state-of-the-art methods. We further demonstrate that scNODE's predictions help cell trajectory inference under the missing timepoint paradigm and the learned latent space is useful for in silico perturbation analysis of relevant genes along a developmental cell path., Availability and Implementation: The data and code are publicly available at https://github.com/rsinghlab/scNODE., (© The Author(s) 2024. Published by Oxford University Press.)

Angelica sinensis, a traditional Chinese medicinal plant, has been primarily reported due to its nutritional value. Pigmentation in this plant is an important appearance trait that directly affects its commercial value. To understand the mechanism controlling purpleness in A. sinensis, hormonal and transcriptomic analyses were performed in three different tissues (leave, root and stem), using two cultivars with contrasting colors. The two-dimensional data set provides dynamic hormonal and gene expression networks underpinning purpleness in A. sinensis. We found abscisic acid as a crucial hormone modulating anthocyanin biosynthesis in A. sinensis. We further identified and validated 7 key genes involved in the anthocyanin biosynthesis pathway and found a specific module containing ANS as a hub gene in WGCNA. Overexpression of a candidate pigment regulatory gene, AsANS (AS08G02092), in transgenic calli of A. sinensis resulted in increased anthocyanin production and caused purpleness. Together, these analyses provide an important understanding of the molecular networks underlying A. sinensis anthocyanin production and its correlation with plant hormones, which can provide an important source for breeding., (© 2024 Scandinavian Plant Physiology Society.)

Background and Research Purpose: Paeoniflorin and albiflorin are monoterpene glycosides that exhibit various medicinal properties in Paeonia species. This study explored the terpene biosynthesis pathway and analyzed the distribution of these compounds in different tissues of two Korean landraces of Paeonia lactiflora to gain insights into the biosynthesis of monoterpene glycosides in P. lactiflora and their potential applications., Materials and Methods: Two Korean landraces, Hongcheon var. and Hwacheon var, of P. lactiflora were used for the analyses. Contents of the paeoniflorin and albiflorin were analyzed using HPLC. RNA was extracted, sequenced, and subjected to transcriptome analysis. Differential gene expression, KEGG, and GO analyses were performed. Paeoniflorin biosynthesis genes were isolated from the transcriptomes using the genes in Euphorbia maculata with the NBLAST program. Phylogenetic analysis of of 1-Deoxy-D-xylulose 5-phosphate synthase (DOXPS), geranyl pyrophosphate synthase (GPPS), and pinene synthase (PS) was carried out with ClustalW and MEGA v5.0., Results and Discussion: Analysis of paeoniflorin and albiflorin content in different tissues of the two P. lactiflora landraces revealed significant variation. Transcriptome analysis yielded 36,602 unigenes, most of which were involved in metabolic processes. The DEG analysis revealed tissue-specific expression patterns with correlations between landraces. The isolation of biosynthetic genes identified 173 candidates. Phylogenetic analysis of the key enzymes in these pathways provides insights into their evolutionary relationships. The sequencing and analysis of DOXPS, GPPS, PS revealed distinct clades and subclades, highlighting their evolutionary divergence and functional conservation. Our findings highlight the roots as the primary sites of paeoniflorin and albiflorin accumulation in P. lactiflora, underscoring the importance of tissue-specific gene expression in their biosynthesis., Conclusion: this study advances our understanding of monoterpene glycoside production and distribution in Paeonia, thereby guiding further plant biochemistry investigations., (© 2024. The Author(s) under exclusive licence to The Genetics Society of Korea.)

Currently, fructooligosaccharides (FOS) are converted from sucrose by purified enzymes or fungal cells, but these methods are costly and time-consuming. Here, the optimal fermentation conditions for strain E326 were determined through fermentation optimization: initial glucose 200 g/L, NaCl 25 g/L, inoculum volume 20 %, dissolved oxygen 20-30 %, pH 3, and glucose feeding concentration 100 g/L, which increased erythritol titer by 1.5 times. The co-expression of HGT1 and APC11 genes alleviated the erythritol synthesis stagnation, shorten the fermentation time by 16.7 %, and increased the erythritol productivity by 17.2 %. The episomal plasmids based on yeast mitochondrial replication origins (mtORIs) were constructed to surface display fructosyltransferase, effectively utilizing waste yeast cells generated during erythritol fermentation. Under the conditions of 60℃ and pH 6, the FOS yield reached 65 %, which to our best of knowledge is so-far the highest yield of FOS obtained. These findings will contribute to the industrial production of erythritol and FOS., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: [Zhiming Rao reports financial support was provided by National Key Research and Development Program of China. Zhiming Rao reports financial support was provided by National Natural Science Foundation of China. Zhiming Rao reports financial support was provided by National First-class Discipline Program of Light Industry Technology and Engineering. Minglong Shao reports financial support was provided by Open Funding Project of Key Laboratory of Industrial Biotechnology. Minglong Shao reports financial support was provided by Fundamental Research Funds for the Central Universities. If there are other authors, they declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.]., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Aging is a gradual, inevitable physiologic process. The organ aging is related to the persistence of chronic inflammation, but the understanding of inflammatory state during renal aging is lacking currently. Single-cell transcriptome sequencing was performed on aging mouse kidney to reveal the molecular phenotype and composition changes of different cell types. In the early stage of aging, immune cells such as T, B cells and mononuclear macrophages increased in kidney. The molecular state of T cells in aging kidney changed and polarized. Among them, we identified a group of GZMK + CD8 + T cells with high expression of Eomes, Pdcd1 and Ifng and a group of Il17a + T cells with high expression of Il17a and Il23r. Moreover, the cytokines and inflammations can aggravate tissue damage eventually. Furthermore, we found the interaction between different types of epithelial cells and T cells increased during the renal aging. These results identify the changes of T cells in the early stage of aging kidney and suggest that GZMK + CD8 + T cells might be a potential target to ameliorate age-associated dysfunctions of kidney(Graphical Abstract)., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Lung neuroendocrine neoplasms (NENs) are a diverse group of tumors characterized by neuroendocrine (NE) differentiation. Among lung NENs, lung carcinoid (LC) is a rare tumor with unique characteristics. Recent research has highlighted the importance of transcription factors (TFs) in establishing gene expression programs in lung NENs such as small cell lung carcinoma. However, the TFs that control the gene expression of LC are largely unknown. In this study, we report the expression and potential function of a TF called Prospero homeobox protein1 (PROX1) in LC. Publicly available transcriptome data suggested that PROX1 was highly expressed in LC tissues, which was confirmed by immunohistochemical analysis on a tissue microarray. Knockdown of PROX1 did not impact the cellular viability of an LC-derived cell line, NCI-H727. Meanwhile, transcriptome analysis revealed that PROX1 knockdown altered the expression of genes involved in NE differentiation. ASCL1, CHGA, CALCA, and LINC00261 were suggested as downstream genes of PROX1. These findings indicate that PROX1 may play an important role in the NE identity of LC by regulating the expression of key target genes., (© 2024. The Author(s) under exclusive licence to Japan Human Cell Society.)

Vicia villosa (VV) and Vicia sativa (VS) are legume forages highly valued for their excellent nitrogen fixation. However, no research has addressed the mechanisms underlying their differences in nitrogen fixation. This study employed physiological, cytological, and comparative transcriptomic approaches to elucidate the disparities in nitrogen fixation between them. Our results showed that the total amount of nitrogen fixed was 60.45% greater in VV than in VS, and the comprehensive nitrogen response performance was 94.19% greater, while the nitrogen fixation efficiency was the same. The infection zone and differentiated bacteroid proportion in mature VV root nodules were 33.76% and 19.35% greater, respectively, than those in VS. The size of the VV genome was 15.16% larger than that of the VS genome, consistent with its greater biomass. A significant enrichment of the flavonoid biosynthetic pathway was found only for VV-specific genes, among which chalcone-flavonone isomerase, caffeoyl-CoA-O-methyltransferase and stilbene synthase were extremely highly expressed. The VV-specific genes also exhibited significant enrichment in symbiotic nodulation; genes related to nodule-specific cysteine-rich peptides (NCRs) comprised 61.11% of the highly expressed genes. qRT‒PCR demonstrated that greater enrichment and expression of the dominant NCR (Unigene0004451) were associated with greater nodule bacteroid differentiation and greater nitrogen fixation in VV. Our findings suggest that the greater total nitrogen fixation of VV was attributed to its larger biomass, leading to a greater nitrogen demand and enhanced fixation physiology. This process is likely achieved by the synergistic effects of high bacteroid differentiation along with high expression of flavonoid and NCR genes., Competing Interests: Declaration of competing interest Authors declare there is no conflict of interest., (Copyright © 2024 Elsevier Masson SAS. All rights reserved.)

Spatholobus suberectus Dunn (Leguminosae) has been used for medicinal purposes for a long period. Flavonoids are the major bioactive components of S. suberectus. However, there is still limited knowledge of the exact method via which transcription factors (TFs) regulate flavonoid biosynthesis. The full-length transcriptome of S. suberectus was analyzed using SMRT sequencing; 61,548 transcripts were identified, including 12,311 new gene loci, 53,336 novel transcripts, 44,636 simple sequence repeats, 36,414 complete coding sequences, 871 long non-coding RNAs and 6781 TFs. The SsMYB158 TF, which is associated with flavonoid biosynthesis, belongs to the R2R3-MYB class and is localized subcellularly to the nucleus. The overexpression of SsMYB158 in Nicotiana benthamiana and the transient overexpression of SsMYB158 in S. suberectus resulted in a substantial enhancement in both flavonoids and catechin levels. In addition, there was a remarkable upregulation in the expression of essential enzyme-coding genes associated with the flavonoid biosynthesis pathways. Our study revealed SsMYB158 as a critical regulator of flavonoid biosynthesis in S. suberectus and laying the foundation for its molecular breeding., Competing Interests: Declaration of competing interest The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Elsevier Masson SAS. All rights reserved.)

Anesthesia serves as a pivotal tool in modern medicine, creating a transient state of sensory deprivation to ensure a pain-free surgical or medical intervention. While proficient in alleviating pain, anesthesia significantly modulates brain dynamics, metabolic processes, and neural signaling, thereby impairing typical cognitive functions. Furthermore, anesthesia can induce notable impacts such as memory impairment, decreased cognitive function, and diminished intelligence, emphasizing the imperative need to explore the concealed repercussions of anesthesia on individuals. In this investigation, we aggregated gene expression profiles (GSE64617, GSE141242, GSE161322, GSE175894, and GSE178995) from public repositories following second-generation sequencing analysis of various anesthetics. Through scrutinizing post-anesthesia brain tissue gene expression utilizing Gene Set Enrichment Analysis (GSEA), Robust Rank Aggregation (RRA), and Weighted Gene Co-expression Network Analysis (WGCNA), this research aims to pinpoint pivotal genes, pathways, and regulatory networks linked to anesthesia. This undertaking not only enhances comprehension of the physiological changes brought about by anesthesia but also lays the groundwork for future investigations, cultivating new insights and innovative perspectives in medical practice., Competing Interests: Declaration of competing interest None Declared., (Copyright © 2024. Published by Elsevier Ltd.)

Dairy cattle with high (HM) versus low muscle (LM) reserves as determined by longissimus dorsi muscle depth (LDD) in late gestation exhibit differential muscle mobilization related to subsequent milk production. Moreover, branched-chain volatile fatty acid (BCVFA) supplementation increased blood glucose levels. We hypothesized that differences in HM and LM reflect distinct muscle metabolism and that BCVFA supplementation altered metabolic pathways. At 42 days before expected calving (BEC), Holstein dairy cows were enrolled in a 2 × 2 factorial study of diet and muscle reserves, by assignment to control (CON)- or BCVFA-supplemented diets and LDD of HM (>4.6 cm) or LM (≤4.6 cm) groups: HM-CON ( n = 13), HM-BCVFA ( n = 10), LM-CON ( n = 9), and LM-BCVFA ( n = 9). Longisumus dorsi muscle was biopsied at 21 days BEC, total RNA was isolated, and protein-coding gene expression was measured with RNA sequencing. Between HM and LM, 713 genes were differentially expressed and 481 between BCVFA and CON ( P < 0.05). Transcriptional signatures indicated differential distribution of type II fibers between groups, with MYH1 greater in LM cattle and MYH2 greater in HM cattle ( P < 0.05). Signatures of LM cattle relative to HM cattle indicated greater activation of autophagy, ubiquitin-proteasome, and Ca 2+ -calpain pathways. HM cattle displayed greater expression of genes that encode extracellular matrix proteins and factors that regulate their proteolysis and turnover. BCVFA modified transcriptomes by increasing expression of genes that regulate fatty acid degradation and flux of carbons into the tricarboxylic acid cycle as acetyl CoA. Molecular signatures support distinct metabolic strategies between LM and HM cattle and that BCVFA supplementation increased substrates for energy generation. NEW & NOTEWORTHY Muscle biopsies of the longissimus dorsi of prepartum dairy cattle indicate that molecular signatures support distinct metabolic strategies between low- and high-muscle cattle and that branched-chain volatile fatty acid supplementation increased substrates for energy generation.

Dendrobium loddigesii, a member of the Orchidaceae family, is a valuable horticultural crop known for its aromatic qualities. However, the mechanisms responsible for the development of its aromatic characteristics remain poorly understood. To elucidate these underlying mechanisms, we assembled the first chromosome-level reference genome of D. loddigesii using PacBio HiFi-reads, Illumina short-reads, and Hi-C data. The assembly comprises 19 pseudochromosomes with N50 contig and N50 scaffold sizes of 55.15 and 89.94 Mb, respectively, estimating the genome size to be 1.68 Gb, larger than that of other sequenced Dendrobium species. During the flowering stages, we conducted a comprehensive analysis combining volatilomics and transcriptomics to understand the characteristics and biosynthetic mechanisms pathways of the floral scent. Our findings emphasize the significant contribution of aromatic terpenoids, especially monoterpenoids, in defining the floral aroma. Furthermore, we identified two crucial terpene synthase (TPS) genes that play a key role in maintaining the aroma during flowering. Through the integration volatilomics data with catalytic assays of DlTPSbs proteins, we identified specific compounds responsible for the aromatic characteristics of D. loddigesii. This integrated analysis of the genome, transcriptome, and volatilome, offers valuable insights into the development and preservation of D. loddigesii's aromatic characteristics, setting the stage for further exploration of the botanical perfumer hypothesis., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Masson SAS. All rights reserved.)

Erythroid cells, the most prevalent cell type in blood, are one of the earliest products and permeate through the entire process of hematopoietic development in the human body, the oxygen-transporting function of which is crucial for maintaining overall health and life support. Previous investigations into erythrocyte differentiation and development have primarily focused on population-level analyses, lacking the single-cell perspective essential for comprehending the intricate pathways of erythroid maturation, differentiation, and the encompassing cellular heterogeneity. The continuous optimization of single-cell transcriptome sequencing technology, or single-cell RNA sequencing (scRNA-seq), provides a powerful tool for life sciences research, which has a particular superiority in the identification of unprecedented cell subgroups, the analyzing of cellular heterogeneity, and the transcriptomic characteristics of individual cells. Over the past decade, remarkable strides have been taken in the realm of single-cell RNA sequencing technology, profoundly enhancing our understanding of erythroid cells. In this review, we systematically summarize the recent developments in single-cell transcriptome sequencing technology and emphasize their substantial impact on the study of erythroid cells, highlighting their contributions, including the exploration of functional heterogeneity within erythroid populations, the identification of novel erythrocyte subgroups, the tracking of different erythroid lineages, and the unveiling of mechanisms governing erythroid fate decisions. These findings not only invigorate erythroid cell research but also offer new perspectives on the management of diseases related to erythroid cells., (© 2024. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Background: A growing body of research has shown that patients with coronavirus disease 2019 (COVID-19) have significantly higher rates of venous thromboembolism (VTE) than healthy. However, the mechanism remains incompletely elucidated. This study aimed to further investigate the molecular mechanisms underlying the development of this complication., Methods: The gene expression profiles of COVID-19 and VTE were downloaded from the Gene Expression Omnibus (GEO) database. After identifying the common differentially expressed genes (DEGs) for COVID-19 and VTE, functional annotation, a protein-protein interactions (PPI) network, module construction, and hub gene identification were performed. Finally, we constructed a transcription factor (TF)-gene regulatory network and a TF-miRNA regulatory network for hub genes., Results: A total of 42 common DEGs were selected for subsequent analyses. Functional analyses showed that biological function and signaling pathways collectively participated in the development and progression of VTE and COVID-19. Finally, 8 significant hub genes were identified using the cytoHubba plugin, including RSL24D1, RPS17, RPS27, HINT1, COX7C, RPL35, RPL34, and NDUFA4, which had preferable values as diagnostic markers for COVID-19 and VTE., Conclusions: Our study revealed the common pathogenesis of COVID-19 and VTE. These common pathways and pivotal genes may provide new ideas for further mechanistic studies., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Mammalian cells are commonly used as hosts in cell culture for biologics production in the pharmaceutical industry. Structured mechanistic models of metabolism have been used to capture complex cellular mechanisms that contribute to varying metabolic shifts in different cell lines. However, little research has focused on the impact of temporal changes in enzyme abundance and activity on the modeling of cell metabolism. In this work, we present a framework for constructing mechanistic models of metabolism that integrate growth-signaling control of enzyme activity and transcript dynamics. The proposed approach is applied to build models for three Chinese hamster ovary (CHO) cell lines using fed-batch culture data and time-series transcript profiles. Leveraging information from the transcriptome data, we develop a parameter estimation approach based on multi-cell-line (MCL) learning, which combines data sets from different cell lines and trains the individual cell-line models jointly to improve model accuracy. The computational results demonstrate the important role of growth signaling and transcript variability in metabolic models as well as the virtue of the MCL approach for constructing cell-line models with a limited amount of data. The resulting models exhibit a high level of accuracy in predicting distinct metabolic behaviors in the different cell lines; these models can potentially be used to accelerate the process and cell-line development for the biomanufacturing of new protein therapeutics., (© 2024 The Author(s). Biotechnology and Bioengineering published by Wiley Periodicals LLC.)

Human embryonic stem cells (hESCs) derived from blastocyst stage embryos present a primed state of pluripotency, whereas mouse ESCs (mESCs) display naïve pluripotency. Their unique characteristics make naïve hESCs more suitable for particular applications in biomedical research. This work aimed to derive hESCs from single blastomeres and determine their pluripotency state, which is currently unclear. We derived hESC lines from single blastomeres of 8-cell embryos and from whole blastocysts, and analysed several naïve pluripotency indicators, their transcriptomic profile and their trilineage differentiation potential. No significant differences were observed between blastomere-derived hESCs (bm-hESCs) and blastocyst-derived hESCs (bc-hESCs) for most naïve pluripotency indicators, including TFE3 localization, mitochondrial activity, and global DNA methylation and hydroxymethylation, nor for their trilineage differentiation potential. Nevertheless, bm-hESCs showed an increased single-cell clonogenicity and a higher expression of naïve pluripotency markers at early passages than bc-hESCs. Furthermore, RNA-seq revealed that bc-hESCs overexpressed a set of genes related to the post-implantational epiblast. Altogether, these results suggest that bm-hESCs, although displaying primed pluripotency, would be slightly closer to the naïve end of the pluripotency continuum than bc-hESCs., Competing Interests: Declaration of competing interest No competing interests declared., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

Recently, numerous studies including various tissues have been carried out on long non-coding RNAs (lncRNAs), but still, its variability has not yet been fully understood. In this study, we characterised the inter-individual variability of lncRNAs in pigs, in the context of number, length and expression. Transcriptomes collected from muscle tissue belonging to six Polish Landrace boars (PL1-PL6), including half-brothers (PL1-PL3), were investigated using bioinformatics (lncRNA identification and functional analysis) and statistical (lncRNA variability) methods. The number of lncRNA ranged from 1289 to 3500 per animal, and the total number of common lncRNAs among all boars was 232. The number, length and expression of lncRNAs significantly varied between individuals, and no consistent pattern has been found between pairs of half-brothers. In detail, PL5 exhibits lower expression than the others, while PL4 has significantly higher expression than PL2-PL3 and PL5-PL6. Noteworthy, comparing the inter-individual variability of lncRNA and mRNA expression, they exhibited concordant patterns. The enrichment analysis for common lncRNA target genes determined a variety of biological processes that play fundamental roles in cell biology, and they were mostly related to whole-body homeostasis maintenance, energy and protein synthesis as well as dynamics of multiple nucleoprotein complexes. The high variability of lncRNA landscape in the porcine genome has been revealed in this study. The inter-individual differences have been found in the context of three aspects: the number, length and expression of lncRNAs, which contribute to a better understanding of its complex nature., (© 2024. The Author(s), under exclusive licence to Institute of Plant Genetics Polish Academy of Sciences.)

Phenotypical innovations during evolution are caused by novel mutations, which are usually heterozygous at the beginning. The gene expressions on two alleles of these mutation sites are not necessarily identical, leading to flexible allele-specific regulation in cell systems. We retrieve the transcriptome data of normal and non-small cell lung cancer (NSCLC) tissues from 47 African Americans (AA) and 50 European Americans (EA). We analyze the differentially expressed genes (DEGs) in NSCLC as well as the tumor-specific mutations. Expression and mutation profiles show convergent evolution in AA and EA populations. The tumor-specific mutations are poorly overlapped, but many of them are located in the same genes, mainly oncogenes and tumor suppressor genes. The DEGs in tumors are majorly caused by the mutated alleles rather than normal alleles. The relative expressions of mutated alleles are highly correlated between AA and EA. The differential expression in NSCLC is predominantly mediated by the mutated alleles on heterozygous sites. This molecular mechanism underlying NSCLC oncogenesis is conserved across different human populations, exhibiting convergent evolution. We present this novel angle that differential expression analysis should be performed separately for different alleles. Our ideas should greatly benefit the cancer community., (© 2023. The Author(s), under exclusive licence to Institute of Plant Genetics Polish Academy of Sciences.)

Background: At present, immunotherapy has become a powerful treatment for advanced gastric cancer (AGC), but not all patients can benefit from it. According to the latest research, the impact of B cell subpopulations on the immune microenvironment of gastric cancer (GC) is unknown. Exploring whether the interaction between B cells and tumor cells in GC affects the effectiveness of immunotherapy has attracted our interest., Methods: This study involved the re-analysis of single-cell RNA (scRNA) and spatial transcriptomics (ST) data from publicly available datasets. The focus was on investigating the subpopulations and differentiation trajectories of B cells in the gastric cancer (GC) tumor immune microenvironment (TIME). Spatial transcriptomics (ST) and multiple immunofluorescence (mIF) revealed a clear co-localization pattern between B cells and tumor cells. Multiple immunotherapy datasets were collected to identify unique immunotherapy biomarkers. The unique immunotherapeutic potential of targeting CCL28 was validated through a mouse gastric cancer model. In addition, flow cytometry revealed changes in the tumor immune microenvironment targeting CCL28., Results: The re-analysis of ST data from multiple cancer types revealed a co-localization pattern between B cells and tumor cells. A significant number of IgA plasma cells were identified in the GC TIME. Five different tumor-infiltrating B cell subpopulations and two unique B cell differentiation trajectories were characterized, along with seven GC-related states. By analyzing the communication between GC cells and B cells, it was further discovered that tumor cells can influence and recruit plasma cells through CCL28-CCR10 signaling. Additionally, there was a crosstalk between GC cells and B cells. Finally, we identified the LAMA/CD44 signaling axis as a potential prognostic marker for immunotherapy through a large amount of immunotherapy data. We also validated through various animal tumor models that targeting CCL28 can significantly promote CD8 + T cell infiltration and function in the TME by regulating B cell and plasma cell functions, and has the ability to synergize immunotherapy., Conclusion: The co-localization and crosstalk between GC cells and B cells significantly affect the efficacy of immunotherapy, and inhibiting the CCL28-CCR10 signal axis is a potential immunotherapy target for GC. Meanwhile, LAMA/CD44 pair may be a potential adverse indicator for immunotherapy and tumor prognosis., (© 2024. The Author(s).)

Background: Localized scleroderma (LoS) is an autoimmune disorder that primarily affects the skin, and is often treated with autologous fat grafting (AFG). Nevertheless, the retention rate of AFG in patients with LoS is typically low. We hypothesize that the low retention rate may be partially attributed to the inherent abnormalities of adipose-derived stem cells (ASCs) from nonlesional sites of patients with LoS., Methods: We performed a comparative analysis of the single-cell transcriptome of the SVF from nonlesional sites of patients with LoS and healthy donors, including cellular compositional analysis, differential expression analysis, and high-dimensional weighted gene coexpression network analysis. Experimental validation with fluorescence-activated cell sorting and bleomycin-induced skin fibrosis mice models were conducted., Results: We found a significant reduction in the relative proportion of CD55 high interstitial progenitors in ASCs under the condition of LoS. Differential expression analysis revealed inherent abnormalities of ASCs from patients with LoS, including enhanced fibrogenesis, reduced anti-inflammatory properties, and increased oxidative stress. Compared with CD55 low ASCs, CD55 high ASCs expressed significantly higher levels of secreted protein genes that had functions related to anti-inflammation and tissue regeneration (such as CD55, MFAP5, and METRNL). Meanwhile, CD55 high ASCs expressed significantly lower levels of secreted protein genes that promote inflammation, such as chemokine and complement protein genes. Furthermore, we provided in vivo experimental evidence that CD55 high ASCs had superior treatment efficacy compared with CD55 low ASCs in bleomycin-induced skin fibrosis mice models., Conclusions: We found that the low retention rate of AFG may be partially ascribed to the reduced pool of interstitial progenitor cells (CD55 high ) present within the ASC population in patients with LoS. We demonstrated the potential for improving the efficacy of AFG in the treatment of LoS by restoring the pool of interstitial progenitors within ASCs. Our study has significant implications for the field of translational regenerative medicine., (© 2024. The Author(s).)

Background: Whole blood host transcript signatures show great potential for diagnosis of infectious and inflammatory illness, with most published signatures performing binary classification tasks. Barriers to clinical implementation include validation studies, and development of strategies that enable simultaneous, multiclass diagnosis of febrile illness based on gene expression., Methods: We validated five distinct diagnostic signatures for paediatric infectious diseases in parallel using a single NanoString nCounter® experiment. We included a novel 3-transcript signature for childhood tuberculosis, and four published signatures which differentiate bacterial infection, viral infection, or Kawasaki disease from other febrile illnesses. Signature performance was assessed using receiver operating characteristic curve statistics. We also explored conceptual frameworks for multiclass diagnostic signatures, including additional transcripts found to be significantly differentially expressed in previous studies. Relaxed, regularised logistic regression models were used to derive two novel multiclass signatures: a mixed One-vs-All model (MOVA), running multiple binomial models in parallel, and a full-multiclass model. In-sample performance of these models was compared using radar-plots and confusion matrix statistics., Results: Samples from 91 children were included in the study: 23 bacterial infections (DB), 20 viral infections (DV), 14 Kawasaki disease (KD), 18 tuberculosis disease (TB), and 16 healthy controls. The five signatures tested demonstrated cross-platform performance similar to their primary discovery-validation cohorts. The signatures could differentiate: KD from other diseases with area under ROC curve (AUC) of 0.897 [95% confidence interval: 0.822-0.972]; DB from DV with AUC of 0.825 [0.691-0.959] (signature-1) and 0.867 [0.753-0.982] (signature-2); TB from other diseases with AUC of 0.882 [0.787-0.977] (novel signature); TB from healthy children with AUC of 0.910 [0.808-1.000]. Application of signatures outside of their designed context reduced performance. In-sample error rates for the multiclass models were 13.3% for the MOVA model and 0.0% for the full-multiclass model. The MOVA model misclassified DB cases most frequently (18.7%) and TB cases least (2.7%)., Conclusions: Our study demonstrates the feasibility of NanoString technology for cross-platform validation of multiple transcriptomic signatures in parallel. This external cohort validated performance of all five signatures, including a novel sparse TB signature. Two exploratory multi-class models showed high potential accuracy across four distinct diagnostic groups., (© 2024. The Author(s).)

Down syndrome (DS), the genetic condition caused by trisomy 21 (T21), is characterized by delayed neurodevelopment, accelerated aging, and increased risk of many co-occurring conditions. Hypoxemia and dysregulated hematopoiesis have been documented in DS, but the underlying mechanisms and clinical consequences remain ill defined. We report an integrative multi-omic analysis of ∼400 research participants showing that people with DS display transcriptomic signatures indicative of elevated heme metabolism and increased hypoxic signaling across the lifespan, along with chronic overproduction of erythropoietin, elevated biomarkers of tissue-specific hypoxia, and hallmarks of stress erythropoiesis. Elevated heme metabolism, transcriptional signatures of hypoxia, and stress erythropoiesis are conserved in a mouse model of DS and associated with overexpression of select triplicated genes. These alterations are independent of the hyperactive interferon signaling characteristic of DS. These results reveal lifelong dysregulation of key oxygen-related processes that could contribute to the developmental and clinical hallmarks of DS., Competing Interests: Declaration of interests J.M.E. has provided consulting services to Eli Lilly and Co., Gilead Sciences Inc., Biohaven Pharmaceuticals, and Perha Pharmaceuticals., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Although RNA structures play important roles in regulating gene expression, the mechanism and function of mRNA folding in plant bacterial pathogens remain elusive. Therefore, we perform dimethyl sulfate sequencing (DMS-seq) on the Pseudomonas syringae under nutrition-rich and -deficient conditions, revealing that the mRNA structure changes substantially in the minimal medium (MM) that tunes global translation efficiency (TE), thereby inducing virulence. This process is led by the increased expression of hfq, which is directly activated by transcription regulators RpoS and CysB. The co-occurrence of Hfq and RpoS in diverse bacteria and the deep conservation of Hfq Y25 is critical for RNA-mediated regulation and implicates the wider biological importance of mRNA structure and feedback loops in the control of global gene expression., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Bacterial gene expression is a complex process involving extensive regulatory mechanisms. Along with growing interests in this field, Nanopore Direct RNA Sequencing (DRS) provides a promising platform for rapid and comprehensive characterization of bacterial RNA biology. However, the DRS of bacterial RNA is currently deficient in the yield of mRNA-mapping reads and has yet to be exploited for transcriptome-wide RNA modification mapping. Here, we showed that pre-processing of bacterial total RNA (size selection followed by ribosomal RNA depletion and polyadenylation) guaranteed high throughputs of sequencing data and considerably increased the amount of mRNA reads. This way, complex transcriptome architectures were reconstructed for Escherichia coli and Staphylococcus aureus and extended the boundaries of 225 known E. coli operons and 89 defined S. aureus operons. Utilizing unmodified in vitro-transcribed (IVT) RNA libraries as a negative control, several Nanopore-based computational tools globally detected putative modification sites in the E. coli and S. aureus transcriptomes. Combined with Next-Generation Sequencing-based N6-methyladenosine (m6A) detection methods, 75 high-confidence m6A candidates were identified in the E. coli protein-coding transcripts, while none were detected in S. aureus. Altogether, we demonstrated the potential of Nanopore DRS in systematic and convenient transcriptome and epitranscriptome analysis., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Accurate detection and quantification of mRNA isoforms from nanopore long-read sequencing remains challenged by technical noise, particularly in single cells. To address this, we introduce Isosceles, a computational toolkit that outperforms other methods in isoform detection sensitivity and quantification accuracy across single-cell, pseudo-bulk and bulk resolution levels, as demonstrated using synthetic and biologically-derived datasets. Here we show Isosceles improves the fidelity of single-cell transcriptome quantification at the isoform-level, and enables flexible downstream analysis. As a case study, we apply Isosceles, uncovering coordinated splicing within and between neuronal differentiation lineages. Isosceles is suitable to be applied in diverse biological systems, facilitating studies of cellular heterogeneity across biomedical research applications., (© 2024. The Author(s).)

Bronchiolitis is the leading cause of infant hospitalization. However, the molecular networks driving bronchiolitis pathobiology remain unknown. Integrative molecular networks, including the transcriptome and metabolome, can identify functional and regulatory pathways contributing to disease severity. Here, we integrated nasopharyngeal transcriptome and metabolome data of 397 infants hospitalized with bronchiolitis in a 17-center prospective cohort study. Using an explainable deep network model, we identified an omics-cluster comprising 401 transcripts and 38 metabolites that distinguishes bronchiolitis severity (test-set AUC, 0.828). This omics-cluster derived a molecular network, where innate immunity-related metabolites (e.g., ceramides) centralized and were characterized by toll-like receptor (TLR) and NF-κB signaling pathways (both FDR < 0.001). The network analyses identified eight modules and 50 existing drug candidates for repurposing, including prostaglandin I 2 analogs (e.g., iloprost), which promote anti-inflammatory effects through TLR signaling. Our approach facilitates not only the identification of molecular networks underlying infant bronchiolitis but the development of pioneering treatment strategies., (© 2024. The Author(s).)

Retinal progenitor cells (RPCs) are a multipotent and highly proliferative population that give rise to all retinal cell types during organogenesis. Defining their molecular signature is a key step towards identifying suitable approaches to treat visual impairments. Here, we performed RNA sequencing of whole eyes from Xenopus at three embryonic stages and used differential expression analysis to define the transcriptomic profiles of optic tissues containing proliferating and differentiating RPCs during retinogenesis. Gene Ontology and KEGG pathway analyses showed that genes associated with developmental pathways (including Wnt and Hedgehog signaling) were upregulated during the period of active RPC proliferation in early retinal development (Nieuwkoop Faber st. 24 and 27). Developing eyes had dynamic expression profiles and shifted to enrichment for metabolic processes and phototransduction during RPC progeny specification and differentiation (st. 35). Furthermore, conserved adult eye regeneration genes were also expressed during early retinal development, including sox2 , pax6 , nrl , and Notch signaling components. The eye transcriptomic profiles presented here span RPC proliferation to retinogenesis and include regrowth-competent stages. Thus, our dataset provides a rich resource to uncover molecular regulators of RPC activity and will allow future studies to address regulators of RPC proliferation during eye repair and regrowth.

One snapshot of the peer review process for "Transcriptome data are insufficient to control false discoveries in regulatory network inference" (Kernfeld et al., 2024). 1 ., Competing Interests: Declaration of interests The author declares no competing interests., (Copyright © 2024. Published by Elsevier Inc.)

Inference of causal transcriptional regulatory networks (TRNs) from transcriptomic data suffers notoriously from false positives. Approaches to control the false discovery rate (FDR), for example, via permutation, bootstrapping, or multivariate Gaussian distributions, suffer from several complications: difficulty in distinguishing direct from indirect regulation, nonlinear effects, and causal structure inference requiring "causal sufficiency," meaning experiments that are free of any unmeasured, confounding variables. Here, we use a recently developed statistical framework, model-X knockoffs, to control the FDR while accounting for indirect effects, nonlinear dose-response, and user-provided covariates. We adjust the procedure to estimate the FDR correctly even when measured against incomplete gold standards. However, benchmarking against chromatin immunoprecipitation (ChIP) and other gold standards reveals higher observed than reported FDR. This indicates that unmeasured confounding is a major driver of FDR in TRN inference. A record of this paper's transparent peer review process is included in the supplemental information., Competing Interests: Declaration of interests A.B. is a stockholder for Alphabet, Inc.; has consulted for Third Rock Ventures; and is a founder of CellCipher, Inc., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: The availability of transcriptomic data for species without a reference genome enables the construction of de novo transcriptome assemblies as alternative reference resources from RNA-Seq data. A transcriptome provides direct information about a species' protein-coding genes under specific experimental conditions. The de novo assembly process produces a unigenes file in FASTA format, subsequently targeted for the annotation. Homology-based annotation, a method to infer the function of sequences by estimating similarity with other sequences in a reference database, is a computationally demanding procedure., Results: To mitigate the computational burden, we introduce HPC-T-Annotator, a tool for de novo transcriptome homology annotation on high performance computing (HPC) infrastructures, designed for straightforward configuration via a Web interface. Once the configuration data are given, the entire parallel computing software for annotation is automatically generated and can be launched on a supercomputer using a simple command line. The output data can then be easily viewed using post-processing utilities in the form of Python notebooks integrated in the proposed software., Conclusions: HPC-T-Annotator expedites homology-based annotation in de novo transcriptome assemblies. Its efficient parallelization strategy on HPC infrastructures significantly reduces computational load and execution times, enabling large-scale transcriptome analysis and comparison projects, while its intuitive graphical interface extends accessibility to users without IT skills., (© 2024. The Author(s).)

Single-cell transcriptomics reveals significant variations in transcriptional activity across cells. Yet, it remains challenging to identify mechanisms of transcription dynamics from static snapshots. It is thus still unknown what drives global transcription dynamics in single cells. We present a stochastic model of gene expression with cell size- and cell cycle-dependent rates in growing and dividing cells that harnesses temporal dimensions of single-cell RNA sequencing through metabolic labeling protocols and cel lcycle reporters. We develop a parallel and highly scalable approximate Bayesian computation method that corrects for technical variation and accurately quantifies absolute burst frequency, burst size, and degradation rate along the cell cycle at a transcriptome-wide scale. Using Bayesian model selection, we reveal scaling between transcription rates and cell size and unveil waves of gene regulation across the cell cycle-dependent transcriptome. Our study shows that stochastic modeling of dynamical correlations identifies global mechanisms of transcription regulation. A record of this paper's transparent peer review process is included in the supplemental information., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

This study introduces a new imaging, spatial transcriptomics (ST), and single-cell RNA-sequencing integration pipeline to characterize neoplastic cell state transitions during tumorigenesis. We applied a semi-supervised analysis pipeline to examine premalignant pancreatic intraepithelial neoplasias (PanINs) that can develop into pancreatic ductal adenocarcinoma (PDAC). Their strict diagnosis on formalin-fixed and paraffin-embedded (FFPE) samples limited the single-cell characterization of human PanINs within their microenvironment. We leverage whole transcriptome FFPE ST to enable the study of a rare cohort of matched low-grade (LG) and high-grade (HG) PanIN lesions to track progression and map cellular phenotypes relative to single-cell PDAC datasets. We demonstrate that cancer-associated fibroblasts (CAFs), including antigen-presenting CAFs, are located close to PanINs. We further observed a transition from CAF-related inflammatory signaling to cellular proliferation during PanIN progression. We validate these findings with single-cell high-dimensional imaging proteomics and transcriptomics technologies. Altogether, our semi-supervised learning framework for spatial multi-omics has broad applicability across cancer types to decipher the spatiotemporal dynamics of carcinogenesis., Competing Interests: Declaration of interests E.M.J. reports other support from Abmeta; personal fees from Genocea; personal fees from Achilles; personal fees from DragonFly; personal fees from Candel Therapeutics; other support from the Parker Institute; grants and other support from Lustgarten; personal fees from Carta; grants and other support from Genentech; grants and other support from AstraZeneca; personal fees from NextCure; and grants and other support from Break Through Cancer outside of the submitted work. E.J.F. is on the Scientific Advisory Board of Viosera Therapeutics/Resistance Bio and is a consultant to Mestag Therapeutics. W.J.H. reports patent royalties from Rodeo/Amgen and speaking/travel honoraria from Exelixis and Standard BioTools., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Glaesserella parasuis (GPS) can cause severe systemic inflammation in pigs, resulting in huge economic losses to the pig industry. At present, no effective method is available for the prevention and control of GPS infection. Molecular breeding for disease resistance is imminent, but disease-resistance genes have not been identified. To study the mechanism of systemic acute inflammation caused by GPS, we established three in vitro infection models (3D4/21 cells, PK15 cells, and PAVEC cells) according to its infection path. There was no significant difference in apoptosis among the three kinds of cells after 12 h of continuous GPS stimulation, while inflammatory factors were significantly upregulated. Subsequent transcriptome analysis revealed 1969, 1207, and 3564 differentially expressed genes (DEGs) in 3D4/21 cells, PK15 cells, and PAVEC cells, respectively, after GPS infection. Many of the DEGs were predicted to be associated with inflammatory responses (C3, CD44, etc.); cell proliferation, growth and apoptosis; gene expression; and protein phosphorylation. Key signaling pathways, including S100 family signaling, bacteria and virus recognition, and pathogen-induced cytokine storm signaling, were enriched based on Ingenuity Pathway Analysis (IPA). Furthermore, a total of three putative transmembrane receptors and two putative G-protein-coupled receptors, namely F3, ICAM1, PLAUR, ACKR3, and GPRC5A, were identified by IPA among the three types of cells. ACKR3 and GPRC5A play pivotal roles in bacterial adhesion, invasion, host immune response and inflammatory response through the S100 family signaling pathway. Our findings provide new insights into the pathological mechanisms underlying systemic inflammation caused by GPS infection in pigs, and they lay a foundation for further research on disease-resistance breeding to GPS.

The cortical microenvironment surrounding malignant glioblastoma is a source of depolarizing crosstalk favoring hyperexcitability, tumor expansion, and immune evasion. Neosynaptogenesis, excess glutamate, and altered intrinsic membrane currents contribute to excitability dyshomeostasis, yet only half of the cases develop seizures, suggesting that tumor and host genomics, along with location, rather than mass effect, play a critical role. We analyzed the spatial contours and expression of 358 clinically validated human epilepsy genes in the human glioblastoma transcriptome compared to non-tumor adult and developing cortex datasets. Nearly half, including dosage-sensitive genes whose expression levels are securely linked to monogenic epilepsy, are strikingly enriched and aberrantly regulated at the leading edge, supporting a complex epistatic basis for peritumoral epileptogenesis. Surround hyperexcitability induced by complex patterns of proepileptic gene expression may explain the limited efficacy of narrowly targeted antiseizure medicines and the persistence of epilepsy following tumor resection and clarify why not all brain tumors provoke seizures., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Background: A correct and stably expressing reference gene is prerequisite for successful quantitative real-time PCR (qRT-PCR). Investigating gene expression profiling during flower development could enhance our understanding of the molecular mechanisms of flower formation and fertility in Lycium., Methods and Results: In this study, 11 candidate reference genes in Lycium flower development were selected from transcriptome sequence data and evaluated with five traditional housekeeping genes from previous studies based on qRT-PCR amplification. Comparing the expression stability result of 16 candidate genes using GeNorm, NormFinder, BestKeeper, and Delta Ct algorithms, Lba04g01649 and Lba12g02820 were validated as the optimal reference genes for the flower development of Lycium., Conclusions: The reference genes identified in this study would improve the accuracy of qRT-PCR quantification of target gene expression in Lycium flower development and facilitate future functional genomics studies on flower development. This research could lay the foundation for the study of the reproduction and development of the Lycium flower., (© 2024. The Author(s), under exclusive licence to Springer Nature B.V.)

Gestational diabetes mellitus (GDM) presents varied manifestations throughout pregnancy and poses a complex clinical challenge. High-depth cell-free DNA (cfDNA) sequencing analysis holds promise in advancing our understanding of GDM pathogenesis and prediction. In 299 women with GDM and 299 matched healthy pregnant women, distinct cfDNA fragment characteristics associated with GDM are identified throughout pregnancy. Integrating cfDNA profiles with lipidomic and single-cell transcriptomic data elucidates functional changes linked to altered lipid metabolism processes in GDM. Transcription start site (TSS) scores in 50 feature genes are used as the cfDNA signature to distinguish GDM cases from controls effectively. Notably, differential coverage of the islet acinar marker gene PRSS1 emerges as a valuable biomarker for GDM. A specialized neural network model is developed, predicting GDM occurrence and validated across two independent cohorts. This research underscores the high-depth cfDNA early prediction and characterization of GDM, offering insights into its molecular underpinnings and potential clinical applications., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Identifying patients with stage II and III colon cancer who will benefit from 5-fluorouracil (5-FU)-based adjuvant chemotherapy is crucial for the advancement of personalized cancer therapy. We employ a semi-supervised machine learning approach to analyze a large dataset with 933 stage II and III colon cancer samples. Our analysis leverages gene regulatory networks to discover an 18-gene prognostic signature and to explore a 10-gene signature that potentially predicts chemotherapy benefits. The 10-gene signature demonstrates strong prognostic power and shows promising potential to predict chemotherapy benefits. We establish a robust clinical assay on the NanoString nCounter platform, validated in a retrospective formalin-fixed paraffin-embedded (FFPE) cohort, which represents an important step toward clinical application. Our study lays the groundwork for improving adjuvant chemotherapy and potentially expanding into immunotherapy decision-making in colon cancer. Future prospective studies are needed to validate and establish the clinical utility of the 10-gene signature in clinical settings., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Proneural transcription factors establish molecular cascades to orchestrate neuronal diversity. One such transcription factor, Atonal homolog 1 (Atoh1), gives rise to cerebellar excitatory neurons and over 30 distinct nuclei in the brainstem critical for hearing, breathing, and balance. Although Atoh1 lineage neurons have been qualitatively described, the transcriptional programs that drive their fate decisions and the full extent of their diversity remain unknown. Here, we analyzed single-cell RNA sequencing and ATOH1 DNA binding in Atoh1 lineage neurons of the developing mouse hindbrain. This high-resolution dataset identified markers for specific brainstem nuclei and demonstrated that transcriptionally heterogeneous progenitors require ATOH1 for proper migration. Moreover, we identified a sizable population of proliferating unipolar brush cell progenitors in the mouse Atoh1 lineage, previously described in humans as the origin of one medulloblastoma subtype. Collectively, our data provide insights into the developing mouse hindbrain and markers for functional assessment of understudied neuronal populations., Competing Interests: Declaration of interests R.S.D. is a paid consultant for AstraZeneca. H.Y.Z. is a member of the Board of Regeneron, co-founder of Cajal Neuroscience, and on the science advisory board of the Column group. This is work in the Zoghbi Lab and has not been licensed anywhere. The work in this paper has no connection to any of the consulting or companies named here., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Huntington's disease (HD) is a neurodegenerative genetic disorder caused by an expansion in the CAG repeat tract of the huntingtin (HTT) gene resulting in behavioural, cognitive, and motor defects. Current knowledge of disease pathogenesis remains incomplete, and no disease course-modifying interventions are in clinical use. We have previously reported the development and characterisation of the OVT73 transgenic sheep model of HD. The 73 polyglutamine repeat is somatically stable and therefore likely captures a prodromal phase of the disease with an absence of motor symptomatology even at 5-years of age and no detectable striatal cell loss. To better understand the disease-initiating events we have undertaken a single nuclei transcriptome study of the striatum of an extensively studied cohort of 5-year-old OVT73 HD sheep and age matched wild-type controls. We have identified transcriptional upregulation of genes encoding N-methyl-D-aspartate (NMDA), α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) and kainate receptors in medium spiny neurons, the cell type preferentially lost early in HD. Further, we observed an upregulation of astrocytic glutamate uptake transporters and medium spiny neuron GABAA receptors, which may maintain glutamate homeostasis. Taken together, these observations support the glutamate excitotoxicity hypothesis as an early neurodegeneration cascade-initiating process but the threshold of toxicity may be regulated by several protective mechanisms. Addressing this biochemical defect early may prevent neuronal loss and avoid the more complex secondary consequences precipitated by cell death., (© The Author(s) 2024. Published by Oxford University Press.)

Glioblastoma (GBM) is a highly prevalent and deadly brain tumor with high mortality rates, especially among adults. Despite extensive research, the underlying mechanisms driving its progression remain poorly understood. Computational analysis offers a powerful approach to explore potential prognostic biomarkers, drug targets, and therapeutic agents for GBM. In this study, we utilized three gene expression datasets from the Gene Expression Omnibus (GEO) database to identify differentially expressed genes (DEGs) associated with GBM progression. Our goal was to uncover key molecular players implicated in GBM pathogenesis and potential avenues for targeted therapy. Analysis of the gene expression datasets revealed a total of 78 common DEGs that are potentially involved in GBM progression. Through further investigation, we identified nine hub DEGs that are highly interconnected in protein-protein interaction (PPI) networks, indicating their central role in GBM biology. Gene Ontology (GO) and pathway enrichment analyses provided insights into the biological processes and immunological pathways influenced by these DEGs. Among the nine identified DEGs, survival analysis demonstrated that increased expression of GMFG correlated with decreased patient survival rates in GBM, suggesting its potential as a prognostic biomarker and preventive target for GBM. Furthermore, molecular docking and ADMET analysis identified two compounds from the NIH clinical collection that showed promising interactions with the GMFG protein. Besides, a 100 nanosecond molecular dynamics (MD) simulation evaluated the conformational changes and the binding strength. Our study highlights the potential of GMFG as both a prognostic biomarker and a therapeutic target for GBM. The identification of GMFG and its associated pathways provides valuable insights into the molecular mechanisms driving GBM progression. Moreover, the identification of candidate compounds with potential interactions with GMFG offers exciting possibilities for targeted therapy development. However, further laboratory experiments are required to validate the role of GMFG in GBM pathogenesis and to assess the efficacy of potential therapeutic agents targeting this molecule., Competing Interests: The authors declare no conflicts of interest., (Copyright © 2024 Md. Wasim Alom et al.)

The homologous E6-AP carboxy-terminal structural domain (HECT) contained in E3 ubiquitin ligases (E3s) is a key factor in protein degradation and maintenance of cellular homeostasis in animals. However, the functional roles and evolutionary aspects of the HECT gene family in bivalve mussels remain unclear and warrant further investigation. In this study, we identified 22 HECT genes within the genome of Mytilus coruscus Gould, all containing a conserved HECT structural domain derived from dispersed repeats, distributed unevenly across 11 chromosomes. Phylogenetic analysis classified M. coruscus HECT genes into six major classes, with amino acid sequences within the same evolutionary clade displaying similar conserved motifs. Homology analysis with HECT genes of four bivalve species revealed that M. coruscus and Mytilus galloprovincialis possessed the largest number of homologous gene pairs, showing a significant correlation between the two in the evolution of the HECT gene family. Homology analysis with HECT genes of four bivalve species revealed that M. coruscus and M. galloprovincialis possessed the largest number of homologous gene pairs, showing a significant correlation between the two in the evolution of the HECT gene family. M. coruscus exhibited pronounced and specific expression in gills and blood tissues. Notably, Mco&#95;UPL3 gene expression was significantly upregulated after 12 h of acute heat stress (33 °C) and 24 h of Vibrio injection (0.4 OD). Gene ontology analysis of the HECT genes in M. coruscus revealed that it is primarily enriched in protein modification and degradation functions. This suggests that HECT genes may play a key role in protein degradation and immunomodulation in M. coruscus . These findings offer valuable insights for the breeding of stress-tolerant traits in M. coruscus . In summary, our data shed light on the potential functions of HECT E3 ligases in response to heat stress and Vibrio infection, providing practical guidance for enhancing resilience through breeding in M. coruscus .