Noise-robust soft clustering of gene expression time-course data

Proteomic characterization of murine hematopoietic stem progenitor cells reveals dynamic fetal-to-adult changes in metabolic-related pathways

Hematopoietic stem progenitor cells (HSPCs) give rise to the hematopoietic system, maintain hematopoiesis throughout the lifespan, and undergo molecular and functional changes during their development and aging. The importance of hematopoietic stem cell (HSC) biology has led to their extensive characterization at genomic and transcriptomic levels. However, the proteomics of HSPCs throughout the murine lifetime still needs to be fully completed. Here, using mass spectrometry (MS)-based quantitative proteomics, we report on the dynamic changes in the proteome of HSPCs from four developmental stages in the fetal liver (FL) and the bone marrow (BM), including E14.5, young (2 months), middle-aged (8 months), and aging (18 months) stages. Proteomics unveils highly dynamic protein kinetics during the development and aging of HSPCs. Our data identify stage-specific developmental features of HSPCs, which can be linked to their functional maturation and senescence. Our proteomic data demonstrated that FL HSPCs depend on aerobic respiration to meet their proliferation and oxygen supply demand, while adult HSPCs prefer glycolysis to preserve the HSC pool. By functional assays, we validated the decreased mitochondrial metabolism, glucose uptake, reactive oxygen species (ROS) production, protein synthesis rate, and increased glutathione S-transferase (GST) activity during HSPC development from fetal to adult. Distinct metabolism pathways and immune-related pathways enriched in different HSPC developmental stages were revealed at the protein level. Our study will have broader implications for understanding the mechanism of stem cell maintenance and fate determination and reversing the HSC aging process.

Quercetin intervention mitigates small intestinal damage and immunologic derangement induced by polystyrene nanoplastics: Insights from multi-omics analysis in mice

Nanoplastics (NPs), which belong to emerging environmental pollutants, threaten environmental sustainability and human health. Despite recent studies have reported that NPs damage the gastrointestinal tract and immune homeostasis, the underlying mechanisms remain unclear. Polyphenols have been found to promote NPs excretion by interacting with intestinal flora (IF). However, the potential mechanisms and action targets of this are still poorly understood. To address these knowledge gaps, we investigated the impact of quercetin and three concentrations of polystyrene nanoplastics (PS-NPs) in mice using an integrated phenotypic and multi-omics analysis. Our findings demonstrated that PS-NPs accumulate within the intestine, resulting in impairments to intestinal tissue and barrier function, as well as disturbing the expression of immune-response small intestinal genes and composition of IF. Exposure to PS-NPs significantly elevate the level of intestinal IgG and CD20+ B cells, while inhibiting T cells activation. Furthermore, PS-NPs could induce systemic immune and serum insulin level disorders. Quercetin might mitigate PS-NPs-induced intestinal damage and immune disorders though reversing IF disorders, gene expression changes, and their interaction.

Mannose-6-Phosphate Isomerase Functional Status Shapes a Rearrangement in the Proteome and Degradome of Mannose-Treated Melanoma Cells

Metabolic reprogramming is a ubiquitous feature of transformed cells, comprising one of the hallmarks of cancer and enabling neoplastic cells to adapt to new environments. Accumulated evidence reports on the failure of some neoplastic cells to convert mannose-6-phosphate into fructose-6-phosphate, thereby impairing tumor growth in cells displaying low levels of mannose-6-phosphate isomerase (MPI). Thus, we performed functional analyses and profiled the proteome landscape and the repertoire of substrates of proteases (degradome) of melanoma cell lines with distinct mutational backgrounds submitted to treatment with mannose. Our results suggest a significant rearrangement in the proteome and degradome of melanoma cell lines upon mannose treatment including the activation of catabolic pathways (such as protein turnover) and differences in protein N-terminal acetylation. Even though MPI protein abundance and gene expression status are not prognostic markers, perturbation in the network caused by an exogenous monosaccharide source (i.e., mannose) significantly affected the downstream interconnected biological circuitry. Therefore, as reported in this study, the proteomic/degradomic mapping of mannose downstream effects due to the metabolic rewiring caused by the functional status of the MPI enzyme could lead to the identification of specific molecular players from affected signaling circuits in melanoma.

Metformin decelerates aging clock in male monkeys

In a rigorous 40-month study, we evaluated the geroprotective effects of metformin on adult male cynomolgus monkeys, addressing a gap in primate aging research. The study encompassed a comprehensive suite of physiological, imaging, histological, and molecular evaluations, substantiating metformin's influence on delaying age-related phenotypes at the organismal level. Specifically, we leveraged pan-tissue transcriptomics, DNA methylomics, plasma proteomics, and metabolomics to develop innovative monkey aging clocks and applied these to gauge metformin's effects on aging. The results highlighted a significant slowing of aging indicators, notably a roughly 6-year regression in brain aging. Metformin exerts a substantial neuroprotective effect, preserving brain structure and enhancing cognitive ability. The geroprotective effects on primate neurons were partially mediated by the activation of Nrf2, a transcription factor with anti-oxidative capabilities. Our research pioneers the systemic reduction of multi-dimensional biological age in primates through metformin, paving the way for advancing pharmaceutical strategies against human aging.

Contrasting the development of larval and adult body plans during the evolution of biphasic lifecycles in sea urchins

Biphasic lifecycles are widespread among animals, but little is known about how the developmental transition between larvae and adults is regulated. Sea urchins are a unique system for studying this phenomenon because of the stark differences between their bilateral larval and pentaradial adult body plans. Here, we use single-cell RNA sequencing to analyze the development of Heliocidaris erythrogramma (He), a sea urchin species with an accelerated, non-feeding mode of larval development. The sequencing time course extends from embryogenesis to roughly a day before the onset of metamorphosis in He larvae, which is a period that has not been covered by previous datasets. We find that the non-feeding developmental strategy of He is associated with several changes in the specification of larval cell types compared to sea urchins with feeding larvae, such as the loss of a larva-specific skeletal cell population. Furthermore, the development of the larval and adult body plans in sea urchins may utilize largely different sets of regulatory genes. These findings lay the groundwork for extending existing developmental gene regulatory networks to cover additional stages of biphasic lifecycles.

Graphene quantum dots disrupt the mitochondrial potential of Trypanosoma brucei by interacting with the p18 subunit of ATP synthase F1 after endocytosis via the VSG recycling pathway

HYPOTHESIS

Trypanosomiasis is one of the main threats to human and animal health in African countries. Trypanosoma brucei can evade the host immune recognition by rapidly altering its variant surface glycoprotein (VSG). The ATP synthase F1 subunit of the parasite exhibits extremely low similarity to that of its mammalian hosts, hypothetically making it an ideal target for the development of novel therapeutics.

EXPERIMENTS

Graphene quantum dots (GQDs) were synthesized, and their adhesion to T. brucei surface and internalization was observed microscopically. The activity of ATP synthase and mitochondrial membrane potential of T. brucei were measured after exposure to GQDs. Proteomics, biolayer interferometry, and molecular dynamic simulations were utilized to evaluate the interaction between GQDs with the target proteins.

FINDINGS

GQDs specifically adhered to the VSG of T. brucei and were conveyed inside the parasite via the VSG internalization pathway. The GQDs promoted intracellular ROS production, interacted with, and inhibited the activity of the p18 subunit of ATP synthase, disrupted parasite mitochondrial membrane potential. Additionally, the GQDs caused a decrease in aminoacyl - tRNA biosynthesis, and upregulated RNA and protein degradation pathways. The findings of this study offer a novel avenue for the target-oriented discovery of anti-trypanosome drugs.

Exploring the antibiotic resistance genes removal dynamics in chicken manure by composting

Prolonged antibiotic usage in livestock farming leads to the accumulation of antibiotic resistance genes in animal manure. Composting has been shown as an effective way of removing antibiotic resistance from manures, but the specific mechanisms remain unclear. This study used time-series sampling and metagenomics to analyse the resistome types and their bacterial hosts in chicken manures. Composting significantly altered the physicochemical properties and microbiome composition, reduced antibiotic resistance genes by 65.71 %, mobile genetic elements by 68.15 % and horizontal gene transfer frequency. Source tracking revealed that Firmicutes, Actinobacteria, and Proteobacteria are the major bacterial hosts involved in the resistome and gene transfer events. Composting reduces the resistome risk by targeting pathogens such as Staphylococcus aureus. Structural equation modelling confirmed that composting reduces resistome risk by changing pH and pathogen abundance. This study demonstrates that composting is an effective strategy for mitigating resistome risk in chicken manure, thereby supporting the One Health initiative.

