The R package Rsubread is easier, faster, cheaper and better for alignment and quantification of RNA sequencing reads

Drought-induced molecular changes in crown of various barley phytohormone mutants

One of the main signal transduction pathways that modulate plant growth and stress responses, including drought, is the action of phytohormones. Recent advances in omics approaches have facilitated the exploration of plant genomes. However, the molecular mechanisms underlying the response in the crown of barley, which plays an essential role in plant performance under stress conditions and regeneration after stress treatment, remain largely unclear. The objective of the present study was the elucidation of drought-induced molecular reactions in the crowns of different barley phytohormone mutants. We verified the hypothesis that defects of gibberellins, brassinosteroids, and strigolactones action affect the transcriptomic, proteomic, and hormonal response of barley crown to the transitory drought influencing plant development under stress. Moreover, we assumed that due to the strong connection between strigolactones and branching the hvdwarf14.d mutant, with dysfunctional receptor of strigolactones, manifests the most abundant alternations in crowns and phenotype under drought. Finally, we expected to identify components underlying the core response to drought which are independent of the genetic background. Large-scale analyses were conducted using gibberellins-biosynthesis, brassinosteroids-signaling, and strigolactones-signaling mutants, as well as reference genotypes. Detailed phenotypic evaluation was also conducted. The obtained results clearly demonstrated that hormonal disorders caused by mutations in the HvGA20ox2, HvBRI1, and HvD14 genes affected the multifaceted reaction of crowns to drought, although the expression of these genes was not induced by stress. The study further detected not only genes and proteins that were involved in the drought response and reacted specifically in mutants compared to the reaction of reference genotypes and vice versa, but also the candidates that may underlie the genotype-universal stress response. Furthermore, candidate genes involved in phytohormonal interactions during the drought response were identified. We also found that the interplay between hormones, especially gibberellins and auxins, as well as strigolactones and cytokinins may be associated with the regulation of branching in crowns exposed to drought. Overall, the present study provides novel insights into the molecular drought-induced responses that occur in barley crowns.

The importance of DNA sequence for nucleosome positioning in transcriptional regulation

Nucleosome positioning is a key factor for transcriptional regulation. Nucleosomes regulate the dynamic accessibility of chromatin and interact with the transcription machinery at every stage. Influences to steer nucleosome positioning are diverse, and the according importance of the DNA sequence in contrast to active chromatin remodeling has been the subject of long discussion. In this study, we evaluate the functional role of DNA sequence for all major elements along the process of transcription. We developed a random forest classifier based on local DNA structure that assesses the sequence-intrinsic support for nucleosome positioning. On this basis, we created a simple data resource that we applied genome-wide to the human genome. In our comprehensive analysis, we found a special role of DNA in mediating the competition of nucleosomes with cis-regulatory elements, in enabling steady transcription, for positioning of stable nucleosomes in exons, and for repelling nucleosomes during transcription termination. In contrast, we relate these findings to concurrent processes that generate strongly positioned nucleosomes in vivo that are not mediated by sequence, such as energy-dependent remodeling of chromatin.

Genome-wide DNA methylation and transcriptomic analysis of liver tissues subjected to early ischemia/reperfusion injury upon human liver transplantation

INTRODUCTION AND OBJECTIVES

Epigenetic changes represent a mechanism connecting external stresses with long-term modifications of gene expression programs. In solid organ transplantation, ischemia-reperfusion injury (IRI) appears to induce epigenomic changes in the graft, although the currently available data are extremely limited. The present study aimed to characterize variations in DNA methylation and their effects on the transcriptome in liver transplantation from brain-dead donors.

PATIENTS AND METHODS

12 liver grafts were evaluated through serial biopsies at different timings in the procurement-transplantation process: T0 (warm procurement, in donor), T1 (bench surgery), and T2 (after reperfusion, in recipient). DNA methylation (DNAm) and transcriptome profiles of biopsies were analyzed using microarrays and RNAseq.

RESULTS

Significant variations in DNAm were identified, particularly between T2 and T0. Functional enrichment of the best 1000 ranked differentially methylated promoters demonstrated that 387 hypermethylated and 613 hypomethylated promoters were involved in spliceosomal assembly and response to biotic stimuli, and inflammatory immune responses, respectively. At the transcriptome level, T2 vs. T0 showed an upregulation of 337 and downregulation of 61 genes, collectively involved in TNF-α, NFKB, and interleukin signaling. Cell enrichment analysis individuates macrophages, monocytes, and neutrophils as the most significant tissue-cell type in the response.

CONCLUSIONS

In the process of liver graft procurement-transplantation, IRI induces significant epigenetic changes that primarily act on the signaling pathways of inflammatory responses dependent on TNF-α, NFKB, and interleukins. Our DNAm datasets are the early IRI methylome literature and will serve as a launch point for studying the impact of epigenetic modification in IRI.

Impact of freshwater rearing history on Atlantic salmon gill response to viral stimulation post seawater transfer

Land-based recirculating aquaculture systems (RAS) have risen in prevalence in recent years for Atlantic salmon production, enabling intensive production which allows increased growth and environmental control, but also having the potential for reducing water use and eutrophication. The Atlantic salmon has an anadromous life history with juvenile stages in freshwater (FW) and on-growing in seawater (SW), enabled by a transformational process known as smoltification. The timing of smoltification and transfer of smolts from FW to SW is critical under commercial production with high mortalities during this period. The impact of FW rearing system on immune function following seawater transfer (SWT) is not well understood. In this study parr were raised in either RAS or a traditional open-LOCH system until smolting and then transferred to a common marine environment. Two-weeks post-SWT fish were immune stimulated with a viral mimic (poly I:C) for 24 h to assess the ability to mount an antiviral immune response, assessed by whole transcriptome analysis of gill tissue, an important immune organ in fish. We show that unstimulated smolts reared in the LOCH had higher immune gene expression than those reared in RAS as determined by functional analysis. However, following stimulation, smolts reared in the RAS mounted a greater magnitude of response with a suite of immune genes displaying higher fold induction of transcription compared to LOCH reared smolts. We suggest RAS smolts have a lower steady state immune-associated transcriptome likely due to an unvarying environment, in terms of environmental factors and lack of exposure to pathogens, which shows a compensatory mechanism following stimulation allowing immune 'catch-up' with those reared in the LOCH. Alternatively, the RAS fish are experiencing an excessive response to the immune stimulation.

Alternative splicing as a driver of natural variation in abscisic acid response

Abscisic acid (ABA) is a crucial player in plant responses to the environment. It accumulates under stress, activating downstream signaling to implement molecular responses that restore homeostasis. Natural variance in ABA sensitivity remains barely understood, and the ABA pathway has been mainly studied at the transcriptional level, despite evidence that posttranscriptional regulation, namely, via alternative splicing, contributes to plant stress tolerance. Here, we identified the Arabidopsis accession Kn-0 as less sensitive to ABA than the reference Col-0, as shown by reduced effects of the hormone on seedling establishment, root branching, and stomatal closure, as well as by decreased induction of ABA marker genes. An in-depth comparative transcriptome analysis of the ABA response in the two variants revealed lower expression changes and fewer genes affected for the least ABA-sensitive ecotype. Notably, Kn-0 exhibited reduced levels of the ABA-signaling SnRK2 protein kinases and lower basal expression of ABA-reactivation genes, consistent with our finding that Kn-0 contains less endogenous ABA than Col-0. ABA also markedly affected alternative splicing, primarily intron retention, with Kn-0 being less responsive regarding both the number and magnitude of alternative splicing events, particularly exon skipping. We find that alternative splicing introduces a more ecotype-specific layer of ABA regulation and identify ABA-responsive splicing changes in key ABA pathway regulators that provide a functional and mechanistic link to the differential sensitivity of the two ecotypes. Our results offer new insight into the natural variation of ABA responses and corroborate a key role for alternative splicing in implementing ABA-mediated stress responses.

Nkd1 functions downstream of Axin2 to attenuate Wnt signaling

Wnt signaling is a crucial developmental pathway involved in early development as well as stem-cell maintenance in adults and its misregulation leads to numerous diseases. Thus, understanding the regulation of this pathway becomes vitally important. Axin2 and Nkd1 are widely utilized negative feedback regulators in Wnt signaling where Axin2 functions to destabilize cytoplasmic β-catenin, and Nkd1 functions to inhibit the nuclear localization of β-catenin. Here, we set out to further understand how Axin2 and Nkd1 regulate Wnt signaling by creating axin2gh1/gh1, nkd1gh2/gh2 single mutants and axin2gh1/gh1;nkd1gh2/gh2 double mutant zebrafish using sgRNA/Cas9. All three Wnt regulator mutants were viable and had impaired heart looping, neuromast migration defects, and behavior abnormalities in common, but there were no signs of synergy in the axin2gh1/gh1;nkd1gh2/gh2 double mutants. Further, Wnt target gene expression by qRT-PCR and RNA-seq, and protein expression by mass spectrometry demonstrated that the double axin2gh1/gh1;nkd1gh2/gh2 mutant resembled the nkd1gh2/gh2 phenotype demonstrating that Nkd1 functions downstream of Axin2. In support of this, the data further demonstrates that Axin2 uniquely alters the properties of β-catenin-dependent transcription having novel readouts of Wnt activity compared with nkd1gh2/gh2 or the axin2gh1/gh1;nkd1gh2/gh2 double mutant. We also investigated the sensitivity of the Wnt regulator mutants to exacerbated Wnt signaling, where the single mutants displayed characteristic heightened Wnt sensitivity, resulting in an eyeless phenotype. Surprisingly, this phenotype was rescued in the double mutant, where we speculate that cross-talk between Wnt/β-catenin and Wnt/Planar Cell Polarity pathways could lead to altered Wnt signaling in some scenarios. Collectively, the data emphasizes both the commonality and the complexity in the feedback regulation of Wnt signaling.

