Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2

Dissemination and persistence of antimicrobial resistance (AMR) along the wastewater-river continuum

Antimicrobial resistance (AMR) is a global health hazard. Although clinical and agricultural environments are well-established contributors to the evolution and dissemination of AMR, research on wastewater treatment works (WwTWs) has highlighted their potential role as disseminators of AMR in freshwater environments. Using metagenomic sequencing and analysis, we investigated the changes in resistomes and associated mobile genetic elements within untreated wastewater influents and treated effluents of five WwTWs, and sediments collected from corresponding river environments in Oxfordshire, UK, across three seasonal periods within a year. Our analysis demonstrated a high diversity and abundance of antimicrobial resistance genes (ARGs) in untreated wastewater influents, reflecting the varied anthropogenic and environmental origins of wastewater. WwTWs effectively reduced AMR in the final effluent, with an average 87 % reduction in normalised ARG abundance and an average 63 % reduction in richness. However, wastewater effluents significantly impacted the antimicrobial resistome of the receiving rivers, with an average 543 % increase in ARG abundance and a 164 % increase in richness from upstream sediments to downstream sediments. The normalised abundance of the human gut-associated bacteriophage crAssphage was highly associated with both ARG abundance and richness. We observed seasonal variation in the resistome of raw influent which was not found in the effluent-receiving sediments. We illustrate the potential of WwTWs as focal points for disseminating ARGs and resistance-selecting chemicals, contributing to the elevation of environmental AMR. Our study emphasises the need for a comprehensive understanding of the anthropogenic impacts on AMR evolution and dissemination in wastewater and river environments, informing efforts to mitigate this growing public health crisis.

Potential Contributions of Human Endogenous Retroviruses in Innate Immune Memory

The phenomenon wherein innate immune cells adopt long-term inflammatory phenotypes following the first stimuli is named trained immunity and can improve host defense against infections. Transcriptional and epigenetic reprogramming are critical mechanisms of trained immunity; however, the regulatory networks are not entirely clear at present. The human endogenous retroviruses (HERVs) provide large amounts of transcriptional regulators in the regulatory pathways. In this study, we analyzed published large omics data to explore the roles of such "dark matter" of the human genome in trained and tolerant macrophages. We collected 80 RNA sequencing data and 62 sequencing data to detect histone modifications and active regulatory regions from nine published studies on trained and tolerant macrophages. By analyzing the characteristics of transcription and epigenetic modification of HERVs, as well as their association with gene expression, we found that 15.3% of HERVs were transcribed nonrandomly from noncoding regions and enriched in specific HERV families and specific chromosomes, such as chromosomes 11, 15, 17, and 19, and they were highly related with the expression of adjacent genes. We found that 295 differentially expressed HERVs are located in 50-kbp flanking regions of 142 differentially expressed genes. We found epigenetic changes of these HERVs and that overlap with predicted enhancers and identified 35 enhancer-like HERVs. The related genes were highly involved in the activation and inflammatory responses, such as the TLR pathway. Other pathways including phosphoinositide signaling and transport of folate and K+ might be also related with trained immunity, which require further study. These results demonstrated that HERVs might play important roles in trained immunity.

YTHDC2 serves a distinct late role in spermatocytes during germ cell differentiation

Posttranscriptional regulation of gene expression by RNA-binding proteins can enhance the speed and robustness of cell state transitions by controlling RNA stability, localization, or if, when, or where mRNAs are translated. The RNA helicase YTHDC2 is required to shut down components of the mitotic program to facilitate a proper switch from mitosis to meiosis in mouse germ cells. Here, we show that YTHDC2 has a second essential role in promoting meiotic progression in late spermatocytes. Inducing conditional knockout of Ythdc2 during the first wave of spermatogenesis, after initiation of meiotic prophase, allowed YTHDC2-deficient germ cells to advance to the pachytene stage and properly express many meiotic markers. However, the YTHDC2-deficient spermatocytes mis-expressed a number of genes, some up-regulated and some down-regulated, failed to transition to the diplotene stage, and then quickly died. Coimmunoprecipitation experiments revealed that YTHDC2 interacts with several RNA-binding proteins in early or late spermatocytes, with many of the interacting proteins, including MEIOC, localizing to granules, similar to YTHDC2. Our findings suggest that YTHDC2 collaborates with other RNA granule components to facilitate proper progression of germ cells through multiple steps of meiosis via mechanisms influencing posttranscriptional regulation of RNAs.

The extension of mammalian pregnancy required taming inflammation: Independent evolution of extended placentation in the tammar wallaby

In the first live-bearing mammals, pregnancy was likely short and ended with a brief period of inflammatory maternal-fetal interaction. This mode of reproduction has been retained in many marsupials. While inflammation is key to successful implantation in eutherians, a key innovation in eutherians is the ability to switch off this inflammation after it has been initiated. This extended period, in which inflammation is suppressed, likely allowed for an extended period of placentation. Extended placentation has evolved independently in one lineage of marsupials, the macropodids (wallabies and kangaroos), with placentation lasting beyond the 2 to 4 d seen in other marsupial taxa, which allows us to investigate the role of inflammation response after attachment in the extension of placentation in mammals. By comparing gene expression changes at attachment in three marsupial species, the tammar wallaby, opossum, and fat-tailed dunnart, we show that inflammatory attachment is an ancestral feature of marsupial implantation. In contrast to eutherians, where attachment-related (quasi-) inflammatory reaction is even involved in epitheliochorial placentation (e.g., pig), this study found no evidence of a distinct attachment-related reaction in wallabies. Instead, only a small number of inflammatory genes are expressed at distinct points of gestation, including IL6 before attachment, LIF throughout placentation, and prostaglandins before birth. During parturition, a more distinct inflammatory reaction is detectable, likely involved in precipitating the parturition cascade similar to eutherians. We suggest that in wallaby, extended gestation became possible by avoiding an inflammatory attachment reaction, which is a different strategy than seen in eutherians.

Tumor microbiota of renal cell carcinoma affects clinical prognosis by influencing the tumor immune microenvironment

Despite reported influences of the intratumoral microbiome on cancer progression, its role in this subtype remains unclear. This study aimed to characterize the microbial landscape and signatures of kidney renal clear cell carcinoma using RNA-Seq data from The Cancer Genome Atlas. Following microbial decontamination, differential microbial analysis was conducted between tumorous and adjacent non-tumorous samples. Compared to non-tumorous samples, tumorous microbiota exhibited reduced α and β diversity and distinct phylum-level communities. Differential microbial analysis between patients exhibiting long and short overall survival revealed ten significant differential microbial genera, with six genera correlating with a positive prognosis (Plasmodium, Babesia, Toxoplasma, Cytobacillus, Alicyclobacillus, Verrucomicrobium) and four with a negative prognosis (Colletotrichum, Leuconostoc, Gluconobacter, and Parabacteroides). Employing Cox regression analysis and support vector machines, a prognosis-related microbiome risk signature was developed, achieving an AUC of 0.809. Based on this risk signature, two microbiome-based subtypes were found to be significantly associated with distinct clinical prognoses and immune microenvironments. These findings were corroborated by significant correlations between prognostic-relevant microorganisms and 30 immune-related differentially expressed genes. Specifically, microbial genera associated with a negative prognosis were linked to a pro-tumor acute inflammatory immune response, whereas genera related to a positive prognosis were associated with an anti-tumor adaptive immune response. In conclusion, microbiome-based subtyping revealed correlations between tumor microbiome, clinical prognosis, and tumor microenvironment, indicating intratumoral microbiota as a promising prognostic biomarker for kidney renal clear cell carcinoma.

Deciphering dynamic interactions between spermatozoa and the ovarian microenvironment through integrated multi-omics approaches in viviparous Sebastes schlegelii

The ovarian microenvironment plays a crucial role in ensuring the reproductive success of viviparous teleosts. However, the molecular mechanism underlying the interaction between spermatozoa and the ovarian microenvironment has remained elusive. This study aimed to contribute to a better understanding of this process in black rockfish (Sebastes schlegelii) using integrated multi-omics approaches. The results demonstrated significant upregulation of ovarian complement-related proteins and pattern recognition receptors, along with remodeling of glycans on the surface of spermatozoa at the early spermatozoa-storage stage (1 month after mating). As spermatozoa were stored over time, ovarian complement proteins were progressively repressed by tryptophan and hippurate, indicating a remarkable adaptation of spermatozoa to the ovarian microenvironment. Before fertilization, a notable upregulation of cellular junction proteins was observed. The study revealed that spermatozoa bind to ZPB2a protein through GSTM3 and that ZPB2a promotes spermatozoa survival and movement in a GSTM3-dependent manner. These findings shed light on a key mechanism that influences the dynamics of spermatozoa in the female reproductive tract, providing valuable insights into the molecular networks regulating spermatozoa adaptation and survival in species with internal fertilization.

Transcriptome analysis reveals key regulatory networks and genes involved in the acquisition of cold stress memory in pepper seedlings