Local protein synthesis at neuromuscular synapses is required for motor functions

Motor neurons are highly polarized, and their axons extend over great distances to form connections with myofibers via neuromuscular junctions (NMJs). Local translation at the NMJs in vivo has not been identified. Here, we utilized motor neuron-labeled RiboTag mice and the TRAP (translating ribosome affinity purification) technique to spatiotemporally profile the translatome at NMJs. We found that mRNAs associated with glucose catabolism, synaptic connection, and protein homeostasis are enriched at presynapses. Local translation at the synapse shifts from the assembly of cytoskeletal components during early developmental stages to energy production in adulthood. The mRNA of neuronal Agrin (Agrn), the key molecule for NMJ assembly, is present at motor axon terminals and locally translated. Disrupting the axonal location of Agrn mRNA causes impairment of synaptic transmission and motor functions in adult mice. Our findings indicate that spatiotemporal regulation of mRNA local translation at NMJs plays critical roles in synaptic transmission and motor functions in vivo.

A polarized FGF8 source specifies frontotemporal signatures in spatially oriented cell populations of cortical assembloids

Organoids generating major cortical cell types in distinct compartments are used to study cortical development, evolution and disorders. However, the lack of morphogen gradients imparting cortical positional information and topography in current systems hinders the investigation of complex phenotypes. Here, we engineer human cortical assembloids by fusing an organizer-like structure expressing fibroblast growth factor 8 (FGF8) with an elongated organoid to enable the controlled modulation of FGF8 signaling along the longitudinal organoid axis. These polarized cortical assembloids mount a position-dependent transcriptional program that in part matches the in vivo rostrocaudal gene expression patterns and that is lost upon mutation in the FGFR3 gene associated with temporal lobe malformations and intellectual disability. By producing spatially oriented cell populations with signatures related to frontal and temporal area identity within individual assembloids, this model recapitulates in part the early transcriptional divergence embedded in the protomap and enables the study of cortical area-relevant alterations underlying human disorders.

Proteomics reveals dynamic metabolic changes in human hematopoietic stem progenitor cells from fetal to adulthood

BACKGROUND

Hematopoietic stem progenitor cells (HSPCs) undergo phenotypical and functional changes during their emergence and development. Although the molecular programs governing the development of human hematopoietic stem cells (HSCs) have been investigated broadly, the relationships between dynamic metabolic alterations and their functions remain poorly characterized.

METHODS

In this study, we comprehensively described the proteomics of HSPCs in the human fetal liver (FL), umbilical cord blood (UCB), and adult bone marrow (aBM). The metabolic state of human HSPCs was assessed via a Seahorse assay, RT‒PCR, and flow cytometry-based metabolic-related analysis. To investigate whether perturbing glutathione metabolism affects reactive oxygen species (ROS) production, the metabolic state, and the expansion of human HSPCs, HSPCs were treated with buthionine sulfoximine (BSO), an inhibitor of glutathione synthetase, and N-acetyl-L-cysteine (NAC).

RESULTS

We investigated the metabolomic landscape of human HSPCs from the fetal, perinatal, and adult developmental stages by in-depth quantitative proteomics and predicted a metabolic switch from the oxidative state to the glycolytic state during human HSPC development. Seahorse assays, mitochondrial activity, ROS level, glucose uptake, and protein synthesis rate analysis supported our findings. In addition, immune-related pathways and antigen presentation were upregulated in UCB or aBM HSPCs, indicating their functional maturation upon development. Glutathione-related metabolic perturbations resulted in distinct responses in human HSPCs and progenitors. Furthermore, the molecular and immunophenotypic differences between human HSPCs at different developmental stages were revealed at the protein level for the first time.

CONCLUSION

The metabolic landscape of human HSPCs at three developmental stages (FL, UCB, and aBM), combined with proteomics and functional validations, substantially extends our understanding of HSC metabolic regulation. These findings provide valuable resources for understanding human HSC function and development during fetal and adult life.

Transcriptome Analysis and Identification of Genes Associated with Cotton Seed Size

Cotton seeds, as the main by-product of cotton, are not only an important raw material for edible oil and feed but also a source of biofuel. The quality of cotton seeds directly affects cotton planting and is closely related to the yield and fiber quality. However, the molecular mechanism governing cotton seed size remains largely unexplored. This study investigates the regulatory mechanisms of cotton seed size by focusing on two cotton genotypes, N10 and N12, which exhibit notable phenotypic variations across multiple environments. Developing seeds were sampled at various stages (5, 20, 30, and 35 DPA) and subjected to RNA-seq. Temporal pattern clustering and WGCNA on differentially expressed genes identified 413 candidate genes, including these related to sugar metabolism that were significantly enriched in transcriptional regulation. A genetic transformation experiment indicated that the overexpression of the GhUXS5 gene encoding UDP-glucuronate decarboxylase 5 significantly increased seed size, suggesting an important role of GhUXS5 in regulating cotton seed size. This discovery provides crucial insights into the molecular mechanisms controlling cotton seed size, helping to unravel the complex regulatory network and offering new strategies and targets for cotton breeding to enhance the economic value of cotton seeds and overall cotton yield.

Epigenetic insights into GABAergic development in Dravet Syndrome iPSC and therapeutic implications

Dravet syndrome (DS) is a devastating early-onset refractory epilepsy syndrome caused by variants in the SCN1A gene. A disturbed GABAergic interneuron function is implicated in the progression to DS but the underlying developmental and pathophysiological mechanisms remain elusive, in particularly at the chromatin level. Induced pluripotent stem cells (iPSCs) derived from DS cases and healthy donors were used to model disease-associated epigenetic abnormalities of GABAergic development. Chromatin accessibility was assessed at multiple time points (Day 0, Day 19, Day 35, and Day 65) of GABAergic differentiation. Additionally, the effects of the commonly used anti-seizure drug valproic acid (VPA) on chromatin accessibility were elucidated in GABAergic cells. The distinct dynamics in the chromatin profile of DS iPSC predicted accelerated early GABAergic development, evident at D19, and diverged further from the pattern in control iPSC with continued differentiation, indicating a disrupted GABAergic maturation. Exposure to VPA at D65 reshaped the chromatin landscape at a variable extent in different iPSC-lines and rescued the observed dysfunctional development of some DS iPSC-GABA. The comprehensive investigation on the chromatin landscape of GABAergic differentiation in DS-patient iPSC offers valuable insights into the epigenetic dysregulations associated with interneuronal dysfunction in DS. Moreover, the detailed analysis of the chromatin changes induced by VPA in iPSC-GABA holds the potential to improve the development of personalized and targeted anti-epileptic therapies.

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

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

Dissecting key regulators of transcriptome kinetics through scalable single-cell RNA profiling of pooled CRISPR screens

We present a combinatorial indexing method, PerturbSci-Kinetics, for capturing whole transcriptomes, nascent transcriptomes and single guide RNA (sgRNA) identities across hundreds of genetic perturbations at the single-cell level. Profiling a pooled CRISPR screen targeting various biological processes, we show the gene expression regulation during RNA synthesis, processing and degradation, miRNA biogenesis and mitochondrial mRNA processing, systematically decoding the genome-wide regulatory network that underlies RNA temporal dynamics at scale.

DUSP6 inhibition overcomes neuregulin/HER3-driven therapy tolerance in HER2+ breast cancer