Mimicking bacterial infection in male mice changes sperm small RNA profiles and multigenerationally alters offspring behavior and physiology

Paternal pre-conceptual exposures, including stress, diet, substance abuse, parasite infection, and viral immune activation via Poly I:C, have been reported to influence the brains and behavior of offspring through sperm epigenetic changes. However, the effects of paternal (F0) pre-conceptual exposure to bacterial-induced immune activation on the behavior and physiology of F1 and F2 generations remain unexplored. We examined this using C57BL/6J mice. Eight-week-old males (F0) received a single intraperitoneal injection of the bacterial mimetic lipopolysaccharide (LPS: 5 mg/kg) or 0.9 % saline (vehicle control) before mating with naïve females at four weeks post-injection. Comprehensive behavioral assessments were conducted to investigate anxiety, social behaviors, depressive-like behaviors and cognition in both the F1 and F2 generations within the age range of 8 to 14 weeks. Results demonstrated that only female offspring of LPS-exposed fathers exhibited reduced anxiety levels in the light/dark box, large open field, and novelty-suppressed feeding test. These F1 female offspring also exhibited heightened sociability in the 3-chambered social interaction test and a reduced preference for saccharin in the saccharin preference test. Additionally, the F1 male offspring of LPS-challenged males demonstrated an increased total distance traveled in the light/dark box and a longer distance covered in the light zone. They also exhibited diminished preference for social novelty in the 3-chambered social interaction test and an elevated novel arm preference index in the Y-maze. In the F2 generation, male descendants of LPS-treated fathers showed reduced latency to feed in the novelty-suppressed feeding test. Additionally, the F2 generation of LPS-challenged fathers, but not the F1 generation, displayed enhanced immune response in both sexes after an acute LPS immune challenge (5 mg/kg). Analysis of sperm small noncoding RNA profiles from LPS-treated F0 mice revealed significant changes at 4 weeks after administration of LPS. These changes included three microRNAs, eight PIWI-interacting RNAs, and two transfer RNAs, exhibiting significant upregulation (mmu-miR-146a-5p, mmu-piR-27082 and mmu-piR-29102) or downregulation (mmu-miR-5110, mmu-miR-467e-3p, mmu-piR-22583, mmu-piR-23548, mmu-piR-36341, mmu-piR-50293, mmu-piR-16583, mmu-piR-36507, Mus_musculus_tRNA-Ile-AAT-2-1 and Mus_musculus_tRNA-Tyr-GTA-1-1). Additionally, we detected 52 upregulated small noncoding RNAs (including 9 miRNAs, 41 piRNAs, and 2 tRNAs) and 7 downregulated small noncoding RNAs (3 miRNAs, 3 piRNAs, and 1 tRNA) in the sperm of F1 offspring from LPS-treated males. These findings provide compelling evidence for the involvement of epigenetic mechanisms in the modulation of brain function and immunity, and associated behavioral and immunological traits, across generations, in response to bacterial infection.

Arabidopsis mRNA export factor MOS11: molecular interactions and role in abiotic stress responses

Transcription and export (TREX) is a multi-subunit complex that links synthesis, processing and export of mRNAs. It interacts with the RNA helicase UAP56 and export factors such as MOS11 and ALYs to facilitate nucleocytosolic transport of mRNAs. Plant MOS11 is a conserved, but sparsely researched RNA-binding export factor, related to yeast Tho1 and mammalian CIP29/SARNP. Using biochemical approaches, the domains of Arabidopsis thaliana MOS11 required for interaction with UAP56 and RNA-binding were identified. Further analyses revealed marked genetic interactions between MOS11 and ALY genes. Cell fractionation in combination with transcript profiling demonstrated that MOS11 is required for export of a subset of mRNAs that are shorter and more GC-rich than MOS11-independent transcripts. The central α-helical domain of MOS11 proved essential for physical interaction with UAP56 and for RNA-binding. MOS11 is involved in the nucleocytosolic transport of mRNAs that are upregulated under stress conditions and accordingly mos11 mutant plants turned out to be sensitive to elevated NaCl concentrations and heat stress. Collectively, our analyses identify functional interaction domains of MOS11. In addition, the results establish that mRNA export is critically involved in the plant response to stress conditions and that MOS11 plays a prominent role at this.

Oxylipins and metabolites from pyroptotic cells act as promoters of tissue repair

Pyroptosis is a lytic cell death mode that helps limit the spread of infections and is also linked to pathology in sterile inflammatory diseases and autoimmune diseases1-4. During pyroptosis, inflammasome activation and the engagement of caspase-1 lead to cell death, along with the maturation and secretion of the inflammatory cytokine interleukin-1β (IL-1β). The dominant effect of IL-1β in promoting tissue inflammation has clouded the potential influence of other factors released from pyroptotic cells. Here, using a system in which macrophages are induced to undergo pyroptosis without IL-1β or IL-1α release (denoted Pyro-1), we identify unexpected beneficial effects of the Pyro-1 secretome. First, we noted that the Pyro-1 supernatants upregulated gene signatures linked to migration, cellular proliferation and wound healing. Consistent with this gene signature, Pyro-1 supernatants boosted migration of primary fibroblasts and macrophages, and promoted faster wound closure in vitro and improved tissue repair in vivo. In mechanistic studies, lipidomics and metabolomics of the Pyro-1 supernatants identified the presence of both oxylipins and metabolites, linking them to pro-wound-healing effects. Focusing specifically on the oxylipin prostaglandin E2 (PGE2), we find that its synthesis is induced de novo during pyroptosis, downstream of caspase-1 activation and cyclooxygenase-2 activity; further, PGE2 synthesis occurs late in pyroptosis, with its release dependent on gasdermin D pores opened during pyroptosis. As for the pyroptotic metabolites, they link to immune cell infiltration into the wounds, and polarization to CD301+ macrophages. Collectively, these data advance the concept that the pyroptotic secretome possesses oxylipins and metabolites with tissue repair properties that may be harnessed therapeutically.

Neuron-specific chromatin disruption at CpG islands and aging-related regions in Kabuki syndrome mice

Many Mendelian developmental disorders caused by coding variants in epigenetic regulators have now been discovered. Epigenetic regulators are broadly expressed, and each of these disorders typically shows phenotypic manifestations from many different organ systems. An open question is whether the chromatin disruption-the root of the pathogenesis-is similar in the different disease-relevant cell types. This is possible in principle, because all these cell types are subject to effects from the same causative gene, which has the same kind of function (e.g., methylates histones) and is disrupted by the same germline variant. We focus on mouse models for Kabuki syndrome types 1 and 2 and find that the chromatin accessibility changes in neurons are mostly distinct from changes in B or T cells. This is not because the neuronal accessibility changes occur at regulatory elements that are only active in neurons. Neurons, but not B or T cells, show preferential chromatin disruption at CpG islands and at regulatory elements linked to aging. A sensitive analysis reveals that regulatory elements disrupted in B/T cells do show chromatin accessibility changes in neurons, but these are very subtle and of uncertain functional significance. Finally, we are able to identify a small set of regulatory elements disrupted in all three cell types. Our findings reveal the cellular-context-specific effect of variants in epigenetic regulators and suggest that blood-derived episignatures, although useful diagnostically, may not be well suited for understanding the mechanistic basis of neurodevelopment in Mendelian disorders of the epigenetic machinery.

Far from Their Origins: A Transcriptomic Investigation on How 2,4-Di-tert-butyl-6-(5-chloro-2H-benzotriazol-2-yl) Phenol Affects Rainbow Trout Alevins

Benzotriazole ultraviolet stabilizers (BUVSs) are a group of widely used chemicals added to a variety of consumer (e.g., plastics) and industrial (e.g., metal coating) goods. Although detected globally as an environmentally persistent pollutant, BUVSs have received relatively little toxicological attention and only recently have been acknowledged to affect development and the endocrine system in vivo. In our previous study, altered behavior, indicative of potential neurotoxicity, was observed among rainbow trout alevins (day 14 posthatching) that were microinjected as embryos with a single environmentally relevant dose of 2,4-di-tert-butyl-6-(5-chloro-2H-benzotriazol-2-yl) phenol (UV-327). In the present follow-up study, we performed whole-transcriptome profiling (RNA sequencing) of newly hatched alevins from the same batch. The primary aim was to identify biomarkers related to behavior and neurology. Dose-specifically, 1 to 176 differentially expressed genes (DEGs) were identified. In the group presenting altered behavior (273.4 ng g-1), 176 DEGs were identified, yet only a fraction was related to neurological functions, including water, calcium, and potassium homeostasis; acetylcholine transmission and signaling; as well insulin and energy metabolism. The second objective was to estimate the transcriptomic point of departure (tPOD) and assess if point estimate(s) are protective of altered behavior. A tPOD was established at 35 to 94 ng UV-327 g-1 egg, making this tPOD protective of behavioral alterations. Holistically, these transcriptomic alterations provide a foundation for future research on how BUVSs can influence rainbow trout alevin development, while providing support to the hypothesis that UV-327 can influence neurogenesis and subsequent behavioral endpoints. The exact structural and functional changes caused by embryonic exposure to UV-327 remain enigmatic and will require extensive investigation before being deciphered and understood toxicologically. Environ Toxicol Chem 2024;00:1-12. © 2024 The Author(s). Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