Temperature is an important limiting factor in the counter-seasonal cultivation of pepper. Currently, there are no studies on transcriptomic analysis of 'cold stress memory' in pepper. In this study, in order to understand the mechanism of 'cold stress memory' in pepper (Capsicum annuum L.), seedlings were subjected to the following treatments: normal temperature treatment (P0), the first cold treatment for 3 days (P3), the recovery temperature treatment for 3 days (R3), and another cold treatment for 3 days (RP3). The results showed that P3 plants wilted the most, RP3 the second and R3 the least. Leaf reactive oxygen species (ROS) and electrolyte leakage were the most in P3, the second in RP3 and the least in R3. In addition, RP3 had the highest accumulation of zeaxanthin, violaxanthin and β-cryptoxanthin, followed by P3, and R3 had the least. These results suggest that pepper seedlings are characterized by 'cold stress memory'. Transcriptomics was used to analyze the key genes and transcription factors involved in the biosynthesis of zeaxanthin, violaxanthin and β-cryptoxanthin during the formation of 'cold stress memory'. This study provides candidate genes and transcription factors for an in-depth study of the cold tolerance mechanism in pepper.

Delivery of a Jagged1-PEG-MAL hydrogel with pediatric human bone cells regenerates critically sized craniofacial bone defects

Current treatments for congenital and acquired craniofacial (CF) bone abnormalities are limited and costly. Conventional methods involve surgical correction, short-term stabilization, and long-term bone grafting, which may include problematic allografts and limited autografts. While bone morphogenetic protein 2 (BMP2) has been used for bone regeneration, it can cause bone overgrowth and life-threatening inflammation. Bone marrow-derived mesenchymal stem cell therapies, though promising, are not Food and Drug Administration approved and are resource intensive. Thus, there is a need for effective, affordable, and less side-effect-prone bone regenerative therapies. Previous research demonstrated that JAGGED1 induces osteoblast commitment in murine cranial neural crest cells through a NOTCH-dependent non-canonical pathway involving JAK2-STAT5. We hypothesize that delivery of JAGGED1 and induction of its downstream NOTCH non-canonical signaling in pediatric human osteoblasts constitutes an effective bone regenerative treatment. Delivering pediatric human bone-derived osteoblast-like cells to an in vivo murine bone loss model of a critically sized cranial defect, we identified that JAGGED1 promotes human pediatric osteoblast commitment and bone formation through p70 S6K phosphorylation. This approach highlights the potential of JAGGED1 and its downstream activators as innovative treatments for pediatric CF bone loss.

Transcriptomic dynamics of ABA response in Brassica napus guard cells

Drought has a significant, negative impact on crop production; and these effects are poised to increase with climate change. Plants acclimate to drought and water stress through diverse physiological responses, primarily mediated by the hormone abscisic acid (ABA). Because plants lose the majority of their water through stomatal pores on aerial surfaces of plants, stomatal closure is one of the rapid responses mediated by ABA to reduce transpirational water loss. The dynamic changes in the transcriptome of stomatal guard cells in response to ABA have been investigated in the model plant Arabidopsis thaliana. However, guard cell transcriptomes have not been analyzed in agronomically valuable crops such as a major oilseed crop, rapeseed. In this study, we investigated the dynamics of ABA-regulated transcriptomes in stomatal guard cells of Brassica napus and conducted comparison analysis with the transcriptomes of A. thaliana. We discovered changes in gene expression indicating alterations in a host of physiological processes, including stomatal movement, metabolic reprogramming, and light responses. Our results suggest the existence of both immediate and delayed responses to ABA in Brassica guard cells. Furthermore, the transcription factors and regulatory networks mediating these responses are compared to those identified in Arabidopsis. Our results imply the continuing evolution of ABA responses in Brassica since its divergence from a common ancestor, involving both protein-coding and non-coding nucleotide sequences. Together, our results will provide a basis for developing strategies for molecular manipulation of drought tolerance in crop plants.

Modulating bacterial function utilizing A knowledge base of transcriptional regulatory modules

Synthetic biology enables the reprogramming of cellular functions for various applications. However, challenges in scalability and predictability persist due to context-dependent performance and complex circuit-host interactions. This study introduces an iModulon-based engineering approach, utilizing machine learning-defined co-regulated gene groups (iModulons) as design parts containing essential genes for specific functions. This approach identifies the necessary components for genetic circuits across different contexts, enhancing genome engineering by improving target selection and predicting module behavior. We demonstrate several distinct uses of iModulons: (i) discovery of unknown iModulons to increase protein productivity, heat tolerance and fructose utilization; (ii) an iModulon boosting approach, which amplifies the activity of specific iModulons, improved cell growth under osmotic stress with minimal host regulation disruption; (iii) an iModulon rebalancing strategy, which adjusts the activity levels of iModulons to balance cellular functions, significantly increased oxidative stress tolerance while minimizing trade-offs and (iv) iModulon-based gene annotation enabled natural competence activation by predictably rewiring iModulons. Comparative experiments with traditional methods showed our approach offers advantages in efficiency and predictability of strain engineering. This study demonstrates the potential of iModulon-based strategies to systematically and predictably reprogram cellular functions, offering refined and adaptable control over complex regulatory networks.

RNAi of nuclear receptor E78 inhibits the cuticle formation in the molting process of spider mite, Tetranychus urticae

BACKGROUND

The two-spotted spider mite, Tetranychus urticae, is an important pest mite in agriculture worldwide. E78, as a member of the nuclear receptor superfamily and a downstream responsive gene of ecdysteroids, plays a crucial role in regulating physiological behaviors such as development and reproduction in insects. However, its function in mites remains unclear. The aim of this study was to explore how E78 functions in the molting process of spider mites.

RESULTS

In this study, reverse transcription quantitative polymerase chain reaction (RT-qPCR) experiments to analyze the expression pattern of TuE78 during the development of Tetranychus urticae, demonstrated that the expression level of TuE78 was higher during the molting state than that after the completion of molting, and it reached a peak expression level when the deutonymph mites entered the molting stage. RNA interference (RNAi)-mediated gene-silencing of TuE78 resulted in 95% deutonymph mite molt failure. A series of analysis under a light microscope, and scanning and transmission electron microscopy revealed that RNAi mites died within the exuvium without ecdysis, and that apolysis had started but the new cuticle was thin and the typical cuticular lamellae were absent, indicating blockage of the post-apolysial processes and explaining molt failure. Hence, transcriptome sequencing confirmed that the expression of cuticle protein and lipid metabolism-related genes was significantly affected after TuE78 silencing.

CONCLUSION

This study demonstrated that TuE78 participates in the molting process of Tetranychus urticae by regulating the post-apolysial processes with the formation of new cuticle and successful ecdysis. This in turn suggests the potential of TuE78 as a target for pest mite control and provides a theoretical basis for further exploration of the molecular regulatory mechanism of spider mite molting. © 2024 Society of Chemical Industry. Published by John Wiley & Sons Ltd.

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

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

Dynamics of diversified A-to-I editing in Streptococcus pyogenes is governed by changes in mRNA stability

Adenosine-to-inosine (A-to-I) RNA editing plays an important role in the post-transcriptional regulation of eukaryotic cell physiology. However, our understanding of the occurrence, function and regulation of A-to-I editing in bacteria remains limited. Bacterial mRNA editing is catalysed by the deaminase TadA, which was originally described to modify a single tRNA in Escherichia coli. Intriguingly, several bacterial species appear to perform A-to-I editing on more than one tRNA. Here, we provide evidence that in the human pathogen Streptococcus pyogenes, tRNA editing has expanded to an additional tRNA substrate. Using RNA sequencing, we identified more than 27 editing sites in the transcriptome of S. pyogenes SF370 and demonstrate that the adaptation of S. pyogenes TadA to a second tRNA substrate has also diversified the sequence context and recoding scope of mRNA editing. Based on the observation that editing is dynamically regulated in response to several infection-relevant stimuli, such as oxidative stress, we further investigated the underlying determinants of editing dynamics and identified mRNA stability as a key modulator of A-to-I editing. Overall, our findings reveal the presence and diversification of A-to-I editing in S. pyogenes and provide novel insights into the plasticity of the editome and its regulation in bacteria.

Mitophagy Defects Exacerbate Inflammation and Aberrant Proliferation in Lymphocytic Thyroiditis