Despite clinical benefits of tyrosine kinase inhibitors (TKIs) in cancer, most tumors can reactivate proliferation under TKI therapy. Here we present transcriptional profiling of HER2+ breast cancer cells transitioning from dormant drug tolerant cells to re-proliferating cells under continuous HER2 inhibitor (HER2i) therapy. Focusing on phosphatases, expression of dual-specificity phosphatase DUSP6 was found inhibited in dormant cells, but strongly induced upon regrowth. DUSP6 expression also selectively associated with poor patient survival in HER2+ breast cancers. DUSP6 overexpression conferred apoptosis resistance, whereas its pharmacological blockade prevented therapy tolerance development under HER2i therapy. DUSP6 targeting also synergized with clinically used HER2i combination therapies. Mechanistically DUSP6 is a positive regulator of HER3 expression, and its impact on HER2i tolerance was mediated by neuregulin-HER3 axis. In vivo, genetic targeting of DUSP6 reduced tumor growth in brain metastasis model, whereas its pharmacological targeting induced synthetic lethal therapeutic effect in combination with HER2i. Collectively this work demonstrates that DUSP6 drives escape from HER2i-induced dormancy, and that DUSP6 is a druggable target to overcome HER3-driven TKI resistance.

Lipidomic analysis identifies long-chain acylcarnitine as a target for ischemic stroke

INTRODUCTION

Lipid metabolism dysfunction is widely involved in the pathological process of acute ischemic stroke (AIS). The coordination of lipid metabolism between neurons and astrocytes is of great significance. However, the full scope of lipid dynamic changes and the function of key lipids during AIS remain unknown. Hence, identifying lipid alterations and characterizing their key roles in AIS is of great importance.

METHODS

Untargeted and targeted lipidomic analyses were applied to profile lipid changes in the ischemic penumbra and peripheral blood of transient middle cerebral artery occlusion (tMCAO) mice as well as the peripheral blood of AIS patients. Infarct volume and neurological deficits were assessed after tMCAO. The cell viability and dendritic complexity of primary neurons were evaluated by CCK8 assay and Sholl analysis. Seahorse, MitoTracker Green, tetramethyl rhodamine methyl ester (TMRM), 2',7'-dichlorodihydrofluorescein diacetate (DCFH-DA) and MitoSOX were used as markers of mitochondrial health. Fluorescent and isotopic free fatty acid (FFA) pulse-chase assays were used to track FFA flux in astrocytes.

RESULTS

Long-chain acylcarnitines (LCACs) were the lipids with the most dramatic changes in the ischemic penumbra and peripheral blood of tMCAO mice. LCACs were significantly elevated on admission in AIS patients and associated with poor outcomes in AIS patients. Increasing LCACs through a bolus administration of palmitoylcarnitine amplified stroke injury, while decreasing LCACs by overexpressing carnitine palmitoyltransferase 2 (CPT2) ameliorated stroke injury. Palmitoylcarnitine aggravated astrocytic mitochondrial damage after OGD/R, while CPT2 overexpression in astrocytes ameliorated cocultured neuron viability. Further study revealed that astrocytes stimulated by OGD/R liberated FFAs from lipid droplets into mitochondria to form LCACs, resulting in mitochondrial damage and lowered astrocytic metabolic support and thereby aggravated neuronal damage.

CONCLUSION

LCACs could accumulate and damage neurons by inducing astrocytic mitochondrial dysfunction in AIS. LCACs play a crucial role in the pathology of AIS and are novel promising diagnostic and prognostic biomarkers for AIS.

Apoptosis and NETotic cell death affect diabetic nephropathy independently: An study integrative study encompassing bioinformatics, machine learning, and experimental validation

OBJECTIVE

Although programmed cell death (PCD) and diabetic nephropathy (DN) are intrinsically conneted, the interplay among various PCD forms remains elusive. In this study, We aimed at identifying independently DN-associated PCD pathways and biomarkers relevant to the related pathogenesis.

METHODS

We acquired DN-related datasets from the GEO database and identified PCDs independently correlated with DN (DN-PCDs) through single-sample Gene Set Enrichment Analysis (ssGSEA) as well as, univariate and multivariate logistic regression analyses. Subsequently, applying differential expression analysis, weighted gene co-expression network analysis (WGCNA), and Mfuzz cluster analysis, we filtered the DN-PCDs pertinent to DN onset and progression. The convergence of various machine learning techniques ultimately spotlighted hub genes, substantiated through dataset meta-analyses and experimental validations, thereby confirming hub genes and related pathways expression consistencies.

RESULTS

We harmonized four DN-related datasets (GSE1009, GSE142025, GSE30528, and GSE30529) post-batch-effect removal for subsequent analyses. Our differential expression analysis yielded 709 differentially expressed genes (DEGs), comprising 446 upregulated and 263 downregulated DEGs. Based on our ssGSEA as well as univariate and multivariate logistic regressions, apoptosis and NETotic cell death were appraised as independent risk factors for DN (Odds Ratio > 1, p < 0.05). Next, we further refined 588 apoptosis- and NETotic cell death-associated genes through WGCNA and Mfuzz analysis, resulting in the identification of 17 DN-PCDs. Integrating protein-protein interaction (PPI) network analyses, network topology, and machine learning, we pinpointed hub genes (e.g., IL33, RPL11, and CX3CR1) as significant DN risk factors with expression corroborating in subsequent meta-analyses and experimental validations. Our GSEA enrichment analysis discerned differential enrichments between DN and control samples within pathways such as IL2/STAT5, IL6/JAK/STAT3, TNF-α via NF-κB, apoptosis, and oxidative phosphorylation, with related proteins such as IL2, IL6, and TNFα, which we subsequently submitted to experimental verification.

CONCLUSION

Innovatively stemming from from PCD interactions, in this study, we discerned PCDs with an independent impact on DN: apoptosis and NETotic cell death. We further screened DN evolution- and progression-related biomarkers, i.e. IL33, RPL11, and CX3CR1, all of which we empirically validated. This study not only poroposes a PCD-centric perspective for DN studies but also provides evidence for PCD-mediated immune cell infiltration exploration in DN regulation. Our results could motivate further exploration of DN pathogenesis, such as how the inflammatory microenvironment mediates NETotic cell death in DN regulation, representing a promising direction for future research.

Dynamic Alternative Polyadenylation during Litopenaeus Vannamei Metamorphosis Development

As an important mechanism in the post-transcriptional regulation of eukaryotic gene expression, alternative polyadenylation (APA) plays a key role in biological processes such as cell proliferation and differentiation. However, the role and dynamic pattern of APA during Litopenaeus vannamei metamorphosis are poorly understood. Here, RNA-seq data covering from the embryo to the maturation (16 time points) of L. vannamei were utilized. We identified 247 differentially expressed APA events between early and adult stages, and through fuzzy mean clustering analysis, we discovered five dynamic APA patterns. Among them, the gradual elongation of the 3'UTR is the major APA pattern that changes over time, and its genes are enriched in the pathways of protein and energy metabolism. Finally, we constructed mRNA-miRNA and PPI networks and detected several central miRNAs that may regulate L. vannamei development. Our results revealed the complex APA mechanisms in L. vannamei metamorphosis, shedding new light on post-transcriptional regulation of crustacean metamorphosis.

Spontaneous fermentation mitigates the frequency of genes encoding antimicrobial resistance spreading from the phyllosphere reservoir to the diet

The phyllosphere microbiome of vegetable products constitutes an important reservoir for multidrug resistant bacteria and Antibiotic Resistance Genes (ARG). Vegetable products including fermented products such as Paocai therefore may serve as a shuttle for extrinsic microorganisms with ARGs into the gut of consumers. Here we study the effect of fermentation on Paocai ARG dissemination by metagenomic analysis. Microbial abundance and diversity of the Paocai microbiome were diminished during fermentation, which correlated with the reduction of abundance in ARGs. Specifically, as fermentation progressed, Enterobacterales overtook Pseudomonadales as the predominant ARG carriers, and Lactobacillales and Enterobacteriales became the determinants of Paocai resistome variation. Moreover, the dual effect of microbes and metal resistance genes (MRGs) was the major contributor driving Paocai resistome dynamics. We recovered several metagenome-assembled genomes (MAGs) carrying acquired ARGs in the phyllosphere microbiome. ARGs of potential clinical and epidemiological relevance such as tet M and emrB-qacA, were mainly hosted by non-dominant bacterial genera. Overall, our study provides evidence that changes in microbial community composition by fermentation aid in constraining ARG dispersal from raw ingredients to the human microbiome but does not eliminate them.