H3.1/3.2 regulate the initial progression of the gene expression program

In mice, transcription from the zygotic genome is initiated at the mid-one-cell stage, and occurs promiscuously in many areas of the genome, including intergenic regions. Regulated transcription from selected genes is established during the two-cell stage. This dramatic change in the gene expression pattern marks the initiation of the gene expression program and is essential for early development. We investigated the involvement of the histone variants H3.1/3.2 in the regulation of changes in gene expression pattern during the two-cell stage. Immunocytochemistry analysis showed low nuclear deposition of H3.1/3.2 in the one-cell stage, followed by a rapid increase in the late two-cell stage. Where chromatin structure is normally closed between the one- and two-cell stages, it remained open until the late two-cell stage when H3.1/3.2 were knocked down by small interfering RNA. Hi-C analysis showed that the formation of the topologically associating domain was disrupted in H3.1/3.2 knockdown (KD) embryos. Promiscuous transcription was also maintained in the late two-cell stage in H3.1/3.2 KD embryos. These results demonstrate that H3.1/3.2 are involved in the initial process of the gene expression program after fertilization, through the formation of a closed chromatin structure to execute regulated gene expression during the two-cell stage.

SMAD4 depletion contributes to endocrine resistance by integrating ER and ERBB signaling in HR + HER2- breast cancer

Endocrine resistance poses a significant clinical challenge for patients with hormone receptor-positive and human epithelial growth factor receptor 2-negative (HR + HER2-) breast cancer. Dysregulation of estrogen receptor (ER) and ERBB signaling pathways is implicated in resistance development; however, the integration of these pathways remains unclear. While SMAD4 is known to play diverse roles in tumorigenesis, its involvement in endocrine resistance is poorly understood. Here, we investigate the role of SMAD4 in acquired endocrine resistance in HR + HER2- breast cancer. Genome-wide CRISPR screening identifies SMAD4 as a regulator of 4-hydroxytamoxifen (OHT) sensitivity in T47D cells. Clinical data analysis reveals downregulated SMAD4 expression in breast cancer tissues, correlating with poor prognosis. Following endocrine therapy, SMAD4 expression is further suppressed. Functional studies demonstrate that SMAD4 depletion induces endocrine resistance in vitro and in vivo by enhancing ER and ERBB signaling. Concomitant inhibition of ER and ERBB signaling leads to aberrant autophagy activation. Simultaneous inhibition of ER, ERBB, and autophagy pathways synergistically impacts SMAD4-depleted cells. Our findings unveil a mechanism whereby endocrine therapy-induced SMAD4 downregulation drives acquired resistance by integrating ER and ERBB signaling and suggest a rational treatment strategy for endocrine-resistant HR + HER2- breast cancer patients.

Amino acid influx via LAT1 regulates iron demand and sensitivity to PPMX-T003 of aggressive natural killer cell leukemia

Aggressive natural killer cell leukemia (ANKL) is a rare hematological malignancy with a fulminant clinical course. Our previous study revealed that ANKL cells proliferate predominantly in the liver sinusoids and strongly depend on transferrin supplementation. In addition, we demonstrated that liver-resident ANKL cells are sensitive to PPMX-T003, an anti-human transferrin receptor 1 inhibitory antibody, whereas spleen-resident ANKL cells are resistant to transferrin receptor 1 inhibition. However, the microenvironmental factors that regulate the iron dependency of ANKL cells remain unclear. In this study, we first revealed that the anti-neoplastic effect of PPMX-T003 was characterized by DNA double-strand breaks in a DNA replication-dependent manner, similar to conventional cytotoxic agents. We also found that the influx of extracellular amino acids via LAT1 stimulated sensitivity to PPMX-T003. Taken together, we discovered that the amount of extracellular amino acid influx through LAT1 was the key environmental factor determining the iron dependency of ANKL cells via adjustment of their mTOR/Myc activity, which provides a good explanation for the different sensitivity to PPMX-T003 between liver- and spleen-resident ANKL cells, as the liver sinusoid contains abundant amino acids absorbed from the gut.

HumanIslets: An integrated platform for human islet data access and analysis

Comprehensive molecular and cellular phenotyping of human islets can enable deep mechanistic insights for diabetes research. We established the Human Islet Data Analysis and Sharing (HI-DAS) consortium to advance goals in accessibility, usability, and integration of data from human islets isolated from donors with and without diabetes at the Alberta Diabetes Institute (ADI) IsletCore. Here we introduce HumanIslets.com , an open resource for the research community. This platform, which presently includes data on 547 human islet donors, allows users to access linked datasets describing molecular profiles, islet function and donor phenotypes, and to perform various statistical and functional analyses at the donor, islet and single-cell levels. As an example of the analytic capacity of this resource we show a dissociation between cell culture effects on transcript and protein expression, and an approach to correct for exocrine contamination found in hand-picked islets. Finally, we provide an example workflow and visualization that highlights links between type 2 diabetes status, SERCA3b Ca 2+ -ATPase levels at the transcript and protein level, insulin secretion and islet cell phenotypes. HumanIslets.com provides a growing and adaptable set of resources and tools to support the metabolism and diabetes research community.

ISAba1 mediated intrinsic chromosomal oxacillinase amplification confers carbapenem resistance in Acinetobacter baumannii

Tandem amplification of carbapenemase genes increases gene copy number and enhances carbapenem resistance. These amplifications are often heterogeneous, transient, and located on plasmids, which also contribute to heteroresistance. Amplification of encoding genes is especially important for enzymes with low hydrolysis activity, which are often overlooked. Here, we reported an intrinsic oxacillinase oxaAb amplification flanked by ISAba1. The amplification is in the chromosome and contains up to twenty-five repeats. We provided genomic, transcriptomic, and proteomic evidence that the amplification resulted in oxacillinase overproduction. Notably, no point mutations of oxaAb were found during the amplification process. Strains of A. baumannii with intrinsic amplified or external transformed ISAba1-oxaAb exhibited higher meropenem hydrolysis activity. Furthermore, the number of repeats in the amplification decreased gradually over a period of 21 days cultured with carbapenem withdrawal. However, upon re-exposure to meropenem, the ISAba1 flanked oxaAb responded rapidly, with repeat numbers reaching or exceeding pre-carbapenem withdrawal levels within 24 hours. Taken together, these findings suggest that ISAba1-mediated gene amplification and overproduction of intrinsic low-activity oxacillinase oxaAb resulted in carbapenem resistance.

Ribosome profiling reveals the role of yeast RNA-binding proteins Cth1 and Cth2 in translational regulation

Iron serves as a cofactor for enzymes involved in several steps of protein translation, but the control of translation during iron limitation is not understood at the molecular level. Here, we report a genome-wide analysis of protein translation in response to iron deficiency in yeast using ribosome profiling. We show that iron depletion affects global protein synthesis and leads to translational repression of multiple genes involved in iron-related processes. Furthermore, we demonstrate that the RNA-binding proteins Cth1 and Cth2 play a central role in this translational regulation by repressing the activity of the iron-dependent Rli1 ribosome recycling factor and inhibiting mitochondrial translation and heme biosynthesis. Additionally, we found that iron deficiency represses MRS3 mRNA translation through increased expression of antisense long non-coding RNA. Together, our results reveal complex gene expression and protein synthesis remodeling in response to low iron, demonstrating how this important metal affects protein translation at multiple levels.

Adult microglial TGFβ1 is required for microglia homeostasis via an autocrine mechanism to maintain cognitive function in mice

While TGF-β signaling is essential for microglial function, the cellular source of TGF-β1 ligand and its spatial regulation remains unclear in the adult CNS. Our data supports that microglia but not astrocytes or neurons are the primary producers of TGF-β1 ligands needed for microglial homeostasis. Microglia-Tgfb1 KO leads to the activation of microglia featuring a dyshomeostatic transcriptome that resembles disease-associated, injury-associated, and aged microglia, suggesting microglial self-produced TGF-β1 ligands are important in the adult CNS. Astrocytes in MG-Tgfb1 inducible (i)KO mice show a transcriptome profile that is closely aligned with an LPS-associated astrocyte profile. Additionally, using sparse mosaic single-cell microglia KO of TGF-β1 ligand we established an autocrine mechanism for signaling. Here we show that MG-Tgfb1 iKO mice present cognitive deficits, supporting that precise spatial regulation of TGF-β1 ligand derived from microglia is required for the maintenance of brain homeostasis and normal cognitive function in the adult brain.