Background: Mitochondrial dysfunction in the thyroid due to defective mitophagy has been observed in lymphocytic thyroiditis (LT). However, the effect of impaired mitophagy on the pathogenesis of LT is not well understood. The aim of this study is to investigate the role of mitophagy dysregulation in the thyroid gland. Methods: We analyzed RNA sequencing data of human thyroid glands with/without LT from Genotype-Tissue Expression (GTEx; n = 653) and performed RNA sequencing in thyroid glands of phosphatase and tensin homolog-induced putative protein kinase 1 (Pink1) knock-out and wild-type mice. We evaluated the phenotypic and histopathologic characteristics of the human (n = 16) and mouse thyroids. Additionally, we assessed cell proliferation, reactive oxygen species (ROS) production, and cytokine secretion of human thyroid epithelial cells (HTori-3) treated with PINK1 siRNA or a mitophagy inhibitor. Results: We found that expression of PINK1, a key regulator of mitophagy, was compromised in human thyroids with LT. Thyroid glands of Pink1-deficient mice exhibited increased inflammatory responses and nodular hyperplasia. Furthermore, mitophagy defects led to the production of pro-inflammatory cytokines and ROS in thyroid cells, resulting in immune cell recruitment. Notably, these mitophagy defects upregulated both the RNA expression and protein secretion of amphiregulin (AREG), an epidermal growth factor receptor (EGFR) ligand, in thyroid cells, while decreasing the protein expression of cAMP response element-binding protein (CREB), a transcription factor that suppresses AREG transcription. Finally, we demonstrated that aberrant cell proliferation in thyroid cells, driven by mitophagy defects, was mitigated after treatment with cetuximab, an EGFR inhibitor. Conclusions: In this study, we observed that mitophagy defects in the thyroid not only intensify inflammation through the accumulation of ROS, cytokine production, and immune cell recruitment but also contribute to hyperplasia via the EGFR pathway, facilitated by increased secretion of AREG from thyroid cells.

Ruminant Trans Fatty Acid Intake Modulates Inflammation Pathways in the Adipose Tissue Transcriptome of C57BL/6 Mice

SCOPE

The study aims to analyze transcriptomic profiles in adipose tissues postconsumption of elaidic acid (EA; trans-18:1n-9) and trans-palmitoleic acid (TPA; trans-16:1n-7), elucidating their different effects on inflammation and glucose metabolism.

METHODS AND RESULTS

Twenty C57BL/6 mice are divided into four groups. Each group receives one of the following formulations in drinking water: lecithin nanovesicles, nanovesicles containing either lecithin with EA or TPA (86:14 w/w), or water (control) for 28 days with a regular fat diet (18% calories from fat). Total RNA is extracted, and paired-end sequencing is performed. TPA intake alters the expression of 351 genes compared to EA intake, including 11 downregulated and 340 upregulated genes (fold change [FC] >1.5, p < 0.05). TPA compares to EA upregulated: Slc5a8, Lcn2, Csf3, Scube1, Mapk13, Bdkrb2, Ctla2a, Slc2a1, Oas3, Cx3cl1, Oas2, Nlrp6, Pycard, Cyba, Ddr1, and Prkab1 and downregulated Fas gene. These genes are related to the NOD-like receptor, lipid and atherosclerosis, IL-17 signaling, TNF, nonalcoholic fatty liver disease, cytokine-cytokine receptor interaction, adipocytokine, glucagon, insulin resistance, and inflammatory mediator regulation of TRP channels signaling.

CONCLUSION

TPA intake has a distinct impact on the regulation of inflammation and diabetes-related pathways in adipose tissue compared to EA.

Evolutionary Conserved and Divergent Responses to Copper Zinc Superoxide Dismutase Inhibition in Plants

After an initial evolution in a reducing environment, life got successively challenged by reactive oxygen species (ROS), especially during the great oxidation event (GOE) that followed the development of photosynthesis. Therefore, ROS are deeply intertwined into the physiological, morphological and transcriptional responses of most present-day organisms. Copper-zinc superoxide dismutases (CuZnSODs) evolved during the GOE and are present in charophytes and extant land plants, but nearly absent from chlorophytes. The chemical inhibitor of CuZnSOD, lung cancer screen 1 (LCS-1), could greatly facilitate the study of SODs in diverse plants. Here, we determined the impact of chemical inhibition of plant CuZnSOD activity, on plant growth, transcription and metabolism. We followed a comparative approach by using different plant species, including Marchantia Polymorpha and Physcomitrium patens, representing bryophytes, the sister lineage to vascular plants, and Arabidopsis thaliana. We show that LCS-1 causes oxidative stress in plants and that the inhibition of CuZnSODs provoked a similar core response that mainly impacted glutathione homoeostasis in all plant species analysed. That said, Physcomitrium and Arabidopsis, which contain multiple CuZnSOD isoforms showed a more complex and exacerbated response. In addition, an untargeted metabolomics approach revealed a specific metabolic signature for each plant species. Our comparative analysis exposes a conserved core response at the physiological and transcriptional level towards LCS-1, while the metabolic response largely varies. These differences correlate with the number and localization of the CuZnSOD isoforms present in each species.

Effects of cadherin mediated contact normalization on oncogenic Src kinase mediated gene expression and protein phosphorylation

Nontransformed cells form heterotypic cadherin junctions with adjacent transformed cells to inhibit tumor cell growth and motility. Transformed cells must override this form of growth control, called "contact normalization", to invade and metastasize during cancer progression. Heterocellular cadherin junctions between transformed and nontransformed cells are needed for this process. However, specific mechanisms downstream of cadherin signaling have not been clearly elucidated. Here, we utilized a β-catenin reporter construct to determine if contact normalization affects Wnt signaling in transformed cells. β-catenin driven GFP expression in Src transformed mouse embryonic cells was decreased when cultured with cadherin competent nontransformed cells compared to transformed cells cultured with themselves, but not when cultured with cadherin deficient nontransformed cells. We also utilized a layered culture system to investigate the effects of oncogenic transformation and contact normalization on gene expression and oncogenic Src kinase mediated phosphorylation events. RNA-Seq analysis found that cadherin dependent contact normalization inhibited the expression of 22 transcripts that were induced by Src transformation, and increased the expression of 78 transcripts that were suppressed by Src transformation. Phosphoproteomic analysis of cells expressing a temperature sensitive Src kinase construct found that contact normalization decreased phosphorylation of 10 proteins on tyrosine residues that were phosphorylated within 1 h of Src kinase activation in transformed cells. Taken together, these results indicate that cadherin dependent contact normalization inhibits Wnt signaling to regulate oncogenic kinase activity and gene expression, particularly PDPN expression, in transformed cells in order to control tumor progression.

Identification and Characterization of Innate Immunity in Actinidia melanandra in Response to Pseudomonas syringae pv. actinidiae

Pseudomonas syringae pv. actinidiae biovar 3 (Psa3) has decimated kiwifruit orchards growing susceptible kiwifruit Actinidia chinensis varieties. Effector loss has occurred recently in Psa3 isolates from resistant kiwifruit germplasm, resulting in strains capable of partially overcoming resistance present in kiwiberry vines (Actinidia arguta, Actinidia polygama, and Actinidia melanandra). Diploid male A. melanandra recognises several effectors, sharing recognition of at least one avirulence effector (HopAW1a) with previously studied tetraploid kiwiberry vines. Sequencing and assembly of the A. melanandra genome enabled the characterisation of the transcriptomic response of this non-host to wild-type and genetic mutants of Psa3. A. melanandra appears to mount a classic effector-triggered immunity (ETI) response to wildtype Psa3 V-13, as expected. Surprisingly, the type III secretion (T3SS) system-lacking Psa3 V-13 ∆hrcC strain did not appear to trigger pattern-triggered immunity (PTI) despite lacking the ability to deliver immunity-suppressing effectors. Contrasting the A. melanandra responses to an effectorless Psa3 V-13 ∆33E strain and to Psa3 V-13 ∆hrcC suggested that PTI triggered by Psa3 V-13 was based on the recognition of the T3SS itself. The characterisation of both ETI and PTI branches of innate immunity responses within A. melanandra further enables breeding for durable resistance in future kiwifruit cultivars.

Loss of TET Activity in the Postnatal Mouse Brain Perturbs Synaptic Gene Expression and Impairs Cognitive Function

Conversion of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) by ten-eleven translocation (TET) family proteins leads to the accumulation of 5hmC in the central nervous system; however, the role of 5hmC in the postnatal brain and how its levels and target genes are regulated by TETs remain elusive. We have generated mice that lack all three Tet genes specifically in postnatal excitatory neurons. These mice exhibit significantly reduced 5hmC levels, altered dendritic spine morphology within brain regions crucial for cognition, and substantially impaired spatial and associative memories. Transcriptome profiling combined with epigenetic mapping reveals that a subset of genes, which display changes in both 5hmC/5mC levels and expression patterns, are involved in synapse-related functions. Our findings provide insight into the role of postnatally accumulated 5hmC in the mouse brain and underscore the impact of 5hmC modification on the expression of genes essential for synapse development and function.

Gene expression differences associated with alcohol use disorder in human brain