Profiling the interplay and coevolution of Microcystis aeruginosa and cyanosiphophage Mic1

The cyanosiphophage Mic1 specifically infects the bloom-forming Microcystis aeruginosa FACHB 1339 from Lake Chaohu, China. Previous genomic analysis showed that its 92,627 bp double-stranded DNA genome consists of 98 putative open reading frames, 63% of which are of unknown function. Here, we investigated the transcriptome dynamics of Mic1 and its host using RNA sequencing. In the early, middle, and late phases of the 10 h lytic cycle, the Mic1 genes are sequentially expressed and could be further temporally grouped into two distinct clusters in each phase. Notably, six early genes, including gp49 that encodes a TnpB-like transposase, immediately reach the highest transcriptional level in half an hour, representing a pioneer cluster that rapidly regulates and redirects host metabolism toward the phage. An in-depth analysis of the host transcriptomic profile in response to Mic1 infection revealed significant upregulation of a polyketide synthase pathway and a type III-B CRISPR system, accompanied by moderate downregulation of the photosynthesis and key metabolism pathways. The constant increase of phage transcripts and relatively low replacement rate over the host transcripts indicated that Mic1 utilizes a unique strategy to gradually take over a small portion of host metabolism pathways after infection. In addition, genomic analysis of a less-infective Mic1 and a Mic1-resistant host strain further confirmed their dynamic interplay and coevolution via the frequent horizontal gene transfer. These findings provide insights into the mutual benefit and symbiosis of the highly polymorphic cyanobacteria M. aeruginosa and cyanophages.

IMPORTANCE

The highly polymorphic Microcystis aeruginosa is one of the predominant bloom-forming cyanobacteria in eutrophic freshwater bodies and is infected by diverse and abundant cyanophages. The presence of a large number of defense systems in M. aeruginosa genome suggests a dynamic interplay and coevolution with the cyanophages. In this study, we investigated the temporal gene expression pattern of Mic1 after infection and the corresponding transcriptional responses of its host. Moreover, the identification of a less-infective Mic1 and a Mic1-resistant host strain provided the evolved genes in the phage-host coevolution during the multiple-generation cultivation in the laboratory. Our findings enrich the knowledge on the interplay and coevolution of M. aeruginosa and its cyanophages and lay the foundation for the future application of cyanophage as a potential eco-friendly and bio-safe agent in controlling the succession of harmful cyanobacterial blooms.

A Dose-Response Study on Functional and Transcriptomic Effects of FSH on Ex Vivo Mouse Folliculogenesis

Follicle-stimulating hormone (FSH) binds to its membrane receptor (FSHR) in granulosa cells to activate various signal transduction pathways and drive the gonadotropin-dependent phase of folliculogenesis. Both FSH insufficiency (due to genetic or nongenetic factors) and FSH excess (as encountered with ovarian stimulation in assisted reproductive technology [ART]) can cause poor female reproductive outcomes, but the underlying molecular mechanisms remain elusive. Herein, we conducted single-follicle and single-oocyte RNA sequencing analysis along with other approaches in an ex vivo mouse folliculogenesis and oogenesis system to investigate the effects of different concentrations of FSH on key follicular events. Our study revealed that a minimum FSH threshold is required for follicle maturation into the high estradiol-secreting preovulatory stage, and such threshold is moderately variable among individual follicles between 5 and 10 mIU/mL. FSH at 5, 10, 20, and 30 mIU/mL induced distinct expression patterns of follicle maturation-related genes, follicular transcriptomics, and follicular cAMP levels. RNA sequencing analysis identified FSH-stimulated activation of G proteins and downstream canonical and novel signaling pathways that may critically regulate follicle maturation, including the cAMP/PKA/CREB, PI3K/AKT/FOXO1, and glycolysis pathways. High FSH at 20 and 30 mIU/mL resulted in noncanonical FSH responses, including premature luteinization, high production of androgen and proinflammatory factors, and reduced expression of energy metabolism-related genes in oocytes. Together, this study improves our understanding of gonadotropin-dependent folliculogenesis and provides crucial insights into how high doses of FSH used in ART may impact follicular health, oocyte quality, pregnancy outcome, and systemic health.

Therapeutic application of circular RNA aptamers in a mouse model of psoriasis

Efforts to advance RNA aptamers as a new therapeutic modality have been limited by their susceptibility to degradation and immunogenicity. In a previous study, we demonstrated synthesized short double-stranded region-containing circular RNAs (ds-cRNAs) with minimal immunogenicity targeted to dsRNA-activated protein kinase R (PKR). Here we test the therapeutic potential of ds-cRNAs in a mouse model of imiquimod-induced psoriasis. We find that genetic supplementation of ds-cRNAs leads to inhibition of PKR, resulting in alleviation of downstream interferon-α and dsRNA signals and attenuation of psoriasis phenotypes. Delivery of ds-cRNAs by lipid nanoparticles to the spleen attenuates PKR activity in examined splenocytes, resulting in reduced epidermal thickness. These findings suggest that ds-cRNAs represent a promising approach to mitigate excessive PKR activation for therapeutic purposes.

Transcriptome analysis of Kluyveromyces marxianus under succinic acid stress and development of robust strains

Kluyveromyces marxianus has become an attractive non-conventional yeast cell factory due to its advantageous properties such as high thermal tolerance and rapid growth. Succinic acid (SA) is an important platform molecule that has been applied in various industries such as food, material, cosmetics, and pharmaceuticals. SA bioproduction may be compromised by its toxicity. Besides, metabolite-responsive promoters are known to be important for dynamic control of gene transcription. Therefore, studies on global gene transcription under various SA concentrations are of great importance. Here, comparative transcriptome changes of K. marxianus exposed to various concentrations of SA were analyzed. Enrichment and analysis of gene clusters revealed repression of the tricarboxylic acid cycle and glyoxylate cycle, also activation of the glycolysis pathway and genes related to ergosterol synthesis. Based on the analyses, potential SA-responsive promoters were investigated, among which the promoter strength of IMTCP2 and KLMA_50231 increased 43.4% and 154.7% in response to 15 g/L SA. In addition, overexpression of the transcription factors Gcr1, Upc2, and Ndt80 significantly increased growth under SA stress. Our results benefit understanding SA toxicity mechanisms and the development of robust yeast for organic acid production. KEY POINTS: • Global gene transcription of K. marxianus is changed by succinic acid (SA) • Promoter activities of IMTCP2 and KLMA_50123 are regulated by SA • Overexpression of Gcr1, Upc2, and Ndt80 enhanced SA tolerance.

One-pot trimodal mapping of unmethylated, hydroxymethylated, and open chromatin sites unveils distinctive 5hmC roles at dynamic chromatin loci

We present a method, named Mx-TOP, for profiling of three epigenetic regulatory layers-chromatin accessibility, general DNA modification, and DNA hydroxymethylation-from a single library. The approach is based on chemo-enzymatic covalent tagging of unmodified CG sites and hydroxymethylated cytosine (5hmC) along with GC sites in chromatin, which are then mapped using tag-selective base-resolution TOP-seq sequencing. Our in-depth validation of the approach revealed its sensitivity and informativity in evaluating chromatin accessibility and DNA modification interactions that drive transcriptional regulation. We employed the technology in a study of chromatin and DNA demethylation dynamics during in vitro neuronal differentiation. The study highlighted the involvement of gene body 5hmC in modulating an extensive decoupling between promoter accessibility and transcription. The importance of 5hmC in chromatin remodeling was further demonstrated by the observed resistance of the developmentally acquired open loci to the global 5hmC erasure in neuronal progenitors.

Mesenchyme instructs growth while epithelium directs branching in the mouse mammary gland