Improved resilience and proteostasis mediate longevity upon DAF-2 degradation in old age

Little is known about the possibility of reversing age-related biological changes when they have already occurred. To explore this, we have characterized the effects of reducing insulin/IGF-1 signaling (IIS) during old age. Reduction of IIS throughout life slows age-related decline in diverse species, most strikingly in the nematode Caenorhabditis elegans. Here we show that even at advanced ages, auxin-induced degradation of DAF-2 in single tissues, including neurons and the intestine, is still able to markedly increase C. elegans lifespan. We describe how reversibility varies among senescent changes. While senescent pathologies that develop in mid-life were not reversed, there was a rejuvenation of the proteostasis network, manifesting as a restoration of the capacity to eliminate otherwise intractable protein aggregates that accumulate with age. Moreover, resistance to several stressors was restored. These results support several new conclusions. (1) Loss of resilience is not solely a consequence of pathologies that develop in earlier life. (2) Restoration of proteostasis and resilience by inhibiting IIS is a plausible cause of the increase in lifespan. And (3), most interestingly, some aspects of the age-related transition from resilience to frailty can be reversed to a certain extent. This raises the possibility that the effect of IIS and related pathways on resilience and frailty during aging in higher animals might possess some degree of reversibility.

Core Transcriptome of Hydrogen Producing Marine Vibrios Reveals Contribution of Glycolysis in Their Efficient Hydrogen Production

Pyruvate (Pyr) is the end product of the glycolysis pathway. Pyr is also renewable and is further metabolized to produce formate, which is the precursor of H2, via pyruvate formate lyase (PFL) under anaerobic conditions. The formate is excluded and re-imported via the formate channel and is then converted to H2 via the formate hydrogenlyase (FHL) complex. In H2 producing marine vibrios, such as Vibrio tritonius and Vibrio porteresiae in the Porteresiae clade of the family Vibrionaceae, apparent but inefficient H2 production from Pyr has been observed. To elucidate the molecular mechanism of why this inefficient H2 production is observed in Pry-metabolized marine vibrio cells and how glycolysis affects those H2 productions of marine vibrios, the "Core Transcriptome" approach to find common gene expressions of those two major H2 producing Vibrio species in Pyr metabolism was first applied. In the Pyr-metabolized vibrio cells, genes for the "Phosphoenolpyruvate (PEP)-Pyruvate-Oxalate (PPO)" node, due to energy saving, and PhoB-, RhaR-, and DeoR-regulons were regulated. Interestingly, a gene responsible for oxalate/formate family antiporter was up-regulated in Pyr-metabolized cells compared to those of Glc-metabolized cells, which provides new insights into the uses of alternative formate exclusion mechanics due to energy deficiencies in Pyr-metabolized marine vibrios cells. We further discuss the contribution of the Embden-Meyerhof-Parnas (EMP) pathway to efficient H2 production in marine vibrios.

The TGFβ Induced MicroRNAome of the Trabecular Meshwork

Primary open-angle glaucoma (POAG) is a progressive optic neuropathy with a complex, multifactorial aetiology. Raised intraocular pressure (IOP) is the most important clinically modifiable risk factor for POAG. All current pharmacological agents target aqueous humour dynamics to lower IOP. Newer therapeutic agents are required as some patients with POAG show a limited therapeutic response or develop ocular and systemic side effects to topical medication. Elevated IOP in POAG results from cellular and molecular changes in the trabecular meshwork driven by increased levels of transforming growth factor β (TGFβ) in the anterior segment of the eye. Understanding how TGFβ affects both the structural and functional changes in the outflow pathway and IOP is required to develop new glaucoma therapies that target the molecular pathology in the trabecular meshwork. In this study, we evaluated the effects of TGF-β1 and -β2 treatment on miRNA expression in cultured human primary trabecular meshwork cells. Our findings are presented in terms of specific miRNAs (miRNA-centric), but given miRNAs work in networks to control cellular pathways and processes, a pathway-centric view of miRNA action is also reported. Evaluating TGFβ-responsive miRNA expression in trabecular meshwork cells will further our understanding of the important pathways and changes involved in the pathogenesis of glaucoma and could lead to the development of miRNAs as new therapeutic modalities in glaucoma.

Reelin marks cocaine-activated striatal ensembles, promotes neuronal excitability, and regulates cocaine reward

Drugs of abuse activate defined neuronal ensembles in brain reward structures such as the nucleus accumbens (NAc), which are thought to promote the enduring synaptic, circuit, and behavioral consequences of drug exposure. While the molecular and cellular effects arising from experience with drugs like cocaine are increasingly well understood, the mechanisms that sculpt NAc ensemble participation are largely unknown. Here, we leveraged unbiased single-nucleus transcriptional profiling to identify expression of the secreted glycoprotein Reelin (encoded by the Reln gene) as a marker of cocaine-activated neuronal ensembles within the rat NAc. Multiplexed in situ detection confirmed selective expression of the immediate early gene Fos in Reln+ neurons after cocaine experience, and also revealed enrichment of Reln mRNA in Drd1 + medium spiny neurons (MSNs) in both the rat and human brain. Using a novel CRISPR interference strategy enabling selective Reln knockdown in the adult NAc, we observed altered expression of genes linked to calcium signaling, emergence of a transcriptional trajectory consistent with loss of cocaine sensitivity, and a striking decrease in MSN intrinsic excitability. At the behavioral level, loss of Reln prevented cocaine locomotor sensitization, abolished cocaine place preference memory, and decreased cocaine self-administration behavior. Together, these results identify Reelin as a critical mechanistic link between ensemble participation and cocaine-induced behavioral adaptations.

Cancer-associated Histone H3 N-terminal arginine mutations disrupt PRC2 activity and impair differentiation

Dysregulated epigenetic states are a hallmark of cancer and often arise from genetic alterations in epigenetic regulators. This includes missense mutations in histones, which, together with associated DNA, form nucleosome core particles. However, the oncogenic mechanisms of most histone mutations are unknown. Here, we demonstrate that cancer-associated histone mutations at arginines in the histone H3 N-terminal tail disrupt repressive chromatin domains, alter gene regulation, and dysregulate differentiation. We find that histone H3R2C and R26C mutants reduce transcriptionally repressive H3K27me3. While H3K27me3 depletion in cells expressing these mutants is exclusively observed on the minor fraction of histone tails harboring the mutations, the same mutants recurrently disrupt broad H3K27me3 domains in the chromatin context, including near developmentally regulated promoters. H3K27me3 loss leads to de-repression of differentiation pathways, with concordant effects between H3R2 and H3R26 mutants despite different proximity to the PRC2 substrate, H3K27. Functionally, H3R26C-expressing mesenchymal progenitor cells and murine embryonic stem cell-derived teratomas demonstrate impaired differentiation. Collectively, these data show that cancer-associated H3 N-terminal arginine mutations reduce PRC2 activity and disrupt chromatin-dependent developmental functions, a cancer-relevant phenotype.

Mapping RNA-protein interactions with subcellular resolution using colocalization CLIP

RNA-binding proteins (RBPs) are essential for RNA metabolism and profoundly impact health and disease. The subcellular organization of RBP interaction networks with target RNAs remains largely unexplored. Here, we develop colocalization CLIP (coCLIP), a method that combines cross-linking and immunoprecipitation (CLIP) with proximity labeling, to explore in-depth the subcellular RNA interactions of the RBP human antigen R (HuR). Using this method, we uncover HuR's dynamic and location-specific interactions with RNA, revealing alterations in sequence preferences and interactions in the nucleus, cytosol, or stress granule (SG) compartments. We uncover HuR's unique binding preferences within SGs during arsenite stress, illuminating intricate interactions that conventional methodologies cannot capture. Overall, coCLIP provides a powerful method for revealing RBP-RNA interactions based on localization and lays the foundation for an advanced understanding of RBP models that incorporate subcellular location as a critical determinant of their functions.

Identification of functional enhancer variants associated with type I diabetes in CD4+ T cells

Type I diabetes is an autoimmune disease mediated by T-cell destruction of β cells in pancreatic islets. Currently, there is no known cure, and treatment consists of daily insulin injections. Genome-wide association studies and twin studies have indicated a strong genetic heritability for type I diabetes and implicated several genes. As most strongly associated variants are noncoding, there is still a lack of identification of functional and, therefore, likely causal variants. Given that many of these genetic variants reside in enhancer elements, we have tested 121 CD4+ T-cell enhancer variants associated with T1D. We found four to be functional through massively parallel reporter assays. Three of the enhancer variants weaken activity, while the fourth strengthens activity. We link these to their cognate genes using 3D genome architecture or eQTL data and validate them using CRISPR editing. Validated target genes include CLEC16A and SOCS1. While these genes have been previously implicated in type 1 diabetes and other autoimmune diseases, we show that enhancers controlling their expression harbor functional variants. These variants, therefore, may act as causal type 1 diabetic variants.

YABBY and diverged KNOX1 genes shape nodes and internodes in the stem