Excessive alcohol consumption is a leading cause of preventable death worldwide. To improve understanding of neurobiological mechanisms associated with alcohol use disorder (AUD) in humans, we compared gene expression data from deceased individuals with and without AUD across two addiction-relevant brain regions: the nucleus accumbens (NAc) and dorsolateral prefrontal cortex (DLPFC). Bulk RNA-seq data from NAc and DLPFC (N ≥50 with AUD, ≥46 non-AUD) were analyzed for differential gene expression using modified negative binomial regression adjusting for technical and biological covariates. The region-level results were meta-analyzed with those from an independent dataset (NNAc = 28 AUD, 29 non-AUD; NPFC = 66 AUD, 77 non-AUD). We further tested for heritability enrichment of AUD-related phenotypes, gene co-expression networks, gene ontology enrichment, and drug repurposing. We identified 176 differentially expressed genes (DEGs; 12 in both regions, 78 in NAc only, 86 in DLPFC only) for AUD in our new dataset. After meta-analyzing with published data, we identified 476 AUD DEGs (25 in both regions, 29 in NAc only, 422 in PFC only). Of these DEGs, 17 were significant when looked up in GWAS of problematic alcohol use or drinks per week. Gene co-expression analysis showed both concordant and unique gene networks across brain regions. We also identified 29 and 436 drug compounds that target DEGs from our meta-analysis in NAc and PFC, respectively. This study identified robust AUD-associated DEGs, contributing novel neurobiological insights into AUD and highlighting genes targeted by known drug compounds, generating opportunity for drug repurposing to treat AUD.

ELiAH: the atlas of E3 ligases in human tissues for targeted protein degradation with reduced off-target effect

The development of therapeutic agents has mainly focused on designing small molecules to modulate target proteins or genes which are conventionally druggable. Therefore, targeted protein degradation (TPD) for undruggable cases has emerged as promising pharmaceutical approach. TPD, often referred PROTACs (PROteolysis TArgeting Chimeras), uses a linker to degrade target proteins by hijacking the ubiquitination system. Therefore, unravel the relationship including reversal and co-expression between E3 ligands and other possible target genes in various human tissues is essential to mitigate off-target effects of TPD. Here, we developed the atlas of E3 ligases in human tissues (ELiAH), to prioritize E3 ligase-target gene pairs for TPD. Leveraging over 2900 of RNA-seq profiles consisting of 11 human tissues from the GTEx (genotype-tissue expression) consortium, users of ELiAH can identify tissue-specific genes and E3 ligases (FDR P-value of Mann-Whitney test < .05). ELiAH unravels 933 830 relationships consisting of 614 E3 ligases and 20 924 of expressed genes considering degree of tissue specificity, which are indispensable for ubiquitination based TPD development. In addition, docking properties of those relationships are also modeled using RosettaDock. Therefore, ELiAH presents comprehensive repertoire of E3 ligases for ubiquitination-based TPD drug development avoiding off-target effects. Database URL: https://eliahdb.org.

The neuroendocrine transition in prostate cancer is dynamic and dependent on ASCL1

Lineage plasticity is a hallmark of cancer progression that impacts therapy outcomes, yet the mechanisms mediating this process remain unclear. Here, we introduce a versatile in vivo platform to interrogate neuroendocrine lineage transformation throughout prostate cancer progression. Transplanted mouse prostate organoids with human-relevant driver mutations (Rb1-/-; Trp53-/-; cMyc+ or Pten-/-; Trp53-/-; cMyc+) develop adenocarcinomas, but only those with Rb1 deletion advance to aggressive, ASCL1+ neuroendocrine prostate cancer (NEPC) resistant to androgen receptor signaling inhibitors. Notably, this transition requires an in vivo microenvironment not replicated by conventional organoid culture. Using multiplexed immunofluorescence and spatial transcriptomics, we reveal that ASCL1+ cells arise from KRT8+ luminal cells, progressing into transcriptionally heterogeneous ASCL1+;KRT8- NEPC. Ascl1 loss in established NEPC causes transient regression followed by recurrence, but its deletion before transplantation abrogates lineage plasticity, resulting in castration-sensitive adenocarcinomas. This dynamic model highlights the importance of therapy timing and offers a platform to identify additional lineage plasticity drivers.

Interactions of Polychlorinated Biphenyls and Their Metabolites with the Brain and Liver Transcriptome of Female Mice