The mammary gland is a unique organ that undergoes dynamic alterations throughout a female's reproductive life, making it an ideal model for developmental, stem cell and cancer biology research. Mammary gland development begins in utero and proceeds via a quiescent bud stage before the initial outgrowth and subsequent branching morphogenesis. How mammary epithelial cells transit from quiescence to an actively proliferating and branching tissue during embryogenesis and, importantly, how the branch pattern is determined remain largely unknown. Here, we provide evidence indicating that epithelial cell proliferation and onset of branching are independent processes, yet partially coordinated by the Eda signaling pathway. Through heterotypic and heterochronic epithelial-mesenchymal recombination experiments between mouse mammary and salivary gland tissues and ex vivo live imaging, we demonstrate that unlike previously concluded, the mode of branching is an intrinsic property of the mammary epithelium whereas the pace of growth and the density of ductal tree are determined by the mesenchyme. Transcriptomic profiling and ex vivo and in vivo functional studies in mice disclose that mesenchymal Wnt/ß-catenin signaling, and in particular IGF-1 downstream of it critically regulate mammary gland growth. These results underscore the general need to carefully deconstruct the different developmental processes producing branched organs.

Investigating HMGB1 as a potential serum biomarker for early diabetic nephropathy monitoring by quantitative proteomics

Current diagnostic methods for diabetic nephropathy (DN) lack precision, especially in early stages and monitoring progression. This study aims to find potential biomarkers for DN progression and evaluate their accuracy. Using serum samples from healthy controls (NC), diabetic patients (DM), early-medium stage DN (DN-EM), and late-stage DN (DN-L), researchers employed quantitative proteomics and Mfuzz clustering analysis revealed 15 proteins showing increased expression during DN progression, hinting at their biomarker potential. Combining Mfuzz clustering with weighted gene co-expression network analysis (WGCNA) highlighted five candidates (HMGB1, CD44, FBLN1, PTPRG, and ADAMTSL4). HMGB1 emerged as a promising biomarker, closely correlated with renal function changes. Experimental validation supported HMGB1's upregulation under high glucose conditions, reinforcing its potential as an early detection biomarker for DN. This research advances DN understanding and identifies five potential biomarkers, notably HMGB1, as a promising early monitoring target. These findings set the stage for future clinical diagnostic applications in DN.

The Mediterranean mussel Mytilus galloprovincialis: a novel model for developmental studies in mollusks

A model organism in developmental biology is defined by its experimental amenability and by resources created for the model system by the scientific community. For the most powerful invertebrate models, the combination of both has already yielded a thorough understanding of developmental processes. However, the number of developmental model systems is still limited, and their phylogenetic distribution heavily biased. Members of one of the largest animal lineages, the Spiralia, for example, have long been neglected. In order to remedy this shortcoming, we have produced a detailed developmental transcriptome for the bivalve mollusk Mytilus galloprovincialis, and have expanded the list of experimental protocols available for this species. Our high-quality transcriptome allowed us to identify transcriptomic signatures of developmental progression and to perform a first comparison with another bivalve mollusk: the Pacific oyster Crassostrea gigas. To allow co-labelling studies, we optimized and combined protocols for immunohistochemistry and hybridization chain reaction to create high-resolution co-expression maps of developmental genes. The resources and protocols described here represent an enormous boost for the establishment of Mytilus galloprovincialis as an alternative model system in developmental biology.

Transcriptomic Analysis Provides Insights into Candidate Genes and Molecular Pathways Involved in Growth of Mytilus coruscus Larvae

Growth is a fundamental aspect of aquaculture breeding programs, pivotal for successful cultivation. Understanding the mechanisms that govern growth and development differences across various stages can significantly boost seedling production of economically valuable species, thereby enhancing aquaculture efficiency and advancing the aquaculture industry. Mytilus coruscus, a commercially vital marine bivalve, underscores this importance. To decipher the intricate molecular mechanisms dictating growth and developmental disparities in marine shellfish, we conducted transcriptome sequencing and meticulously analyzed gene expression variations and molecular pathways linked to growth traits in M. coruscus. This study delved into the molecular and gene expression variations across five larval development stages, with a specific focus on scrutinizing the differential expression patterns of growth-associated genes using RNA sequencing and quantitative real-time PCR analysis. A substantial number of genes-36,044 differentially expressed genes (DEGs)-exhibited significant differential expression between consecutive developmental stages. These DEGs were then categorized into multiple pathways (Q value < 0.05), including crucial pathways such as the spliceosome, vascular smooth muscle contraction, DNA replication, and apoptosis, among others. In addition, we identified two pivotal signaling pathways-the Hedgehog (Hh) signaling pathway and the TGF-beta (TGF-β) signaling pathway-associated with the growth and development of M. coruscus larvae. Ten key growth-related genes were pinpointed, each playing crucial roles in molecular function and the regulation of growth traits in M. coruscus. These genes and pathways associated with growth provide deep insights into the molecular basis of physiological adaptation, metabolic processes, and growth variability in marine bivalves.

Fibronectin mediates activin A-promoted human trophoblast migration and acquisition of endothelial-like phenotype

BACKGROUND

During human early placentation, a proportion of extravillous trophoblasts (EVTs) migrate to the maternal decidua, differentiating into endovascular EVTs to remodel spiral arteries and ensure the establishment of blood circulation at the maternal-fetal interface. Inadequate EVT migration and endovascular differentiation are closely associated with adverse pregnancy outcomes such as miscarriage. Activin A and fibronectin are both secretory molecules abundantly expressed at the maternal-fetal interface. Activin A has been reported to regulate EVT biological functions. However, whether fibronectin mediates activin A-promoted EVT migration and acquisition of endothelial-like phenotype as well as the underlying molecular mechanisms remain unknown. Additionally, the role of fibronectin in pregnancy establishment and maintenance warrants further investigation.

METHODS

Primary and immortalized (HTR8/SVneo) human EVTs were used as in vitro study models. Cultured human first-trimester chorionic villous explants were utilized for ex vivo validation. A local fibronectin knockdown model in ICR mouse uteri, achieved by nonviral in vivo transfection with small interfering RNA (siRNA) targeting fibronectin 1 (si-Fn1), was employed to explore the roles of fibronectin in the establishment and maintenance of early pregnancy.

RESULTS

Our results showed that activin A treatment significantly induced fibronectin 1 (FN1) mRNA expression and fibronectin protein production, which is essential for human trophoblast migration and endothelial-like tube formation. Both basal and activin A-upregulated fibronectin expression were abolished by the TGF-β type I receptor inhibitor SB431542 or siRNA-mediated knockdown of activin receptor-like kinase (ALK4) or SMAD4. Moreover, activin A-increased trophoblast migration and endothelial-like tube formation were attenuated following the depletion of fibronectin. Fibronectin knockdown via intrauterine siRNA administration reduced CD31 and cytokeratin 8 (CK8) expression at the maternal-fetal interface, resulting in a decrease in the number of implantation sites and embryos.

CONCLUSIONS

Our study demonstrates that activin A promotes trophoblast cell migration and acquisition of endothelial-like phenotype via ALK4-SMAD2/3-SMAD4-mediated fibronectin upregulation. Furthermore, through a local fibronectin knockdown model in mouse uteri, we found that the absence of fibronectin at the maternal-fetal interface impedes endovascular migration of trophoblasts and decidual vascularization, thereby interfering with early embryo implantation and the maintenance of pregnancy. These findings provide novel insights into placental development during early pregnancy establishment and contribute to the advancement of therapeutic approaches for managing pregnancy complications related to trophoblast dysfunction.

Longitudinal single cell atlas identifies complex temporal relationship between type I interferon response and COVID-19 severity

Due to the paucity of longitudinal molecular studies of COVID-19, particularly those covering the early stages of infection (Days 1-8 symptom onset), our understanding of host response over the disease course is limited. We perform longitudinal single cell RNA-seq on 286 blood samples from 108 age- and sex-matched COVID-19 patients, including 73 with early samples. We examine discrete cell subtypes and continuous cell states longitudinally, and we identify upregulation of type I IFN-stimulated genes (ISGs) as the predominant early signature of subsequent worsening of symptoms, which we validate in an independent cohort and corroborate by plasma markers. However, ISG expression is dynamic in progressors, spiking early and then rapidly receding to the level of severity-matched non-progressors. In contrast, cross-sectional analysis shows that ISG expression is deficient and IFN suppressors such as SOCS3 are upregulated in severe and critical COVID-19. We validate the latter in four independent cohorts, and SOCS3 inhibition reduces SARS-CoV-2 replication in vitro. In summary, we identify complexity in type I IFN response to COVID-19, as well as a potential avenue for host-directed therapy.