Plant stems comprise nodes and internodes that specialize in solute exchange and elongation. However, their boundaries are not well defined, and how these basic units arise remains elusive. In rice with clear nodes and internodes, we found that one subclade of class I knotted1-like homeobox (KNOX1) genes for shoot meristem indeterminacy restricts node differentiation and allows internode formation by repressing YABBY genes for leaf development and genes from another node-specific KNOX1 subclade. YABBYs promote nodal vascular differentiation and limit stem elongation. YABBY and node-specific KNOX1 genes specify the pulvinus, which further elaborates the nodal structure for gravitropism. Notably, this KNOX1 subclade organization is specific to seed plants. We propose that nodes and internodes are distinct domains specified by YABBY-KNOX1 cross-regulation that diverged in early seed plants.

Endocrine disruptor-induced epimutagenesis in vitro : Insight into molecular mechanisms

Endocrine disrupting chemicals (EDCs) such as bisphenol S (BPS) are xenobiotic compounds that can disrupt endocrine signaling following exposure due to steric similarities to endogenous hormones within the body. EDCs have been shown to induce disruptions in normal epigenetic programming (epimutations) that accompany dysregulation of normal gene expression patterns that appear to predispose disease states. Most interestingly, the prevalence of epimutations following exposure to many different EDCs often persists over multiple subsequent generations, even with no further exposure to the causative EDC. Many previous studies have described both the direct and prolonged effects of EDC exposure in animal models, but many questions remain about molecular mechanisms by which EDCs initially induce epimutations or contribute to the propagation of EDC-induced epimutations either within the exposed generation or to subsequent generations. Additional questions remain regarding the extent to which there may be differences in cell-type specific susceptibilities to various EDCs, and whether this susceptibility is correlative with expression of relevant hormone receptors and/or the location of relevant hormone response elements (HREs) in the genome. To address these questions, we exposed cultured mouse pluripotent (induced pluripotent stem [iPS]), somatic (Sertoli and granulosa), and germ (primordial germ cell like [PGCLC]) cells to BPS and measured changes in DNA methylation levels at the epigenomic level and gene expression at the transcriptomic level. We found that there was indeed a difference in cell-type specific susceptibility to EDC-induced epimutagenesis and that this susceptibility correlated with differential expression of relevant hormone receptors and, in many cases, tended to generate epimutations near relevant HREs within the genome. Additionally, however, we also found that BPS can induce epimutations in a cell type that does not express relevant receptors and in genomic regions that do not contain relevant HREs, suggesting that both canonical and non-canonical signaling mechanisms can be disrupted by BPS exposure. Most interestingly, we found that when iPS cells were exposed to BPS and then induced to differentiate into PGCLCs, the prevalence of epimutations and differentially expressed genes (DEGs) initially induced in the iPSCs was largely retained in the resulting PGCLCs, however, >90% of the specific epimutations and DEGs were not conserved but were rather replaced by novel epimutations and DEGs following the iPSC to PGCLC transition. These results are consistent with a unique concept that many EDC-induced epimutations may normally be corrected by germline and/or embryonic epigenetic reprogramming but that due to disruption of the underlying chromatin architecture induced by the EDC exposure, many novel epimutations may emerge during the reprogramming process as well. Thus, it appears that following exposure to a disruptive agent such as an EDC, a prevalence of epimutations may transcend epigenetic reprogramming even though most individual epimutations are not conserved during this process.

Human trophectoderm becomes multi-layered by internalization at the polar region

To implant in the uterus, mammalian embryos form blastocysts comprising trophectoderm (TE) surrounding an inner cell mass (ICM), confined to the polar region by the expanding blastocoel. The mode of implantation varies between species. Murine embryos maintain a single layered TE until they implant in the characteristic thick deciduum, whereas human blastocysts attach via polar TE directly to the uterine wall. Using immunofluorescence (IF) of rapidly isolated ICMs, blockade of RNA and protein synthesis in whole embryos, or 3D visualization of immunostained embryos, we provide evidence of multi-layering in human polar TE before implantation. This may be required for rapid uterine invasion to secure the developing human embryo and initiate formation of the placenta. Using sequential fluorescent labeling, we demonstrate that the majority of inner TE in human blastocysts arises from existing outer cells, with no evidence of conversion from the ICM in the context of the intact embryo.

Autotrophic adaptive laboratory evolution of the acetogen Clostridium autoethanogenum delivers the gas-fermenting strain LAbrini with superior growth, products, and robustness

Microbes able to convert gaseous one-carbon (C1) waste feedstocks are increasingly important to transition to the sustainable production of renewable chemicals and fuels. Acetogens are interesting biocatalysts since gas fermentation using Clostridium autoethanogenum has been commercialised. However, most acetogen strains need complex nutrients, display slow growth, and are not robust for bioreactor fermentations. In this work, we used three different and independent adaptive laboratory evolution (ALE) strategies to evolve the wild-type C. autoethanogenum to grow faster, without yeast extract and to be robust in operating continuous bioreactor cultures. Multiple evolved strains with improved phenotypes were isolated on minimal media with one strain, named "LAbrini", exhibiting superior performance regarding the maximum specific growth rate, product profile, and robustness in continuous cultures. Whole-genome sequencing of the evolved strains identified 25 mutations. Of particular interest are two genes that acquired seven different mutations across the three ALE strategies, potentially as a result of convergent evolution. Reverse genetic engineering of mutations in potentially sporulation-related genes CLAU_3129 (spo0A) and CLAU_1957 recovered all three superior features of our ALE strains through triggering significant proteomic rearrangements. This work provides a robust C. autoethanogenum strain "LAbrini" to accelerate phenotyping and genetic engineering and to better understand acetogen metabolism.

Time-course swRNA-seq uncovers a hierarchical gene regulatory network in controlling the response-repair-remodeling after wounding

Wounding initiates intricate responses crucial for tissue repair and regeneration. Yet, the gene regulatory networks governing wound healing remain poorly understood. Here, employing single-worm RNA sequencing (swRNA-seq) across 12 time-points, we delineated a three-stage wound repair process in C. elegans: response, repair, and remodeling. Integrating diverse datasets, we constructed a dynamic regulatory network comprising 241 transcription regulators and their inferred targets. We identified potentially seven autoregulatory TFs and five cross-autoregulatory loops involving pqm-1 and jun-1. We revealed that TFs might interact with chromatin factors and form TF-TF combinatory modules via intrinsically disordered regions to enhance response robustness. We experimentally validated six regulators functioning in transcriptional and translocation-dependent manners. Notably, nhr-76, daf-16, nhr-84, and oef-1 are potentially required for efficient repair, while elt-2 may act as an inhibitor. These findings elucidate transcriptional responses and hierarchical regulatory networks during C. elegans wound repair, shedding light on mechanisms underlying tissue repair and regeneration.

Circuit-Wide Gene Network Analysis Reveals Sex-Specific Roles for Phosphodiesterase 1b in Cocaine Addiction

Cocaine use disorder is a significant public health issue without an effective pharmacological treatment. Successful treatments are hindered in part by an incomplete understanding of the molecular mechanisms that underlie long-lasting maladaptive plasticity and addiction-like behaviors. Here, we leverage a large RNA sequencing dataset to generate gene coexpression networks across six interconnected regions of the brain's reward circuitry from mice that underwent saline or cocaine self-administration. We identify phosphodiesterase 1b (Pde1b), a Ca2+/calmodulin-dependent enzyme that increases cAMP and cGMP hydrolysis, as a central hub gene within a nucleus accumbens (NAc) gene module that was bioinformatically associated with addiction-like behavior. Chronic cocaine exposure increases Pde1b expression in NAc D2 medium spiny neurons (MSNs) in male but not female mice. Viral-mediated Pde1b overexpression in NAc reduces cocaine self-administration in female rats but increases seeking in both sexes. In female mice, overexpressing Pde1b in D1 MSNs attenuates the locomotor response to cocaine, with the opposite effect in D2 MSNs. Overexpressing Pde1b in D1/D2 MSNs had no effect on the locomotor response to cocaine in male mice. At the electrophysiological level, Pde1b overexpression reduces sEPSC frequency in D1 MSNs and regulates the excitability of NAc MSNs. Lastly, Pde1b overexpression significantly reduced the number of differentially expressed genes (DEGs) in NAc following chronic cocaine, with discordant effects on gene transcription between sexes. Together, we identify novel gene modules across the brain's reward circuitry associated with addiction-like behavior and explore the role of Pde1b in regulating the molecular, cellular, and behavioral responses to cocaine.

Myocardin-Related Transcription Factor Mediates Epithelial Fibrogenesis in Polycystic Kidney Disease