Exposure to polychlorinated biphenyls (PCBs) is linked to neurotoxic effects. This study aims to close knowledge gaps regarding the specific modes of action of PCBs in female C57BL/6J mice (>6 weeks) orally exposed for 7 weeks to a human-relevant PCB mixture (MARBLES mix) at 0, 0.1, 1, and 6 mg/kg body weight/day. PCB and hydroxylated PCB (OH-PCBs) levels were quantified in the brain, liver, and serum; RNA sequencing was performed in the striatum, prefrontal cortex, and liver, and metabolomic analyses were performed in the striatum. Profiles of PCBs but not their hydroxylated metabolites were similar in all tissues. In the prefrontal cortex, PCB exposure activated the oxidative phosphorylation respiration pathways, while suppressing the axon guidance pathway. PCB exposure significantly changed the expression of genes associated with neurodevelopmental and neurodegenerative diseases in the striatum, impacting pathways like growth hormone synthesis and dendrite development. PCBs did not affect the striatal metabolome. In contrast to the liver, which showed activation of metabolic processes following PCB exposure and the induction of cytochrome P450 enzymes, the expression of xenobiotic processing genes was not altered by PCB exposure in either brain region. Network analysis revealed complex interactions between individual PCBs (e.g., PCB28 [2,4,4'-trichlorobiphenyl]) and their hydroxylated metabolites and specific differentially expressed genes (DEGs), underscoring the need to characterize the association between specific PCBs and DEGs. These findings enhance the understanding of PCB neurotoxic mechanisms and their potential implications for human health.

Four Core Genotypes mice harbour a 3.2MB X-Y translocation that perturbs Tlr7 dosage

The Four Core Genotypes (FCG) is a mouse model system used to disentangle the function of sex chromosomes and hormones. We report that a copy of a 3.2 MB region of the X chromosome has translocated to the YSry- chromosome and thus increased the expression of X-linked genes including the single-stranded RNA sensor and autoimmune disease mediator Tlr7. This previously-unreported X-Y translocation complicates the interpretation of studies reliant on C57BL/6J FCG mice.

Bioavailable testosterone and androgen receptor activation, but not total testosterone, are associated with muscle mass and strength in females

Testosterone, the major androgen, influences the reproductive and non-reproductive systems in males and females via binding to the androgen receptor (AR). Both circulating endogenous testosterone and muscle AR protein content are positively associated with muscle mass and strength in males, but there is no such evidence in females. Here, we tested whether circulating testosterone levels were associated with muscle mass, function, or the muscle anabolic response to resistance training in pre-menopausal females. Twenty-seven pre-menopausal, untrained females (aged 23.5 ± 4.8 years) underwent a 12-week resistance training programme. Muscle strength, size, power, and plasma and urine androgen hormone levels were measured. Skeletal muscle biopsies were collected before and after the training programme to quantify the effect of resistance training on AR content and nuclear localisation. Primary muscle cell lines were cultured from a subset (n = 6) of the participants' biopsies and treated with testosterone to investigate its effect on myotube diameter, markers of muscle protein synthesis and AR cellular localisation. Physiological levels of total testosterone were not associated with muscle mass or strength at baseline or with the changes in muscle mass and strength that occurred in response to resistance training in our cohort of pre-menopausal females. In contrast, bioavailable testosterone and the proportion of nuclear-localised AR were positively associated with skeletal muscle mass and strength in pre-menopausal females. In vitro, supra-physiological doses of testosterone increased myocyte diameter, but this did not occur via the Akt/mTOR pathway as previously suggested. Instead, we show a marked increase in AR nuclear localisation with testosterone administration in vitro. KEY POINTS: Total circulating testosterone was not related to muscle mass or strength before or after resistance training in pre-menopausal females. Bioavailable testosterone was positively related to exercise-induced muscle hypertrophy in pre-menopausal females. In vivo nuclear localisation of the androgen receptor was positively related to muscle mass in pre-menopausal females at baseline, but not to resistance training-induced hypertrophy. Testosterone treatment induced androgen receptor nuclear translocation but did not induce mTOR signalling in primary skeletal myocytes cultured from pre-menopausal female muscle.

Transcriptomic analyses of host-virus interactions during in vitro infection with wild-type and glycoprotein g-deficient (ΔgG) strains of ILTV in primary and continuous cell cultures

Infectious laryngotracheitis (ILT) remains a significant concern for the poultry industry worldwide due to its impact on animal welfare and its substantial economic consequences. The disease is caused by the alphaherpesvirus, infectious laryngotracheitis virus (ILTV). This study investigated in vitro host-virus interactions of a glycoprotein G (gG) deletion mutant vaccine strain of ILTV (ΔgG ILTV), and its parent wild-type strain (CSW-1 ILTV). Inoculations were performed separately for the two strains of ILTV using both a primary (chicken embryonic kidney, CEK) and a continuous culture (leghorn male hepatoma, LMH) of chicken cells. Transcriptome analysis was performed at 12 hours post infection. Each cell-type displayed distinct effects on host and viral gene transcription, with a greater number of viral and host genes differentially transcribed in CEK cells and LMH cells, respectively. Both cell-types infected with either strain demonstrated enrichment of pathways related to signalling, and gene ontologies (GO) associated with chemotaxis. Infection with either strain upregulated both SOCS proteins and certain proto-oncogenes, which may contribute to prolonged viral persistence by promoting immunosuppression and preventing apoptosis, respectively. Patterns of gene transcription related to cytokines, chemokines, endosomal TLRs, and interferon responses, as well as pathways associated with histone acetylation, transport, and extracellular matrix organization were similar within each cell type, regardless of the viral strain. In CEK cells, GO terms and pathways were downregulated uniquely after CSW-1 ILTV infection, indicating a viral-strain specific effect in this cell-type. Overall, this study highlights that the observed differences in host and ILTV gene transcription in vitro were more strongly influenced by the cell-types used rather than the presence or absence of gG. This underscores the importance of cell-line selection in studying host-virus interactions and interpreting experimental results.

LINE1 modulate human T cell function by regulating protein synthesis during the life span

The molecular mechanisms responsible for the heightened reactivity of quiescent T cells in human early life remain largely elusive. Our previous research identified that quiescent adult naïve CD4+ T cells express LINE1 (long interspersed nuclear elements 1) spliced in previously unknown isoforms, and their down-regulation marks the transition to activation. Here, we unveil that neonatal naïve T cell quiescence is characterized by enhanced energy production and protein synthesis. This phenotype is associated with the absence of LINE1 expression attributed to tonic T cell receptor/mTOR complex 1 (mTORC1) signaling and (polypyrimidine tract-binding protein 1 (PTBP1)-mediated LINE1 splicing suppression. The absence of LINE1 expression primes these cells for rapid execution of the activation program by directly regulating protein synthesis. LINE1 expression progressively increases in childhood and adults, peaking in elderly individuals, and, by decreasing protein synthesis, contributes to immune senescence in aging. Our study proposes LINE1 as a critical player of human T cell function across the human life span.

Best practices for differential accessibility analysis in single-cell epigenomics

Differential accessibility (DA) analysis of single-cell epigenomics data enables the discovery of regulatory programs that establish cell type identity and steer responses to physiological and pathophysiological perturbations. While many statistical methods to identify DA regions have been developed, the principles that determine the performance of these methods remain unclear. As a result, there is no consensus on the most appropriate statistical methods for DA analysis of single-cell epigenomics data. Here, we present a systematic evaluation of statistical methods that have been applied to identify DA regions in single-cell ATAC-seq (scATAC-seq) data. We leverage a compendium of scATAC-seq experiments with matching bulk ATAC-seq or scRNA-seq in order to assess the accuracy, bias, robustness, and scalability of each statistical method. The structure of our experiments also provides the opportunity to define best practices for the analysis of scATAC-seq data beyond DA itself. We leverage this understanding to develop an R package implementing these best practices.

Chronic Activation of Tubulin Tyrosination Improves Heart Function

BACKGROUND

Hypertrophic cardiomyopathy (HCM) is the most common cardiac genetic disorder caused by sarcomeric gene variants and associated with left ventricular hypertrophy and diastolic dysfunction. The role of the microtubule network has recently gained interest with the findings that microtubule detyrosination (dTyr-MT) is markedly elevated in heart failure. Acute reduction of dTyr-MT by inhibition of the detyrosinase (VASH [vasohibin]/SVBP [small VASH-binding protein] complex) or activation of the tyrosinase (TTL [tubulin tyrosine ligase]) markedly improved contractility and reduced stiffness in human failing cardiomyocytes and thus posed a new perspective for HCM treatment. In this study, we tested the impact of chronic tubulin tyrosination in an HCM mouse model (Mybpc3 knock-in), in human HCM cardiomyocytes, and in SVBP-deficient human engineered heart tissues (EHTs).

METHODS

Adeno-associated virus serotype 9-mediated TTL transfer was applied in neonatal wild-type rodents, in 3-week-old knock-in mice, and in HCM human induced pluripotent stem cell-derived cardiomyocytes.

RESULTS

We show (1) TTL for 6 weeks dose dependently reduced dTyr-MT and improved contractility without affecting cytosolic calcium transients in wild-type cardiomyocytes; (2) TTL for 12 weeks reduced the abundance of dTyr-MT in the myocardium, improved diastolic filling, compliance, cardiac output, and stroke volume in knock-in mice; (3) TTL for 10 days normalized cell area in HCM human induced pluripotent stem cell-derived cardiomyocytes; (4) TTL overexpression activated transcription of tubulins and other cytoskeleton components but did not significantly impact the proteome in knock-in mice; (5) SVBP-deficient EHTs exhibited reduced dTyr-MT levels, higher force, and faster relaxation than TTL-deficient and wild-type EHTs. RNA sequencing and mass spectrometry analysis revealed distinct enrichment of cardiomyocyte components and pathways in SVBP-deficient versus TTL-deficient EHTs.

CONCLUSIONS

This study provides the first proof of concept that chronic activation of tubulin tyrosination in HCM mice and in human EHTs improves heart function and holds promise for targeting the nonsarcomeric cytoskeleton in heart disease.

Gene expression profiling of Cacopsylla pyricola (Hemiptera: Psyllidae) infected with Ca. Phytoplasma pyri (Acholeplasmatales: Acholeplasmataceae) reveals candidate effectors and mechanisms of infection

Phytoplasmas can negatively or positively alter vector host fitness. "Candidatus Phytoplasma pyri," is the causal agent of pear decline in commercial pear (Pyrus communis L.; Rosales: Rosaceae) and peach yellow leafroll in peach [Prunus persica (L.); Rosaceae]. This plant pathogen is transmitted by several species of pear psyllids (Cacopsylla spp. Hemiptera: Psyllidae). We sought to explore the relationship between the pear decline phytoplasma and its US vector, Cacopsylla pyricola (Förster), at the molecular genetic level through transcriptomic analysis using RNA-sequencing methodology. We also focused on phytoplasma and insect effectors, which are secreted proteins that can modulate interactions within a pathosystem. In this study, we identified 30 differentially expressed genes, 14 candidate insect effector genes, and 8 Ca. Phytoplasma pyri candidate effectors. Two strains of Ca. Phytoplasma pyri were identified based on immunodominant membrane protein sequence analysis from C. pyricola collected in the Pacific Northwest agricultural region. Here, we present a first genetic look at the pear decline pathosystem and report gene candidates for further exploration of infection mechanisms and potential tools for integrated pest management.

Construction and validation of a prognostic model of angiogenesis-related genes in multiple myeloma

BACKGROUND

Angiogenesis is associated with tumour growth, infiltration, and metastasis. This study aimed to detect the mechanisms of angiogenesis-related genes (ARGs) in multiple myeloma (MM) and to construct a new prognostic model.

METHODS

MM research foundation (MMRF)-CoMMpass cohort, GSE47552, GSE57317, and ARGs were sourced from public databases. Differentially expressed genes (DEGs) in the tumour and control cohorts in GSE47552 were determined through differential expression analysis and were enriched with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. Weighted gene coexpression network analysis (WGCNA) was applied to derive modules linked to the ARG scores and obtain module genes in GSE47552. Differentially expressed ARGs (DE-ARGs) were selected for subsequent analyses by overlapping DEGs and module genes. Furthermore, prognostic genes were selected using univariate Cox and least absolute shrinkage and selection operator (LASSO) regression analyses. Depending on the prognostic genes, a risk model was constructed, and risk scores were determined. Moreover, MM samples from MMRF-CoMMpass were sorted into high- and low-risk teams on account of the median risk score. Additionally, correlations among clinical characteristics, gene set variation analysis (GSVA), gene set enrichment analysis (GSEA), immune analysis, immunotherapy predictions and the mRNA‒miRNA‒lncRNA network were carried out.

RESULTS

A total of 898 DEGs, 211 module genes, 24 DE-ARGs and three prognostic genes (AKAP12, C11orf80 and EMP1) were selected for this study. Enrichment analysis revealed that the DEGs were related to 86 GO terms, such as 'cytoplasmic translation', and 41 KEGG pathways, such as 'small cell lung cancer'. A prognostic gene-based risk model was created in MMRF-CoMMpass and confirmed with the GSE57317 dataset. Moreover, a nomogram was established on the basis of independent prognostic factors that have proven to be good predictors. In addition, the immune cell infiltration results suggested that memory B cells were enriched in the high-risk group and that immature B cells were enriched in the low-risk group. Finally, the mRNA‒miRNA‒lncRNA network demonstrated that hsa-miR-508-5p was tightly associated with EMP1 and AKAP12. RT‒qPCR was used to validate the expression of the genes associated with prognosis.

CONCLUSION

A new prognostic model of MM associated with ARGs was created and validated, providing a new perspective for exploring the connection between ARGs and MM.

A tumor cornification and immune-infiltration-based scheme for anti-PD-1 plus chemotherapy response in advanced squamous cell lung carcinoma

BACKGROUND

Anti-PD-1 immunotherapy plus chemotherapy (combo) exhibits significantly prolonged survival for squamous cell lung cancer (LUSC). An exploration of predictive biomarkers is still needed.

METHODS

High-throughput RNA sequencing (RNA-seq) of 349 LUSC samples from the randomized, multi-center, phase 3 trial ORIENT-12 (ClinicalTrials.gov: NCT03629925) was conducted for biomarker discovery, followed by flow cytometry and multiplex immunohistochemistry (mIHC) in additional clinical cohorts, and in vitro experiments were performed for verification.

RESULTS

A high abundance of activated CD8+ T and CD56bright natural killer (NK) cells benefited patients' outcomes (progression-free survival [PFS]; overall survival [OS]) with combo treatment. Tumor cornification level remarkably affected the infiltration of the two crucial immune cells. Thus, a novel scheme of LUSC immune infiltration and cornification characterization-based classification (LICC) was established for combo efficacy prediction. Patients who received combo treatment achieved significant PFS improvements in LICC1 (hazard ratio [HR] = 0.43, 95% confidence interval [CI]: 0.25-0.75, p = 0.0029) and LICC2 (HR = 0.32, 95% CI: 0.17-0.58, p = 0.0002) subtypes but not in the LICC3 subtype (HR = 0.86, 95% CI: 0.60-1.23, p = 0.4053). Via single-cell RNA-seq analysis, the tumor cornification signal was mainly mapped to SPRR3+ tumor cells, whose relationships with activated CD8+ T or CD56bright NK cells were verified using flow cytometry and mIHC. Our data suggest that SPRR3+ tumor cells might evade immune surveillance via the CD24-SIGLEC10 (M2 macrophage) axis to maintain a suppressive tumor microenvironment.

CONCLUSIONS

Tumor cornification greatly impacts immune infiltration, and the LICC scheme may guide clinical medication of anti-PD-1+chemo treatment in patients with LUSC.

FUNDING

The study was funded by the National Key R&D Program of China, the National Natural Science Foundation of China, Shanghia Multidisplinary Cooperation Building Project for Diagnosis and Treatment of Major Disease, and Innovent Biologics, Inc.

DNA cytosine methylation suppresses meiotic recombination at the sex-determining region

Meiotic recombination between homologous chromosomes is vital for maximizing genetic variation among offspring. However, sex-determining regions are often rearranged and blocked from recombination. It remains unclear whether rearrangements or other mechanisms might be responsible for recombination suppression. Here, we uncover that the deficiency of the DNA cytosine methyltransferase DNMT1 in the green alga Chlamydomonas reinhardtii causes anomalous meiotic recombination at the mating-type locus (MT), generating haploid progeny containing both plus and minus mating-type markers due to crossovers within MT. The deficiency of a histone methyltransferase for H3K9 methylation does not lead to anomalous recombination. These findings suggest that DNA methylation, rather than rearrangements or histone methylation, suppresses meiotic recombination, revealing an unappreciated biological function for DNA methylation in eukaryotes.

Placental Transcriptomic Signatures of Prenatal Phthalate Exposure and Identification of Placenta-Brain Genes Associated with the Effects of Phthalate Exposure on Neurodevelopment in Children

Prenatal exposure to phthalates may affect placental function and fetal development, but the underlying mechanisms are unclear. The aim of our study was to explore the alterations in the placental transcriptome associated with prenatal phthalate exposure and to further analyze whether the placental-brain axis (PBA) genes play a mediating role in the association between prenatal phthalate exposure and children's neurodevelopment. We included 172 participants from the Ma'anshan Birth Cohort and collected data on seven phthalate metabolites in urine during pregnancy, placental tissue RNA-seq, and neurodevelopment of offspring. Bioinformatics analysis revealed that aberrant regulation of the placental transcriptome was associated with prenatal phthalate exposure. Exposure to phthalates during pregnancy was found to be associated with neurodevelopmental delay in children aged 6, 18, and 48 months using the multiple linear regression model. Meanwhile, employing mediation analysis, nine PBA genes were identified that mediate the association between exposure to phthalates during pregnancy and the neurodevelopment of children. Our study will provide a basis for potential mechanisms by which prenatal exposure to phthalates affects placental function and children's neurodevelopment.

Genome-wide mapping of native co-localized G4s and R-loops in living cells

The interplay between G4s and R-loops are emerging in regulating DNA repair, replication, and transcription. A comprehensive picture of native co-localized G4s and R-loops in living cells is currently lacking. Here, we describe the development of HepG4-seq and an optimized HBD-seq methods, which robustly capture native G4s and R-loops, respectively, in living cells. We successfully employed these methods to establish comprehensive maps of native co-localized G4s and R-loops in human HEK293 cells and mouse embryonic stem cells (mESCs). We discovered that co-localized G4s and R-loops are dynamically altered in a cell type-dependent manner and are largely localized at active promoters and enhancers of transcriptional active genes. We further demonstrated the helicase Dhx9 as a direct and major regulator that modulates the formation and resolution of co-localized G4s and R-loops. Depletion of Dhx9 impaired the self-renewal and differentiation capacities of mESCs by altering the transcription of co-localized G4s and R-loops -associated genes. Taken together, our work established that the endogenous co-localized G4s and R-loops are prevalently persisted in the regulatory regions of active genes and are involved in the transcriptional regulation of their linked genes, opening the door for exploring broader roles of co-localized G4s and R-loops in development and disease.

Methotrexate promotes the release of granulocyte-macrophage colony-stimulating factor from rheumatoid arthritis fibroblast-like synoviocytes via autocrine interleukin-1 signaling

BACKGROUND

Activated fibroblast-like synoviocytes (FLS) are drivers of synovitis and structural joint damage in rheumatoid arthritis (RA). Despite the use of disease-modifying drugs, only about 50% of RA patients reach remission in real-world settings. We used an unbiased approach to investigate the effects of standard-of-care methotrexate (MTX) and a Janus kinase inhibitor, tofacitinib (TOFA), on gene expression in RA-FLS, in order to identify untargeted disease mediators.

METHODS

Primary RA-FLS were activated by stimulation with interleukin-1β (IL-1β) or platelet-derived growth factor + IL-1β in the presence or absence of MTX or TOFA, with or without additional inhibitors. Co-cultures of synovial cells were performed in direct and indirect systems. Cells were collected for RNA sequencing or qPCR, and supernatants were analyzed for protein concentrations.

RESULTS

Six thousand three hundred fifty genes were differentially expressed, the majority being upregulated, in MTX-treated activated RA-FLS and 970 genes, the majority being downregulated, in TOFA-treated samples. Pathway analysis showed that MTX had largest effects on 'Molecular mechanisms of cancer' and TOFA on 'Interferon signaling'. Targeted analysis of disease-associated genes revealed that MTX increased the expression of cell cycle-regulating genes but also of pro-inflammatory mediators like IL-1α (IL1A) and granulocyte-macrophage colony-stimulating factor, GM-CSF (CSF2). The MTX-promoted expression of CSF2 in activated RA-FLS peaked at 48 h, could be mediated via either NF-κB or AP-1 transcription factors, and was abrogated by IL-1 inhibitors (IRAK4 inhibitor and anakinra). In a co-culture setting, MTX-treatment of activated RA-FLS induced IL1B expression in macrophages.

CONCLUSIONS

MTX treatment induces secretion of IL-1 from activated RA-FLS which by autocrine signaling augments their release of GM-CSF. This unexpected effect of MTX might contribute to the persistence of synovitis.

Gene regulatory network inference from CRISPR perturbations in primary CD4+ T cells elucidates the genomic basis of immune disease

The effects of genetic variation on complex traits act mainly through changes in gene regulation. Although many genetic variants have been linked to target genes in cis, the trans-regulatory cascade mediating their effects remains largely uncharacterized. Mapping trans-regulators based on natural genetic variation has been challenging due to small effects, but experimental perturbations offer a complementary approach. Using CRISPR, we knocked out 84 genes in primary CD4+ T cells, targeting inborn error of immunity (IEI) disease transcription factors (TFs) and TFs without immune disease association. We developed a novel gene network inference method called linear latent causal Bayes (LLCB) to estimate the network from perturbation data and observed 211 regulatory connections between genes. We characterized programs affected by the TFs, which we associated with immune genome-wide association study (GWAS) genes, finding that JAK-STAT family members are regulated by KMT2A, an epigenetic regulator. These analyses reveal the trans-regulatory cascades linking GWAS genes to signaling pathways.

Inhibition of PRMT5 moderately suppresses prostate cancer growth in vivo but enhances its response to immunotherapy

Protein arginine methylation is a common post-translational modification (PTM) catalyzed by nine protein arginine methyltransferases (PRMTs). As the major symmetric arginine methyltransferase that methylates both histone and non-histone substrates, PRMT5 plays key roles in a number of biological processes critical for development and tumorigenesis. PRMT5 overexpression has been reported in multiple cancer types including prostate cancer (PCa), but the exact biological and mechanistic understanding of PRMT5 in aggressive PCa remains ill-defined. Here, we show that PRMT5 is upregulated in PCa, correlates with worse patient survival, promotes corrupted RNA splicing, and functionally cooperates with an array of pro-tumorigenic pathways to enhance oncogenesis. PRMT5 inhibition via either genetic knockdown or pharmacological inhibition reduces stemness with paralleled differentiation and arrests cell cycle progression without causing appreciable apoptosis. Strikingly, the severity of antitumor effect of PRMT5 inhibition correlates with disease aggressiveness, with AR+ PCa being less affected. Molecular characterization pinpoints MYC, but not (or at least to a lesser degree) AR, as the main partner of PRMT5 to form a positive feedback loop to exacerbate malignancy in both AR+ and AR- PCa cells. Inspired by the surprising finding that PRMT5 negatively correlates with tumor immune infiltration and transcriptionally suppresses an immune-gene program, we further show that although PRMT5 inhibitor (PRMT5i) EPZ015666 or anti-PD-1 immunotherapy alone exhibits limited antitumor effects, combination of PRMT5i with anti-PD-1 displays superior efficacy in inhibiting castration-resistant PCa (CRPC) in vivo. Finally, to expand the potential use of PRMT5i through a synthetic lethality concept, we also perform a global CRISPR/Cas9 knockout screen to unravel that many clinical-grade drugs of known oncogenic pathways can be repurposed to target CRPC when used in combination with PRMT5i at low doses. Collectively, our findings establish a rationale to exploit PRMT5i in combination with immunotherapy or other targeted therapies to treat aggressive PCa.

Synergy between Arbuscular Mycorrhizal Fungi and Rhizosphere Bacterial Communities Increases the Utilization of Insoluble Phosphorus and Potassium in the Soil by Maize

Arbuscular mycorrhizal (AM) fungi can enhance plant uptake of phosphorus (P) and potassium (K), but it is not yet clear whether rhizosphere bacteria can enhance the ability of AM fungi to acquire insoluble P and K from the soil. Here, pot experiments confirmed that AM fungus-promoted insoluble P and K uptake by plants requires rhizosphere bacteria. The changes of rhizosphere bacterial communities associated with AM fungi were explored by 16S rRNA amplicon sequencing and metagenomic sequencing. Five core bacteria genera identified were involved in P and K cycles. Synthetic community (SynCom) inoculation revealed that SynCom increased soil available P and K and its coinoculation with AM fungi increased P and K concentration in the plants. This study revealed that AM fungi interact with rhizosphere bacteria and promote insoluble P and K acquisition, which provided a foundation for the application of AM fungal-bacterial biofertilizers and was beneficial for the sustainable development of agriculture.

Virome Characterization of Native Wild-Rice Plants Discovers a Novel Pathogenic Rice Polerovirus With World-Wide Circulation

Pandemics originating from zoonotic viruses have posed significant threats to human health and agriculture. Recent discoveries have revealed that wild-rice plants also harbour viral pathogens capable of severely impacting rice production, a cornerstone food crop. In this study, we conducted virome analysis on ~1000 wild-rice individual colonies and discovered a novel single-strand positive-sense RNA virus prevalent in these plants. Through comprehensive genomic characterization and comparative sequence analysis, this virus was classified as a new species in the genus Polerovirus, designated Rice less tiller virus (RLTV). Our investigations elucidated that RLTV could be transmitted from wild rice to cultivated rice via a specific insect vector, the aphid Rhopalosiphum padi, causing less tiller disease symptoms in rice plants. We generated an infectious cDNA clone for RLTV and demonstrated systemic infection of rice cultivars and induction of severe disease symptoms following mechanical inoculation or stable genetic transformation. We further illustrated transmission of RLTV from stable transgenic lines to healthy rice plants by the aphid vector, leading to the development of disease symptoms. Notably, our database searches showed that RLTV and another polerovirus isolated from a wild plant species are widely circulating not only in wild rice but also cultivated rice around the world. Our findings provide strong evidence for a wild plant origin for rice viruses and underscore the imminent threat posed by aphid-transmitted rice Polerovirus to rice cultivar.

Critical role of G protein-coupled receptor 40 in B cell response and the pathogenesis of rheumatoid arthritis in mice and patients

Rheumatoid arthritis (RA) is marked by joint damage and inflammation, with B cells playing a key role by generating autoantibodies. This study shows that G protein-coupled receptor 40 (GPR40) deficiency in B cells leads to increased activation, proliferation, antibody production, germinal center formation, and class switch recombination. GPR40 regulates Plcγ2 phosphorylation and intracellular calcium flux downstream of the B cell receptor by binding to the Gαq protein. In GPR40-deficient mice, susceptibility to collagen-induced arthritis was higher. GPR40 agonists showed potential as therapeutic agents, and their reduced expression in patients with RA correlated with disease onset, suggesting GPR40 as a potential therapeutic target and diagnostic marker.

Drought-responsive dynamics of H3K9ac-marked 3D chromatin interactions are integrated by OsbZIP23-associated super-enhancer-like promoter regions in rice

BACKGROUND

In response to drought stress (DS), plants undergo complex processes that entail significant transcriptome reprogramming. However, the intricate relationship between the dynamic alterations in the three-dimensional (3D) genome and the modulation of gene co-expression in drought responses remains a relatively unexplored area.

RESULTS

In this study, we reconstruct high-resolution 3D genome maps based on genomic regions marked by H3K9ac, an active histone modification that dynamically responds to soil water variations in rice. We discover a genome-wide disconnection of 3D genome contact upon DS with over 10,000 chromatin loops lost, which are partially recovered in the subsequent re-watering. Loops integrating promoter-promoter interactions (PPI) contribute to gene expression in addition to basal H3K9ac modifications. Moreover, H3K9ac-marked promoter regions with high affinities in mediating PPIs, termed as super-promoter regions (SPRs), integrate spatially clustered PPIs in a super-enhancer-like manner. Interestingly, the knockout mutation of OsbZIP23, a well-defined DS-responsive transcription factor, leads to the disassociation of over 80% DS-specific PPIs and decreased expression of the corresponding genes under DS. As a case study, we show how OsbZIP23 integrates the PPI cluster formation and the co-expression of four dehydrin genes, RAB16A-D, through targeting the RAB16C SPR in a stress signaling-dependent manner.

CONCLUSIONS

Our high-resolution 3D genome maps unveil the principles and details of dynamic genome folding in response to water supply variations and illustrate OsbZIP23 as an indispensable integrator of the yet unique 3D genome organization that is essential for gene co-expression under DS in rice.

Chromosomal integrons are genetically and functionally isolated units of genomes

Integrons are genetic elements that increase the evolvability of bacteria by capturing new genes and stockpiling them in arrays. Sedentary chromosomal integrons (SCIs) can be massive and highly stabilized structures encoding hundreds of genes, whose function remains generally unknown. SCIs have co-evolved with the host for aeons and are highly intertwined with their physiology from a mechanistic point of view. But, paradoxically, other aspects, like their variable content and location within the genome, suggest a high genetic and functional independence. In this work, we have explored the connection of SCIs to their host genome using as a model the Superintegron (SI), a 179-cassette long SCI in the genome of Vibrio cholerae N16961. We have relocated and deleted the SI using SeqDelTA, a novel method that allows to counteract the strong stabilization conferred by toxin-antitoxin systems within the array. We have characterized in depth the impact in V. cholerae's physiology, measuring fitness, chromosome replication dynamics, persistence, transcriptomics, phenomics, natural competence, virulence and resistance against protist grazing. The deletion of the SI did not produce detectable effects in any condition, proving that-despite millions of years of co-evolution-SCIs are genetically and functionally isolated units of genomes.

Therapeutic effect of bloodletting on bone deterioration induced by hypobaric hypoxia in young rats

OBJECTIVES

High-altitude regions, comprising hypoxic conditions, are associated with different altitude-induced pathologies, including a reduction in bone density. Elucidating the mechanisms underlying bone degradation in such environments and developing targeted interventions and therapeutics is important. Bloodletting therapy has promising clinical applications, but its effects on the skeletal system and bone homeostasis are not well understood. The aim of this study was to investigate the effects of a hypobaric hypoxia environment on specific femoral morphological and structural properties, including bone volume, cortical thickness, and trabecular microarchitecture, in juvenile Sprague-Dawley (SD) rats, and to explore the potential modulating effects of a bloodletting intervention on these parameters.

METHODS

Male SD rats, 6 weeks of age, were subjected to a simulated hypobaric hypoxia environment, replicating a 5000-m altitude, for 12 weeks. For the bloodletting intervention group, rats were subjected to a weekly 500 μL tail vein blood withdrawal. Micro-CT technology, hematoxylin and eosin staining, and tartrate-resistant acid phosphatase staining were employed to comprehensively assess the femoral microstructure, tissue architecture, and cellular morphology. Additionally, immunofluorescence analysis was conducted to quantify the expression of key proteins, and transcriptome analysis was performed to identify differentially expressed genes.

RESULTS

Exposure of rats to hypobaric hypoxia led to a significant reduction in bone mineral content, trabecular bone number, and cortical bone thickness, suggesting a deterioration of bone microstructure. Additionally, the hypoxic environment upregulated the expression of RANKL and HIF-1α, while downregulating RUNX2. Notably, although bloodletting intervention did not significantly reverse these bone structural changes, transcriptome analysis revealed its regulatory influence on the expression of key genes, particularly Mmp2, Fosl2, and URS0000B2A65A, which are implicated in pathways governing the hypoxic response, osteoclast differentiation, and PI3K-Akt signaling.

CONCLUSION

This study highlights the detrimental effect of hypobaric hypoxia on the bone microstructure of juvenile rats and underscores the therapeutic potential of bloodletting to ameliorate this condition. Additionally, our study on the regulatory mechanisms mediating bloodletting's effects on gene expression offers fresh perspectives on bone alterations. It suggests promising avenues for the development of novel preventative measures and targeted therapies to address the challenges posed by related bone disorders.

Systematic perturbations of SETD2, NSD1, NSD2, NSD3, and ASH1L reveal their distinct contributions to H3K36 methylation

BACKGROUND

Methylation of histone 3 lysine 36 (H3K36me) has emerged as an essential epigenetic component for the faithful regulation of gene expression. Despite its importance in development and disease, how the molecular agents collectively shape the H3K36me landscape is unclear.

RESULTS

We use mouse mesenchymal stem cells to perturb the H3K36me methyltransferases (K36MTs) and infer the activities of the five most prominent enzymes: SETD2, NSD1, NSD2, NSD3, and ASH1L. We find that H3K36me2 is the most abundant of the three methylation states and is predominantly deposited at intergenic regions by NSD1, and partly by NSD2. In contrast, H3K36me1/3 are most abundant within exons and are positively correlated with gene expression. We demonstrate that while SETD2 deposits most H3K36me3, it may also deposit H3K36me2 within transcribed genes. Additionally, loss of SETD2 results in an increase of exonic H3K36me1, suggesting other (K36MTs) prime gene bodies with lower methylation states ahead of transcription. While NSD1/2 establish broad intergenic H3K36me2 domains, NSD3 deposits H3K36me2 peaks on active promoters and enhancers. Meanwhile, the activity of ASH1L is restricted to the regulatory elements of developmentally relevant genes, and our analyses implicate PBX2 as a potential recruitment factor.

CONCLUSIONS

Within genes, SETD2 primarily deposits H3K36me3, while the other K36MTs deposit H3K36me1/2 independently of SETD2 activity. For the deposition of H3K36me1/2, we find a hierarchy of K36MT activities where NSD1 > NSD2 > NSD3 > ASH1L. While NSD1 and NSD2 are responsible for most genome-wide propagation of H3K36me2, the activities of NSD3 and ASH1L are confined to active regulatory elements.

Discovery of prognostic lncRNAs in colorectal cancer using spatial transcriptomics

Colorectal cancer (CRC) exhibits significant genetic and epigenetic diversity, evolving into sub-clonal populations with varied metastatic potentials and treatment responses. Predicting metastatic disease in CRC patients remains challenging, underscoring the need for reliable biomarkers. While most research on therapeutic targets and biomarkers has focused on proteins, non-coding RNAs such as long non-coding RNAs (lncRNAs) comprise most of the transcriptome and demonstrate superior tissue- and cancer-specific expression. We utilised spatial transcriptomics to investigate lncRNAs in CRC tumours, offering more precise cell-type-specific expression data compared to bulk RNA sequencing. Our analysis identified 301 lncRNAs linked to malignant CRC regions, which we validated with public data. Further validation using RNA-FISH revealed three lncRNAs (LINC01978, PLAC4, and LINC01303) that are detectable in stage II tumours but not in normal epithelium and are upregulated in metastatic tissues. These lncRNAs hold potential as biomarkers for early risk assessment of metastatic disease.

Type I interferon signaling and peroxisomal dysfunction contribute to enhanced inflammatory cytokine production in Irgm1-deficient macrophages

The human IRGM gene has been linked to inflammatory diseases including sepsis and Crohn's disease. Decreased expression of human IRGM, or of the mouse orthologues Irgm1 and Irgm2, leads to increased production of a number of inflammatory chemokines and cytokines in vivo and/or in cultured macrophages. Prior work has indicated that increased cytokine production is instigated by metabolic alterations and by changes in mitochondrial homeostasis; however, a comprehensive mechanism has not been elucidated. In the studies presented here, RNA deep sequencing and quantitative PCR were used to show that increases in cytokine production, as well as most changes in the transcriptional profile of Irgm1-/- bone marrow-derived macrophages (BMM), are dependent on increased type I IFN production seen in those cells. Metabolic alterations that drive increased cytokines in Irgm1-/- BMM - specifically increases in glycolysis and increased accumulation of acyl-carnitines - were unaffected by quenching type I IFN signaling. Dysregulation of peroxisomal homeostasis was identified as a novel upstream pathway that governs type I IFN production and inflammatory cytokine production. Collectively, these results enhance our understanding of the complex biochemical changes that are triggered by lack of Irgm1 and contribute to inflammatory disease seen with Irgm1-deficiency.

Stress granules are not present in Kras mutant cancers and do not control tumor growth

Stress granules (SG) are membraneless ribonucleoprotein-based cytoplasmic organelles that assemble in response to stress. Their formation is often associated with an almost global suppression of translation, and the aberrant assembly or disassembly of these granules has pathological implications in neurodegeneration and cancer. In cancer, and particularly in the presence of oncogenic KRAS mutations, in vivo studies concluded that SG increase the resistance of cancer cells to stress. Hence, SG have recently been considered a promising target for therapy. Here, starting from our observations that genes coding for SG proteins are stimulated during development of pancreatic ductal adenocarcinoma, we analyze the formation of SG during tumorigenesis. We resort to in vitro, in vivo and in silico approaches, using mouse models, human samples and human data. Our analyses do not support that SG are formed during tumorigenesis of KRAS-driven cancers, at least that their presence is not universal, leading us to propose that caution is required before considering SG as therapeutic targets.

Single cell atlas reveals multilayered metabolic heterogeneity across tumour types

BACKGROUND

Metabolic reprogramming plays a pivotal role in cancer progression, contributing to substantial intratumour heterogeneity and influencing tumour behaviour. However, a systematic characterization of metabolic heterogeneity across multiple cancer types at the single-cell level remains limited.

METHODS

We integrated 296 tumour and normal samples spanning six common cancer types to construct a single-cell compendium of metabolic gene expression profiles and identify cell type-specific metabolic properties and reprogramming patterns. A computational approach based on non-negative matrix factorization (NMF) was utilised to identify metabolic meta-programs (MMPs) showing intratumour heterogeneity. In-vitro cell experiments were conducted to confirm the associations between MMPs and chemotherapy resistance, as well as the function of key metabolic regulators. Survival analyses were performed to assess clinical relevance of cellular metabolic properties.

FINDINGS

Our analysis revealed shared glycolysis upregulation and divergent regulation of citric acid cycle across different cell types. In malignant cells, we identified a colorectal cancer-specific MMP associated with resistance to the cuproptosis inducer elesclomol, validated through in-vitro cell experiments. Furthermore, our findings enabled the stratification of patients into distinct prognostic subtypes based on metabolic properties of specific cell types, such as myeloid cells.

INTERPRETATION

This study presents a nuanced understanding of multilayered metabolic heterogeneity, offering valuable insights into potential personalized therapies targeting tumour metabolism.

FUNDING

National Key Research and Development Program of China (2021YFA1300601). National Natural Science Foundation of China (key grants 82030081 and 81874235). The Shenzhen High-level Hospital Construction Fund and Shenzhen Basic Research Key Project (JCYJ20220818102811024). The Lam Chung Nin Foundation for Systems Biomedicine.

kallisto, bustools and kb-python for quantifying bulk, single-cell and single-nucleus RNA-seq

The term 'RNA-seq' refers to a collection of assays based on sequencing experiments that involve quantifying RNA species from bulk tissue, single cells or single nuclei. The kallisto, bustools and kb-python programs are free, open-source software tools for performing this analysis that together can produce gene expression quantification from raw sequencing reads. The quantifications can be individualized for multiple cells, multiple samples or both. Additionally, these tools allow gene expression values to be classified as originating from nascent RNA species or mature RNA species, making this workflow amenable to both cell-based and nucleus-based assays. This protocol describes in detail how to use kallisto and bustools in conjunction with a wrapper, kb-python, to preprocess RNA-seq data. Execution of this protocol requires basic familiarity with a command line environment. With this protocol, quantification of a moderately sized RNA-seq dataset can be completed within minutes.