Cuscuta campestris fine-tunes gene expression during haustoriogenesis as an adaptation to different hosts

The Cuscuta genus comprises obligate parasitic plants that have an unusually wide host range. Whether Cuscuta uses different infection strategies for different hosts or whether the infection strategy is mechanistically and enzymatically conserved remains unknown. To address this, we investigated molecular events during the interaction between field dodder (Cuscuta campestris) and two host species of the Solanum genus that are known to react differently to parasitic infection. We found that host gene induction, particularly of cell wall fortifying genes, coincided with a differential induction of genes for cell wall degradation in the parasite in the cultivated tomato (Solanum lycopersicum) but not in a wild relative (Solanum pennellii). This indicates that the parasite can adjust its gene expression in response to its host. This idea was supported by the increased expression of C. campestris genes encoding an endo-β-1,4-mannanase in response to exposure of the parasite to purified mono- and polysaccharides in a host-independent infection system. Our results suggest multiple key roles of the host cell wall in determining the outcome of an infection attempt.

JSI-124 inhibits cell proliferation and tumor growth by inducing autophagy and apoptosis in murine malignant mesothelioma

Malignant pleural mesothelioma (MPM), mainly caused by asbestos exposure, has a poor prognosis and lacks effective treatment compared with other cancer types. The intracellular transcription factor signal transducer and activator of transcription 3 (STAT3) is overexpressed and hyperactivated in most human cancers. In this study, the role of STAT3 in murine MPM was examined. Inhibition of the Janus kinase 2 (JAK2)/STAT3 pathway with the selective inhibitor JSI-124 has an antitumor effect in murine MPM. Specifically, we demonstrated that JSI-124 inhibits murine MPM cell growth and induces apoptotic and autophagic cell death. Exposure of RN5 and AB12 cells to JSI-124 resulted in apoptosis via the Bcl-2 family of proteins. JSI-124 triggered autophagosome formation, accumulation, and conversion of LC3I to LC3II. Autophagy inhibitors, Chloroquine (CQ) and Bafilomycin A1 (Baf-A1), inhibited autophagy and sensitized RN5 and AB12 cells to JSI-124-induced apoptosis. Our data indicate that JSI-124 is a promising therapeutic agent for MPM treatment.

On-capillary alkylation micro-reactor: a facile strategy for proteo-metabolome profiling in the same single cells

Single-cell multi-omics analysis can provide comprehensive insights to study cell-to-cell heterogeneity in normal and disease physiology. However, due to the lack of amplification technique, the measurement of proteome and metabolome in the same cell is challenging. Herein, a novel on-capillary alkylation micro-reactor (OCAM) was developed to achieve proteo-metabolome profiling in the same single cells, by which proteins were first covalently bound to an iodoacetic acid functionalized open-tubular capillary micro-reactor via sulfhydryl alkylation reaction, and metabolites were rapidly eluted, followed by on-column digestion of captured proteins. Compared with existing methods for low-input proteome sample preparation, OCAM exhibited improved efficiency, anti-interference ability and recovery, enabling the identification of an average of 1509 protein groups in single HeLa cells. This strategy was applied to single-cell proteo-metabolome analysis of mouse oocytes at different stages, 3457 protein groups and 171 metabolites were identified in single oocytes, which is the deepest coverage of proteome and metabolome from single mouse oocytes to date, achieving complementary characterization of metabolic patterns during oocyte maturation.

A transcriptional atlas identifies key regulators and networks for the development of spike tissues in barley

Grain number and size determine grain yield in crops and are closely associated with spikelet fertility and grain filling in barley (Hordeum vulgare). Abortion of spikelet primordia within individual barley spikes causes a 30%-50% loss in the potential number of grains during development from the awn primordium stage to the tipping stage, after that grain filling is the primary factor regulating grain size. To identify transcriptional signatures associated with spike development, we use a six-rowed barley cultivar (Morex) to develop a spatiotemporal transcriptome atlas containing 255 samples covering 17 stages and 5 positions along the spike. We identify several fundamental regulatory networks, in addition to key regulators of spike development and morphology. Specifically, we show HvGELP96, encoding a GDSL domain-containing protein, as a regulator of spikelet fertility and grain number. Our transcriptional atlas offers a powerful resource to answer fundamental questions in spikelet development and degeneration in barley.

The regulatory landscape of 5' UTRs in translational control during zebrafish embryogenesis

The 5' UTRs of mRNAs are critical for translation regulation, but their in vivo regulatory features are poorly characterized. Here, we report the regulatory landscape of 5' UTRs during early zebrafish embryogenesis using a massively parallel reporter assay of 18,154 sequences coupled to polysome profiling. We found that the 5' UTR is sufficient to confer temporal dynamics to translation initiation, and identified 86 motifs enriched in 5' UTRs with distinct ribosome recruitment capabilities. A quantitative deep learning model, DaniO5P, revealed a combined role for 5' UTR length, translation initiation site context, upstream AUGs and sequence motifs on in vivo ribosome recruitment. DaniO5P predicts the activities of 5' UTR isoforms and indicates that modulating 5' UTR length and motif grammar contributes to translation initiation dynamics. This study provides a first quantitative model of 5' UTR-based translation regulation in early vertebrate development and lays the foundation for identifying the underlying molecular effectors.

Achieving quantitative reproducibility in label-free multisite DIA experiments through multirun alignment

DIA is a mainstream method for quantitative proteomics, but consistent quantification across multiple LC-MS/MS instruments remains a bottleneck in parallelizing data acquisition. One reason for this inconsistency and missing quantification is the retention time shift which current software does not adequately address for runs from multiple sites. We present multirun chromatogram alignment strategies to map peaks across columns, including the traditional reference-based Star method, and two novel approaches: MST and Progressive alignment. These reference-free strategies produce a quantitatively accurate data-matrix, even from heterogeneous multi-column studies. Progressive alignment also generates merged chromatograms from all runs which has not been previously achieved for LC-MS/MS data. First, we demonstrate the effectiveness of multirun alignment strategies on a gold-standard annotated dataset, resulting in a threefold reduction in quantitation error-rate compared to non-aligned DIA results. Subsequently, on a multi-species dataset that DIAlignR effectively controls the quantitative error rate, improves precision in protein measurements, and exhibits conservative peak alignment. We next show that the MST alignment reduces cross-site CV by 50% for highly abundant proteins when applied to a dataset from 11 different LC-MS/MS setups. Finally, the reanalysis of 949 plasma runs with multirun alignment revealed a more than 50% increase in insulin resistance (IR) and respiratory viral infection (RVI) proteins, identifying 11 and 13 proteins respectively, compared to prior analysis without it. The three strategies are implemented in our DIAlignR workflow (>2.3) and can be combined with linear, non-linear, or hybrid pairwise alignment.

Isoform-resolved transcriptome of the human preimplantation embryo

Human preimplantation development involves extensive remodeling of RNA expression and splicing. However, its transcriptome has been compiled using short-read sequencing data, which fails to capture most full-length mRNAs. Here, we generate an isoform-resolved transcriptome of early human development by performing long- and short-read RNA sequencing on 73 embryos spanning the zygote to blastocyst stages. We identify 110,212 unannotated isoforms transcribed from known genes, including highly conserved protein-coding loci and key developmental regulators. We further identify 17,964 isoforms from 5,239 unannotated genes, which are largely non-coding, primate-specific, and highly associated with transposable elements. These isoforms are widely supported by the integration of published multi-omics datasets, including single-cell 8CLC and blastoid studies. Alternative splicing and gene co-expression network analyses further reveal that embryonic genome activation is associated with splicing disruption and transient upregulation of gene modules. Together, these findings show that the human embryo transcriptome is far more complex than currently known, and will act as a valuable resource to empower future studies exploring development.

RBPome identification in egg-cell like callus of Arabidopsis