Polycystic kidney disease (PKD) is characterized by extensive cyst formation and progressive fibrosis. However, the molecular mechanisms whereby the loss/loss-of-function of Polycystin 1 or 2 (PC1/2) provokes fibrosis are largely unknown. The small GTPase RhoA has been recently implicated in cystogenesis, and we identified the RhoA/cytoskeleton/myocardin-related transcription factor (MRTF) pathway as an emerging mediator of epithelium-induced fibrogenesis. Therefore, we hypothesized that MRTF is activated by PC1/2 loss and plays a critical role in the fibrogenic reprogramming of the epithelium. The loss of PC1 or PC2, induced by siRNA in vitro, activated RhoA and caused cytoskeletal remodeling and robust nuclear MRTF translocation and overexpression. These phenomena were also manifested in PKD1 (RC/RC) and PKD2 (WS25/-) mice, with MRTF translocation and overexpression occurring predominantly in dilated tubules and the cyst-lining epithelium, respectively. In epithelial cells, a large cohort of PC1/PC2 downregulation-induced genes was MRTF-dependent, including cytoskeletal, integrin-related, and matricellular/fibrogenic proteins. Epithelial MRTF was necessary for the paracrine priming of the fibroblast-myofibroblast transition. Thus, MRTF acts as a prime inducer of epithelial fibrogenesis in PKD. We propose that RhoA is a common upstream inducer of both histological hallmarks of PKD: cystogenesis and fibrosis.

The Neo-Open Reading Frame Peptides That Comprise the Tumor Framome Are a Rich Source of Neoantigens for Cancer Immunotherapy

Identification of immunogenic cancer neoantigens as targets for therapy is challenging. Here, we integrate the whole-genome and long-read transcript sequencing of cancers to identify the collection of neo-open reading frame peptides (NOP) expressed in tumors. We termed this collection of NOPs the tumor framome. NOPs represent tumor-specific peptides that are different from wild-type proteins and may be strongly immunogenic. We describe a class of hidden NOPs that derive from structural genomic variants involving an upstream protein coding gene driving expression and translation of noncoding regions of the genome downstream of a rearrangement breakpoint, i.e., where no gene annotation or evidence for transcription exists. The entire collection of NOPs represents a vast number of possible neoantigens particularly in tumors with many structural genomic variants and a low number of missense mutations. We show that NOPs are immunogenic and epitopes derived from NOPs can bind to MHC class I molecules. Finally, we provide evidence for the presence of memory T cells specific for hidden NOPs in peripheral blood from a patient with lung cancer. This work highlights NOPs as a major source of possible neoantigens for personalized cancer immunotherapy and provides a rationale for analyzing the complete cancer genome and transcriptome as a basis for the detection of NOPs.

TSC2 loss in neural progenitor cells suppresses translation of ASD/NDD-associated transcripts in an mTORC1- and MNK1/2-reversible fashion

Tuberous sclerosis complex (TSC), an inherited neurodevelopmental (ND) disorder with frequent manifestations of epilepsy and autism spectrum disorder (ASD). TSC is caused by mutations in TSC1 or TSC2 tumor suppressor genes, with encoded proteins hamartin/TSC1 and tuberin/TSC2 forming a functional complex inhibiting mechanistic target of rapamycin complex-1 (mTORC1) signaling, leading to FDA-approved allosteric mTORC1-selective rapamycin analogs for TSC tumors. Rapalogs are effective for TSC-associated hamartomas, however, they are not effective for treating ND manifestations. mTORC1 signaling plays an essential role in protein synthesis through mTORC1-eIF4F and MNK-eIF4E-mediated mRNA translation. Further, the effects on mRNA translation by specific mTORC1 and MNK inhibitors such as RMC-6272 and eFT-508 in TSC have never been explored. Here, employing CRISPR-modified, isogenic TSC2 patient-derived neural progenitor cells (NPCs), we have examined mRNA translation upon loss of TSC2 . Our results reveal dysregulated translation in TSC2 -Null NPCs, which significantly overlap with the translatome from TSC1 -Null NPCs, which we reported recently. Most notably, numerous non-monogenic ASD-NDD- and epilepsy-associated genes identified in patients harboring putative loss-of-function mutations, including protein truncating, or damaging missense variants, were translationally suppressed in TSC2 -Null NPCs, and their translation were reversed upon RMC-6272 or eFT-508 treatment. Our study here establishes the importance of mTORC1-eIF4F and MNK-eIF4E-mediated mRNA translation in TSC, ASD and other neurodevelopmental disorders and lay the groundwork for evaluating drugs in clinical development that target these pathways as a treatment strategy for TSC as well as ASD/NDD.

Nextflow vs. plain bash: different approaches to the parallelization of SNP calling from the whole genome sequence data

This study compared computational approaches to parallelization of an SNP calling workflow. The data comprised DNA from five Holstein-Friesian cows sequenced with the Illumina platform. The pipeline consisted of quality control, alignment to the reference genome, post-alignment, and SNP calling. Three approaches to parallelization were compared: (i) a plain Bash script in which a pipeline for each cow was executed as separate processes invoked at the same time, (ii) a Bash script wrapped in a single Nextflow process and (iii) a Nextflow script with each component of the pipeline defined as a separate process. The results demonstrated that on average, the multi-process Nextflow script performed 15-27% faster depending on the number of assigned threads, with the biggest execution time advantage over the plain Bash approach observed with 10 threads. In terms of RAM usage, the most substantial variation was observed for the multi-process Nextflow, for which it increased with the number of assigned threads, while RAM consumption of the other setups did not depend much on the number of threads assigned for computations. Due to intermediate and log files generated, disk usage was markedly higher for the multi-process Nextflow than for the plain Bash and for the single-process Nextflow.

Transcriptomic analysis and fusion gene identifications of midbody remnants released from colorectal cancer cells reveals they are molecularly distinct from exosomes and microparticles

Previously, we reported that human primary (SW480) and metastatic (SW620) colorectal (CRC) cells release three classes of membrane-encapsulated extracellular vesicles (EVs); midbody remnants (MBRs), exosomes (Exos), and microparticles (MPs). We reported that MBRs were molecularly distinct at the protein level. To gain further biochemical insights into MBRs, Exos, and MPs and their emerging role in CRC, we performed, and report here, for the first time, a comprehensive transcriptome and long noncoding RNA sequencing analysis and fusion gene identification of these three EV classes using the next-generation RNA sequencing technique. Differential transcript expression analysis revealed that MBRs have a distinct transcriptomic profile compared to Exos and MPs with a high enrichment of mitochondrial transcripts lncRNA/pseudogene transcripts that are predicted to bind to ribonucleoprotein complexes, spliceosome, and RNA/stress granule proteins. A salient finding from this study is a high enrichment of several fusion genes in MBRs compared to Exos, MPs, and cell lysates from their parental cells such as MSH2 (gene encoded DNA mismatch repair protein MSH2). This suggests potential EV-liquid biopsy targets for cancer detection. Importantly, the expression of cancer progression-related transcripts found in EV classes derived from SW480 (EGFR) and SW620 (MET and MACCA1) cell lines reflects their parental cell types. Our study is the report of RNA and fusion gene compositions within MBRs (including Exos and MPs) that could have an impact on EV functionality in cancer progression and detection using EV-based RNA/ fusion gene candidates for cancer biomarkers.

Defining the relationship between cellular and extracellular vesicle (EV) content in breast cancer via an integrative multi-omic analysis

Much recent research has been dedicated to exploring the utility of extracellular vesicles (EVs) as circulating disease biomarkers. Underpinning this work is the assumption that the molecular cargo of EVs directly reflects the originating cell. Few attempts have been made, however, to empirically validate this on the -omic level. To this end, we have performed an integrative multi-omic analysis of a panel of breast cancer cell lines and corresponding EVs. Whole transcriptome analysis validated that the cellular transcriptome remained stable when cultured cells are transitioned to low serum or serum-free medium for EV collection. Transcriptomic profiling of the isolated EVs indicated a positive correlation between transcript levels in cells and EVs, including disease-associated transcripts. Analysis of the EV proteome verified that HER2 protein is present in EVs, however neither the estrogen (ER) nor progesterone (PR) receptor proteins are detected regardless of cellular expression. Using multivariate analysis, we derived an EV protein signature to infer cellular patterns of ER and HER2 expression, though the ER protein could not be directly detected. Integrative analyses affirmed that the EV proteome and transcriptome captured key phenotypic hallmarks of the originating cells, supporting the potential of EVs for non-invasive monitoring of breast cancers.

Differential impact of sex on regulation of skeletal muscle mitochondrial function and protein homeostasis by hypoxia-inducible factor-1α in normoxia