RNA binding proteins (RBPs) have multiple and essential roles in transcriptional and posttranscriptional regulation of gene expression in all living organisms. Their biochemical identification in the proteome of a given cell or tissue requires significant protein amounts, which limits studies in rare and highly specialized cells. As a consequence, we know almost nothing about the role(s) of RBPs in reproductive processes such as egg cell development, fertilization and early embryogenesis in flowering plants. To systematically identify the RBPome of egg cells in the model plant Arabidopsis, we performed RNA interactome capture (RIC) experiments using the egg cell-like RKD2-callus and were able to identify 728 proteins associated with poly(A+)-RNA. Transcripts for 97 % of identified proteins could be verified in the egg cell transcriptome. 46 % of identified proteins can be associated with the RNA life cycle. Proteins involved in mRNA binding, RNA processing and metabolism are highly enriched. Compared with the few available RBPome datasets of vegetative plant tissues, we identified 475 egg cell-enriched RBPs, which will now serve as a resource to study RBP function(s) during egg cell development, fertilization and early embryogenesis. First candidates were already identified showing an egg cell-specific expression pattern in ovules.

Genomic and transcriptomic analyses reveal polygenic architecture for ecologically important traits in aspen (Populus tremuloides Michx.)

Intraspecific genetic variation in foundation species such as aspen (Populus tremuloides Michx.) shapes their impact on forest structure and function. Identifying genes underlying ecologically important traits is key to understanding that impact. Previous studies, using single-locus genome-wide association (GWA) analyses to identify candidate genes, have identified fewer genes than anticipated for highly heritable quantitative traits. Mounting evidence suggests that polygenic control of quantitative traits is largely responsible for this "missing heritability" phenomenon. Our research characterized the genetic architecture of 30 ecologically important traits using a common garden of aspen through genomic and transcriptomic analyses. A multilocus association model revealed that most traits displayed a highly polygenic architecture, with most variation explained by loci with small effects (likely below the detection levels of single-locus GWA methods). Consistent with a polygenic architecture, our single-locus GWA analyses found only 38 significant SNPs in 22 genes across 15 traits. Next, we used differential expression analysis on a subset of aspen genets with divergent concentrations of salicinoid phenolic glycosides (key defense traits). This complementary method to traditional GWA discovered 1243 differentially expressed genes for a polygenic trait. Soft clustering analysis revealed three gene clusters (241 candidate genes) involved in secondary metabolite biosynthesis and regulation. Our work reveals that ecologically important traits governing higher-order community- and ecosystem-level attributes of a foundation forest tree species have complex underlying genetic structures and will require methods beyond traditional GWA analyses to unravel.

An 11-point time course midgut transcriptome across 72 h after bloodfeeding provides detailed temporal resolution of transcript expression in the arbovirus vector, Aedes aegypti

As the major vector for dengue, Zika, yellow fever, and chikungunya viruses, the mosquito Aedes aegypti is one of the most important insects in public health. These viruses are transmitted by bloodfeeding, which is also necessary for the reproduction of the mosquito. Thus, the midgut plays an essential role in mosquito physiology as the center for bloodmeal digestion and as an organ that serves as the first line of defense against viruses. Despite its importance, transcriptomic dynamics with fine temporal resolution across the entire digestion cycle have not yet been reported. To fill this gap, we conducted a transcriptomic analysis of A. aegypti female midguts across a 72-h bloodmeal digestion cycle for 11 time points, with a particular focus on the first 24 h. PCA analysis confirmed that 72 h is indeed a complete digestion cycle. Cluster and GO enrichment analysis showed the orchestrated modulation of thousands of genes to accomplish the midgut's role as the center for digestion, as well as nutrient transport with a clear progression with sequential emphasis on transcription, translation, energy production, nutrient metabolism, transport, and finally, autophagy by 24-36 h. We further determined that many serine proteases are robustly expressed as if to prepare for unexpected physiological challenges. This study provides a powerful resource for the analysis of genomic features that coordinate the rapid and complex transcriptional program induced by mosquito bloodfeeding.

Genome-wide association study of multiyear dynamic growth traits in hybrid Liriodendron identifies robust genetic loci associated with growth trajectories

The genetic factors underlying growth traits differ over time points or stages. However, most current studies of phenotypes at single time points do not capture all loci or explain the genetic differences underlying growth trajectories. Hybrid Liriodendron exhibits obvious heterosis and is widely cultivated, although its complex genetic mechanism underlying growth traits remains unknown. A genome-wide association study (GWAS) is an effective method for elucidating the genetic architecture by identifying genetic loci underlying complex quantitative traits. In the present study, using a GWAS, we identified robust loci associated with growth trajectories in hybrid Liriodendron populations. We selected 233 hybrid progenies derived from 25 crosses for resequencing, and measured their tree height (H) and diameter at breast height (DBH) for 11 consecutive years; 192 972 high-quality single nucleotide polymorphisms (SNPs) were obtained. The dynamics of the multiyear single-trait GWAS showed that year-specific SNPs predominated, and only five robust SNPs for DBH were identified in at least three different years. Multitrait GWAS analysis with model parameters as latent variables also revealed 62 SNPs for H and 52 for DBH associated with the growth trajectory, displaying different biomass accumulation patterns, among which four SNPs exerted pleiotropic effects. All identified SNPs also exhibited temporal variations in effect sizes and inheritance patterns potentially related to different growth and developmental stages. The haplotypes resulting from these significant SNPs might pyramid favorable loci, benefitting the selection of superior genotypes. The present study provides insights into the genetic architecture of dynamic growth traits and lays a basis for future molecular-assisted breeding.

A surgically optimized intraoperative poly(I:C)-releasing hydrogel prevents cancer recurrence

Recurrences frequently occur following surgical removal of primary tumors. In many cancers, adjuvant therapies have limited efficacy. Surgery provides access to the tumor microenvironment, creating an opportunity for local therapy, in particular immunotherapy, which can induce local and systemic anti-cancer effects. Here, we develop a surgically optimized biodegradable hyaluronic acid-based hydrogel for sustained intraoperative delivery of Toll-like receptor 3 agonist poly(I:C) and demonstrate that it significantly reduces tumor recurrence after surgery in multiple mouse models. Mechanistically, poly(I:C) induces a transient interferon alpha (IFNα) response, reshaping the tumor/wound microenvironment by attracting inflammatory monocytes and depleting regulatory T cells. We demonstrate that a pre-existing IFN signature predicts response to the poly(I:C) hydrogel, which sensitizes tumors to immune checkpoint therapy. The safety, immunogenicity, and surgical feasibility are confirmed in a veterinary trial in canine soft tissue tumors. The surgically optimized poly(I:C)-loaded hydrogel provides a safe and effective approach to prevent cancer recurrence.

Improved RNA stability estimation indicates that transcriptional interference is frequent in diverse bacteria

We used stochastic simulations and experimental data from E. coli, K. aerogenes, Synechococcus PCC 7002 and Synechocystis PCC 6803 to provide evidence that transcriptional interference via the collision mechanism is likely a prevalent mechanism for bacterial gene regulation. Rifampicin time-series data can be used to globally monitor and quantify collision between sense and antisense transcription-complexes. Our findings also highlight that transcriptional events, such as differential RNA decay, partial termination, and internal transcriptional start sites often deviate from gene annotations. Consequently, within a single gene annotation, there exist transcript segments with varying half-lives and transcriptional properties. To address these complexities, we introduce 'rifi', an R-package that analyzes transcriptomic data from rifampicin time series. 'rifi' employs a dynamic programming-based segmentation approach to identify individual transcripts, enabling accurate assessment of RNA stability and detection of diverse transcriptional events.

3D chromatin reprogramming primes human memory TH2 cells for rapid recall and pathogenic dysfunction

Memory T cells provide long-lasting defense responses through their ability to rapidly reactivate, but how they efficiently "recall" an inflammatory transcriptional program remains unclear. Here, we show that human CD4⁺ memory T helper 2 (T_H2) cells carry a chromatin landscape synergistically reprogrammed at both one-dimensional (1D) and 3D levels to accommodate recall responses, which is absent in naive T cells. In memory T_H2 cells, recall genes were epigenetically primed through the maintenance of transcription-permissive chromatin at distal (super)enhancers organized in long-range 3D chromatin hubs. Precise transcriptional control of key recall genes occurred inside dedicated topologically associating domains ("memory TADs"), in which activation-associated promoter-enhancer interactions were preformed and exploited by AP-1 transcription factors to promote rapid transcriptional induction. Resting memory T_H2 cells from patients with asthma showed premature activation of primed recall circuits, linking aberrant transcriptional control of recall responses to chronic inflammation. Together, our results implicate stable multiscale reprogramming of chromatin organization as a key mechanism underlying immunological memory and dysfunction in T cells.

Crucial role of iron in epigenetic rewriting during adipocyte differentiation mediated by JMJD1A and TET2 activity

Iron metabolism is closely associated with the pathogenesis of obesity. However, the mechanism of the iron-dependent regulation of adipocyte differentiation remains unclear. Here, we show that iron is essential for rewriting of epigenetic marks during adipocyte differentiation. Iron supply through lysosome-mediated ferritinophagy was found to be crucial during the early stage of adipocyte differentiation, and iron deficiency during this period suppressed subsequent terminal differentiation. This was associated with demethylation of both repressive histone marks and DNA in the genomic regions of adipocyte differentiation-associated genes, 　including Pparg, which encodes PPARγ, the master regulator of adipocyte differentiation. In addition, we identified several epigenetic demethylases to be responsible for iron-dependent adipocyte differentiation, with the histone demethylase jumonji domain-containing 1A and the DNA demethylase ten-eleven translocation 2 as the major enzymes. The interrelationship between repressive histone marks and DNA methylation was indicated by an integrated genome-wide association analysis, and was also supported by the findings that both histone and DNA demethylation were suppressed by either the inhibition of lysosomal ferritin flux or the knockdown of iron chaperone poly(rC)-binding protein 2. In summary, epigenetic regulations through iron-dependent control of epigenetic enzyme activities play an important role in the organized gene expression mechanisms of adipogenesis.

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.

Chemical Inhibition of Bromodomain Proteins in Insect-Stage African Trypanosomes Perturbs Silencing of the Variant Surface Glycoprotein Repertoire and Results in Widespread Changes in the Transcriptome

The eukaryotic protozoan parasite Trypanosoma brucei is transmitted by the tsetse fly to both humans and animals, where it causes a fatal disease called African trypanosomiasis. While the parasite lacks canonical DNA sequence-specific transcription factors, it does possess histones, histone modifications, and proteins that write, erase, and read histone marks. Chemical inhibition of chromatin-interacting bromodomain proteins has previously been shown to perturb bloodstream specific trypanosome processes, including silencing of the variant surface glycoprotein (VSG) genes and immune evasion. Transcriptomic changes that occur in bromodomain-inhibited bloodstream parasites mirror many of the changes that occur as parasites developmentally progress from the bloodstream to the insect stage. We performed transcriptome sequencing (RNA-seq) time courses to determine the effects of chemical bromodomain inhibition in insect-stage parasites using the compound I-BET151. We found that treatment with I-BET151 causes large changes in the transcriptome of insect-stage parasites and also perturbs silencing of VSG genes. The transcriptomes of bromodomain-inhibited parasites share some features with early metacyclic-stage parasites in the fly salivary gland, implicating bromodomain proteins as important for regulating transcript levels for developmentally relevant genes. However, the downregulation of surface procyclin protein that typically accompanies developmental progression is absent in bromodomain-inhibited insect-stage parasites. We conclude that chemical modulation of bromodomain proteins causes widespread transcriptomic changes in multiple trypanosome life cycle stages. Understanding the gene-regulatory processes that facilitate transcriptome remodeling in this highly diverged eukaryote may shed light on how these mechanisms evolved. IMPORTANCE The disease African trypanosomiasis imposes a severe human and economic burden for communities in sub-Saharan Africa. The parasite that causes the disease is transmitted to the bloodstream of a human or ungulate via the tsetse fly. Because the environments of the fly and the bloodstream differ, the parasite modulates the expression of its genes to accommodate two different lifestyles in these disparate niches. Perturbation of bromodomain proteins that interact with histone proteins around which DNA is wrapped (chromatin) causes profound changes in gene expression in bloodstream-stage parasites. This paper reports that gene expression is also affected by chemical bromodomain inhibition in insect-stage parasites but that the genes affected differ depending on life cycle stage. Because trypanosomes diverged early from model eukaryotes, an understanding of how trypanosomes regulate gene expression may lend insight into how gene-regulatory mechanisms evolved. This could also be leveraged to generate new therapeutic strategies.

A cytosolic surveillance mechanism activates the mitochondrial UPR

The mitochondrial unfolded protein response (UPR^mt) is essential to safeguard mitochondria from proteotoxic damage by activating a dedicated transcriptional response in the nucleus to restore proteostasis^1,2. Yet, it remains unclear how the information on mitochondria misfolding stress (MMS) is signalled to the nucleus as part of the human UPR^mt (refs. ^3,4). Here, we show that UPR^mt signalling is driven by the release of two individual signals in the cytosol-mitochondrial reactive oxygen species (mtROS) and accumulation of mitochondrial protein precursors in the cytosol (c-mtProt). Combining proteomics and genetic approaches, we identified that MMS causes the release of mtROS into the cytosol. In parallel, MMS leads to mitochondrial protein import defects causing c-mtProt accumulation. Both signals integrate to activate the UPR^mt; released mtROS oxidize the cytosolic HSP40 protein DNAJA1, which leads to enhanced recruitment of cytosolic HSP70 to c-mtProt. Consequently, HSP70 releases HSF1, which translocates to the nucleus and activates transcription of UPR^mt genes. Together, we identify a highly controlled cytosolic surveillance mechanism that integrates independent mitochondrial stress signals to initiate the UPR^mt. These observations reveal a link between mitochondrial and cytosolic proteostasis and provide molecular insight into UPR^mt signalling in human cells.

TWEAK Signaling-Induced ID1 Expression Drives Malignant Transformation of Hepatic Progenitor Cells During Hepatocarcinogenesis

The malignant transformation of hepatic progenitor cells (HPCs) in the inflammatory microenvironment is the root cause of hepatocarcinogenesis. However, the potential molecular mechanisms are still elusive. The HPCs subgroup is identified by single-cell RNA (scRNA) sequencing and the phenotype of HPCs is investigated in the primary HCC model. Bulk RNA sequencing (RNA-seq) and proteomic analyses are also performed on HPC-derived organoids. It is found that tumors are formed from HPCs in peritumor tissue at the 16th week in a HCC model. Furthermore, it is confirmed that the macrophage-derived TWEAK/Fn14 promoted the expression of inhibitor of differentiation-1 (ID1) in HPCs via NF-κB signaling and a high level of ID1 induced aberrant differentiation of HPCs. Mechanistically, ID1 suppressed differentiation and promoted proliferation in HPCs through the inhibition of HNF4α and Rap1GAP transcriptions. Finally, scRNA sequencing of HCC patients and investigation of clinical specimens also verified that the expression of ID1 is correlated with aberrant differentiation of HPCs into cancer stem cells, patients with high levels of ID1 in HPCs showed a poorer prognosis. This study provides important intervention targets and a theoretical basis for the clinical diagnosis and treatment of HCC.

Seiðr: Efficient calculation of robust ensemble gene networks

Gene regulatory and gene co-expression networks are powerful research tools for identifying biological signal within high-dimensional gene expression data. In recent years, research has focused on addressing shortcomings of these techniques with regard to the low signal-to-noise ratio, non-linear interactions and dataset dependent biases of published methods. Furthermore, it has been shown that aggregating networks from multiple methods provides improved results. Despite this, few useable and scalable software tools have been implemented to perform such best-practice analyses. Here, we present Seidr (stylized Seiðr), a software toolkit designed to assist scientists in gene regulatory and gene co-expression network inference. Seidr creates community networks to reduce algorithmic bias and utilizes noise corrected network backboning to prune noisy edges in the networks. Using benchmarks in real-world conditions across three eukaryotic model organisms, Saccharomyces cerevisiae, Drosophila melanogaster, and Arabidopsis thaliana, we show that individual algorithms are biased toward functional evidence for certain gene-gene interactions. We further demonstrate that the community network is less biased, providing robust performance across different standards and comparisons for the model organisms. Finally, we apply Seidr to a network of drought stress in Norway spruce (Picea abies (L.) H. Krast) as an example application in a non-model species. We demonstrate the use of a network inferred using Seidr for identifying key components, communities and suggesting gene function for non-annotated genes.