Hypoxia-inducible factor (HIF)-1α is continuously synthesized and degraded in normoxia. During hypoxia, HIF1α stabilization restricts cellular/mitochondrial oxygen utilization. Cellular stressors can stabilize HIF1α even during normoxia. However, less is known about HIF1α function(s) and sex-specific effects during normoxia in the basal state. Since skeletal muscle is the largest protein store in mammals and protein homeostasis has high energy demands, we determined HIF1α function at baseline during normoxia in skeletal muscle. Untargeted multiomics data analyses were followed by experimental validation in differentiated murine myotubes with loss/gain of function and skeletal muscle from mice without/with post-natal muscle-specific Hif1a deletion (Hif1amsd). Mitochondrial oxygen consumption studies using substrate, uncoupler, inhibitor, titration protocols; targeted metabolite quantification by gas chromatography-mass spectrometry; and post-mitotic senescence markers using biochemical assays were performed. Multiomics analyses showed enrichment in mitochondrial and cell cycle regulatory pathways in Hif1a deleted cells/tissue. Experimentally, mitochondrial oxidative functions and ATP content were higher with less mitochondrial free radical generation with Hif1a deletion. Deletion of Hif1a also resulted in higher concentrations of TCA cycle intermediates and HIF2α proteins in myotubes. Overall responses to Hif1amsd were similar in male and female mice, but changes in complex II function, maximum respiration, Sirt3 and HIF1β protein expression and muscle fibre diameter were sex-dependent. Adaptive responses to hypoxia are mediated by stabilization of constantly synthesized HIF1α. Despite rapid degradation, the presence of HIF1α during normoxia contributes to lower mitochondrial oxidative efficiency and greater post-mitotic senescence in skeletal muscle. In vivo responses to HIF1α in skeletal muscle were differentially impacted by sex. KEY POINTS: Hypoxia-inducible factor -1α (HIF1α), a critical transcription factor, undergoes continuous synthesis and proteolysis, enabling rapid adaptive responses to hypoxia by reducing mitochondrial oxygen consumption. In mammals, skeletal muscle is the largest protein store which is determined by a balance between protein synthesis and breakdown and is sensitive to mitochondrial oxidative function. To investigate the functional consequences of transient HIF1α expression during normoxia in the basal state, myotubes and skeletal muscle from male and female mice with HIF1α knockout were studied using complementary multiomics, biochemical and metabolite assays. HIF1α knockout altered the electron transport chain, mitochondrial oxidative function, signalling molecules for protein homeostasis, and post-mitotic senescence markers, some of which were differentially impacted by sex. The cost of rapid adaptive responses mediated by HIF1α is lower mitochondrial oxidative efficiency and post-mitotic senescence during normoxia.

Grainyhead-like-2, an epithelial master programmer, promotes interferon induction and suppresses breast cancer recurrence

Type-I and -III interferons play a central role in immune rejection of pathogens and tumors, thus promoting immunogenicity and suppressing tumor recurrence. Double strand RNA is an important ligand that stimulates tumor immunity via interferon responses. Differentiation of embryonic stem cells to pluripotent epithelial cells activates the interferon response during development, raising the question of whether epithelial vs. mesenchymal gene signatures in cancer potentially regulate the interferon pathway as well. Here, using genomics and signaling approaches, we show that Grainyhead-like-2 (GRHL2), a master programmer of epithelial cell identity, promotes type-I and -III interferon responses to double-strand RNA. GRHL2 enhanced the activation of IRF3 and relA/NF-kB and the expression of IRF1; a functional GRHL2 binding site in the IFNL1 promoter was also identified. Moreover, time to recurrence in breast cancer correlated positively with GRHL2 protein expression, indicating that GRHL2 is a tumor recurrence suppressor, consistent with its enhancement of interferon responses. These observations demonstrate that epithelial cell identity supports interferon responses in the context of cancer.

Differential in vivo roles of Mpl cytoplasmic tyrosine residues in murine hematopoiesis and myeloproliferative disease

Thrombopoietin (Tpo), which binds to its specific receptor, the Mpl protein, is the major cytokine regulator of megakaryopoiesis and circulating platelet number. Tpo binding to Mpl triggers activation of Janus kinase 2 (Jak2) and phosphorylation of the receptor, as well as activation of several intracellular signalling cascades that mediate cellular responses. Three tyrosine (Y) residues in the C-terminal region of the Mpl intracellular domain have been implicated as sites of phosphorylation required for regulation of major Tpo-stimulated signalling pathways: Mpl-Y565, Mpl-Y599 and Mpl-Y604. Here, we have introduced mutations in the mouse germline and report a consistent physiological requirement for Mpl-Y599, mutation of which resulted in thrombocytopenia, deficient megakaryopoiesis, low hematopoietic stem cell (HSC) number and function, and attenuated responses to myelosuppression. We further show that in models of myeloproliferative neoplasms (MPN), where Mpl is required for pathogenesis, thrombocytosis was dependent on intact Mpl-Y599. In contrast, Mpl-Y565 was required for negative regulation of Tpo responses; mutation of this residue resulted in excess megakaryopoiesis at steady-state and in response to myelosuppression, and exacerbated thrombocytosis associated with MPN.

A WOX homolog disrupted by a transposon led to the loss of spines and contributed to the domestication of lettuce

The loss of spines is one of the most important domestication traits for lettuce (Lactuca sativa). However, the genetics and regulation of spine development in lettuce remain unclear. We examined the genetics of spines in lettuce using a segregating population derived from a cross between cultivated and wild lettuce (Lactuca serriola). A gene encoding WUSCHEL-related homeobox transcription factor, named as WOX-SPINE1 (WS1), was identified as the candidate gene controlling the spine development in lettuce, and its function on spines was verified. A CACTA transposon was found to be inserted into the first exon of the ws1 allele, knocking out its function and leading to the lack of spines in cultivated lettuce. All lettuce cultivars investigated have the nonfunctional ws1 gene, and a selection sweep was found at the WS1 locus, suggesting its important role in lettuce domestication. The expression levels of WS1 were associated with the density of spines among different accessions of wild lettuce. At least two independent loss-of-function mutations in the ws1 gene caused the loss of spines in wild lettuce. These findings provide new insights into the development of spines and facilitate the exploitation of wild genetic resources in future lettuce breeding programs.

Dynamic Changes in Lymphocyte Populations Establish Zebrafish as a Thymic Involution Model

The thymus is the site of T lymphocyte development and T cell education to recognize foreign, but not self, Ags. B cells also reside and develop in the thymus, although their functions are less clear. During "thymic involution," a process of lymphoid atrophy and adipose replacement linked to sexual maturation, thymocytes decline. However, thymic B cells decrease far less than T cells, such that B cells comprise ∼1% of human neonatal thymocytes but up to ∼10% in adults. All jawed vertebrates possess a thymus, and we and others have shown zebrafish (Danio rerio) also have thymic B cells. In this article, we investigated the precise identities of zebrafish thymic T and B cells and how they change with involution. We assessed the timing and specific details of zebrafish thymic involution using multiple lymphocyte-specific, fluorophore-labeled transgenic lines, quantifying the changes in thymic T- and B-lymphocytes pre- versus postinvolution. Our results prove that, as in humans, zebrafish thymic B cells increase relative to T cells postinvolution. We also performed RNA sequencing on D. rerio thymic and marrow lymphocytes of four novel double-transgenic lines, identifying distinct populations of immature T and B cells. Collectively, this is, to our knowledge, the first comprehensive analysis of zebrafish thymic involution, demonstrating its similarity to human involution and establishing the highly genetically manipulatable zebrafish model as a template for involution studies.

Multiomic analyses uncover immunological signatures in acute and chronic coronary syndromes

Acute and chronic coronary syndromes (ACS and CCS) are leading causes of mortality. Inflammation is considered a key pathogenic driver of these diseases, but the underlying immune states and their clinical implications remain poorly understood. Multiomic factor analysis (MOFA) allows unsupervised data exploration across multiple data types, identifying major axes of variation and associating these with underlying molecular processes. We hypothesized that applying MOFA to multiomic data obtained from blood might uncover hidden sources of variance and provide pathophysiological insights linked to clinical needs. Here we compile a longitudinal multiomic dataset of the systemic immune landscape in both ACS and CCS (n = 62 patients in total, n = 15 women and n = 47 men) and validate this in an external cohort (n = 55 patients in total, n = 11 women and n = 44 men). MOFA reveals multicellular immune signatures characterized by distinct monocyte, natural killer and T cell substates and immune-communication pathways that explain a large proportion of inter-patient variance. We also identify specific factors that reflect disease state or associate with treatment outcome in ACS as measured using left ventricular ejection fraction. Hence, this study provides proof-of-concept evidence for the ability of MOFA to uncover multicellular immune programs in cardiovascular disease, opening new directions for mechanistic, biomarker and therapeutic studies.

m6Aexpress-enet: Predicting the regulatory expression m6A sites by an enet-regularization negative binomial regression model

As the most abundant mRNA modification, m6A controls and influences many aspects of mRNA metabolism including the mRNA stability and degradation. However, the role of specific m6A sites in regulating gene expression still remains unclear. In additional, the multicollinearity problem caused by the correlation of methylation level of multiple m6A sites in each gene could influence the prediction performance. To address the above challenges, we propose an elastic-net regularized negative binomial regression model (called m6Aexpress-enet) to predict which m6A site could potentially regulate its gene expression. Comprehensive evaluations on simulated datasets demonstrate that m6Aexpress-enet could achieve the top prediction performance. Applying m6Aexpress-enet on real MeRIP-seq data from human lymphoblastoid cell lines, we have uncovered the complex regulatory pattern of predicted m6A sites and their unique enrichment pathway of the constructed co-methylation modules. m6Aexpress-enet proves itself as a powerful tool to enable biologists to discover the mechanism of m6A regulatory gene expression. Furthermore, the source code and the step-by-step implementation of m6Aexpress-enet is freely accessed at https://github.com/tengzhangs/m6Aexpress-enet.

Development and validation of a radiosensitivity model to evaluate radiotherapy benefits in pan-cancer

Radiotherapy, one of the most fundamental cancer treatments, is confronted with the dilemma of treatment failure due to radioresistance. To predict the radiosensitivity and improve tumor treatment efficiency in pan-cancer, we developed a model called Radiation Intrinsic Sensitivity Evaluation (RISE). The RISE model was built using cell line-based mRNA sequencing data from five tumor types with varying radiation sensitivity. Through four cell-derived datasets, two public tissue-derived cohorts, and one local cohort of 42 nasopharyngeal carcinoma patients, we demonstrated that RISE could effectively predict the level of radiation sensitivity (area under the ROC curve [AUC] from 0.666 to 1 across different datasets). After the verification by the colony formation assay and flow cytometric analysis of apoptosis, our four well-established radioresistant cell models successfully proved higher RISE values in radioresistant cells by RT-qPCR experiments. We also explored the prognostic value of RISE in five independent TCGA cohorts consisting of 1137 patients who received radiation therapy and found that RISE was an independent adverse prognostic factor (pooled multivariate Cox regression hazard ratio [HR]: 1.84, 95% CI 1.39-2.42; p < 0.01). RISE showed a promising ability to evaluate the radiotherapy benefit while predicting the prognosis of cancer patients, enabling clinicians to make individualized radiotherapy strategies in the future and improve the success rate of radiotherapy.

A GGC-repeat expansion in ZFHX3 encoding polyglycine causes spinocerebellar ataxia type 4 and impairs autophagy

Despite linkage to chromosome 16q in 1996, the mutation causing spinocerebellar ataxia type 4 (SCA4), a late-onset sensory and cerebellar ataxia, remained unknown. Here, using long-read single-strand whole-genome sequencing (LR-GS), we identified a heterozygous GGC-repeat expansion in a large Utah pedigree encoding polyglycine (polyG) in zinc finger homeobox protein 3 (ZFHX3), also known as AT-binding transcription factor 1 (ATBF1). We queried 6,495 genome sequencing datasets and identified the repeat expansion in seven additional pedigrees. Ultrarare DNA variants near the repeat expansion indicate a common distant founder event in Sweden. Intranuclear ZFHX3-p62-ubiquitin aggregates were abundant in SCA4 basis pontis neurons. In fibroblasts and induced pluripotent stem cells, the GGC expansion led to increased ZFHX3 protein levels and abnormal autophagy, which were normalized with small interfering RNA-mediated ZFHX3 knockdown in both cell types. Improving autophagy points to a therapeutic avenue for this novel polyG disease. The coding GGC-repeat expansion in an extremely G+C-rich region was not detectable by short-read whole-exome sequencing, which demonstrates the power of LR-GS for variant discovery.

M2 macrophage-derived TGF-β induces age-associated loss of adipogenesis through progenitor cell senescence

OBJECTIVES

Adipose tissue is an endocrine and energy storage organ composed of several different cell types, including mature adipocytes, stromal cells, endothelial cells, and a variety of immune cells. Adipose tissue aging contributes to the pathogenesis of metabolic dysfunction and is likely induced by crosstalk between adipose progenitor cells (APCs) and immune cells, but the underlying molecular mechanisms remain largely unknown. In this study, we revealed the biological role of p16high senescent APCs, and investigated the crosstalk between each cell type in the aged white adipose tissue.

METHODS

We performed the single-cell RNA sequencing (scRNA-seq) analysis on the p16high adipose cells sorted from aged p16-CreERT2/Rosa26-LSL-tdTomato mice. We also performed the time serial analysis on the age-dependent bulk RNA-seq datasets of human and mouse white adipose tissues to infer the transcriptome alteration of adipogenic potential within aging.

RESULTS

We show that M2 macrophage-derived TGF-β induces APCs senescence which impairs adipogenesis in vivo. p16high senescent APCs increase with age and show loss of adipogenic potential. The ligand-receptor interaction analysis reveals that M2 macrophages are the donors for TGF-β and the senescent APCs are the recipients. Indeed, treatment of APCs with TGF-β1 induces senescent phenotypes through mitochondrial ROS-mediated DNA damage in vitro. TGF-β1 injection into gonadal white adipose tissue (gWAT) suppresses adipogenic potential and induces fibrotic genes as well as p16 in APCs. A gWAT atrophy is observed in cancer cachexia by APCs senescence, whose induction appeared to be independent of TGF-β induction.

CONCLUSIONS

Our results suggest that M2 macrophage-derived TGF-β induces age-related lipodystrophy by APCs senescence. The TGF-β treatment induced DNA damage, mitochondrial ROS, and finally cellular senescence in APCs.

Biosensor capability of the endometrium is mediated in part, by altered miRNA cargo from conceptus-derived extracellular vesicles

We tested the hypothesis that the biosensor capability of the endometrium is mediated in part, by the effect of different cargo contained in the extracellular vesicles secreted by the conceptus during the peri-implantation period of pregnancy. We transferred Bos taurus taurus embryos of different origin, in vivo (high developmental potential (IV)), in vitro (intermediate developmental potential (IVF)), or cloned (low developmental potential (NT)), into Bos taurus indicus recipients. Extracellular vesicles (EVs) recovered from Day 16 conceptus-conditioned medium were characterized and their microRNA (miRNA) cargo sequenced alongside RNA sequencing of their respective endometria. There were substantial differences in the endometrial response to in vivo versus in vitro and in vivo versus cloned conceptuses (1153 and 334DEGs respectively) with limited differences between in vitro Vs cloned conceptuses (36 DEGs). The miRNA cargo contained in conceptus-derived EVs was similar between all three groups (426 miRNA in common). Only 8 miRNAs were different between in vivo and cloned conceptuses, while only 6 miRNAs were different between in vivo and in vitro-derived conceptuses. Treatment of endometrial epithelial cells with mimic or inhibitors for miR-128 and miR-1298 changed the proteomic content of target cells (96 and 85, respectively) of which mRNAs are altered in the endometrium in vivo (PLXDC2, COPG1, HSPA12A, MCM5, TBL1XR1, and TTF). In conclusion, we have determined that the biosensor capability of the endometrium is mediated in part, by its response to different EVs miRNA cargo produced by the conceptus during the peri-implantation period of pregnancy.

Vaginal microbes alter epithelial transcriptome and induce epigenomic modifications providing insight into mechanisms for susceptibility to adverse reproductive outcomes

The cervicovaginal microbiome is highly associated with women's health, with microbial communities dominated by Lactobacillus species considered optimal. Conversely, a lack of lactobacilli and a high abundance of strict and facultative anaerobes, including Gardnerella vaginalis, have been associated with adverse reproductive outcomes. However, how host-microbial interactions alter specific molecular pathways and impact cervical and vaginal epithelial function remains unclear. Using RNA-sequencing, we characterized the in vitro cervicovaginal epithelial transcriptional response to different vaginal bacteria and their culture supernatants. We showed that G. vaginalis upregulates genes associated with an activated innate immune response. Unexpectedly, G. vaginalis specifically induced inflammasome pathways through activation of NLRP3-mediated increases in caspase-1, IL-1β and cell death, while live L. crispatus had minimal transcriptomic changes on epithelial cells. L. crispatus culture supernatants resulted in a shift in the epigenomic landscape of cervical epithelial cells that was confirmed by ATAC-sequencing showing reduced chromatin accessibility. This study reveals new insights into host-microbe interactions in the lower reproductive tract and suggests potential therapeutic strategies leveraging the vaginal microbiome to improve reproductive health.

BMP9 is a key player in endothelial identity and its loss is sufficient to induce arteriovenous malformations

AIMS

BMP9 is a high affinity ligand of ALK1 and endoglin receptors that are mutated in the rare genetic vascular disorder hereditary hemorrhagic telangiectasia (HHT). We have previously shown that loss of Bmp9 in the 129/Ola genetic background leads to spontaneous liver fibrosis via capillarization of liver sinusoidal endothelial cells (LSEC) and kidney lesions. We aimed to decipher the molecular mechanisms downstream of BMP9 to better characterize its role in vascular homeostasis in different organs.

METHODS AND RESULTS

For this, we performed an RNA-seq analysis on LSEC from adult WT and Bmp9-KO mice and identified over 2000 differentially expressed genes. Gene ontology analysis showed that Bmp9 deletion led to a decrease in BMP and Notch signalling, but also LSEC capillary identity while increasing their cell cycle. The gene ontology term 'glomerulus development' was also negatively enriched in Bmp9-KO mice vs. WT supporting a role for BMP9 in kidney vascularization. Through different imaging approaches (electron microscopy, immunostainings), we found that loss of Bmp9 led to vascular enlargement of the glomeruli capillaries associated with alteration of podocytes. Importantly, we also showed for the first time that the loss of Bmp9 led to spontaneous arteriovenous malformations (AVMs) in the liver, gastrointestinal tract, and uterus.

CONCLUSION

Altogether, these results demonstrate that BMP9 plays an important role in vascular quiescence both locally in the liver by regulating endothelial capillary differentiation markers and cell cycle but also at distance in many organs via its presence in the circulation. It also reveals that loss of Bmp9 is sufficient to induce spontaneous AVMs, supporting a key role for BMP9 in the pathogenesis of HHT.